朱 吉 ，劉 建 ，孫建幫 ，楊仕柳 ，李靜如 ，謝菊蘭 ，譚建坤 ，彭英林 ，※

（1.湖南省畜牧獸醫(yī)研究所，湖南長沙 410131；2.湘村高科農(nóng)業(yè)股份有限公司，湖南婁底417000；3.湖南省畜牧水產(chǎn)局，湖南長沙 410006）

隨著育種技術的進步和發(fā)展，豬的育種目標不再單純注重降低背膘厚、提高生長率，逐步發(fā)展為提高母豬的窩產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)等繁殖性狀，因此探索豬繁殖性狀中各種效應的作用和規(guī)律，并將其運用到養(yǎng)豬育種實踐中顯得越來越重要。

育種方案的制定必須以準確、可靠的遺傳參數(shù)估計為前提條件。國內(nèi)外的大量研究證明，MTDFRE mL方法是目前動物育種中遺傳參數(shù)估計較理想的方法，且得到了廣泛應用。本研究旨在用MTDFRmL方法估計湘村黑豬某些繁殖性狀的遺傳力、母體效應和遺傳相關，為選擇提供必要的遺傳參數(shù)[1～5]。

1 材料與方法

1.1 材料

本文所用資料取自2005～2010年湘村高科農(nóng)業(yè)發(fā)展有限公司黑豬原種場進行的湘村黑豬選育群育種資料。本文分析的性狀主要有：產(chǎn)仔數(shù)(TNB)、產(chǎn)活仔數(shù)(NBA)、初生窩重(LWB)和 21日齡窩重（LW21)等繁殖性狀，總記錄數(shù)為823胎次。

1.2 數(shù)據(jù)整理

用EXCEL對系譜資料和數(shù)據(jù)資料進行整理,然后利用SAS統(tǒng)計分析比較固定效應不同水平間的差異顯著性，然后將差異不顯著的水平合并從而減少水平數(shù)，根據(jù)結(jié)果重新合并固定效應、確定水平；隨機效應中，母豬個體號計算個體的動物效應，系譜中的母豬號用于計算母體遺傳效應、永久環(huán)境效應，母豬出生所在窩號計算窩效應。

1.3 混和模型：

Y—觀察值向量；b—固定效應向量（世代、季節(jié)、胎次）；u—隨機動物效應向量；

m—隨機母體遺傳效應向量；pe—影響同一母豬永久環(huán)境效應向量；t—母豬出生所在窩的窩效應向量；e—隨機殘差效應向量；X、Z、W、S1、S2－分別對應固定效應、隨機效應、母體效應、永久環(huán)境效應、窩效應的結(jié)構(gòu)矩陣。

2 結(jié)果與分析

2.1 表型參數(shù)

所分析的性狀有效樣本數(shù)、平均數(shù)和標準差見表1。由表1可見，經(jīng)過五個世代的選育，湘村黑豬的各項繁殖性狀基本上處于穩(wěn)步上升的趨勢，表明選育方法是適當?shù)摹?010年與2005年比較，產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)、初生窩重、21日齡窩重分別提高了 0.03頭、0.94頭、1.26 kg、5.51 kg。

?

?

2.2 繁殖性狀的遺傳力

表2列出了湘村黑豬各繁殖性狀的遺傳力、母體效應率、窩效應率估計值。繁殖性狀估計遺傳力的變化范圍為：0.07～0.19，TNB(頭)、NBA(頭)、LWB（kg)、LW21(kg)的遺傳力估計值分別是 0.19、0.15、0.07、0.11。母體效應率、窩效應率的變化范圍分別是 0.07～0.11、0.14～0.34。

2.3 繁殖性狀間的相關分析

表3可見各繁殖性狀間的相關分析結(jié)果。兩兩性狀間遺傳相關、永久環(huán)境相關、殘差相關、表型相關范圍分別為 0.40～0.91、0.68～0.97、0.33～0.62、0.31～0.87，遺傳相關和表型相關均表現(xiàn)為TN B/NBA的最高（0.91，0.87），而 TN B/LW21最低（0.40，0.31），總的來說，性狀間相關達到了中等偏上水平。

?

3 討論

3.1 繁殖性狀遺傳力

繁殖性狀的遺傳力一般都很低。因為窩產(chǎn)仔數(shù)、窩重等窩性狀由公豬、母豬及胚胎的基因型等因素互作決定。本試驗結(jié)果，TNB、NBA、LWB、LW21的遺傳力估值分別為0.19、0.15、0.07、0.11，略高于國外部分學者的報道結(jié)果[6-10]。分析差異產(chǎn)生的原因可能如下：①國內(nèi)一般樣本含量小，取樣誤差大，而國外樣本含量大，取樣誤差??；②采用方法的不同；③估計遺傳參數(shù)時，模型的選擇和設計不同。

3.2 母體效應率和窩效應率

許多性狀受母體環(huán)境影響，特別是對那些哺育早期度量的性狀。當母性效應被忽視時，其它效應的方差組分就有可能被過高地估計。但目前對母性效應影響的估計結(jié)果相差較大，Preze[11]等結(jié)果發(fā)現(xiàn)母性效應不顯著或沒有，而Ferraz[12]等發(fā)現(xiàn)母性方差占表型方差的比例為0.01～0.13。本研究四個性狀的母體效應率為0.08～0.10、窩效應率的變化范圍分別為0.10～0.35，表明母體遺傳效應、窩效應均占有一定比例，對繁殖性狀產(chǎn)生一定的影響。

3.3 性狀相關

兩兩性狀間遺傳相關、永久環(huán)境相關、殘差相關、表型相關范圍分別為0.40～0.91、0.68～0.97、0.33～0.62、0.31～0.87，基本處于中等偏上的正相關，符合一般規(guī)律。R．Roehe和B．W．Kennedy[13]報道了Yorkshire TNB與NBA的遺傳相關為0.993，張偉峰[14]等報道北京花豬TNB和NBA的遺傳相關為0.91，LBW與TNB、NBA的遺傳相關為0.71、0.75，與本試驗的結(jié)果基本相近。 □

[1] 張 沅，等.畜禽育種中的線型模型[M]，北京：北京農(nóng)業(yè)大學出版社，1993.

[2] 張 勤.家畜育種值和遺傳參數(shù)估計方法的發(fā)展及現(xiàn)狀[J]，草食家畜,1990,(6)：1-4，46.

[3] Monchai等，動物育種中解決混合模型復雜計算問題的軟件應用[J]，華南農(nóng)業(yè)大學學報，2005，26(S1)：83-78.

[4] 張勝利.遺傳參數(shù)不同估計方法的比較[J]，遺傳，1995，17(6)：21-23.

[5] 鄭友民.用MTDFREML方法估計大白豬遺傳參數(shù)[J]，養(yǎng)豬，1999，（1）：32-33.

[6] Cho.Y C，et al,Estimation of genetic parameters for major economic traits in swine[J],Korean journal of animal sicence,1998,40(5)：447-454.

[7] Bidanel J P,et al,Use of a BLUP animal model for gnetic evaluation of French Large white and French Landrace Pig for Prolificacy[J],Pig news and information,1998,18(l)：4.

[8] Estany J,et al.Estimation of genetic Parameters for liter size in Danish Landrace and yorkshire Pigs[J],Anim.sci,1995,60(2)：315-324.

[9] Irgang R.,et al.,Genetic parameters for liter size of different paritied in Durco,Landrace and Large white sows[J].J of Anim.sci,1994,72(9)：2237-2246.

[10] Rosa V C．Revista-Brasileria de genetica[J]．1995，18(2)：199-205.

[11] Misztel I,Perez-Enciso M.Sparse matrix inversion for REML estimation of variance components by exepctation maximization[J].J.diary Sci.1993,76(5)：1479-1483.

[12] FerrnaJ B S,et al,Animal model estimation of genetic parameters and response to selcetion for litter size and weight，gorwth and backfatin c losed seedstock porductions of Large white and Landrace swine[J].J of Anim sci.1993,71(4)：850-858.

[13] Roehe R,et al,estimation of genetic paramaters for little size in canadian Yorkshire and Landrace swine with each parity of forrowing treated different trait[J]．J Anim Sci，1995，73：2959-2970.

[14] 張偉峰，等.應用動物模型REML估計豬的遺傳參數(shù)的研究[J]．養(yǎng)豬，1999，(3)：31-32.