趙開(kāi)森，王 杰，孫文強(qiáng)，賈先波，高 揚(yáng)，曾 偉，譚 熊，李 健，賴(lài)松家*

（1. 四川農(nóng)業(yè)大學(xué)動(dòng)物科技學(xué)院，成都 611130；2. 四川省陽(yáng)平種牛場(chǎng)，四川 眉山 620360）

產(chǎn)奶量是泌乳性能最直觀的體現(xiàn)，是決定牧場(chǎng)經(jīng)濟(jì)效益的關(guān)鍵指標(biāo)[1]。大多數(shù)指定的西門(mén)塔爾牛的綜合選擇指數(shù)都涉及了體形的選擇[2-3]?；I(yè)浩等[4]研究體現(xiàn)了乳房線性指標(biāo)與西門(mén)塔爾牛產(chǎn)奶量的關(guān)聯(lián)，吐?tīng)栠d江·吾木爾艾力等[5]研究表明馬的產(chǎn)奶量與體高、胸圍、體長(zhǎng)等體尺性狀顯著相關(guān)。西門(mén)塔爾牛關(guān)于體尺性狀的研究主要集中于生長(zhǎng)發(fā)育[6-8]，而各體尺性狀與產(chǎn)奶量的相關(guān)性研究還比較缺乏。研究發(fā)現(xiàn)，體高、胸圍等傳統(tǒng)體尺性狀對(duì)荷斯坦牛的產(chǎn)奶性能均有不同的影響[9-13]，肛殖距作為一個(gè)新興的研究熱點(diǎn)，與荷斯坦奶牛的繁殖性能和產(chǎn)奶性能均有密切的相關(guān)性[9,14]，但關(guān)于肛殖距對(duì)西門(mén)塔爾牛繁殖和產(chǎn)奶性能的影響還未見(jiàn)報(bào)道。因此有必要將西門(mén)塔爾牛的各體尺性狀與產(chǎn)奶量的關(guān)系進(jìn)行系統(tǒng)分析。本研究通過(guò)對(duì)四川地區(qū)的西門(mén)塔爾牛群體的肛殖距、體高、胸圍、體斜長(zhǎng)和管?chē)润w尺指標(biāo)進(jìn)行分析，闡明各體尺性狀與產(chǎn)奶量的相關(guān)關(guān)系，同時(shí)構(gòu)建西門(mén)塔爾牛產(chǎn)奶量的回歸方程可為實(shí)際的產(chǎn)奶量測(cè)定提供理論參考。

1 材料和方法

1.1 實(shí)驗(yàn)動(dòng)物

試驗(yàn)于四川省陽(yáng)平種牛場(chǎng)進(jìn)行，選擇系譜資料完整的西門(mén)塔爾牛成年母牛378頭。

1.2 體尺指標(biāo)數(shù)據(jù)測(cè)定

按照常規(guī)方法對(duì)選擇的成年母牛的肛殖距、體高、胸圍、體斜長(zhǎng)和管?chē)润w尺指標(biāo)進(jìn)行測(cè)量并記錄。

1.3 產(chǎn)奶量數(shù)據(jù)

378頭西門(mén)塔爾牛成年母牛的產(chǎn)奶量數(shù)據(jù)均來(lái)自該牛場(chǎng)的檔案記錄，中國(guó)荷斯坦奶牛產(chǎn)奶量及奶成分?jǐn)?shù)據(jù)來(lái)源于四川雪寶鴻豐牧場(chǎng)[9]、常玲玲等[15]和余詩(shī)強(qiáng)等[16]文獻(xiàn)記錄。

1.4 數(shù)據(jù)統(tǒng)計(jì)與分析

利用Microsoft Excel 2016軟件對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行初步整理，剔除各個(gè)指標(biāo)的異常值，然后使用SPSS 27.0進(jìn)行統(tǒng)計(jì)分析，采用t檢驗(yàn)進(jìn)行組間差異分析，并以各體形性狀為自變量，產(chǎn)奶量為因變量，采用多元回歸分析的逐步回歸法建立中國(guó)西門(mén)塔爾牛產(chǎn)奶量與體形指標(biāo)的最佳回歸方程，剖析各性狀對(duì)產(chǎn)奶量的影響作用。

2 結(jié)果與分析

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

經(jīng)過(guò)對(duì)數(shù)據(jù)初步處理后，西門(mén)塔爾牛7 個(gè)體形數(shù)據(jù)和3 個(gè)產(chǎn)奶性能相關(guān)數(shù)據(jù)的描述性統(tǒng)計(jì)見(jiàn)表1。體高、胸圍、腹圍等6個(gè)體尺性狀的變異系數(shù)在4.43%～8.04%之間，肛殖距、產(chǎn)奶量的變異系數(shù)均高于10.00%，乳蛋白率和乳脂率的變異系數(shù)較低（均低于2.00%）。說(shuō)明在群體選育水平對(duì)肛殖距和產(chǎn)奶量的選育水平較低，性狀在不同的西門(mén)塔爾牛個(gè)體之間的差異較大，特別是肛殖距這一性狀，因?yàn)榉植驾^寬可能導(dǎo)致數(shù)據(jù)的離散性較大。將西門(mén)塔爾牛的產(chǎn)奶量、乳脂率、乳蛋白率與中國(guó)荷斯坦牛[15-16]進(jìn)行比較發(fā)現(xiàn)，西門(mén)塔爾牛的產(chǎn)奶量明顯低于中國(guó)荷斯坦奶牛，而乳脂率和乳蛋白率略高于中國(guó)荷斯坦奶牛（圖1）。

圖1 西門(mén)塔爾牛與中國(guó)荷斯坦奶牛產(chǎn)奶量、乳蛋白率和乳脂率的差異Figure 1 Differences in milk production, milk protein percentage, and milk fat percentage between Simmental cattle and Chinese Holstein cattle

表1 西門(mén)塔爾牛體尺性狀及產(chǎn)奶性狀的描述性統(tǒng)計(jì)Table 1 Descriptive statistics of body size traits and milk production traits in Simmental cattle

2.2 肛殖距、體尺及產(chǎn)奶量相關(guān)性分析

相關(guān)性分析表明（見(jiàn)表2），產(chǎn)奶量與胸圍之間呈現(xiàn)中等相關(guān)性(r=0.300 5,P＜0.01）；而胸圍與腹圍、肛殖距都存在顯著的相關(guān)關(guān)系(P＜0.05)，其中胸圍與腹圍相關(guān)性最高(r=0.718 4)，與肛殖距的相關(guān)關(guān)系次之；此外，我們發(fā)現(xiàn)胸圍與體高，體高與體斜長(zhǎng)、體直長(zhǎng)以及管?chē)南嚓P(guān)性達(dá)到顯著水平(P＜0.05）。由此可知，胸圍可以直接反映西門(mén)塔爾牛的產(chǎn)奶量，而體高、腹圍、體斜長(zhǎng)及管?chē)? 個(gè)體尺指標(biāo)直接存在顯著相關(guān)關(guān)系，可以間接評(píng)價(jià)西門(mén)塔爾牛的產(chǎn)奶量，可以用來(lái)間接反映西門(mén)塔爾牛產(chǎn)奶量的重要指標(biāo)，此外，在之前關(guān)于中國(guó)荷斯坦牛產(chǎn)奶量的研究中，發(fā)現(xiàn)肛殖距可作為影響產(chǎn)奶量的候選性狀，因此在本研究中也將這一指標(biāo)考慮在內(nèi)。在遺傳育種工作中可以通過(guò)考慮各性狀之間的相關(guān)性來(lái)選育高產(chǎn)的中國(guó)西門(mén)塔爾牛。

表2 西門(mén)塔爾牛各性狀之間的相關(guān)性Table 2 Correlation between various traits in Simmental cattle

2.3 西門(mén)塔爾牛產(chǎn)奶量與肛殖距、各體尺性狀的回歸分析

以西門(mén)塔爾牛的肛殖距、各體尺性狀作為自變量，產(chǎn)奶量作為因變量，建立關(guān)于產(chǎn)奶量的回歸方程：Y=4 152.38-22.78X1+16.29X2-2.82X3+35.63X4-30.66X5-1.15X6，其中X1為體高，X2為胸圍，X3為體斜長(zhǎng)，X4為體直長(zhǎng)，X5為管?chē)?，X6為肛殖距。

為了驗(yàn)證以上回歸方程的準(zhǔn)確性，選用該牛場(chǎng)191頭成年母牛的體尺數(shù)據(jù)代入回歸方程預(yù)測(cè)其產(chǎn)奶量，并利用SPSS 27.0軟件的配對(duì)樣本t檢驗(yàn)(置信區(qū)間95%)來(lái)分析實(shí)際產(chǎn)奶量和預(yù)測(cè)產(chǎn)奶量之間的差異見(jiàn)表3，結(jié)果發(fā)現(xiàn)實(shí)際產(chǎn)奶量和預(yù)測(cè)產(chǎn)奶之間沒(méi)有顯著差異(P=0.96)，表明本試驗(yàn)建立的關(guān)于中國(guó)西門(mén)塔爾牛產(chǎn)奶量回歸預(yù)測(cè)模型具有統(tǒng)計(jì)學(xué)意義，模型建立比較成功。

表3 西門(mén)塔爾牛實(shí)際產(chǎn)奶量和回歸方程預(yù)測(cè)產(chǎn)奶量樣本配對(duì)t檢驗(yàn)Table 3 Paired t-test for actual milk production and regression equation predicted milk production samples of Simmental cattle

3 討論

牛的體尺性狀不僅反映牛的生長(zhǎng)發(fā)育狀況，而且也直接體現(xiàn)了牛不同生理階段飼養(yǎng)管理水平及遺傳因素影響，進(jìn)而間接影響其產(chǎn)奶性能的發(fā)揮[17-19]。程月等[20]通過(guò)對(duì)奶牛體尺、體重和產(chǎn)奶量的測(cè)量和分析發(fā)現(xiàn)體斜長(zhǎng)、胸寬、胸深、腰高、胸圍、腹圍對(duì)高峰奶量有極顯著影響，但體高對(duì)于產(chǎn)奶量的影響較小，未達(dá)到顯著水平。然而徐啟明[21]通過(guò)對(duì)無(wú)為縣76 頭后備母牛體尺和產(chǎn)奶量的相關(guān)性分析發(fā)現(xiàn)體高、胸圍、管?chē)c產(chǎn)奶量存在正相關(guān)關(guān)系，但相關(guān)性的強(qiáng)弱、是否顯著可能與牛的品種、飼養(yǎng)地區(qū)有關(guān)。本研究中西門(mén)塔爾牛的產(chǎn)奶量與胸圍存在顯著的正相關(guān)關(guān)系，與體高、體斜長(zhǎng)、管?chē)鹊南嚓P(guān)關(guān)系與荷斯坦奶牛相比較弱且沒(méi)有達(dá)到顯著水平，與程郁昕等[22]和吳錦波等[23]研究結(jié)果相似，與吐?tīng)栠d江·吾木爾艾力等[24]報(bào)道的產(chǎn)奶量與體高、管?chē)蕵O顯著正相關(guān)存在差異。在西門(mén)塔爾牛的選育中，由于各表型性狀存在一定的相關(guān)關(guān)系，在改良產(chǎn)奶量的同時(shí)不能忽略體尺性狀的改良。

多元回歸分析利用與目標(biāo)性狀具有較強(qiáng)相關(guān)性的其他性狀建立回歸模型，具有準(zhǔn)確性高、實(shí)用性強(qiáng)等特點(diǎn)，是重要的目標(biāo)性狀預(yù)測(cè)方法[25]。經(jīng)過(guò)分析，最終確定體高（X1）、胸圍(X2)、體斜長(zhǎng)(X3)、體直長(zhǎng)(X4)、管?chē)?X5)和肛殖距(X6) 6個(gè)指標(biāo)作為自變量建立回歸模型：Y=4 152.38-22.78X1+16.29X2-2.82X3+35.63X4-30.66X5-1.15X6。相比劉舜才等[26]建立的黑白花奶牛產(chǎn)奶量回歸方程考慮了肛殖距和管?chē)鷮?duì)產(chǎn)奶量的影響，預(yù)測(cè)的結(jié)果更加接近實(shí)際產(chǎn)奶量，同時(shí)也完善了虎業(yè)浩等[4]研究缺乏體尺性狀數(shù)據(jù)的不足。結(jié)合實(shí)際生產(chǎn)，發(fā)現(xiàn)西門(mén)塔爾牛的預(yù)測(cè)產(chǎn)奶量與實(shí)際產(chǎn)奶量無(wú)顯著差異，表明本模型具有統(tǒng)計(jì)學(xué)意義，并且所選的6 個(gè)自變量，測(cè)量方便易行，可應(yīng)用于生產(chǎn)之中。

4 結(jié)論

本研究發(fā)現(xiàn)西門(mén)塔爾牛的產(chǎn)奶量與胸圍存在顯著的相關(guān)關(guān)系，構(gòu)建的基于肛殖距和體尺性狀的產(chǎn)奶量預(yù)測(cè)模型：Y=4 152.38-22.78X1+16.29X2-2.82X3+35.63X4-30.66X5-1.15X6，可應(yīng)用于估測(cè)西門(mén)塔爾牛產(chǎn)奶量的實(shí)際生產(chǎn)中，為提高四川地區(qū)西門(mén)塔爾牛的產(chǎn)奶性能以及我國(guó)西門(mén)塔爾牛的選育提供理論支撐，對(duì)我國(guó)的奶業(yè)發(fā)展具有一定意義。