廖圓圓,王燕新,劉翊中,李海健,徐紅偉,楊具田,阿依木古麗,蔡 勇?

(1 西北民族大學(xué)生命科學(xué)與工程學(xué)院;2 西北民族大學(xué)實驗教學(xué)部;3 西北民族大學(xué)生物醫(yī)藥研究中心,蘭州730030)

寧夏鹽池灘羊?qū)儆诿F的裘皮用綿羊品種,是從蒙古羊中分化出來的,在長期的自然選擇及牧民精心培育下養(yǎng)成。 鹽池灘羊生產(chǎn)的二毛裘皮是寧夏“五寶”之一[1-2],也是鹽池“三寶”之首。 鹽池縣因其得天獨厚的地理環(huán)境,培養(yǎng)出的灘羊具有獨特的品質(zhì)。 近年來,關(guān)于寧夏鹽池灘羊的生產(chǎn)性能、布魯氏菌病檢測的研究以及灘羊的遺傳特性的報道較多[3-5],但鮮見寧夏鹽池灘羊其表型性狀中的體重與各個體尺數(shù)據(jù)的相關(guān)性研究。 體重、體尺不僅反映羊的體格大小,而且能看出其營養(yǎng)狀況。 因此,本研究擬對寧夏鹽池灘羊的體重與各體尺指標(biāo)進(jìn)行分析,建立寧夏鹽池灘羊體重與各體尺間的最優(yōu)回歸方程,以期為寧夏鹽池灘羊后續(xù)的選種、選育工作提供參考。

1 材料與方法

1.1 實驗動物

實驗動物為隨機(jī)選擇健康、生長發(fā)育正常的寧夏鹽池灘羊745 只(1—1.5 歲,母羊400 只,公羊345 只,由寧夏鹽池縣鑫海食品有限公司提供)。 參考家畜育種學(xué)教材中對體尺指標(biāo)的測量方法[6],對745 只羊的體斜長(X1)、體高(X2)、胸寬(X3)、胸深(X4)、胸圍(X5)、管圍(X6)、十字部高(X7)、體重(Y)進(jìn)行測定。 灘羊的體重使用精度為0.01 kg 的電子秤進(jìn)行測量,各個體尺數(shù)據(jù)則使用測杖、卷尺。

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

所測的數(shù)據(jù)先使用Excel 2016 初步進(jìn)行整理,之后采用SPSS 22.0 統(tǒng)計分析的軟件進(jìn)行各種數(shù)據(jù)的分析,同時對各個相關(guān)的系數(shù)進(jìn)行分解,最終建立起寧夏鹽池灘羊的體重與各體尺性狀之間的最優(yōu)回歸方程。

2 結(jié)果與分析

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

由表1 可知,體重(Y)的變異系數(shù)最大,公、母羊分別為18.05%和17.24%;其次為母灘羊的胸寬(X3)和胸深(X4),分別為11.53%和10.52%,為中等變異;其他體尺指標(biāo)變異系數(shù)均小于10%,為弱變異。

表1 寧夏鹽池灘羊體重與體尺的表型統(tǒng)計量Table 1 Phenotype statistics of body weight and body size of Ningxia Yanchi tan sheep

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

采用SPSS 22.0 統(tǒng)計分析軟件對數(shù)據(jù)進(jìn)行分析,得出體重與各體尺性狀之間的相關(guān)系數(shù)。 根據(jù)寧夏鹽池灘羊體重與體尺的表型相關(guān)系數(shù)(表2)可知,體重(Y)與各體尺性狀之間呈極顯著正相關(guān)。 母灘羊的相關(guān)程度由小到大為:胸寬(X3) ＜管圍(X6) ＜體斜長(X1) ＜體高(X2) ＜十字部高(X7) ＜胸深(X4) ＜胸圍(X5);公灘羊的相關(guān)程度由小到大為:胸深(X4) ＜體斜長(X1) ＜體高(X2) ＜胸寬(X3) ＜管圍(X6) ＜十字部高(X7) ＜胸圍(X5)。 體重(Y)與胸圍(X5)的表型相關(guān)系數(shù)最大,分別為0.776和0.745。

表2 寧夏鹽池灘羊體重與體尺的表型相關(guān)系數(shù)Table 2 Phenotypic correlation coefficient between body weight and body size of Ningxia Yanchi tan sheep

2.3 體重與體尺的通徑分析

根據(jù)表1、表2 的初步分析后,對寧夏鹽池灘羊的體重與各個體尺的數(shù)據(jù)進(jìn)行通徑分析,最終確立以通徑系數(shù)(P)為變量的多元一次方程組。

母灘羊的方程組為:

對多元一次方程組進(jìn)行求解,得到通徑系數(shù)。

母灘羊:P1=0.041,P2= -0.169,P3=0.036,P4=0.139,P5=0.518,P6= -0.009,P7=0.370。

公灘羊:P1=0.212,P2= -0.130,P3= -0.101,P4= -0.157,P5=0.703,P6=0.010,P7=0.327。

2.4 體尺指標(biāo)對體重相關(guān)系數(shù)的分解

根據(jù)各個相關(guān)系數(shù)間的組合效應(yīng),可將寧夏鹽池灘羊各體尺性狀與體重(Y)的表型相關(guān)系數(shù)分解,一是各個體尺性狀對體重產(chǎn)生直接作用,二是各個體尺性狀通過作用其他體尺性狀從而對體重產(chǎn)生間接的影響。 根據(jù)表3、表4 可知,母灘羊各體尺指標(biāo)對體重的直接作用從小到大為:體高(X2) ＜管圍(X6) ＜胸寬(X3) ＜體斜長(X1) ＜胸深(X4) ＜十字部高(X7) ＜胸圍(X5)。 各體尺性狀對體重的間接影響從小到大為:胸圍(X5) ＜十字部高(X7) ＜胸寬(X3) ＜體斜長(X1) ＜管圍(X6) ＜胸深(X4) ＜體高(X2)。 公灘羊各體尺指標(biāo)對體重的直接作用從小到大為:胸深(X4) ＜體高(X2) ＜胸寬(X3) ＜管圍(X6) ＜體斜長(X1) ＜十字部高(X7) ＜胸圍(X5)。 各體尺指標(biāo)對體重的間接影響從小到大為:胸圍(X5) ＜體斜長(X1) ＜十字部高(X7) ＜管圍(X6) ＜胸深(X4) ＜體高(X2) ＜胸寬(X3)。 根據(jù)以上結(jié)果可以看出,胸圍(X5)、十字部高(X7)對體重有較強(qiáng)的直接作用;而體高(X2)、胸深(X4)、胸寬(X3)主要是通過間接作用對體重產(chǎn)生影響,對體重的直接作用較弱。

表3 寧夏鹽池母灘羊體尺對體重相關(guān)系數(shù)的分解Table 3 Decomposition of body size to body weight correlation coefficient of Ningxia Yanchi female tan sheep

表4 寧夏鹽池公灘羊體尺對體重相關(guān)系數(shù)的分解Table 4 Decomposition of body size to body weight correlation coefficient of Ningxia Yanchi male tan sheep

2.5 體重和體尺的逐步回歸分析

運用SPSS 22.0 軟件對數(shù)據(jù)進(jìn)行多元逐步回歸分析。 數(shù)據(jù)的因變量是體重(Y),自變量則分別是體斜長(X1)、體高(X2)、胸寬(X3)、胸深(X4)、胸圍(X5)、管圍(X6)、十字部高(X7)。 從表5—7 可以看出,管圍(X6)、胸寬(X3)對體重的作用效果不是很大,故將這兩個體尺性狀剔除。 母灘羊所得到的體重和各體尺性狀之間的最優(yōu)回歸方程為Y= -48.167 +0.549X5+0.444X7+0.380X4(R=0.832,P＜0.01),而方差分析所的F值為297.520,為極顯著水平(P＜0.01)(表7)。 公灘羊所得到的體重和各體尺性狀之間的最優(yōu)回歸方程為Y= -63.541 +0.89X5+0.246X7+0.178X1(R=0.784,P＜0.01),方差分析的F值為181.126,為極顯著水平(P＜0.01)(表7)。

表5 寧夏鹽池灘羊體尺指標(biāo)對體重回歸模型匯總Table 5 Body size indexes to body weight regression model summary of Ningxia Yanchi tan sheep

表6 寧夏鹽池灘羊多元回歸系數(shù)的參數(shù)估計值檢驗Table 6 Parameter estimation test of multivariate regression coefficient of Ningxia Yanchi tan sheep

表7 寧夏鹽池母灘羊多元回歸方差分析Table 7 Analysis of multivariate regression variance of Ningxia Yanchi tan sheep

3 討論與結(jié)論

體重、體尺作為綿羊的外貌體況評定的重要因素,與綿羊的各個生產(chǎn)性能均有極為密切的聯(lián)系。 變異系數(shù)反映的是其相對變化程度的指標(biāo)[7],當(dāng)其小于10%時為弱變異,大于10%而小于100%時為中等變異,研究結(jié)果顯示:鹽池灘羊體重(Y)的變異系數(shù)(公羊:18.05%,母羊:17.24%)最大,其次是母灘羊的胸寬(X3)、胸深(X4),分別為11.53%、10.52%。 其變異幅度與體重的相比則比較小,說明寧夏鹽池灘羊的體重(Y)具有較大的提高潛力。

皮爾遜(Pearson)相關(guān)系數(shù)的大小可以用來判斷變量間相關(guān)關(guān)系的密切程度[8]。 本研究通過相關(guān)性分析表明,體重(Y)與其他體尺均有顯著的相關(guān)性。 寧夏鹽池母灘羊體高(X2)對體重(Y)的綜合間接影響最大,為0.686,其次為胸深(X4),胸圍(X5)對體重(Y)的綜合間接影響最小,為0.258。 寧夏鹽池公灘羊胸寬(X3)對體重(Y)的綜合間接影響最大,為0.458,其次為體高(X2),胸圍(X5)對體重(Y)的綜合間接影響最小,為0.041。 根據(jù)上述結(jié)果可以看出:胸圍、十字部高這兩個體尺對體重的綜合作用效果較大,一是表現(xiàn)為直接作用,二是通過輔助其他性狀從而對灘羊體重產(chǎn)生的間接影響,其一或者二者皆有。 因此在后續(xù)的灘羊選育生產(chǎn)中,需進(jìn)行全面分析,以達(dá)到最優(yōu)的效果。

通徑系數(shù)表示相關(guān)變量之間的因果關(guān)系。 本研究表明,胸圍、十字部高對體重的直接作用效果明顯,而間接的影響效果較小。 母灘羊的胸深對體重的直接作用明顯,同時也有一些間接影響。 體高(X2)對體重(Y)的綜合間接影響最大,為0.686;胸圍(X5)對體重(Y)的綜合間接影響最小,為0.258。 公灘羊的體斜長對體重具有明顯的直接作用,又有一定間接影響。 胸寬(X3)對體重(Y)的綜合間接影響最大,為0.458;胸圍(X5)對體重(Y)的綜合間接作用最小,為0.041。 根據(jù)上述結(jié)果可以看出,胸圍(X5)、十字部高(X7)對體重的作用效果較強(qiáng);而體高(X2)、胸深(X4)則主要是通過間接影響從而對體重產(chǎn)生一定的作用,對體重的直接影響較小。 故在選種過程中,除了加強(qiáng)對體重影響直接作用強(qiáng)的因素(如寧夏鹽池灘羊胸圍、十字部高)的選擇,同時也應(yīng)關(guān)注其他體尺性狀對體重的影響。

通過體重和體尺的逐步回歸分析,得到母灘羊各體尺性狀對體重的最優(yōu)回歸方程為:Y= -48.167 +0.549X5+0.444X7+0.380X4(R=0.832,P＜0.01)。 公灘羊各體尺對體重的最優(yōu)回歸方程為:Y=-63.541 +0.89X5+0.246X7+0.178X1(R=0.784,P＜0.01)。R作為體重和各體尺指標(biāo)之間的復(fù)相關(guān)系數(shù),R2為體重和各體尺指標(biāo)之間的相關(guān)指數(shù)[9]。 回歸方程中R2分別為0.693、0.614,說明該回歸方程的回歸平方和占總變異量的69.3%及61.4%,所選的體尺指標(biāo)對體重具有較大的決定作用,具有一定的參考價值。 康建兵等[10]對貴州白山羊的最優(yōu)回歸模型入選了體高、體斜長、胸圍和管圍4 個性狀。 李勝開等[11]對隆林山羊的最優(yōu)回歸模型入選了體高、體長和胸圍3 個性狀。 王志明等[12]對甘肅高山細(xì)毛羊(母)的最優(yōu)回歸模型入選了胸圍和體長兩個性狀。 索效軍等[13]對烏骨山羊(母)的最優(yōu)回歸模型入選了體高、體長和胸圍3 個性狀。 楊燕等[14]對沂蒙黑山公羊的最優(yōu)回歸模型入選了胸圍、胸深和管圍3 個性狀。 白雅琴等[15]對岷縣黑裘皮公羊的最優(yōu)回歸模型入選了體高、體長、胸圍和胸寬4 個性狀,母羊的最優(yōu)回歸模型入選了體高、體長、胸圍3 個性狀,王志明[16]對甘肅高山細(xì)毛的最優(yōu)回歸模型入選了體高、體長2 個性狀,王偉峰[17]對白薩?？搜虻淖顑?yōu)回歸模型入選了胸圍、體斜長2 個性狀,陸會寧等[18]對蘭州大尾羊×小尾寒羊F1代的最優(yōu)回歸模型入選了體長和胸圍2 個性狀,吳錦波等[19]對三江牛的最優(yōu)回歸模型入選了體高、體長、胸圍3 個性狀。 本試驗的母灘羊最優(yōu)回歸模型入選了胸圍、管圍、體斜長和體高4 個性狀,公灘羊最優(yōu)回歸模型入選了胸圍、十字部高和體斜長3 個性狀。 結(jié)果存在一定的差異,可能是因為品種、飼養(yǎng)方式、人員的管理以及各自生活的環(huán)境等種種不同因素造成的。 母灘羊剩余因子e=1 -R2=0.307,公灘羊的剩余因子為0.386,說明對寧夏鹽池灘羊體重產(chǎn)生影響的自變量不止胸圍、十字部高等,還有其他影響因素需要進(jìn)一步分析。