王浩，徐新明，石剛，張小偉，魏趁，關(guān)鳴軒，何軍敏，楊雪梅，杜建文，曾偉丹，黃錫霞*，田可川*

(1. 新疆農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院，新疆 烏魯木齊 830052；2. 新疆畜牧科學(xué)院畜牧研究所，新疆 烏魯木齊 830011)

我國(guó)綿羊存欄數(shù)量較大，根據(jù)國(guó)家統(tǒng)計(jì)局?jǐn)?shù)據(jù)顯示，截至2018年底我國(guó)綿羊存欄量達(dá)到1.61億只[1]。近年來，綿羊絨毛的售價(jià)出現(xiàn)了小范圍的上漲，但在2018年以來，活體羊和羊肉價(jià)格的增加范圍較絨毛更大，羊肉的性價(jià)比較羊絨毛更好[2]。對(duì)于綿羊生產(chǎn)性能相關(guān)性狀的選育一直以來都是養(yǎng)羊業(yè)關(guān)注的重點(diǎn)問題。初生重和斷奶重是研究羔羊健康狀況和生產(chǎn)性能的重要性狀，受到遺傳、環(huán)境和母體效應(yīng)等多種因素的影響[3]。Bahreini等[4]研究發(fā)現(xiàn)新生羔羊初生重與斷乳重、6月齡體重、9月齡體重、周歲重之間存在著高度遺傳相關(guān)；Afolayan等[5]的研究表明，對(duì)羔羊的初生重進(jìn)行選擇，可以有效地改善羊群的生產(chǎn)性能和繁殖性能。

綿羊不僅是人類需要的主要肉類食品和輔助畜產(chǎn)品的來源，也是人類農(nóng)業(yè)文明傳播的主要組成部分[6]。因此，本研究通過分析非遺傳因素對(duì)多胎細(xì)毛羊羔羊初生重和斷奶重的影響，探討環(huán)境因素的影響，為今后遺傳分析提供基礎(chǔ)和依據(jù)。

1 材料與方法

1.1 數(shù)據(jù)來源

以新疆科創(chuàng)畜牧繁育中心多胎細(xì)毛羊?yàn)樵囼?yàn)動(dòng)物，使用最小二乘方差分析(GLM)對(duì)多胎細(xì)毛羊的初生重和斷奶重的因素進(jìn)行分析。收集該中心2009—2019年多胎細(xì)毛羊的產(chǎn)羔記錄數(shù)據(jù)共1 567條，包括：多胎細(xì)毛羊的個(gè)體號(hào)、出生月份、出生年份、性別、初生重、斷奶重、死亡日期和死亡原因。

1.2 飼養(yǎng)管理方式

新疆科創(chuàng)畜牧繁育中心的細(xì)毛羊以舍飼和放牧相結(jié)合的方式進(jìn)行飼養(yǎng)，每年6月至9月進(jìn)行放牧，10月至來年5月對(duì)羊群進(jìn)行舍飼管理，舍飼時(shí)早晚各飼喂1次，飼料以青貯飼料為主。母羊配種一般集中在秋季10—11月，2—4月產(chǎn)羔。

1.3 性狀及影響因素

研究性狀為多胎細(xì)毛羊的初生重和斷奶重，對(duì)這2個(gè)性狀運(yùn)用統(tǒng)計(jì)軟件進(jìn)行分析：以出生年份、出生月份、性別和母羊年齡為固定效應(yīng)(非遺傳因素)。效應(yīng)水平劃分的方法：將出生年份劃分為11個(gè)水平(2009—2019年)；將出生月份按羔羊出生的日期劃分為5個(gè)水平(1—5月，分別用1～5表示)；將母羊的年齡劃分為6個(gè)水平(2～7歲，7歲以上的合并至7歲)；將性別按羔羊出生的自然性別劃分為2個(gè)水平(公羔、母羔，分別用1、2表示)[7]。

1.4 表型數(shù)據(jù)整理

將收集的數(shù)據(jù)錄入到Excel表格中進(jìn)行整理，為保證數(shù)據(jù)分析的準(zhǔn)確性，剔除無號(hào)、補(bǔ)號(hào)區(qū)分不出母羊年齡的個(gè)體數(shù)據(jù)，并運(yùn)用平均值加減3倍標(biāo)準(zhǔn)差的方式確定初生重和斷奶重的范圍，剔除極端值和異常值等不合理的記錄數(shù)據(jù)。

1.5 統(tǒng)計(jì)分析

采用SAS 9.2統(tǒng)計(jì)軟件GLM程序?qū)Ω鱾€(gè)固定效應(yīng)進(jìn)行方差分析，求得最小二乘均值，并進(jìn)行多重比較，其模型如下：

Yijklm=μ+ai+nj+mk+sl+eijklm，

式中，Yijklm表示性狀的觀察值；μ是群體均值；ai母羊年齡第i個(gè)水平效應(yīng)；nj是出生年份第j個(gè)水平效應(yīng)；mk是出生月份第k個(gè)水平效應(yīng)；sl是性別第l個(gè)水平效應(yīng)；eijklm是隨機(jī)殘差效應(yīng)。

2 結(jié)果與分析

2.1 描述性統(tǒng)計(jì)

不同年度多胎細(xì)毛羊的生產(chǎn)資料所記錄的2個(gè)性狀(初生重和斷奶重)的數(shù)據(jù)有部分缺失，是由于不同出生類型多胎細(xì)毛羊身體健康狀況有差異，由于營(yíng)養(yǎng)不良體弱等原因?qū)е虏糠侄嗵ゼ?xì)毛羊個(gè)體出現(xiàn)死胎或未斷奶便已經(jīng)死亡，或者部分多胎細(xì)毛羊個(gè)體未出生便已經(jīng)死亡，未收集到2018年和2019年羔羊的斷奶記錄數(shù)據(jù)，所以造成初生重和斷奶重的數(shù)據(jù)不均衡。因此，初生重和斷奶重性狀的樣本量并非一致。由表1可知，初生重和斷奶重的實(shí)際樣本量分別為1 404條和926條，公羔的初生重和斷奶重均值均高于母羊。

表1 初生重和斷奶重的描述性統(tǒng)計(jì)

2.2 產(chǎn)羔率分析

由圖1可以看出，從2009—2019年多胎細(xì)毛羊的產(chǎn)羔率呈現(xiàn)波動(dòng)上升趨勢(shì)。其中，2019年產(chǎn)羔率最高，達(dá)到184.51%，產(chǎn)多胎的母羊數(shù)目最多，達(dá)到90只。說明多胎細(xì)毛羊的選育有明顯的成效。

圖1 2009—2019年多胎細(xì)毛羊產(chǎn)羔率統(tǒng)計(jì)分析

2.3 不同因素對(duì)初生重、斷奶重的影響

由表2可知，出生年份和性別對(duì)各性狀均具有極顯著影響(P<0.01)，公羔的初生重和斷奶重高于母羔；出生月份和母羊年齡對(duì)初生重和斷奶重均無顯著影響(P>0.05)。

表2 不同因素對(duì)初生重和斷奶重的最小二乘方差分析

2.4 各水平的多重比較

由表3可見，出生年份和性別對(duì)初生重和斷奶重有極顯著影響(P<0.01)，出生月份和母羊年齡對(duì)初生重和斷奶重沒有顯著影響(P>0.05)。羔羊出生年份不同，初生重和斷奶重也有一定的區(qū)別，出生在2017年和2019年初生重?zé)o顯著差異(P>0.05)，但均顯著低于其他年份，出生在2014年的羔羊斷奶重最大，顯著大于其他年份(P<0.05)。對(duì)于出生月份，初生重在3月最大，出生在2月和3月的羔羊無顯著差異(P>0.05)，1月、2月、4月和5月出生的羔羊無顯著差異(P>0.05)，出生在5月的羔羊斷奶重最低，與其他月份差異顯著(P<0.05)。不同性別羔羊的初生重和斷奶重差異極顯著(P<0.01)，母羊的初生重和斷奶重低于公羊。不同年齡母羊所產(chǎn)的羔羊初生重差異不顯著(P>0.05)，當(dāng)母羊年齡為6歲時(shí)，羔羊的斷奶重最小，顯著低于其他年份(P>0.05)。

表3 初生重和斷奶重的最小二乘均值(LSM)及標(biāo)準(zhǔn)誤(SE)

續(xù)表3

3 討論

在全球眾多的綿羊品種中，除了少量的高產(chǎn)、性成熟比較早或者是常年發(fā)情的品種以外，大部分品種都是季節(jié)性發(fā)情的單產(chǎn)類型[8]。產(chǎn)羔率是細(xì)毛羊育種中的一個(gè)重要經(jīng)濟(jì)性狀，在細(xì)毛羊的育種過程中應(yīng)當(dāng)依據(jù)經(jīng)濟(jì)性狀的重要性進(jìn)行選育。

平衡育種的理念一直是科研育種工作關(guān)注的焦點(diǎn)[9]，隨著豐富的統(tǒng)計(jì)方法及其軟件的開發(fā)，使整個(gè)過程變得更易操作，在育種過程中盡可能地考慮非遺傳因素和遺傳因素的影響，對(duì)原始數(shù)據(jù)進(jìn)行系統(tǒng)全面地統(tǒng)計(jì)分析，為提高多胎細(xì)毛羊遺傳參數(shù)估計(jì)的準(zhǔn)確性提供更為可靠的依據(jù)，為生產(chǎn)實(shí)踐工作提供更加科學(xué)的指導(dǎo)，進(jìn)而提高育種的效率[10]。梅步俊等[11]的報(bào)道顯示山羊產(chǎn)羔數(shù)具有低遺傳力，高產(chǎn)羔能力可能受到飼養(yǎng)管理方式和環(huán)境等因素的多重影響。

3.1 出生年份對(duì)初生重和斷奶重的影響

出生年份對(duì)初生重和斷奶重有極顯著影響，出生在2017年和2019年的初生重?zé)o顯著差異，但均顯著低于其他年份，出生在2018年的初生重與其他年份差異顯著；出生在2014年的羔羊斷奶重最大，顯著大于其他年份。本研究所得結(jié)果與黃錫霞等[12]、楊潔等[13]對(duì)相關(guān)性狀的研究結(jié)果一致，表明集自然生態(tài)和飼養(yǎng)管理等2個(gè)條件于一體的年份是一個(gè)復(fù)雜的綜合效應(yīng)，不能人為控制。由于細(xì)毛羊選種工作每年都要進(jìn)行，因此，年份效應(yīng)對(duì)選種本身的影響是不大的，但是，可能會(huì)影響到育種值的估計(jì)和遺傳參數(shù)的準(zhǔn)確性，選擇效果會(huì)因此而受到影響。由此看出，在遺傳參數(shù)育種值估計(jì)方面，年份是一個(gè)值得關(guān)注的因素。

3.2 出生月份對(duì)初生重和斷奶重的影響

出生月份既可以反映羊羔的出生時(shí)間，也可以提供細(xì)毛羊的配種繁育時(shí)間，而這些時(shí)間也是與環(huán)境空間緊密聯(lián)系的。本研究中，出生月份對(duì)初生重和斷奶重?zé)o顯著影響，但有一定差異，初生重在3月份最大，出生在5月份的羔羊斷奶重最低。這與一些文獻(xiàn)報(bào)道的出生月份對(duì)綿羊有顯著影響不一致。李秀麗等[14-15]研究表明，出生月份對(duì)細(xì)毛羊的斷奶重有極顯著影響。分析可能是因?yàn)槠渌芯康难驁?chǎng)母羊的配種時(shí)間有差異, 但所有羔羊同期斷奶, 從而使得不同出生月份羔羊的斷奶重不同。而本研究的細(xì)毛羊群體沒有強(qiáng)制羔羊在同一天斷乳，這可能是結(jié)果存在差異的原因。張劍搏等[16]研究表明，出生月份、群別及其一些非遺傳因素間的互作效應(yīng)對(duì)高山美利奴羊的初生重和斷奶重均有極顯著影響。范景勝等[17]報(bào)道出生月份對(duì)簡(jiǎn)陽(yáng)大耳羊影響顯著。細(xì)毛羊發(fā)情主要集中在每年10到11月份，所以生產(chǎn)中的最佳配種期應(yīng)該為10到11月份。因此，多胎細(xì)毛羊在2、3和4月份產(chǎn)羔都是比較合適的[18]。本研究結(jié)果與上述研究結(jié)果不同可能是由于品種、地理環(huán)境的不同，或者是沒有考慮互作效應(yīng)造成的。

3.3 性別對(duì)初生重和斷奶重的影響

本研究中，性別對(duì)多胎細(xì)毛母羊的初生重和斷奶重有極顯著影響，母羊的初生重和斷奶重低于公羊。研究結(jié)果與陶衛(wèi)東等[19]報(bào)道的性別對(duì)南疆絨山羊羔羊初生重影響不顯著的結(jié)果不一致。而俞春山等[20]報(bào)道的性別對(duì)中衛(wèi)山羊羔羊初生重有顯著影響，這與本研究結(jié)果是一致的。這些差異可能是由于新疆獨(dú)特的地理和氣候環(huán)境或所研究品種不同所引起的。公羊和母羊的生理結(jié)構(gòu)、激素以及日增重不同，所以初生重和斷奶重會(huì)在公羔與母羔之間產(chǎn)生一定的差異。

3.4 母羊年齡對(duì)初生重和斷奶重的影響

本研究中，母羊年齡效應(yīng)對(duì)羔羊的初生重和斷奶重?zé)o顯著影響，不同年齡母羊所產(chǎn)的羔羊初生重差異不顯著，當(dāng)母羊年齡為6歲時(shí)，羔羊的斷奶重最小。研究結(jié)果與羅保華等[21]報(bào)道的母羊年齡對(duì)內(nèi)蒙古絨山羊的體重性狀影響極顯著的結(jié)果不一致?？赡苡捎诒狙芯康娜后w較小，導(dǎo)致差異不顯著。

4 結(jié)論

本研究利用最小二乘方差分析的方法對(duì)新疆科創(chuàng)畜牧繁育中心2009—2019年多胎細(xì)毛羊的初生重和斷奶重性狀進(jìn)行分析，證實(shí)出生年份和性別對(duì)多胎細(xì)毛羊的初生重和斷奶重性狀有顯著影響，而出生月份和母羊年齡對(duì)初生重和斷奶重性狀均沒有顯著影響。分析結(jié)果為改善多胎細(xì)毛羊的生產(chǎn)性能和提高成活率等研究提供理論依據(jù)，為今后培育優(yōu)質(zhì)的多胎細(xì)毛羊奠定基礎(chǔ)。