田云方，葉瑩瑩，吳常文，傅澤欽

（浙江海洋大學(xué)海洋科學(xué)與技術(shù)學(xué)院，國(guó)家海洋設(shè)施養(yǎng)殖工程中心，浙江舟山 316022）

基于線粒體COI基因序列的麗文蛤群體遺傳多樣性和遺傳結(jié)構(gòu)分析

田云方，葉瑩瑩，吳常文，傅澤欽

（浙江海洋大學(xué)海洋科學(xué)與技術(shù)學(xué)院，國(guó)家海洋設(shè)施養(yǎng)殖工程中心，浙江舟山 316022）

為了評(píng)估麗文蛤的遺傳多樣性，為麗文蛤的保護(hù)與開(kāi)發(fā)提供基礎(chǔ)遺傳信息，本研究利用線粒體COI基因片段序列比較分析了我國(guó)沿海麗文蛤4個(gè)地理群體（漳州、珠海、?？诤捅焙＃┑倪z傳多樣性和遺傳結(jié)構(gòu)。4個(gè)群體共54個(gè)樣本的線粒體COI基因部分序列經(jīng)處理得到長(zhǎng)度均為623 bp的核苷酸片段，檢測(cè)到19個(gè)單倍型。群體遺傳多樣性分析顯示，4個(gè)群體的單倍型多樣性為0.464～0.889，核苷酸多樣性為0.000 80～0.011 04。單倍型多樣性較高的是漳州群體和?？谌后w，北海群體單倍型多樣性相對(duì)較低。核苷酸多樣性中漳州群體最高，而北海群體同樣顯示最低，北海群體多樣性較低的原因可能是樣本量小。分子變異分析（AMOVA）結(jié)果顯示96.20%的變異來(lái)自群體內(nèi)。群體間遺傳分化指數(shù)顯示漳州與珠海和北海群體之間存在較大的遺傳分化（FST值分別為0.162 49和0.117 31），且均達(dá)到顯著水平；而漳州與?？谌后w之間遺傳分化小且分化不顯著。珠海、?？诤捅焙Ｈ后w之間遺傳分化?。‵ST值0.070 99～0.094 30）且分化不顯著。聚類(lèi)分析和單倍型網(wǎng)絡(luò)圖均顯示麗文蛤漳州群體與珠海和北海群體分化明顯，而珠海、?？诤捅焙Ｈ后w之間分化不明顯，其結(jié)果與FST值相一致。

麗文蛤；線粒體COI；遺傳多樣性；遺傳結(jié)構(gòu)

麗文蛤Meretrix lusoria屬瓣鰓綱，簾蛤科，文蛤?qū)?，生活在潮間帶和淺海沙質(zhì)海底。它內(nèi)含豐富的蛋白質(zhì)，具有高營(yíng)養(yǎng)、易養(yǎng)殖的特點(diǎn)，是我國(guó)沿海灘涂養(yǎng)殖的主要品種之一。麗文蛤外觀與文蛤極為相似，莊啟謙[1]提出文蛤?qū)僭谖覈?guó)僅有3種，即文蛤M.meretrix、麗文蛤和斧文蛤M.lamarckii。關(guān)于文蛤和斧文蛤的研究有很多，但是對(duì)于麗文蛤的研究?jī)H限于形態(tài)學(xué)[2]及外緣凝集素[3]和貝殼采集[4]的相關(guān)研究。而關(guān)于麗文蛤遺傳差異的研究?jī)H有陳愛(ài)輝等[5]基于線粒體COI基因序列的文蛤?qū)伲ㄜ涹w動(dòng)物門(mén)：簾蛤科）系統(tǒng)發(fā)育關(guān)系。針對(duì)于麗文蛤線粒體DNA標(biāo)記研究麗文蛤的遺傳結(jié)構(gòu)的文章，未有報(bào)道。

由于近年來(lái)過(guò)度捕撈和比較嚴(yán)重的海洋污染，文蛤?qū)俑魑锓N的群體多樣性面臨極大挑戰(zhàn)，由于麗文蛤在外形上與文蛤存在著諸多相似，一直被忽視，游離在研究之外，尚未從分子水平加以研究分析。潘寶平等[6]首先利用分子生物學(xué)手段對(duì)麗文蛤和文蛤、斧文蛤的親緣性進(jìn)行了研究。陳愛(ài)輝等[7]基于線粒體COI基因序列的文蛤?qū)俚南到y(tǒng)發(fā)育研究。TORRI，et al[8]中國(guó)，韓國(guó)和印度的麗文蛤的研究，標(biāo)明是關(guān)于麗文蛤的形態(tài)及遺傳特征的研究，但是文章內(nèi)容著重研究形態(tài)學(xué)及麗文蛤的分類(lèi)，并沒(méi)有采用分子手段對(duì)麗文蛤群體遺傳多樣性進(jìn)行研究。而國(guó)內(nèi)尚未有采用分子手段對(duì)麗文蛤的群體多樣性進(jìn)行研究的報(bào)道。

細(xì)胞色素C氧化酶亞基I（COI）是線粒體DNA中一個(gè)重要基因，基因變異性較大、進(jìn)化速率比較快，在采用分子手段研究貝類(lèi)[9-10]和物種群體的遺傳多樣性及結(jié)構(gòu)分析[11-13]等方向應(yīng)用廣泛。本文采用線粒體COI基因進(jìn)行標(biāo)記，對(duì)麗文蛤四個(gè)不同沿海地理種群的分子遺傳學(xué)研究分析，為我國(guó)麗文蛤種質(zhì)資源的研究和麗文蛤及文蛤?qū)傥锓N多樣性研究提供理論依據(jù)。

1 材料與方法

1.1 樣品采集和DNA提取

本研究所用麗文蛤分別來(lái)源于4個(gè)地理群體，為廣西北海（BH）、海南海口（HK）、廣東珠海（ZH）、福建漳州（ZZ）如圖1。每個(gè)地點(diǎn)隨機(jī)采樣活體解剖取閉殼肌(表1)，保存于無(wú)水乙醇中，用于DNA提取。采用改良鹽析法[14]提取總DNA，利用1%的瓊脂糖電泳進(jìn)行檢測(cè)，無(wú)菌超純水稀釋至50 ng/μL備用。

1.2 線粒體COI序列擴(kuò)增

麗文蛤線粒體COI擴(kuò)增引物來(lái)源于通用引物L(fēng)CO1490和HCO2198[15]，COI-F 5'-GGTCAACAAATCATAAAGATATTGG-3'COI-R 5'-TAAACTTCAGGGTGACCAAAAAATCA-3'PCR反應(yīng)體系為 25.0 μL，內(nèi)含 Taq 酶 mix 12.5 μL，正、反向引物各 1.0 μL（10 μmol/L），模板 DNA 1.0 μL，ddH2O 9.5 μL。PCR 反應(yīng)程序?yàn)?94 ℃預(yù)變性為 2 min，94 ℃變性時(shí)間為30 s，52℃退火時(shí)間為30 s，72℃延伸時(shí)間為1 min，35個(gè)循環(huán)反應(yīng)后，72 ℃延伸時(shí)間為10 min。PCR產(chǎn)物經(jīng)1%～2%瓊脂糖凝膠電泳檢測(cè)后，PCR產(chǎn)物最后送杭州擎科（杭州擎科梓熙生物技術(shù)有限公司）進(jìn)行雙向測(cè)序。

圖1 本研究麗文蛤地理群體采集點(diǎn)地圖Fig.1 Collection locations of M.lusoria geographical populations in this study

表1 麗文蛤4個(gè)群體單倍型多樣性及核苷酸多樣性Tab.1 Haplotype diversity and nucleotide diversity of four populations of M.lusoria

1.3 數(shù)據(jù)統(tǒng)計(jì)與分析

實(shí)驗(yàn)測(cè)序全部采用雙向測(cè)序，為了確認(rèn)樣品分類(lèi)的準(zhǔn)確性，序列結(jié)果均在NCBI（https：//www.ncbi.nlm.nih.gov/）上進(jìn)行Blast對(duì)比，統(tǒng)計(jì)好單倍型并上傳NCBI（https：//www.ncbi.nlm.nih.gov/）數(shù)據(jù)庫(kù)，并獲得登錄號(hào)（表2）。利用Mega版本5.0[16]對(duì)獲得的COI序列進(jìn)行整理編輯，結(jié)合人工對(duì)編輯結(jié)果進(jìn)行矯正。利用DNASP[17]統(tǒng)計(jì)單倍型，計(jì)算每個(gè)實(shí)驗(yàn)群體的單倍型多樣性（h）和核苷酸多樣性（π）。利用Arlequin 3.5[18]軟件計(jì)算群體間及群體內(nèi)的遺傳分化指數(shù)（FST），并進(jìn)行P檢驗(yàn)，通過(guò)10 000次重復(fù)置換檢驗(yàn)（Permutation test）分析其顯著性。同時(shí)分析群體內(nèi)和群體間的分子方差（Analysis of Molecular Variance，AMOVA），以10 000次重復(fù)抽樣檢驗(yàn)其顯著性，并進(jìn)行Fu`s Fs檢驗(yàn)。利用Network 5.0.0.1以核苷酸進(jìn)化模型計(jì)算各單倍型之間的遺傳距離并構(gòu)建單倍型網(wǎng)絡(luò)圖，利用 Mega 5.0[15]以Kimura 2-parameter模型建立4個(gè)群體的Unweighted Pair Group Method with Arithmatic Mean（UPGMA）單倍型系統(tǒng)進(jìn)化樹(shù)，抽樣次數(shù)設(shè)置為1 000次。

2 結(jié)果

2.1 群體內(nèi)遺傳變異

對(duì)4個(gè)麗文蛤地理群體的DNA進(jìn)行PCR擴(kuò)增，獲得所需的COI基因片段，經(jīng)BLAST比對(duì)分析，確認(rèn)所得片段為麗文蛤COI序列，經(jīng)Mega同源排序，去除部分端部序列，獲得54條623 bp的COI基因片段。序列組成顯示平均 A、T、G、C 堿基含量分別為 A：11.43%、T：41.54%、C：23.92%、G：23.11%。其中 A+T：52.97%、C+G：47.03%。A+T含量大于G+C的含量，具有明顯的堿基組成偏倚性。共檢測(cè)出單倍型19個(gè)（表2）（登錄號(hào)MF285209-MF285227），其中H1被漳州、珠海、?？诤捅焙Ｈ后w共享，H3被漳州和?？谌后w共享，Hap9被珠海和?？谌后w共享。H2、H5、H6、H7、H8為漳州群體特有單倍型，H10-19為?？谌后w特有單倍型。

麗文蛤4個(gè)地理群體的遺傳學(xué)參數(shù)見(jiàn)表2，平均單倍型及核苷酸多樣性處于中等水平（h=0.851，π=0.040 01）。單倍型多樣性最高的是漳州群體（h=0.945）和?？谌后w（h=0.917），北海群體單倍型多樣性相對(duì)較低（h=0.750）。核苷酸多樣性中漳州群體最高（π=0.179 7），而北海群體核苷酸多樣性顯示最低（π=0.014 7）。北海群體和珠海群體的單倍型多樣性和核苷酸多樣性偏低，這可能是由于樣本量比較小，但是我們不排除樣品本身的原因造成珠海和北海本身麗文蛤的遺傳多樣性低于其他地區(qū)。

2.2 群體遺傳結(jié)構(gòu)研究

用MEGA5.0構(gòu)建4個(gè)群體單倍型系統(tǒng)發(fā)育樹(shù)如圖2，系統(tǒng)發(fā)育樹(shù)可以看出，單倍型9以下沒(méi)有珠海和?？谌后w出現(xiàn)。漳州和海口群體大部分群體都聚集在下方，且沒(méi)有明顯分化。用Network5.0.1.1構(gòu)建單倍型網(wǎng)絡(luò)圖如圖3，結(jié)果顯示單倍型網(wǎng)絡(luò)圖右側(cè)沒(méi)有珠海和北海群體出現(xiàn)，這表明漳州與珠海和北海群體發(fā)生了顯著分化。這與系統(tǒng)發(fā)育樹(shù)和FST值結(jié)果一致。

麗文蛤4個(gè)群體COI序列的分子變異分析（AMOVA）結(jié)果見(jiàn)表4，群體間遺傳變異占5.38%，群體內(nèi)遺傳變異占94.60%，群體內(nèi)的遺傳變異大于群體間，群體間分化系數(shù)P值小于0.05，達(dá)到顯著水平，表明群體間發(fā)生顯著的遺傳分化（P＜0.05）。

表2 單倍型序列比對(duì)Tab.2 Sequence comparison of haplotypes

表3 麗文蛤4個(gè)群體的COI序列選擇中性檢驗(yàn)Tab.3 Selective neutrality test analysis of 4 populations of M.lusoria based on COI sequences

圖2 麗文蛤4個(gè)群體的單倍型UPGMA系統(tǒng)進(jìn)化樹(shù)Fig.2 Haplotype UPGMA tree of 4 populations of M.lusiria

圖3 麗文蛤4個(gè)群體的單倍型網(wǎng)絡(luò)圖Fig.3 Haplotype network of 4 populations of M.lusiria

表4 麗文蛤COI序列遺傳差異的分子方差分析（AMOVA）Tab.4 Analysis of molecular variance(AMOVA)of COI gene of M.lusiria

麗文蛤群體間的FST顯著性檢測(cè)見(jiàn)表5，兩兩之間的差異數(shù)值介于0.018 16～0.162 49之間，其中珠海群體和漳州群體（FST=0.162 49）（P＜0.05）、北海群體和漳州群體（FST=0.117 31）（P＜0.05），說(shuō)明珠海群體和漳州群體，北海群體和漳州群體發(fā)生了顯著分化。從地理位置來(lái)看，珠海所在地理位置位于珠江三角洲，海洋生物類(lèi)易于堆積。而北海所在位置與漳州之間隔有陸地，這可能就是阻擋基因交流的重要原因。漳州與?？谥g沒(méi)有阻隔，斧文蛤群體以及幼體均可以隨海流和潮流四處移動(dòng)，基因間交流較大，分化不明顯。中性檢驗(yàn)結(jié)果顯示，平均Fu`s Fs值為-0.501 04（P＞0.05）（表3），F(xiàn)u`s Fs值未達(dá)到顯著水平，表明麗文蛤未偏離中性選擇，麗文蛤群體沒(méi)有經(jīng)歷過(guò)群體擴(kuò)張。

表5 基于線粒體COI序列的麗文蛤群體遺傳分化Tab.5 Estimates of population genetic differentiation based on COI genefor the M.lusiria

3 討論

線粒體COI基因具有高度的保守性，一般多可以用于物種的多樣性比較。本研究檢測(cè)到麗文蛤線粒體COI序列的單倍型多樣性為0.790，核苷酸多樣性為0.007 27，與我國(guó)其他海水貝類(lèi)[19-21]相比，遺傳多樣性不高。這可能與我們實(shí)驗(yàn)樣本量不足引起的，但也有可能我們實(shí)驗(yàn)所采集的樣本中有效群體小、群體之間親本性別比例不平衡等[22]造成，從而導(dǎo)致麗文蛤群體得COI基因序列遺傳多樣性偏低。從單倍型網(wǎng)絡(luò)圖來(lái)看，各群體單倍型2～10個(gè)不等。本實(shí)驗(yàn)研究樣品來(lái)源于我國(guó)主要的文蛤養(yǎng)殖海域，由于麗文蛤與文蛤在形態(tài)學(xué)上難以鑒定，文蛤各群體的養(yǎng)殖規(guī)模大，我們也可以說(shuō)是麗文蛤養(yǎng)殖規(guī)模比較大。養(yǎng)殖群體的群體遺傳多樣性偏低，放殖到野生環(huán)境中與野生群體想混合生殖，導(dǎo)致野生群體的遺傳多樣性整體偏低。

群體遺傳中的遺傳分化系數(shù)是衡量物種分化的重要指標(biāo)，本文4個(gè)群體麗文蛤群體間的FST為0.053 8（P＜0.05），表明群體間的相似性不大，群體間發(fā)生了顯著分化，而群體內(nèi)的遺傳變異達(dá)到0.946 2（P＜0.05），群體內(nèi)分化顯著。其中群體間珠海群體和漳州群體（FST=0.162 49）（P＜0.05）、北海群體和漳州群體（FST=0.117 31）（P＜0.05），漳州與珠海和北海群體間分化顯著。某些雙殼類(lèi)必要時(shí)可憑借貝殼的急劇開(kāi)合和外套膜觸手的作用在海中進(jìn)行蝶式游泳。大部分水生貝類(lèi)營(yíng)底棲生活，或在水底匍匐、爬行，或在底質(zhì)中挖穴隱居，或附著在其他外物上生活。所以貝類(lèi)會(huì)隨著海流和潮流進(jìn)行移動(dòng)[23]，本研究中的四個(gè)地理位置中珠海位于珠江三角洲，海洋生物易于堆積，而不易于遷移。北海與漳州之間被海峽阻隔，基因交流不大。相對(duì)來(lái)說(shuō)珠海、海口與北海之間相距較近，基因交流的機(jī)會(huì)較大，所以群體間沒(méi)有明顯的遺傳分化。漳州與?？谥g雖然相距甚遠(yuǎn)但是地理位置間沒(méi)有阻隔，隨著海流與潮流，基因交流的可能性較大，因此群體間分化不明顯。實(shí)驗(yàn)所采用的樣品為野生樣品，流動(dòng)性較大，這可能與野生群體的異地養(yǎng)殖和增殖放流有關(guān)。李宏俊等[24]也認(rèn)為貝類(lèi)的養(yǎng)殖可能對(duì)當(dāng)?shù)匾吧后w的遺傳結(jié)構(gòu)造成影響。單倍型系統(tǒng)發(fā)育樹(shù)顯示，漳州群體和珠海和北海群體發(fā)生了顯著分化，系統(tǒng)發(fā)育樹(shù)H9以下沒(méi)有珠海和北海群體出現(xiàn)，單倍型網(wǎng)絡(luò)圖顯示右側(cè)沒(méi)有珠海和北海群體出現(xiàn)，表明珠海群體和北海群體與漳州群體發(fā)生明顯的分化，這與地理環(huán)境和海流有關(guān)[23]，也可能物種本身經(jīng)過(guò)長(zhǎng)時(shí)間的自身遺傳進(jìn)化所引起的。

我們衡量物種進(jìn)化潛力的依據(jù)是分析物種的遺傳多樣性和遺傳結(jié)構(gòu)，通常情況下，遺傳多樣性的高低代表物種適應(yīng)環(huán)境的能力。本研究表明麗文蛤群體在經(jīng)歷長(zhǎng)時(shí)間的進(jìn)化，已經(jīng)產(chǎn)生了明顯的分化，這表明麗文蛤自身的適應(yīng)能力在逐步提升。但是由于近年來(lái)海洋污染的加劇，海洋生物面臨極大的挑戰(zhàn)，很多物種面臨滅絕的危險(xiǎn)。文蛤?qū)賰H有三種貝類(lèi)，文蛤及斧文蛤的研究明顯對(duì)于麗文蛤，麗文蛤的分子研究也一直游離于研究之外，由于麗文蛤與文蛤在外形上極為相似，不易分辨，所以麗文蛤的養(yǎng)殖及野生群體的總體數(shù)量并不是很清楚，這很不利于麗文蛤種資資源的保護(hù)與開(kāi)發(fā)。我們的研究結(jié)果將為麗文蛤的深入研究及其開(kāi)發(fā)保護(hù)提供幫助。

4 小結(jié)

本文采用4個(gè)野生麗文蛤?yàn)闃颖緦?duì)麗文蛤的遺傳差異進(jìn)行線粒體COI分子研究，結(jié)果顯示群體間發(fā)生了明顯的分化，群體內(nèi)分化較為明顯，這可能與文蛤群體的野生異地養(yǎng)殖和增值放流有關(guān)，但不排除麗文蛤本身的基因特點(diǎn)。麗文蛤作為文蛤?qū)倨渲幸环N，研究麗文蛤的遺傳差異對(duì)文蛤?qū)俚恼w物種多樣性具有重要意義。目前關(guān)于麗文蛤的研究相對(duì)較少，這對(duì)其本身的發(fā)展及保護(hù)不利。麗文蛤和文蛤相類(lèi)似，比較容易養(yǎng)殖，而且味道鮮美，所以今后麗文蛤的養(yǎng)殖和野生馴養(yǎng)以及野生群體和養(yǎng)殖群體的差異性比較等都是值得研究的課題。

[1]莊啟謙.中國(guó)動(dòng)物志軟體動(dòng)物門(mén)雙殼綱簾蛤科[M].北京：科學(xué)出版社,2001：229-236.

[2]陳愛(ài)華,吳楊平,姚國(guó)興,等.麗文蛤與文蛤4個(gè)地理群體的形態(tài)差異分析[J].海洋漁業(yè),2010,30(2)：125-131

[3]朱越雄,曹廣力.麗文蛤血清中外源凝集素的凝集性能[J].海洋科學(xué),2000(10)：55-56.

[4]小池裕子,李慧冬,王 青.麗文蛤在推斷繩文時(shí)期貝類(lèi)采集季節(jié)中的作用[J].南方文物,2008(3)：130-133.

[5]陳愛(ài)輝,李朝霞,封功能.基于線粒體 COI基因序列的文蛤?qū)?軟體動(dòng)物門(mén)：簾蛤科)系統(tǒng)發(fā)育關(guān)系[J].動(dòng)物學(xué)研究,2009,30(3)：233-239.

[6]潘寶平,吳 琪,張素萍,等.文蛤?qū)?Meretrix)16SrRNA基因及ITS1序列的系統(tǒng)學(xué)分析[J].海洋與湖沼,2006,37(4)：342-347.

[7]陳愛(ài)輝,李朝霞,封功能.基于線粒體 COI基因序列的文蛤?qū)?軟體動(dòng)物門(mén)：簾蛤科)系統(tǒng)發(fā)育關(guān)系[J].動(dòng)物學(xué)研究,2009,30(3)：233-239.

[8]TORII H S,SATO S C,et al.The comparison of shell morphology and geneticrelationship between Meretrix lusoria and M.petechialis in Japan and Korea.Plankton and Benthos Research[J].Plankton&Benthos Research,2010(5)：231-241.

[9]程漢良,夏德全,吳婷婷,等.6種簾蛤科貝類(lèi)及4個(gè)地理種群文蛤線粒體COI基因片段序列分析[J].海洋學(xué)報(bào),2007,29(5)：109-116.

[10]程漢良,彭永興,董志國(guó),等.基于線粒體細(xì)胞色素c氧化酶亞基I基因序列的簾蛤科貝類(lèi)分子系統(tǒng)發(fā)育研究[J].生態(tài)學(xué)報(bào),2013,33(9)：2 744-2 753.

[11]牛東紅,李家樂(lè),沈和定,等.縊蟶六群體線粒體DNA-COI基因序列變異及群體遺傳結(jié)構(gòu)分析[J].海洋學(xué)報(bào),2008,30(3)：109-116.

[12]沈玉幫,張俊彬,馮冰冰,等.基于線粒體COI序列分析對(duì)紫貽貝群體遺傳多樣性的研究分析[J].海洋通報(bào),2011,30(4)：435-440.

[13]ZHAN A B,PEREPELIZIN P V,GHABOOLI S,et al.Scale-dependent post-establishment spread and genetic diversity in an invading mollusc in South America[J].Diversity and Distributions,2012,18(10)：1 042-1 055.

[14]ALJANABI S M,MARTINEZ I.Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques[J].Nucleic acids Res,1997(25)：4 692-4 693.

[15]FOLMER O,BLACK M,HOEH W,et al.DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates[J].Molecular Marine Biology Biotechnology,1994,3(5)：294-299.

[16]TAMURA K,DUDLEY J,NEI M,et al.MEGA 4：molecular evolutionary genetics analysis(MEGA)software version 4.0[J].

Molecular Biology and Evolution,2007,24(8)：1 596-1 599.

[17]ROZAS J,SáNCHEA-DELBARRIO J C,MESSEGUER X,et al.DNA polymorphism analyses by the coalescentand other methods[J].Bioinformatics,2003,19(18)：2 496-2 497.

[18]EXCOFFIER L,LISCHER H E L.Arlequin suite ver 3.5：a new series of programs to perform population genetics analyses under Linux and Windows[J].Molecular Ecology Resources,2010,10(3)：564-567.

[19]王 超,陳愛(ài)華,曹 奕,等.6個(gè)不同海域文蛤地理群體的親緣關(guān)系分析[J].海洋漁業(yè),2016,38(3)：262-272.

[20]鄭文娟,朱世華,沈錫權(quán),等基于線粒體COI基因序列探討泥蚶的遺傳分化[J].動(dòng)物學(xué)研究,2009,30(1)：17-23.

[21]MAO Y L,GAO T X,YANAGIMOTO T,et al.Molecular phylogeography of Ruditapes philippinarum in the Northwestern Pacific Ocean based on COI gene[J].Journal of Experimental Marine Biology and Ecology,2011,407(2)：171-181.

[22]AVISE J C,ARNOLD J,Ball R M,et al.Intraspecificphylogeography：the mitochondrial DNA bridge between population genetics and systematics[J].Annual Review of Ecology and Systematics,1987,18：489-522.

[23]馬 超.黑潮對(duì)東中國(guó)海主要流系的影響[D].青島：中國(guó)海洋大學(xué),2009.

[24]李宏俊,張晶晶,袁秀堂,等.利用線粒體COI和微衛(wèi)星標(biāo)記分析文蛤7個(gè)地理群體的遺傳變異[J].生態(tài)學(xué)報(bào),2016,36(2)：499-507.

Genetic Diversity and Genetic Structure of Meretrix lusoria Population based on Mitochondrial COI Gene Sequences

TIAN Yun-fang,YE Ying-ying,WU Chang-wen,et al
(School of Ocean Science and Technology of Zhejiang Ocean University,National Engineering Research Center For Marine Aquaculture,Zhoushan 316022,China)

In order to evaluate the genetic diversity of Meretrix lusoria,make a contribution to protect and develope the M.lusoria.In this study,the genetic diversity of four geographical populations(Zhangzhou,Zhuhai,Haihou and Beihai)of M.lusoriawas calculated based on the mitochondrial cytochrome c oxidase subunit I(COI)gene sequences.A 623 bp nucleotide fragment was obtained in the 54 individuals in all 4 populations and 19 haplotypes were detected in them.Population genetic diversity analysis showed that the haplotype diversity of each population was 0.464-0.889.Zhangzhou shows the highest haplotype diversity in the four populations,and the haplotype diversity of Beihai is relatively low.This may be caused by a small sample size.Molecular variation analysis(AMOVA)results showed 96.20%variation from within the population.The genetic differentiation indexes between Zhangzhou/Zhuhai and Zhangzhou/Beihai populations showed a large genetic differentiation (FSTvalues is 0.162 49 and 0.117 31,respectively),and reached a significant level.The genetic differentiation between Zhangzhou and Haikou was not significant.The genetic differentiation between Zhuhai/Beihai and Haikou/Beihai populations were not significant.Cluster analysis and haplotype network map showed that the differentiation between Zhangzhou and Zhangzhou populations was significant,while the differentiation among Zhuhai,Haikou and Beihai populations werenot obvious,and the results were consistent with the FSTvalues.

Meretrix lusoria;MT COI;genetic diversity;genetic structure

Q959.215+.4

A

2096-4730(2017)04-0320-06

2017-05-10

浙江海洋大學(xué)“水產(chǎn)”省一流學(xué)科開(kāi)放課題（20160002）

田云方（1990-），女，河南上蔡人，碩士研究生，研究方向：海洋生物學(xué).E-mail：xiaotian081215@126.com

葉瑩瑩，研究方向：分子遺傳育種.E-mail：Yeyingying5559@163.com