邱炬維　刁宏山　賴勝琪　尹聰　栗志民

摘要 利用形態(tài)特征分析、方差分析、主成分分析、判別分析和聚類分析5種多元統(tǒng)計(jì)分析方法，對(duì)北部灣沿岸東興（DX）、北海（BH）和防城港（FCG）3個(gè)野生鈍綴錦蛤（Tapes dorsatus）群體進(jìn)行形態(tài)學(xué)比較。結(jié)果表明：東興群體體型最大，貝殼較為厚重；北海群體軟體部最為飽滿；防城港群體殼型“凸”形明顯，貝殼隆起程度最高。通過主成分分析共獲得3個(gè)主成分，主成分1、主成分2、主成分3的貢獻(xiàn)率分別為35.369%、21.282%和14.968%，累計(jì)貢獻(xiàn)率為71.619%。通過逐步判別法選出8個(gè)貢獻(xiàn)率較大的性狀進(jìn)行判別分析，建立了3個(gè)種群的判別函數(shù)，判別準(zhǔn)確率P1為50.00％～90.00％，P2為60.98%～78.95%，3個(gè)種群的綜合判別率為70.37％。聚類分析結(jié)果顯示，北海群體與防城港群體的鈍綴錦蛤具有更為接近的形態(tài)特征，而東興群體由于殼型較扁，在形態(tài)上與另外2個(gè)群體存在明顯差異。北部灣地區(qū)鈍綴錦蛤不同群體間的形態(tài)學(xué)差異是種間雜交和環(huán)境共同作用的結(jié)果。該研究結(jié)果可為北部灣鈍綴錦蛤群體的人工養(yǎng)殖與種質(zhì)資源保護(hù)提供依據(jù)。

關(guān)鍵詞 鈍綴錦蛤；形態(tài)學(xué)差異；主成分分析；判別分析；聚類分析

中圖分類號(hào) S917.4 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2023）06-0088-05

doi：10.3969/j.issn.0517-6611.2023.06.022

Study on the Morphological Difference of Different Geographical Populations of Tapes dorsatus along the Coast of Beibu Gulf

QIU Ju－wei， DIAO Hong－shan， LAI Sheng－qi et al

（Fisheries College， Guangdong Ocean University， Zhanjiang，Guangdong 524000）

Abstract By using five multivariate statistical analysis methods（morphological characteristics analysis， variance analysis， principal component analysis， discriminant analysis and cluster analysis）， the morphological traits among three wild populations of Tapes dorsatus along the coast of Beibu Gulf， including Dongxing （DX）， Beihai （BH）， Fangchenggang （FCG）， were compared. The results showed that Dongxing population had the largest body type and thick shell. The soft part of Beihai population was the most full. Fangchenggang population had obvious convex shell type and? the highest shell uplift. Three principal components were obtained? by the principal component analysis method， the contributory rates of principal component 1，2，3 were 35.369%， 21.282% and 14.968%， respectively， the cumulative contributory rate was 71.619%. The stepwise discriminant method was used to select 8 characters with higher contribution rate for discriminant analysis. The discrimination functions of the three populations were set up， and the identification accuracy P1 and P2 were 50.00%-90.00% and 60.98%-78.95% respectively. The total discriminate accuracy was 70.37%.? The results of cluster analysis indicated that the morphological characters were similar between the populations of BH and FCG. Dongxing population was different from the other two populations in morphology because of its flat shell type. The morphological variations of T.dorsatus in Beibu Gulf region among different populations were resulted from the joint effects of interspecific cross and environment. The research results could provide basis for the conservation of germ plasma resources and artificial breeding of T.dorsatus in Beibu Gulf region.

Key words Tapes dorsatus;Morphological difference;Principal component analysis;Discrimination analysis;Cluster analysis

鈍綴錦蛤（Tapes dorsatus）隸屬軟體動(dòng)物門（Mollusca）瓣鰓綱（Bivalvia）簾蛤目（Veneroida）簾蛤科（Veneridae）綴錦蛤?qū)伲═apes），是一種熱帶、亞熱帶海洋貝類，廣泛分布于西南太平洋地區(qū)。在日本、菲律賓、印度、印度尼西亞、澳大利亞、新西蘭等國(guó)家的沿海淺灘和泥灘地帶均發(fā)現(xiàn)野生鈍綴錦蛤群體［1］，而在我國(guó)則常見于東海、南海、北部灣地區(qū)與臺(tái)灣島海域。鈍綴錦蛤生活在潮間帶中、低潮區(qū)至潮下帶沙泥和泥沙底質(zhì)中，營(yíng)埋棲生活，一般棲息深度為3～15 cm，喜碎石泥沙環(huán)境，在潮間帶下部珊瑚塊處也能找到野生集群［2-3］。鈍綴錦蛤個(gè)體較大，具有生長(zhǎng)速度快、適應(yīng)性廣、抗污能力強(qiáng)等優(yōu)點(diǎn)，是一種品質(zhì)優(yōu)良的可食用貝類。目前鈍綴錦蛤已經(jīng)發(fā)展成為一種新的海水貝類養(yǎng)殖品種，具有較高的經(jīng)濟(jì)價(jià)值及較好的養(yǎng)殖發(fā)展前景［4］。近年來(lái)，國(guó)外有關(guān)鈍綴錦蛤的研究報(bào)道較少。Nell等［5］報(bào)道了鈍綴錦蛤二倍體和三倍體的人工孵化，認(rèn)為它是一種具有養(yǎng)殖潛力的新品種。國(guó)內(nèi)科研人員從多方面對(duì)該種貝類開展了研究，包括海洋生物源肝素開發(fā)［6］、形態(tài)分類［7］、生態(tài)習(xí)性［2］、分子學(xué)分類［8］、繁殖發(fā)育［4］以及人工育苗［9］和育種［10］等。

目前，鈍綴錦蛤的人工繁育和養(yǎng)殖受到科研人員的廣泛關(guān)注。2019年，廣東海洋大學(xué)劉志剛教授團(tuán)隊(duì)參與了廣西創(chuàng)新驅(qū)動(dòng)發(fā)展專項(xiàng)，開展了鈍綴錦蛤的人工繁育，并開展了鈍綴錦蛤良種選育，2021年11—12月培育苗種（殼長(zhǎng)2～3 mm）約1.9億粒。然而，鈍綴錦蛤的人工繁育和養(yǎng)殖仍處于試驗(yàn)階段，培育的人工苗種還很少。隨著市場(chǎng)對(duì)鈍綴錦蛤需求量的日益增大，市場(chǎng)價(jià)格不斷攀升。人工捕撈天然鈍綴錦蛤不僅無(wú)法滿足市場(chǎng)的需求，而且酷漁濫捕造成自然資源量逐年銳減。為了保護(hù)北部灣海域鈍綴錦蛤群體，對(duì)北部灣海域該貝類種質(zhì)資源現(xiàn)狀進(jìn)行評(píng)估十分必要。貝類不同地理群體形態(tài)學(xué)差異分析可作為評(píng)估貝類種質(zhì)資源的一種重要手段。目前，科研人員已對(duì)多種貝類開展了相關(guān)研究，例如毛蚶［11］、泥蚶［12］、菲律賓蛤仔［13］、青蛤［14］、真曲巴非蛤［15］等。筆者利用5種多元統(tǒng)計(jì)分析方法比較北部灣海域鈍綴錦蛤3個(gè)不同地理群體的形態(tài)學(xué)差異，并對(duì)該貝類的種質(zhì)資源進(jìn)行評(píng)估，旨在為保護(hù)北部灣地區(qū)野生鈍綴錦蛤群體資源、研究該貝類的遺傳育種技術(shù)提供科學(xué)依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

2020年8月，選取北部灣海區(qū)3個(gè)不同地理群體野生鈍綴錦蛤2齡貝，采捕海區(qū)分別為廣西東興（107°97′E，21°54′N）、廣西北海（109°47′E，21°55′N）、廣西防城港（108°20′E，21°37′N）。

1.2 試驗(yàn)方法

1.2.1 形態(tài)學(xué)性狀測(cè)量。使用游標(biāo)卡尺對(duì)殼長(zhǎng)（LS）、殼寬（WS）、殼高（HS）、前緣長(zhǎng)（LE）、后緣長(zhǎng)（LT）、外韌帶長(zhǎng)（LY）、外韌帶槽寬（WY）、殼厚（TS）共8個(gè)殼形態(tài)學(xué)性狀（圖1）進(jìn)行測(cè)量，精確到0.01 mm。總質(zhì)量（WT）、軟體質(zhì)量（PW）、殼質(zhì)量（SW）使用電子天平稱量，精確到0.01 g。該研究使用的是各個(gè)性狀的比值，以消除鈍綴錦蛤個(gè)體規(guī)格不同所帶來(lái)的影響［16］。比例性狀包括殼高/殼長(zhǎng)（HS/LS）、殼寬/殼長(zhǎng)（WS/LS）、殼寬/殼高（WS/HS）、前緣長(zhǎng)/后緣長(zhǎng)（LE/LT）、外韌帶槽寬/外韌帶長(zhǎng)（WY/LY）、殼厚/殼長(zhǎng)（TS/LS）、總質(zhì)量/殼長(zhǎng)（WT/LS）、軟體質(zhì)量/殼長(zhǎng)（PW/LS）、殼質(zhì)量/殼長(zhǎng)（SW/LS）9個(gè)指標(biāo)。

利用Excel 2007軟件分析3個(gè)群體鈍綴錦蛤的形態(tài)學(xué)性狀，統(tǒng)計(jì)各比例性狀的平均值（）、標(biāo)準(zhǔn)差（SD）和變異系數(shù)（CV）。

1.2.2 方差分析。

采用SPSS 20.0統(tǒng)計(jì)軟件對(duì)鈍綴錦蛤不同群體間形態(tài)學(xué)性狀進(jìn)行單因素方差（One-way ANOVA）分析。根據(jù)P值來(lái)判斷9個(gè)比例性狀在不同群體間是否存在顯著差異。P＜0.01表示差異極顯著，P＜0.05表示差異顯著，P＞0.05表示無(wú)顯著差異。

1.2.3 主成分分析。

使用SPSS 20.0統(tǒng)計(jì)軟件對(duì)3個(gè)群體鈍

綴錦蛤的形態(tài)學(xué)性狀進(jìn)行主成分分析，從9個(gè)比例性狀中選出3個(gè)綜合性指標(biāo)，即互不關(guān)聯(lián)的3個(gè)主成分，并分別計(jì)算各主成分的貢獻(xiàn)率。根據(jù)得到的3個(gè)主成分貢獻(xiàn)率，按從大到小的順序進(jìn)行排序，繪制主成分散布圖。

1.2.4 判別分析。

對(duì)所有樣本進(jìn)行逐步判別分析。判別準(zhǔn)確率的計(jì)算公式如下：

式中，實(shí)測(cè)樣本數(shù)是指該種群的采樣個(gè)體數(shù)，判別樣本數(shù)是指被判為該種群的個(gè)體數(shù)量；Ai為第i個(gè)種群中判別正確的個(gè)體數(shù)，Bi為第i個(gè)種群的樣本數(shù)量，k為種群數(shù)。

1.2.5 聚類分析。

利用分層聚類法對(duì)3個(gè)種群的9個(gè)比例性狀進(jìn)行分析，采用NCSS 2007統(tǒng)計(jì)軟件求出各個(gè)比例性狀的平均值，并進(jìn)行系統(tǒng)聚類，最后用樹狀圖顯示3個(gè)群體間的親疏程度。

2 結(jié)果與分析

2.1 各群體形態(tài)比例性狀統(tǒng)計(jì)

3個(gè)不同群體9個(gè)比例性狀的統(tǒng)計(jì)結(jié)果見表1。從表1可以看出，，3個(gè)群體鈍綴錦蛤的HS/LS平均值均小于1，其中防城港群體（FCG）鈍綴錦蛤的HS/LS平均值最大（0.657 3），說(shuō)明該群體鈍綴錦蛤殼型較圓；東興群體（DX）鈍綴錦蛤HS/LS平均值最小，為0.652 6，說(shuō)明該群體鈍綴錦蛤殼型較扁。從WS/LS來(lái)看，防城港群體（FCG）WS/LS平均值最大（0.425 9），說(shuō)明該群體鈍綴錦蛤貝殼隆起程度最高；北海群體（BH）WS/LS平均值最?。?.410 7），說(shuō)明該群體鈍綴錦蛤貝殼隆起程度最低。從TS/LS來(lái)看，防城港群體（FCG）TS/LS平均值最大，為0.304 2；北海群體（BH）TS/LS平均值最小，僅為0.254 1。防城港群體（FCG）除PW/LS、SW/LS以外的其他指標(biāo)均最大，而北海群體這2個(gè)指標(biāo)最大。在3個(gè)群體中，TS/LS、LE/LT以及WY/LY的變異系數(shù)均超過0.100 0。根據(jù)統(tǒng)計(jì)學(xué)原理［17］可知，這3個(gè)性狀在各群體間表現(xiàn)出廣泛的變異。從SW/LS來(lái)看，東興群體與北海群體鈍綴錦蛤SW/LS的變異系數(shù)均超過0.100 0，表明東興群體與北海群體這一性狀出現(xiàn)了變異。

2.2 方差分析

3個(gè)群體鈍綴錦蛤9個(gè)比例性狀的單因素方差分析結(jié)果如表2所示。從表2可以看出，不同群體鈍綴錦蛤除HS/LS、TS/LS、WY/LY外，其他比例性狀均在不同群體間差異極顯著（P＜0.01）。該研究結(jié)果揭示出不同性狀在3個(gè)野生群體間存在一定的差異，為不同群體鈍綴錦蛤的鑒別提供了可能。

2.3 主成分分析

為明確鈍綴錦蛤群體形態(tài)學(xué)差異的主成分因子，利用SPSS 20.0統(tǒng)計(jì)軟件進(jìn)行綜合運(yùn)算，以經(jīng)過校正的形態(tài)學(xué)指標(biāo)作為基礎(chǔ)數(shù)據(jù)進(jìn)行分析，得到了3個(gè)互不關(guān)聯(lián)的主成分，結(jié)果如表3所示。

由表3可知，主成分1的貢獻(xiàn)率為35.369%，主成分2的貢獻(xiàn)率為21.282%，主成分3的貢獻(xiàn)率為14.968%，累計(jì)貢獻(xiàn)率為71.619%。在主成分1中，SW/LS的影響最大，其貢獻(xiàn)率為81.2%；其次為WT/LS，其貢獻(xiàn)率為78.6%。在主成分2中，TS/LS的影響最大，其貢獻(xiàn)率為52.2%；其次為WS/LS，其貢獻(xiàn)率為49.5%。在主成分3中，HS/LS的影響最大，其貢獻(xiàn)率為48.4%；其次為L(zhǎng)E/LT，其貢獻(xiàn)率為46.7%。

根據(jù)表3中各比例性狀與主成分載荷矩陣，計(jì)算出每個(gè)樣品的主成分值，以3個(gè)主成分大小的中間值（即主成分2）作為基本關(guān)聯(lián)量，分別以主成分1和主成分3作為散布圖的橫坐標(biāo)軸，以主成分2作為散布圖的縱坐標(biāo)軸，繪制鈍綴錦蛤3個(gè)地理群體的散布圖，如圖2所示。圖2中主成分1、主成分2、主成分3均為相對(duì)值，無(wú)量綱。在圖2A中，從主成分1和主成分2的關(guān)聯(lián)性可以看出，北海群體與防城港群體明顯交疊在一起，說(shuō)明這2個(gè)群體主成分1和主成分2的關(guān)聯(lián)性具有較高的相似性；東興群體與另外2個(gè)群體間無(wú)重疊區(qū)域，且分布較為零散，說(shuō)明東興群體與其他2個(gè)群體間的差異較大。在圖2B中，3個(gè)群體之間無(wú)重疊區(qū)域，但在主成分2軸上北海群體與防城港群體具有較高的相似度，而東興群體則明顯不同；在主成分3軸上，3個(gè)群體分成3組，證實(shí)了該項(xiàng)指標(biāo)在不同群體間存在一定的差異。

2.4 判別分析

采用逐步判別分析對(duì)矯正后的鈍綴錦蛤性狀特征值進(jìn)行篩選，分別以X1、X2、X3、X4、X5、X6、X7、X8代表殼高/殼長(zhǎng)（HS/LS）、殼寬/殼長(zhǎng)（WS/LS）、殼厚/殼長(zhǎng)（TS/LS）、總質(zhì)量/殼長(zhǎng)（WT/LS）、軟體質(zhì)量/殼長(zhǎng)（PW/LS）、殼質(zhì)量/殼長(zhǎng)（SW/LS）、前緣長(zhǎng)/后緣長(zhǎng)（LE/LT）、外韌帶槽寬/外韌帶長(zhǎng)（WY/LY），建立3個(gè)鈍綴錦蛤群體的判別函數(shù)

，得到3個(gè)群體鈍綴錦蛤形態(tài)比例性狀判別函數(shù)的各項(xiàng)系數(shù)及常數(shù)項(xiàng)，如表4所示。

3個(gè)群體的判別函數(shù)如下：

YDX=1 595.532X1+393.058X2-3 121.808X3+68.007X4-56.243X5-71.163X6+110.875X7-85.193X8-603.123

YBH=1 629.373X1+351.312X2-2 796.137X3+70.943X4-86.095X5-126.497X6+131.571X7-57.262X8-611.873

YFCG=1 591.711X1+421.394X2-2 960.636X3+53.773X4-19.729X5-64.730X6+102.767X7-105.057X8-610.642

判別分析結(jié)果如表5所示，3個(gè)群體鈍綴錦蛤判別準(zhǔn)確率P1為50.00%～90.00％，判別分析準(zhǔn)確率P2為60.98%～78.95%，綜合判別率為70.37%。

2.5 聚類分析

3個(gè)群體鈍綴錦蛤的聚類分析表明，北海與防城港群體的形態(tài)學(xué)差異不大，較為接近；東興群體的趨異程度最大，與其他2個(gè)群體的形態(tài)學(xué)差異最大（圖3）。

3 討論

3.1 不同群體的形態(tài)學(xué)差異分析

3個(gè)鈍綴錦蛤群體中，防城港群體的殼高/殼長(zhǎng)（HS/LS）平均值最大，殼型較圓；東興群體HS/LS平均值最小，殼型較扁。從殼寬/殼長(zhǎng)（WS/LS）來(lái)看，防城港群體WS/LS平均值最大，說(shuō)明該群體鈍綴錦蛤貝殼隆起程度最高；北海群體WS/LS平均值最小，說(shuō)明該群體鈍綴錦蛤貝殼隆起程度最低。在3個(gè)鈍綴錦蛤群體中，殼厚/殼長(zhǎng)、前緣長(zhǎng)/后緣長(zhǎng)以及外韌帶槽寬/外韌帶長(zhǎng)的變異系數(shù)均超過0.100 0，說(shuō)明這3個(gè)性狀在同一群體內(nèi)不同個(gè)體間具有廣泛的變異。從殼質(zhì)量/殼長(zhǎng)來(lái)看，東興與北海群體SW/LS的變異系數(shù)均超過了0.100 0，表明東興與北海群體這一性狀出現(xiàn)了變異。牛泓博等［13］對(duì)遼寧沿海菲律賓蛤仔（Ruditapes philippinarum）8 個(gè)野生群體進(jìn)行形態(tài)學(xué)比較，結(jié)果發(fā)現(xiàn)由于各海域生態(tài)環(huán)境存在一定差異，不同海域菲律賓蛤仔群體形態(tài)表現(xiàn)出一定的差異。劉童等［18］以韓國(guó)釜山、中國(guó)榮成、中國(guó)即墨、中國(guó)長(zhǎng)島４個(gè)地理群體魁蚶為研究對(duì)象，結(jié)果發(fā)現(xiàn)4個(gè)群體魁蚶形態(tài)存在一定差異，筆者認(rèn)為魁蚶形態(tài)不同可能是對(duì)不同環(huán)境適應(yīng)的結(jié)果。環(huán)境底質(zhì)是決定形態(tài)變異的主要因素之一，這種表型特征的地理變異也是長(zhǎng)期演變的結(jié)果。該研究中北部灣海域鈍綴錦蛤3個(gè)野生群體（東興、北海、防城港）在形態(tài)上存在一定的差異，可能與北部灣海域3個(gè)群體所在的海域環(huán)境有關(guān)。

主成分分析是一種常用的簡(jiǎn)化數(shù)據(jù)集方法，能夠?qū)⒃械亩鄠€(gè)指標(biāo)轉(zhuǎn)化為少數(shù)幾個(gè)代表性較好的綜合性指標(biāo)，且保持各個(gè)指標(biāo)之間相對(duì)獨(dú)立，以此來(lái)反映原來(lái)指標(biāo)大部分的信息。判別分析是根據(jù)事物各特征的變量值和它所屬類別得出判別函數(shù)，以此為依據(jù)對(duì)研究對(duì)象進(jìn)行所屬類判別的方法，這2種方法已在貝類形態(tài)學(xué)差異分析中廣泛應(yīng)用［19］。古恒光等［14］分析了4個(gè)野生青蛤群體的7個(gè)比例性狀，得到3個(gè)主成分，累計(jì)貢獻(xiàn)率為63.428%；將7個(gè)比例性狀全部引入判別公式，綜合判別率為64.3%。劉達(dá)博等［20］對(duì)我國(guó)沿海7個(gè)地理群體縊蟶的12個(gè)形態(tài)學(xué)性狀進(jìn)行了研究，得到3個(gè)主成分，累計(jì)貢獻(xiàn)率為71.23%；建立7個(gè)群體縊蟶的判別函數(shù)，綜合判別率為77.93%。紀(jì)燕如等［21］對(duì)我國(guó)東南沿海5個(gè)織錦巴非蛤群體11個(gè)形態(tài)學(xué)性狀進(jìn)行比較研究，通過主成分分析構(gòu)建了3個(gè)主成分，累計(jì)貢獻(xiàn)率為66.005％；建立了5個(gè)群體的判別函數(shù)，綜合判別率為52.60％。該研究結(jié)果表明，通過主成分分析共獲得3個(gè)主成分，累計(jì)貢獻(xiàn)率為71.619%；利用逐步判別法選擇出8個(gè)貢獻(xiàn)率較大的性狀進(jìn)行判別分析，建立了3個(gè)種群的判別函數(shù)，綜合判別率為70.37％。

3.2 環(huán)境條件對(duì)不同鈍綴錦蛤群體形態(tài)的影響

在貝類生長(zhǎng)發(fā)育過程中，不同的地理環(huán)境會(huì)顯著影響貝類的形態(tài)學(xué)性狀，造成群體與群體之間、個(gè)體與個(gè)體之間都會(huì)存在一定的差異。這些影響因素主要包括溫度、光照、鹽度、餌料等外界條件。大多數(shù)學(xué)者認(rèn)為，貝類的形態(tài)變異是一種生態(tài)適應(yīng)現(xiàn)象，是對(duì)棲息地生活環(huán)境長(zhǎng)期適應(yīng)的表現(xiàn)，是遺傳因子與環(huán)境因子共同作用的結(jié)果。馮建彬等［22］對(duì)我國(guó)4個(gè)海區(qū)的野生文蛤群體進(jìn)行形態(tài)學(xué)分析，結(jié)果發(fā)現(xiàn)造成文蛤不同群體間形態(tài)學(xué)性狀差異的主要因素為溫度、鹽度和底質(zhì)營(yíng)養(yǎng)。宋菲菲等［11］認(rèn)為地理距離限制是造成毛蚶群體變異程度不同的主要原因。Laudien等［23］認(rèn)為浪蛤的形態(tài)學(xué)差異是由地理環(huán)境差異決定的，而不同群體間的遺傳相似性則是由基因交流、選擇壓力和進(jìn)化造成的。該研究通過形態(tài)特征分析、方差分析、主成分分析、判別分析與聚類分析發(fā)現(xiàn)，鈍綴錦蛤3個(gè)群體在不同地理環(huán)境下出現(xiàn)了形態(tài)學(xué)變異。北部灣沿海海岸線曲折，受東南亞熱帶季風(fēng)的影響，容易產(chǎn)生小氣候區(qū)［24］，形成獨(dú)特的生態(tài)環(huán)境。鈍綴錦蛤活動(dòng)能力較差，無(wú)法進(jìn)行長(zhǎng)距離快速移動(dòng)，極易受到棲息地環(huán)境的影響。研究表明，東興市6—8月平均氣溫約29.2 ℃［25］，而防城港市與北海市氣溫相差不大，6—8月平均氣溫約30.1 ℃［26］。鹽度調(diào)查表明，防城港海域平均海水鹽度為28.0［27］，北海海域平均海水鹽度為29.28，均高于東興海域平均海水鹽度［28］。陳宜展等［29］對(duì)北部灣北部季節(jié)環(huán)流的研究表明，防城港以東至北海的外海海域呈現(xiàn)出夏季順時(shí)針環(huán)流、冬季逆時(shí)針環(huán)流的特殊環(huán)境使得兩地海域環(huán)境出現(xiàn)極大的相似性。東興市地處中國(guó)與越南邊境接壤處，受全年逆時(shí)針環(huán)流的影響，使得該地區(qū)氣溫與海水鹽度更容易受到東南亞海洋季風(fēng)的影響，呈現(xiàn)出不同的生態(tài)環(huán)境。這些環(huán)境差異很可能影響鈍綴錦蛤的生長(zhǎng)發(fā)育速度，從而反映為外部形態(tài)特征的差異。

3.3 鈍綴錦蛤種質(zhì)資源保護(hù)與利用

近年來(lái)，如何保護(hù)和利用好貝類種質(zhì)資源、如何確保我國(guó)養(yǎng)殖海水貝類產(chǎn)量穩(wěn)步增加等問題層出不窮，掀起了我國(guó)水產(chǎn)動(dòng)物研究和海水貝類養(yǎng)殖的新熱潮。種質(zhì)資源的完好是一個(gè)物種能夠延續(xù)至今的重要法寶，也是所有水產(chǎn)科學(xué)研究前沿工作者的重點(diǎn)研究?jī)?nèi)容。每一個(gè)物種都具有獨(dú)特的基因多樣性［30］，因此物種是組成大自然基因多樣性的基礎(chǔ)，而自然生物資源評(píng)價(jià)則很大程度取決于大自然的基因多樣性［31］。由于多年來(lái)人們對(duì)北部灣海域鈍綴錦蛤的酷漁濫捕，該種貝類的野生群體已經(jīng)遭到嚴(yán)重破壞，自然資源量銳減，因此開展該貝類種質(zhì)資源保護(hù)與利用就顯得十分必要。種質(zhì)資源的研究對(duì)于了解貝類種源的適應(yīng)性、起源、基因資源分布、基因資源保護(hù)以及種群改良等方面都具有重要價(jià)值，已經(jīng)越來(lái)越為廣大育種和養(yǎng)殖人員所重視。鈍綴錦蛤形態(tài)學(xué)分析可為揭示其種群遺傳學(xué)特征提供重要依據(jù)。鈍綴錦蛤在我國(guó)北部灣海域廣泛分布，近年來(lái)該種貝類已經(jīng)成為這一海域重要的養(yǎng)殖品種，目前已經(jīng)有少量苗種供應(yīng)養(yǎng)殖，然而隨著養(yǎng)殖規(guī)模的擴(kuò)大，當(dāng)?shù)孛绶N供給不夠充足，無(wú)法滿足養(yǎng)殖的需要，因此不得不將來(lái)自不同海域的鈍綴錦蛤苗種引進(jìn)北部灣海域進(jìn)行人工養(yǎng)殖。由于不同地理種群的交流，容易造成雜交污染，使北部灣海域土著鈍綴錦蛤種群遺傳多樣性降低。該研究中筆者采用各個(gè)性狀的比值作為分析指標(biāo)，很好地消除了不同個(gè)體間因規(guī)格不同所帶來(lái)的誤差。

采用方差分析、主成分分析、判別分析和聚類分析等多種統(tǒng)計(jì)方法，分析各群體鈍綴錦蛤的形態(tài)學(xué)指標(biāo)，從不同角度分析了鈍綴錦蛤不同群體間的形態(tài)學(xué)差異。主成分分析和聚類分析顯示，防城港和北海群體相似性較大，可能與地理距離較近，這兩個(gè)群體的基因交流頻繁有關(guān)。因此，只用貝類的形態(tài)學(xué)指標(biāo)不能完全揭示各群體的特性。在今后的研究中應(yīng)結(jié)合形態(tài)學(xué)分析開展分子遺傳學(xué)的相關(guān)研究，有利于全面證實(shí)該種貝類不同群體間的遺傳關(guān)系。

參考文獻(xiàn)

［1］ 蔡英亞，莊啟謙.關(guān)于拉漢軟體動(dòng)物名稱一些問題的商榷［J］.湛江海洋大學(xué)學(xué)報(bào)，2001，21（4）：84.

［2］ 黃洋，杜濤，楊世平.鈍綴錦蛤生態(tài)習(xí)性的初步研究［J］.水產(chǎn)科學(xué)，2008，27（4）：175-178.

［3］ 巫旗生，曾志南，寧岳，等.鈍綴錦蛤形態(tài)性狀對(duì)活體質(zhì)量的影響［J］.水產(chǎn)科學(xué)，2018，37（1）：110-114.

［4］ 巫旗生，文宇，曾志南，等.鈍綴錦蛤繁殖周期和胚胎發(fā)育［J］.中國(guó)水產(chǎn)科學(xué)，2017，24（3）：488-496.

［5］ NELL J A，O’CONNOR W A，HAND R E，et al.Hatchery production of diploid and triploid clams，Tapes dorsatus（Lamarck 1818）：A potential new species for aquaculture［J］.Aquaculture，1995，130（4）：389-394.

［6］ 杜振興，周斯儀，鐘賽意，等.不同海洋生物源肝素的理化性質(zhì)及抗凝血活性［J］.食品科學(xué)，2019，40（17）：134-140.

［7］ 蔡英亞，謝紹河.廣東的海貝［M］.修訂版.汕頭：汕頭大學(xué)出版社，2006：336.

［8］ 陳道海，李洪英，吳秋穎，等.10種簾蛤科貝類COI基因序列分析及系統(tǒng)發(fā)育研究［J］.生物資源，2018，40（3）：277-284.

［9］ 彭銀輝.鈍綴錦蛤的人工育苗方法［R］.廣西壯族自治區(qū)，廣西壯族自治區(qū)海洋研究所，2014-08-06.

［10］ 聶振平，彭慧婧，鄒杰，等.鈍綴錦蛤選育群體F2生長(zhǎng)性狀相關(guān)性及遺傳力分析［J］.廣西科學(xué)，2020，27（3）：241-247.

［11］ 宋菲菲，王其翔，劉名，等.不同地區(qū)毛蚶群體的形態(tài)差異的比較研究［J］.中國(guó)海洋大學(xué)學(xué)報(bào)（自然科學(xué)版），2012，42（11）：40-45.

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

［13］ 牛泓博，聶鴻濤，趙力強(qiáng)，等.遼寧沿海菲律賓蛤仔不同地理群體形態(tài)差異研究［J］.海洋科學(xué)，2015，39（11）：54-60.

［14］ 古恒光，劉建勇，楊國(guó)宏.我國(guó)南海4個(gè)青蛤野生群體的形態(tài)差異分析［J］.廣東海洋大學(xué)學(xué)報(bào)，2010，30（3）：1-5.

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

［16］ 錢榮華，李家樂，董志國(guó)，等.中國(guó)五大湖三角帆蚌形態(tài)差異分析［J］.海洋與湖沼，2003，34（4）：436-443.

［17］ 蔡一林，岳永生.水產(chǎn)生物統(tǒng)計(jì)［M］.北京：中國(guó)農(nóng)業(yè)出版社，2004：11-21.

［18］ 劉童，鄒琰，王英俊，等.魁蚶4個(gè)群體形態(tài)差異比較分析［J］.水產(chǎn)科學(xué)，2019，38（5）：676-681.

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

［20］ 劉達(dá)博，牛東紅，姜志勇，等.縊蟶7群體的形態(tài)差異與判別分析［J］.海洋漁業(yè)，2009，31（4）：363-368.

［21］ 紀(jì)燕如，林志華，伍榮聰，等.我國(guó)東南沿海5個(gè)織錦巴非蛤地理群體的形態(tài)差異分析［J］.熱帶生物學(xué)報(bào)，2011，2（3）：219-225.

［22］ 馮建彬，李家樂，王美珍，等.我國(guó)四海區(qū)不同群體文蛤形態(tài)差異與判別分析［J］.浙江海洋學(xué)院學(xué)報(bào)（自然科學(xué)版），2005，24（4）：318-323.

［23］ LAUDIEN J，F(xiàn)LINT N S，VAN DER BANK F H，et al.Genetic and morphological variation in four populations of the surf clam Donax serra（Rding）from southern African sandy beaches［J］.Biochemical systematics and ecology，2003，31（7）：751-772.

［24］ 符華兒，尹業(yè)民，羅海鵬，等.季風(fēng)對(duì)北部灣北岸海平面的影響［J］.廣西科學(xué)院學(xué)報(bào)，1983（1）：68-71.

［25］ 夏華永，古萬(wàn)才.廣西沿海海洋站觀測(cè)海水溫度的統(tǒng)計(jì)分析［J］.海洋通報(bào)，2000，19（4）：15-21.

［26］ 黎彥佐，李妍君.防城港市氣溫與降水南北氣候差異分析［J］.氣象研究與應(yīng)用，2018，39（3）：30-34.

［27］ 曹振軼，鮑敏，管衛(wèi)兵，等.北部灣東北部水團(tuán)分布及季節(jié)變化分析［J］.海洋與湖沼，2019，50（3）：532-542.

［28］ 楊芳.三個(gè)細(xì)基江蘺繁枝變種種群的生長(zhǎng)和氮營(yíng)養(yǎng)特征研究［D］.南寧：廣西大學(xué)，2014：1-2.

［29］ 陳宜展，楊威，曹永港，等.北部灣北部環(huán)流的季節(jié)特征［J］.廣東海洋大學(xué)學(xué)報(bào)，2020，40（4）：68-74.

［30］ PRIMACK R B.Tropical community dynamics and conservation biology［J］.BioScience，1992，42（11）：818-821.

［31］ 孫振興，常林瑞.貝類種質(zhì)資源保護(hù)研究進(jìn)展［J］.海洋湖沼通報(bào)，2005（3）：103-108.

基金項(xiàng)目 2019年廣西創(chuàng)新驅(qū)動(dòng)發(fā)展專項(xiàng)（桂科AA19254032-3）。

作者簡(jiǎn)介 邱炬維（2000—），男，云南曲靖人，研究方向：水產(chǎn)養(yǎng)殖。*通信作者，教授，從事貝類遺傳育種研究。

收稿日期 2022-04-01