張哲 楊章武 鄭雅友 李正良 葛輝 楊志宏 許智海

摘要 為研究凡納濱對蝦選育群體自交和雜交F1的形態(tài)親緣關系，采用方差分析和多元統計分析方法，通過12個測量性狀比較F1各組和普通苗組的外部形態(tài)特征。方差分析結果顯示，F1各組之間、與普通苗組之間在部分形態(tài)特征上表現出明顯差異，主要體現在F1各組的平均鮮出肉率比普通苗組提高了6.12%，普通苗組表現出較高的體長變異系數和較低的存活率等特征。通過主成分分析共構建了4個主成分，累計貢獻率為92.33%。第一主成分主要反映頭胸甲長（X3）、頭胸甲寬（X4）、第一腹節(jié)長（X5）、第二腹節(jié)長（X6）、第六腹節(jié)長（X7）的特征，方差貢獻率為40.14%；第二主成分主要反映肥滿度指數Ⅰ（X9）、肥滿度指數Ⅱ（X10）、重長指數Ⅰ（X11）、重長指數Ⅱ（X12）的特征，方差貢獻率為25.34%。逐步判別分析顯示，所建Fisher分類函數方程組可清晰區(qū)分普通苗組、BB（B♀×B♂）、AB（A♀×B♂）和BA（B♀×A♂）各組對蝦群體，其判別準確率分別為96%、75%、84%和97%，但AA組（A♀×A♂）的判別準確率較低，僅為13%。聚類分析結果顯示，普通苗組與選育群體F1代明顯分為2支，F1代各組中自交組間親緣關系最近，正交組與自交組間親緣關系較近，反交組與各組親緣關系較遠。

關鍵詞 凡納濱對蝦；雜交選育；形質差異；多元統計分析

中圖分類號 S968? 文獻標識碼 A? 文章編號 0517-6611（2023）13-0081-05

doi：10.3969/j.issn.0517-6611.2023.13.019

Morphological Variation Analysis of Self-crossing and Hybrid F1 from Two Breeding Populations of Litopenaeus vannamei

ZHANG Zhe，YANG Zhang-wu，ZHENG Ya-you et al

（Key Laboratory of Cultivation and High-value Utilization of Marine Organisms in Fujian Province/Fisheries Research Institute of Fujian，Xiamen，Fujian 361013）

Abstract In order to study the morphological relationship between self-crossing and hybrid F1 of Litopenaeus vannamei breeding population，the external morphological characteristics of F1 group and common seedling group were compared by 12 measured traits by using variance analysis and multivariate statistical analysis.The variance analysis results showed that there were significant differences in some morphological characteristics between F1 groups and ordinary seedling groups，which mainly reflected that the average fresh meat rate of F1 group was 6.12% higher than that of ordinary seedling group.And the ordinary seedling group showed higher variation coefficient of body length and lower survival rate.Through the principal component analysis，four principal components were constructed，the cumulative contribution rate was 92.33%.The first principal component mainly reflected the characteristics of carapace length （X3），carapace width （X4），the length of the first pleon segment （X5），the length of the second pleon segment （X6），the length of the sixth pleon segment length （X7） ，the variance contribution rate was? 40.14%.The second principal component mainly reflected the characteristics of fatness index Ⅰ （X9），fatness index Ⅱ （X10），weight-length index I （X11） and weight-length index Ⅱ （X12），the variance contribution rate was 25.34%. Step discriminant analysis showed that the Fisher classification function equations could clearly distinguish the groups of common seedling group，BB （B♀×B♂），AB （A♀×B♂） and BA （B♀×A♂），with the discriminant accuracy of 96%，75%，84% and 97%，respectively.But the discriminant accuracy of AA （A♀×A♂） was lower，only 13%.The results of cluster analysis showed that the common seedling group and the breeding population F1 generation were divided into two branches.Among the F1 generation groups，the self-crossing group had the closest relationship，the relationship between the direct cross group and the reciprocal cross group was closer，the reciprocal cross had farther relationship with each group.

Key words Litopenaeus vannamei;Hybridization;Morphological variation;Multivariate statistical analysis

基金項目 福建省種業(yè)工程項目（2017FJSCZY02）；2021年福建省海洋與漁業(yè)結構調整專項（2021HYJG16）。

作者簡介 張哲（1985—），男，安徽宣城人，助理研究員，碩士，從事對蝦選育與增養(yǎng)殖研究。通信作者，研究員，從事海水養(yǎng)殖與水產育種研究。

收稿日期 2022-06-12

凡納濱對蝦（Litopenaeus vannamei）又名南美白對蝦，俗稱白蝦，原產于南美洲太平洋沿岸的暖水水域，具有養(yǎng)殖產量高，適溫、鹽范圍廣，生長速度快等特點，是世界三大養(yǎng)殖對蝦中單產量最高的蝦種，也是我國對蝦養(yǎng)殖最主要的種類。2020年我國凡納濱對蝦海水養(yǎng)殖119.8萬t、淡水養(yǎng)殖66.5萬t，凡納濱對蝦養(yǎng)殖產量占當年對蝦養(yǎng)殖產量的86.5%。福建省也是對蝦養(yǎng)殖的重要省份之一，2020年對蝦養(yǎng)殖產量22.6萬t，其中凡納濱對蝦20.3萬t，占對蝦養(yǎng)殖產量的89.8%。凡納濱對蝦蝦苗數量4 269億尾，占全國對蝦養(yǎng)殖數量的27.3%［1］。苗種是水產增養(yǎng)殖產業(yè)的源頭，對蝦種苗問題是對蝦養(yǎng)殖業(yè)發(fā)展的核心問題。近年來，國內各大研究機構開展了凡納濱對蝦人工育種研究，已成功培育出9個凡納濱對蝦新品種［2-5］。其中，有關凡納濱對蝦的遺傳多樣性研究，國內外主要集中在SSR［6-8］、SNP分子標記［9-10］等方面對選育群體或家系的遺傳結構進行分析，但對傳統的形態(tài)特征研究較少。目前利用多元統計分析方法對種群形態(tài)特征的分析已在魚類［11］、蟹類［12］、貝類［13］和蝦類［14］中廣泛應用。然而，利用多元統計分析方法分析凡納濱對蝦不同選育群體形質差異的研究還較少。目前國家大力提倡種業(yè)創(chuàng)新，鼓勵水產良種自主可控，因此不同選育世代遺傳特性研究在對蝦選育過程中至關重要，是實現良種創(chuàng)新的有效途徑。筆者培育凡納濱對蝦選育后代，通過形態(tài)學測定，采用主成分分析、判別分析和聚類分析3種多元統計分析方法分析了2個凡納濱對蝦選育群體自交和雜交F1的形態(tài)質量特征差異，旨在為后續(xù)良種本地化選育和親緣關系研究提供理論依據和方法參考。

1 材料與方法

1.1 試驗材料

A組泰國選育系和B組海南選育系2個凡納濱對蝦群體均通過商業(yè)渠道獲得。特征性狀均為速長型，A組抗逆性較強。每個群體經過2代以上選擇并自繁純化，于2年后進行親本群體內自交和群體間雜交，親本產卵孵化后進行苗種培育，獲得4個不同交配組合的F1。為消除不同環(huán)境對凡納濱對蝦生長的影響，對蝦幼體在培育至體長2～3 cm后，轉入同一水泥池，分不同網箱養(yǎng)成，對凡納濱對蝦形態(tài)性狀進行對比測試。將大小相當的普通苗作為對照，普通苗為未經選育的本地苗，俗稱“土苗”，進行為期200 d的養(yǎng)殖，對12個形態(tài)和質量指標進行對比分析。試驗用凡納濱對蝦的來源和樣本信息見表1。

1.2 形質指標的測量

使用游標卡尺（精度為0.05 mm）和電子秤（精度為0.01 g）隨機抽樣，不同組合F1和普通苗組均隨機抽取100尾樣本，測定各組形態(tài)性狀和質量性狀，包括體長（BL）、頭胸甲高（CH）、頭胸甲長（CL）、頭胸甲寬（CW）、第一腹節(jié)長（FSL）、第二腹節(jié)長（SSL）、第六腹節(jié)長（SISL）以及體重（BW）、頭重（HW）和腹重（AW），共10個形質指標。凡納濱對蝦形態(tài)測量示意圖見圖1。

1.3 數據處理

1.3.1 數據標準化。

為消除凡納濱對蝦個體大小對形態(tài)學參數的影響，需要將每尾對蝦的所有參數做如下處理，所測7個形態(tài)指標（BL、CH、CL、CW、FSL、SSL、SISL）通過以下公式進行轉化［15-17］，分別得出7個形態(tài)特征指標（以X1～X7表示）。

Ms=Mo（Ls/Lo）b（1）

b=log10Mo/log10Lo（2）

式中，Ms為各性狀標準化數據；Mo為各性狀測量數據；Ls為頭胸甲長測量值平均值；Lo為頭胸甲長實際測量值。

對3個質量性狀參數通過公式計算進行轉化，得出5個質量評價指標（以X8～X12表示），分別為X8（鮮出肉率=腹重/體重×100%）、X9［肥滿度指數Ⅰ= 體重/（體長）3×100%］、X10［肥滿度指數Ⅱ=肉重/（體長）3×100%］、X11（重長指數Ⅰ=體重/體長）、X12（重長指數Ⅱ=肉重/體長）。

1.3.2 數據統計分析。

采用Excel 2013和SPSS 25.0軟件統計凡納濱對蝦群體各形質評價指標的平均值和標準差，對12組標準化數據進行多元統計分析。其中，通過主成分分析計算出各主成分的得分、方差貢獻率和累計方差貢獻率，并根據貢獻較大的主成分得分，繪制主成分散布圖［18-19］。判別分析采用逐步判別分析方法，每步選擇Willks的λ統計量最小的變量進入判別函數［20-21］。聚類分析根據5個群體間的歐式距離進行聚類［22］。

2 結果與分析

2.1 方差分析

方差分析結果（表2～3）顯示，普通苗組與AA組除X2、X4、X9外的其余9個性狀均有顯著差異（P<0.05）；普通苗組與BB組除X9外的其余11個性狀均有顯著差異（P<0.05）；普通苗組與AB組除X9外的其余11個性狀均有顯著差異（P<0.05）；普通苗組與BA組除X7、X9、X11、X12外的其余8個性狀均有顯著差異（P<0.05）；AA與BB組除X1、X5、X11、X12外的其余8個性狀均有顯著差異（P<0.05）；BB、AB、BA組間除X2、X4、X6、X8～X10的6個性狀均有顯著差異（P<0.05）。

2.2 主成分分析 通過Bartlett球形檢驗和KMO適合度檢驗發(fā)現，形質評價指標相關系數矩陣與單位矩陣有顯著差異（P<0.05）且適合度尚可（KMO=0.785>0.700），表明該研究所涉及形質評價指標適合進行主成分分析。對所有樣本的12個形質評價指標進行主成分分析，共獲得4個主成分，其方差貢獻率分別為40.14%、25.34%、15.98%和10.87%，累計貢獻率為92.33%。對于第一主成分，X3、X4、X5、X6、X7的影響較大，其負荷值分別為0.678、0.607、0.968、0.925和0.907；對于第二主成分，X9、X10、X11、X12的影響較大，其負荷值分別為0.783、0.831、0.747和0.797；對于第三主成分，X2、X4、X8的影響較大，其負荷值分別為0.691、0.654、0.632；對于第四個主成分，X8、X9、X3的影響較大，其負荷值分別為-0.597、0.537和0.415（表4）。不同對蝦群體第一、第二主成分散布圖見圖2。從圖2可以看出，除了普通苗組和AA組可區(qū)分外，BB、AB和BA組主成分落點相互重疊，不能清晰區(qū)分開來。

2.3 判別分析

通過判別分析共獲得4個典型判別函數，前3個典型判別函數的累積方差為99.3%（表5）。測量性狀與判別函數的相關系數見表6，從12個特征性狀中篩選出有顯著貢獻的7個變量，即X2、X4、X5、X7、X8、X11、X12。第一和第二判別函數散布圖見圖3。從圖3可以看出，普通苗組與各F1組區(qū)分明顯。在不同交配組合的F1各組中，BA組除了與AA組存在部分重疊外，與BB組、AB組區(qū)分明顯，AA組、BB組與AB組存在較多重疊。經過判別分析（表7）發(fā)現，普通苗組、BB、AB、BA各組的正確判別率較高，分別為96%、75%、84%和97%。AA組的錯判率較高，其判別準確率僅13%。判別結果進一步表明，普通苗組與F1代各組在形質特征性狀上存在顯著差異，各選育群體自交和雜交F1代各組間也存在差異，其中雜交組差異較為明顯。

2.4 聚類分析

對所有樣本的12個測量性狀進行聚類分析，得出5個凡納濱對蝦群體的聚類圖（圖4）。從圖4可以看出，普通苗組與F1代各組明顯聚為2支，F1代各組中自交組AA和BB親緣關系最近，雜交組中正交組（AB）與自交組親緣關系較近，反交組（BA）與各組親緣關系較遠。這進一步驗證了判別分析的結果，說明凡納濱對蝦通過雜交和自交選育的后代與普通苗產生了較大的形質特征差異，形成了具有單獨特征的新種群。

3 討論

在研究對蝦不同群體特征性狀時，樣本量是非常重要的前提條件。若樣本量太少，則不能充分反映群體的變異特征，從而會錯誤判斷不同群體間的差異。據文獻顯示，在進行主成分分析和判別分析時，樣本量與所測量性狀數的比值必須大于3［23］。該研究中不同對蝦群體的樣本量均為100尾，該比值均超過8，說明該研究結果具備統計學意義。

3.1 不同群體凡納濱對蝦形質特征方差分析

方差分析表明，經過項目組選育的凡納濱對蝦自交和雜交F1各群體與普通苗在12個指標性狀上存在不同程度的顯著差異，F1各組的平均鮮出肉率比普通苗組提高了6.12%，普通苗組表現出頭胸甲較細長的特征。F1各組體長的變異系數均明顯低于普通苗組，表現出選育后代更為整齊的特征。F1各組存活率明顯高于普通苗組，多代選育后具有對水泥池養(yǎng)殖環(huán)境更好的適應性。

3.2 不同群體凡納濱對蝦雜種優(yōu)勢分析

凡納濱對蝦屬于外來蝦種，我國海域沒有天然資源。近年來，國內科研機構和企業(yè)已培育出多個凡納濱對蝦新品種，且培育方法多以雜交育種為主。2016年通過審定的“廣泰1號”凡納濱對蝦新品種采用品系繁育結合正反交測試選育方法，在產業(yè)上得到較好應用。該研究結果進一步表明，雜交組合的子一代在各生長性狀上表現出一定的雜種優(yōu)勢。雜交組合中正交和反交子一代的性狀存在一定差異。反交BA組在各性狀中除了X8、X9和X10與正交AB組沒有顯著差異外，其他各性狀均顯著高于正交AB組，這種正反交雜種優(yōu)勢存在差異的現象也普遍存在于其他水產動物中。陳炳霖等［24］對不同種群羅非魚正反交子代奧尼羅非魚的生長性狀進行了對比分析，發(fā)現正交奧尼羅非魚群體全長、體高和尾柄高對體質量的影響較大，而反交群體全長和體高對體質量的影響較大。李永等［25］對2個地理種群斑節(jié)對蝦雜交子代的生長性能進行研究，結果表明正交組合的雜種優(yōu)勢均高于反交組合。張躍環(huán)等［26］對不同地理群體牡蠣開展種間雜交研究，結果表明正交組生長和存活性狀表現出積極的單親雜種優(yōu)勢，而反交組則表現出明顯的遠交衰退現象，變態(tài)率極低。該研究中正反交F1代在生長性狀上均表現出一定的母系遺傳效應［27］。

3.3 不同群體凡納濱對蝦形質特征差異的原因分析

在判別分析中，AA組判別準確率較低，僅13%。該研究結果表明，AA組中有31%被判別為BB組、19%被判別為AB組、37%被判別為BA組，說明BB、AB、BA組均與AA組在形態(tài)上存在不同程度的相似性，難以正確區(qū)分AA組。根據A組和B組選育群體的來源，分析原因可能是因為來自海南的選育群體B組最初也來源于泰國的選育系，二者之間可能存在一定的基因交流或混雜現象。泰國與我國海南的地理位置相近，貿易往來頻繁。海南地區(qū)常引進泰國的選育系作為對蝦種蝦，產出子一代銷往全國，也就是俗稱的“一代苗”，因此這種基因交流的可能性較大［28-29］。判別結果中BA組和普通苗組的判別準確率較高，說明其與其他各組在形質特征上已存在明顯差異。其中，普通苗組為未經選育但多代生長于本地的群體，已逐漸形成適應本地環(huán)境的新種群。從不同群體歐氏距離聚類圖也可以看出，普通苗組與選育群體子代組的親緣關系已較遠。BA組通過選育群體雜交而獲得，其在形質特征上與親本存在一定差異，表現出一定的雜種優(yōu)勢，且親緣關系較正交子代更遠。

凡納濱對蝦的生長周期一般為1～2年，商品蝦通常養(yǎng)殖90～150 d即可出塘。該研究過程中可能由于養(yǎng)殖網箱的限制，試驗對蝦生長普遍偏慢，生長數據低于正常養(yǎng)殖水平，因此項目組將最終測定時間延長至200 d，但由于生長較慢，選育群體的生長優(yōu)勢未能充分體現出來，與普通苗的大小差異沒有達到預期效果。今后將改進設施條件，開展進一步研究。

參考文獻

［1］ 農業(yè)農村部漁業(yè)漁政管理局，全國水產技術推廣總站，中國水產學會.2021年中國漁業(yè)統計年鑒［M］.北京：中國農業(yè)出版社，2021：26-34.

［2］ 代平，孔杰，欒生.我國凡納濱對蝦種質資源引進與分析［J］.科學養(yǎng)魚，2018（1）：3-5.

［3］ 趙永貞.凡納濱對蝦 “桂海 1號”［J］.海洋與漁業(yè)，2013（3）：54.

［4］ 孔杰，何建國，江謝武，等.凡納濱對蝦 “海興農 2 號”［J］.中國水產，2017（11）：67-71.

［5］ 葉寧.凡納濱對蝦新品種 “興海 1 號” 的選育技術［J］.海洋與漁業(yè)，2018（7）：70-71.

［6］ 黃小帥，徐煜，胡曉娟，等.利用微衛(wèi)星標記分析7個凡納濱對蝦引進群體子一代的遺傳多樣性［J］.南方水產科學，2019，15（1）：54-62.

［7］ 彭敏，陳慧芳，李強勇，等.凡納濱對蝦連續(xù)3個世代選育群體的遺傳多樣性分析［J］.南方農業(yè)學報，2020，51（6）：1442-1450.

［8］ 方振朋，孟憲紅，李旭鵬，等.基于微衛(wèi)星分子標記的凡納濱對蝦商業(yè)苗種遺傳多樣性分析［J］.漁業(yè)科學進展，2020，41（5）：101-109.

［9］ 申淑慧，戴習林.基于生長和抗逆功能基因SNP分子標記的凡納濱對蝦野生及選育群體遺傳多樣性分析［J］.南方農業(yè)學報，2020，51（11）：2836-2845.

［10］ 劉峻宇，劉均輝，孔杰，等.凡納濱對蝦SNP標記開發(fā)與家系親緣關系驗證分析［J］.漁業(yè)科學進展，2021，42（1）：108-116.

［11］ 李振通，成美玲，田永勝，等.雜交種“云龍石斑魚”與親本的形態(tài)差異分析［J］.漁業(yè)科學進展，2019，40（4）：73-83.

［12］ 唐永成，彭姣，趙運林，等.池養(yǎng)中華絨螯蟹不同性別形態(tài)及質量差異分析［J］.漁業(yè)科學進展，2019，40（6）：114-120.

［13］ 滕爽爽，林興管，方軍，等.縊蟶不同種群數量形態(tài)特征的比較分析［J］.海洋湖沼通報，2020（4）：127-137.

［14］ 黃光華，雷燕，馬華威，等.羅氏沼蝦不同地理種群親本形態(tài)及繁育性能比較［J］.水產科技情報，2019，46（4）：187-191.

［15］ ANASTASIADOU C，LIASKO R，LEONARDOS I D.Biometric analysis of lacustrine and riverine populations of Palemonetes antennarius （H.Milne-Edwards，1837） （Crustacea，Decapoda，Palaemonidae） from north-western Greece［J］.Limnologica，2009，39（3）：244-254.

［16］ HURLBUT T，CLAY D.Morphometric and meristic differences between shallow-and deep-water populations of white hake （Urophycis tenuis） in the southern Gulf of St Lawrence［J］.Canadian journal of fisheries and aquatic sciences，1998，55（10）：2274-2282.

［17］ FERRITO V，MANNINO M C，PAPPALARDO A M，et a1.Morphological variation among populations of Aphanius fasciatus Nardo，1827 （Teleostei，Cyprinodontidae）from the Mediterranean［J］.Fish biology，2007，70（1）：1-20.

［18］ SMITH L I.A tutorial on principal components analysis［M］.New York：Cornell University，2002：51-52.

［19］ HINTON P R，MCMURRAY I，BROWNLOW C.SPSS explained［M］.London：Routledge，2014：21-22.

［20］ CACOULLOS T.Discriminant analysis and applications［M］.New York：Academic Press，2014：1-2.

［21］ 李秀啟，安麗，王亞楠，等.南水北調東線南四湖湖鱭種群的形態(tài)判別與地理親緣關系［J］.水生態(tài)學雜志，2015，36（2）：81-87.

［22］ 余來寧，夏小平，楊東，等.草魚♀ × 鳡♂雜交F1及其親本的形態(tài)差異分析［J］.安徽農業(yè)科學，2014，42（35）：12483-12486.

［23］ KOCOVSKY P M，ADAMS J V，BRONTE C R.The effect of sample size on the stability of principal components analysis of truss-based fish morphometrics［J］.Transactions of the American fisheries society，2009，138（3）：487-496.

［24］ 陳炳霖，肖煒，鄒芝英，等.正反交奧尼羅非魚形態(tài)性狀差異及對體質量的影響［J］.漁業(yè)科學進展，2020，41（3）：94-102.

［25］ 李永，黃建華，楊其彬，等.斑節(jié)對蝦2個地理種群自交與雜交F1的生長性能［J］.中國水產科學，2012，19（5）：784-789.

［26］ 張躍環(huán)，王昭萍，閆喜武，等.太平洋牡蠣與近江牡蠣的種間雜交［J］.水產學報，2012，36（8）：1215-1224.

［27］ HEATH D D，FOX C W，HEATH J W.Maternal effects on offspring size：Variation through early development of chinook salmon［J］.Evolution，1999，53（5）：1605-1611.

［28］ 周勁松，曹哲明，楊國梁，等.羅氏沼蝦緬甸引進種和浙江本地種及其雜交種的生長性能與SRAP分析［J］.中國水產科學，2006，13（4）：667-673.

［29］ 林紅軍，沈琪，張呂平，等.凡納濱對蝦生長性狀的雙列雜交分析［J］.熱帶海洋學報，2010，29（6）：51-56.