利用平均信息最大似然法估計(jì)中國美利奴羊羊毛性狀遺傳參數(shù)

張艷花，于麗娟，蔣曉梅，瑪爾孜亞·亞森，阿米尼古麗·阿不來孜，吳偉偉，古麗努爾·馬哈提，瑪依肯·沙力，拉扎提·艾尼瓦爾，田可川

(1.新疆畜牧科學(xué)院，烏魯木齊 830011；2.中國農(nóng)業(yè)大學(xué)，北京 100193)

利用平均信息最大似然法估計(jì)中國美利奴羊羊毛性狀遺傳參數(shù)

張艷花1，2，于麗娟1，蔣曉梅1，瑪爾孜亞·亞森1，阿米尼古麗·阿不來孜1，吳偉偉1，古麗努爾·馬哈提1，瑪依肯·沙力1，拉扎提·艾尼瓦爾1，田可川1

(1.新疆畜牧科學(xué)院，烏魯木齊 830011；2.中國農(nóng)業(yè)大學(xué)，北京 100193)

【目的】估計(jì)鞏乃斯種羊場中國美利奴羊羊毛性狀遺傳力及性狀間的遺傳相關(guān)，為種羊選擇及育種方案的制定提供依據(jù)?！痉椒ā繙y定鞏乃斯種羊場1 351只中國美利奴羊周歲母羊羊毛自然長度、剪毛量、羊毛纖維直徑、羊毛纖維直徑變異系數(shù)等性狀，利用DMUv6軟件及SAS v8.2，對(duì)測定數(shù)據(jù)進(jìn)行分析，采用多性狀動(dòng)物模型和平均信息最大似然法(AIREML)分析各性狀的方差組分，估計(jì)遺傳參數(shù)?！窘Y(jié)果】研究估計(jì)出中國美利奴羊羊毛自然長度、羊毛纖維直徑、羊毛纖維直徑變異系數(shù)、剪毛量的遺傳力分別為0.39、0.52、0.30、0.37。羊毛纖維直徑及羊毛纖維直徑變異系數(shù)間遺傳相關(guān)為-0.73，羊毛自然長度與剪毛量、羊毛纖維直徑與羊毛自然長度之間遺傳相關(guān)為0.24、0.33?！窘Y(jié)論】研究估計(jì)出中國美利奴羊羊毛自然長度、剪毛量、羊毛纖維直徑、羊毛纖維直徑變異系數(shù)以及剪毛量的遺傳力，可為今后制定中國美利奴羊育種方案、種羊遺傳評(píng)估提供數(shù)據(jù)支持。

羊毛纖維直徑；羊毛纖維直徑變異系數(shù)；剪毛量；羊毛自然長度；遺傳參數(shù)；中國美利奴羊

0 引 言

【研究意義】中國美利奴羊是我國主要的細(xì)毛羊品種，從1972年起，由國家農(nóng)牧漁業(yè)部牽頭，組織新疆(包括兵團(tuán))、內(nèi)蒙古、吉林等三省(區(qū))四方聯(lián)合攻關(guān)，于1985年育成。中國美利奴羊育成以后，從鞏乃斯種羊場推廣到了國內(nèi)主要細(xì)毛羊產(chǎn)區(qū)，在我國細(xì)羊毛生產(chǎn)中發(fā)揮了重要作用[1]。此后的幾十年，我國的細(xì)毛羊育種工作者對(duì)中國美利奴羊不斷選育提高，使其各項(xiàng)生產(chǎn)性能都得到了提高。估計(jì)中國美利奴羊的毛用性狀遺傳力及遺傳相關(guān)，找到其遺傳規(guī)律，為提高中國美利奴羊的羊毛質(zhì)量及產(chǎn)量性狀的育種目標(biāo)提供數(shù)據(jù)支持?！厩叭搜芯窟M(jìn)展】遺傳參數(shù)估計(jì)是家畜育種中的一項(xiàng)基本任務(wù)，在進(jìn)行種羊選擇、制定育種規(guī)劃、預(yù)測選擇效果中具有重要作用[2]。近20年來，國外細(xì)毛羊育種工作者，在羊的遺傳參數(shù)估計(jì)方面做了大量的工作，主要采取的方法有混合模型約束最大似然法(REML，Restricted maximum likelihood)，還有部分采取了貝葉斯方法(Bayesian)，估計(jì)的性狀主要為產(chǎn)毛性狀、產(chǎn)肉性狀以及產(chǎn)奶性狀，其中羊毛性狀包括原毛量、凈毛重、羊毛纖維直徑、羊毛纖維直徑變異系數(shù)、羊毛纖維直徑標(biāo)準(zhǔn)差、羊毛自然長度、羊毛纖維強(qiáng)度、剪毛量以及卷曲率等[3-17]。我國育種工作者針對(duì)不同地方、不同品種的細(xì)毛羊估計(jì)了遺傳參數(shù)，商數(shù)岐等[18]估計(jì)了東北細(xì)毛羊羊毛自然長度、細(xì)度、密度、腹毛覆蓋度等性狀的遺傳參數(shù)；齊侃虎[19]估計(jì)了良種細(xì)毛羊羊毛細(xì)度、剪毛量的遺傳參數(shù)；李俊年等[20]估計(jì)了中國美利奴(新疆型)細(xì)毛羊羊毛自然長度、凈毛率、污毛重、羊毛纖維直徑等性狀的遺傳參數(shù)；李文輝等[21]估計(jì)了甘肅高山細(xì)毛羊斷奶及1.5歲時(shí)毛密度、羊毛自然長度、油汗評(píng)分、彎曲評(píng)分等性狀的遺傳參數(shù)；梅花[22]估計(jì)了敖漢細(xì)毛羊斷乳毛細(xì)度、周歲及2歲羊羊毛自然長度、細(xì)度、剪毛量等性狀的遺傳參數(shù)；黃錫霞等[23]估計(jì)了優(yōu)質(zhì)細(xì)毛羊羊毛自然長度、細(xì)度、剪毛量等性狀的遺傳參數(shù)；張亞軍等[24]估計(jì)了中國美利奴羊(新疆型)羊毛自然長度、羊毛纖維直徑、剪毛量等性狀的遺傳參數(shù)；李莉等[25]計(jì)算了敖漢細(xì)毛羊剪毛量、羊毛自然長度、羊毛纖維直徑的遺傳力；徐新明等[26]估計(jì)了中國美利奴羊初生重等性狀的遺傳參數(shù)；國內(nèi)研究者主要采取的方法有動(dòng)物模型非求導(dǎo)約束最大似然法(MTDFREML)、半同胞組內(nèi)相關(guān)法。在遺傳參數(shù)估計(jì)方法方面，1995年約翰遜等開發(fā)了平均信息最大似然法(AIREML， Average Information Restricted Maximum Likelihood )，估計(jì)多性狀遺傳參數(shù)及方差組分，它比EM算法(Expectation Maximization Algorithm，最大期望值算法)效率提高了2～3倍[27]?！颈狙芯壳腥朦c(diǎn)】中國美利奴羊育成后，30年來鞏乃斯種羊場的選育工作從未中斷過，各性狀的生產(chǎn)性能與育成初期相比已有明顯變化。隨著計(jì)算技術(shù)的發(fā)展，新的計(jì)算方法、計(jì)算軟件不斷出現(xiàn)，使遺傳參數(shù)估計(jì)的方法更加的高效、準(zhǔn)確。研究利用丹麥農(nóng)業(yè)科學(xué)院的Per Madsen與Just Jensen開發(fā)的DMU(2008年4月最新發(fā)布的DMU v6 Release 4.7)軟件分析數(shù)據(jù)，采用AIREML法估計(jì)中國美利奴羊遺傳參數(shù)?！緮M解決的關(guān)鍵問題】研究估計(jì)中國美利奴羊毛用性狀的遺傳參數(shù)，找到其遺傳規(guī)律，為中國美利奴羊種羊選擇、育種方案的制定提供數(shù)據(jù)支持。

1 材料與方法

1.1 材 料

1.1.1 羊群

試驗(yàn)羊群為新疆鞏乃斯種羊場2007年度出生的1 351只中國美利奴羊周歲母羊，試驗(yàn)羊只具有系譜信息、育種記錄檔案。實(shí)驗(yàn)羊群采用暖季放牧、冷季舍飼的飼養(yǎng)方式，由牧戶承包飼養(yǎng)，羊場統(tǒng)一管理。試驗(yàn)測定了周歲羊羊毛自然長度、剪毛量，并采集毛樣分析羊毛纖維直徑及羊毛纖維直徑變異系數(shù)。

1.2 方 法

1.2.1 羊毛性狀及測定

研究分析的性狀包括羊毛自然長度、剪毛量、羊毛纖維直徑、羊毛纖維直徑變異系數(shù)等4個(gè)毛用性狀。羊毛自然長度為周歲母羊體左側(cè)中線，肩胛骨后緣一掌處，順毛叢方向以直尺測量的毛叢自然狀態(tài)的長度(以cm為單位)。在周歲母羊身體左側(cè)，肩胛骨后緣一掌處，采集毛樣，利用光學(xué)纖維直徑分析儀OFA100測定羊毛纖維直徑，測定結(jié)果為平均羊毛纖維直徑(以μm為單位)；羊毛纖維直徑變異系數(shù)為測定平均羊毛纖維直徑時(shí)獲得的羊毛纖維直徑變異系數(shù)(以%表示)。剪毛量為周歲羊剪毛后的羊毛稱重(以kg為單位)。

1.2.2 數(shù)據(jù)統(tǒng)計(jì)

1.2.2.1 羊毛性狀統(tǒng)計(jì)

將周歲羊鑒定記錄、剪毛稱重記錄、羊毛纖維直徑測定記錄等數(shù)據(jù)輸入Excel數(shù)據(jù)庫，導(dǎo)入SASv8.2軟件系統(tǒng)，作描述統(tǒng)計(jì)分析，剔除不可信的極值數(shù)據(jù)。各性狀可用于分析的有效記錄數(shù)及其簡單描述性統(tǒng)計(jì)分析結(jié)果。將要計(jì)算的性狀信息排序、整理、合并，按照DMU要求的格式整理成文本型數(shù)據(jù)文件、系譜文件。表1

表1 中國美利奴羊周歲母羊毛用性狀基本統(tǒng)計(jì)量

Table 1 Data Instruction on Different Trait of Chinese Merino

性狀Trait記錄數(shù)Recode平均值A(chǔ)verage標(biāo)準(zhǔn)誤Standarderror羊毛自然長度Woollength(cm)135110.211.16羊毛纖維直徑Woolfibrediameter(μm)135118.261.72羊毛纖維直徑變異系數(shù)Coefficientofvariationoffibrediameter(%)135121.632.75剪毛量Greasyfleeceweight(kg)13513.590.70

1.2.2.2 正態(tài)性檢驗(yàn)

樣本數(shù)小于2 000，所以采用W檢驗(yàn)統(tǒng)計(jì)量進(jìn)行正態(tài)性檢驗(yàn)，利用SAS軟件INSIGHT功能模塊對(duì)中國美利奴羊羊毛自然長度、羊毛纖維直徑變異系數(shù)、羊毛纖維直徑等測定數(shù)據(jù)進(jìn)行正態(tài)分布檢驗(yàn)。

1.2.2.3 統(tǒng)計(jì)模型

采用多狀動(dòng)物模型對(duì)中國美利奴羊周歲母羊4個(gè)毛用性狀進(jìn)行方差組分及遺傳力估計(jì)，并分析協(xié)方差組分及遺傳相關(guān)。由于研究群體為2008年度鞏乃斯羊場周歲母羊，性別、年齡、季節(jié)等因素一致，因此在對(duì)中國美利奴羊周歲母羊主要生產(chǎn)性能遺傳參數(shù)估計(jì)時(shí)，模型中只考慮羊群效應(yīng)。根據(jù)分析，對(duì)羊毛自然長度、剪毛量、羊毛纖維直徑、羊毛纖維直徑變異系數(shù)等4個(gè)性狀采用公式1所示統(tǒng)計(jì)模型：

trait=μ+Xβ+Zαα+e.

(1)

利用DMU V6軟件AI模塊計(jì)算各性狀方差組分，再根據(jù)公式2、3、4，計(jì)算各性狀遺傳力、遺傳相關(guān)、表型相關(guān)：

(2)

(3)

(4)

2 結(jié)果與分析

2.1 測定數(shù)據(jù)的正態(tài)性檢驗(yàn)

用SAS軟件進(jìn)行正態(tài)分布檢驗(yàn)，測定數(shù)據(jù)的正態(tài)分布檢驗(yàn)結(jié)果，可以看出羊毛纖維直徑、羊毛纖維直徑變異系數(shù)、剪毛量、羊毛自然長度等性狀P值大于0.05，均服從正態(tài)分布。羊毛纖維直徑、羊毛纖維直徑測定服從正態(tài)分布。圖1，圖2，表2

表2 正態(tài)分布檢驗(yàn)結(jié)果

Table 2 Normal test results

變量VariableW檢驗(yàn)Shapiro-WilkP值羊毛纖維直徑(μm)Woolfibrediameter0.993194>0.05羊毛纖維直徑變異系數(shù)(%)Coefficientofvariationoffibrediameter0.986881>0.05剪毛量(kg)Greasyfleeceweight0.990823>0.05羊毛自然長度(cm)Woollength0.969678>0.05

2.2 四個(gè)羊毛性狀的遺傳力

研究表明，中國美利奴羊周歲羊羊毛自然長度、羊毛纖維直徑、羊毛纖維直徑變異系數(shù)、剪毛量等4個(gè)性狀的加性方差組分及標(biāo)準(zhǔn)誤估計(jì)值分別為0.36±0.13、1.5±0.65、2.17±1.14、0.16±0.09；殘差方差及及標(biāo)準(zhǔn)誤估計(jì)值分別為0.56±0.12、1.38±0.58、5.16±1.07、0.27±0.09；表型方差及及標(biāo)準(zhǔn)誤估計(jì)值分別為0.92±0.96、2.88±1.69、7.34±2.71、0.44±0.66。羊毛自然長度的遺傳力為0.39±0.02，羊毛纖維直徑的遺傳力為0.52±0.04，羊毛纖維直徑變異系數(shù)的遺傳力為0.30±0.02，剪毛量的遺傳力為0.37±0.04，可見4個(gè)毛用性狀的遺傳力均大于0.3。表3～6

圖1 羊毛纖維直徑數(shù)據(jù)分布

Fig.1 Data distribution of wool fiber diameter

圖2 羊毛纖維直徑變異系數(shù)數(shù)據(jù)分布

Fig.2 Data distribution of CVFD

表3 羊毛性狀的加性方差組分估計(jì)值及標(biāo)準(zhǔn)誤

表4 羊毛性狀的殘差方差組分估計(jì)值分及標(biāo)準(zhǔn)誤

表5 羊毛性狀的表型方差組分估計(jì)值分及標(biāo)準(zhǔn)誤

Table 5 Estimates of Phenotypic variance Components(δ2p)& Standard error (SE)of wool Traits

表6 羊毛性狀的遺傳力估計(jì)值分及標(biāo)準(zhǔn)誤

Table 6 Estimates of heritability & Standard error (SE)of wool Traits

項(xiàng)目Items遺傳力(?2)標(biāo)準(zhǔn)誤(SE)羊毛自然長度Woollength0.390.02羊毛纖維直徑Woolfibrediameter0.520.04羊毛纖維直徑變異系數(shù)Coefficientofvariationoffibrediameter0.300.02剪毛量Greasyfleeceweight0.370.04

2.3 四個(gè)羊毛性狀間的遺傳相關(guān)及標(biāo)準(zhǔn)誤

中國美利奴羊周歲羊4個(gè)羊毛性狀間的遺傳相關(guān)及標(biāo)準(zhǔn)誤表明，4個(gè)性狀間的遺傳相關(guān)值范圍為-0.73～0.33。其中羊毛纖維直徑性狀與羊毛自然長度之間遺傳相關(guān)為0.33，為中等程度的正相關(guān)；羊毛纖維直徑與羊毛纖維直徑變異系數(shù)之間為-0.73，二者之間存在較高的負(fù)遺傳相關(guān)。表7

表7 四個(gè)羊毛性狀間遺傳相關(guān)及標(biāo)準(zhǔn)誤

Table 7 Estimates of genetic correlations and standard errors among four wool traits of Chinese Merino

項(xiàng)目Items羊毛纖維直徑Woolfibrediameter羊毛纖維直徑變異系數(shù)Coefficientofvariationoffibrediamete剪毛量Greasyfleeceweight羊毛自然長度(Woollength)0.33(0.3)0.07(0.3)0.24(0.33)羊毛纖維直徑(Woolfibrediameter)-0.73(0.31)0.04(0.80)羊毛纖維直徑變異系數(shù)(Coefficientofvariationoffibrediameter)0.4(0.37)

2.4 四個(gè)羊毛性狀間的表型相關(guān)及顯著性

中國美利奴羊周歲羊4個(gè)羊毛性狀間的表型相關(guān)及顯著性可以看出，4個(gè)羊毛性狀間表型相關(guān)值較小，羊毛直徑變異系數(shù)及羊毛纖維直徑性狀間的表型相關(guān)不顯著。表8

表8 四個(gè)羊毛性狀間的表型相關(guān)及顯著性

Table 8 Estimates of phenotypic correlations with significance among four wool traits of Chinese Merino

注：*表示P<0.05，**表示P<0.01，ns表示不顯著Note: * indicates P<0.05, ** indicates P<0.01, ns indicates not significant

3 討 論

3.1 羊毛性狀的遺傳力

研究估計(jì)出中國美利奴羊羊毛自然長度遺傳力為0.39，高于與張亞軍等[24]估計(jì)的中國美利奴羊(新疆型)羊羊毛自然長度遺傳力0.13，低于黃錫霞等[23]估計(jì)的優(yōu)質(zhì)細(xì)毛羊羊毛自然長度遺傳力0.56，低于左北瑤等[28]估計(jì)出的優(yōu)質(zhì)細(xì)毛羊(無角型)羊毛自然長度遺傳力0.45，造成這種差異的原因可能與估計(jì)的群體遺傳結(jié)構(gòu)不同，計(jì)算方法的差異有關(guān)。Safari等[29]2005年對(duì)已報(bào)導(dǎo)的毛用羊的羊毛自然長度的遺傳力估計(jì)值進(jìn)行統(tǒng)計(jì)分析，得出羊毛自然長度遺傳力平均估計(jì)值為0.46±0.4，可見研究與國外美利奴羊羊毛自然長度遺傳力估計(jì)值平均數(shù)接近。這可能與中國美利奴羊育種過程中持續(xù)引入澳大利亞美利奴羊改良提高有關(guān)。

研究估計(jì)出中國美利奴羊羊毛纖維直徑遺傳力為0.52，高于李莉等[25]2006年用MTDFREML法估計(jì)的敖漢細(xì)毛羊羊毛細(xì)度(目測支數(shù))遺傳力值0.12、高于張亞君等[24]2008年用MTDFREML估計(jì)的新疆鞏乃斯種羊場中國美利奴羊羊毛纖維直徑遺傳力值0.36，高于黃錫霞等[23]2005年用MTDFREML估計(jì)的新疆塔城地區(qū)種羊場優(yōu)質(zhì)細(xì)毛羊羊毛細(xì)度(目測支數(shù))遺傳力值0.33，與左北瑤[28]2006年估計(jì)的新疆巴州種畜場優(yōu)質(zhì)細(xì)毛羊(無角型)羊毛纖維直徑遺傳力0.62接近。接近于Safari等[29]統(tǒng)計(jì)分析的33個(gè)已報(bào)導(dǎo)的毛用羊的羊毛纖維直徑遺傳力平均值0.59±0.02。從研究及國內(nèi)外的報(bào)導(dǎo)可以看出羊毛纖維直徑性狀屬于高遺傳力性狀，因此在細(xì)毛羊育種中，針對(duì)羊毛纖維直徑性狀可以進(jìn)行個(gè)體選擇。另外通過比較發(fā)現(xiàn)，對(duì)羊毛纖維直徑性狀的度量方法不同影響性狀遺傳力的估計(jì)值，采用羊毛細(xì)度客觀檢測數(shù)據(jù)估計(jì)遺傳力的值略高于采用目測支數(shù)數(shù)據(jù)估計(jì)值。此外，不同研究者采用的統(tǒng)計(jì)模型、計(jì)算方法以及數(shù)據(jù)結(jié)構(gòu)的差異，都是造成遺傳力估計(jì)值存在差異的重要原因。

研究估計(jì)出中國美利奴羊羊毛纖維直徑變異系數(shù)遺傳力為0.30，低于Safari等[29]統(tǒng)計(jì)分析的14個(gè)已報(bào)導(dǎo)的毛用羊的羊毛纖維直徑變異系數(shù)平均遺傳力0.52±0.04。國內(nèi)關(guān)于羊毛纖維直徑變異系數(shù)遺傳力的報(bào)導(dǎo)很少，但羊毛纖維直徑變異系數(shù)是衡量羊毛細(xì)度勻度的重要指標(biāo)，對(duì)羊毛纖維直徑變異系數(shù)遺傳力的估計(jì)將對(duì)細(xì)毛羊育種中羊毛細(xì)度勻度的選擇具有指導(dǎo)作用。

在國內(nèi)學(xué)者對(duì)中國美利奴羊、東北細(xì)毛羊、敖漢細(xì)毛羊、青海細(xì)毛羊、新吉細(xì)毛羊等細(xì)毛羊品種的剪毛量遺傳力都進(jìn)行了估計(jì)，估計(jì)值范圍為0.12～0.68[18, 22, 23, 30]。研究估計(jì)出中國美利奴羊剪毛量的遺傳力為0.37，在國內(nèi)報(bào)導(dǎo)的細(xì)毛羊剪毛量遺傳力估計(jì)范圍之內(nèi)[22-24]，與Safari[29]統(tǒng)計(jì)分析20個(gè)已報(bào)導(dǎo)的毛用羊剪毛量平均遺傳力0.37±0.02一致。

3.2 羊毛性狀間的相關(guān)分析

研究結(jié)果表明，中國美利奴羊4個(gè)羊毛性狀間的表型相關(guān)較低，但羊毛纖維直徑與羊毛自然長度、羊毛纖維直徑與羊毛纖維直徑變異系數(shù)間有較強(qiáng)的遺傳相關(guān)。羊毛纖維直徑性狀與羊毛自然長度之間遺傳相關(guān)為0.33，可在實(shí)際育種中加以運(yùn)用，針對(duì)羊毛纖維直徑選擇時(shí)，羊毛纖維直徑越細(xì)則羊毛自然長度變短，針對(duì)超細(xì)毛羊育種時(shí)要注意控制由于對(duì)細(xì)度的過度選擇造成羊毛自然長度變短。羊毛纖維直徑與羊毛纖維直徑變異系數(shù)之間的遺傳相關(guān)為-0.73，表明羊毛纖維直徑越細(xì)，則羊毛纖維直徑變異系數(shù)越大，在針對(duì)細(xì)度性狀的選擇中要注意避免因?yàn)檠蛎w維直徑變細(xì)造成細(xì)度勻度方面的損失。研究估計(jì)出的羊毛自然長度與剪毛量的遺傳相關(guān)為0.24，表明對(duì)細(xì)毛羊羊毛自然長度的正向選擇可以增加剪毛量。研究估計(jì)出的羊毛纖維直徑與剪毛量之間遺傳相關(guān)與表型相關(guān)分別為0.04、0.01，低于Safari等[29]的統(tǒng)計(jì)分析均值0.36與0.31，造成這種明顯差異的原因，這可能與數(shù)據(jù)結(jié)構(gòu)、分析方法等因素有關(guān)。

4 結(jié) 論

4.1 研究估計(jì)出中國美利奴羊羊毛自然長度、羊毛纖維直徑、羊毛纖維直徑變異系數(shù)、剪毛量的遺傳力分別為0.39、0.52、0.30、0.37。研究估計(jì)出羊毛纖維直徑及羊毛纖維直徑變異系數(shù)的遺傳相關(guān)為-0.73，羊毛自然長度與剪毛量、羊毛纖維直徑與羊毛自然長度之間遺傳相關(guān)為0.24、0.33，均為中等以上的遺傳相關(guān)。可為今后中國美利奴羊育種方案的制定、種羊遺傳評(píng)估提供數(shù)據(jù)支持。

4.2 中國美利奴羊羊毛纖維直徑性狀遺傳力較高，因此通過個(gè)體選擇，對(duì)以后世代或現(xiàn)有羊群進(jìn)行選擇，都能使羊毛纖維直徑變細(xì)，但在針對(duì)羊毛纖維直徑性狀進(jìn)行選擇時(shí)，所選的羊必須年齡相同，而且是在同一環(huán)境中飼養(yǎng)。在育種及生產(chǎn)過程中可根據(jù)個(gè)體生產(chǎn)性能進(jìn)行選配。

4.3 中國美利奴羊羊毛纖維直徑與直徑變異系數(shù)間為較高的負(fù)遺傳相關(guān)，因此在超細(xì)型細(xì)毛羊育種中要注意采取方法，降低因羊毛直徑變細(xì)造成細(xì)度勻度的損失。羊毛自然長度與剪毛量有中等程度的正遺傳相關(guān)，因此可同時(shí)對(duì)剪毛量、羊毛自然長度進(jìn)行正向選擇。

4.4 中國美利奴羊是在引進(jìn)澳大利亞美利奴羊的基礎(chǔ)上育成的，多年來中國美利奴羊不斷向各細(xì)毛羊主產(chǎn)區(qū)推廣，對(duì)細(xì)毛羊主產(chǎn)區(qū)羊毛品質(zhì)的提高起了積極作用[31]。國內(nèi)主產(chǎn)區(qū)細(xì)毛羊群體之間存在一定的遺傳聯(lián)系。中國美利奴羊周歲母羊羊毛纖維直徑遺傳力估計(jì)值對(duì)國內(nèi)各場細(xì)毛羊育種具有借鑒意義。

References)

[1] 田可川.新疆的細(xì)毛羊育種回顧與展望[C]//. 第五屆中國羊業(yè)發(fā)展大會(huì); 2008; 中國遼寧蓋州.

TIAN Ke-chuan. (2008).ReviewandProspectofFineWoolSheepBreedinginXinjiang[C]//. 2008ChinaSheepIndustryProgressProceeding.(in Chinese)

[2] 張沅編著. 現(xiàn)代動(dòng)物育種原理和方法[M]. 北京:北京農(nóng)業(yè)大學(xué)出版社,1989:81.

ZHANG Yuan. (1989).PrinciplesandMethodsofModernAnimalBreeding[M].Beijing. China Agriculture Press: 81. (in Chinese)

[3] Cloete, S., Wyk, J. V., & Neser, F. (2004). Estimates of genetic and environmental (co)variances for live weight and fleece traits in yearling south african mutton merino sheep.SouthAfricanJournalofAnimalScience, 34(1):37-43.

[4] Wyk, J. B. V., Swanepoel, J. W., & Delport, G. J. (2006). Estimates of genetic parameters for yearling body weight and fleece traits in the South African Dohne Merino population.Proceedingsofthe8thWorldCongressonGeneticsAppliedtoLivestockProduction,BeloHorizonte,MinasGerais,Brazil: 13-18 August, 2006.

[5] Di, J., Zhang, Y., Tian, K. C., Lazate, Liu, J. F., & Xu, X. M., et al. (2011). Estimation of (co)variance components and genetic parameters for growth and wool traits of chinese superfine merino sheep with the use of a multi-trait animal model.LivestockScience, 138(1-3):278-288.

[6] Ciappesoni, C. G., Ravagnolo, O., Gimeno, D., Montossi, F., & Barbieri, I. D. (2006). Estimation of genetic parameters and genetic trends for wool production and quality for the Uruguayan Merino.WorldCongressonGeneticsAppliedToLivestockProduction.

[7] Dixit, S. P., Singh, G., & Dhillon, J. S. (2011). Genetic and environmental factors affecting fleece traits in bharat merino sheep.IndianJournalofAnimalSciences, 81(1):80-83.

[8] Swp, C., Gilmour, A. R., Olivier, J. J., & van Wyk, J. B. (2004). Genetic and phenotypic trends and parameters in reproduction, greasy fleece weight and liveweight in merino lines divergently selected for multiple rearing ability.AustralianJournalofExperimentalAgriculture, 44(8):745-754.

[9] Preston, J. W. V., Hatcher, S., & Villalobos, N. L. (2013). Genetic estimates for along and across fibre diameter variation and its use to improve staple strength in Merino sheep.ProceedingsoftheTwentiethConferenceoftheAssociationfortheAdvancementofAnimalBreedingandGenetics,TranslatingScienceintoAction,Napier,NewZealand, 20th-23rd October 2013.

[10] Valera, M., Arrebola, F., Juárez, M., & Molina, A. (2009). Genetic improvement of wool production in spanish merino sheep: genetic parameters and simulation of selection strategies.AnimalProductionScience, 49(1):43-47.

[11] Bromley, C. M., Snowder, G. D., & Vleck, L. D. V. (2000). Genetic parameters among weight, prolificacy, and wool traits of columbia, polypay, rambouillet, and targhee sheep.JournalofAnimalScience,78(78):846-858.

[12] Huisman, A. E., & Brown, D. J. (2008). Genetic parameters for bodyweight, wool, and disease resistance and reproduction traits in merino sheep. 2. genetic relationships between bodyweight traits and other traits.AustralianJournalofExperimentalAgriculture, 48(9):1,186-1,193.

[13] Rose, M., & Pepper, P. M. (2000). Genetic parameters for staple length and staple strength of merino wool produced in central and north west queensland.Aaap&Asap, (13):87-90.

[14] Greeff, J.C., Schlink, A.C. (2013).Geneticparametersforstaplestrengthandcoefficientofvariationoffibrediameterinmerinowoolofdifferentstaplelength[C]//.Proceedings of the Twentieth Conference of the Association for the Advancement of Animal Breeding and Genetics, Translating Science into Action, Napier, New Zealand : 393-396.

[15] Sherlock, R. G., Lopezvillalobos, N., & Garrick, D. J. (2003). Genetic parameters for wool traits in ultra-fine New Zealand Merinos. 50yearsofDNA:ProceedingsoftheFifteenthConference,AssociationfortheAdvancementofAnimalBreedingandGenetics,Melbourne,Australia, 7-11 July 2003. (Vol.15) :277-280.

[16] Swan, A. A., Purvis, I. W., & Piper, L. R. (2008). Genetic parameters for yearling wool production, wool quality and bodyweight traits in fine wool merino sheep.AnimalProductionScience, 48(9):1,168-1,176.

[17] Safari, A., & Fogarty, N. M. (2003). Genetic parameters for sheep production traits.

[18] 商樹岐, 董維官, 焦驊, 等. 東北細(xì)毛羊遺傳力估測的討論[J]. 中國畜牧雜志,1980,(5):14-17.

SHANG Shu-qi, DONG Wei-guan, JIAO Hua, et al. (1980). The discussion of estimating Northeast fine-wool sheep heritability [J].ChineseJournalofAnimalScience, (5):14-17. (in Chinese)

[19] 齊侃虎. 良種細(xì)毛羊主要經(jīng)濟(jì)性狀的遺傳參數(shù)及優(yōu)化選種方法的研究[J]. 北京農(nóng)業(yè)大學(xué)學(xué)報(bào),1986,(2):219-229.

QI Kan-hu. (1986). The study of estimate main economic traits genetic parameters of fine-wool sheep and optimization method of breeding selection [J].JournalofBeijingAgriculturalUniversity, (2):219-229. (in Chinese)

[20] 李俊年, 阿扎提. 中國美利奴(新疆型)細(xì)毛羊遺傳參數(shù)的估計(jì)[J]. 草食家畜,1999,(1):21-23.

LI Jun-nian, Azhati. (1999). Estimate genetic parameters of Chinese Merino (Xinjiang type) wool sheep [J].Grass-FeedingLivestock, (1):21-23. (in Chinese)

[21] 李文輝, 蘇文娟, 保國俊. 甘肅高山細(xì)毛羊遺傳參數(shù)估計(jì)[J]. 中國畜牧雜志,1997,(6):24-25.

LI Wen-hui, SU Wen-juan, BAO Guo-jun. (1997). Estimate genetic parameters of Gansu Alpine Merino [J].ChineseJournalofAnimalScience, (6):24-25. (in Chinese)

[22] 梅花. 敖漢細(xì)毛羊主要數(shù)量性狀的遺傳參數(shù)估計(jì)[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學(xué)碩士論文,2005.

MEI Hua. (2005).EstimategeneticparametersmajorquantitativetraitsofAohanfine-woolsheep[D]. Master Dissertation. Inner Mongolia Agricultural University, Hohhot. (in Chinese)

[23] 黃錫霞, 吳常信, 田可川, 等. 用MTDFREML方法估計(jì)優(yōu)質(zhì)細(xì)毛羊主要經(jīng)濟(jì)性狀的遺傳參數(shù)[J]. 畜牧獸醫(yī)學(xué)報(bào),2005,(10):1 087-1 090.

HUANG Xi-xia, WU Chang-xin, TIAN Ke-chuan, et al. (20905). Estimate genetic parameters of high-quality fine-wool sheep main economic traits with MTDFREML method [J].AnimalHusbandryandVeterinaryScience, (10):1,087-1,090. (in Chinese)

[24] 張亞軍, 趙文生. 用MTDFREML估計(jì)中國美利奴羊(新疆型)產(chǎn)毛性狀的遺傳參數(shù)[J]. 中國草食動(dòng)物,2008,(1):8-10.

ZHANG Ya-jun, ZHAO Wen-sheng. (2008). Estimate genetic parameters Chinese Merino (Xinjiang type) wool traits with MTDFREML [J].ChineseHerbivore, (1):8-10. (in Chinese)

[25] 李莉. 用BLUP和MTDFREML對(duì)敖漢細(xì)毛羊遺傳評(píng)定和遺傳參數(shù)估計(jì)[D]. 呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué)碩士論文,2006.

LI Li. (2006).EstimategeneticparametersofAohanfine-woolsheepwithBLUPandMTDFREML[D]. Master Dissertation. Inner Mongolia Agricultural University; Hohhot. (in Chinese)

[26] 徐新明, 李彥飛, 付雪峰,等. 中國美利奴羊(新疆型)初生重的遺傳力估計(jì)及非遺傳因素分析[J]. 中國畜牧獸醫(yī),2014,(6):168-171.

XU Xin-ming, LI Yan-fei, FU Xue-feng,et al. (2014). Heritability Estimates and Non-genetic Factors Analysis of Birth Weight of Chinese Merino (Xinjiang type) [J].ChinaAnimalHusbandry&VeterinaryMedicine, (6):168-171. (in Chinese)

[27] Gilmour, A., & Johnson, D. (1995). Recent advances in genetic parameter estimation.

[28] 左北瑤. 優(yōu)質(zhì)細(xì)毛羊(無角類型)選育 [D]. 北京：中國農(nóng)業(yè)科學(xué)院碩士論文,2006.

ZUO Bei-yao. (2006).SelectionoftheGoodQualityFineSheepofPollType[D].Master Dissertation. Chinese Academy of Agricultural Sciences, Beijing. (in Chinese)

[29] Safari, E., Fogarty, N. M., & Gilmour, A. R. (2005). A review of genetic parameter estimates for wool, growth, meat and reproduction traits in sheep.LivestockProductionScience, 92(3):271-289.

[30] 柳楠. 良種細(xì)毛羊經(jīng)濟(jì)性狀遺傳參數(shù)和選種方法的研究[J]. 吉林農(nóng)業(yè)科學(xué),1985,(4):54-61.

LIU Nan. (1985). The Study of Economic Trait Genetic Parameters and Selection Methods on Fine-wool Sheep [J].JilinAgriculturalSciences, (4):54-61. (in Chinese)

[31] 蔣英. 中國美利奴羊由國家命名驗(yàn)收[J]. 草與畜雜志,1986,(2):41.

JIANG Ying. (1986).Chinese Merino Sheep Named by the national inspection [J].ChinaHerbivoreScience, (2):41. (in Chinese)

Fund project:Supported by National technology system construction for wool sheep industry (CARS-40); The special funds for basic science and technology research of non-profit research institutions of Xinjiang Uygur Autonomous Region (KY2015 011) and Special funds for innovation and development of scientific research institutions of Xinjiang Uygur Autonomous Region (2012027)

Using AIREML to Estimate Genetic Parameter of Chinese Merino Wool Traits

ZHANG Yan-hua1,2, YU Li-juan1, JIANG Xiao-mei1, Marziya Yasen1, Gulinur MAHATI1,Aminiguli Abulaizi1, WU Wei-wei1, Mayiken Sali1, Lazhati Aniwaier1, TIAN Kechuan1

(1.XinjiangAcademyofAnimalSciences，Urumqi830011，China；2.ChinaAgriculturalUniversity，Beijing100193，China)

【Objective】 To estimate the heritability and correlation of Chinese Merino wool traits on Kunas sheep stud farm in the hope of providing genetic parameters to develop sheep breeding and selection programs.【Method】Collecting wool length, wool fiber diameter, coefficient of variation of wool fiber diameter and greasy fleece weight of 1,351 Chinese Merino hogget ewes on Kunas sheep stud farm. Using DMU v6 software and SAS v8.2 to analyze data, employing multi-trait animal model and Average Information Restricted Maximum Likelihood (AIREML) to estimate traits variance component and genetic parameters of each trait.【Result】The result showed that heritability estimates for the Chinese Merino hogget wool length, wool fiber diameter, coefficient of variation of wool fiber diameter and greasy fleece weight were 0.39,0.52 and 0.30,0.37, respectively. Genetic correlation between wool fiber diameter and coefficient of variation of fiber diameter was -0.73, wool length and greasy fleece weight, wool length and wool fiber diameter were 0.24 and 0.33, respectively.【Conclusion】The result can provide the heritability of Chinese Merino hogget wool length, wool fibere diameter, coefficient of variation of wool fibre diameter and greasy fleece weight which can be used to develop Chinese Merino sheep breeding programs and evaluate Chinese Merino sire.

wool fiber；coefficient of variation of fiber diameter；greasy fleece weight；wool length；genetic parameters；Chinese Merino

2016-08-12

國家絨毛用羊產(chǎn)業(yè)技術(shù)體系建設(shè)(CARS-40)；自治區(qū)公益性科研院所基本科研業(yè)務(wù)費(fèi)(KY2015011)；自治區(qū)科研機(jī)構(gòu)創(chuàng)新發(fā)展專項(xiàng)資金(2012027)

張艷花(1976-)，女，甘肅人，副研究員，博士，研究方向?yàn)檠虻倪z傳育種，(E-mail)zhangyanhuav@163.com

田可川(1963-)，男，寧夏人，研究員，博士，研究方向?yàn)榻q毛用羊遺傳育種與繁殖，(E-mail)tiankechuan@163.com

10.6048/j.issn.1001-4330.2016.12.024

S826；S188

：A

：1001-4330(2016)12-2344-09