許藝蘭，陳 忠，覃俊奇，文露婷，潘賢輝，周康奇，黃 姻，林 勇，杜雪松**

(1.廣西壯族自治區(qū)水產科學研究院,廣西水產遺傳育種與健康養(yǎng)殖重點實驗室，廣西南寧 530021；2.廣西壯族自治區(qū)獸醫(yī)研究所,廣西獸醫(yī)生物技術重點實驗室，廣西南寧 530001)

全州禾花鯉(Cyprinuscarpiovar.Quanzhouensis)又名禾花魚、禾花烏鯉，因長期在稻田養(yǎng)殖、食禾花而得名，分布于廣西北部的全州縣及其周邊地區(qū)。與其他鯉魚養(yǎng)殖品種相比，其典型特征是體型粗短、體背部烏褐色、鰓蓋和腹部呈半透明，是一個極具廣西特色的稻漁綜合種養(yǎng)優(yōu)良品種[1]，但其生長速度慢的缺點也較為明顯，早期測定顯示1齡魚平均體質量為4.0-13.5 g[2]。近年來，隨著稻田養(yǎng)殖產業(yè)的迅速發(fā)展，各地陸續(xù)引進生長性狀更優(yōu)的鯉魚品種開展稻田養(yǎng)殖，有意或無意的雜交導致本地禾花鯉種群逐漸出現種質資源混亂、特色性狀退化等問題[3]，不利于全州禾花鯉種質資源的長期保存和利用。因此，對全州禾花鯉原種群開展生長速度改良選育是目前亟待開展的研究工作之一。

遺傳參數估計作為育種選育的基本參考數據，遺傳參數和育種值的評估對育種方案的制定及優(yōu)化具有重要的指導意義。自1974年Henderson[4]提出最佳線性無偏預測(Best Linear Unbiased Prediction,BLUP)法后，統(tǒng)計程序逐漸發(fā)展起來，動物模型開始被應用于遺傳參數估計的研究工作中。這種方法可充分利用系譜等信息，全面、快速地評估群體的遺傳參數和遺傳進展，有利于加快遺傳選育效率，減少成本消耗。因此，動物模型BLUP法在國內外已成為水產動物選擇育種的主要評估方法，在對蝦(Penaeusorientalis)[5]、牙鲆(Paralichthysolivaceus)[6]、鳙(Hypophthalmichthysnobilis)[7]、長牡蠣(Crassostreagigas)[8]、鯉魚(Cyprinuscarpio)[9]、尼羅羅非魚(Oreochromisniloticus)[10]等的報道中均獲得了較為理想的選育效果。本研究采用全州禾花鯉性成熟個體為親本，通過人工授精的方法構建家系，采用BLUP法首次對全州禾花鯉體質量、體長、體寬、日增重以及代表其體型粗短特點的體長/體寬等遺傳參數和育種值進行分析，比對依據育種值選擇和依據表型值選擇的差異性，以期為全州禾花鯉經濟性狀的遺傳改良提供參考。

1 材料與方法

1.1 全州禾花鯉親本

全州禾花鯉親本為廣西桂林市綠淼生態(tài)農業(yè)有限公司在全州縣才灣鎮(zhèn)所收集的養(yǎng)殖群體，在廣西農業(yè)良種海南南繁育種基地養(yǎng)殖至性成熟。配對所用父本平均體質量(422.6±47.8) g，母本平均體質量(524.6±36.8) g。

1.2 家系構建與魚苗培育

選擇性腺發(fā)育成熟的全州禾花鯉40尾(父本20尾、母本20尾)，采用人工催產授精、嵌套交配的方法，利用9尾父本和8尾母本構建全同胞家系23個(表1)。催產藥物采用促黃體素釋放激素類似物(LRH-A)，使用計量為母本10 μg/kg、父本5 μg/kg，每個家系各挑取適量的受精卵，置于孵化箱內孵化。待仔魚孵化后，每個家系隨機選取1 000尾魚苗放入網箱中養(yǎng)殖。

1.3 家系培育與性狀測量

待網箱中魚苗生長至6 g以上，每個家系隨機選取40尾，采用PIT魚類芯片(1.35 mm×8 mm)標記后，在3口面積為200 m2水池中混合養(yǎng)殖。待全州禾花鯉養(yǎng)殖180日齡左右時，捕獲池塘中70%以上的個體，測量體質量、體長、體寬，計算日增重和體長/體寬。測量時，將其放入濃度為100 mg/L的MS-222麻醉劑中快速麻醉，用電子天平測定體質量，精確度為0.1 g。采用游標卡尺測量體長和體寬，精確度為0.1 mm。

表1 全州禾花鯉家系構建

1.4 數據處理與分析

1.4.1 描述性統(tǒng)計分析

測量數據采用SPSS 26.0軟件進行描述性統(tǒng)計分析和正態(tài)性檢驗，采用R語言Hmisc包進行性狀間相關性分析和顯著性檢驗。

1.4.2 遺傳參數估計

采用ASReml-R軟件評估全州禾花鯉各個生長性狀的方差組分，選擇約束極大似然(Restricted Maximum Likelihood，REML)法和最佳線性無偏預測(Best Linear Unbiased Prediction，BLUP)法計算遺傳參數和育種值。其中，體質量、體長、體寬及體長/體寬采用的分析模型[11]如下：

yijklmn=μ+Pondi+Agej+Sexk+Pondi×

Agej+Agej×Sexk+Pondi×Sexk+al+fm+eijklmn，

(1)

式中：yijklmn表示個體的性狀測定值，μ表示總體平均值，Pondi表示養(yǎng)殖池效應，Agej表示養(yǎng)殖日齡效益，Sexk表示性別效應，Pondi×Agej、Agej×Sexk、Pondi×Sexk分別表示池塘、日齡、性別之間的互作效應，al是加性遺傳效應，fm是家系在混合養(yǎng)殖之前單獨飼養(yǎng)時的隨機效應，eijklmn是個體隨機殘差。

日增重采用的遺傳分析模型[12]如下：

yijklm=μ+Pondi+Sexj+Pondi×Sexj+ak+

fl+eijklm，

(2)

式中：Sexj表示性別效應，Pondi×Sexj表示池塘與性別之間的互作效應，ak是加性遺傳效應，fl是家系在混合養(yǎng)殖之前單獨飼養(yǎng)時的隨機效應。

2 結果與分析

2.1 全州禾花鯉家系的生長性狀及相關性

全州禾花鯉的體質量、體長、體寬、體長/體寬、日增重的描述性統(tǒng)計結果如表2。體質量和日增重指標變異系數較高，分別為0.52和0.51，表明全州禾花鯉種群體質量性狀分化程度較大，個體間生長速度存在較大差異。體寬性狀變異系數僅有0.10，顯示出較高的一致性。由于體型粗短是全州禾花鯉的典型特征之一，所以本研究采用體長/體寬衡量該表型，參數越小說明體型粗短特征越顯著。統(tǒng)計顯示，其平均值為5.53，變異系數為0.18，在所有性狀中處于中等水平。正態(tài)性檢驗體質量、體長、體寬、體長/體寬、日增重性狀均不符合標準正態(tài)分布，其中體寬(偏度0.31)、體長/體寬(偏度0.27)最接近標準正態(tài)分布，可能是由于養(yǎng)殖條件導致的種內競爭和樣本量過小。性狀相關性分析顯示，日增重與體長相關性最大(R2=0.96，P<0.01)，其次為日增重與體寬(R2=0.89，P<0.01)，均達到極顯著水平，日增重與體長/體寬相關性最小(R2=0.26，P>0.05)(圖1)。

表2 全州禾花鯉生長性狀的描述性統(tǒng)計

** means extreme significance (P<0.01)，* means significance (P<0.05)

2.2 方差分量和遺傳力

利用ASReml-R軟件統(tǒng)計分析測量的表型數據，采用單性動物模型獲得了全州禾花鯉體質量、體長、體寬、體長/體寬、日增重的方差組分和遺傳力(表3)。結果表明體質量、體長、體寬和日增重為高遺傳力，分別是0.45(P<0.01)、0.42(P<0.01)、0.36(P<0.01)和0.39(P<0.05)，體長/體寬為低遺傳力(0.12，P>0.05)，說明全州禾花鯉體型粗短性狀在選育中更難獲得遺傳進展。

表3 全州禾花鯉生長性狀的方差分量和遺傳力

2.3 表型相關和遺傳相關

采用兩性動物模型估計性狀之間的表型相關與遺傳相關(表4)，結果表明體質量、體寬、體長、日增重性狀之間均有較高的表型相關與遺傳相關(0.76-0.99)，體質量與體長/體寬、日增重與體長/體寬之間的表型相關與遺傳相關較低，分別是0.62，0.67和0.55，0.58。表明全州禾花鯉體型粗短特征和生長性狀的關聯性較弱，在選育過程中僅開展針對生長性狀的選育很難獲得在體型方面的遺傳進展，應采用多性狀復合選育策略，以達到保障全州禾花鯉體型特色的同時提高其生長速度。

表4 利用動物模型估計的全州禾花鯉表型相關和遺傳相關

2.4 個體育種值與個體選擇比較

對全州禾花鯉各性狀育種值進行相關性分析發(fā)現，日增重與體長、日增重與體寬、體長與體寬育種值之間相關性極為顯著(P<0.01)，其他性狀之間相關性較弱，如體質量與體寬之間的相關系數僅為0.13(P<0.01)(圖2)。另外，根據性狀表型值和育種值分別對各性狀按10%的留種率進行留種選擇，統(tǒng)計2種留種方法入選個體的家系數量、平均表型值、平均育種值以及同時入選的個體相同率，結果表明按照表型值選擇前10%的個體來自18-23個家系，按照育種值選擇的前10%個體來自14-18個家系。從個體水平分析，2種留種方法相比，上述5種性狀排在前10%的個體相同率分別為23.53%、29.41%、44.12%、55.88%和63.50%(表5)。可見采用表型值選擇和采用育種值選擇，在選留家系和個體層面均存在明顯的差異。

** means extreme significance (P<0.01),* means significance (P<0.05)

表5 基于表型值和育種值全州禾花鯉個體的選擇比較

3 討論

3.1 全州禾花鯉生長性狀的特點

生長性狀是魚類最重要的經濟性狀之一，提高魚類生長速度可縮短養(yǎng)殖周期，有利于節(jié)約養(yǎng)殖空間和降低養(yǎng)殖成本。全州禾花鯉作為廣西特色水產品種，具有高品質的同時也存在生長速度慢的突出缺點。本研究測定6月齡全州禾花鯉的平均體質量為(27.38±14.01) g，平均日增重為0.16 g/d，遠低于松浦鏡鯉(C.carpiovar.Songpu,1.98 g/d)、建鯉(C.carpiovar.Jian,1.23 g/d)的生長速度[13]。因此，提高全州禾花鯉的生長速度是今后種質改良的主要目標。在體質量、體長、體寬、體長/體寬、日增重5個性狀中，變異系數分別為0.52，0.17，0.10，0.18和0.51，其中體質量的變異系數最大。在魚類中，體質量性狀變異系數較大且較為常見，如盧薛等[14]發(fā)現，不同日齡的鱖魚(Sinipercachuatsi)體質量、體長、體高3個生長性狀中，體質量性狀的變異系數為0.43-0.55，與本研究測定的結果接近。李祥孔等[15]對不同時期牙鲆的生長性狀分析變異系數，結果均屬體質量性狀的變異系數(0.33-0.38)最大，類似的研究結果也出現在銀大麻哈魚(Oncorhynchuskisutch)、甌江彩鯉(C.carpiovar.Color)等物種中[16,17]。但與目前大規(guī)模養(yǎng)殖的建鯉[18]等相比，全州禾花鯉體質量變異系數較高，個體差異大，性狀遺傳不穩(wěn)定，這不利于全州禾花鯉規(guī)模養(yǎng)殖效益的提高，也是今后選育工作需要重點關注的特征，另一方面也可能說明該群體遺傳多樣性較高，性狀改良潛力大。

3.2 遺傳力、表型相關及遺傳相關

為了提高選育效率，參考目標性狀的遺傳參數以制定適宜的選育策略是常用方法之一。本研究采用動物模型對全州禾花鯉生長性狀遺傳力估計發(fā)現，體質量和體長性狀為高遺傳力(h2>0.4)，體寬和日增重性狀為中等遺傳力(0.2

遺傳相關和表型相關是描述不同性狀關聯性的重要參數，利用性狀關聯性進行間接選育，對于改良度量難度較大的性狀極為有用[30]。在選擇育種中，有時對目標性狀的選擇出現一定難度時，通過遺傳相關來確定某一性狀與目標性狀的相關程度，可間接對目標性狀進行相關選擇[31]。如韋信鍵等[32]研究大黃魚(Larimichthyscrocea)早期不同階段生長性狀遺傳參數時發(fā)現，1月齡大黃魚體質量與全長遺傳相關為0.83-0.90，6月齡時則為0.97-0.98，并提出可以利用6月齡時2個性狀高度的遺傳相關進行性狀間的關聯選擇。在本研究的全州禾花鯉5個性狀的遺傳相關和表型相關分析結果中，體質量、體寬、體長、日增重性狀之間具有顯著相關性，其中日增重與體長性狀之間的遺傳相關(0.98)和表型相關(0.96)最高，表明可利用體長性狀對日增重進行間接選育。

3.3 育種值與個體選擇

育種值作為個體性狀中可遺傳的部分，依靠育種值進行優(yōu)良選育相對于依靠表型值選育更穩(wěn)妥、更有效率[33]。鄭靜靜等[34]對60日齡、105日齡、150日齡不同生長階段的凡納濱對蝦(Litopenaeusvannamei)選取體質量性狀的育種值和表型值前10%個體計算平均育種值，前者比后者分別高出約20.06%、26.03%、12.01%，說明依靠育種值選擇比表型值選擇更有優(yōu)勢。本研究中，采用表型值和育種值分別對各性狀按10%的留種率進行留種選擇，分析發(fā)現采用育種值篩選后，保留個體所在的家系數量比表型值選擇減少，被兩種方法同時選留的個體占所有保留個體的23.53%-63.50%，表明在全州禾花鯉選育中采用育種值篩選方法的有效性。全州禾花鯉是地方性特色品種，除生長性狀之外，其體型和體色也是決定其特色和經濟價值的重要性狀。綜上所述，對全州禾花鯉的性狀改良應采取多性狀綜合選育策略，以生長速度、體型和體色性狀的同步改進和提升為選育目標開展家系選育，從而獲得優(yōu)質高產的禾花鯉新品種。

4 結論

利用全州禾花鯉家系的生長性狀指標對其遺傳參數進行估計，主要結論包括：(1)全州禾花鯉體質量性狀變異系數高，規(guī)格差異大。(2)全州禾花鯉生長性狀(體質量、體長、體寬和日增重)為高遺傳力性狀，體型性狀(體長/體寬)為低遺傳力性狀。(3)體質量、體寬、全長、日增重性狀相互之間有較高的表型相關與遺傳相關，體質量與體長/體寬、日增重與體長/體寬之間表型相關與遺傳相關較低。(4)采用育種值開展個體選擇與表型值選擇，結果存在明顯差異。(5)全州禾花鯉種質改良應以獲得優(yōu)質高產禾花鯉新品種為目標，采用多性狀綜合選育策略對生長速度、體型和體色性狀進行同步選擇。