安 麗，孟慶磊，張龍崗，董學颯，李 嫻，李 寧，朱樹人

(山東省淡水漁業(yè)研究院，山東省淡水水產(chǎn)遺傳育種重點實驗室，濟南 250117)

翹嘴紅鲌雌雄個體的形態(tài)差異分析

安 麗，孟慶磊，張龍崗，董學颯，李 嫻，李 寧，朱樹人

(山東省淡水漁業(yè)研究院，山東省淡水水產(chǎn)遺傳育種重點實驗室，濟南 250117)

為研究翹嘴紅鲌(Eryghroculterilishaeformis)主要形態(tài)指標和雌雄個體的形態(tài)差異，2016年3月采集了67尾翹嘴紅鲌，對其全長/體長、體高/體長、體寬/體長、尾柄長/體長、尾柄高/體長、眼間距/頭長、頭高/頭長、頭寬/頭長、吻長/頭長和豐滿度10項標準化比例性狀進行主成分分析、R-聚類分析和逐步判別分析，并建立了雌雄判別方程式。主成分分析和R-聚類分析結果顯示，翹嘴紅鲌主要形態(tài)指標可分為肥瘦程度特征、軀干部特征、尾部特征和頭部特征4個方面。散布圖結果顯示雌雄個體在尾部形態(tài)上存在差異，同判別分析結果一致。通過逐步判別法從67尾翹嘴紅鲌的10項標準化性狀中篩選出體高/體長、體寬/體長、尾柄高/體長和豐滿度4個比例性狀，建立性別判別方程，對群體識別正確率達91.0%。經(jīng)T檢驗，除體高/體長外，其他3個比例性狀在兩性群體間均存在顯著或極顯著差異，表明雌性相較于雄性體型較寬、較豐滿，但尾柄較矮。上述差異性狀及雌雄判別方程式可為翹嘴紅鲌性別的鑒定提供方法和理論依據(jù)。

翹嘴紅鲌(Eryghroculterilishaeformis); 形態(tài)差異; 雌雄判別

在魚類的人工繁殖育種過程中，要實現(xiàn)雌雄個體分養(yǎng)、保證催產(chǎn)合理性比、提高人工繁殖及選育成效等都需要對性別的準確鑒別[1、2]。對于兩性異形明顯的魚類[3-6]，可直接用肉眼區(qū)分雌性和雄性。但對于無性征的魚類則很難通過外型進行鑒定，需通過解剖觀察、制作性腺切片的方法鑒別，但這種方法需解剖魚體，既不利于稀有動物的保護，也難以實現(xiàn)性別鑒別后的再利用。近年來隨著分子育種技術的發(fā)展，通過分子遺傳標記可準確的鑒別出雌雄個體[7-9]，但這種方法需早期開發(fā)性別特異標記，在實際操作中應用性較差，代價較高。因此在實際生產(chǎn)和操作中，通過魚類傳統(tǒng)形態(tài)指標的測量和多元統(tǒng)計分析成為了一種簡單易行、實用有效的鑒別雌雄個體的方法。

翹嘴紅鲌(Eryghroculterilishaeformis)，又稱白條、和順、太湖白魚，是“太湖三白”之一，屬于鯉形目鯉科鲌亞科紅鲌屬，是一種大型經(jīng)濟淡水魚類，分布廣泛，在我國黑龍江、黃河、長江、珠江等水系均有分布 。該種具有個體大、生長迅速、抗病力強、肉質(zhì)白嫩、味道鮮美、營養(yǎng)價值豐富等特點，深受廣大百姓的歡迎和喜愛[10-12]。本研究通過測量翹嘴紅鲌形態(tài)指標，采用多元統(tǒng)計分析尋找雌雄形態(tài)差異，建立判別雌雄個體的數(shù)學模型，旨在為以后翹嘴紅鲌繁殖、選育過程中提供鑒別雌雄個體的方法和理論依據(jù)。

1 材料與方法

1.1 實驗材料

2016年3月，實驗用翹嘴紅鲌取自山東省淡水漁業(yè)研究院構建的選育群系，隨機撈取67尾翹嘴紅鲌，其中雄魚31尾，雌魚36尾，18月齡。雌雄性別通過解剖觀察其性腺發(fā)育特征來確定。

1.2 參數(shù)測量

使用60 mg/L的丁香油將翹嘴紅鲌麻醉，吸干體表多余的水分，測量體重及其他形態(tài)性狀。用直尺測量全長(TL)和體長(BL)，精確到0.1 cm；用數(shù)顯游標卡尺測量體高(BH)、體寬(BB)、眼間距(ED)、頭長(HL)、頭高(HH)、頭寬(HB)、吻長(SL)、尾柄長(CL)和尾柄高(CH)等形態(tài)性狀指標，精確到0.01 cm；使用電子天平秤取體重(BW)，精確到0.1 g。豐滿度(F)通過BW/BL3計算得到。

1.3 數(shù)據(jù)處理

為消除個體大小對測量參數(shù)的影響，將軀干部的測量性狀除以體長、頭部的測量性狀除以頭長，轉化為標準化性狀，連同豐滿度共計10項性狀(表1)，進行主成分分析、R-聚類分析和逐步判別分析，篩選出雌雄群體間具有顯著差異的性狀建立翹嘴紅鲌雌雄個體判別方程式。應用Excle 2010和SPSS 16.0統(tǒng)計軟件對數(shù)據(jù)進行處理和分析。

表1 標準化性狀參數(shù)

2 結果

2.1 雌雄形態(tài)性狀分析

翹嘴紅鲌各形態(tài)性狀及其變異系數(shù)見表2。由表2可知，實驗魚雄魚個體全長為230.0～360.0 mm，體重為89.8～328.2 g；雌性個體全長為245.0～425.0 mm，體重為125.8～496.4 g。測量的12個性狀中，其中7個性狀雌魚個體的變異大于雄魚，雌雄個體變異最大的均為體重，變異系數(shù)分別為0.339和0.400，雄魚相較于雌魚略顯均勻。

表2 翹嘴紅鲌各形態(tài)性狀統(tǒng)計量

2.2 主成分分析

對翹嘴紅鲌雌雄個體的10個標準化性狀進行主成分分析，共獲得4個主成分，累積貢獻率為70.17%，4個主成分的貢獻率以及各性狀在這4個主成分上的負荷量見表3。其中，主成分1的貢獻率為27.29%，主要包括體高/體長、體寬/體長、頭寬/頭長和豐滿度4個性狀，反映了翹嘴紅鲌體型的肥瘦程度特征；主成分2的貢獻率為16.05%，主要包括全長/體長和眼間距/頭長2個性狀，反映了翹嘴紅鲌的軀干及眼間距特征；主成分3的貢獻率為14.10%，主要包括尾柄長/體長和尾柄高/體長2個性狀，反映了翹嘴紅鲌的尾部特征；主成分4的貢獻率為12.73%，主要包括頭高/頭長和吻長/頭長2個性狀，反映了翹嘴紅鲌的頭部特征。

根據(jù)主成分分析的得分，分別對主成分1與主成分2及主成分1與主成分3之間的雌雄個體繪制散布圖，見圖1和圖2。圖1中雌雄個體分布在主成分2軸上表現(xiàn)出稍明顯的分離態(tài)勢，通過得分可知差異主要由眼間距造成，圖2中雌雄個體在主成分3上稍有分離，比較尾柄高和尾柄長的得分，兩者影響相當，散布圖顯示雌雄個體可能在眼間距和尾部形態(tài)上稍有差異。

表3 4個主成分貢獻率及各性狀負荷量

圖1 翹嘴紅鲌雌雄個體在主成分1和主成分2的分布Fig.1 The distribution of principal component 1 and 2 of E.ilishaeformis

圖2 翹嘴紅鲌雌雄個體在主成分1和主成分3的分布Fig.2 The distribution of principal component 1 and 3 of E.ilishaeformis

2.3 R-聚類分析

對67尾翹嘴紅鲌的10個標準化形態(tài)性狀進行R-聚類分析，結果見圖3。由圖可見，以上性狀分為兩類，第一類為X3、X10、X2、X8、X1、X4和X5，主要反映了體型的特征參數(shù)；第二類為X6、X9和X7，主要反映了頭部的特征參數(shù)。第一類又細分為2小類，X3、X10、X2、X8和X1組成第1小類，X4和X5組成第2小類。

R-聚類分析中第一類中的第1小類與主成分分析所得第1主成分和第2主成分相對應，反映軀干形態(tài)特征；第一類中的第2小類與第3主成分相對應，反映了尾部形態(tài)特征；第二類與第4主成分相對應，反映了頭部形態(tài)特征，兩種分析結果大體一致。

圖3 翹嘴紅鲌標準化形態(tài)性狀R-聚類分析圖Fig.3 R-Cluster dendrogram of standardized morphological indexes for E.ilishaeformis

2.4 判別分析

對67尾翹嘴紅鲌個體進行逐步判別分析，根據(jù)各性狀對模型的貢獻大小，逐步剔除不相關的性狀，最終篩選出4個變量，即X2、X3、X5和X10。建立判別方程式為：

雄：F1=1116X2+1335X3+6148X5-86.525X10-389.941

雌：F2=865.707X2+1831X3+5765X5-70.369X10-369.042

判別函數(shù)的顯著性檢驗顯示，P<0.01，表明判別函數(shù)達到極顯著水平。計算出雌雄個體的判別分數(shù)值，分別得到雌雄個體的頻布圖，見圖4和圖5。由圖可見，只有少數(shù)的個體出現(xiàn)錯判的現(xiàn)象，該模型可以有效地區(qū)分翹嘴紅鲌的性別。

圖4 鲌雌性個體判別分數(shù)的頻布圖Fig.4 Frequency distributions of discriminant scores obtained from female E.ilishaeformis

圖5 雄性個體判別分數(shù)頻布圖Fig.5 Frequency distributions of discriminant scores obtained from male E.ilishaeformis

利用所建立的判別方程式對翹嘴紅鲌67尾個體進行性別判別，將每尾個體的4個性狀參數(shù)值帶入方程式計算F1和F2，若F1>F2，判定為雄魚，反之則為雌魚。判別結果見表4，雄性判別準確率達93.5%，雌性判別準確率達88.9%，綜合判別成功率達91.0%，67尾翹嘴紅鲌僅6尾被錯判為其它性別。

表4 翹嘴紅鲌雌雄判別結果

對篩選出的用于建立判別函數(shù)的4個標準化性狀進行T檢驗，結果見表5。結果顯示體寬/體長和尾柄高/體長的P<0.01，差異極顯著，豐滿度的P<0.05，差異顯著，體高/體長的P>0.05，差異不顯著。以上數(shù)據(jù)表明，雌性個體的體型較寬較豐滿，但尾柄比較矮。

表5 翹嘴紅鲌性狀的T檢驗及雌雄形態(tài)差異

3 討論

3.1 多元統(tǒng)計分析在雌雄差異分析上的應用

隨著計算機技術的發(fā)展及各種統(tǒng)計軟件的開發(fā)，多元統(tǒng)計分析方法也廣泛應用于水產(chǎn)科學領域。主成分分析是在不損失或很少損失原有信息的基礎上，將多個實測變量簡化為較少變量，更能反映事物本質(zhì)的一種分析方法。聚類分析和判別分析都是研究事物分類的基本方法，但聚類分析是根據(jù)事物本身的特性研究個體分類的方法，而判別分析是先求出判別函數(shù)，然后根據(jù)判別函數(shù)對未知所屬類別的事物進行分類的一種分析方法[13]。近年來，許多學者將主成分分析、R-聚類分析和判別分析等多元分析的方法應用于水產(chǎn)動物雌雄個體形態(tài)差異的分析上。郭弘藝等[14]對153尾銀色鰻的16個直接測量性狀進行主成分分析，17項標準化比例性狀進行R-聚類分析，還建立了雌雄個體的判別模型，得出結果，長江口銀色鰻的雌性個體比雄性更為豐滿，體型較高但眼較小。吳波等[15]對178尾金錢魚形態(tài)性狀進行主成分分析和R-聚類分析，建立了金錢魚雌雄判別方程，結果表明雌性個體相比于雄性具有體形較高、頭較大且吻較長的特征。諶微等[16]分析了554尾大黃魚12項計量性狀及11項標準化性狀，雌魚個體間變異較雄魚更大，雌雄魚的體長/體厚差異最大。岳亮等[17]采用主成分分析、逐步判別分析和T檢驗的方法對紅鰭東方鲀雌雄個體的形態(tài)差異進行了分析研究，結果表明紅鰭東方鲀雄性與雌性相比，體型較寬，體周長2較長。

3.2 翹嘴紅鲌雌雄差異分析

本研究對67尾翹嘴紅鲌10項標準化性狀進行多元統(tǒng)計分析，得到的4個主成分分別代表了翹嘴紅鲌肥瘦程度特征、軀干部特征、尾部特征和頭部特征，通過對這些指標的測定和分析可以基本反映出翹嘴紅鲌的外部形態(tài)。聚類分析中第一類中的第1小類與主成分分析所得第1主成分和第2主成分相對應，反映軀干形態(tài)特征；第一類中的第2小類與第3主成分相對應，反映了尾部形態(tài)特征；第二類與第4主成分相對應，反映了頭部形態(tài)特征，除眼間距分析結果不同以外其他分析結果相一致。在主成分分析的散布圖可見雌雄個體在主成分2軸上表現(xiàn)出稍明顯的分離態(tài)勢，說明雌雄個體在主成分2的指標有差異。根據(jù)得分情況可知應該是眼間距對此分離的作用較大，分析原因可能是在實際操作中眼間距的測量值較小，人為操作時產(chǎn)生的誤差造成的影響更顯突出，也因此造成兩種分析在眼間距上有差異。由散布圖可見在主成分3上雌雄個體也稍有分離，說明雌雄個體在尾部形態(tài)上存在差異，這一結果在下面的判別分析中也可以得到驗證。運用逐步判別分析的T檢驗結果表明翹嘴紅鲌雌雄差異主要集中體型肥瘦程度和尾部形態(tài)上，即雌性相較于雄性體型顯著較寬、較豐滿，但尾柄較矮。同時建立了雌雄判別方程式，最終保留了體高/體長、體寬/體長、尾柄高/體長和豐滿度4個變量，綜合判別成功率達91.0%，說明通過以上4個性狀指標可以快速鑒別出翹嘴紅鲌性別。

3.3 存在的問題

水產(chǎn)動物中雌雄個體間的差異主要包括雌雄形態(tài)差異和生長差異兩部分，隨著年齡的增長雌雄間的差異會有不同的變化，建立的判別雌雄模型或只適用于特定的體尺度范圍之內(nèi)或特定的體重范圍之內(nèi)。吳波等[15]認為所篩選出的5個外形特征只可以判別體重在100 g以上的金錢魚個體的性別，100 g以下的還需要進一步驗證。岳亮等[17]驗證了新選用的兩個東方紅鰭鲀雌雄群體均不適用于已建立的雌雄判別模型。本研究通過測量、分析翹嘴紅鲌外部形態(tài)指標，建立了判別雌雄的數(shù)學模型，進而得出雌性相較于雄性較寬較豐滿，尾柄較矮，但陳世通等[18]研究中認為翹嘴紅鲌個體大小并無顯著的兩性差異，局部形態(tài)特征方面也不存在顯著的兩性差異。說明本研究建立的判別模型只適用于體長為20～37 cm的翹嘴紅鲌，對作者陳世通等[18]選用的體長為7.1～19.4 cm的魚體不適用，因此尋找不同雌雄差異的體尺度適用范圍或體重的適用范圍成為下一步需要研究的方向。

自然水域中的翹嘴紅鲌是一種大型兇猛肉食性魚類，其生長速度、個體大小、雌雄差異等會受到自身性選擇、水生環(huán)境因子(食物的豐度、水中物種的分布)、年齡差別等因素的影響。本研究選用了翹嘴紅鲌的選育群系，雖然選擇強度較小但其生長的環(huán)境、攝食模式等發(fā)生較大變化，因此野生群體和選育群體可能存在一定的差別，需進一步比較研究。

[1]湯亞斌,馬達文,易 翀,等.大鯢親本性別鑒定的幾種方法[J].中國水產(chǎn),2012,(10):67.

[2]張 濤,顏世偉,莊 平,等.西伯利亞鱘性別及性腺發(fā)育的超聲波鑒定[J].水產(chǎn)學報,2010,34(1):72-79.

[3]李潤潼.魚類的性別[J].生物學通報,1993,28(11):8-9.

[4]Echeverria T W.Sexual dimorphism in four species of rockfish genus Sebastes (Scorpaenidae) [J].Environ Biol Fish, 1986,15(3):181-190.

[6]倪海兒.東海鰳魚雌雄性狀比較及其鑒別[J].浙江水產(chǎn)學院學報,1994,13(2):93-98.

[7]Chen S L, Li J,Deng S P,et al. Isolation of female-specific AFLP markers and molecular identification of genetic sex in half-smooth tongue sole (Cynoglossussemilaevis)[J].Mar Biotechnol,2007,9(2):273-280.

[8]楊 玲,孟慶磊,李 嫻,等. 一個奧里亞羅非魚雌性特異性SCAR標記的建立[J].漁業(yè)科學進展,2011,32(4):34-40.

[9]劉改艷,陳昆慈,鄭光明,等.SSR-BSA技術對烏鱧性別差異標記的初步篩選[J].水產(chǎn)學報,2011,35(2):170-175.

[10]孟慶磊,朱樹人,李 寧,等.黃河翹嘴鲌染色體核型研究[J].齊魯漁業(yè),2016,33(6):5-7.

[11]張永正.翹嘴紅鲌人工繁殖技術研究[D].浙江：浙江海洋學院水產(chǎn)養(yǎng)殖系,2005.

[12]楊 玲,李 寧,孟慶磊.黃河翹嘴紅鲌同工酶組織特異性研究[J].山東師范大學學報(自然科學版),2015,30(4):99-103.

[13]張慶利.SPSS寶典 [M]，第2版.北京：電子工業(yè)出版社,2011:327-403.

[14]郭弘藝,魏 凱,謝正麗,等.長江口銀色鰻的形態(tài)指標體系及其雌雄鑒別[J].水產(chǎn)學報,2011,35(1):1-9.

[15]吳 波,張敏智,鄧思平,等.金錢魚雌雄個體的形態(tài)差異分析[J].上海海洋大學學報,2014,23(1):64-68.

[16]諶 微,王盼盼,肖世俊,等.大黃魚形態(tài)指標體系及雌雄差異分析[J].集美大學學報(自然科學版),2014,19(6):401-408.

[17]岳 亮,王新安,馬愛軍,等.紅鰭東方鲀(Takifugurubripes)雌、雄個體的形態(tài)特征比較[J].漁業(yè)科學進展,2016,37(1):30-35.

[18]陳世通,趙仁友,張亞平,等.翹嘴紅鲌的形態(tài)兩性異性和雌性生育力[J].麗水學院學報,2013,35(2):11-15.

(責任編輯：鄧 薇)

Analysis of morphological indexes and discriminationof male and femaleEryghroculterilishaeformis

AN Li, MENG Qing-lei, ZHANG Long-gang, DONG Xue-sa, LI Xian, LI Ning, ZHU Shu-ren

(ShandongFreshwaterFisheriesResearchInstitute,ShandongProvincialKeyLaboratoryofFreshwaterGeneticsandBreeding,Ji’nan250117,Shandong,China)

To study the main morphological indexes ofEryghroculterilishaeformisand morphological variations between male and female, a total of 12 quantitative characters of 67 individuals collected in March 2016 were measured. 10 characteristics(total length/body length, body height/body length, body breadth/body length, caudal peduncle length/body length, caudal peduncle height/body length, eye distance/head length, head height/head length, head breadth/head length, snout length/head length and fatness) were normalized and further used for Principal Component Analysis, R-cluster Analysis and Stepwise Discriminant Analysis. Principal Component Analysis and R-cluster Analysis indicated that the morphological index system ofE.ilishaeformiscould be classified into the level of obesity, the body feature, the caudal peduncle and head feature. The distribution of principal component 1 and 3 showed thatE.ilishaeformiscould be classified by caudal peduncle feature, which was in keeping with Stepwise Discriminant Analysis. Four normalized traits involving body height/body length, body breadth/body length, caudal peduncle height/body length and fatness were screened out by Stepwise Discriminant Analysis based on the contribution of variables to the discriminant function. On the basis of the four normalized traits, sex discriminant equations were established and the accuracy rate for identifying 67 samples was 91.0%.Further, the T-text results of four normalized traits between male and female were significantly different (P<0.05 orP<0.01), except for the feature of body hepth/body length. Therefore, we concluded that the femaleE.ilishaeformiswere superior to male in the respect of body breadth and fatness but having smaller caudal peduncle height than that of male. In conclusion, the differential characters and the discriminant equations for male and female could be a simple and efficient method for sex identification ofE.ilishaeformis.

E.ilishaeformis; morphological variation; sex discrimination

2016-10-10；

2017-01-17

山東省農(nóng)業(yè)良種工程項目“翹嘴紅鲌優(yōu)異基因資源發(fā)掘與創(chuàng)新種質(zhì)培育”；山東省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術體系魚類創(chuàng)新團隊(SDAIT-12-01)；國家科技基礎條件平臺“水產(chǎn)種質(zhì)資源平臺(2016DKA30470)”子課題“山東省淡水水產(chǎn)種質(zhì)資源整理、整合與共享服務”

安 麗(1982-)，女，碩士，助理研究員，主要從事水產(chǎn)動物遺傳育種研究。E-mail: ali325@163.com 通訊作者：孟慶磊。E-mail:qingleimeng@126.com

S931

A

1000-6907-(2017)02-0036-06