徐其建 黃鎮(zhèn) 林崗

摘?要：初步探究凹線仙女蜆的系統(tǒng)發(fā)育地位，為凹線仙女蜆的分子遺傳學(xué)研究提供前期數(shù)據(jù)。采用PCR方法擴增得到凹線仙女蜆Cyrenobatissa subsulcata線粒體細胞色素b（Cytb）基因片段序列，序列長為624 bp，并在NCBI上使用NT數(shù)據(jù)庫進行Blast比對分析，進而構(gòu)建系統(tǒng)進化樹，探討凹線仙女蜆的分類地位。結(jié)果顯示NT數(shù)據(jù)庫中未見凹線仙女蜆的Cytb基因信息。研究擴增獲得的序列與蜆屬其他種的Cytb基因片段有88%的相似性，與其他簾蛤目貝類的相似度77.02%～78.53%，表明擴增出的片段是凹線仙女蜆所特有的Cytb基因片段。基于Kimura′s 2-Parameter模型計算得出，3個樣本個體的遺傳距離僅為0～0.0016，遠小于蜆屬內(nèi)種間遺傳距離（0.002～0.114），與蜆屬的遺傳距離為0.127～0.144、與簾蛤目簾蛤科種類遺傳距離為0.282～0.335。進化樹顯示，凹線仙女蜆樣品獨立聚為一支，與蜆屬聚為一大支，支持度為99，表明仙女蜆屬與蜆屬的親緣關(guān)系較近，與簾蛤科的遺傳距離較遠。

關(guān)鍵詞：凹線仙女蜆;Cytb基因片段;進化樹;遺傳距離

中圖分類號：S917.4?文獻標志碼：A?文章編號：0253-2301（2020）10-0001-09

DOI： 10.13651/j.cnki.fjnykj.2020.10.001

Abstract： The phylogenetic status of Cyrenobatissa subsulcata was preliminarily explored in this study to provide preliminary data for the molecular genetics research of Cyrenobatissa subsulcata. The PCR method was used to amplify the mitochondrial cytochrome b （Cytb） gene fragment sequence of Cyrenobatissa subsulcata. The sequence length was 624 bp. The Blast comparison analysis was performed on the NCBI by using the NT database to construct a phylogenetic tree and discuss the taxonomic status of Cyrenobatissa subsulcata. The results showed that there was no Cytb gene information of Cyrenobatissa subsulcata in the database. The sequence amplified in this study had 88% similarity with the Cytb gene fragments of other species of the Corbicula， and had 77.02%-78.53% similarity with other clam shellfish， indicating that the fragments amplified by PCR were from the Cytb gene fragment of Cyrenobatissa subsulcata. Based on Kimura′s 2-parameter model， the genetic distance between the three individuals of Cyrenobatissa subsulcata was only 0-0.0016， which was far less than that between the species in Corbicula （0.002-0.114） and with Corbicula was 0.127-0.144， and the genetic distance between the three individuals of Corbicula and the species of Veneridae was 0.282-0.335. The phylogenetic tree showed that the samples of Cyrenobatissa subsulcata were clustered into one branch independently and clustered with Corbicula into a large branch with a support degree of 99， which indicated that Cyrena was closely related to Corbicula， while had a long genetic distance from Veneridae.

Key words： Cyrenobatissa subsulcata;Cytb gene fragment;Phylogenetic tree;Genetic distance

凹線仙女蜆Cyrenobatissa subsulcata隸屬于軟體動物門Mollusca、雙殼綱Bivalvia、異齒亞綱Heterodonta、簾蛤目Veneroida、蜆科Corbiculidae、仙女蜆屬Cyrenobatissa。凹線仙女蜆是中國的地方性物種，分布具有局限性，主要棲息于底質(zhì)為沙或沙泥的咸淡水交匯處，主要分布于廣東?。ㄉ穷^韓江、吳川鑒江）、廣西?。ǚ莱牵?、海南?。ㄅ鸥郏┖团_灣省（淡水、蘇澳、安平）[1]。蜆科種類在沿海河口半咸水區(qū)域也有較多分布，硬殼蜆屬的紅樹蜆Polymesoda erosa分布于中國、菲律賓等地區(qū)，我國主要分布于廣東省、廣西省、海南省、臺灣地區(qū)等[2];花蜆屬的花蜆Cyrenodonax formosana分布于河口砂質(zhì)區(qū)，產(chǎn)于我國臺灣地區(qū)（淡水）、福建?。ㄔ葡觯┖蛷V西?。ǚ莱牵3];蜆屬的河蜆廣泛分布于江河、湖泊、溝渠及咸淡水交匯處等水域，在福建省的閩江、烏龍江地區(qū)有百年的養(yǎng)殖歷史，是河口底棲生物的重要組成部分和重要的經(jīng)濟貝類[4-5]。

凹線仙女蜆的研究主要在分類學(xué)和生態(tài)學(xué)[6]，對河蜆Corbicula fluminea等蜆科物種研究較多，包括生態(tài)特性、遺傳差異分析、生長與繁殖等方面[7-9]。凹線仙女蜆是蜆科中個體最大的物種，具有較高的經(jīng)濟價值[6]。蜆的軟體部含有較低的脂肪、較高的蛋白質(zhì)，蜆肉酶解產(chǎn)物具有解酒護肝的功效，可食用也可作為天然調(diào)味品，具有較高的營養(yǎng)價值[10]。由于過度捕撈等一系列對自然環(huán)境的不合理利用問題，造成河流生態(tài)系統(tǒng)結(jié)構(gòu)破壞、生物多樣性驟減和生態(tài)系統(tǒng)生產(chǎn)力下降，有些地區(qū)的蜆類已處于瀕危狀態(tài)，制約了蜆類產(chǎn)業(yè)的發(fā)展[11-13]。因此，對凹線仙女蜆種質(zhì)資源的保護、利用研究已顯得十分緊迫。物種的分子鑒定是種質(zhì)資源研究的重要基礎(chǔ)，細胞色素b基因（Cytb）作為一種線粒體DNA分子標記，其結(jié)構(gòu)和功能是線粒體DNA的13個蛋白質(zhì)編碼基因中了解最清楚的基因，已廣泛用于魁蚶Scapharca broughtonii、鰱Hypophthalmichthys molitrix、松江鱸Trachidermus fasciatus和河蜆Corbicula fluminea等水生物種的分子鑒定、種間變異程度和親緣關(guān)系、不同地理群體遺傳多樣性和遺傳結(jié)構(gòu)等研究[14-19]，但目前尚未見凹線仙女蜆線粒體基因研究的相關(guān)報道。

本研究以分布于福州市長樂區(qū)凹線仙女蜆群體為研究對象，對其Cytb基因片段進行擴增分析，探討Cytb基因片段作為凹線仙女蜆種類鑒定的分子標記的可行性，可為凹線仙女蜆的遺傳多樣性、種質(zhì)資源研究提供基礎(chǔ)資料。

1?材料與方法

1.1?試驗材料

1.1.1?凹線仙女蜆樣本信息?凹線仙女蜆Cyrenobatissa subsulcata個體全部取自福州市長樂區(qū)文武砂鎮(zhèn)規(guī)劃中的濱海新區(qū)濕地公園（圖1）。

1.1.2?試驗儀器?PCR儀、凝膠成像分析系統(tǒng)、Nanodrop分光光度計、體視熒光顯微鏡、電泳儀、解剖鏡。

1.1.3?試驗試劑?蛋白酶K、Quick-DNA Universal Kit 試劑盒（Zymo Research生物科技公司）、瓊脂糖（Sigma-Aldrich）、引物（福州擎科生物技術(shù)有限公司合成）、 DNA marker（Thermo Scientific）、 PCR mixture（北京全式金生物技術(shù)有限公司）。

1.2?試驗方法

1.2.1?形態(tài)學(xué)性狀測定?從采集來的凹線仙女蜆中隨機挑選30個，將貝殼表面的水分吸干后，用游標卡尺測量仙女蜆的殼長、殼高、殼寬，精確到0.01 mm。用電子天平測量活體的軟體重、殼重、總重精確到0.01 g。凹線仙女蜆的外部形態(tài)見圖2。

1.2.2?凹線仙女蜆基因組DNA的提取?采用Quick-DNA Universal Kit試劑盒提取基因片段，在所選的30個凹線仙女蜆中隨機挑選3個個體，分別提取，每個個體取0.1 g的肌肉組織放入EP管，剪碎后加入95 μL無菌水、95 μL Solid Tissue Buffer（blue）和10 μL蛋白酶K。將EP管放入55℃的水浴鍋中水浴1.5 h。當(dāng)組織溶解后，12000 r·min-1離心10 min，取上層液150 μL，移入新的EP管，加300 μL的Genomic Binding Buffer溶液混合均勻?；旌弦恨D(zhuǎn)入Zymo-spin IIC-xl小管中，12000 r·min-1離心1 min，棄掉下管的液體。再加 400 μL的DNA pre-wash Buffer到小管，12000 r·min-1離心1 min。然后加700 μL的g-DNA wash Buffer液體到小管，12000 r·min-1速度下離心1 min，棄掉下管的液體。最后加入37 μL雙蒸水。12000 r·min-1離心1 min后得到DNA溶液，之后進行電泳，并用Nanodrop測定DNA完整度和純度。

1.2.3?凹線仙女蜆Cytb基因片段的克隆和序列測定?Cytb基因引物設(shè)計參考Yamada等[20]，序列擴增通用引物為CBF6（5′GCAAGCTTTTGATTCTGTTGTTCACATTG3′）和CBR6（5′GTAGGAATCCTACGCAAAATGAGAATAAGCG3′ ）。據(jù)梯度PCR確定最適溫度57℃，在此溫度下進行PCR產(chǎn)物的大體系擴增。PCR 100 μL體系如下：模板DNA 2 μL、上下游引物1 μL、2×PCR mixture 50 μL和無菌水46 μL。分別以3只凹線仙女蜆基因組DNA作模板，放入PCR儀擴增。擴增程序為：預(yù)變性94℃ 7 min，變性94℃ 1 min，復(fù)性57℃ 45 s，72℃延伸2 min，擴增35個循環(huán)結(jié)束后，72℃再延伸7 min。預(yù)期擴增片段大小640 bp，將PCR產(chǎn)物進行1%的凝膠電泳，選擇亮度完好、有特異性條帶的PCR產(chǎn)物送到尚亞生物公司測序。

1.3?數(shù)據(jù)分析

運用Excel對形態(tài)學(xué)數(shù)據(jù)進行初步處理，計算變異系數(shù)。利用Mega7軟件將本研究得到的基因片段在NCBI上使用NT數(shù)據(jù)庫進行Blast比對和序列信息分析。運用Mega7軟件內(nèi)Clustal W進行多序列比對，并輔以人工校對，計算出凹線仙女蜆Cytb基因片段序列長度以及堿基組成;通過Kimura′s 2-Parameter模型計算出凹線仙女蜆與蜆屬的其他種Cytb基因片段序列的遺傳距離;通過Maximum Likelihood方法對凹線仙女蜆Cytb基因片段進行分子進化樹的構(gòu)建。

2?結(jié)果與分析

2.1?凹線仙女蜆形態(tài)學(xué)性狀分析

30個樣品各性狀表型統(tǒng)計量見表1。各性狀的變異系數(shù)差異較大，其中殼長、殼高、殼寬形態(tài)性狀的相對較小，而軟體重、殼重和總重性狀的相對較大，其中軟體重的變異系數(shù)最大，為79.35%。

2.2?凹線仙女蜆Cytb基因片段分析

2.2.1?凹線仙女蜆基因組DNA的提取?電泳結(jié)果（圖3）顯示，獲得的DNA基因片段大小均在2000 bp以上。所提取的DNA大小完整，條帶清晰，說明試劑盒能夠有效地提取凹線仙女蜆DNA。進一步對DNA的純度進行檢測，濃度均在19 ng·μL-1左右，A260/280均在1.80左右。說明提取的DNA純度高，無RNA、蛋白質(zhì)污染，可以進行下一步PCR 擴增反應(yīng)。

2.2.2?凹線仙女蜆Cytb基因片段擴增?用所設(shè)計的引物在57℃退火溫度下進行PCR擴增，PCR產(chǎn)物的電泳結(jié)果見圖4。結(jié)果表明，所擴增的3個凹線仙女蜆Cytb基因的PCR擴增產(chǎn)物，大小在640 bp左右，測定的序列峰圖完整，沒有重疊峰的出現(xiàn)，得到的PCR產(chǎn)物均為單一條帶，說明成功擴增出凹線仙女蜆Cytb基因片段。將3個PCR擴增產(chǎn)物進行測序，獲得序列。

2.2.3?凹線仙女蜆Cytb基因片段比對分析?經(jīng)查詢，在NCBI數(shù)據(jù)庫尚無凹線仙女蜆Cytb基因序列。將測序得到的凹線仙女蜆Cytb基因片段序列在NCBI上使用NT數(shù)據(jù)庫進行Blast比對及相似性檢索，得到100條比對信息，結(jié)果見表2。擴增得到凹線仙女蜆的Cytb基因片段序列與蜆屬的河蜆、日本蜆、環(huán)紋蜆的Cytb基因片段比對的相似度在86.49%～88.25%;與簾蛤目貝類的相似度77.02%～78.53%;與凡氏費氏繭蜂、雙斑蟳、三疣梭子蟹的Cytb基因片段比對的相似度在74.12%～74.55%，表明凹線仙女蜆Cytb基因片段序列區(qū)分度顯著。

2.2.4?凹線仙女蜆Cytb基因片段序列?經(jīng)測序獲得的3條Cytb基因序列中，最長為663 bp，最短為640 bp。運用Mega7軟件內(nèi)Clustal W進行多序列比對，并輔以人工校對，凹線仙女蜆Cytb基因片段的核酸序列為624 bp（圖5），該段序列共編碼208個氨基酸。從圖5可以看出，對3個樣品的PCR產(chǎn)物進行比對，序列幾乎一致，只在1個堿基上存在突變。

2.2.5?凹線仙女蜆Cytb基因片段序列特征分析?運用Mega7軟件內(nèi)Clustal W進行多序列比對，分析其中624 bp的序列，計算其堿基組成，結(jié)果（表3）表明：T（U）、C、A、G的平均堿基組成為45.3%、12.4%、21.5%、20.8%。保守位點623個（99.84%），變異位點1個（0.16%），顛換數(shù)為0。其中胸腺嘧啶的比例較高，胞嘧啶的比例偏低;堿基含量T>A>G>C，其中，A+T的含量（66.8%）大于C+G的含量（33.2%），符合線粒體DNA的特征。

2.2.6?凹線仙女蜆種群內(nèi)遺傳距離?通過Mega7軟件中的Kimura′s 2-Parameter模型計算試驗樣品的Cytb基因片段序列的遺傳距離。結(jié)果如表4所示：凹線仙女蜆種群內(nèi)的遺傳距離在0～0.00160，平均遺傳距離為0.00053，說明種群內(nèi)差異小。

2.2.7?凹線仙女蜆與參比種的屬間遺傳距離?通過Kimura′s 2-Parameter模型計算出凹線仙女蜆與蜆屬的其他種Cytb基因片段序列的遺傳距離（表5）。3個樣品個體的平均遺傳距離為0.00053，遠遠小于蜆屬內(nèi)種間遺傳距離（0.002～0.114）、與蜆屬的遺傳距離為0.127～0.144、與簾蛤目的簾蛤科種類遺傳距離為0.282～0.335，表明仙女蜆屬與蜆屬能明確區(qū)分開，也表明采集的樣品為同一物種。

2.2.8?凹線仙女蜆的系統(tǒng)進化樹?基于Maximum Likelihood 方法對Cytb基因片段序列進行對比和分析，構(gòu)建了包括凹線仙女蜆在內(nèi)的8個雙殼類物種的系統(tǒng)進化樹（圖6）。結(jié)果顯示，凹線仙女蜆樣品獨立聚為一支，蜆屬的河蜆與環(huán)紋蜆聚集再與仙臺蜆聚類，最后與日本蜆聚集為一支。仙女蜆屬與蜆屬又聚為一大支，支持度為99，表明仙女蜆屬和蜆屬的親緣關(guān)系較近，與簾蛤科的親緣關(guān)系較遠。

3?討論

環(huán)境因素是造成雙殼類動物個體差異的主要因素之一。蜆科軟體動物有因棲息地環(huán)境因素形成高度形態(tài)變異的特性，僅利用殼體形態(tài)特征劃分河蜆種類比較困難[21-27]。基于DNA分子標記技術(shù)的物種分子鑒定可以彌補形態(tài)學(xué)分類上存在的樣本量多、易形成誤差、不能分析基因?qū)哟蔚奈⒂^變化情況等缺點，有利于群體遺傳育種工作的研究，提高可靠性[28]。而線粒體DNA中Cytb基因分子標記更被大量運用在物種鑒定、系統(tǒng)發(fā)育和群體遺傳學(xué)研究[29-31]。

凹線仙女蜆的殼長、殼寬、殼高的變異系數(shù)差異較小，軟體重的變異系數(shù)最大。這與張俊杰等[32]報道紅樹蜆、王慶恒等[33]報道波紋巴非蛤的研究結(jié)果基本一致。

凹線仙女蜆Cytb基因片段A+T的含量為66.8%，明顯大于C+G的含量（33.2%），這一結(jié)果與李瑤瑤等[16]在雙殼貝類Cytb基因中觀察到的結(jié)果一致。凹線仙女蜆Cytb基因片段序列與蜆屬的河蜆、日本蜆、環(huán)紋蜆Cytb基因片段Blast比對結(jié)果，結(jié)合系統(tǒng)進化樹分析表明仙女蜆屬與蜆屬能明確區(qū)分開，也從分子水平上支持凹線仙女蜆為蜆科仙女蜆屬而異于蜆屬[34]。

保護優(yōu)質(zhì)的貝類種質(zhì)資源成為亟待解決的問題，本研究獲得的Cytb基因片段序列可用于凹線仙女蜆的分子鑒定，可為凹線仙女蜆的種質(zhì)資源研究、蜆科物種間系統(tǒng)進化關(guān)系的研究、凹線仙女蜆的生物資源保護和開發(fā)利用提供基礎(chǔ)資料。

參考文獻：

[1]楊文，蔡英亞， 鄺雪梅.中國南海經(jīng)濟貝類原色圖譜[M].北京：中國農(nóng)業(yè)出版社，2013.

[2]蔡英亞，黃翔鵠，吳洞科.紅樹蜆的生態(tài)觀察[J].熱帶海洋，1995（1）：94-98.

[3]徐鳳山，張素萍.中國海產(chǎn)雙殼類圖志[M].北京：科學(xué)出版社，2008.

[4]鄧勇輝，付輝云，陶志英，等.我國河蜆研究進展與展望[J].江西水產(chǎn)科技，2016（2）：33-37.

[5]劉月英，張文珍，王躍先，等.中國經(jīng)濟動物志——淡水軟體動物[M].北京：科學(xué)出版社，1979.

[6]ZHANG M K.Catalogue of Fresh Water Shells of Taiwan[J].The Malacological Society of ROC，1991（16）：86-96.

[7]杞桑，林美心.珠江廣州河段河蜆種群的若干生態(tài)學(xué)研究[J].生態(tài)學(xué)報，1987（2）：161-169.

[8]聞海波.我國主要淡水經(jīng)濟貝類的種質(zhì)資源研究[D].南京：南京農(nóng)業(yè)大學(xué)，2009.

[9]吳叢迪.洪澤湖河蜆繁殖特性及黑黃殼色遺傳差異分析[D].上海：上海海洋大學(xué)，2019.

[10]張帥，章超樺，高加龍，等.河蜆肉酶解產(chǎn)物解酒護肝功效[J].食品與機械，2019，35（3）：166-171.

[11]嚴維輝，潘元潮，郝忱，等.洪澤湖底棲生物調(diào)查報告[J].水利漁業(yè)，2007（3）：65-66.

[12]袁永滸，王興元，陳安來，等.洪澤湖螺蜆資源調(diào)查報告[J].水產(chǎn)養(yǎng)殖，1994（6）：15-16.

[13]畢婷婷，侯剛，張勝宇，等.基于殼長頻率分析的洪澤湖河蜆漁獲物年齡結(jié)構(gòu)、生長與死亡參數(shù)估計[J].水生生物學(xué)報，2014，38（4）：797-800.

[14]ARIF I A， KHAN H A，BAHKALI A H，et al.DNA marker technology for wildlife conservation[J]. Saudi J Biol，2011，18（3）：219-225.

[15]CHEN J，SONG Z，RONG M，et al.The association analysis between Cytb polymorphism and growth traits in three Chinese donkey breeds[J]. Livestock ence，2009，126（1-3）：306-309.

[16]李瑤瑤，劉云國，劉凌霄，等.基于線粒體基因ND1、COX1和Cytb的雙殼貝類分子系統(tǒng)進化研究[J].福建農(nóng)林大學(xué)學(xué)報（自然科學(xué)版），2018，47（2）：188-197.

[17]王劍平.洞庭湖蜆屬的系統(tǒng)分類與生殖特征研究[D].長沙：湖南農(nóng)業(yè)大學(xué)，2018.

[18]魏建軍，尤宏爭，夏蘇東，等.松江鱸魚線粒體Cytb基因分析及其在種質(zhì)鑒定中的應(yīng)用[J].河北漁業(yè)，2019（8）：21-25，61.

[19]楊習(xí)文，劉熠，薛向平，等.基于微衛(wèi)星標記的長江江蘇段鰱（Hypophthalmichthys molitrix）增殖放流資源貢獻率的評估[J].湖泊科學(xué)，2020，32（4）：1154-1164.

[20]YAMADA M，ISHIBASHI R，KAWAMURA K，et al.Interrelationships of the freshwater clams Corbicula leana Prime，1864. and C. fluminea （Müller， 1774） distributed in Japan inferred from shell type and mitochondrial DNA Cytb region[J]. Nippon Suisan Gakkaishi，2010，76（5）：926-932.

[21]HILLIS D M， PATTON J C. Morphological and Electrophoretic Evidence for Two Species of Corbicula （Bivalvia： Corbiculidae） in North America[J]. The American Midland Naturalist，1982，108（1）：74-80.

[22]MORTON B. Polymorphism in Corbicula fluminea（Bivalvia：Corbiculoidea）from HongKong[J].Malacological Review，1987，20：105-127.

[23]ALYAKRINSKAYA I O. Functional Significance and Weight Properties of the Shell in Some Mollusks[J].Biology Bulletin，2005，32（4）：397-418.

[24]BRITTON J C，MORTON B.Polymorphism in Corbicula fluminea（Bivalvia：Corbiculoidea）from North America[J].Malacological Review，1986，19：1-44.

[25]RENARD E，BACHMANN V，CARIOU M L，et al.Morphological and molecular differentiation of invasive freshwater species of the genus Corbicula（Bivalvia， Corbiculidea） suggest the presence of three taxa in French rivers[J].Molecular Ecology，2000，9：2009-2016.

[26]BAGATINI Y M，PANARARI R S.Morphological and molecular characterization of Corbicula（Mollusca， Bilvalvia）at Rosana Resevoir，Brazil[J].Acta Scientiarum Biological Sciences，2005，27：397-404

[27]SOUSA R，F(xiàn)REIRE R，RUFINO M，et al.Genetic and shell morphological variability of the invasive bivalve Corbicula fluminea（Mǜller， 1774） in two Portuguese estuaries[J].Estuarine，Coastal and Shelf Science，2007，74：166-174.

[28]劉瑩.雙齒圍沙蠶的分子系統(tǒng)地理學(xué)研究[D].青島：中國海洋大學(xué)，2014.

[29]辛一.線粒體COI、COII和CYTB基因在鮑屬物種鑒定中的適用性分析[J].海洋科學(xué)，2011，35（11）：58-62.

[30]于濤，吳彪，楊愛國，等.櫛孔扇貝、蝦夷扇貝及其雜交子代線粒體COI和Cytb基因遺傳多樣性分析[J].海洋科學(xué)，2016，40（3）：1-9.

[31]位正鵬.翡翠貽貝（Perna viridis）線粒體COI片段序列特性和香港巨牡蠣（Crassostrea hongkongensis）線粒體基因組全序列研究[D].青島：中國海洋大學(xué)， 2008.

[32]張俊杰，周浩郎，邢永澤，等.紅樹蜆表型性狀對秋季軟體部重量的影響效果分析[J].廣西科學(xué)，2012，19（4）：384-386.

[33]王慶恒，鄧岳文，杜曉東.波紋巴非蛤Paphia undulata表型性狀通徑和回歸分析[J].熱帶海洋學(xué)報，2010，29（5）：132-135.

[34]徐鳳山.中國海雙殼類軟體動物[M].北京：科學(xué)出版社，1997.

（責(zé)任編輯：柯文輝）