基于16S和COⅠ基因探討巖扇貝在扇貝科中的分類地位

曹善茂 朱麗潔 聶鴻濤

摘要：采用PCR擴增和序列測定等技術(shù)，對巖扇貝（Crassadoma gigantea）線粒體DNA 16S rRNA和CO Ⅰ基因片段進行初步研究，得到16S rRNA基因片段的大小在613.0～634.0 bp之間，A、T（U）、G、C的平均含量分別為26.5%、29.7%、24.5%、19.3%;CO Ⅰ基因片段的大小為660 bp，堿基A、T、G、C的平均含量分別為17.2%、40.3%、27.2%、15.3%。16S基因在巖扇貝中偏好使用GUU、UUG、AUU、ACG、UAU、CAA、GAG、UUU、UCU、GAU、GGU、GAG、AGG和GGA等25個密碼子。COⅠ系統(tǒng)進化樹表明，與巖扇貝親緣關(guān)系最近的是柿孔扇貝（Azumapecten farreri），其次是冰島扇貝（Chlamys islandica）。16S系統(tǒng)進化樹親緣關(guān)系由近到遠依次是北美扇貝（Patinopecten caurinus）、蝦夷扇貝（Mizuhopecten yessoensis）、Chlamys islandica。

關(guān)鍵詞：巖扇貝;CO Ⅰ基因;16S rRNA;密碼子偏好性;系統(tǒng)進化樹分析

中圖分類號： S917.4 ?文獻標志碼： A ?文章編號：1002-1302（2019）07-0028-04

巖扇貝（Crassadoma gigantea）隸屬于軟體動物門、瓣鰓綱、珍珠貝目、扇貝科。原產(chǎn)于北美洲太平洋沿岸海域，廣泛分布于美國的阿拉斯加，加利福尼亞半島以及墨西哥灣附近。該貝具有閉殼肌大、肉質(zhì)鮮美、風味獨特、生長速度快、對環(huán)境適應能力強等優(yōu)點，2013年大連海洋大學科研人員從加拿大溫哥華島引進巖扇貝成貝，在遼寧省沿海各地進行人工育苗和人工養(yǎng)殖，并對其進行了人工育苗技術(shù)和幼貝攝食生理等研究[1-3]。國外學者對巖扇貝的營養(yǎng)成分、口味、冷藏肉質(zhì)的穩(wěn)定性、野生種群的生長速率、排精產(chǎn)卵機制、幼貝的餌料配比、附著生理、胚胎發(fā)育等方面進行了研究[4-9]。日常巖扇貝通過巖扇貝殼型、顏色等形態(tài)特征進行分類，比較而言應用分子生物學的分析方法對巖扇貝的分類及鑒定較少。Bernard建立了包含必須以固著而生的巖扇貝的巖扇貝屬（Crassadoma）[10]。Waller擴展了巖扇貝屬的范圍，包含其他固著和非固著的Hinnites種，并創(chuàng)造了包括Caribachlamys屬在內(nèi)的Crassadomini[11]。可見巖扇貝在扇貝科下的親緣關(guān)系并不明確，Saavedra等研究了巖扇貝與太平洋部分扇貝和加勒比海扇貝的進化關(guān)系[12]。本研究在其基礎(chǔ)上增加國內(nèi)經(jīng)濟貝類，如櫛孔扇貝、蝦夷扇貝等，應用DNA條形碼技術(shù)進行進一步的研究。DNA條形碼是利用足夠變異的標準化短基因片段構(gòu)建物種的鑒別體系，具有對物種進行快速、準確鑒定的新生物身份識別的優(yōu)勢[13-15]。本研究采用PCR擴增和序列測定技術(shù)，對巖扇貝線粒體DNA 16S rRNA和CO Ⅰ基因片段進行初步研究，進一步證實巖扇貝在扇貝科中的分類地位及親緣關(guān)系。

1 材料與方法

1.1 試驗材料

巖扇貝于2016年11月取自遼寧省旅順河口養(yǎng)殖基地，鑒定后活體解剖取扇貝閉殼肌于液氮冷凍后轉(zhuǎn)-80 ℃超低溫冰箱保存?zhèn)溆谩?/p>

1.2 方法

1.2.1 基因組DNA提取和檢測 隨機選取巖扇貝2個，取大約30 mg閉殼肌組織，盡量剪碎，用DNA提取試劑盒提取基因組DNA，并溶解于100 μL Tris-EDTA緩沖buffer（TE溶液）中，用紫外分光光度計測定其濃度和純度，然后將基因組DNA稀釋至50 ng/μL，4 ℃ 下保存?zhèn)溆谩?/p>

1.2.2 PCR擴增及測序 本試驗使用CO Ⅰ和16S引物序列進行PCR擴增（表1），所有PCR反應均采用預混合酶試劑（premix），配制 10 μL 體系，在PTC-100型PCR儀上進行反應。PCR反應程序為94 ℃預變性3 min;94 ℃變性30 s，16S變性1 min，CO Ⅰ和16S分別48 ℃退火1 min，72 ℃延伸 1 min，共運行35個循環(huán);最后1個循環(huán)結(jié)束后72 ℃再延伸 5 min。PCR產(chǎn)物用含有溴化乙錠的1%瓊脂糖凝膠電泳檢測，凝膠成像系統(tǒng)觀察、照相。用UNIQ-10柱式PCR產(chǎn)物純化試劑盒[生工生物工程（上海）股份有限公司生產(chǎn)]純化后，進行雙向測序。

1.2.4 數(shù)據(jù)處理 本試驗選取2個巖扇貝個體進行測序，所獲原始序列首先與其對應峰圖對比檢查，以確認序列質(zhì)量;每個個體測得正反向序列并用BioEdit軟件進行拼接，結(jié)合人工校正;在NCBI數(shù)據(jù)庫中進行Blast比對分析，確定序列準確性。下載扇貝科下物種序列[16]，應用MEGA 7軟件檢查平均一致度估測聯(lián)配結(jié)果的可信度;計算16S的種內(nèi)遺傳距離;基于CO Ⅰ和16S rDNA，應用鄰接樹（neighbor joining，簡稱NJ）法分別重建系統(tǒng)發(fā)生樹;運用MEGA 7軟件Codonusage功能計算同義密碼子出現(xiàn)的次數(shù)及使用度[17]，16S基因密碼子偏好性分析，同義密碼子使用度（RSCU）計算公式為

RSCUij=xijx=xij1ni∑nij=1xij。

式中：xij表示編碼第i個氨基酸的第j個密碼子的出現(xiàn)次數(shù);ni表示編碼第i個氨基酸的同義密碼子數(shù)量（其值為1～6）。當密碼子的RSCU值等于1時，表示該密碼子沒有偏好性;當密碼子的RSCU值小于1時，表示該密碼子使用較少;當密碼子的RSCU值大于1時，表示該密碼子使用較多[18]。

2 結(jié)果與分析

2.1 巖扇貝CO Ⅰ和16S的DNA序列組成

本研究對巖扇貝的CO Ⅰ和16S基因進行PCR擴增，得到清晰的基因片段擴增產(chǎn)物，純化后直接進行雙向序列測定，經(jīng)過比對和校正后提交GenBank數(shù)據(jù)庫。測序結(jié)果表明，CO Ⅰ 基因片斷大小為660 bp，序列的A、T、G、C平均含量分別為17.2%、40.3%、27.2%、15.3%，其中A+T的含量（57%）明顯高于 G+C的含量（42%）;16S基因片段大小為613.0～634.0 bp（不包含引物），序列A、T（U）、G、C平均含量分別為26.5%、297%、24.5%、19.3%（表2），其中A+T（U）的含量（56.2%）高于G+C的含量（43.8%）。CO Ⅰ和16S堿基均出現(xiàn)偏倚性，符合線粒體堿基組成的特點。將研究中測定及下載的所有16S序列進行整理，并裁剪為同一長度529 bp進行分析，其保守位點153個，變異位點369個，簡約信息位點315個，單獨位點54個。

2.2 系統(tǒng)發(fā)育樹的物種鑒定

本試驗選取珍珠貝科馬氏珠母貝、牡蠣科長牡蠣、鶯蛤科白斑珍珠蛤作為外群，在CO Ⅰ和16S進化樹中C. gigantea單獨聚為1支，馬氏珠母貝（Pinctada martensi），白斑珍珠蛤（Pinctada maculata）聚為1支。CO Ⅰ進化樹的p-distance<08，16S序列的Jukes-Cantor分析值為0.246，這2個值分別說明CO Ⅰ和16S所選序列適合構(gòu)建NJ樹。CO Ⅰ的NJ樹表明，C. gigantea進化依次接近于櫛孔扇貝（Azumapecten farreri），冰島扇貝（Chlamys islandica），北美扇貝（Patinopecten caurinus），蝦夷扇貝（Mizuhopecten yessoensis），其中后三者聚為1支（圖1）。16S的NJ樹表明，巖扇貝5個序列形成具有較高支持度的單分支且依次接近于北美扇貝（Patinopecten caurinus）、蝦夷扇貝（Mizuhopecten yessoensis）、Chlamys islandica（圖2）。扇貝科的5個屬，其中類櫛孔扇貝屬的華貴櫛孔扇貝（Mimachlamys nobilis）單獨聚為1支;海灣扇貝屬的海灣扇貝（Argopecten irradians）、紫扇貝（A. purpuratus）、A. ventricosus聚為1支;櫛孔扇貝屬的C. farreri、A. farreri和S. livida為1支;掌扇貝屬的平瀨掌扇貝（Volachlamys hirasei）、新加坡掌扇貝（V. singaporina）聚為1支;盤扇貝屬的M. yessoensis、P. caurinus聚為1支，聚類結(jié)果基本與傳統(tǒng)的分類結(jié)果一致。

2.3 巖扇貝中16S基因的密碼子偏好性

16S在巖扇貝中擴增偏向使用以U或G結(jié)尾的密碼子，總使用度分別為12.76、8.02。其中GUU、UUG、AUU、ACG、UAU、CAA、GAG、UUU、UCU、GAU、GGU、GAG、AGG和GGA等25個密碼子的RSCU值均大于1，其中GUU的RSCU值大于2，為16S在巖扇貝中偏好性使用次數(shù)最多的（表3）。

3 討論

孔曉瑜等做了櫛孔扇貝和海灣扇貝線粒體16S rRNA基因片段的比較研究，A+T的含量分別為54.73%、55.17%[18]。胡麗萍等在紫扇貝和海灣扇貝中擴增16S，獲得基因片段長度均為542 bp，紫扇貝和海灣扇貝的堿基組成 A+T所占比例分別為54.2%、55.1%[19]。朱立靜等對四角蛤蜊CO Ⅰ基因的研究表明，四角蛤蜊的G+C含量明顯低于A+T含量[20];程漢良等對簾蛤目中6種貝類進行CO Ⅰ序列的擴增結(jié)果也表明，A+T含量明顯高于G+C含量[21]。這與本研究中CO Ⅰ在巖扇貝中的擴增結(jié)果相似。一般認為，

在線粒體基因組中，16S rRNA基因進化速率低，比較保守，而CO Ⅰ基因為變異性較大的區(qū)域，在不少無脊椎動物中檢測到了較大的變異[22-23]。從本試驗結(jié)果可以看出，在剪裁為 529 bp 的16S rRNA片段中，其保守位點153個，變異位點369個，變異位點比例達到69.8%，即巖扇貝種內(nèi)變異率較高，原因可能是16S存在較高的堿基轉(zhuǎn)換和顛換，這部分有待后續(xù)進行進一步討論。

Thompson等認為成對聯(lián)配的氨基酸序列的平均一致性百分比太低時，多序列聯(lián)配的結(jié)果準確性會下降，一致性百分比等于1減去p-distance[24]。 當p-distance﹤0.8，聯(lián)配結(jié)果是可接受的;當p-distance≥0.8則是不可靠的。本試驗CO Ⅰ的p-distance為0.746，對應的一致性為25.4%，介于20%和30%之間的臨界區(qū)域，氨基酸殘基的聯(lián)配正確率在80%左右，另一項研究表明，氨基酸聯(lián)配的精確性大于50%時，它對系統(tǒng)進化樹準確性的影響就微乎其微了，即該COⅠ系統(tǒng)進化樹是可信且可用的[25]。同理，16S序列的Jukes-Cantor分析值為0.246，說明16S的NJ樹可信且可用。

Saavedr等分析表明，C. gigantea進化上接近C. islandica，其次是巖扇貝屬的其他成員如Crassadoma multistriata、Mimachlamys varia[12]。在本試驗16S的NJ樹中，C. gigantea更接近P. caurinus和M. yessoensis，其次是C. islandica，在CO Ⅰ的NJ樹中更接近于A. farreri，2個NJ樹不一致的原因可能是2個基因片段剪裁序列的長度不同導致聯(lián)配一致度不同。Feng 等研究表明，在CO Ⅰ進化樹中芬香海扇蛤（D. plica）、大海扇蛤（P. maximus）、A. pleuronectes、A.irradians聚為1支，在16S進化樹中Pecten maximus、Amusium pleuronectes、Argopecten irradians和其他扇貝聚為1支而Decatopecten plica單獨聚為1支[16]。

參考文獻：

[1]曹善茂，汪 健，王 謙，等. 巖扇貝人工育苗的初步研究[J]. 大連海洋大學學報，2017，32（1）：1-6.

[2]曹善茂，王 昊，陳 煒，等. 巖扇貝閉殼肌營養(yǎng)成分的分析及與中國 3 種扇貝的比較[J]. 大連海洋大學學報，2016，31（5）：544-550.

[3]曹善茂，梁偉鋒，汪 健，等. 巖扇貝幼貝濾食率的基礎(chǔ)研究[J].?大連海洋大學學報，2016，31（6）：612-617.

[4]Phleger C F，Holtz R B，Grimes P W，et al.Chemical and sensory analysis of the purple-hinge rock scallo Hinnites multirugosus Gale[J]. Journal of Food Science，1978，43（6）：1793-1796.

[5]Maxwell-Miller G，Josephson R V，Spindler A A，et al.Chilled（5 ℃）and frozen（-18 ℃） storage stability of the purple-hinge rock scallop，Hinnites multirugosus Gale[J]. Journal of Food Science，1982，47（5）：1654-1661.

[6]Monical James B. Comparative studies on growth of the purple-hinge rock scallop Hinnites multirugosus （Gale） in southern California[J]. Journal of the World Aquaculture Society，1979，10（1/2/3/4）：648.

[7]Laurén D J.Oogenesis and protandry in the purple-hinge rock scallop，Hinnites giganteus，in upper Puget Sound，Washington，U.S.A.[J]. Canadian Journal of Zoology，1982，60（10）：2333-2336.

[8]Cary S C，Leighton D L，Phleger C F.Food and feeding strategies in culture of larval and early juvenile purple -hinge rock scallops，Hinnites multirugosus （Gale）[J]. Journal of the World Aquaculture Society，1981，12（1）：156-169.

[9]Culver C S，Richards J B，Page H M.Plasticity of attachment in the purple-hinge rock scallop，Crassadoma gigantea：implications for commercial culture[J]. Aquaculture，2006，254（1/2/3/4）：361-369.

[10]Bernard F R. Crassadoma gen. nov. for “Hinnites”giganteus （Gray，1825） from the North-eastern Pacific Ocean （Bivalvia：Pectinidae）[J]. Venus 1986，45：70-74

[11]Waller T R. The evolution of “Chlamys” （Mollusca：Bivalvia：Pectinidae） in the tropical western Atlantic and eastern Pacific[J]. Am Malacol Bull，1993，10：195-249.

[12]Saavedra C，Pea J B. Phylogenetics of American scallops （Bivalvia：Pectinidae） based on partial 16S and 12S ribosomal RNA gene sequences[J]. Marine Biology，2006，150：111-119.

[13]Hebert P D N，Ratnasingham S，de Waard J R . Barcoding animal life：cytochrome c oxidase subunit 1 divergences among closely related species[J]. Proceedings of the Royal Society B：Biological Sciences，2003，270（S1）：96-99.

[14]Ball S L，Hebert P D N，Burian S K，et al. Biological identifications of mayflies （Ephemeroptera） using DNA barcodes[J]. Journal of the North American Benthological Society，2005，24（3）：508-524.

[15]程希婷，王愛民，顧志峰，等. DNA條形碼研究進展[J]. 基因組學與應用生物學，2011，30（6）：748-758.

[16]Feng Y W，Li Q，Kong L F，et al. DNA barcoding and phylogenetic analysis of Pectinidae （Mollusca：Bivalvia） based on mitochondrial CO Ⅰ and 16S rRNA genes[J]. Molecular Biology Reports，2011，38（1）：291-299.

[17]宋喬喬，柴志欣，鐘金城. 禽流感病毒基因的密碼子偏好性及聚類分析[J]. 生物技術(shù)，2014，24（2）：48-53.

[18]孔曉瑜，劉亞軍，喻子牛，等. 櫛孔扇貝和海灣扇貝線粒體16S rRNA基因片段的比較研究[C]// 中國動物學會、中國海洋湖沼學會貝類學分會代表大會暨學術(shù)討論會，1999.

[19]胡麗萍，姜黎明，黃曉婷，等. 基于16S rDNA基因序列探討引進物種紫扇（Argopecten purpuratus）在海灣扇貝屬（Argopecten）中的分類地位[J]. 海洋與湖沼，2016，47（6）：1149-1157.

[20]朱立靜，陳淑吟，許曉風，等. 四角蛤蜊江蘇群體線粒體CO Ⅰ基因片段序列研究[J]. 江蘇農(nóng)業(yè)科學，2010，38（4）：33-35.

[21]程漢良，夏德全，吳婷婷，等. 6種簾蛤科貝類及4個地理種群文蛤線粒體CO Ⅰ基因片段序列分析[J]. 海洋學報，2007，29（5）：109-116.

[22]Howland D E，Hewitt G M. Phylogeny of the colepters based on mitochondrial cytochrome oxidase Ⅰ sequence data[J]. Insect Mol Bio，1995，4：203-205.

[23]Spicer G S . Phylogenetic utility of the mitochondrial cytochrome oxidase gene：molecular evolution of the Drosophila buzzatiispecies complex[J]. Journal of Molecular Evolution，1996，41（6）：749-759.

[24]Thompson J D，Plewniak F，Poch O . A comprehensive comparison of multiple sequence alignment programs[J]. Nucleic Acids Research，1999，27（13）：2682-2690.

[25]Ogdenw T H，Rosenberg M S. Multiple sequence alignment accuracy and phylogenetic inference[J]. Systematic Biology，2006，55（2）：314-328.劉冬梅，婁喜艷，吳 狄，等. 基于棉花轉(zhuǎn)錄組測序的SSR分子標記的開發(fā)[J]. 江蘇農(nóng)業(yè)科學，2019，47（7）：32-35.