程雪艷, 滕爽爽 , 肖國(guó)強(qiáng) , 邵艷卿 , 張炯明 , 方 軍 , 柴雪良

(1. 溫州醫(yī)科大學(xué) 檢驗(yàn)醫(yī)學(xué)院 生命科學(xué)學(xué)院, 浙江 溫州 325035; 2. 浙江省海洋水產(chǎn)養(yǎng)殖研究所, 浙江 溫州 325005; 3. 浙江省近岸水域生物資源開(kāi)發(fā)和保護(hù)重點(diǎn)實(shí)驗(yàn)室 浙江 溫州 325001)

泥蚶(Tegillarca granosa)屬軟體動(dòng)物門(mén)(Mollusca)、雙殼綱(Bivalvia)、蚶目(Arcoida)、蚶科(Arcidae), 是一種棲息于沿海灘涂的廣溫廣鹽性雙殼貝類, 在我國(guó)主要分布于山東以南沿海[1], 是我國(guó)主要的養(yǎng)殖經(jīng)濟(jì)貝類之一, 具有很高的經(jīng)濟(jì)價(jià)值。但近年來(lái), 過(guò)度釆捕和生態(tài)環(huán)境的改變對(duì)泥蚶的自然資源造成嚴(yán)重的破壞。人工養(yǎng)殖的泥蚶種質(zhì)大多是未經(jīng)遺傳改良的野生型, 在經(jīng)過(guò)多年累代養(yǎng)殖后,已出現(xiàn)種質(zhì)退化等制約泥蚶養(yǎng)殖業(yè)持續(xù)健康發(fā)展的現(xiàn)象[2]。因此, 有必要對(duì)泥蚶遺傳背景進(jìn)行深入調(diào)查研究, 為泥蚶自然資源的恢復(fù)和遺傳改良奠定理論基礎(chǔ)。

目前, 對(duì)于泥蚶群體的遺傳多樣性的研究涉及形態(tài)學(xué)水平[1,3]、同工酶水平[4-5], DNA分子水平[3,6-8]等方面。微衛(wèi)星標(biāo)記(microsatellite)廣泛分布于真核生物基因組中, 由 2～6個(gè)核苷酸的串聯(lián)重復(fù)片段構(gòu)成。又稱為簡(jiǎn)單重復(fù)序列(simple sequence repeats)或短串聯(lián)重復(fù)序列(simple tandem repeats, STRs)[9]。由于具有數(shù)量豐富、多態(tài)性豐富、共顯性標(biāo)記、實(shí)驗(yàn)重復(fù)性好以及結(jié)果可靠等優(yōu)點(diǎn)[10], 被廣泛應(yīng)用于物種的種群遺傳結(jié)構(gòu)和遺傳多樣性分析。在水產(chǎn)動(dòng)物種群遺傳結(jié)構(gòu)的分析中已被廣泛使用, 如海灣扇貝(Argopecten irradias)[11]、蝦夷扇貝(Patinopecten yessoensis)[12]、文蛤(Meretrix meretrix)[13]、長(zhǎng)牡蠣(Crassostrea gigas)[14]、縊蟶(Sinonovacula constricta)[15-16]、馬氏珠母貝(Pinctada fucata)[17]、斧文蛤(Meretrix lamarckii)[18]等。

本研究采集了韓國(guó)以及我國(guó)沿海地區(qū) 5個(gè)地理群體泥蚶, 利用本實(shí)驗(yàn)室開(kāi)發(fā)的泥蚶 SSR引物對(duì) 5個(gè)群體進(jìn)行比較分析, 對(duì)泥蚶的遺傳多樣性水平進(jìn)行評(píng)估, 為更好地開(kāi)發(fā)利用該物種生物資源提供理論基礎(chǔ)。

1 材料與方法

1.1 樣品采集

泥蚶(T. granosa)的 5個(gè)地理群體分別取自浙江樂(lè)清(ZJ), 廣西企沙(GX), 海南海口(HN), 山東日照(SD),韓國(guó)釜山(KR)等地海域, 5個(gè)地理群體均為自然群體,數(shù)量分別為 24、23、23、26、24個(gè)。樣品采集后活體運(yùn)回實(shí)驗(yàn)室, –70℃超低溫冰箱保存?zhèn)溆谩?/p>

1.2 基因組DNA提取

將所采集樣品分別解剖, 取閉殼肌組織?；蚪MDNA提取采用酚-氯仿法抽提[19], 對(duì)所提取DNA進(jìn)行純度、濃度以及完整性檢測(cè)后, 4℃保存?zhèn)溆谩?/p>

1.3 PCR擴(kuò)增與電泳檢測(cè)

本研究所用引物由本課題組開(kāi)發(fā)設(shè)計(jì)而成[20],引物序列及退火溫度見(jiàn)表1, 引物由上海生工生物工程有限公司合成。在Eppendorf擴(kuò)增儀上進(jìn)行PCR反應(yīng), 反應(yīng)體系: 10×PCR buffer(Mg2+plus)1.5 μl,2.5 mmol/L dNTP 1.2 μl, 10 μmol/L 上下游引物各 1.5 μl,Taq 酶(5 U/μl)0.075 μl, DNA 模板 100 ng, 滅菌水補(bǔ)足至15 μl。PCR反應(yīng)條件設(shè)置: 94℃預(yù)變性5 min;94℃變性30 s, 50℃退火30 s, 72℃延伸45 s, 30個(gè)循環(huán); 最后72℃延伸7 min。PCR結(jié)束后取4 μl擴(kuò)于2%瓊脂糖凝膠電泳檢測(cè)(1×TBE, 5 V/cm 恒壓),DuGreen染色, 凝膠成像系統(tǒng)觀察和拍照。PCR 擴(kuò)增結(jié)束后將產(chǎn)物于 4℃保存, 擴(kuò)增產(chǎn)物用 6%的變性聚丙烯酰胺凝膠電泳檢測(cè)。

表1 泥蚶17個(gè)微衛(wèi)星位點(diǎn)的特性Tab.1 Characteristics of 17 microsatellites for T. granosa

1.4 數(shù)據(jù)統(tǒng)計(jì)

2 結(jié)果

2.1 泥蚶5個(gè)地理群體遺傳參數(shù)分析

使用17對(duì)微衛(wèi)星引物對(duì)泥蚶5個(gè)不同地理群體共120個(gè)樣品進(jìn)行了DNA檢測(cè), 17對(duì)引物在5個(gè)泥蚶群體中分別檢測(cè)到 102(山東)、99(廣西)、91(海南)、86(韓國(guó))、99個(gè)(浙江)等位基因, 由于不同群體檢出的等位基因不一致, 共獲得115個(gè)不同等位基因。不同位點(diǎn)檢測(cè)的等位基因數(shù)為2～12個(gè)不等, 其中位點(diǎn)Teg39等位基因數(shù)最多, 為12個(gè)等位基因; Teg12次之, 為11個(gè)等位基因; Teg18最少, 為2個(gè)等位基因。其中Teg28在其他群體中為2～3個(gè)等位基因, 在韓國(guó)群體中為1個(gè)等位基因。

通過(guò)基因產(chǎn)生的頻率計(jì)算17個(gè)位點(diǎn)的PIC(表2),17個(gè)位點(diǎn)的平均多態(tài)性信息含量為0.142～0.857。其中Teg18最低, 為0.142; Teg28、Teg24、Teg8的多態(tài)性信息含量為0.176、0.273、0.377, 10個(gè)位點(diǎn)的多態(tài)性信息含量在0.5以上, 表現(xiàn)出高度多態(tài)性。5個(gè)群體的PIC為0.525～0.608, 說(shuō)明泥蚶5個(gè)群體的遺傳多樣性處于中等偏上水平。

5個(gè)地理群體的平均觀測(cè)雜合度 Ho為 0.430～0.516, 平均期望雜合度He為0.573～0.656, 其中山東群體的Ho和He均為最高。

在85個(gè)群體位點(diǎn)中, Teg28在韓國(guó)群體中的位點(diǎn)為單態(tài)的, 對(duì)84個(gè)多態(tài)性的群體位點(diǎn)進(jìn)行分析發(fā)現(xiàn),41個(gè)群體位點(diǎn)符合Hardy-Weinberg平衡(P>0.05), 而其余 43個(gè)(51.2%)群體位點(diǎn)由于雜合子缺失導(dǎo)致偏離Hardy-Weinberg平衡(P<0.05)。在 17個(gè)微衛(wèi)星位點(diǎn)中,Teg12、Teg17、Teg18、Teg19四個(gè)位點(diǎn)在5個(gè)群體中均沒(méi)有偏離Hardy-Weinberg平衡, Teg29、Teg37、Teg39三個(gè)位點(diǎn)在5個(gè)群體中均偏離Hardy- Weinberg平衡。

2.2 不同群體間遺傳結(jié)構(gòu)

使用Popgene 32(Version 1.31)軟件分析泥蚶5個(gè)地理群體遺傳距離及遺傳相似性指數(shù)。各群體間遺傳相似性指數(shù)以廣西群體和山東群體之間最大(0.936 5), 其次為浙江群體和韓國(guó)群體之間(0.015 9),廣西群體和韓國(guó)群體之間的遺傳相似性最小(0.842 8);相應(yīng)的, 各群體之間的遺傳距離以廣西群體和韓國(guó)群體之間最大(0.171 1), 以廣西群體和山東群體(0.065 6)、浙江群體和韓國(guó)群體之間(0.087 9)相對(duì)較小(表3)。利用MEGA6.0軟件對(duì)泥蚶 5個(gè)群體進(jìn)行聚類分析中, 獲得UPGMA聚類圖(圖1)。浙江群體和韓國(guó)群體親緣關(guān)系最近首先聚在一起, 然后與海南群體聚合; 廣西群體和山東群體親緣關(guān)系也較近,聚為一支, 然后與上面3個(gè)群體聚合。

圖1 基于泥蚶5個(gè)群體遺傳距離構(gòu)建的UPGMA聚類圖Fig.1 UPGMA dendrogram based on genetic distances of five populations of T. granosa

對(duì)5個(gè)泥蚶群體進(jìn)行了F-統(tǒng)計(jì)量以及基因流分析(表 4)。5個(gè)泥蚶群體中 Teg17、Teg19、Teg34位點(diǎn)的近交系數(shù) FIS為負(fù)值, 其余位點(diǎn)的 FIS為正值;FIS值為–0.009 9～0.738 2, 平均值為0.217 4。各位點(diǎn)的總?cè)后w近交系數(shù)為 0.000 7～0.774 9, 平均值為0.264 2。度量群體間遺傳差異程度的FST為0.014 2～0.164 4, 平均值為 0.059 8, 基因流 Nm為 1.271～17.343 9, 平均值為3.927 5, 不同位點(diǎn)的基因流值變化波動(dòng)較大。

利用Arlequin (version 3.1)軟件對(duì)泥蚶5個(gè)群體之間的FST分析, 結(jié)果如表5。泥蚶5個(gè)群體間的FST值在 0.012和 0.062之間, 其中韓國(guó)群體-山東群體(0.062)的FST值最高, 韓國(guó)群體-廣西群體(0.060), 海南群體-山東群體(0.051)的FST值僅次于韓國(guó)群體-山東群體。韓國(guó)群體-海南群體(0.012)、山東群體-廣西群體(0.016)的 FST值最低, 韓國(guó)群體-海南群體與山東群體-廣西群體之外的其他群體間的 FST值均達(dá)到顯著性差異。對(duì) 5個(gè)地理種群的泥蚶進(jìn)行 AMOVA分析, 結(jié)果見(jiàn)表 6, 發(fā)現(xiàn)種群間變異占 3.59%, 種群個(gè)體間變異占 96.41%, 結(jié)果表明引起種群總體變異的主要因素還是種群內(nèi)的遺傳變異。

表2 17個(gè)微衛(wèi)星位點(diǎn)在泥蚶5個(gè)地理群體上的遺傳多樣性分析Tab.2 Parameters of genetic diversity analyzed by seventeen loci in T. granosa five geographic populations

表3 泥蚶5個(gè)地理群體間的Nei’s遺傳相似性(上三角)和遺傳距離(下三角)Tab.3 Nei’s genetic identity(above diagonal) and genetic distance(below diagonal) among five populations of T. granosa

表4 17個(gè)微衛(wèi)星位點(diǎn)的固定指數(shù)及基因流Tab.4 Fixation index and gene flow of 17 polymorphic loci

表5 泥蚶5個(gè)群體樣本的遺傳分化系數(shù)FSTTab.5 The FST value among 5 populations of T. granosa

表6 泥蚶群體之間的AMOVA分析結(jié)果Tab.6 Analysis of AMOVA of T. granosa populations

用 IBDWS 3.23對(duì)泥蚶群體間的地理距離和遺傳距離的關(guān)系進(jìn)行分析, 結(jié)果如圖2所示, 發(fā)現(xiàn)泥蚶5個(gè)群體間的遺傳分化程度和地理距離之間沒(méi)有顯著的相關(guān)性。

3 討論

3.1 泥蚶不同地理群體的遺傳多樣性分析

圖2 地理距離和遺傳距離的關(guān)系圖Fig.2 Diagram of the relationship between geographic distance and genetic distance

遺傳多樣性是物種和某個(gè)群體長(zhǎng)期進(jìn)化的產(chǎn)物,是衡量一個(gè)種群種質(zhì)資源質(zhì)量的重要標(biāo)準(zhǔn)[22]。一個(gè)種群的遺傳多樣性越高, 其對(duì)生存環(huán)境的適應(yīng)能力就越強(qiáng), 進(jìn)化的潛力也越大[23]。

目前, 對(duì)于泥蚶群體的遺傳多樣性已有較多研究, 喻子牛等[24]采用水平淀粉凝膠電泳技術(shù)研究了三個(gè)泥蚶群體的等位基因酶遺傳變異, 在12種等位基因酶中共檢測(cè)到了27個(gè)基因位點(diǎn), 平均觀察雜合度為0.075, 平均有效等位基因數(shù)為1.372; 李太武等[6]采用隨機(jī)擴(kuò)增多態(tài)性DNA(RAPD)技術(shù)對(duì)5個(gè)泥蚶群體的遺傳多樣性進(jìn)行研究, 其平均期望雜合度為0.260 7; 呂振明等[4]采用垂直聚丙烯酰胺同工酶電泳技術(shù)對(duì) 7個(gè)泥蚶群體的遺傳多樣性進(jìn)行分析, 不同群體間位點(diǎn)數(shù)為 21～28, 平均期望雜合度為0.086～0.186, 平均觀測(cè)雜合度為 0.035～0.083。本實(shí)驗(yàn)采用微衛(wèi)星標(biāo)記技術(shù)對(duì) 5個(gè)泥蚶群體的遺傳多樣性進(jìn)行分析, 共檢測(cè)到115個(gè)不同等位基因, 5個(gè)群體的He為0.573～0.656, Ho為0.430～0.516, 明顯高于其他標(biāo)記的等位基因數(shù)和雜合度, 這與姚韓韓等[3]采用相同技術(shù)對(duì) 2個(gè)泥蚶群體的遺傳多樣性的分析結(jié)果相一致。較高的等位基因數(shù)和雜合度, 說(shuō)明微衛(wèi)星標(biāo)記技術(shù)相對(duì)于其他標(biāo)記技術(shù)可以揭示出較高水平的遺傳多樣性, 是研究物種遺傳多樣性的理想工具, 也有可能是因?yàn)楸緦?shí)驗(yàn)中選擇的17個(gè)位點(diǎn)均有較高的的理論多態(tài)信息含量。

17對(duì)微衛(wèi)星引物的多態(tài)信息含量為0.142～0.857,依據(jù)Botstein等[21]的劃分標(biāo)準(zhǔn), 本研究的17個(gè)位點(diǎn)中, 有5個(gè)位點(diǎn)為中度多態(tài)位點(diǎn), 10個(gè)位點(diǎn)為高度多態(tài)位點(diǎn), 這些微衛(wèi)星位點(diǎn)在泥蚶遺傳多樣性分析中可提供有效的信息。Teg28和Teg18的多態(tài)信息含量最低, 其中Teg28在韓國(guó)群體中為單態(tài)位點(diǎn), 分析原因可能是樣品數(shù)量不夠充足或者位點(diǎn)本身多態(tài)信息含量低, 可能在韓國(guó)群體中為保守位點(diǎn)。5個(gè)泥蚶群體的平均多態(tài)信息含量介于 0.525～0.608, 其中山東群體的平均多態(tài)信息含量最高, 為0.608; 5個(gè)群體的平均多態(tài)信息含量均處于中度偏上水平, 具有較高的遺傳多樣性水平。

群體內(nèi)Ho和He之間的平衡關(guān)系是由FIS反映的。當(dāng)FIS<0時(shí), 說(shuō)明雜合子過(guò)度; FIS>0時(shí), 說(shuō)明雜合子缺失, FIS越接近0, 說(shuō)明基因型的分布越接近平衡狀態(tài)[25]。本實(shí)驗(yàn)中, FIS除了 3個(gè)微衛(wèi)星位點(diǎn)外, 其余14個(gè)位點(diǎn)均大于0, FIS的平均值為0.217 4, 說(shuō)明了泥蚶群體中雜合子缺失, 海洋貝類的雜合子缺失現(xiàn)象在很多貝類中都有相關(guān)報(bào)道[26,27]。由于雜合子缺失導(dǎo)致 43個(gè)(51.2%)群體位點(diǎn)偏離 Hardy-Weinberg平衡, Yu等[27]認(rèn)為主要原因可能是無(wú)效等位基因的存在; 喻子牛等[24]認(rèn)為可能是由于近交衰退、Wahlund效應(yīng)、電泳圖譜的閱讀、自然選擇等原因造成的。

3.2 泥蚶群體的遺傳分化

遺傳距離是用來(lái)估計(jì)不同種群之間遺傳分化程度的一個(gè)指標(biāo), 根據(jù)其大小可以確定群體間親緣關(guān)系的遠(yuǎn)近。群體間的遺傳距離越大, 則基因型差異越大, 表明群體間的親緣關(guān)系越遠(yuǎn), 遺傳相似性越小;反之亦然[28]。本研究采用17個(gè)微衛(wèi)星位點(diǎn)對(duì)5個(gè)泥蚶群體進(jìn)行遺傳多樣性分析, 根據(jù)Nei's指數(shù)法對(duì)泥蚶 5個(gè)群體的遺傳結(jié)構(gòu)進(jìn)行遺傳距離和遺傳相似性統(tǒng)計(jì)分析, 結(jié)果表明廣西-山東群體遺傳相似度最大(0.936 5), 遺傳距離最小(0.065), 聚為一支。浙江-韓國(guó)群體遺傳相似度(0.915 9)僅次于廣西-山東群體,也聚為一支。海南群體與上述兩支群體間的遺傳相似度比較接近, 所以處于兩支群體之間。其中浙江-韓國(guó)群體遺傳距離比較接近, 在很多其他的研究中均有此結(jié)論[6,29], 5個(gè)群體間遺傳分化的水平的高低并不完全符合地理分布的特征, 分析原因可能是由于在泥蚶大規(guī)模養(yǎng)殖育種過(guò)程中, 不同地理群體之間互相引種[1,6], 造成不同地理群體之間交流頻繁,尤其是泥蚶養(yǎng)殖業(yè)比較發(fā)達(dá)的地區(qū), 致使群體間的遺傳分化程度并不完全符合地理分布特征。

遺傳分化系數(shù) FST(F-statistic)和基因流(gene flow)是反映種群間的遺傳分化程度和各種群間基因交流情況的重要指標(biāo)[30]。當(dāng)0

鄭文娟等[7]利用 COI標(biāo)記對(duì)泥蚶地理群體進(jìn)行種群結(jié)構(gòu)分析, 得到了我國(guó)沿海泥蚶已分化形成福建以南和以北兩大類群的結(jié)論; Ni等[8]利用COI、ITS標(biāo)記對(duì)中國(guó)東部海域和南海泥蚶類群進(jìn)行種群遺傳結(jié)構(gòu)分析, 提出了兩個(gè)進(jìn)化顯著單元的觀點(diǎn); 邵艷卿等[32]利用 AFLP對(duì)泥蚶種群結(jié)構(gòu)進(jìn)行分析, 提出未發(fā)現(xiàn)有效的南北間基因流屏障的結(jié)論; 劉春芳[33]利用 COI標(biāo)記對(duì)泥蚶 7個(gè)地理群體進(jìn)行分析, 發(fā)現(xiàn)泥蚶群體的遺傳分化主要來(lái)自于群體內(nèi)部, 與之前學(xué)者提出的南北分化且群體間的遺傳距離大于群體內(nèi)遺傳距離的觀點(diǎn)不甚相同。所以, 目前關(guān)于我國(guó)沿海泥蚶的南北區(qū)系分化的說(shuō)法依然存在爭(zhēng)議, 分析原因可能是由于取樣的時(shí)間不同, 而且取樣數(shù)量均較少, 不能代表整個(gè)地理群體的種質(zhì)資源情況; 而且不同類型的標(biāo)記方法特點(diǎn)不同, 都具有各自的優(yōu)勢(shì)和局限性, 檢測(cè)基因組的部位也不同, 可能會(huì)導(dǎo)致結(jié)果不一樣。

本研究利用17種微衛(wèi)星標(biāo)記對(duì)泥蚶5個(gè)地理群體進(jìn)行遺傳多樣性分析, 結(jié)果顯示泥蚶 5個(gè)群體的遺傳多樣性水平較高, 種質(zhì)資源良好, 對(duì)泥蚶養(yǎng)殖產(chǎn)業(yè)的良種選育、種質(zhì)資源保護(hù)等工作具有指導(dǎo)意義。但是, 本研究中采集的泥蚶的群體數(shù)和個(gè)體數(shù)有限, 為了更好地保護(hù)泥蚶的種質(zhì)資源, 我們應(yīng)加強(qiáng)對(duì)泥蚶群體資源的遺傳學(xué)研究, 制定一套健康科學(xué)的保護(hù)措施, 保證泥蚶種質(zhì)資源得到合理和可持續(xù)性的利用。

[1] 張永普, 林志華, 應(yīng)雪萍, 等.不同地理種群泥蚶的形態(tài)差異與判別分析[J]. 水產(chǎn)學(xué)報(bào), 2004, 28(3):339-342.

[2] 姚韓韓, 董迎輝, 林志華, 等. 泥蚶 4個(gè)快速生長(zhǎng)家系的遺傳變異分析[J]. 水產(chǎn)學(xué)報(bào), 2011, 35(3):340-347.

[3] 姚韓韓, 林志華, 周小龍, 等. 泥蚶奉化種群與韓國(guó)種群形態(tài)特征和遺傳結(jié)構(gòu)的差異分析[J]. 江蘇農(nóng)業(yè)科學(xué), 2014, 42(3): 183-186.

[4] 呂振明, 李太武, 蘇秀榕, 等. 泥蚶遺傳多樣性的研究[J]. 高科技通訊, 2005, 15(12): 104-110.

[5] 王日昕, 李太武, 呂振明, 等. 泥蚶不同地理居群同工酶變異及遺傳分化的研究[J]. 海洋與湖沼, 2005,36(3): 228-234.

[6] 李太武, 李成華, 宋林生, 等. 5個(gè)泥蚶群體遺傳多樣性的RAPD分析[J].生物多樣性, 2003, 11(2): 118-124.

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

[8] Ni G, Li Q, Kong LF, et al. Phylogeography of the bivalve Tegillarca granosa in coastal China:implications for management and conservation[J].Marine Ecology Progress Series, 2012, 452: 119-130.

[9] 劉華, 賈繼增. 指紋圖譜在作物品種鑒定中的應(yīng)用[J]. 作物品種資源, 1997, (2): 45-48.

[10] Powell W, Machray G C, Provan J. Polymorphism revealed by simple sequence repeats[J]. Trends Plant Sci, 1996, 1: 215-222.

[11] 李宏俊, 劉曉, 杜雪地, 等. 海灣扇貝微衛(wèi)星標(biāo)記開(kāi)發(fā)及其分離方式分析[J]. 海洋科學(xué), 2009, 33(12):4-8.

[12] 韓斐斐, 張繼彪, 李莉, 等. 蝦夷扇貝養(yǎng)殖群體及其子代的遺傳多樣性分析[J]. 海洋科學(xué), 2012, 36(9):1-8.

[13] Lu X, Wang H X, Liu B Z, et al. Microsatellite-based genetic and growth analysis for a diallel mating design of two stocks of the clam, Meretrix meretrix[J].Aquaculture research, 2012, 43(2): 260-270.

[14] 張琳琳, 李莉, 張國(guó)范. 長(zhǎng)牡蠣(Crassostrea gigas)EST 串聯(lián)重復(fù)序列的組成和分布[J]. 海洋科學(xué),2011, 35(4): 9-14.

[15] 劉博, 邵艷卿, 滕爽爽, 等. 樂(lè)清灣養(yǎng)殖縊蟶群體遺傳結(jié)構(gòu)的微衛(wèi)星標(biāo)記分析[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2012,28(2): 69-73.

[16] 劉博, 邵艷卿, 王侃, 等. 4 個(gè)縊蟶群體遺傳多樣性和系統(tǒng)發(fā)生關(guān)系的微衛(wèi)星分析[J]. 海洋科學(xué), 2013,37(8): 96-102.

[17] 湯健, 管云雁, 劉文廣, 等. 馬氏珠母貝家系遺傳結(jié)構(gòu)的微衛(wèi)星分析[J]. 海洋科學(xué), 2013, 37(8): 35-41.

[18] 滕爽爽, 方軍, 金文輝, 等. 4個(gè)斧文蛤群體微衛(wèi)星標(biāo)記的遺傳多樣性分析[J]. 中國(guó)水產(chǎn)科學(xué), 2015, 22(3):580-587.

[19] Sambrook J, Fritsch E, Maniatis T. Molecular Cloning:a Laboratory Manual[M]. third ed. New York : Cold Spring Harbor Laboratory Press, 2001.

[20] Liu B, Teng S S, Shao Y Q, et al. Isolation and characterization of 39 novel polymorphic EST-SSR loci for the blood clam, Tegillarca granosa[J]. Conservation Genet Resour, 2012, 4(2): 375-378.

[21] Botstein D, White R L, Skolnick M, et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. Am J Hum Genet,1980, 32(3): 314-331.

[22] 姚韓韓. 泥蚶野生群體與選育群體的分子遺傳基礎(chǔ)研究[D].上海: 上海海洋大學(xué), 2011.

[23] 鄧春興. 基于 COI基因序列的中國(guó)東南沿海四指馬鲅屬魚(yú)類的遺傳多樣性分析[D]. 廣州: 暨南大學(xué),2014.

[24] 喻子牛, 孔曉瑜, 楊銳, 等. 泥蚶等位基因酶遺傳變異研究[J]. 中國(guó)水產(chǎn)科學(xué), 1997, 4(5): 5-21.

[25] Weir B S, Cockerham C C. Estimation F-statistics for the analysis of population structure[J].Evolution, 1984,38: 1358-1370.

[26] 李剛, 金啟增, 姜衛(wèi)國(guó), 等. 合浦珠母貝和長(zhǎng)耳珠母貝的生化遺傳變異[J]. 遺傳學(xué)報(bào), 1985, 12(3):204-212.

[27] Yu H, Li Q. Genetic variation of wild and hatchery populations of the Pacific oyster Crassostrea gigas assessed by microsatellite markers[J]. Journal of Genetics and Genomics 2007, 34(12): 1114-1122.

[28] 張愛(ài)兵, 王正軍, 譚聲江, 等. 分子生態(tài)學(xué)重要概念-遺傳距離及其測(cè)定的理論研究概況[J]. 生態(tài)學(xué)報(bào),2002, 22(6): 943-949.

[29] 李成華, 李太武, 宋林生, 等. 福建南北泥蚶種內(nèi)分化的RAPD分析[J]. 動(dòng)物學(xué)研究, 2003, 24(5): 362-366.

[30] Allendorf F W. Isolation, gene flow and genetic differentiation among populations[J]. Genetic and Conservation, 1983, 5l-65.

[31] Weight S. Evolution and the genetics of population[M],Volume 4: variability within and among Natural Population Chicago. Chicago: university of Chicago press, 1987.

[32] 邵艷卿, 柴雪良, 董迎輝等.地理距離對(duì)基因流的影響: 6個(gè)泥蚶群體的AFLP分析[C]//中國(guó)動(dòng)物學(xué)會(huì)、中國(guó)海洋湖沼. 中國(guó)海洋湖沼學(xué)會(huì)貝類學(xué)會(huì)分會(huì)第十四次學(xué)會(huì)研討會(huì)論文集.青島: 中國(guó)海洋湖沼學(xué)會(huì),2009: 95.

[33] 劉春芳.中國(guó)近海蚶總科(Arcoidea)貝類的系統(tǒng)發(fā)育及泥蚶(Tegillarca granosa)不同群體遺傳多樣性研究[D]. 青島: 中國(guó)科學(xué)院海洋研究所, 2014.