祁劍飛，寧 岳，巫旗生，郭 香，鄭雅友，曾志南，江澤貴

(福建省水產(chǎn)研究所、福建省海洋生物增養(yǎng)殖與高值化利用重點(diǎn)實(shí)驗(yàn)室，福建 廈門 361013)

菲律賓蛤仔(Ruditapesphilippinarum)是一種廣溫、廣鹽、廣分布的海洋經(jīng)濟(jì)貝類，2017年全國(guó)蛤類產(chǎn)量達(dá)417.79萬噸，占全國(guó)海水養(yǎng)殖產(chǎn)量的20.88%[1]，其中菲律賓蛤仔產(chǎn)量超過300萬噸，是我國(guó)單品種養(yǎng)殖產(chǎn)量最高的海洋生物物種，其總產(chǎn)量在所有貝類物種中也居世界第二[2]。遍及全國(guó)沿海，特別是遼寧和山東的大規(guī)模菲律賓蛤仔養(yǎng)殖活動(dòng)得益于福建省大規(guī)模的苗種供應(yīng)[3]，據(jù)業(yè)內(nèi)估計(jì)，其苗種產(chǎn)量占全國(guó)的90%以上。但是這種長(zhǎng)期的異地苗種底播養(yǎng)殖方式會(huì)對(duì)當(dāng)?shù)匾吧N群的遺傳結(jié)構(gòu)產(chǎn)生一定影響[4]。因此對(duì)菲律賓蛤仔野生種群資源的調(diào)查和評(píng)估很有必要，既有利于保護(hù)野生種群資源，也可為雜交育種提供參考。而外部形態(tài)差異則是不同野生地理群體基于形態(tài)測(cè)量學(xué)的判定條件。

關(guān)于菲律賓蛤仔遺傳差異的報(bào)道[5-9]較多，而關(guān)于其不同地理群體形態(tài)差異的研究，除鄒琰(2010)和牛泓博等(2015)分別研究了山東和遼寧沿海菲律賓蛤仔形態(tài)差異[10-11]外，鮮有報(bào)道。本研究測(cè)定了采自不同海域的5個(gè)菲律賓蛤仔野生群體的形態(tài)參數(shù)，分析了其形態(tài)差異，可為今后的種質(zhì)資源保護(hù)和遺傳選育提供參考。

1 材料與方法

1.1 樣品采集

樣品于2018年9—12月，采自海州(HZ: 37°59′26.7″N，125°44′45.8″E)、南浦(NP: 39°27′49.5″N，125°12′40.2″E )、萊州(LZ: 37°21′24.2″N，119°50′18.9″E)、莆田(PT: 25°10′29.3″N，119°15′43.8″E)和汕尾(SW: 22°40′31.7″N，115°28′17.2″E)自然海區(qū)。通過低溫?；罘侄嗯芜\(yùn)至廈門小嶝水產(chǎn)科技有限公司實(shí)驗(yàn)基地，流水充氣暫養(yǎng)。樣本信息見表1。

表1 菲律賓蛤仔5個(gè)群體樣本信息Tab.1 Information of 5 populations of Manila clam

1.2 形態(tài)指標(biāo)及測(cè)定

每個(gè)群體測(cè)定50個(gè)樣本，外部形態(tài)測(cè)量指標(biāo)參考劉建勇等(2010)[12]，采用游標(biāo)卡尺(精確到0.01 mm)測(cè)量了殼長(zhǎng)(AB)、殼寬(GH)、殼高(OC)、韌帶長(zhǎng)(OD)、楯面長(zhǎng)(OE)、楯面寬(MN)、小月面長(zhǎng)(OF)、前端到腹緣(BC)、后端到腹緣(AC)9個(gè)形態(tài)性狀(圖1)，共獲取了2 250個(gè)數(shù)據(jù)。

圖1 菲律賓蛤仔形態(tài)測(cè)量位點(diǎn)Fig.1 Landmark points for morphological measurement of Manila clam

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

為消除不同群體規(guī)格大小對(duì)數(shù)據(jù)的影響，每個(gè)個(gè)體的數(shù)據(jù)除以殼長(zhǎng)進(jìn)行校正。應(yīng)用SPSS 17.0以校正后的數(shù)據(jù)對(duì)5個(gè)群體形態(tài)差異進(jìn)行單因素方差分析、聚類分析、主成分分析和判別分析。

1.3.1 單因素方差分析 采用Tukey s-b(K)法進(jìn)行多重比較，并劃分出同類子集，其顯著性水平設(shè)置為0.05。

1.3.2 聚類分析 利用各參數(shù)校正值的均值，采用平方歐氏距離法進(jìn)行系統(tǒng)聚類分析。

1.3.3 主成分分析 以校正后的8個(gè)參數(shù)值為變量，250個(gè)個(gè)體數(shù)據(jù)整體進(jìn)行主成分分析，選取特征值大于1的主成分進(jìn)行結(jié)果分析。

1.3.4 判別分析 采用逐步判別分析法，判別準(zhǔn)確率公式如下：

P1=n/N×100

(1)

P2=n/Nd×100

(2)

式(1～3)中：P1、P2為判別準(zhǔn)確率(%)，P為綜合判別率(%)，n為判別正確的菲律賓蛤仔數(shù)(個(gè))，N為實(shí)測(cè)菲律賓蛤仔數(shù)(個(gè))； Nd為判別菲律賓蛤仔數(shù)(個(gè))； Ai為第i個(gè)群體判別正確的菲律賓蛤仔數(shù)(個(gè))，Bi為第i個(gè)群體實(shí)際判別的菲律賓蛤仔數(shù)(個(gè))，k為群體數(shù)(個(gè))。

2 結(jié)果與討論

2.1 單因素方差分析

菲律賓蛤仔HZ組GH/AB、OC/AB、OD/AB和MN/AB顯著大于其他組(p<0.05，表2)，這表明HZ群體殼寬較厚，殼高較高，殼形接近球形，與其他組形態(tài)差異明顯。其他組之間各指標(biāo)互有差異，整體殼形差別不明顯。

表2 菲律賓蛤仔5個(gè)群體形態(tài)的方差分析Tab.2 ANOVA analysis of morphological parameters on 5 populations of Manila clam

注：同列數(shù)據(jù)上標(biāo)字母“a、b、c”完全不同表示有顯著差異(p<0.05)。

2.2 聚類分析

5個(gè)群體形態(tài)性狀的聚類分析(圖2)可以看出NP、LZ和PT等3組形態(tài)差異不大，而NP與LZ群體形態(tài)尤其接近；整體來看除HZ外的4組之間差異較小，而其整體與HZ差異很大，這說明HZ組形態(tài)與其它4組差異明顯。

2.3 主成分分析

按照特征值大于1，共提取2個(gè)主成分(表3)。主成分1解釋的總方差貢獻(xiàn)率為28.28%，主成分2為20.77%，累計(jì)貢獻(xiàn)率為49.05%，即2個(gè)主成分可解釋5個(gè)群體間形態(tài)差異的49.05%。在主成分1中GH/AB影響最大，其負(fù)荷值為0.802，其次為MN/AB，負(fù)荷值為0.751，這說明主成分1代表了“寬度”對(duì)殼型的影響；在主成分2中前端到BC/AB影響最大，其負(fù)荷值為0.794，其次為OF/AB，負(fù)荷值為0.618，這說明主成分2代表了“長(zhǎng)度”對(duì)殼型的影響。

圖2 菲律賓蛤仔5個(gè)群體形態(tài)性狀聚類分析Fig.2 Clustering dendrogram by morphological traits of 5 populations of Manila clam

表3 菲律賓蛤仔5個(gè)群體形態(tài)的2個(gè)主成分矩陣負(fù)荷值和貢獻(xiàn)率

Tab.3 Contribution ratios and loading of 2 principal components for morphological parameters of 5 populations of Manila clam

形態(tài)性狀主成分1主成分2GH/AB0.8020.239OC/AB0.520.406OD/AB0.4850.156OE/AB0.3720.275MN/AB0.751-0.012OF/AB0.050.618BC/AB-0.3340.794AC/AB0.591-0.572貢獻(xiàn)率28.28%20.77%

從圖3可見看出，除HZ群體外，其它4組混雜在一起，難以區(qū)分，這和聚類分析的結(jié)果較一致。

圖3 菲律賓蛤仔5個(gè)群體形態(tài)主成分散布圖Fig.3 Scattering diagram of the principal components for morphological parameters of 5 populations of Manila clam

2.4 判別分析

采用逐步判別分析，篩選出X1(GH/AB)、X2(OE/AB)、X3(MN/AB)和X4(AC/AB)共4個(gè)性狀特征值建立了5個(gè)群體的判別函數(shù)：

YHZ=814.485X1+220.653X2+331.141X3+

472.528X4-451.783

YNP=745.604X1+222.846X2+150.492X3+

448.072X4-386.301

YLZ=714.053X1+214.790X2+219.063X3+

438.213X4-367.209

YPT=755.939X1+235.811X2+179.968X3+

451.513X4-402.781

YSW=727.249X1+237.870X2+192.687X3+

466.209X4-402.165

表4為5個(gè)判別函數(shù)的各項(xiàng)系數(shù)和常數(shù)，表5為判別分析結(jié)果，判別準(zhǔn)確率P1和P2均在36.0%～80.0%之間；5個(gè)群體的綜合判別率為53.6%。HZ群體的判別率最高，為80.0%，這說明該群體特征明顯，易于判別。

表4 5個(gè)菲律賓蛤仔群體的Fisher線性判別式的各項(xiàng)系數(shù)及常數(shù)Tab.4 Coefficients and constants of Fisher discriminant functions for 5 populations of Manila clam

表5 5個(gè)菲律賓蛤仔群體判別分析結(jié)果Tab.5 Fisher discriminant results of 5 populations of Manila clam

2.5 討論

2.5.1 菲律賓蛤仔的形態(tài)差異 貝類形態(tài)變異受遺傳因子和環(huán)境因子共同影響，是對(duì)棲息環(huán)境長(zhǎng)期適應(yīng)的結(jié)果。本研究通過幾種多元分析方法對(duì)菲律賓蛤仔8個(gè)形態(tài)校正參數(shù)比較分析，結(jié)果表明HZ群體與NP、LZ、PT和SW群體在外部形態(tài)上差異顯著。雖然聚類分析將NP、LZ、PT和SW聚在一起，顯示出從北到南差異逐漸增大的趨勢(shì)，但是HZ與其它組的差異顯然不是地理隔離造成的，因?yàn)榫嚯x相近的NP與之差異巨大而與另外3個(gè)群體相似。其他研究中也存在類似的結(jié)果，牛泓博等認(rèn)為遼寧沿海菲律賓蛤仔不同群體間的形態(tài)變異程度與地理距離有一定關(guān)系，但主要受棲息地環(huán)境影響[11]。宋菲菲等(2012)認(rèn)為毛蚶(Scapharcasubcrenata)群體變異程度可能與地理距離有一定關(guān)系，同時(shí)受到棲息地環(huán)境的影響[13]。貝類均具有較強(qiáng)的表型環(huán)境可塑性，其表型變異并不總是由遺傳導(dǎo)致，更多是受棲息地環(huán)境影響[14-15]。Kwon等(1999)將韓國(guó)沿海外殼形態(tài)差異明顯的3組菲律賓蛤仔轉(zhuǎn)移到同一養(yǎng)殖區(qū)域，13個(gè)月后發(fā)現(xiàn)它們外形差異變得不顯著，證明棲息地環(huán)境是其外殼形態(tài)的主要影響因素[16]。Watanabe等(2010)對(duì)日本4個(gè)海灣菲律賓蛤仔形態(tài)學(xué)研究發(fā)現(xiàn)，棲息于貧營(yíng)養(yǎng)環(huán)境中的個(gè)體較營(yíng)養(yǎng)條件良好的個(gè)體殼寬更大[17]；Caill-Milly等(2012)研究表明，群體密度較低及褐色肌肉病會(huì)導(dǎo)致菲律賓蛤仔個(gè)體呈現(xiàn)更高的GH/AB比值，使外形接近于球形[18]。然而HZ組殼寬較大的形態(tài)差異是否由海州呈貧營(yíng)養(yǎng)化海區(qū)環(huán)境和較低種群密度造成，尚需對(duì)采樣地生態(tài)條件開展進(jìn)一步研究。

2.5.2 蛤類形態(tài)差異判別方法 傳統(tǒng)貝類分類鑒定依據(jù)是殼長(zhǎng)、殼高、殼寬、殼色和紋理等幾個(gè)指標(biāo)，主觀性較強(qiáng)，誤判率較高。因此基于形態(tài)測(cè)量學(xué)的聚類分析、主成分分析、判別分析等多元分析方法在蛤類不同地理群體的形態(tài)差異分析中廣泛應(yīng)用[19-21]，這些方法雖然基于測(cè)量數(shù)據(jù)，依靠數(shù)學(xué)統(tǒng)計(jì)分析得出區(qū)分結(jié)果，消除了人的主觀影響，但是由于測(cè)量的形態(tài)學(xué)指標(biāo)較少，在分析形態(tài)相似的地理群體間差異時(shí)仍是差強(qiáng)人意。也有學(xué)者應(yīng)用廣義線性模型和傳統(tǒng)形態(tài)測(cè)量學(xué)結(jié)合的方法對(duì)菲律賓蛤仔地理群體進(jìn)行判別分析[22]，但該方法篩選后采用的預(yù)測(cè)變量組合只有1種，沒有突破傳統(tǒng)形態(tài)測(cè)量學(xué)方法的短板，且由于要求響應(yīng)變量Y服從二項(xiàng)分布，導(dǎo)致研究群體只能為2個(gè)，更加局限了該方法的應(yīng)用。為消除形態(tài)學(xué)指標(biāo)少的缺陷，近年來基于橢圓傅里葉分析的輪廓線法[23]和基于空間坐標(biāo)系的地標(biāo)點(diǎn)法[24]等幾何形態(tài)測(cè)量學(xué)方法得到一定應(yīng)用，這兩種方法依靠計(jì)算機(jī)圖形處理技術(shù)采集更多更精確的形態(tài)幾何數(shù)據(jù)并做更準(zhǔn)確的定量分析。但幾何測(cè)量學(xué)方法需要專門的圖像采集硬件和專業(yè)的圖形處理軟件，成本較高，操作繁瑣，限制了它們的應(yīng)用。

在貝類地理群體差異分析中，形態(tài)差異分析是一個(gè)重要方面，然而大多數(shù)雙殼類尤其是蛤類外殼近似圓形，形態(tài)差異分析具有局限性，將其與生理生化、分子生物學(xué)指標(biāo)結(jié)合將會(huì)使結(jié)果更加準(zhǔn)確。

2.5.3 菲律賓蛤仔種質(zhì)資源保護(hù)和遺傳育種 我國(guó)是世界上進(jìn)行菲律賓蛤仔人工育苗規(guī)模最大、數(shù)量最多的國(guó)家。僅在福建省墾區(qū)育苗生產(chǎn)中，面積為86.6 hm2的池塘就可培育2 mm的稚貝604.9億粒[25]。大規(guī)模的育苗生產(chǎn)需要大規(guī)模的親貝供應(yīng)，但業(yè)者在選擇菲律賓蛤仔親貝時(shí)僅考慮性腺成熟度，而忽略其他性狀，使用一齡性腺飽滿的親貝育苗已是普遍現(xiàn)象。從育種學(xué)的角度看，這些菲律賓蛤仔親貝大多是未經(jīng)系統(tǒng)選育或改良的人工養(yǎng)殖群體，經(jīng)過多年反復(fù)養(yǎng)殖再育苗，存在近親繁殖的風(fēng)險(xiǎn)，商品蛤經(jīng)濟(jì)性狀日趨退化。因此，保護(hù)原始的優(yōu)良品種，并利用這些優(yōu)良的種質(zhì)資源培育新品種、新品系，已成為貝類養(yǎng)殖業(yè)亟待解決的重要課題[26]。本研究測(cè)量了多個(gè)形態(tài)性狀參數(shù)，通過形態(tài)性狀的多元分析評(píng)估了5個(gè)地理群體菲律賓蛤仔的種質(zhì)資源狀況，可為其種質(zhì)資源保護(hù)和遺傳育種工作提供依據(jù)。

3 結(jié)論

應(yīng)用多元分析方法對(duì)不同海域的5個(gè)菲律賓蛤仔野生群體形態(tài)數(shù)據(jù)進(jìn)行比較分析。結(jié)果表明HZ群體與其它4組形態(tài)差異顯著，GH/AB明顯較大，但這種差異是由遺傳差異還是環(huán)境差異造成，需進(jìn)一步研究。本研究可為菲律賓蛤仔種質(zhì)資源保護(hù)和遺傳育種提供基礎(chǔ)資料。