遼寧地區(qū)西門塔爾種公牛體重變化與生產(chǎn)性能之間的關(guān)聯(lián)分析

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

（1.遼寧省牧經(jīng)種牛繁育中心有限公司，遼寧 沈陽(yáng) 110032；2.遼寧省畜牧業(yè)經(jīng)濟(jì)管理站，遼寧 沈陽(yáng) 110032；3.遼寧省鐵嶺市種畜禽監(jiān)督管理站，遼寧 鐵嶺 112000）

公牛的生產(chǎn)性能受多方因素的影響，到目前為止，一種可靠而又準(zhǔn)確的關(guān)于體重變化在生產(chǎn)性能上的研究還未見(jiàn)報(bào)道。在牛生產(chǎn)過(guò)程中，不僅體重時(shí)刻變化著，而且外界環(huán)境氣候也同時(shí)在變。加上種公牛資源相對(duì)稀少，嚴(yán)格地限定其他條件而去研究體重是很難做到的。因而一些相關(guān)研究只是從不同側(cè)面進(jìn)行分析[1]，并不能準(zhǔn)確地指導(dǎo)生產(chǎn)工作。

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

為分析西門塔爾種公牛體重變化與生產(chǎn)性能之間的關(guān)聯(lián)，選取相同飼養(yǎng)管理?xiàng)l件下，開(kāi)始采精前2年內(nèi)一個(gè)生產(chǎn)周期的公牛生產(chǎn)數(shù)據(jù)，形成年群體生產(chǎn)性能變化曲線。將不同體重種公牛的生產(chǎn)數(shù)據(jù)與之相比較，分析體重變化與產(chǎn)精量、精子密度之間的關(guān)系，以期得出可靠的結(jié)論用以指導(dǎo)相關(guān)的生產(chǎn)與育種工作。

1 材料與方法

1.1 數(shù)據(jù)來(lái)源遼寧省牧經(jīng)種牛繁育中心有限公司的相關(guān)生產(chǎn)記錄。體重值為開(kāi)始采精后每季度最后1個(gè)月中旬前后的測(cè)量值（空腹），時(shí)間為2016年6月至2017年6月。

1.2 分析設(shè)計(jì)

1.2.1 體重與生產(chǎn)性能之間關(guān)聯(lián)的建立由于無(wú)法剔除時(shí)間因素對(duì)體重與生產(chǎn)性能的共同干擾，體重與生產(chǎn)性能之間的關(guān)系只能進(jìn)行間接性分析。在成年種公牛群體中體重隨時(shí)間表現(xiàn)為上下波動(dòng)式變化，生產(chǎn)性能隨時(shí)間表現(xiàn)為規(guī)律性變化。即可建立一個(gè)相對(duì)穩(wěn)定且波動(dòng)不劇烈的群體生產(chǎn)性能變化曲線，用個(gè)體體重變化下的生產(chǎn)曲線與群體的作對(duì)比，若差異顯著則表明體重的變化干擾到了生產(chǎn)曲線，即體重與生產(chǎn)性能之間存在關(guān)聯(lián)。

1.2.2 群體生產(chǎn)性能曲線模型的建立P=G＋E＋I(xiàn)GE（P：數(shù)量性狀表型值；G：基因型值；E：環(huán)境效應(yīng)值；IGE：基因型與環(huán)境的互作偏差效應(yīng)值）；E又分為持久性環(huán)境效應(yīng)（EP，環(huán)境因素對(duì)表型值的影響是單向的持久的，在群體變化規(guī)律中表現(xiàn)同正或同負(fù)）與暫時(shí)性環(huán)境效應(yīng)（ET，環(huán)境因素對(duì)表型值的影響是雙向的，在群體變化規(guī)律中可以相互抵消），即P=G＋EP＋ET＋I(xiàn)GE。

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

1.2.4 所有分析個(gè)體中U一致性的處理由于U在不同年齡或有無(wú)采精經(jīng)歷上所產(chǎn)生的效應(yīng)大小及方向是不同的，而且無(wú)法用擴(kuò)大群體規(guī)模的方式來(lái)消除，因而需選取U方向一致且變化差異不大的同一生產(chǎn)時(shí)間段，來(lái)盡可能地消除U對(duì)群體曲線波動(dòng)劇烈程度的干擾。有文獻(xiàn)[5]報(bào)導(dǎo)，西門塔爾種公牛群體在開(kāi)始采精前2年次采精量與次精子密度均無(wú)差異性變化，而且均為初采牛只，與體重測(cè)量值在同一年齡段，因而選取開(kāi)始采精前2年內(nèi)相關(guān)公牛的生產(chǎn)數(shù)據(jù)進(jìn)行群體生產(chǎn)曲線的建立準(zhǔn)確性比較高。

1.2.5 群體生產(chǎn)性能曲線的建立選取所有飼養(yǎng)管理?xiàng)l件相同的同齡牛，16月齡初次采精，前2年采精期間任一個(gè)生產(chǎn)周期的相關(guān)生產(chǎn)數(shù)據(jù)，進(jìn)行變化曲線的建立。

1.2.6 體重變化與生產(chǎn)性能關(guān)聯(lián)的判斷分析牛只的生產(chǎn)性能值采用體重測(cè)量日60 d（2個(gè)月）后的月次平均值，若月次平均值的變化規(guī)律符合綜合曲線的變化規(guī)律，則該牛只體重變化對(duì)生產(chǎn)性能無(wú)影響。曲線的高低則由U決定，與體重?zé)o關(guān)。

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

2 結(jié)果與分析

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

2.2 相鄰季度體重變化下生產(chǎn)性能的變化是否服從綜合曲線變化的卡方檢驗(yàn)經(jīng)卡方檢驗(yàn)均不顯著，即測(cè)量體重期間相鄰季度的生產(chǎn)性能變化是服從群體生產(chǎn)性能曲線變化的。說(shuō)明在采精前2年內(nèi)，分析牛只的體重變化在季度上對(duì)生產(chǎn)性能不產(chǎn)生影響。

圖1 西門塔爾群體開(kāi)始采精前兩年內(nèi)一個(gè)生產(chǎn)周期的綜合曲線。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 分析牛只生產(chǎn)性能的整體變化是否服從群體生產(chǎn)性能曲線變化的卡方檢驗(yàn)結(jié)果經(jīng)卡方檢驗(yàn)均不顯著，即測(cè)量體重期間生產(chǎn)性能的整體變化是服從群體生產(chǎn)性能曲線變化的。說(shuō)明在采精前2年內(nèi)，分析牛只的體重變化對(duì)生產(chǎn)性能不產(chǎn)生影響。

2.4 測(cè)量值與群體生產(chǎn)性能曲線平均值之間的顯著性檢驗(yàn)及變化趨勢(shì)結(jié)果由圖2可以看出月平均次采精量E牛只的最好，F(xiàn)牛只的最差。B、C、D、A曲線依次下降，G、H、I曲線依次上升，但下降與上升相互交錯(cuò)。結(jié)合表3體重?cái)?shù)據(jù)可以看出，在該體重分析范圍內(nèi)，牛只初采體重與月平均次采精量無(wú)關(guān)，而且體重增加與否對(duì)次季度月平均次采精量不產(chǎn)生影響。

由圖2可以看出月平均次精子密度A、G、H牛只的較好，C牛只的最差，B牛和E牛的變化趨勢(shì)很相似。結(jié)合表3體重?cái)?shù)據(jù)可以看出，在該體重分析范圍內(nèi)，牛只初采體重與月平均次精子密度無(wú)關(guān)，而且體重增加與否對(duì)次季度月平均次精子密度也不產(chǎn)生影響。

由圖2可以分析出，在相同飼養(yǎng)管理?xiàng)l件下，

同齡牛在同一個(gè)生產(chǎn)周期內(nèi)，其生產(chǎn)性能曲線基本保持同一趨勢(shì)。即在該前提條件下，可以通過(guò)任一月平均生產(chǎn)性能值的大小來(lái)評(píng)定不同牛只的好壞，而且該評(píng)定與體重?zé)o關(guān)。在此分析結(jié)果中，月平均次采精量與月平均次精子密度之間似乎存在一定的關(guān)聯(lián)，而且該關(guān)聯(lián)與體重也無(wú)關(guān)，但需大量的研究分析進(jìn)一步論證。

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

表 3 各分析值與綜合值之間的t檢驗(yàn)及對(duì)應(yīng)體重測(cè)量值Table 3 The T test of each analysis value and the comprehensive value and the corresponding weight measurement value

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

3 討論

環(huán)境因素在各月對(duì)同一品種牛的影響基本是相同的[6]，而特定地區(qū)的環(huán)境因素具有周期性，因而在其他條件相似的前提下，不同年但同月的環(huán)境影響方向可以看作是相同的。在這個(gè)前提下，對(duì)于飼養(yǎng)管理相同的同一品種的初采公牛，在不同年份的同一年齡月份下，其生產(chǎn)數(shù)據(jù)具有很高的可參考性。又由于初次采精對(duì)西門塔爾種公牛前2年的生產(chǎn)性能無(wú)顯著影響[7]，因而對(duì)西門塔爾種公牛來(lái)說(shuō)，開(kāi)采前2年的生產(chǎn)性能具有很強(qiáng)的規(guī)律性，可以建立群體生產(chǎn)性能曲線進(jìn)行指導(dǎo)生產(chǎn)與育種工作。

分析顯示，西門塔爾種公牛在開(kāi)始采精前2年內(nèi)的群體生產(chǎn)性能曲線具有很高的規(guī)律性，其月平均次采精量呈拋物線型，在4月左右有個(gè)最大值；月平均次精子密度呈V型，在8月左右有個(gè)最低值[8]。符合當(dāng)?shù)丨h(huán)境影響效果及生產(chǎn)規(guī)律，對(duì)檢測(cè)其他牛只生產(chǎn)性能的好壞與體重影響因素具有一定的可參考性。

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

綜上分析可發(fā)現(xiàn)，自然飼養(yǎng)管理下，初采精量及精子密度與體重的大小關(guān)系不大，而且體重在100 kg之間變化不會(huì)影響相鄰季度的相關(guān)生產(chǎn)性能，在200 kg之間變化對(duì)整體生產(chǎn)性能不會(huì)產(chǎn)生影響。一頭牛生產(chǎn)性能的好壞主要受遺傳因素的影響，與體重?zé)o關(guān)。在確定牛只基因的好壞及體態(tài)優(yōu)劣的前提下，可以一定范圍內(nèi)減少體重來(lái)節(jié)約飼養(yǎng)成本及預(yù)防相關(guān)疾病的發(fā)生[9-10]。