朱丹麗，王曉清*，曾志南，秦溱，熊鋼，曾丹，嚴(yán)璐琪

（1.湖南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)技術(shù)學(xué)院，湖南 長(zhǎng)沙 410128；2.水產(chǎn)高效健康生產(chǎn)湖南省協(xié)同創(chuàng)新中心，湖南 常德415000；3.福建省水產(chǎn)研究所，福建 廈門 361000；4.湖南生物機(jī)電職業(yè)技術(shù)學(xué)院，湖南 長(zhǎng)沙 410127）

?

3個(gè)群體方斑東風(fēng)螺線粒體COI基因的遺傳多樣性分析

朱丹麗1，2，王曉清1，2*，曾志南3*，秦溱1，2，熊鋼1，4，曾丹1，2，嚴(yán)璐琪1，2

（1.湖南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)技術(shù)學(xué)院，湖南 長(zhǎng)沙 410128；2.水產(chǎn)高效健康生產(chǎn)湖南省協(xié)同創(chuàng)新中心，湖南 常德415000；3.福建省水產(chǎn)研究所，福建 廈門 361000；4.湖南生物機(jī)電職業(yè)技術(shù)學(xué)院，湖南 長(zhǎng)沙 410127）

采用 PCR產(chǎn)物直接測(cè)序法對(duì)來(lái)自福建東山、廣東湛江和廣西北海的 3個(gè)群體方斑東風(fēng)螺(Babylonia areolata）的COI基因進(jìn)行比較分析。結(jié)果表明：COI基因片段長(zhǎng)度為680 bp，堿基組成顯示A+T比例較高，約為62.42%；種群間多態(tài)位點(diǎn)數(shù)為32，單倍型多樣性指數(shù)為0.707±0.092，核苷酸多態(tài)性指數(shù)為0.001 85±0.000 39，平均核苷酸差異數(shù)為1.252。以泥東風(fēng)螺（Babylonia lutosa）和織紋螺（Nassarius）為外群，構(gòu)建了NJ系統(tǒng)發(fā)育樹(shù)，其拓?fù)浣Y(jié)構(gòu)顯示，3個(gè)群體的方斑東風(fēng)螺首先聚類到一起，然后與泥東風(fēng)螺聚為一支，與織紋螺形成不同的分支。

方斑東風(fēng)螺；COI基因；單倍型；遺傳多樣性

投稿網(wǎng)址：http：//xb.ijournal.cn

方斑東風(fēng)螺（Babylonia areolata）簡(jiǎn)稱東風(fēng)螺，俗稱花螺，隸屬于腹足綱，新腹足目（Neogastropoda），蛾螺科（Buccinidae），東風(fēng)螺屬（Babylonia）［1］。在中國(guó)主要分布于福建、廣東、廣西和海南等省的沿海地區(qū)［2］，是中國(guó)境內(nèi)分布的幾種東風(fēng)螺中個(gè)體較大、生長(zhǎng)較快、抗病能力強(qiáng)的一種優(yōu)質(zhì)貝類。近年來(lái)，由于市場(chǎng)需求量增加，過(guò)度捕撈造成方斑東風(fēng)螺的自然資源銳減［3］。方斑東風(fēng)螺現(xiàn)已發(fā)展成為新興的海水養(yǎng)殖對(duì)象，在中國(guó)的海南、廣東和福建等省，方斑東風(fēng)螺養(yǎng)殖行業(yè)呈上升趨勢(shì)。

線粒體DNA（mitochondrial DNA，mtDNA）是承載線粒體遺傳密碼的物質(zhì)，主要呈雙鏈環(huán)狀。mtDNA表現(xiàn)為母系遺傳，進(jìn)化速度快且不易發(fā)生重組現(xiàn)象。動(dòng)物體內(nèi)的線粒體基因（mtDNA）是唯一的核外遺傳信息載體，嚴(yán)格遵守母系遺傳，已被廣泛應(yīng)用于多種動(dòng)物類群的系統(tǒng)進(jìn)化研究。細(xì)胞色素氧化酶I亞基（cytochrome oxidasesubunit I，COI）是線粒體氧化呼吸鏈的重要成員。由于COI基因變異較大［4］，其序列已被廣泛應(yīng)用于貝類的種質(zhì)鑒定［5-6］、種群遺傳結(jié)構(gòu)分析［7］和分子系統(tǒng)發(fā)育分析［8］等領(lǐng)域。

筆者對(duì)福建東山、廣東湛江和廣西北海3個(gè)群體方斑東風(fēng)螺的線粒體COI基因進(jìn)行克隆和分析，并對(duì)其遺傳變異進(jìn)行比較研究，分析其序列組成特性、系統(tǒng)關(guān)系以及序列多態(tài)性，旨在為進(jìn)一步研究其遺傳分析、種質(zhì)鑒定和系統(tǒng)進(jìn)化等提供參考依據(jù)。

1　材料與方法

1.1材料

1.1.1試驗(yàn)動(dòng)物

試驗(yàn)用方斑東風(fēng)螺采自福建東山（DSYF）、廣東湛江（ZJYF）和廣西北海（BHYF）3個(gè)地區(qū)，采集數(shù)量分別為 27、25、22個(gè)，且均為養(yǎng)殖群體。測(cè)定采集樣品的形態(tài)指標(biāo)，然后進(jìn)行編號(hào)，放入-20 ℃冰箱保存。

1.1.2試劑

總DNA提取采用索萊寶柱式動(dòng)物DNA OUT試劑盒，PCR擴(kuò)增用Mix為全式金提供的2*Easy Taq PCR Supermix（+dye），切膠回收為全式金提供的EasyPure Quick Gel Extraction Kit瓊脂糖凝膠回收試劑盒。

1.1.3引物

COI基因引物序列為L(zhǎng)1490（5'-GGTCAACAA ATCATAAAGATATTGG-3'）和H2198（5'-TAAAC TTCAGGGTGACCAAAAAATCA-3'）［9］，均由上海鉑尚生物技術(shù)有限公司合成。

1.2方法

1.2.1基因組DNA的提取

取方斑東風(fēng)螺腹足約100 mg，按照索萊寶柱式動(dòng)物 DNA OUT 試劑盒提示的操作步驟提取總DNA?？侱NA 用1% 瓊脂糖凝膠電泳檢測(cè)，-20 ℃保存，備用。

1.2.2目的片段的PCR擴(kuò)增和測(cè)序

PCR反應(yīng)體系總體積25 μL，其中DNA模板1 μL，正反引物各0.5 μL，Mix 12.5 μL，加水至25 μL。

PCR擴(kuò)增反應(yīng)條件：94 ℃預(yù)變性5 min；94 ℃變性45 s，52 ℃退火1 min，72 ℃延伸1 min，35個(gè)循環(huán)；72 ℃延伸5 min。

PCR產(chǎn)物經(jīng)1% 瓊脂糖電泳、切膠和回收純化后，送上海鉑尚生物技術(shù)有限公司進(jìn)行雙向測(cè)序。

1.2.3序列分析

測(cè)序獲得的片段先用BLAST在線比對(duì)，確定片段的可靠性，然后由 Clustral［10］進(jìn)行多重比對(duì)分析，并用DNAman軟件分析序列長(zhǎng)度和計(jì)算堿基組成，用DnaSp 4.0［11］軟件分析多態(tài)數(shù)位點(diǎn)（singleton polymorphic sites）、單倍型多樣性（haplotype diversity，Hd） 、核苷酸多樣性（nucleotide diversity，Pi） 等遺傳多樣性參數(shù)，最后以泥東風(fēng)螺（Babylonia lutosa）和織紋螺（Nassarius）為外群，用MEGA6.0［12］軟件構(gòu)建NJ系統(tǒng)發(fā)育樹(shù)，并進(jìn)行同源性分析和遺傳距離分析。

2　結(jié)果與分析

2.13個(gè)群體方斑東風(fēng)螺COI基因片段的堿基組成

通過(guò)PCR擴(kuò)增，經(jīng)過(guò)雙向直接測(cè)序，測(cè)得COI基因的序列長(zhǎng)度約為710 bp。通過(guò)在線BLAST比對(duì)，證實(shí)所得序列為方斑東風(fēng)螺COI基因片段。通過(guò)Clustal X多重序列比對(duì)分析后，除去引物及部分雜合序列，最終獲得大小約為680 bp的COI基因片段。

由表1可知，3個(gè)群體方斑東風(fēng)螺COI基因序列中的堿基組成均顯示A+T所占比例較高，存在著明顯的A/T堿基偏向，并且C的含量最低。

表1　3個(gè)群體方斑東風(fēng)螺COI基因片段的堿基組成Table 1 The base composition of COI gene fragments from three different groups of Babylonia areolata

2.2序列多態(tài)性分析結(jié)果

遺傳多樣性檢測(cè)結(jié)果如表2所示。3個(gè)群體間存在的多態(tài)位點(diǎn)數(shù)為 7～18，單倍型多態(tài)性指數(shù)為（0.459±0.119）～（0.840±0.078），核苷酸多態(tài)性指數(shù)為（0.000 76±0.000 23）～（0.002 52±0.000 49），平均核苷酸差異數(shù)為 0.519～1.710。遺傳多樣性分析結(jié)果顯示：3個(gè)群體方斑東風(fēng)螺的核苷酸多樣性指數(shù)以東山方斑東風(fēng)螺的最低，以北海方斑東風(fēng)螺的最高。北海養(yǎng)殖群體的遺傳多樣性最高，其次是湛江養(yǎng)殖群體，東山養(yǎng)殖群體的遺傳多樣性最低。

表2　3個(gè)群體方斑東風(fēng)螺COI基因片段的遺傳多樣性參數(shù)Table 2 Genetic diversity parameters of COI gene fragments among three different groups of Babylonia areolata

由表3可知，在3個(gè)群體方斑東風(fēng)螺的COI基因序列中檢測(cè)出了28個(gè)單倍型，不計(jì)插入和缺失，共發(fā)現(xiàn)了32個(gè)變異位點(diǎn)，其中簡(jiǎn)約信息位點(diǎn)5個(gè)，轉(zhuǎn)換位點(diǎn)數(shù)30個(gè)，顛換位點(diǎn)數(shù)2個(gè)，平均轉(zhuǎn)顛換比為15。轉(zhuǎn)換方式以A-G（12個(gè)）為主，C-T（8個(gè)）其次，然后是T-C（6個(gè)），G-A（4個(gè)）最少；顛換方式只有G-C（2個(gè)）。由表4可知，具有單倍型Hap1的個(gè)體有40個(gè)，為優(yōu)勢(shì)單倍型；具有單倍型Hap26的個(gè)體有4個(gè)；具有單倍型Hap10和Hap16的個(gè)體均有3個(gè)；其余單倍型均為1個(gè)個(gè)體。從方斑東風(fēng)螺COI基因序列的組成來(lái)看，群體間的單倍型交流較少。東山群體的樣本數(shù)雖然最多，但是其單倍型數(shù)最少；北海群體樣本數(shù)雖然最少，但其單倍型數(shù)最多。

表3　3個(gè)群體方斑東風(fēng)螺COI基因的單倍型及變異位點(diǎn)位置Table 3 Haplotypes and variable loci of COI gene in three different groups of Babylonia areolata

表3　(續(xù))

表4　3個(gè)群體方斑東風(fēng)螺COI基因的共享單倍型個(gè)體數(shù)Table 4 Number of shared haplotypes of COI gene in three different groups of Babylonia areolata

2.3遺傳距離與聚類分析結(jié)果

由表5可知，東山群體內(nèi)遺傳距離為0.001，湛江群體內(nèi)遺傳距離為 0.002，北海群體內(nèi)遺傳距離為0.003，3個(gè)群體間的遺傳距離為0.002～0.005。3個(gè)群體方斑東風(fēng)螺與泥東風(fēng)螺的遺傳距離為0.127～0.129，與織紋螺的遺傳距離為0.181～0.183。

表5　3個(gè)群體方斑東風(fēng)螺及其外群之間COI基因的遺傳距離Table 5 Distance matrix among three different groups of Babylonia areolata and outgroups on COI gene fragments

拓?fù)浣Y(jié)構(gòu)顯示3個(gè)養(yǎng)殖群體的方斑東風(fēng)螺首先聚類到一起，然后與泥東風(fēng)螺聚為一支，與織紋螺形成不同的分支。3個(gè)群體之間無(wú)明顯分化現(xiàn)象。本結(jié)果與遺傳多樣性指數(shù)分析結(jié)果相一致。貝葉斯推斷樹(shù)相應(yīng)節(jié)點(diǎn)的后驗(yàn)概率值為0.001～0.099。

圖1　3個(gè)群體方斑東風(fēng)螺COI基因的NJ系統(tǒng)發(fā)育樹(shù)Fig.1 NJ phylogenetic trees based on COI gene fragments from three different groups of Babylonia areolata and outgroups

3　結(jié)論與討論

對(duì)于任何一個(gè)物種而言，遺傳多樣性越豐富，說(shuō)明該物種對(duì)環(huán)境變化的適應(yīng)能力和擴(kuò)展種群分布范圍的能力越強(qiáng)［13］。3個(gè)養(yǎng)殖群體方斑東風(fēng)螺的COI基因序列堿基組成均顯示較高的A+T比例，這與蘇天鳳等［5］對(duì)來(lái)自粵東和粵西方斑東風(fēng)螺和來(lái)自臺(tái)灣東風(fēng)螺的線粒體16S rRNA和COI基因片段的序列分析結(jié)果基本一致，且與其他貝類COI基因的分析結(jié)果也一致［13-17］。Folmer等［9］認(rèn)為，A+T含量高的線粒體基因在進(jìn)化中更具優(yōu)勢(shì)。本研究中，3個(gè)群體方斑東風(fēng)螺的A+T占比為62.42 %，說(shuō)明該物種具有較高的遺傳變異能力和進(jìn)化優(yōu)勢(shì)，但其遺傳多樣性指數(shù)分析結(jié)果顯示其具有較低的遺傳多樣性。導(dǎo)致這一現(xiàn)象的原因：一方面是由于長(zhǎng)期人工繁殖導(dǎo)致親螺的親緣關(guān)系越來(lái)越接近，減少了地理隔離造成的遺傳分化；另一方面是因?yàn)榉桨邧|風(fēng)螺自身的奠基者效應(yīng)或近期發(fā)生的瓶頸效應(yīng)。此外，人類活動(dòng)對(duì)方斑東風(fēng)螺生存環(huán)境的影響也不容忽視，如對(duì)方斑東風(fēng)螺棲息環(huán)境的破壞和近年來(lái)的過(guò)度捕撈，使野生種群日趨縮?。?8］，導(dǎo)致方斑東風(fēng)螺遺傳多樣性水平較低。

3個(gè)群體方斑東風(fēng)螺的核苷酸多樣性指數(shù)為0.001 85±0.000 39。當(dāng)核苷酸多樣性指數(shù)為0.001～0.004 7時(shí)，表明該物種種群的遺傳多樣性較低［19］。從遺傳距離來(lái)看，3個(gè)群體間遺傳距離的差異不明顯。聚類分析結(jié)果顯示，東山群體與湛江群體首先聚類到一起，然后再與北海群體聚到一起，且三者之間無(wú)明顯分化。從地理位置看，福建東山臨東海，而廣東湛江和廣西北海地處南海，盡管這與遺傳多樣性指數(shù)和聚類分析結(jié)果不相符，但從整體而言，東海與南海以臺(tái)灣島南端和閩粵兩省交界處的連線為界，兩海其實(shí)南北相連，都屬北太平洋西部的陸緣海。隨著方斑東風(fēng)螺苗種繁育技術(shù)和成螺養(yǎng)殖技術(shù)的日趨成熟，長(zhǎng)期人工養(yǎng)殖導(dǎo)致方斑東風(fēng)螺群體內(nèi)的遺傳分化較少，因此，體現(xiàn)出較小的遺傳差異。本研究中的3個(gè)群體均為養(yǎng)殖群體，其遺傳差異較小也是符合常理的。

綜上所述，開(kāi)展方斑東風(fēng)螺種群遺傳多樣性的研究和資源保護(hù)刻不容緩。首先，要加強(qiáng)漁政執(zhí)法，規(guī)定起捕規(guī)格，控制捕撈強(qiáng)度，實(shí)行嚴(yán)格的禁漁期制度，以逐步恢復(fù)方斑東風(fēng)螺的野生資源；其次，應(yīng)開(kāi)展方斑東風(fēng)螺種群遺傳多樣性研究，以便為其種群保護(hù)和遺傳多樣性監(jiān)測(cè)提供依據(jù)；最后，建議每年多進(jìn)行一些野外放流，在苗種養(yǎng)殖場(chǎng)進(jìn)行人工繁育時(shí)盡量采用遺傳距離較遠(yuǎn)的親本，避免近交繁育帶來(lái)的遺傳衰退，以利于保護(hù)其遺傳多樣性。

［1］ 王如才．中國(guó)水生貝類原色圖鑒［M］．杭州：浙江科學(xué)技術(shù)出版社，1988：92．

［2］ 張漢華，吳進(jìn)鋒，陳利雄，等．東風(fēng)螺人工育苗、養(yǎng)殖及產(chǎn)業(yè)化發(fā)展前景［J］．南方水產(chǎn)科學(xué)，2004（11）：2-5．

［3］ 陳利雄，吳進(jìn)鋒．東風(fēng)螺的增養(yǎng)殖技術(shù)及產(chǎn)業(yè)化前景［J］．齊魯漁業(yè)，2004，21（10）：9-11．

［4］ Hoeh WR，Stewart DT，Sutherland BW，et al. Cytochrome c oxidase sequence comparisons suggest an unusually high rate of mitochondrial DNA evolution in Mytilus （Mollusca：Bivalvia）［J］．Mol Biol Evol，1996，13（2）：418-421．DOI：10.1093/oxfordjournals.molbev.a025600．

［5］ 蘇天鳳，黃建華，吳進(jìn)鋒，等．2種東風(fēng)螺線粒體基因序列的多態(tài)性研究［J］．中國(guó)水產(chǎn)科學(xué)，2007，14（3）：369-376．DOI：10.3321/j.issn：1005-8737.2007.03.004．

［6］ 陳麗梅，孔曉瑜，喻子牛，等．3種蟶類線粒體16S rRNA 和 COI基因片段的序列比較及其系統(tǒng)學(xué)初步研究［J］．海洋科學(xué)，2005，29（8）：27-32．DOI：10.3969/j.issn. 1000-3096.2005.08.008．

［7］ King TL，Eackles MS，Gjetvaj B，et al．Intraspecific phylogeography of Lasmigona subviridis （Bivalvia：Unionidae）： conservation implications of range discontinuity［J］．Mol Ecol，1999，8（12 Suppl 1）：65-78．DOI：10.1046/j.1365-294x.1999.00784.x．

［8］ An HS，Jee YJ，Min KS，et al．Phylogenetic analysis of six species of Pacific abalone （Haliotidae） based on DNA sequences of 16s rRNA and cytochrome c oxidase subunit I mitochondrial genes［J］．Mar Biotechnol，2005，7（4）：373-380．DOI：10．1007/s10126-004-4405-2．

［9］ Folmer O，Black M，Hoeh W，et al．DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates［J］．Mol Marine Biol Biotechnol，1994，3（5）：294-299．

［10］ Thompson JD，Gibson TJ，Plewniak F，et al．The CLUSTAL_X Windows interface：flexible strategies for multiple sequence alignment aided by quality analysis tools［J］． Nucleic Acids Res， 1997， 25（24）：4876-4882．DOI：10．1093/nar/25．24．4876．

［11］ Librado P，Rozas J．DnaSP v5：a software for comprehensive analysis of DNA polymorphism data［J］．Bioinformatics，2009，25（11）：1451-1452．DOI：10.1093/bioinformatics/btp187.［PubMed］

［12］ Kumar S，Tamura K，Jakobsen IB，et al．MEGA2：molecular evolutionary genetics analysis software［J］. Bioinformatics，2001，17（12）：1244-1245．

［13］ 陳靈芝．中國(guó)的生物多樣性現(xiàn)狀及其保護(hù)對(duì)策［M］．北京：科學(xué)出版社，1993：99-113．

［14］ 蘇天鳳，江世貴，朱彩艷，等．廣西欽洲灣養(yǎng)殖牡蠣線粒體16SrRNA基因片段序列變異分析［J］．中國(guó)水產(chǎn)科學(xué)，2005，12（1）：1-4．

［15］ 王琳楠，閆喜武，秦艷杰，等．中國(guó)簾蛤目16種經(jīng)濟(jì)貝類DNA條形碼及分子系統(tǒng)發(fā)育的研究［J］．大連海洋大學(xué)學(xué)報(bào)，2013，28（5）：431-437．DOI：10.3969/j.issn. 2095-1388.2013.05.004．

［16］ 劉亞軍，喻子牛，姜艷艷，等．櫛孔扇貝 16S rRNA基因片段序列的多態(tài)性研究［J］．海洋與湖沼，2002，33（5）：477-483．DOI：10．3321/j．issn：0029-814X ．2002 ．05．004．

［17］ Zheng Xiaodong，Wang Rucai，Xiao Shu，et al．Genetic diversity in populations of sepiella maindroni using 16S rRNA gene sequence analysis［J］ ．High Technology Letters，2003，9（1）：1-4．

［18］ 董津茂，蔡如星．浙江動(dòng)物志（軟體動(dòng)物）［M］．杭州：浙江科學(xué)技術(shù)出版社，1991：17-115．

［19］ Lan H，Shi L．The origin and genetic differentiation of native breeds of pigs in southwest China：an approach from mitochondrial DNA polymorphism［J］．Biochem Genet，1993，31（1-2）：51-60．DOI：10.1007/bf02399819．

責(zé)任編輯：王賽群

英文編輯：王 庫(kù)

Genetic polymorphisms of three different groups of Babylonia areolata based on COI gene

Zhu Danli1，2， Wang Xiaoqing1，2*， Zeng Zhinan3*， Qin Qin1，2， Xiong Gang1，4， Zeng Dan1，2， Yan Luqi1，2
（1.College of Animal Science and Technology， Hunan Agricultural University， Changsha 410128， China； 2.Collaborative Innovation Center for Aquatic Efficient Health Production of Fisheries in Hunan Province， Changde， Hunan 415000，China； 3.Fisheries Research Institute of Fujian， Xiamen， Fujian 361000， China； 4.Hunan Biological Electromechanical Vocational and Technical College， Changsha， 410127， China）

We analyzed and compared the polymorphism of mitochondrial COI gene sequence derived from three different groups of Babylonia areolata in Dongshan， Zhanjiang and Beihai by means of PCR. The results showed that the nucleotide sequence of COI gene was 680 bp with a high percentage of A+T base composition （about 62.42%）； there were 32 singleton poylmorphic sites， 0.707±0.092 haplotype diversity， 0.001 85±0.000 39 nucleotide polymorphisms index， and 1.252 average number of nucleotide differences. Taking Babylonia lutosa and Nassarius as outgroups， we built NJ phylogenetic tree， its topology suggested that all the Babylonia areolata clustered together first， then mixed with Babylonia lutosa for one， and finally Nassarius formed independent.

Babylonia areolata; COI gene； phyletic evolution； genetic polymorphism

朱丹麗（1991—），女，上海人，碩士研究生，主要從事水產(chǎn)養(yǎng)殖研究，463609105@qq.com；*通信作者，王曉清，博士，教授，主要從事水產(chǎn)動(dòng)物育種研究，wangxiao8258@126.com；*通信作者，曾志南，博士，研究員，主要從事海水養(yǎng)殖和貝類遺傳育種研究，xmzzn@sina.com

Q959.212

A

1007-1032（2016）04-0429-06

2015-11-02 修回日期：2016-06-20

海洋公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)項(xiàng)目（201205021-1）；國(guó)家自然科學(xué)基金項(xiàng)目（31372530）