王琳楠,閆喜武,秦艷杰,聶鴻濤,牛泓博,張國范

(1.大連海洋大學水產(chǎn)與生命學院遼寧省貝類良種繁育工程技術(shù)研究中心,遼寧大連116023;2.中國科學院海洋研究所海洋生物技術(shù)研發(fā)中心,山東青島266071)

中國簾蛤目16種經(jīng)濟貝類DNA條形碼及分子系統(tǒng)發(fā)育的研究

王琳楠1,閆喜武1,秦艷杰1,聶鴻濤1,牛泓博1,張國范2

(1.大連海洋大學水產(chǎn)與生命學院遼寧省貝類良種繁育工程技術(shù)研究中心,遼寧大連116023;2.中國科學院海洋研究所海洋生物技術(shù)研發(fā)中心,山東青島266071)

對中國簾蛤目主要經(jīng)濟貝類的分子系統(tǒng)發(fā)育進行了研究,應(yīng)用DNA條形碼通用引物擴增了6種中國近海簾蛤目經(jīng)濟貝類共計60個個體的COI基因片段,與GenBank收錄的10種簾蛤目貝類50條同源序列進行比對。結(jié)果表明:簾蛤目貝類COI基因存在堿基插入和缺失現(xiàn)象,在16個物種中有6個物種存在103個插入和缺失位點,其中雜色蛤仔Ruditapes variegata、裂紋哥特蛤Katelysia hiantina插入和缺失位點均為30個,大竹蟶Solen grandis為27個;堿基的組成出現(xiàn)偏倚現(xiàn)象,A+T含量 (64.2%)明顯高于G+C含量(35.8%);基于K2P模型的計算,16個物種的種內(nèi)平均遺傳距離為0.010 6,種間平均遺傳距離為0.388 4,后者是前者的21.34倍;系統(tǒng)發(fā)育樹聚類分析表明,COI基因在科、屬、種水平上的鑒定及其系統(tǒng)進化關(guān)系重構(gòu)方面與傳統(tǒng)的形態(tài)學分類一致性較高。研究表明,線粒體COI基因作為簾蛤目貝類DNA條形碼在物種鑒定的適用性上提供了一定的依據(jù),同時也為形態(tài)分類系統(tǒng)提供了必要補充。

簾蛤目;COI基因;系統(tǒng)發(fā)育;分類鑒定;DNA條形碼

簾蛤目Veneroida隸屬于軟體動物門Mollusca、雙殼綱Bivalvia,是世界上分布較廣的貝類[1]。中國的簾蛤目貝類種類豐富,已有記錄的近千種,若僅僅依據(jù)殼型、放射肋、小月面等外部形態(tài)特征進行種屬水平的形態(tài)學分類難度很大[2-3],需要以基于分子生物學的分析方法進行補充。DNA條形碼技術(shù) (DNA barcoding)是利用COI基因構(gòu)建物種的鑒別體系[4-5],現(xiàn)已成為物種鑒定的一門新興技術(shù)。DNA條形碼的兩大主要目的是鑒定已知種和發(fā)現(xiàn)新種,它相對于傳統(tǒng)的生物鑒定的優(yōu)勢就是可以揭示隱存種。近年來,COI基因序列已經(jīng)廣泛應(yīng)用于貝類的種質(zhì)鑒定[6-10]、 種群遺傳結(jié)構(gòu)分析[11-13]、分子系統(tǒng)發(fā)生研究[14-21]等領(lǐng)域。 研究證實,COI基因序列適用種內(nèi)、種間和屬以上階元的遺傳多樣性分析[22-25]。本研究中以中國黃渤海、東海、南海16種重要的簾蛤目經(jīng)濟貝類為材料,由于傳統(tǒng)的形態(tài)學鑒定一般都是根據(jù)殼形、鉸合部和閉殼肌痕等形態(tài)特征來進行鑒定,而實際上殼形的多變以及發(fā)育中存在的中間形態(tài)和中間生境,使得形態(tài)學分類比較困難,為此本研究基于COI序列的DNA條形碼技術(shù)對這些貝類進行分析,從而證明DNA條形碼在簾蛤目貝類輔助物種鑒定中的適用性,旨在為全球貝類DNA條形碼數(shù)據(jù)庫提供可靠的資料。

1 材料與方法

1.1 材料

試驗用60個簾蛤目樣品均采自中國沿海地區(qū),參照 《中國海產(chǎn)雙殼類圖志》進行形態(tài)學鑒定,這些樣品隸屬于4科、5屬、6種,分別為櫻蛤科亮櫻蛤?qū)?種、竹蟶科竹蟶屬2種、簾蛤科畸心蛤?qū)?種、簾蛤科鏡蛤?qū)?種、綠螂科綠螂屬1種。

1.2 方法

1.2.1 DNA的提取 從GenBank 中篩選出已經(jīng)發(fā)表的簾蛤目COI基因序列信息50條,分別屬于簾蛤目3科、9屬、10種,共選取簾蛤目16個種, 110個個體進行研究。選取長牡蠣為外類群 (表 1)。試驗時每個種取10～12個,取其斧足、閉殼肌或者外套膜約50 mg,采用酚氯仿法抽提DNA。DNA溶于TE中,于-20℃下保存?zhèn)溆谩?/p>

表1 16種簾蛤目貝類的COI基因信息Tab.1 Information of COI genes in 16 species in Veneroida

1.2.2 PCR擴增及產(chǎn)物測序 COI片段采用無脊椎動物通用的COI引物序列LCO1490(5＇GGTCAACAAATCATAAAGATATTGG 3＇)和 HC02198 (5＇TAAACTTCAGGGTGACCAAAAAATCA 3＇)[26]。引物由英濰捷基 (上海)貿(mào)易有限公司合成。擴增產(chǎn)物經(jīng)20 g/L瓊脂糖凝膠電泳檢測后,由Invitrogen(上海)有限公司進行雙向測序。

1.3 數(shù)據(jù)處理

利用 Clustal X 1.81軟件進行序列比對排序[27],經(jīng)DnaSP 4.0[28]單倍型分析后,將單倍型序列提交GenBank數(shù)據(jù)庫。利用Mega 4.0軟件統(tǒng)計不同樣本的核苷酸組成、突變位點、遺傳距離、轉(zhuǎn)換/顛換值[18,29];利用DnaSP 4.0計算單倍型數(shù)量、單倍型多態(tài)性、核苷酸多樣性指數(shù)、多態(tài)性位點數(shù),以及插入和缺失位點數(shù);使用Mega 4.0中的Kimura雙參數(shù)法(K2P)計算種內(nèi)和種間的遺傳距離,并構(gòu)建基于K2P的鄰接樹(neighbor joining, NJ)和最大簡約樹(maximum parsimony,MP)[30]。

2 結(jié)果

2.1 簾蛤目COI序列分析

本研究中,對6種簾蛤目的COI基因進行PCR擴增,挑選較好的PCR產(chǎn)物進行測序,經(jīng)過比對和校正后提交GenBank數(shù)據(jù)庫,得到的登錄號見表1。將研究中測定及下載的所有COI序列進行整理,并裁剪為同一長度628 bp進行分析。統(tǒng)計結(jié)果表明,在研究的110個個體中共產(chǎn)生了56個單倍型,單倍型多態(tài)性指數(shù)為0.533 0～1.000 0;核苷酸多態(tài)性指數(shù)為0.001 07～0.088 04;在16個物種 (GP1～GP16)中有6個物種存在堿基插入和缺失現(xiàn)象,占到了37.50%,共存在103個插入和缺失位點 (表2)。

表2 16種簾蛤目貝類的COI基因序列特征Tab.2 The characters of COI sequences from 16 species in Veneroida

110條簾蛤目貝類COI序列的平均堿基組成: T為 41.4%,C為 14.3%,A為 22.8%,G為21.4%,其中A+T的含量 (64.2%)明顯高于G+ C的含量 (35.8%),堿基出現(xiàn)明顯的偏倚性,符合線粒體堿基組成的特點 (表3)。

表3 16種簾蛤目貝類的COI序列堿基組成的頻率Tab.3 Average nucleotide frequencies of COI sequences in 16 species in Veneroida %

本研究中的110條簾蛤目COI序列的核苷酸變異情況如表4所示。

2.2 種內(nèi)和種間的遺傳距離

基于K2P算法計算出110個個體的簾蛤目COI基因序列的種內(nèi)和種間遺傳距離如表5所示。16個物種的種內(nèi)平均遺傳距離為0.010 6,56.25%的種內(nèi)遺傳距離小于0.010 0,其中裂紋哥特蛤的種內(nèi)遺傳距離最大,為0.018 2,并未超出DNA條形碼COI基因種類遺傳距離2%的界限。

表4 16種簾蛤目貝類的COI部分序列密碼子的變異Tab.4 The variation in some codons of COI gene sequences of 16 species in Veneroida

2.3 簾蛤目的分類及系統(tǒng)發(fā)育關(guān)系

從基于NJ法和MP法構(gòu)建的簾蛤目系統(tǒng)發(fā)育樹(圖1)可見,用兩種方法構(gòu)建的系統(tǒng)進化樹得到相似的拓撲結(jié)構(gòu),但對不同分類階元而言,分子系統(tǒng)關(guān)系與傳統(tǒng)形態(tài)學分類結(jié)果不盡相同。在種水平上,NJ樹和MP樹的分類一致性均為100%;在屬水平上,分子系統(tǒng)與形態(tài)學的分類結(jié)果略有差異。

表5 16種簾蛤目貝類的種間和種內(nèi)遺傳距離Tab.5 Genetic distances in intraspecific and interspecies in 16 species in Veneroida

圖1 基于NJ法和MP法構(gòu)建的簾蛤目系統(tǒng)發(fā)育樹Fig.1 NJ and MP trees based on COI gene data from Veneroida

3 討論

3.1 簾蛤目線粒體COI基因的結(jié)構(gòu)特征

本研究結(jié)果表明,在簾蛤目16個物種中有6個物種COI基因存在堿基插入和缺失現(xiàn)象,占到了37.50%。共存在103個插入和缺失位點,其中雜色蛤仔、裂紋哥特蛤插入和缺失位點均多達30個,大竹蟶插入和缺失位點為27個。堿基的插入和缺失現(xiàn)象在軟體動物中較為常見[31-36],這可能與軟體動物進化過程中線粒體拷貝數(shù)的變化有一定關(guān)聯(lián)。簾蛤目貝類COI基因堿基組成中也存在著很普遍的偏倚現(xiàn)象。朱立靜等[37]對四角蛤蜊COI基因的研究表明,四角蛤蜊的A+T含量明顯高于G+C含量;程漢良等[38]對簾蛤目中6種貝類進行COI序列的擴增結(jié)果也表明,A+T含量明顯高于G+C含量。這與本研究結(jié)果相似。線粒體堿基組成的偏倚現(xiàn)象符合無脊椎動物線粒體基因組的特點,其他無脊椎動物如蝦、蟹、頭足類等均存在這種相似的現(xiàn)象[39]。

COI基因存在著大量序列變異,適于用其進行群體遺傳學多樣性分析。在線粒體序列變異中,變異位點的轉(zhuǎn)換較易在近緣種間頻繁發(fā)生,而顛換在較遠緣種間頻繁發(fā)生;在同種動物中,轉(zhuǎn)換往往在數(shù)量上遠超過顛換。本研究中的物種分別屬于不同種屬關(guān)系,故顛換多于轉(zhuǎn)換。

3.2 DNA條形碼與形態(tài)學方法在簾蛤目貝類物種鑒定中的異同

DNA條形碼對物種的鑒定可以彌補傳統(tǒng)形態(tài)學鑒定存在的不足,為形態(tài)學鑒定提供充足的理論依據(jù)[26]。徐鳳山[1]用形態(tài)學鑒定方法,將本研究中的16個種分成了6科,其中簾蛤科有菲律賓蛤仔、雜色蛤仔、畸心蛤、江戶布目蛤、日本鏡蛤、文蛤、青蛤、綴錦蛤、巴非蛤和裂紋哥特蛤;蛤蜊科有四角蛤蜊;櫻蛤科有虹光亮櫻蛤;綠螂科有中國綠螂;竹蟶科有大竹蟶和長竹蟶;截蟶科有縊蟶。本研究中構(gòu)建的NJ系統(tǒng)樹在科水平上與用形態(tài)學方法鑒定的一致性為83.3%,其中虹光亮櫻蛤與中國綠螂分別屬于不同的科,卻聚為一支,而在MP樹中則不存在這個問題;在屬水平上,相同屬的種類可以很好地聚為一支,一致性為100%。本研究中發(fā)現(xiàn),青蛤與文蛤聚為一支,綴錦蛤、巴非蛤與裂紋哥特蛤聚為一支,大竹蟶、長竹蟶與縊蟶聚為一支,這可以說明他們之間的親緣關(guān)系較近,這與形態(tài)學分類結(jié)果相符合[1],這樣就為以往的形態(tài)學分類增加了可靠的依據(jù)。

3.3 DNA條形碼在簾蛤目貝類物種鑒定中的適用性和局限性

本研究中所開發(fā)的全部條形碼序列逐一與生物條形碼數(shù)據(jù)庫 (BoLD)中已經(jīng)公布的生物條形碼進行比對,匹配率高達100%,并與形態(tài)學鑒定相吻合,表明DNA條形碼在簾蛤目貝類鑒定方面具有較高的適用性,完全可以通過COI基因序列的對比完成物種的準確鑒定。同時Hebert等[26]對動物界11個門13 320個物種的COI基因進行分析,發(fā)現(xiàn)其序列間的差異能夠很好地區(qū)分所有研究物種,并認為在動物界中COI基因是合適的DNA條形碼標準基因。Vences等[40]于 2005年探討了DNA條形碼對兩棲爬行類物種的鑒定能力,結(jié)果顯示選用的COI基因能夠準確辨別各物種。此外Vences等[40]選用16S rRNA基因進行條形碼分析,發(fā)現(xiàn)在兩棲爬行類中16S rRNA基因具有較強的引物適應(yīng)性,并且序列相對保守。對于國際上分類混亂的兩棲爬行類,使用DNA條形碼技術(shù),不僅能夠準確地鑒別出各物種,而且還能準確無誤地識別出其不同地理類群的變異。對北美哥斯達黎加地區(qū)1 000多種鱗翅目昆蟲生物多樣性的調(diào)查中,Hajibabaei等[41]利用DNA條形碼的COI基因序列建立了一個分類學平臺,能夠鑒定物種以及發(fā)現(xiàn)這一區(qū)域存在的模糊種,并且能準確鑒定97%的物種,區(qū)分不能被確認的物種是親緣關(guān)系極近的物種或形態(tài)定義模糊的物種,由此認為,DNA條形碼能夠很好地協(xié)助進行一個地區(qū)的生物多樣性調(diào)查。Hebert等[26]利用COI基因序列對北美260余種鳥類進行了有效區(qū)分。利用DNA條形碼也能夠?qū)︳~類以及熱帶鱗翅目昆蟲進行有效的物種鑒別。

Hebert等[26]指出,利用COI序列有效地進行物種鑒別的關(guān)鍵是種間的遺傳距離必須大于種內(nèi)的遺傳距離。Hebert等[26]在對動物界11個門13 320個物種進行研究后得出,物種種內(nèi)的遺傳距離很少有大于2%的,大部分物種種內(nèi)遺傳距離小于1%。本研究結(jié)果表明,16種簾蛤目貝類種間的平均遺傳距離為0.388 4,種內(nèi)平均遺傳距離為0.010 6,種間平均遺傳距離為種內(nèi)平均遺傳距離的21.34倍;而種內(nèi)遺傳距離介于0.002 3～0.018 2,均小于2%,與Herbert等提出的觀點相符。

COI適用于物種鑒定的原因在于COI基因既保證了足夠多變異的同時又很容易被通用引物進行擴增,而且發(fā)生插入和缺失的也是很小的一部分序列,甚至不會發(fā)生插入和缺失現(xiàn)象。本研究中發(fā)現(xiàn)簾蛤目的貝類有一些堿基插入和缺失現(xiàn)象,孔曉瑜等[42]的研究發(fā)現(xiàn),魁蚶也存在堿基的插入和缺失現(xiàn)象。而且這種插入和缺失的現(xiàn)象普遍存在于其他各種動物體中。

雖然基于COI基因的DNA條形碼可以較好地實現(xiàn)物種的鑒定,但DNA條形碼仍存在一定的缺陷,不可以完全取代形態(tài)學方法在分類學上的地位。DNA介導的物種鑒定重點在于種內(nèi)和種間遺傳差異的限定。種內(nèi)和種間變異范圍都不明確,而且不同物種的變異范圍也不相同。就線粒體DNA而言,種內(nèi)差異通常小于種間差異,但也有一些特例[34]。由于不同的種系間線粒體基因的水平轉(zhuǎn)移,近緣物種共享線粒體多態(tài)性以及同屬物種間的線粒體DNA修剪機制,都將造成快速分化物種與通過雜交產(chǎn)生的新物種在鑒定上的困難。本研究結(jié)果表明,DNA條形碼在簾蛤目貝類的鑒定中有一定的適用性,但仍有許多問題需要解決。

[1] 徐鳳山.中國海洋雙殼類軟體動物[M].北京:科學出版社, 1997:100-123.

[2] 齊鐘彥.中國經(jīng)濟軟體動物[M]北京:中國農(nóng)業(yè)出版社,1996: 218-220.

[3] 齊鐘彥.黃渤海的軟體動物[M]北京:農(nóng)業(yè)出版社,1987:146-150.

[4] Hebert PD N,Ratnasingham S,DeWard JR.Cytochrome c oxidase subunit Idivergences among closely related species[J].Barcoding Animal Life,2003,270(1):96-99.

[5] Ball S L,Hebert P D N,Burian S K.Biological identifications of mayflies Ephemeroptera using DNA barcodes[J].Journal of the North American Benthological Society,2005,24(3):508-524.

[6] Maynard B T,Kerr L J,Mckiernan JM,et al.Mitochondrial DNA sequence and gene organization in the Australian blacklip abaloneHaliotis rubra(Leach)[J].Mar Biotechnol,2005,7(1):645-658.

[7] Milbury C A,Gaffney PM.Completemitochondrial DNA sequence of the eastern oysterCrassostrea virginica[J].Mar Biotechnol, 2005,7(2):697-712.

[8] Passamonti M,Boore J L,Scali V.Molecular evolution and recombination in gender associated mitochondrial DNAs of the Manila clamTapes philippinarun[J].Genetics,2003,164(2):603-611.

[9] 陳麗梅,孔曉瑜,喻子牛,等.3種蟶類線粒體16S rRNA和COI基因片段的序列比較及其系統(tǒng)學初步研究[J].海洋科學, 2005,29(8):27-32.

[10] Park JK,Kim W.TwoCorbicula(Corbieulidae:Bivalvia)mitochondrial lineages arewidely distributed in Asian freshwater environment[J].Molecular Phylogenetics and Evolution,2003,29 (3):529-539.

[11] Pfenninger K,Reinhardt F,Streit B.Evidence for cryptic hybridization between different evolutionary lineages of the invasive clam genusCorbicula(Veneroida,Bivalvia)[J].Evol Biol, 2002,15(1):818-829.

[12] Smith P J,Mcveagh SM,Won Y,et al.Genetic heterogeneity among New Zealand speciesof hydrothermal ventmussels(Mytilidae:Bathymodiolus)[J].Marine Biology,2004,144(1):537-545.

[13] Therriault TW,Docker M F,Rlova M L,et al.Molecular resolution of the family Dreissenidae(Mollusca:Bivalvia)with emphasis on Ponto-Caspian species,including first report ofMytilopsis leucophaeatain the Black Sea basin[J].Molecular Phylogenetics and Evolution,2004,30(1):479-489.

[14] King TL,EackiesM S,Gjetvaj B,etal.Intraspecific phylogeography ofLasmigona subviridis(Bivalvia:Unionidae):conservation implications of range discontinuity[J].Molecular Ecology,1999, 8:S65-S78.

[15] Nagashima K,Sato M,Kawamata K,et al.Genetic structure of Japanese scallop population in Hokkaido analyzed by mitochondrial haplotype distribution[J].Mar Biotechnol,2005,7:1-10.

[16] Stepien C A,Morton B,Dabrowska K A,et al.Genetic diversity and evolutionary relationships of the troglodytic living fossilCongeri akusceri(Bivlvia:Dreissenidae)[J].Molecular Ecology, 2001,10(1):1873-1879.

[17] An H S,Jee Y J,Min K S,et al.Phylogenetic analysis of six species of Pacific abalone(Haliotidae)based on DNA sequences of 16S rRNA and cytochrome oxidase subunit Imitochondrial genes [J].Mar Biotechnol,2005,7(1):373-380.

[18] Carlini D B,Young R E,Vecchione M.A molecular phylogeny of the Octopoda(Mollusca:Cephalopoda)evaluated in light ofmorphological evidence[J].Molecular Phylogenetics and Evolution, 2001,21(3):388-397.

[19] Donald KM,KennedyM,Spencer H G.The phylogeny and taxonomy of Austral monodontine topshells(Mollusca:Gastropoda: Trochidae),inferred from DNA sequences[J].Molecular Phylogenetics and Evolution,2005,37(2):474-483.

[20] Goffredi SK,Hurtado L A,Hallam S,et al.Evolutionary relationships of deepsea ventand cold seep clams(Mollusca:Vesicomyidae)of the species complex[J].Marine Biology,2003,142(2): 311-320.

[21] Lin X,Zheng X,Xiao S,etal.Phylogeny of the cuttlefishes(Mollusca:Cephalopoda)based on mitochondrial COIand 16S rRNA gene sequence data[J].Acta Oceanologica Sinica,2004,23(4): 699-707.

[22] Machordom A,Araujo R,Erpenbeck D,et al.Phylogeography and conservation genetics of endangered European Margaritiferidae (Bivalvia:Unionoidea)[J].Biological J of the Linnean Society, 2003,78(2):235-252.

[23] Minton R L,Lydeard C.Phylogeny,taxonomy,geneticsand global heritage ranks of an imperiled freshwater snail genusLithasia(Pleuroceridae)[J].Molecular Ecology,2003,12(2):75-87.

[24] Nikula R,V In L R.Phylogeography ofCerastodermaglaucum(Bivlvia:Cardiidae)across Europe:amajor break in the EasternMediterranean[J].Marine Biology,2003,143(2):339-350.

[25] Folmer O,Black M,Hoeh W,et al.DNA primers for amplification ofmitochondrial cytochrome c oxidase subunit I from divercemetazoan invertebrates[J].Mol Mar Biol Biotech,1994,3(3):294-299.

[26] Hebert PD,Cywinska A,Ball SL,et al.Biological identifications through DNA barcodes[J].Proc Biol Sci,2003,270:313-321.

[27] Kumar S,Tamura K,NeiM.MEGA 3:Integrated software formolecular evolutionary genetics analysis and sequence alignment [J].Briefings in Bioinformatics,2004,5(2):150-163.

[28] Rozas J,Sánchez-DelBarrio JC,Messeguer X.DnaSPDNA polymorphism analyses by the coalescent and other methods[J]. Bioinformatics,2003,19(2):2496-2497.

[29] 徐鋼春,魏廣蓮,李建林,等.基于線粒體DNA D-loop序列分析養(yǎng)殖刀鱭與湖鱭的遺傳多樣性[J].大連海洋大學學報, 2012,27(5):448-452.

[30] Boore JL.Complete DNA sequence of themitochondrial genome of the black chiton,Katharina tunicate[J].Genetics,1994,138: 423-443.

[31] Boore J L.Animal mitochondrial genomes[J].Nucleic Acids Res,1999,27(8):1767-1780.

[32] Lindgren R A.Molecular inference of phylogenetic relationships among Decapodiformes(Mollusca:Cephalopoda)with special focus on the squid order Oegopsida[J].Mol Phylogen Evol,2010, 56(1):77-90.

[33] Akasaki T,Nikaido M,Tsuchiya K,et al.Extensivemitochondrial gene arrangements in coleoid Cephalopoda and their phylogenetic implications[J].Mol Phylogen Evol,2006,38(3):648-658.

[34] Brown W M,George JM,Wilson A C.Rapid evolution of animal mitochondrial DNA[J].Proc Natl Acad Sci USA,1979,76(3): 1967-1971.

[35] Masaaki Y,Kazuhiko T,Hisetaka F.Phylogeny of selected Sepiidae(Mollusca,Cephalopoda)based on 12S,16S,and COI sequence,with comments on the taxonomy reliability of severalmorphological characters[J].Zoolog Sci,2006,23:341-351.

[36] Zheng X D,Xiao S,Chen B,et al.Study on COI and 16S rRNA gene sequences of the common Chinese cuttlefishSepiella maindroni[J].Transaction of the Chinese Society of Malacology, 2001,9:64-69.

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

[38] 程漢良,吳婷婷,夏德全,等.6種簾蛤科貝類及4個地理種群文蛤線粒體COI基因片段序列分析[J].海洋學報,2007 (5):109-116.

[39] Yancheng.Relationships of the genusMeretrix(Mollusca:Veneridae)based on mitochondrial COI gene sequences[J].Zoology Res,2009,30(3):233-239.

[40] Vences M,Thomas M,Bonett R M,et al.Deciphering amphibian diversity through DNA barcoding:chances and challenges[J]. Phil Trans R Soc B,2005,360:1859-1868.

[41] HajibabaeiM,Janzen D H,Burns JM,et al.DNA barcodes distinguish species of tropical lepidoptera[J].Proc Natl Acad Sci USA,2006,103(4):968-971.

[42] 孔曉瑜,姜艷艷,相建海,等.魁蚶線粒體16S rRNA和COI基因片段序列測定及其應(yīng)用前景[J].海洋科學,2001,25(12): 46-48.

DNA barcoding and molecular phylogeny of 16 species of economical shellfish in Veneroida along China Coast

WANG Lin-nan1,YAN Xi-wu1,QIN Yan-jie1,NIE Hong-tao1,NIU Hong-bo1,ZHANG Guo-fan2
(1.Engineering Research Center of Shellfish Culture and Breeding in Liaoning Province,College of Fisheries and Life Science,Dalian Ocean University,Dalian 116023,China;2.Research and Development Center of Marine Biotechnology,Institute of Oceanology,Chinese Academy of Sciences,Qingdao 266071,China)

Mitochondrial COI genes were sequenced from 60 individuals of 6 common species in Veneroida along China coast using the universal barcoding primers,and the 60 sequences were compared with the other 50 homologous sequences from 10 shellfish species in Veneroida derived from the GenBank.The results showed that there were Insert-Deletion(indel)sites in many COI sequences in Veneroida species,total indel siteswere 103,including COI sequence indels of up to 30 for Manila clamRuditapes variegateandKatelysia hiantina,and 27 forSolengrandis.On average,the content of A+T(64.2%)was found significantly higher than that of G+C(35.8%). The Kimera-2-parameter model indicated that there was mean distances of 0.388 4 in pairwise-species and 0.010 6 in species.According tomaximum parsimony and neighbor joining trees for all 110 sequences of Veneroida,there was very high consistency between molecular analysis and morphological classification on family,genus, and species levels.Consequently,the findings provided essential revision and supplement to morphology taxonomy.

Veneroida;COI gene;phylogeny;taxonomy;DNA barcoding

S917.4

A

2095-1388(2013)05-0431-07

2013-02-21

現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項資金資助項目 (CARS-48)

王琳楠 (1988-),女,碩士研究生。E-mail:xlypt@126.com

閆喜武 (1962-),男,博士,教授。E-mail:Yanxiwu2002t@163.com