波紋巴非蛤5個(gè)地理群體的遺傳多樣性分析

駱　軒,徐小偉,張鵬飛,虞晉晉,黃妙琴,柯才煥

(廈門大學(xué) 海洋與地球?qū)W院,福建省海洋生物資源開發(fā)利用協(xié)同創(chuàng)新中心,福建廈門361102)

?

波紋巴非蛤5個(gè)地理群體的遺傳多樣性分析

駱軒,徐小偉,張鵬飛,虞晉晉,黃妙琴,柯才煥*

(廈門大學(xué) 海洋與地球?qū)W院,福建省海洋生物資源開發(fā)利用協(xié)同創(chuàng)新中心,福建廈門361102)

采用形態(tài)特征和微衛(wèi)星分子標(biāo)記對(duì)采自菲律賓(PH)、泰國(guó)(TH)、福建云霄(YX)、海南儋州(HN)和廣西北海(BH)的5個(gè)地理群體的波紋巴非蛤(Paphiaundulata)進(jìn)行遺傳多樣性分析．形態(tài)差異分析結(jié)果顯示,國(guó)外群體和國(guó)內(nèi)群體形態(tài)差異明顯:2個(gè)國(guó)外群體(PH和TH)形態(tài)接近,聚為一支;國(guó)內(nèi)支系中,YX和HN群體的形態(tài)較近,先聚為一支,而BH群體的外部形態(tài)與前兩者差異較大.Hardy-Weinberg 平衡檢驗(yàn)檢測(cè)到所有群體在5個(gè)位點(diǎn)上均出現(xiàn)不同程度地偏離平衡.5個(gè)群體間遺傳分化指數(shù)(FST，0.000 5～0.183 1)、基因流值(Nm,2.230 46～1 066.23)以及分子方差分析(AMOVA,群體間的變異貢獻(xiàn)率為10.21%，p<0.000 1)結(jié)果表明,群體間存在中等程度的分化.基于遺傳距離的非加權(quán)組平均法(UPGMA)聚類分析顯示:2個(gè)國(guó)外群體和3個(gè)國(guó)內(nèi)群體分別聚為一支,再聚為一個(gè)整體;國(guó)內(nèi)支系中,YX和BH群體先聚為一體,再與HN群體聚類．研究結(jié)果揭示了波紋巴非蛤較高的遺傳多樣性水平和中等程度的遺傳分化水平,同時(shí)也暗示其現(xiàn)有的遺傳分化程度受到人類生產(chǎn)活動(dòng)的干擾．

波紋巴非蛤;遺傳多樣性;形態(tài)特征;微衛(wèi)星;地理群體

遺傳多樣性是生物多樣性的主要形式之一,多指種內(nèi)遺傳多樣性[1],是生命進(jìn)化和物種分化的基礎(chǔ),更是評(píng)價(jià)自然生物資源的重要依據(jù)[2]．遺傳多樣性的高低決定了一個(gè)物種抵御惡劣環(huán)境的能力和自身進(jìn)化的潛力．

波紋巴非蛤(Paphiaundulata)隸屬簾蛤科(Veneridae)[3],是我國(guó)東南沿海一種重要的經(jīng)濟(jì)貝類,在我國(guó)自然分布于浙江南部、福建、廣東、海南和廣西沿岸海域,俗稱“油蛤”、“花蚶”和“花甲螺”[4]．由于營(yíng)養(yǎng)價(jià)值高且經(jīng)濟(jì)效益好,波紋巴非蛤現(xiàn)已成為一個(gè)優(yōu)良的海水養(yǎng)殖品種．然而,近年來隨著需求的增加和出口量的增大,波紋巴非蛤所面臨的采捕壓力也逐年增大,其海區(qū)自然資源受到了嚴(yán)重破壞．目前,國(guó)內(nèi)外對(duì)波紋巴非蛤的研究還主要集中于生物學(xué)和生態(tài)學(xué)特性[4-7]、增養(yǎng)殖技術(shù)[8-9]、食品應(yīng)用[10]以及微衛(wèi)星標(biāo)記開發(fā)[11]等方面,有關(guān)其遺傳多樣性的研究鮮有報(bào)道．

通過對(duì)我國(guó)東南沿海各個(gè)波紋巴非蛤產(chǎn)區(qū)的資源調(diào)查發(fā)現(xiàn),市售波紋巴非蛤來源混雜,其中廣西北海(BH)、海南儋州(HN)及福建云霄(YX)幾處的波紋巴非蛤幾乎占據(jù)了整個(gè)波紋巴非蛤的市場(chǎng),并且相互間存在著交換和苗種貿(mào)易，可見人類水產(chǎn)活動(dòng)對(duì)其自然狀態(tài)的干擾強(qiáng)度日漸增大．這種非自然的交流對(duì)種質(zhì)保持是不利的,長(zhǎng)此以往將影響我國(guó)東南沿海波紋巴非蛤的遺傳多樣性水平．因此，充分認(rèn)識(shí)波紋巴非蛤群體的遺傳多樣性現(xiàn)狀對(duì)其種質(zhì)資源的合理開發(fā)、利用和保護(hù)具有重要意義．

貝類的形態(tài)學(xué)差異建立在一定的遺傳基礎(chǔ)之上,因而不同地理種群的表型差異也能在一定程度上體現(xiàn)其遺傳多樣性．雖然已有關(guān)于波紋巴非蛤形態(tài)差異分析的報(bào)道[12],但結(jié)合分子標(biāo)記分析能更充分反映群體內(nèi)和群體間的遺傳變異水平．微衛(wèi)星DNA標(biāo)記具有多態(tài)性高、重復(fù)性好、共顯性等優(yōu)勢(shì),已廣泛應(yīng)用于長(zhǎng)牡蠣(Crassostreagigas)[13]、蝦夷扇貝(Patinopectenyessoensis)[14]、中國(guó)蛤蜊(Mactrachinensis)[15]、櫛江珧(Atrinapectinata)[16]、香港巨牡蠣(Crassostreahongkongensis)[17]和魁蚶(Scapharcabroughtonii)[18]等海洋雙殼貝類的遺傳多樣性研究．本研究基于形態(tài)特征與微衛(wèi)星標(biāo)記,分析了波紋巴非蛤的遺傳多樣性,以期為其種質(zhì)資源的利用和保護(hù)提供基礎(chǔ)數(shù)據(jù)和科學(xué)依據(jù)．

1　材料與方法

1.1樣品采集與處理

實(shí)驗(yàn)樣品于2011 年12月至 2013年12月分別采自泰國(guó)(TH)、菲律賓(PH)、YX、HN和BH自然海區(qū),各采樣點(diǎn)地理位置如圖1所示．每個(gè)群體包含的個(gè)體數(shù)不少于30,活體樣本經(jīng)冰塊保鮮運(yùn)回實(shí)驗(yàn)室后,隨機(jī)選取外殼完好的個(gè)體,取其斧足及閉殼肌肌肉樣保存于95%(體積分?jǐn)?shù))乙醇中以備DNA的提取,殼樣經(jīng)洗凈擦干后以備測(cè)量．

圖1　采樣點(diǎn)的地理位置Fig.1Locations of the sampling stations

1.2表觀形態(tài)參數(shù)的測(cè)量

采用游標(biāo)卡尺作為測(cè)量工具(精確到0.01 mm),獲得殼長(zhǎng)(LS,前端到后端的距離)、殼高(HS,殼頂?shù)礁咕壍木嚯x)、殼寬(WS,左右兩殼面間最大的距離)、韌帶長(zhǎng)(LH,鉸合部背面褐色幾丁質(zhì)韌帶的長(zhǎng)度)、前端-腹緣距(AF,前端到腹緣的距離)、后端-腹緣距(AP,后端到腹緣的距離)6個(gè)數(shù)量性狀,各測(cè)量部位見圖2．參考劉建勇等[12]用電子天平(精確到0.01 g)稱得殼質(zhì)量(mS,去除軟體部分并吸干表面水分后的殼質(zhì)量)．

圖2　波紋巴非蛤形態(tài)學(xué)測(cè)量Fig.2Morphological measurements of P. undulata

1.3基因組DNA的提取

采用北京天根生化科技有限公司的海洋動(dòng)物組織基因組DNA提取試劑盒提取波紋巴非蛤的基因組DNA,用1%(質(zhì)量分?jǐn)?shù))的瓊脂糖凝膠電泳檢測(cè)其純度,用微量紫外分光光度計(jì)(ND 2000)檢測(cè)其濃度,隨后將其稀釋到適合PCR的最佳濃度,置于4 ℃冰箱中保存?zhèn)溆茫?/p>

1.4微衛(wèi)星引物的合成

本研究采用8對(duì)能夠同時(shí)在5個(gè)波紋巴非蛤群體中進(jìn)行有效擴(kuò)增、結(jié)果穩(wěn)定且重復(fù)性高的微衛(wèi)星引物(表1),并由上海生工生物工程股份有限公司合成．

1.5PCR擴(kuò)增及產(chǎn)物檢測(cè)

每個(gè)群體隨機(jī)選取32個(gè)個(gè)體用于PCR反應(yīng)．PCR體系為20 μL,內(nèi)含10×PCR Buffer(含Mg2+) 2 μL,Taq酶(5 U/μL) 0.16 μL,正、反向引物(10 mmol/L)各0.8 μL,dNTPs(10 mmol/L) 1.6 μL,用ddH2O補(bǔ)足至20 μL．PCR反應(yīng)條件為:94 ℃預(yù)變性5 min;94 ℃變性30 s,退火30 s,72 ℃延伸30 s,30個(gè)循環(huán);72 ℃延伸7 min;10 ℃保存．將PCR產(chǎn)物送至上海生工生物工程股份有限公司,使用ABI公司的3730測(cè)序儀進(jìn)行毛細(xì)管電泳測(cè)序分析．

1.6數(shù)據(jù)分析

1.6.1形態(tài)數(shù)據(jù)分析

表1　8對(duì)微衛(wèi)星引物的特征

Tab.1　Characterization of 8 microsatellite primer pairs

引物名稱引物序列(5'→3')退火溫度/℃核心重復(fù)序列目的片段長(zhǎng)度/bpP15[11]F:TATGGTGTGGTGTGGTGTGG56(TGGCG)694～146R:CACTCGAACACGTCAACGCP22[11]F:TGTAAATAATGTATGAATGTATGTGCG54(GT)8264～356R:AGTGCTCTTCTGGTGATGCCP34[11]F:GGTATGGTGTGGTGTGGTGT57(TGGCG)690～144R:ACTCGAACATGTCAACGCCP48[11]F:TCAGGAGCATTTTGGTTTTATG50(GT)8(GTGC)9N(GC)6N(GT)22202～294R:TGTGCTACCTACCTACCGACCP51[11]F:TGGATTGGTGTGTCGGTAGA54(GT)8(GTGC)9N(GC)6N(GT)23180～276R:TTGGCTATCGCTGTAGTTATCAP67[11]F:TGTGTGTCGGTAGGTCGGTA54(AT)7262～288R:TCCATGTTATGGTTCAAGGAAGP79[19]F:GCCAAGGGACACTACACAGAGGG60(CA)3180～219R:TGGGGATGGGACAACTGAAAATGP8R[20]F:ACAAGGTGATGTGAGGTG55(AT)5139～168R:TCTGTCATCAGTGAAGGC

為消除不同群體波紋巴非蛤規(guī)格大小對(duì)形態(tài)特征造成的影響,選用比例性狀作為形態(tài)特征校正值來衡量各群體間的差異,共選用5個(gè)比例性狀(分別為HS/LS，WS/LS、LH/LS、mS、LS和AF/AP),采用SPSS 20.0軟件進(jìn)行分析．利用聚類分析、主成分分析和判別分析等多種多元統(tǒng)計(jì)分析方法對(duì)5個(gè)波紋巴非蛤群體的形態(tài)差異進(jìn)行研究．

1.6.2微衛(wèi)星數(shù)據(jù)分析

利用POPGENE32[21]計(jì)算各位點(diǎn)的Shannon′s多態(tài)信息指數(shù)(PIC)、等位基因數(shù)(Na)、有效等位基因數(shù)(Ne)、觀測(cè)雜合度(Ho)、期望雜合度(He)、固定指數(shù)(Fis)、Hardy-Weinberg遺傳偏離指數(shù)(HWD)、群體間的遺傳分化系數(shù)(FST)、群體間的遺傳距離(GD)和遺傳相似性系數(shù)(GI)．采用MEGA 6.0[22]軟件根據(jù)已獲得的GD,對(duì)5個(gè)群體進(jìn)行聚類分析構(gòu)建親緣關(guān)系圖．

利用Arlequin 3.5軟件[23]對(duì)5個(gè)群體間和群體內(nèi)的遺傳變異進(jìn)行分子方差分析(analysis of molecular variance,AMOVA),探討群體間的遺傳結(jié)構(gòu)．

2　結(jié)果與分析

2.1形態(tài)特征分析

2.1.1表觀形態(tài)參數(shù)測(cè)量

5個(gè)波紋巴非蛤群體的表觀形態(tài)特征參數(shù)測(cè)量結(jié)果如表2所示．

2.1.2聚類分析

聚類分析(圖3)表明：國(guó)內(nèi)3個(gè)群體和國(guó)外2個(gè)群體差異明顯,各自聚成一支；其中,YX群體與HN群體距離最短,TH群體與PH群體次之,說明兩組形態(tài)都比較接近；而BH群體與YX、HN群體趨異程度較大．

2.1.3主成分分析

對(duì)所有樣本的5個(gè)比例性狀進(jìn)行主成分分析后,得到2個(gè)主成分PC1和PC2,貢獻(xiàn)率分別為38.114%和29.221%,累計(jì)貢獻(xiàn)率為67.335%,即2個(gè)主成分可以解釋不同群體之間形態(tài)差異的67.335%(表3)．在PC1中,HS/LS的影響最大,貢獻(xiàn)率為 87.0%;WS/LS次之,貢獻(xiàn)率為76.6%．PC2中mS/LS影響最大,貢獻(xiàn)率為79.4%．

PC1和PC2的散布圖見圖4．由圖可見5個(gè)群體大體分成2簇：左上方主要為2個(gè)國(guó)外群體的重疊區(qū),幾乎難以區(qū)分；右下方則為3個(gè)國(guó)內(nèi)群體的重疊區(qū),其中YX與HN重疊區(qū)明顯,而BH群體與前兩者趨異程度更大.該分布規(guī)律與聚類分析結(jié)果一致．

表2　波紋巴非蛤5個(gè)地理群體樣本數(shù)及形態(tài)測(cè)量指標(biāo)

Tab.2　The sample size and morphological indices of five P. undulata populations

群體nLS/mmHS/mmWS/mmLH/mmAF/mmAP/mmmS/gTH5045.05±2.0525.57±1.2013.80±0.8113.19±0.9722.06±1.5629.31±1.664.39±0.87PH5046.20±2.3826.17±1.2814.04±0.7713.75±0.9723.26±1.3929.90±1.964.27±0.63YX3537.99±1.9322.03±1.1612.58±0.6311.34±0.8419.74±1.5823.89±1.442.70±0.45HN3540.11±1.8422.90±1.0313.32±0.6611.81±0.7520.85±1.3925.50±1.633.35±0.50BH3643.12±2.5724.11±1.5813.65±0.8611.56±0.9922.34±1.5827.40±2.143.15±0.70

2.1.4判別分析

通過逐步判別分析,篩選出4個(gè)比例性狀的特征值建立5個(gè)波紋巴非蛤群體的判別函數(shù)．

X1、X2、X3、X4分別代表HS/LS、WS/LS、LH/LS、mS/LS4個(gè)比例性狀的特征值．

5個(gè)群體的費(fèi)歇爾判別公式如下:

圖3　波紋巴非蛤5個(gè)地理群體聚類分析圖Fig.3Diagram of cluster analysis of five P. undulata populations

YPH=3 289.824X1+123.88X2+448.522X3-

652.233X4-988.107，

YTH=3 288.987X1+110.671X2+482.855X3-

692.221X4-989.962，

YYX=3 267.434X1+456.939X2+465.905X3-

968.773X4-1 060.37，

YHN=3 176.063X1+485.317X2+442.205X3-

858.268X4-1 018.68，

YBH=3 186.31X1+409.991X2+346.032X3-

868.656X4-972.148．

表3　波紋巴非蛤形態(tài)特征的主成分的負(fù)荷值和貢獻(xiàn)率

圖4　波紋巴非蛤5個(gè)地理種群的第一、二主成分的散布圖Fig.4Scatter diagram for PC1 and PC2 of five P. undulata populations

判別分析結(jié)果如表4所示：判別準(zhǔn)確率P1為54.0%～80.0%,判別準(zhǔn)確率P2為18.0%～52.0%，5個(gè)群體的綜合判別準(zhǔn)確率為65.0%．

表4　波紋巴非蛤5個(gè)群體的判別分析結(jié)果

2.2微衛(wèi)星數(shù)據(jù)分析

2.2.1群體內(nèi)的遺傳多樣性

本研究采用8個(gè)微衛(wèi)星分子標(biāo)記對(duì)5個(gè)波紋巴非蛤群體共160個(gè)個(gè)體進(jìn)行遺傳多樣性分析．由表5可知:5個(gè)群體的遺傳差異較小,平均Na為7.500～12.252,平均Ne為3.169～3.588,平均Ho為 0.963～0.986,平均He為0.670～0.706,平均PIC為1.391～1.648．綜合各遺傳多樣性指數(shù)看,5個(gè)群體中TH群體的遺傳多樣性最高(平均He為 0.706,平均PIC為1.648),HN群體的遺傳多樣性最低(平均He為 0.679,平均PIC為1.391)．

Hardy-Weinberg平衡檢測(cè)得到的HWD表明,所有群體除了在P67、P79、P8R 3個(gè)位點(diǎn)上表現(xiàn)為平衡狀態(tài)外,在其他5個(gè)位點(diǎn)均表現(xiàn)為偏離平衡,表明5個(gè)波紋巴非蛤地理群體均不同程度地偏離平衡狀態(tài)．

表5　8個(gè)微衛(wèi)星位點(diǎn)在5個(gè)波紋巴非蛤群體上的遺傳多樣性及Hardy-Weinberg平衡檢驗(yàn)

2.2.2群體間的遺傳分化及基因流(Nm)分析

AMOVA分析結(jié)果(表6)表明,10.21%的遺傳差異來源于群體間,且達(dá)到極顯著水平(p<0.000 1),17.77%來自于群體內(nèi),個(gè)體間的遺傳變異占總遺傳變異的72.02%．5個(gè)群體間兩兩分化的FST和Nm見表7．FST介于0.000 5～0.183 1之間,除YX與BH群體外,其余兩兩群體間的遺傳分化均達(dá)到極顯著水平(p<0.000 1)．3個(gè)國(guó)內(nèi)群體(YX、HN和BH)和2個(gè)國(guó)外群體(PH和TH)兩兩比較的FST比國(guó)內(nèi)群體間兩兩比較的FST更大,表明國(guó)內(nèi)和國(guó)外的波紋巴非蛤群體遺傳分化顯著．YX與BH群體的FST最小,且差異不顯著．5個(gè)群體間兩兩的Nm均為正值,YX與BH群體間的估測(cè)值高達(dá)1 066.23,其次為HN與BH群體,Nm為9.955 93,表明5個(gè)群體間均存在較頻繁的基因交流,分化水平有限(FST最大僅為0.183 1)．

表6　波紋巴非蛤5個(gè)地理群體的AMOVA

注：Va,Vb,Vc表示差異顯著；a表示p< 0.000 1．

2.2.3遺傳距離與聚類分析

按照1978年Nei[24]提出的方法計(jì)算出5個(gè)群體之間的GD和GI,結(jié)果(表7)表明:YX群體和BH群體的GD最近(0.008 0),遺傳一致性最高(GI為0.992 0);PH群體和HN群體的GD最遠(yuǎn)(0.115 3),遺傳一致性最低(GI為0.891 1)．

表7　波紋巴非蛤5個(gè)地理群體間的遺傳分化系數(shù)、基因流、遺傳相似性指數(shù)和遺傳距離

注：**表示p< 0.000 1.

根據(jù)各群體之間的GD構(gòu)建UPGMA 樹,聚類分析結(jié)果(圖5)表明：國(guó)外的PH與TH群體GD較近聚成一支；國(guó)內(nèi)的BH與YX群體GD較近,首先聚成一支,再與HN群體聚為一類．

3　討　論

3.15個(gè)地理群體的形態(tài)差異

本研究中,波紋巴非蛤5個(gè)地理群體形態(tài)特征的主成分分析與聚類分析結(jié)果基本一致,5個(gè)群體大致分成兩簇,簇內(nèi)群體間的主成分分別存在不同程度的重疊區(qū)域．國(guó)外2個(gè)群體(PH與TH)重疊率較高,形態(tài)難以區(qū)分，國(guó)內(nèi)3個(gè)群體中YX與HN的重疊明顯,而BH群體與二者趨異較大,該結(jié)果與劉建勇等[12]對(duì)波紋巴非蛤的群體分析結(jié)果一致,說明幾個(gè)地理群體之間在形態(tài)上既相似又存在明顯的差異．然而,本研究構(gòu)建的2個(gè)主成分的累計(jì)貢獻(xiàn)率僅為67.335%,不符合累計(jì)貢獻(xiàn)率大于或等于85%的要求,說明用幾個(gè)相互獨(dú)立的因子來概括波紋巴非蛤不同群體間的形態(tài)差異并不適宜．

本研究的判別分析中,HS/LS、WS/LS、LH/LS和mS/LS4個(gè)比例性狀進(jìn)入判別式,具有較大的判別貢獻(xiàn)．對(duì)5個(gè)群體的判別準(zhǔn)確率P1為54.0%～80.0%,判別準(zhǔn)確率P2為18.0%～52.0%,其綜合判別率為65.0%,低于翡翠貽貝(Pernavirdis(Linnaeus))[25]、毛蚶(Scapharcasubcrenata(Lischke))[26]、真曲巴非蛤(Paphiaeuglypta)[27]和縊蟶(Sinonovaculaconstricta)[28].但判別準(zhǔn)確率也受判別種群跨度和差異復(fù)雜程度影響,低判別準(zhǔn)確率可能與種內(nèi)判別有關(guān)．判別結(jié)果與實(shí)際結(jié)果基本相符,說明所構(gòu)建的判別函數(shù)是基本可靠的．

3.25個(gè)地理群體的遺傳多樣性

本研究所選用的8個(gè)微衛(wèi)星位點(diǎn)的PIC為0.831～2.229,平均值為1.391～1.648,根據(jù)Botstein等[29]提出的標(biāo)準(zhǔn),屬于高多態(tài)性(PIC>0.5),表明所用標(biāo)記可提供豐富的遺傳信息,適用于波紋巴非蛤遺傳多樣性的分析；He和Na能體現(xiàn)出群體的遺傳變異程度,He越高則遺傳變異越大,遺傳信息含量豐富．波紋巴非蛤5個(gè)地理群體的平均He為0.670～0.706,平均Na為7.500～12.252,高于魁蚶[30]和縊蟶[31]群體,與皺紋盤鮑(Haliotisdiscushannai)[32]和合浦珠母貝(Pinctadamartensii)[33-34]相似,表明其遺傳多樣性仍處于較高水平．

3.3形態(tài)聚類與微衛(wèi)星數(shù)據(jù)聚類的比較

一般認(rèn)為,在聚類分析中,聚類的先后順序反映了物種或種群間親緣關(guān)系的遠(yuǎn)近．本研究中,形態(tài)聚類和微衛(wèi)星數(shù)據(jù)聚類的結(jié)果具有一定的相似性,均分為2個(gè)獨(dú)立的分支,分別對(duì)應(yīng)國(guó)外群體和國(guó)內(nèi)群體,這說明形態(tài)學(xué)表型的差異是建立在一定的遺傳基礎(chǔ)之上的．然而,二者也有細(xì)微差別:在國(guó)內(nèi)群體分支中,形態(tài)聚類的結(jié)果是YX群體先與HN群體聚為一支,然后再與BH群體聚為一支；而微衛(wèi)星數(shù)據(jù)聚類的結(jié)果則是YX群體先與BH群體聚為一支,再與HN群體聚為一支．波紋巴非蛤?yàn)闇\海底棲貝類,其生活史要經(jīng)歷幼蟲浮游期,BH和HN同處北部灣內(nèi),BH群體與HN群體間基因交流的可能性較大．因此,不論是形態(tài)聚類還是微衛(wèi)星數(shù)據(jù)聚類,理論上,都應(yīng)該是BH群體與HN群體間先聚為一支,而后再與YX群體聚為一支．然而,形態(tài)聚類和微衛(wèi)星數(shù)據(jù)聚類的結(jié)果均未與理論情況相符合．究其原因,可能是由于波紋巴非蛤在我國(guó)廣泛分布于福建淺海內(nèi)灣、廣東沿岸、雷州半島西面、海南淺海和北部灣海域,是一種具有較高市場(chǎng)地位的經(jīng)濟(jì)貝類[4],而BH、HN和YX等作為海水貝類養(yǎng)殖發(fā)達(dá)區(qū),相互之間存在頻繁的人為引種、貿(mào)易及增養(yǎng)放流等活動(dòng)[34-35]．因此,波紋巴非蛤國(guó)內(nèi)群體的遺傳特點(diǎn)極有可能受到了人為活動(dòng)的影響．

雖然波紋巴非蛤5個(gè)地理群體間存在一定程度的遺傳分化,但從兩兩群體間的FST值(0.000 5～0.183 1)以及群體間的Nm(YX與BH群體間的估測(cè)值高達(dá)1 066.23)可知,波紋巴非蛤各群體間存在著充分的基因交流．而YX和BH群體間的GD最近(0.008 0),GI最高(0.992 0),更是從側(cè)面說明廣西和福建兩省的水產(chǎn)貿(mào)易活動(dòng)較為頻繁,以至于波紋巴非蛤的BH和YX群體的GD相對(duì)更近．

從形態(tài)聚類和微衛(wèi)星數(shù)據(jù)聚類特點(diǎn)出發(fā),本研究推斷波紋巴非蛤的遺傳多樣性特點(diǎn)在一定程度上受到了人為活動(dòng)的干擾．盡管本研究監(jiān)測(cè)到的波紋巴非蛤遺傳多樣性仍處于較高水平,但考慮到人類活動(dòng)對(duì)海洋生物遺傳多多樣性的影響潛力,建議在在今后資源管理和利用中,應(yīng)加強(qiáng)對(duì)波紋巴非蛤的遺傳資源的保護(hù),維持原種特性,避免因人為的基因滲透導(dǎo)致遺傳多樣性的丟失,以保證波紋巴非蛤種質(zhì)資源的可持續(xù)利用．

[1]GRAY J S.Marine biodiversity:patterns,threats and conservation needs[J].Biodiversity & Conservation,1997,6(1):153-175.

[2]崔朝霞,張恒,宋林生,等.中國(guó)重要海洋動(dòng)物遺傳多樣性的研究進(jìn)展[J].生物多樣性,2011,19(6):815-833.

[3]莊啟謙.中國(guó)動(dòng)物志:軟體動(dòng)物門[M].北京:科學(xué)出版社,2001:278.

[4]左江鵬,黃旭忠,黃雪芬,等.波紋巴非蛤的生物學(xué)及養(yǎng)殖技術(shù)[J].水產(chǎn)科技,2008,35(4):16-18.

[5]陳堅(jiān),柯愛英,范景水,等.波紋巴非蛤生物學(xué)性狀及生態(tài)習(xí)性的初步觀察[J].浙江海洋學(xué)院學(xué)報(bào)(自然科學(xué)版),2008,26(3):343-346.

[6]LEETHOCHAVALIT S,CHALERMWAT K,UPA-THAME S,et al.Occurrence ofPerkinsussp. in undu-lated surf clamsPaphiaundulatafrom the Gulf of Thailand[J].Diseases of Aquatic Organisms,2004,60(2):165-171.

[7]馬慶濤,林清,李春曉,等.不同因素對(duì)深澳灣波紋巴非蛤?yàn)V水率的影響[J].熱帶生物學(xué)報(bào),2015,6(1):5.

[8]王萬東.波紋巴非蛤人工育苗技術(shù)的初步研究[J].養(yǎng)殖與飼料,2010(8):14-17.

[9]黃松木,方火順,吳和平,等.云霄礁美海區(qū)波紋巴非蛤生物學(xué)特性及增殖研究[J].福建水產(chǎn),1984(3):14-19.

[10]關(guān)志強(qiáng),李敏,宋小勇,等.凍藏條件對(duì)文蛤和波紋巴非蛤組織結(jié)構(gòu)影響的實(shí)驗(yàn)研究[J].食品與發(fā)酵工業(yè),2007,33(1):147-153.

[11]郭昱嵩,謝子強(qiáng),鄧岳文,等.利用 FIASCO 技術(shù)進(jìn)行波紋巴非蛤微衛(wèi)星 DNA 標(biāo)記分離與篩選的研究[J].中國(guó)農(nóng)學(xué)通報(bào),2012,28(17):160-164.

[12]劉建勇,吳繼興,孫成波.我國(guó)東南沿海5個(gè)波紋巴非蛤 (Paphiaundulata)地理群體的形態(tài)差異分析[J].海洋與湖沼,2010,41(1):114-120.

[13]LI Q,YU H,YU R.Genetic variability assessed by mi-crosatellites in cultured populations of the Pacific oyster (Crassostreagigas) in China[J].Aquaculture,2006,259(1):95-102.

[14]AN H S,LEE Y G,PARK J Y,et al.Genetic characte-rization of four East Asian giant scallop (Mizuhopectenyessoensis) populations using microsatellite markers[J].Aquaculture Research,2009,40(5):619-624.

[15]NI L H,LI Q,KONG L F.Microsatellites reveal fine-scale genetic structure of the Chinese surf clamMactrachinensis(Mollusca,Bivalvia,Mactridae) in Northern China[J].Marine Ecology,2011,32(4):488-497.

[16]AN H S,LEE J W,DONG C M.Population genetic structure of Korean pen shell (Atrinapectinata) in Korea inferred from microsatellite marker analysis[J].Genes & Genomics,2012,34(6):681-688.

[17]LI L,WU X Y,YU Z N.Genetic diversity and substantial population differentiation inCrassostreahongkongensisrevealed by mitochondrial DNA[J].Marine Genomics,2013,11:31-37.

[18]YU H,GAO S,CHEN A L,et al.Genetic diversity and population structure of the ark shellScapharcabrough-toniialong the coast of China based on microsatellites[J].Biochemical Systematics and Ecology,2015，58:235-241.

[19]孫宏,李佳凱,武祥偉,等.基于轉(zhuǎn)錄組測(cè)序的波紋巴非蛤微衛(wèi)星標(biāo)記研究[J].集美大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,20(4):260-265.

[20]閆喜武,虞志飛,秦艷杰,等.菲律賓蛤仔EST-SSRs標(biāo)記開發(fā)及不同地理群體遺傳多樣性[J].生態(tài)學(xué)報(bào),2011,31(15):4190-4198.

[21]YEH F C,YANG R,BOYLE T J， et al.POPGENE 32,Microsoft Windows-based freeware for population genetic analysis[M].Alberta:Molecular Biology and Biotechnology Centre,University of Alberta,2000:153-156.

[22]TAMURA K,STECHER G,PETERSON D,et al.MEGA6:molecular evolutionary genetics analysis version 6.0[J].Molecular Biology and Evolution,2013,30(12):2725-2729.

[23]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.

[24]NEI M.Estimation of average heterozygosity and genetic distance from a small number of individuals[J].Genetics,1978,89(3):583-590.

[25]王慶恒,師尚麗,杜曉東,等.翡翠貽貝三個(gè)野生種群形態(tài)性狀的差異分析[J].廣東海洋大學(xué)學(xué)報(bào),2009,29(6):7-11.

[26]宋菲菲.山東近岸毛蚶群體的形態(tài)學(xué)和遺傳學(xué)研究[D].青島:中國(guó)海洋大學(xué),2013:44-58.

[27]劉建勇,孫成波,李文,等.我國(guó)南海4個(gè)真曲巴非蛤自然種群的形態(tài)差異分析[J].熱帶生物學(xué)報(bào),2010,1(1):72-77.

[28]劉達(dá)博,牛東紅,姜志勇,等.縊蟶7群體的形態(tài)差異與判別分析[J].海洋漁業(yè),2009,31(4):363-368.

[29]BOTSTEIN D,WHITE R L,SKOLNICK M,et al.Construction of a genetic-linkage map in man using restriction fragment length polymorphisms[J].American Journal of Human Genetics,1980,32(3):314-331.

[30]田吉騰,劉志鴻,楊愛國(guó),等.魁蚶4個(gè)地理群體遺傳多樣性微衛(wèi)星分析[J].漁業(yè)科學(xué)進(jìn)展,2013,34(6):59-67.

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

[32]李莉,孫振興,楊樹德,等.用微衛(wèi)星標(biāo)記分析皺紋盤鮑群體的遺傳變異[J].遺傳,2006,28(12):1549-1554.

[33]曲妮妮.合浦珠母貝微衛(wèi)星DNA標(biāo)記分離與遺傳多樣性研究[D].武漢：華中農(nóng)業(yè)大學(xué),2008：42-52.

[34]彭敏,陳秀荔,蔣偉明,等.企鵝珍珠貝不同地理群體遺傳多樣性的fAFLP分析[J].水生生物學(xué)報(bào),2012,36(1):102-108.

[35]嚴(yán)加坤,楊愛國(guó),周麗青,等.基于線粒體16S rRNA基因研究5個(gè)櫛江珧野生群體的遺傳多樣性[J].海洋科學(xué),2013,37(2):36-42.

Genetic Diversity Analysis of Five Populations of Paphia undulata

LUO Xuan,XU Xiaowei,ZHANG Pengfei,YU Jinjin,HUANG Miaoqin,KE Caihuan*

(Fujian Collaborative Innovation Center for Exploitation and Utilization of Marine Biological Resources, College of Ocean & Earth Sciences,Xiamen University,Xiamen 361102,China)

Genetic diversity of five populations ofPaphiaundulatasampled from Philippines (PH),Thailand (TH),Fujian Yunxiao (YX),Hainan Danzhou (HN) and Guangxi Beihai (BH) was studied using the morphological variation analysis and 8 micro-satellite markers.The morphological variations showed significant differences between the Chinese group and the exotic group.Within the Chinese group,the samples from YX and HN populations were rather similar in morphology,whereas BH population differed greatly from the above two populations.Meanwhile,significant departures from Hardy-Weinberg equilibrium were observed in 5 of the 8 microsatellite loci in all the 5 populations.The lines ofFST(0.000 5 to 0.183 1),Nm(2.230 46 to 1 066.23) and AMOVA analysis (10.21% of genetic variation was distributed within groups with the significant value,p<0.000 1) across all populations indicated that there was medium level of divergence among the 5P.undulatapopulations.The UPGMA clustering tree based on genetic distance demonstrated that the populations of PH and TH clustered into one group,and the other three Chinese populations clustered into one other group.Within the Chinese group,the HN population clustered into one group,while the populations of YX and BH clustered together,which indicated that the YX population was more closely related to the BH population.In this study,high levels of genetic diversity within populations and moderate levels of genetic differentiation among populations were presented forP.undulata.Our data also revealed that the genetic pattern ofP.undulatawas likely to be disturbed by human-mediated passive dispersal via aquaculture activities.These results would provide scientific basis for the conservation and reasonable utilization of natural resources ofP.undulata.

Paphiaundulata;genetic diversity;morphological traits;microsatellites;population

10.6043/j.issn.0438-0479.201605101農(nóng)業(yè)生產(chǎn)專題

2016-05-04錄用日期：2016-06-23

海洋公益性行業(yè)科研專項(xiàng)(201205021-2);國(guó)家貝類產(chǎn)業(yè)技術(shù)體系項(xiàng)目(CARS-48)

chke@xmu.edu.cn

駱軒,徐小偉,張鵬飛，等.波紋巴非蛤5個(gè)地理群體的遺傳多樣性分析[J].廈門大學(xué)學(xué)報(bào)(自然科學(xué)版)，2016,55(5)：637-645.

LUO X,XU X W,ZHANG P F,et al．Genetic diversity analysis of five populations ofPaphiaundulata[J].Journal of Xiamen University(Natural Science)，2016,55(5)：637-645．(in Chinese)

S 968.3

A

0438-0479(2016)05-0637-09