雷寧 曾建剛 沈小明 葉金云 朱俊杰

摘要 運(yùn)用3種多元統(tǒng)計(jì)分析方法對(duì)采自太湖、錢塘江桐廬段的野生花鰩骨以及由太湖野生群體經(jīng)人工繁育3代的花鰩骨這三個(gè)群體的30個(gè)形態(tài)參數(shù)進(jìn)行了比較。結(jié)果顯示，太湖群體和東林養(yǎng)殖群體的形態(tài)較為接近，而與桐廬群體的差異較大。利用主成分分析構(gòu)建了3個(gè)主成分，其貢獻(xiàn)率分別為主成分1 55.08%、主成分2 14.01%、主成分3 6.87%，它們的累積貢獻(xiàn)率為75.97%。逐步判別選入3個(gè)貢獻(xiàn)率較大的性狀進(jìn)行判別分析，通過建立的3個(gè)群體判別函數(shù)，其判別準(zhǔn)確率P1為75%～100%，P2 為77%～100%，綜合判別率為86%，結(jié)果表明3個(gè)群體形態(tài)差異顯著（P<0.01）。該研究為花鰩骨不同群體的鑒別積累了基礎(chǔ)資料以及為花鰩骨的增殖放流、優(yōu)良品種的培育和種質(zhì)資源保護(hù)提供理論依據(jù)。

關(guān)鍵詞 花鰩骨;形態(tài)特征;聚類分析;主成分分析;判別分析

中圖分類號(hào) Q954? 文獻(xiàn)標(biāo)識(shí)碼 A

文章編號(hào) 0517-6611（2021）04-0088-04

doi：10.3969/j.issn.0517-6611.2021.04.023

Comparison of Morphological Characteristics? of Hemibarbus maculates under Different Growth Environments

LEI Ning1，ZENG Jian-gang2，SHEN Xiao-ming3 et al （1.College of Life Sciences，Huzhou University，Huzhou，Zhejiang 313001;2.Bureau of Agriculture and Rural Affairs in Deqing County，Huzhou，Zhejiang 313100;3.Deqing Haoyuan Aquatic Seed Industry Co.，Ltd.，Huzhou，Zhejiang 313200）

Abstract Three multivariate analysis methods were used to compare 30 morphological characteristics parameters of? three populations of Hemibarbus maculatus，including? two wild populations collected from Taihu Lake and Qiantangjiang River（Tonglu section），Donglin breeding population after 3 generations of artificial breeding.The results showed that the morphological characteristics of Taihu population and Donglin breeding population were similar，which had greater difference with those of Tonglu population.Through principal component analysis，three principal components were established，their contribution rates were 55.08%，14.01%，6.87%，respectively，and the cumulative contribution rate was 75.97%.Three traits with higher contribution rate were selected for discriminant analysis by stepwise discrimination.The discrimination functions of the three populations were set up.The discrimination accuracy P1 was 75%-100%，P2 was 77%-100%，comprehensive discrimination rate was 86%.The stepwise discrimination analysis revealed that the morphological characteristics were significantly different among three populations（P<0.01）.This study accumulated basic data for the identification of different populations of H.maculatus，and provided theoretical basis for the proliferation and release of H.maculatus， the cultivation of godd breed and the protection? of germplasm resources.

Key words Hemibarbus maculates;Morphological characteristics;Cluster analysis;Principal component analysis;Discriminant analysis

基金項(xiàng)目 浙江省重大科技專項(xiàng)（2016C020055-3-3）;國(guó)家水產(chǎn)種質(zhì)資源共享平臺(tái)項(xiàng)目（2018DKA30470）;浙江省重點(diǎn)研發(fā)科研項(xiàng)目（2019C02082）。

作者簡(jiǎn)介 雷寧（1996—），女，廣西柳州人，碩士研究生，研究方向：水產(chǎn)種質(zhì)資源。通信作者，副研究員，博士，碩士生導(dǎo)師，從事水產(chǎn)動(dòng)物遺傳育種研究。

收稿日期 2020-07-02

花鰩骨（Hemibarbus maculates Bleeker） 隸屬鯉形目鯉科鰩骨亞科鰩骨屬，俗稱“季魚”“雞骨郎”“季郎魚”，通常分布于亞洲境內(nèi)的中國(guó)、朝鮮和日本等淡水水域以及俄羅斯境內(nèi)的部分黑龍江水系[1]，該魚肉鮮味美，頗受消費(fèi)者的青睞，經(jīng)濟(jì)價(jià)值較高。隨著人工繁殖和苗種培育技術(shù)的持續(xù)突破[2-3]，花鰩骨在江、浙、滬地區(qū)已成為優(yōu)良的養(yǎng)殖品種。野生花鰩骨自然群體數(shù)量卻不斷削減，全因捕撈過度、水質(zhì)污染及富營(yíng)養(yǎng)化嚴(yán)重和競(jìng)爭(zhēng)性魚類放流等生態(tài)因子變化的影響。為保護(hù)花鰩骨野生資源，浙江省自2015年以來(lái)連續(xù)4年在自然水域中增殖放流花鰩骨魚種，但放流的花鰩骨魚種多來(lái)自于人工繁育群體，其對(duì)自然花鰩骨群體是否會(huì)造成種質(zhì)上的沖擊有待深入探索。

魚類形態(tài)特征是其內(nèi)在遺傳特性的外在表現(xiàn)。從形態(tài)學(xué)和表型性狀上研究遺傳變異是經(jīng)典遺傳學(xué)研究的主要內(nèi)容，是資源鑒定、物種分類、生物起源和進(jìn)化系統(tǒng)樹等研究的基礎(chǔ)[2]。生態(tài)環(huán)境適應(yīng)性是指生物的基礎(chǔ)能力，生物一直處于不斷適應(yīng)其生態(tài)環(huán)境的過程中，并演化產(chǎn)生不同的形態(tài)特征[4]。野生和養(yǎng)殖條件就是2種不同環(huán)境，這2種生態(tài)環(huán)境下花鰩骨的形態(tài)特征尚未見到相關(guān)報(bào)道。筆者對(duì)花鰩骨太湖野生群體、桐廬野生群體和東林養(yǎng)殖群體的形態(tài)參數(shù)進(jìn)行了系統(tǒng)對(duì)比，將傳統(tǒng)形態(tài)度量學(xué)方法與統(tǒng)計(jì)學(xué)比較分析相結(jié)合，以期為花鰩骨不同群體的劃分、良種選育和種質(zhì)資源保護(hù)提供參考。

1 材料與方法

1.1 試驗(yàn)材料 3個(gè)花鰩骨群體的樣本，分別于2019年11月采自浙江錢塘江流域桐廬段，于2018年11月采自浙江湖州太湖流域，于2018年11月采自浙江湖州東林某養(yǎng)殖基地。其中，人工養(yǎng)殖群體是由太湖野生群體經(jīng)人工繁育3代的養(yǎng)殖群體。樣本采集情況見表1。

表1 花鰩骨采集時(shí)間、地點(diǎn)與樣本數(shù)

Table 1 Sampling time，site and sample number of H.maculates

群體Population 采集地 Sampling site 采集時(shí)間Sampling time 樣本數(shù) Sample number? 平均全長(zhǎng)Average whole length∥cm

東林養(yǎng)殖群體Donglin breeding population（DL）湖州東林2018-111931.31±2.27

太湖野生群體Taihu wild population（TH）湖州太湖2018-111921.50±3.01

桐廬放流群體Tonglu released population（TL）桐廬富春江2019-111623.25±4.39

1.2 數(shù)據(jù)測(cè)量 共30項(xiàng)形態(tài)數(shù)據(jù)指標(biāo)，含傳統(tǒng)可量性狀及框架參數(shù)。參照李思發(fā)[5] 進(jìn)行框架測(cè)量坐標(biāo)點(diǎn)的選擇，具體見圖1和表2。魚體外形的21個(gè)框架參數(shù)，以10個(gè)定位點(diǎn)之間的距離構(gòu)建而成，例如D2-4 表示定位點(diǎn)2 與 4 之間的距離。

1.3 分析方法

1.3.1 數(shù)據(jù)預(yù)處理。

每個(gè)采集樣本，測(cè)量包括9個(gè)可量性狀與 21個(gè)框架參數(shù)在內(nèi)的30個(gè)指標(biāo)?？傆?jì)測(cè)量54個(gè)樣本，共得到1 620個(gè)數(shù)據(jù)。將各形態(tài)參數(shù)值除以樣本體長(zhǎng)值再進(jìn)行校正，以消除花鰩骨樣本個(gè)體差異帶來(lái)的誤差或影響。校正后得到29項(xiàng)形態(tài)比例參數(shù)，使用SPSS 23.0統(tǒng)計(jì)軟件進(jìn)行數(shù)據(jù)統(tǒng)計(jì)與分析[6]。

1.3.2 主成分分析。利用主成分分析考察各指標(biāo)間的相互關(guān)系，利用正交變換降維將29項(xiàng)參數(shù)轉(zhuǎn)換為較少且互不相關(guān)的幾個(gè)綜合指標(biāo)，使得后續(xù)研究更加簡(jiǎn)便、直觀。

1.3.3 聚類分析。

用3個(gè)群體29項(xiàng)形態(tài)特征比例參數(shù)的平均值進(jìn)行聚類分析。采用的聚類方法為歐式距離法。

1.3.4 判別分析及判別準(zhǔn)確率的計(jì)算。采用SPSS 23.0統(tǒng)計(jì)軟件的逐步判別法開展判別分析，選取對(duì)群體形態(tài)差異貢獻(xiàn)較大的參數(shù)建立判別公式。判別準(zhǔn)確率的計(jì)算公式如下：

判別準(zhǔn)確率P1=判別正確的魚類樣本數(shù)/實(shí)際魚類樣本數(shù)（1）

P2=判別正確的魚類樣本數(shù)/ 判別魚類樣本數(shù)（2）

綜合判別率=ni=1Ai/ni=1Bi（3）

式中，實(shí)際魚類樣本數(shù)是指該魚類樣本的采樣數(shù)目，判別魚類樣本數(shù)是指被判為該群體的魚類樣本數(shù);Ai 為第 i 個(gè)群體中判別正確的魚類樣本數(shù)，Bi為第 i 個(gè)群體的實(shí)際魚類樣本數(shù)，n 為魚類樣本群體數(shù)。

2 結(jié)果與分析

2.1 主成分分析 在主成分分析中，提取了方差貢獻(xiàn)率較大的3個(gè)主成分，它們的累積貢獻(xiàn)率為75.97%。這說明這3種主成分包含了其總變異的絕大部分，成功實(shí)現(xiàn)了數(shù)據(jù)降維目標(biāo)，可利用這3個(gè)相互獨(dú)立的因子來(lái)分析不同種群間花

鰩骨的形態(tài)差異。這3個(gè)主成分的方差貢獻(xiàn)率分別為主成分1 55.088%、主成分2 14.01%、主成分3 6.87%。

對(duì)主成分分析求得的初始變量相關(guān)矩陣進(jìn)行方差最大化正交旋轉(zhuǎn)，旋轉(zhuǎn)后的因子載荷矩陣見表3。由表3可知，對(duì)總方差貢獻(xiàn)率最大的主成分1主要取決于尾柄長(zhǎng)/體長(zhǎng)、D2-3/體長(zhǎng)、D2-4/體長(zhǎng)、D3-4/體長(zhǎng)、D4-5/體長(zhǎng)、 D5-6/體長(zhǎng)、D7-8/體長(zhǎng)、D8-10/體長(zhǎng)、D9-10/體長(zhǎng)這9個(gè)比例參數(shù)，反映了頭部及尾部等形態(tài)特征，由此可見花鰩骨 3個(gè)群體的個(gè)體在形態(tài)上的差異很大部分是由這9個(gè)形態(tài)差異引起的。

花鰩骨主成分2對(duì)主成分1繪圖，結(jié)果如圖2所示，主成分1、主成分2為相對(duì)值。從圖2可以看出，太湖野生群體和東林養(yǎng)殖群體重疊較多，難以完全區(qū)分開，表明二者差異較小，它們與桐廬放流群體差異較大。

2.2 系統(tǒng)聚類分析

采用系統(tǒng)聚類分析，東林、太湖、桐廬3個(gè)花鰩骨群體間的類平均距離見表4，聚類結(jié)果見圖3。聚類分析將源于同一水域環(huán)境（太湖）的太湖野生群體和人工養(yǎng)殖群體歸于一類，而將不同水域的桐廬群體歸于另一類。以上結(jié)果表明3個(gè)花鰩骨群體中，人工養(yǎng)殖群體和太湖野生群體較為相似，其分析結(jié)果與主成分分析相一致。

2.3 判別分析

采用SPSS 23.0統(tǒng)計(jì)軟件中的逐步判別法進(jìn)行判別分析，對(duì)29項(xiàng)形態(tài)參數(shù)進(jìn)行分析，判別準(zhǔn)確率為95%。為提高公式的實(shí)用性和簡(jiǎn)便性，從29項(xiàng)形態(tài)參數(shù)中選出對(duì)判別貢獻(xiàn)較大的12項(xiàng)參數(shù)，再進(jìn)行判別，判別準(zhǔn)確率P1為75%～100%，判別準(zhǔn)確率P2為77%～100%，3個(gè)花鰩骨群體樣本的綜合判別率為86%（表5） 。

利用挑選的 12個(gè)參數(shù)建立3個(gè)群體的判別公式：

Y1=1 630X1+1 436X2+1 197X3+255X4+37X5+1 741X6+324X7 +1 010X8-414X9+934X10-92X11-407X12-535

Y2=1 573X1+1 396X2+1 257X3+299X4+38X5+1 782X6+263X7 +1 020X8-391X9+989X10-87X11-421X12-539

Y3=1 160X1+2 283X2+1 183X3+375X4+34X5+2 078X6+259X7 +1 246X8-476X9+1 037X10-171X11-479X12-508

式中，X1、X2、 X3、 X4、 X5、 X6、 X7、X8、X9、X10、X11、X12分別為頭長(zhǎng)/體長(zhǎng)、眼徑/體長(zhǎng)、尾柄長(zhǎng)/體長(zhǎng)、吻長(zhǎng)/體長(zhǎng)、尾柄高/體長(zhǎng)、眼間距/體長(zhǎng)、D5-8、D6-8、D1-2、D7-8、D3-6、D5-6。

根據(jù)所得的判別公式，可準(zhǔn)確判別出未知來(lái)源花鰩骨樣本屬于哪個(gè)群體。判別方法如下：測(cè)量出未知來(lái)源的花鰩骨樣本的12個(gè)形態(tài)比例參數(shù)（頭長(zhǎng)/體長(zhǎng)、眼徑/體長(zhǎng)、尾柄長(zhǎng)/體長(zhǎng)、吻長(zhǎng)/體長(zhǎng)、尾柄高/體長(zhǎng)、眼間距/體長(zhǎng)、D5-8、D6-8、D1-2、D7-8、D3-6、D5-6），將所測(cè)比例參數(shù)代入每個(gè)判別函數(shù)，則未知樣本屬于判別函數(shù)值大的群體。

3 討論

3.1 主成分分析、判別分析和聚類分析的應(yīng)用

判別分析廣泛應(yīng)用于各個(gè)領(lǐng)域，它是在分類已知的條件下尋找可以表征事物特點(diǎn)的判據(jù)，通過建立判別函數(shù)對(duì)事物進(jìn)行分類。目前，逐步判別是判別分析中最常用的手段。

判別分析是魚類群體鑒定常用方法之一[7-8] ，也常用于鑒別不同水環(huán)境中河蟹[9-11]和羅非魚的不同品系[12-13] 。目前，判別分析在水產(chǎn)領(lǐng)域的應(yīng)用已取得了一些成果，比如建立了不同水系河蟹的判別公式，利用判別函數(shù)區(qū)分不同水系的河蟹，對(duì)于避免河蟹種質(zhì)混雜上有直接作用。

該研究采用可量性狀與框架數(shù)據(jù)相結(jié)合的方法，對(duì)花鰩骨的3個(gè)群體進(jìn)行了判別分析。源于不同水域環(huán)境（太湖、桐廬）的群體判別準(zhǔn)確率100%，判別效果顯著。這說明判別分析對(duì)于親緣關(guān)系很近的群體的鑒別是很有效的，而該研究所建立的判別公式可用于鑒別不同水域的花鰩骨群體，在花鰩骨種質(zhì)資源保護(hù)上具有一定的借鑒意義。

該研究通過聚類分析、判別分析和主成分分析3種多元分析方法進(jìn)行統(tǒng)計(jì)分析，分析結(jié)果均具有相似性。主成分分析是將多種指標(biāo)歸納為少量具有代表性的因子來(lái)說明不同群體之間的差異大小，并根據(jù)不同種群的主成分值找出其各自差異較大的參數(shù)。該研究利用主成分分析得到的第一、二主成分繪制散點(diǎn)圖，直觀地分析了不同群體的區(qū)分程度，為后續(xù)數(shù)據(jù)分析提供了依據(jù)。聚類分析是參照遺傳距離的大小進(jìn)行劃分，遺傳距離較小的2個(gè)品種歸為一類，遺傳距離較大的則分在不同的類。該研究的聚類分析的結(jié)果表明同一水域環(huán)境（太湖）的太湖野生群體和人工養(yǎng)殖群體歸于一類，而將不同水域的桐廬群體歸于另一類，是較為合理可信。這個(gè)分類結(jié)果也在判別分析中得到進(jìn)一步的驗(yàn)證，即源于不同水域環(huán)境（太湖、桐廬）的群體判別準(zhǔn)確率達(dá)100%，判別效果顯著，而源于同一水域的太湖野生和東林人工養(yǎng)殖判別準(zhǔn)確率分別為94.94%和100%，表明判別分析與聚類分析的趨同一致性。

3.2 ?野生與人工養(yǎng)殖的生存環(huán)境對(duì)花鰩骨形態(tài)特征的影響

有關(guān)野生與人工養(yǎng)殖的2種不同生存環(huán)境對(duì)魚類生長(zhǎng)的影響的研究，以往主要集中在肌肉營(yíng)養(yǎng)成分的差異分析，其結(jié)果表明不同的生存環(huán)境和餌料使得野生和養(yǎng)殖群體的肌肉成分和組成存在顯著差異[14-16]。

形態(tài)學(xué)的特征是受遺傳因子和環(huán)境因子的共同影響。該研究結(jié)果表明野生和養(yǎng)殖環(huán)境的差異對(duì)花鰩骨形態(tài)存在影響，二者形態(tài)差異主要取決于尾柄長(zhǎng)/體長(zhǎng)、D2-3/體長(zhǎng)、D2-4/體長(zhǎng)、D3-4/體長(zhǎng)、D4-5/體長(zhǎng)、D5-6/體長(zhǎng)、D7-8/體長(zhǎng)、D8-10/體長(zhǎng)、D9-10/體長(zhǎng)這9個(gè)比例參數(shù)。這可能是因?yàn)樵撗芯恐械娜斯ゐB(yǎng)殖群體是由太湖野生群體經(jīng)人工馴養(yǎng)而來(lái)，僅經(jīng)過數(shù)代的繁育，生存環(huán)境的變化對(duì)該群體的影響時(shí)間較短，并未導(dǎo)致明顯的形態(tài)特征差異，所以源于同一水域環(huán)境（太湖）的太湖野生群體和人工養(yǎng)殖群體更為相似。

基因組遺傳多樣性分析結(jié)果表明不科學(xué)的人工繁育會(huì)降低繁育群體的遺傳多樣性水平[17-22]。源于太湖野生群體F2代以內(nèi)的人工繁育苗種可用于太湖流域的花鰩骨增殖放流，但不建議用于錢塘江水系花鰩骨的增殖放流。

參考文獻(xiàn)

[1]中國(guó)科學(xué)院中國(guó)動(dòng)物志編輯委員會(huì).中國(guó)動(dòng)物志：硬骨魚綱鯉形目（中卷）[M].北京：科學(xué)出版社，1998.

[2] 黃邦星.花鰩骨人工繁殖及苗種培育試驗(yàn)[J].漁業(yè)致富指南，2016（11）：48-52.

[3] 余連渭，夏前征，王英.花（鰩骨）水花魚苗培育夏花魚種試驗(yàn)[J].江西水產(chǎn)科技，2015（2）：26-27.

[4]庹云.袁河花（鰩骨）形態(tài)生物學(xué)研究[J].安徽農(nóng)業(yè)科學(xué)，2012，40（36）：17605-17607.

[5] 李思發(fā).中國(guó)淡水主要養(yǎng)殖魚類種質(zhì)研究[M].上海：上?？茖W(xué)技術(shù)出版社，1998：181-210.

[6] 盧紋岱.SPSS for Windows統(tǒng)計(jì)分析[M].3版.北京：電子工業(yè)出版社，2006.

[7] MAYR E，GORTON LINSLEY E，USINGER R L.Methods and principles of systematic zoology[M].New York：McGraw Hill，1953：23-39，125-154.

[8] 李思發(fā)，吳力釗，王強(qiáng)，等.長(zhǎng)江、珠江、黑龍江鰱、鳙、草魚種質(zhì)資源研究[M].上海：上?？茖W(xué)技術(shù)出版社，1990：145-162.

[9] 趙金良，李思發(fā).中國(guó)大陸沿海六水系絨螯蟹（中華絨螯蟹和日本絨螯蟹） 群體親緣關(guān)系：生化遺傳差異分析[J].水產(chǎn)學(xué)報(bào)，1999，23（4）：331-336.

[10] 許加武，任明榮，李思發(fā). 長(zhǎng)江、遼河、甌江中華絨鰲蟹種群的形態(tài)判別[J].水產(chǎn)學(xué)報(bào)，1997，21（3）：269-274.

[11] 李晨虹，李思發(fā).中國(guó)大陸沿海六水系絨螯蟹（中華絨螯蟹和日本絨螯蟹）群體親緣關(guān)系：形態(tài)判別分析[J].水產(chǎn)學(xué)報(bào)，1999，23（4）：337-342.

[12] 李思發(fā)，李晨虹，李家樂.尼羅羅非魚品系間形態(tài)差異分析[J].動(dòng)物學(xué)報(bào)，1998，44（4）：450-457.

[13] 李家樂，李思發(fā)，李勇，等.尼奧魚[尼羅羅非魚 （♀）×奧利亞羅非魚 （♂）同其親本的形態(tài)和判別[J].水產(chǎn)學(xué)報(bào)，1999，23（3）：261-265.

[14] 陳建明，葉金云，沈斌乾，等.野生和池塘養(yǎng)殖花（鰩骨）肌肉營(yíng)養(yǎng)組成的比較分析[J].上海水產(chǎn)大學(xué)學(xué)報(bào)，2007，16（1）：87-91.

[15] 宋超，莊平，章龍珍，等.野生及人工養(yǎng)殖中華鱘幼魚肌肉營(yíng)養(yǎng)成分的比較[J].動(dòng)物學(xué)報(bào)，2007，53（3）：502-510.

[16] 王偉，陳立僑，顧志敏，等.六個(gè)群體翹嘴紅魚白肌肉生化組成的比較[J].水產(chǎn)學(xué)報(bào)，2007，31（S1）：92-99.

[17] 孟憲紅，孔杰，莊志猛，等.真鯛自然群體和人工繁殖群體的遺傳多樣性[J].生物多樣性，2000，8（3）：248-252.

[18] 王軍，全成干，蘇永全，等.宮井洋大黃魚遺傳多樣性的RAPD 分析[J].海洋學(xué)報(bào)，2001，23（3）：87-91.

[19] 舒琥，蒙子寧，易祖盛，等.唐魚野生與養(yǎng)殖群體遺傳多樣性的隨機(jī)擴(kuò)增多態(tài)DNA（RAPD）分析[J].中山大學(xué)學(xué)報(bào)（自然科學(xué)版），2006，45（1）：77-81.

[20] 陳雪峰，楊國(guó)梁，孔杰，等.人工養(yǎng)殖與選育對(duì)羅氏沼蝦遺傳多樣性的影響[J].水生生物學(xué)報(bào)，2012，36（5）：866-873.

[21] 常玉梅，王文文，徐萬(wàn)土，等.人工繁育大黃魚（Pseudosciaena crocea）群體F2及F3遺傳差異分析[J].海洋與湖沼，2009，40（4）：414-422.

[22] 林茂，黃景，李正，等.大鯢野生親代與人工繁育子二代的隨機(jī)擴(kuò)增多態(tài)DNA分析（英文）[J].上海水產(chǎn)大學(xué)學(xué)報(bào)，2003，12（S1）：20-23.