貴州黑山羊體重與體尺性狀的相關(guān)性研究

吳 萍，郭振剛*，吳 瑛，白洋洋，藍賢勇，宋德榮，龔婷婷，吳道全

（1.畢節(jié)市畜牧獸醫(yī)科學(xué)研究所，貴州畢節(jié) 551700；2.西北農(nóng)林科技大學(xué)動物科技學(xué)院，陜西西安 712100；3.貴州新烏蒙生態(tài)牧業(yè)發(fā)展有限公司，貴州赫章 553200；4.畢節(jié)市老科學(xué)技術(shù)工作者協(xié)會，貴州畢節(jié) 551700）

貴州黑山羊?qū)偃庥眯蜕窖?，是貴州省的地方特色種質(zhì)資源，已被列為貴州省地方優(yōu)良品種保護名錄。畢節(jié)市屬典型高寒山區(qū)，平均海拔為1950m，是貴州黑山羊的主產(chǎn)區(qū)。貴州黑山羊為適應(yīng)畢節(jié)高寒氣候環(huán)境，具有肉質(zhì)好、體質(zhì)結(jié)實、耐粗飼、抗逆性強等優(yōu)良特性[1]。研究表明，動物體重是動物生長發(fā)育、肉品質(zhì)和經(jīng)濟性能最直接的反映，因此，體重發(fā)育規(guī)律及影響因素已逐漸被人們重視。朱俊紅等[2]對努比亞山羊體重與體尺指標(biāo)的相關(guān)關(guān)系進行分析，建立了體重與體尺指標(biāo)之間的最優(yōu)回歸方程，Y=0.900X3（R2=0.811，P＜0.01）。齊麗娜等[3]采 用Von Bertalanffy、Gompertz、Logistic 3 種非線性模型對0～12 月齡湘東黑山羊的體重進行生長曲線擬合，發(fā)現(xiàn)3 種方法均可預(yù)測湘東黑山羊的生長發(fā)育規(guī)律。陶林等[4]采用通徑分析云上黑山羊體尺性狀對體重的影響，發(fā)現(xiàn)胸圍是影響云上黑山羊母羊體重的主要體尺因素；李麗娟等[5]通過對不同年齡階段貴州黑山羊的生長曲線擬合發(fā)現(xiàn)體重變化拐點，同時還利用三種不同動物模型計算了貴州黑山羊生長性狀的遺傳參數(shù)并進行遺傳力估計[6]。但作為與體重直接相關(guān)的體尺性狀，貴州黑山羊體重與體尺性狀之間的相關(guān)性分析尚未見報道。因此，本研究對貴州黑山羊體重與體尺性狀進行相關(guān)分析，并構(gòu)建體重與各體尺性狀之間的最優(yōu)回歸方程，以期為貴州黑山羊后續(xù)的選育工作提供理論依據(jù)。

1 材料與方法

1.1 試驗動物與數(shù)據(jù)收集

試驗羊1729 只，其中育成羊927 只，均為10～12 月齡，成年羊792 只，均為2～3 歲，羊只全部來自赫章縣貴州黑山羊保種區(qū)養(yǎng)殖戶，數(shù)據(jù)測定包括體重和體尺（體高、體長、胸深、胸寬、胸圍、管圍）。試驗羊均生長發(fā)育正常、體質(zhì)健康、生長環(huán)境相同、飼養(yǎng)管理一致，且試驗數(shù)據(jù)由同一個技術(shù)人員按照動物育種學(xué)中講述的體尺指標(biāo)測量方法測定，以減小數(shù)據(jù)誤差。

1.2 數(shù)據(jù)分析

利用Excel 2019 整理收集到的1729 只貴州黑山羊的體重與體尺數(shù)據(jù)，分析群體表型統(tǒng)計量。利用Pearson 法分析體重與各體尺之間的相關(guān)系數(shù)。按照年齡大小對群體進行分類，以體重性狀（Y）為因變量，體高（X1）、體長（X2）、胸深（X3）、胸寬（X4）、胸圍（X5）、管圍性狀（X6）為自變量進行通徑分析，并分解相關(guān)系數(shù)。最后，利用逐步回歸分析構(gòu)建回歸模型，建立貴州黑山羊體重和各體尺性狀間的最優(yōu)回歸方程。

2 結(jié)果與分析

2.1 體重與體尺性狀的表型統(tǒng)計量分析

由表1 可知，在體尺性狀方面，育成羊胸深的變異系數(shù)為10.96%，成年羊胸寬的變異系數(shù)為11.14%，均為中等變異，而其余體尺性狀均為弱變異；在體重上，育成羊的體重變異系數(shù)為18.02%，成年羊的體重變異系數(shù)為12.75%，均為中等變異。

表1 貴州黑山羊體重與體尺的表型統(tǒng)計量

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

利用Pearson 法對貴州黑山羊育成羊和成年羊體重和體尺性狀進行相關(guān)性分析。由表2 可知，貴州黑山羊成年羊與育成羊的體重（Y）與各體尺性狀之間均呈現(xiàn)顯著正相關(guān)的關(guān)系（P＜0.05）。在育成羊中，與體重的相關(guān)系數(shù)由大到小排列為體長（X2）＞胸圍（X5）＞胸深（X3）＞體高（X1）＞管圍（X6）＞胸寬（X4），體長（X2）和胸圍（X5）與體重（Y）的相關(guān)系數(shù)最高，為0.889 和0.882。由表3 可知，在成年羊中，與體重的相關(guān)系數(shù)由大到小排列為胸圍（X5）＞體長（X2）＞胸深（X3）＞管圍（X6）＞體高（X1）＞胸寬（X4），胸圍（X5）與體重（Y）的相關(guān)系數(shù)最高，為0.864。

表2 貴州黑山羊育成羊體重與體尺的表型相關(guān)系數(shù)

表3 貴州黑山羊成年羊體重與體尺的表型相關(guān)系數(shù)

2.3 體重與體尺性狀的通徑分析

為了更好地揭示不同年齡階段貴州黑山羊的體重與體尺性狀之間的關(guān)系，本研究將育成羊和成年羊群體分開，分別研究體重與體尺性狀間的關(guān)系。首先，利用體重與各體尺性狀的表型相關(guān)系數(shù)建立以通徑系數(shù)（P）為變量的多元一次方程組，計算直接通徑系數(shù)。

貴州黑山羊育成羊體尺性狀方程組如下：

貴州黑山羊成年羊體尺性狀方程組如下：

對多元一次方程組進行求解，得到如下通徑系數(shù)：

育成羊：P1=0.169；P2=0.541；P3=-0.240；P4=-0.015；P5=0.518；P6=0.034；

成年羊：P1=-0.088；P2=0.348；P3=0.203；P4=-0.085；P5=0.547；P6=0.049

2.4 體重與體尺性狀相關(guān)系數(shù)的分解

根據(jù)相關(guān)系數(shù)的組成效應(yīng)可以將相關(guān)系數(shù)分為直接作用與間接作用，為探究不同年齡階段貴州黑山羊各體尺對體重的直接與間接影響，本研究將體重與體尺的相關(guān)系數(shù)進行分解。

由表4 可知，在貴州黑山羊育成羊群體中，體長（X2）和胸圍（X5）對體重的相關(guān)系數(shù)最高且對體重的直接作用最強（P=0.541；P=0.518），而其余體尺性狀的直接作用較弱。其中，胸深（X3）和胸寬性狀（X4）對體重的直接作用最小（P=-0.240；P=-0.015），主要通過間接作用影響體重。各體尺性狀對體重的直接作用由大到小為：體長（X2）＞胸圍（X5）＞體高（X1）＞管圍（X6）＞胸寬（X4）＞胸深（X3），間接作用由大到小為：胸深（X3）＞體高（X1）＞胸寬（X4）＞管圍（X6）＞胸圍（X5）＞體長（X2）。以上結(jié)果表明，體長（X2）和胸圍性狀（X5）通過較強的直接作用影響體重，胸深（X3）、胸寬（X4）、體高（X1）和管圍（X6）主要通過與其他體尺性狀之間的間接作用影響體重。

表4 貴州黑山羊育成羊體重與體尺相關(guān)系數(shù)的分解

在貴州黑山羊成年羊群體中，胸圍（X5）對體重的相關(guān)系數(shù)最高且對體重的直接作用最強（P=0.547），而其余體尺性狀的直接作用較弱（表5）。其中，體高（X1）和胸寬（X4）性狀對體重的直接作用最小（P=-0.088；P=-0.049），主要通過間接作用影響體重。各體尺性狀對體重的直接作用由大到小為：胸圍（X5）＞體長（X2）＞胸深（X3）＞管圍（X6）＞胸寬（X4）＞體高（X1），間接作用由大到小為：胸深（X3）＞體高（X1）＞體長（X2）＞管圍（X6）＞胸寬（X4）＞胸圍（X5）（表5）。以上結(jié)果表明，胸圍性狀（X5）通過較強的直接作用影響體重，胸深（X3）和體高（X1）主要通過與其他體尺性狀之間的間接作用影響體重，而其他性狀對體重的影響較為平均。

表5 貴州黑山羊成年羊體重與體尺相關(guān)系數(shù)的分解

2.5 體尺性狀對體重的逐步回歸分析

運用SPSS 26.0 對貴州黑山羊體重與體尺性狀進行多元逐步回歸分析（表6～9）。試驗將貴州黑山羊育成羊和成年羊群體分離，分別計算最優(yōu)回歸方程。以體重（Y）為因變量，將體高（X1）、體長（X2）、胸圍（X3）、胸寬（X4）、胸圍（X5）和管圍（X6）6 個體尺性狀作為自變量逐步引入方程，采用步進的方法，選取顯著性強的性狀依次引入回歸方程，不顯著的性狀剔除，以確保引入回歸方程的變量都為顯著性變量且結(jié)果最優(yōu)。通過逐步回歸分析，貴州黑山羊育成羊群體的體重與各體尺之間的最優(yōu)回歸方程為：Y=0.606 X2+0.469X5-42.814（R2=0.931）（表6～7），貴州黑山羊成年羊群體的體重與各體尺之間的最優(yōu)回歸方程為：Y=0.745X5+0.409X2+0.459X3-65.322（R2=0.900），兩個回歸方程均達到極顯著水平（P＜0.01），說明上述模型具有統(tǒng)計學(xué)意義（表8～9）。

表6 貴州黑山羊育成羊體尺指標(biāo)對體重回歸模型匯總

表7 貴州黑山羊育成羊逐步回歸系數(shù)的參數(shù)估計值檢驗

表8 貴州黑山羊成年羊體尺指標(biāo)對體重回歸模型匯總

表9 貴州黑山羊成年羊逐步回歸系數(shù)的參數(shù)估計值檢驗

3 討論

貴州黑山羊作為貴州省重要的肉山羊品種，其優(yōu)良肉用性能廣泛受到人們青睞。但貴州黑山羊以放牧飼養(yǎng)為主，其體重信息較難測定，因此，研究體尺對體重的影響，對貴州黑山羊的體重估計及選育提高等具有重要意義。本研究將貴州黑山羊群體分為育成羊和成年羊進行平均值、標(biāo)準(zhǔn)差及變異系數(shù)等計算，發(fā)現(xiàn)體重性狀在育成羊（18.02%）和成年羊（12.75%）中均為中等變異，說明在不同年齡階段，山羊體重變化的相對程度更高，具有較大的選擇潛力，而其余體尺性狀則多為弱變異，提示體重性狀的選育潛力及程度往往更高。朱俊紅等[2]研究努比亞山羊體重與體尺的相關(guān)性分析，結(jié)果體重的變異系數(shù)最大為33.44%，其余體尺的變異系數(shù)較小，揭示體重具有較大的選擇空間；李勝開等[7]在隆林山羊品種上的研究也表明，體重的變異系數(shù)最大，為22.43%；上述研究與本研究結(jié)論一致，但是變異系數(shù)差異較大，這可能與品種差異、測量工具不同以及群體數(shù)量等有關(guān)，具體原因還需要進一步驗證。

Pearson 相關(guān)系數(shù)可以反映兩個線性變量間的相關(guān)程度，相關(guān)性系數(shù)越高的變量其相關(guān)程度更加密切[8]。本研究，利用Pearson 法對體重與體尺之間的相關(guān)系數(shù)進行檢驗，全部體尺與體重的相關(guān)性結(jié)果說明體尺性狀與體重存在顯著的線性正相關(guān)關(guān)系，側(cè)面反映了體尺性狀對體重的增長密不可分，其中，胸圍和胸深性狀與體重的相關(guān)系數(shù)最高，推測胸圍和胸深性狀更高的山羊具有更強的體重優(yōu)勢。李勝開等[7]研究表明，隆林山羊12月齡的體重與體高、體長、胸圍和管圍均呈極顯著正相關(guān)，與體重的相關(guān)性以體高（0.748）最大，該結(jié)論與本研究結(jié)果稍有差異，但是朱俊紅等研究，努比亞山羊的胸圍與體重的相關(guān)系數(shù)最大，達0.900，與本研究結(jié)論基本一致。

通徑系數(shù)是標(biāo)準(zhǔn)化的偏回歸系數(shù)，可以直接反映兩個變量之間的相關(guān)程度，分為直接通徑系數(shù)與間接通徑系數(shù)[9]。本研究通過分解貴州黑山羊育成羊與成年羊的相關(guān)系數(shù)，發(fā)現(xiàn)育成羊群體中體長與胸圍性狀對體重的直接作用較強，胸深的間接作用較強，加強體長和胸圍性狀的選育可能較好的影響育成體重；而成年羊群體中只有胸圍性狀對體重的直接作用較強，間接作用性狀分布較平均，說明加強胸圍性狀的選育可能較好的影響成年羊體重。通過對比貴州黑山羊育成羊和成年羊影響體重的直接因素和間接因素，胸圍性狀作為不同年齡階段影響體重的較大直接因素，對胸圍性狀的選育更有助于提高不同年齡階段山羊體重，但山羊生長過程中各體尺對體重的影響逐漸趨平均與穩(wěn)定，說明在選育體重性狀時還需考慮山羊其他間接影響因素的作用情況，以達到更好的選育效果。

通過對育成羊和成年羊體重與體尺性狀進行逐步回歸分析，得到育成羊體重與體尺的最優(yōu)回歸方程：Y=0.606X2+0.469X5-42.814（R2=0.931），方差分析F 值為111.298，達到極顯著水平。該方程引入胸圍和體長性狀進入回歸方程，說明胸圍與體長性狀更加顯著影響體重性狀，且體長的系數(shù)最高，說明體長性狀可能是影響育成山羊體重性狀的主要因素；成年體重與體尺最優(yōu)回歸方程：Y=0.745X5+0.409X2+0.459X3-65.322（R2=0.900），方差分析F 值為124.573，達到極顯著水平。該方程引入胸圍、體長和胸深性狀進入回歸方程，說明上述性狀更加顯著影響體重性狀，且體胸圍系數(shù)最高，說明胸圍性狀可能是影響成年山羊體重性狀的主要因素。綜上，育成羊與成年羊均引入體長與胸圍性狀進入回歸方程，說明體長與胸圍性狀對山羊生長中體重的變化有著至關(guān)重要的作用。在楊燕等[10]、白雅琴等[11]、陸會寧等[12]的研究中，均篩選出胸圍性狀為影響山羊體重的重要體尺因素，揭示了胸圍對體重影響的重要性。本研究提出的回歸方程R2均大于0.8，說明此回歸方程的回歸平方和占總變異量的80%以上，引入回歸模型的變量有較強的作用，預(yù)測的準(zhǔn)確值較大，可為實際生產(chǎn)中貴州黑山羊的體重估計提供基礎(chǔ)模型，也為地方良種的選育與改良工作提供實踐指導(dǎo)。

4 結(jié)論

貴州黑山羊體重與各體尺之間存在極顯著正相關(guān)的關(guān)系。按照不同年齡階段構(gòu)建回歸方程，育成羊的體長和胸圍性狀對體重的直接作用較強，其體重與各體尺間的最優(yōu)回歸方程為：Y=0.606X2+0.469X5-42.814（R2=0.931）；成年羊的胸圍性狀對體重的直接作用較強，其體重與各體尺間的最優(yōu)回歸方程為：Y=0.745X5+0.409X2+0.459X3-65.322（R2=0.900），為貴州黑山羊的選育工作提供理論指導(dǎo)及實踐依據(jù)。