莊蕾，閻明毅，余忠祥，張強(qiáng)龍，唐燕花，多杰才讓，唐燕青，李永欽，吳森

(1.青海大學(xué)畜牧獸醫(yī)科學(xué)院 青海省高原家畜遺傳資源保護(hù)與創(chuàng)新利用重點(diǎn)實(shí)驗(yàn)室，西寧810016;2.青海省河南蒙古族自治縣生態(tài)有機(jī)畜牧業(yè)技術(shù)服務(wù)中心，河南縣811500;3.青海省河南蒙古族自治縣有機(jī)產(chǎn)業(yè)科技園區(qū)管理委員會(huì)，河南縣811500;4.青海省河南蒙古族自治縣托葉瑪鄉(xiāng)農(nóng)牧業(yè)技術(shù)服務(wù)站，河南縣811500;5.青海省河南蒙古族自治縣獸醫(yī)站，河南縣811500)

歐拉羊是青藏高原藏系綿羊的一個(gè)重要類群，曾被稱為“歐拉型藏羊”，2018年被列入國(guó)家畜禽資源，更名為“歐拉羊”，因其主產(chǎn)區(qū)位于歐拉山而得名。歐拉羊具有發(fā)育快、體格結(jié)實(shí)、肢體高大、生長(zhǎng)快、繁殖性能好、適應(yīng)高寒氣候、耐粗飼、肉質(zhì)好等特點(diǎn)［1，2］，其成年羊體重在60 kg以上，遠(yuǎn)大于其他藏羊。另外，歐拉羊主產(chǎn)區(qū)位于“中華水塔”三江源腹地的河南蒙古族自治縣，擁有國(guó)家認(rèn)證的有機(jī)牧場(chǎng)，其肉質(zhì)獨(dú)特，被世人所關(guān)注［3］。歐拉羊具備較高的開(kāi)發(fā)利用價(jià)值，近年來(lái)，在本品種不斷選育的基礎(chǔ)上，現(xiàn)選育出大群無(wú)角型歐拉羊，性格溫順，易于管理，且符合現(xiàn)代規(guī)模化、集約化養(yǎng)殖趨勢(shì)及高效畜牧業(yè)生產(chǎn)理念，歐拉羊是良好的現(xiàn)代化藏羊品種育成資源。

羔羊期是綿羊生長(zhǎng)的關(guān)鍵時(shí)期，容易受代謝印記影響后期發(fā)育。羔羊期生長(zhǎng)發(fā)育相關(guān)研究有助于群體早期選種及生長(zhǎng)發(fā)育評(píng)估，特別是公羔羊。公羔羊的體況評(píng)估研究有利于種公羊的早期選留，減少不必要的養(yǎng)殖負(fù)擔(dān)，實(shí)現(xiàn)種群的快速優(yōu)化。本研究對(duì)無(wú)角型歐拉羊公羔羊的體重與體尺性狀進(jìn)行了線性回歸分析、通徑分析、逐步回歸分析，建立了最優(yōu)回歸方程，旨在為無(wú)角型歐拉羊的早期選種選育、優(yōu)良種公羊培育等工作提供技術(shù)和理論參考。

1 材料與方法

1.1 試驗(yàn)動(dòng)物

本研究以無(wú)角型歐拉羊?yàn)檠芯繉?duì)象，測(cè)定收集青海省黃南藏族自治州河南蒙古族自治縣優(yōu)干寧鎮(zhèn)經(jīng)多年本品種選育的無(wú)角型歐拉羊公羔羊(8～10月齡)653只，所有測(cè)定羊均發(fā)育正常、生長(zhǎng)健康，實(shí)行常規(guī)免疫，自然放牧。

1.2 測(cè)定指標(biāo)與方法

剪毛后對(duì)其體重(Y)和主要體尺性狀，即體高(X1)、體長(zhǎng)(X2)、胸寬(X3)、胸深(X4)、尻寬(X5)、胸圍(X6)進(jìn)行測(cè)定，生產(chǎn)性能測(cè)定主要參照張沅主編《家畜育種學(xué)》的第三章第七節(jié)體型鑒定［4］及《綿山羊生產(chǎn)性能測(cè)定技術(shù)規(guī)范》行業(yè)標(biāo)準(zhǔn)［5］進(jìn)行，由固定測(cè)量專員對(duì)體高、體長(zhǎng)、胸寬、胸深、尻寬、胸圍進(jìn)行測(cè)定，每個(gè)部位測(cè)定3次，取平均值。

1.3 逐步回歸分析

建立最優(yōu)回歸方程的方法很多，常用的是逐步回歸法(Stepwise)［6］，其原理是將自變量逐個(gè)引入，對(duì)全部的自變量X1，X2，...，Xp，按其對(duì)因變量Y影響的大小進(jìn)行比較，通過(guò)F檢驗(yàn)，選擇偏回歸平方和顯著的變量進(jìn)入回歸方程。同時(shí)，每引入一個(gè)新變量后，對(duì)前面已引入回歸模型的變量逐個(gè)進(jìn)行檢驗(yàn)，將經(jīng)檢驗(yàn)后偏回歸平方和不顯著的變量刪除，以保證所得自變量子集中每一個(gè)變量都是顯著的，經(jīng)過(guò)多次循環(huán)，直到不能再引入新變量為止。此時(shí)的多元線性回歸方程即為最優(yōu)多元線性回歸方程［7，8］。

1.4 數(shù)據(jù)處理

用Excel進(jìn)行數(shù)據(jù)基本整理，用SPSS 23.0統(tǒng)計(jì)軟件進(jìn)行統(tǒng)計(jì)分析，對(duì)相關(guān)系數(shù)進(jìn)行分解，建立體重與各體尺性狀之間的最優(yōu)回歸方程。

2 結(jié)果與分析

2.1 無(wú)角型歐拉羊體重與體尺性狀的基本統(tǒng)計(jì)量

由表1可知，無(wú)角型歐拉羊體高變異系數(shù)是7.87%，為弱變異(CV＜10%)。其余表型均為中等變異(10%＜CV＜100%)，胸寬變異系數(shù)為57.36%，接下來(lái)依次是體長(zhǎng)43.66%、尻寬36.63%、體重23.93%、胸深12.36%、胸圍11.98%。

表1 無(wú)角型歐拉羊體重與體尺性狀的基本統(tǒng)計(jì)量

2.2 無(wú)角型歐拉羊體重與體尺性狀的相關(guān)性分析

由圖1可知，體重與體高、胸深、胸圍之間存在極顯著正相關(guān)(P＜0.01)，相關(guān)程度由大到小為體高＞胸圍＞胸深;體重與體長(zhǎng)、胸寬、尻寬之間存在顯著正相關(guān)(P＜0.05)，相關(guān)程度由大到小為體長(zhǎng)=胸寬＞尻寬;其余各性狀兩兩之間也存在正相關(guān)。

圖1 無(wú)角型歐拉羊體重與體尺性狀的相關(guān)性分析(Pearson)

2.3 無(wú)角型歐拉羊體尺指標(biāo)和體重相關(guān)系數(shù)的分解

對(duì)無(wú)角型歐拉羊的各體尺性狀與體重?cái)?shù)據(jù)進(jìn)行逐步回歸分析。根據(jù)相關(guān)系數(shù)的組成效應(yīng)，可將無(wú)角型歐拉羊羔羊的各體尺性狀(Xi)對(duì)體重(Y)的相關(guān)系數(shù)劃分為直接作用(Pi)和通過(guò)性狀(Xj)的間接作用(∑RijPj)兩部分，即rxiy=Pi+∑RijPj。由表2可知，體高對(duì)體重的直接作用最大，為0.421;胸深對(duì)體重的間接作用最大，為0.2034。

表2 無(wú)角型歐拉羊體高、體長(zhǎng)、胸寬、胸深、尻寬、胸圍和體重相關(guān)系數(shù)的分解

2.4 最優(yōu)回歸方程的建立

回歸方程1引入的體尺性狀是體高，回歸方程2引入的體尺性狀是體高和胸圍，回歸方程3引入的體尺性狀是體高、胸圍、胸深。由表3可知，體高、胸圍、胸深被依次引入回歸方程后，其R值逐漸變大，估計(jì)標(biāo)準(zhǔn)誤逐漸變小。

表3 無(wú)角型歐拉羊體重與體尺的逐步回歸模型與參數(shù)估計(jì)值檢驗(yàn)

由上述分析可知，回歸方程3的決定系數(shù)最大，所以回歸方程3為最優(yōu)回歸方程，復(fù)相關(guān)系數(shù)R為0.616，校正的決定系數(shù)R2為0.376?；貧w方程3中對(duì)體重影響最顯著的體尺性狀是體高、胸圍、胸深，各性狀標(biāo)準(zhǔn)偏回歸系數(shù)為0.425、0.253、0.122。經(jīng)T檢驗(yàn)，體高、胸圍、胸深均達(dá)到極顯著水平(P＜0.01)。綜上所述，逐步回歸方程Y=0.425X1+0.253X6+0.122X4為最優(yōu)回歸方程，但Y=0.544X1可作為野外無(wú)稱重條件下快速預(yù)測(cè)歐拉羊體重的最優(yōu)回歸方程。

3 討論

3.1 體重與體尺的表型分析

體重和體尺性狀是綿羊外貌體況評(píng)定的重要指標(biāo)，不僅反映綿羊品種自身的體型特征，同時(shí)與綿羊的生產(chǎn)性能有密切聯(lián)系［9，10］。變異系數(shù)可以將不同單位觀測(cè)數(shù)據(jù)的相對(duì)變化程度直觀的表現(xiàn)出來(lái)，變異系數(shù)低于10%時(shí)稱為弱變異，在10%～100%之間稱為中等變異［7］。在本研究中，無(wú)角型歐拉羊的體高變異系數(shù)是7.87%，為弱變異;體重及其余體尺性狀變異系數(shù)均處于10%～100%之間，為中等變異，且各性狀變異系數(shù)均較大，說(shuō)明無(wú)角型歐拉羊公羔羊在體重、體尺性狀方面具有較大的選育潛力。

3.2 體重與體尺性狀間的相關(guān)性分析

體重是重要的經(jīng)濟(jì)與選種指標(biāo)，主要受到遺傳和環(huán)境的共同影響，因此在同品種內(nèi)和同環(huán)境下，研究體重與各體尺性狀之間的關(guān)系具有一定的代表性［11］。無(wú)角型歐拉羊的體高、胸深、胸圍與體重之間存在極顯著正相關(guān)，體長(zhǎng)、胸寬、尻寬與體重之間存在顯著正相關(guān)，且體高、胸深、胸圍對(duì)體重均有直接或間接作用，說(shuō)明這3個(gè)體尺性狀對(duì)無(wú)角型歐拉羊在體重增長(zhǎng)方面影響較大。在對(duì)無(wú)角型歐拉羊體重增長(zhǎng)影響最大的3個(gè)體尺性狀中，由相關(guān)系數(shù)可知，體高對(duì)體重的影響最大，為0.544;其次是胸圍，為0.429;胸深對(duì)體重的影響最小，為0.324。體重與體高、胸圍、胸深之間的相關(guān)是由于體重與動(dòng)物的肌肉、骨骼、內(nèi)臟的生長(zhǎng)發(fā)育緊密相關(guān)［12］。在梁小鵬等［13］對(duì)新疆也木勒白羊、許鑫等［14］對(duì)特克塞爾羊的研究中發(fā)現(xiàn)，不僅胸圍對(duì)體重會(huì)有較大影響，還有體長(zhǎng)等體尺性狀也會(huì)對(duì)體重有較大影響，結(jié)果與本研究不同，可能是本文測(cè)定的無(wú)角型歐拉羊的月齡較小和品種差異造成的。白俊艷等［15］對(duì)大尾寒羊的研究中發(fā)現(xiàn)，體高和胸圍不僅對(duì)體重有較強(qiáng)的直接作用，還輔助其他性狀產(chǎn)生較大的間接作用，與本研究結(jié)果相似。在李江文等［16］對(duì)呼倫貝爾短尾羊、海占忠等［17］對(duì)加什科羊、楊卓等［18］對(duì)青海加什科羊等的研究中也有相似的結(jié)果。因此，今后在無(wú)角型歐拉羊的品種選育、提高生產(chǎn)性能等研究中，應(yīng)以體重為主，兼顧體高、胸圍、胸深，以取得較好的選育效果。

3.3 最優(yōu)回歸方程的建立

體尺性狀不僅反映其外貌的特征，同時(shí)也可以為體重的預(yù)估提供參考［17］。本研究以體重為因變量，體尺性狀為自變量建立逐步回歸方程:Y=0.425X1+0.253X6+0.122X4，因?yàn)轶w高(X1)、胸深(X4)、胸圍(X6)對(duì)體重影響較大，所以被選入回歸方程，擬合優(yōu)度的決定系數(shù)為0.379，校正R2值為0.376，回歸關(guān)系達(dá)到極顯著水平(P＜0.01)，表明應(yīng)用該回歸方程能夠較準(zhǔn)確地評(píng)估無(wú)角型歐拉羊的體重變化。

綜上所述，體高、胸深、胸圍是影響無(wú)角型歐拉羊公羔羊體重的主要因素，無(wú)角型歐拉羊公羔羊的體重-體尺最優(yōu)回歸方程為Y=0.425X1+0.253X6+0.122X4為最優(yōu)回歸方程，Y=0.544X1可作為野外無(wú)稱重條件下快速預(yù)測(cè)歐拉羊體重的最優(yōu)回歸方程。

4 結(jié)論

無(wú)角型歐拉羊公羔羊各體尺性狀兩兩之間存在正相關(guān)，體高、胸深、胸圍與體重之間存在極顯著正相關(guān)，這3個(gè)體尺性狀是影響體重的主要因素，因此在對(duì)無(wú)角型歐拉羊進(jìn)行選育時(shí)，要重點(diǎn)關(guān)注體高、胸深、胸圍3個(gè)性狀，以獲得更好的選育結(jié)果。