貴州黑山羊體重與體尺性狀的相關性研究

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

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

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

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

利用Excel 2019 整理收集到的1729 只貴州黑山羊的體重與體尺數(shù)據(jù)，分析群體表型統(tǒng)計量。利用Pearson 法分析體重與各體尺之間的相關系數(shù)。按照年齡大小對群體進行分類，以體重性狀（Y）為因變量，體高（X1）、體長（X2）、胸深（X3）、胸寬（X4）、胸圍（X5）、管圍性狀（X6）為自變量進行通徑分析，并分解相關系數(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 體重與體尺性狀的相關性分析

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

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

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

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

為了更好地揭示不同年齡階段貴州黑山羊的體重與體尺性狀之間的關系，本研究將育成羊和成年羊群體分開，分別研究體重與體尺性狀間的關系。首先，利用體重與各體尺性狀的表型相關系數(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 體重與體尺性狀相關系數(shù)的分解

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

由表4 可知，在貴州黑山羊育成羊群體中，體長（X2）和胸圍（X5）對體重的相關系數(shù)最高且對體重的直接作用最強（P=0.541；P=0.518），而其余體尺性狀的直接作用較弱。其中，胸深（X3）和胸寬性狀（X4）對體重的直接作用最?。≒=-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 貴州黑山羊育成羊體重與體尺相關系數(shù)的分解

在貴州黑山羊成年羊群體中，胸圍（X5）對體重的相關系數(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 貴州黑山羊成年羊體重與體尺相關系數(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)計學意義（表8～9）。

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

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

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

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

3 討論

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

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

通徑系數(shù)是標準化的偏回歸系數(shù)，可以直接反映兩個變量之間的相關程度，分為直接通徑系數(shù)與間接通徑系數(shù)[9]。本研究通過分解貴州黑山羊育成羊與成年羊的相關系數(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ù)最高，說明胸圍性狀可能是影響成年山羊體重性狀的主要因素。綜上，育成羊與成年羊均引入體長與胸圍性狀進入回歸方程，說明體長與胸圍性狀對山羊生長中體重的變化有著至關重要的作用。在楊燕等[10]、白雅琴等[11]、陸會寧等[12]的研究中，均篩選出胸圍性狀為影響山羊體重的重要體尺因素，揭示了胸圍對體重影響的重要性。本研究提出的回歸方程R2均大于0.8，說明此回歸方程的回歸平方和占總變異量的80%以上，引入回歸模型的變量有較強的作用，預測的準確值較大，可為實際生產(chǎn)中貴州黑山羊的體重估計提供基礎模型，也為地方良種的選育與改良工作提供實踐指導。

4 結(jié)論

貴州黑山羊體重與各體尺之間存在極顯著正相關的關系。按照不同年齡階段構(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），為貴州黑山羊的選育工作提供理論指導及實踐依據(jù)。