遼寧地區西門塔爾種公牛體重變化與生產性能之間的關聯分析

錢漢超，黃承俊，王大鵬，李紅宇

（1.遼寧省牧經種牛繁育中心有限公司，遼寧 沈陽 110032；2.遼寧省畜牧業經濟管理站，遼寧 沈陽 110032；3.遼寧省鐵嶺市種畜禽監督管理站，遼寧 鐵嶺 112000）

公牛的生產性能受多方因素的影響，到目前為止，一種可靠而又準確的關于體重變化在生產性能上的研究還未見報道。在牛生產過程中，不僅體重時刻變化著，而且外界環境氣候也同時在變。加上種公牛資源相對稀少，嚴格地限定其他條件而去研究體重是很難做到的。因而一些相關研究只是從不同側面進行分析[1]，并不能準確地指導生產工作。

公牛的生產不僅具有周期性，而且其生產性能在生產中也具有一定規律性，這些規律性變化都是一個品種在一定環境條件下長期互作的結果。在夏季高溫地區，牛生產性能會有所降低，在年生產曲線中表現為低點[2-3]；隨著年齡的增加，相同品種生產性能基本上按著同一變化發展[4]。這些都是研究分析體重與生產性能關系的前提條件。將體重的變化放在一定規律變化中研究能更準確地找到其與生產性能之間的聯系。

為分析西門塔爾種公牛體重變化與生產性能之間的關聯，選取相同飼養管理條件下，開始采精前2年內一個生產周期的公牛生產數據，形成年群體生產性能變化曲線。將不同體重種公牛的生產數據與之相比較，分析體重變化與產精量、精子密度之間的關系，以期得出可靠的結論用以指導相關的生產與育種工作。

1 材料與方法

1.1 數據來源遼寧省牧經種牛繁育中心有限公司的相關生產記錄。體重值為開始采精后每季度最后1個月中旬前后的測量值（空腹），時間為2016年6月至2017年6月。

1.2 分析設計

1.2.1 體重與生產性能之間關聯的建立由于無法剔除時間因素對體重與生產性能的共同干擾，體重與生產性能之間的關系只能進行間接性分析。在成年種公牛群體中體重隨時間表現為上下波動式變化，生產性能隨時間表現為規律性變化。即可建立一個相對穩定且波動不劇烈的群體生產性能變化曲線，用個體體重變化下的生產曲線與群體的作對比，若差異顯著則表明體重的變化干擾到了生產曲線，即體重與生產性能之間存在關聯。

1.2.2 群體生產性能曲線模型的建立P=G＋E＋IGE（P：數量性狀表型值；G：基因型值；E：環境效應值；IGE：基因型與環境的互作偏差效應值）；E又分為持久性環境效應（EP，環境因素對表型值的影響是單向的持久的，在群體變化規律中表現同正或同負）與暫時性環境效應（ET，環境因素對表型值的影響是雙向的，在群體變化規律中可以相互抵消），即P=G＋EP＋ET＋IGE。

1.2.3 條件假設 假設基因型效應與環境效應之間無互作，或互作很小，即IGE=0；假設固定場所下持久性環境效應對表型值的影響在同一月都是單向的；假設固定飼養管理模式下EP均相同，ET只表現在體重（T）的上下波動上。由于單一品種條件下基因型效應是單向的，將G、EP、IGE合并成U（單向影響因素），則模型可看成P=U＋T。

1.2.4 所有分析個體中U一致性的處理由于U在不同年齡或有無采精經歷上所產生的效應大小及方向是不同的，而且無法用擴大群體規模的方式來消除，因而需選取U方向一致且變化差異不大的同一生產時間段，來盡可能地消除U對群體曲線波動劇烈程度的干擾。有文獻[5]報導，西門塔爾種公牛群體在開始采精前2年次采精量與次精子密度均無差異性變化，而且均為初采牛只，與體重測量值在同一年齡段，因而選取開始采精前2年內相關公牛的生產數據進行群體生產曲線的建立準確性比較高。

1.2.5 群體生產性能曲線的建立選取所有飼養管理條件相同的同齡牛，16月齡初次采精，前2年采精期間任一個生產周期的相關生產數據，進行變化曲線的建立。

1.2.6 體重變化與生產性能關聯的判斷分析牛只的生產性能值采用體重測量日60 d（2個月）后的月次平均值，若月次平均值的變化規律符合綜合曲線的變化規律，則該牛只體重變化對生產性能無影響。曲線的高低則由U決定，與體重無關。

1.3 統計分析用Spss 17.0統計軟件進行數據的處理，與Excel 10.1辦公軟件進行相關圖的繪制。數據采用單因素方差分析，齊性分析選用LSD方法，非齊性分析選用Games-Howell方法，顯著性采用Dunnett雙尾檢驗，由Duncan方法進行描述，數值以平均值±標準差（SD）的形式表示，P＜0.05為差異有統計學意義。測量值與綜合曲線值（平均值）之間的差異采用t檢驗。適合性檢驗為卡方檢驗，符合比例為相應的群體生產性能曲線值的比。

2 結果與分析

2.1 群體生產性能曲線的建立結果圖1顯示，西門塔爾種公牛在開始采精前2年內，在一個生產周期內，其月平均次采精量呈拋物線型，在4月左右有個最大值；月平均次精子密度呈V型，在8月左右有個最低值。

2.2 相鄰季度體重變化下生產性能的變化是否服從綜合曲線變化的卡方檢驗經卡方檢驗均不顯著，即測量體重期間相鄰季度的生產性能變化是服從群體生產性能曲線變化的。說明在采精前2年內，分析牛只的體重變化在季度上對生產性能不產生影響。

圖1 西門塔爾群體開始采精前兩年內一個生產周期的綜合曲線。Fig.1 The comprehensive trend curve of the production cycle in the first two years after the Simmental bull began to collect the sperm

2.3 分析牛只生產性能的整體變化是否服從群體生產性能曲線變化的卡方檢驗結果經卡方檢驗均不顯著，即測量體重期間生產性能的整體變化是服從群體生產性能曲線變化的。說明在采精前2年內，分析牛只的體重變化對生產性能不產生影響。

2.4 測量值與群體生產性能曲線平均值之間的顯著性檢驗及變化趨勢結果由圖2可以看出月平均次采精量E牛只的最好，F牛只的最差。B、C、D、A曲線依次下降，G、H、I曲線依次上升，但下降與上升相互交錯。結合表3體重數據可以看出，在該體重分析范圍內，牛只初采體重與月平均次采精量無關，而且體重增加與否對次季度月平均次采精量不產生影響。

由圖2可以看出月平均次精子密度A、G、H牛只的較好，C牛只的最差，B牛和E牛的變化趨勢很相似。結合表3體重數據可以看出，在該體重分析范圍內，牛只初采體重與月平均次精子密度無關，而且體重增加與否對次季度月平均次精子密度也不產生影響。

由圖2可以分析出，在相同飼養管理條件下，

同齡牛在同一個生產周期內，其生產性能曲線基本保持同一趨勢。即在該前提條件下，可以通過任一月平均生產性能值的大小來評定不同牛只的好壞，而且該評定與體重無關。在此分析結果中，月平均次采精量與月平均次精子密度之間似乎存在一定的關聯，而且該關聯與體重也無關，但需大量的研究分析進一步論證。

表 2 整體分析值的卡方檢驗Table 2 The chi-square of analysis values overall

表 3 各分析值與綜合值之間的t檢驗及對應體重測量值Table 3 The T test of each analysis value and the comprehensive value and the corresponding weight measurement value

圖 2 分析牛只生產性能的走勢情況與群體生產性能變化曲線Fig.2 The trend curve of production performance of the bull analyzed and the comprehensive trend curve

3 討論

環境因素在各月對同一品種牛的影響基本是相同的[6]，而特定地區的環境因素具有周期性，因而在其他條件相似的前提下，不同年但同月的環境影響方向可以看作是相同的。在這個前提下，對于飼養管理相同的同一品種的初采公牛，在不同年份的同一年齡月份下，其生產數據具有很高的可參考性。又由于初次采精對西門塔爾種公牛前2年的生產性能無顯著影響[7]，因而對西門塔爾種公牛來說，開采前2年的生產性能具有很強的規律性，可以建立群體生產性能曲線進行指導生產與育種工作。

分析顯示，西門塔爾種公牛在開始采精前2年內的群體生產性能曲線具有很高的規律性，其月平均次采精量呈拋物線型，在4月左右有個最大值；月平均次精子密度呈V型，在8月左右有個最低值[8]。符合當地環境影響效果及生產規律，對檢測其他牛只生產性能的好壞與體重影響因素具有一定的可參考性。

卡方檢驗結果顯示，所有分析牛只的生產性能整體變化服從群體生產性能曲線，在采精前2年內分析牛只的體重變化對生產性能不產生影響。而且相鄰季度的生產性能變化也服從群體生產性能曲線，在采精前2年內，分析牛只的體重變化在季度上對生產性能不產生影響。通過進一步分析顯示，在分析的體重范圍內，牛只初采體重與月平均次采精量和月平均次精子密度無關，而且體重增加與否對次季度月平均次采精量和月平均次精子密度不產生影響。該分析結果與一些報道不同[1]，可能由于分析方法的不同而產生。

綜上分析可發現，自然飼養管理下，初采精量及精子密度與體重的大小關系不大，而且體重在100 kg之間變化不會影響相鄰季度的相關生產性能，在200 kg之間變化對整體生產性能不會產生影響。一頭牛生產性能的好壞主要受遺傳因素的影響，與體重無關。在確定牛只基因的好壞及體態優劣的前提下，可以一定范圍內減少體重來節約飼養成本及預防相關疾病的發生[9-10]。