葛叢叢,劉華云,李寒妹,付艷芳,付志新,劉錚鑄,鞏元芳*
(1.河北科技師范學(xué)院 動(dòng)物科技學(xué)院,河北省特色動(dòng)物種質(zhì)資源挖掘與創(chuàng)新重點(diǎn)實(shí)驗(yàn)室,河北 秦皇島 066004;2.河北省畜牧總站,河北 石家莊 050035)
家畜的體重和體尺指標(biāo)與其生產(chǎn)性能、繁殖性能等密切相關(guān),并隨生態(tài)環(huán)境和選育方式的變更而發(fā)生改變[1]。絨山羊是經(jīng)過(guò)長(zhǎng)期自然選擇和人工選育而形成的一類絨毛用山羊品種,因其采食消化能力強(qiáng),纖維消化率高,可在貧瘠的山區(qū)放牧下生存[2],受到養(yǎng)殖者的喜愛(ài)。我國(guó)絨山羊主要分布在華北、西北、東北及青藏高原等地區(qū),主要包括遼寧絨山羊、內(nèi)蒙古絨山羊和河西絨山羊等品種[3]。體重、產(chǎn)絨量是直接影響絨山羊經(jīng)濟(jì)效益的重要指標(biāo)[4-5]。燕山絨山羊主產(chǎn)區(qū)在河北燕山山區(qū),多在秦皇島的青龍縣、承德的寬城縣等地飼養(yǎng)[6],是當(dāng)?shù)剞r(nóng)民以本地羊?yàn)榛A(chǔ),經(jīng)過(guò)數(shù)十年不太系統(tǒng)的選育和雜交,形成的一個(gè)絨肉兼用型的絨山羊新類群[7],截止到目前為止,成年羊群體內(nèi)仍存在體型大小、產(chǎn)絨量及絨品質(zhì)等性能參差不齊的問(wèn)題,嚴(yán)重影響其經(jīng)濟(jì)效益。因此,加快燕山山區(qū)優(yōu)良絨山羊的選育,提高產(chǎn)絨量和羊絨品質(zhì)勢(shì)在必行。本研究擬對(duì)燕山絨山羊的體尺指標(biāo)、體重、絨厚與產(chǎn)絨量等指標(biāo)進(jìn)行測(cè)定與分析,建立最優(yōu)回歸方程,旨在找出燕山絨山羊各體尺指標(biāo)、體重、絨厚與產(chǎn)絨量等的相互關(guān)系,為燕山絨山羊的選育和進(jìn)一步開(kāi)發(fā)、利用提供參考依據(jù)。
2019年3~4月份,在河北省秦皇島市青龍滿族自治縣某絨山羊養(yǎng)殖場(chǎng)選取年齡在1~5歲健康無(wú)病的燕山絨山羊母羊162只,采集體尺、抓絨后體重、絨厚和產(chǎn)絨量等相關(guān)數(shù)據(jù)。
利用皮尺、鋼尺、測(cè)杖、電子秤等工具對(duì)受試山羊的產(chǎn)絨量(Y)、體長(zhǎng)(X1)、體高(X2)、胸圍(X3)、管圍(X4)、絨厚(X5)和抓絨后的體重(X6)進(jìn)行了測(cè)定,具體操作方法見(jiàn)文獻(xiàn)[8]。
利用Excel和SPSS20.0軟件進(jìn)行簡(jiǎn)單歸納處理,剔除無(wú)效數(shù)據(jù)后,對(duì)各指標(biāo)進(jìn)行相關(guān)系數(shù)、回歸方程、通徑系數(shù)等參數(shù)的分析。
利用SPSS20.0軟件對(duì)所得燕山絨山羊體尺、抓絨后體重、絨厚和產(chǎn)絨量的數(shù)據(jù)進(jìn)行了描述統(tǒng)計(jì),結(jié)果見(jiàn)表1。由表1可知,燕山絨山羊體長(zhǎng)和體高的變異系數(shù)較小,為弱變異(CV<10%)[9],而胸圍、管圍、絨厚、抓絨后體重和產(chǎn)絨量為中等變異(10% 表1 燕山絨山羊各體尺指標(biāo)、抓絨后體重、絨厚與產(chǎn)絨量的描述性統(tǒng)計(jì)Table 1 Descriptive statistics of body size,body weight after fleece,cashmere thickness and cashmere yield of Yanshan Cashmere Goat 相關(guān)分析可以確定兩性狀之間的相關(guān)程度,用相關(guān)系數(shù)r來(lái)表示,r為正且值越接近1,表示兩性狀之間的正相關(guān)程度越強(qiáng)。利用SPSS20.0軟件統(tǒng)計(jì)分析了燕山絨山羊各體尺指標(biāo)、抓絨后體重、絨厚與產(chǎn)絨量之間的相關(guān)程度,獲得的相關(guān)系數(shù)見(jiàn)表2。 由表2可知,體長(zhǎng)、體高、胸圍、絨厚、抓絨后體重與產(chǎn)絨量均呈極顯著正相關(guān)(P<0.01),其中絨厚與產(chǎn)絨量的相關(guān)系數(shù)最大,為0.728;體長(zhǎng)、胸圍、抓絨后體重與產(chǎn)絨量相關(guān)系數(shù)較大,分別為0.449、0.639和0.558;體高與產(chǎn)絨量相關(guān)系數(shù)較小,為0.274。管圍與產(chǎn)絨量相關(guān)不顯著(P>0.05),相關(guān)系數(shù)僅為0.071。體高、管圍與絨厚相關(guān)性不顯著(P>0.05),其余各性狀之間均為極顯著正相關(guān)(P<0.01)。 表2 燕山絨山羊體尺、抓絨后體重、絨厚和產(chǎn)絨量的相關(guān)分析結(jié)果Table 2 Correlation analysis results of body size,body weight after fleece,cashmere thickness and cashmere yield of Yanshan Cashmere Goat 根據(jù)以上各指標(biāo)相關(guān)系數(shù)的組成效應(yīng),可將燕山絨山羊各體尺指標(biāo)、抓絨后體重、絨厚中某一性狀(Xi)與產(chǎn)絨量(Y)的相關(guān)系數(shù)剖分為某一性狀(Xi)對(duì)產(chǎn)絨量的直接影響(Piy)和某一性狀(Xi)通過(guò)另一性狀(Xj)對(duì)產(chǎn)絨量的間接影響(∑rijPj)兩部分,即相關(guān)系數(shù)riy=Xi與Y的直接通徑系數(shù)(Piy)+所有Xi與Y的間接通徑系數(shù)(∑rijPj),Xi與Y的間接通徑系數(shù)=相關(guān)系數(shù)(rij)×通徑系數(shù)(Pjy)[10]。利用SPSS20.0軟件對(duì)各指標(biāo)與產(chǎn)絨量之間的關(guān)系進(jìn)行建模,獲得相關(guān)系數(shù)和通徑系數(shù),所得結(jié)果見(jiàn)表3。 由表3各指標(biāo)相關(guān)系數(shù)的分析可知,體長(zhǎng)、體高、胸圍、管圍、抓絨后體重、絨厚對(duì)產(chǎn)絨量的直接影響中,絨厚對(duì)產(chǎn)絨量的直接作用最大,為0.547,其次是胸圍、體長(zhǎng)對(duì)產(chǎn)絨量的直接作用,分別為0.194和0.138,而體高、管圍、抓絨后體重對(duì)產(chǎn)絨量的直接作用均較小。間接影響中,抓絨后體重對(duì)產(chǎn)絨量的間接作用最大,為0.461,其次是胸圍,為0.447,管圍對(duì)產(chǎn)絨量的間接作用最小,為0.150。綜上,絨厚、胸圍、體長(zhǎng)主要通過(guò)直接作用影響產(chǎn)絨量,體高、管圍、抓絨后體重主要通過(guò)胸圍對(duì)產(chǎn)絨量產(chǎn)生間接影響,其中抓絨后體重產(chǎn)生的間接影響最大。 表3 燕山絨山羊各體尺指標(biāo)、抓絨后體重、絨厚與產(chǎn)絨量相關(guān)系數(shù)的分析Table 3 Analysis of correlation coefficients of body size,body weight after fleece,cashmere thickness and cashmere yield of Yanshan Cashmere Goat 利用逐步法進(jìn)行回歸分析,得到燕山絨山羊各體尺指標(biāo)、抓絨后體重、絨厚與產(chǎn)絨量關(guān)系的回歸方程。根據(jù)顯著性檢驗(yàn),選擇回歸關(guān)系極顯著(P<0.01)的最優(yōu)回歸方程為:Y=-598.036+4.527X1+4.81X3+66.282X5,將此回歸方程標(biāo)準(zhǔn)化后為:Y=0.161X1+0.244X3+0.555X5。 標(biāo)準(zhǔn)化回歸方程內(nèi)部各指標(biāo)經(jīng)顯著性分析可知,胸圍(X3)、絨厚(X5)與產(chǎn)絨量(Y)的回歸關(guān)系呈極顯著水平(P<0.01),體長(zhǎng)(X1)與產(chǎn)絨量(Y)的回歸關(guān)系呈顯著水平(P<0.05),說(shuō)明這些性狀對(duì)產(chǎn)絨量的關(guān)系是真實(shí)存在的。 絨山羊的產(chǎn)絨量和體重是衡量其經(jīng)濟(jì)效益的重要指標(biāo),也是養(yǎng)殖者選育的首要指標(biāo),二者受遺傳、營(yíng)養(yǎng)和環(huán)境等多方面的影響[11]。在選育過(guò)程中,既要考慮產(chǎn)絨量和體重等重要經(jīng)濟(jì)性狀,更要考慮各體尺指標(biāo)對(duì)產(chǎn)絨量的影響,達(dá)到精確選育的目的?!吨袊?guó)畜禽遺傳資源志-羊志》記錄遼寧絨山羊母羊的產(chǎn)絨量是641 g[12],明顯高于本研究所得燕山絨山羊的產(chǎn)絨量(548.58 g),說(shuō)明燕山絨山羊的產(chǎn)絨量仍有提高的空間。另外,本研究測(cè)得燕山絨山羊抓絨后體重和產(chǎn)絨量的變異系數(shù)分別為25.20 %和33.30 %,說(shuō)明采集的個(gè)體間指標(biāo)數(shù)據(jù)變動(dòng)幅度較大,群體內(nèi)絨山羊個(gè)體之間生長(zhǎng)發(fā)育差異較大。其中,產(chǎn)絨量的變異系數(shù)大于各體尺指標(biāo)的變異系數(shù),這與依明·蘇來(lái)曼等[13]對(duì)青格里絨山羊的研究和李景玉等[14]對(duì)遼寧絨山羊的研究結(jié)果基本一致,說(shuō)明燕山絨山羊產(chǎn)絨量的群體增長(zhǎng)趨勢(shì)與其他絨山羊一致。燕山絨山羊作為一種絨肉兼用型品種,抓絨后體重、產(chǎn)絨量變異系數(shù)較大,也為其群體內(nèi)選育提供可能。 通過(guò)相關(guān)分析和通徑分析發(fā)現(xiàn),燕山絨山羊的體長(zhǎng)、體高、胸圍、絨厚、抓絨后體重與產(chǎn)絨量均呈極顯著正相關(guān)(P<0.01),這與前人對(duì)陜北絨山羊[15-16]和遼寧絨山羊[17]的研究結(jié)果一致。本研究發(fā)現(xiàn),燕山絨山羊絨厚(0.547)、胸圍(0.194)、體長(zhǎng)(0.138)主要通過(guò)直接作用對(duì)產(chǎn)絨量產(chǎn)生影響;間接作用最大的是抓絨后體重(0.461),抓絨后體重主要通過(guò)影響胸圍(0.158)和絨厚(0.222)等指標(biāo)間接影響產(chǎn)絨量。管圍對(duì)產(chǎn)絨量的直接作用和間接作用均較小。綜上,燕山絨山羊的體長(zhǎng)、胸圍和絨厚對(duì)產(chǎn)絨量的直接作用和間接作用都較強(qiáng)。雖然尚未有研究證實(shí)絨厚是產(chǎn)絨量的第一影響因子,但胸圍對(duì)產(chǎn)絨量的直接作用在藏山羊[18-19]、青格里絨山羊等[13]多個(gè)絨山羊品種上均得到證實(shí)。說(shuō)明在選育高產(chǎn)絨量的燕山絨山羊時(shí),除了參考絨厚這一指標(biāo)外,胸圍也應(yīng)作為一個(gè)重要的參考指標(biāo)。 由回歸分析可知,雖然體高、管圍、抓絨后體重也通過(guò)絨厚、胸圍和體長(zhǎng)對(duì)產(chǎn)絨量產(chǎn)生了間接影響,但在建立最優(yōu)線性回歸方程時(shí)發(fā)現(xiàn),僅絨厚、胸圍與產(chǎn)絨量的回歸關(guān)系呈極顯著水平(P<0.01),體長(zhǎng)與產(chǎn)絨量的回歸關(guān)系呈顯著水平(P<0.05),體高、管圍、抓絨后體重與產(chǎn)絨量關(guān)系不顯著(P>0.05),因此并未計(jì)入回歸方程中。體高未計(jì)入回歸方程,與賀曉宏等[17]對(duì)遼寧絨山羊的研究結(jié)果一致;抓絨后體重、管圍未計(jì)入回歸方程,與李景玉等[14]和閆海龍等[20]對(duì)遼寧絨山羊和陜北白絨山羊上的研究結(jié)果不一致,絨山羊的體尺指標(biāo)只有胸圍和體長(zhǎng)進(jìn)入了回歸方程,導(dǎo)致這一結(jié)果可能是絨山羊品種之間存在差異,或者是與各指標(biāo)測(cè)量的樣本量不足夠多有關(guān)。后續(xù)將繼續(xù)擴(kuò)大樣本量,多方面尋找影響燕山絨山羊經(jīng)濟(jì)性狀的影響因素,加快其選育的進(jìn)程。 本研究結(jié)果表明,燕山絨山羊體長(zhǎng)、體高、胸圍、絨厚和抓絨后體重與產(chǎn)絨量均呈極顯著正相關(guān)。體長(zhǎng)、胸圍和絨厚對(duì)燕山絨山羊的直接作用和間接作用都較強(qiáng)。燕山絨山羊體尺與產(chǎn)絨量的最優(yōu)線性回歸方程標(biāo)準(zhǔn)化后為Y=0.161X1(體長(zhǎng))+0.244X3(胸圍)+0.555X5(絨厚)。2.2 燕山絨山羊各指標(biāo)的相關(guān)分析
2.3 燕山絨山羊各指標(biāo)的通徑分析
2.4 燕山絨山羊各指標(biāo)與產(chǎn)絨量的回歸分析
3 討 論
4 結(jié) 論