施海娜， 劉雨田， 李世恩， 張金霞， 高 鈺， 梁萬鵬， 朱正生， 陳兆勇， 李 凱， 徐振飛*

(1.慶陽市農(nóng)業(yè)科學(xué)研究院，甘肅 慶陽 745000；2.甘肅中盛農(nóng)牧發(fā)展有限公司，甘肅 慶陽 745000)

“澳洲白”是澳大利亞第1個(gè)利用現(xiàn)代基因測定手段培育的品種。該品種集成了白杜泊綿羊、萬瑞綿羊、無解道賽特羊和特克賽爾羊等品種基因，通過對多個(gè)品種羊特定肌肉生長基因標(biāo)記和抗寄生蟲基因標(biāo)記的選擇(MyoMAX，LoinMAX，WormSTAR)，培育而成的專門用于與杜泊綿羊配套的，屬于粗毛型的中、大型肉羊品種，于2009年10月在澳大利亞注冊。其特點(diǎn)是體型大、生長快、成熟早、全年發(fā)情，有很好的自動(dòng)換毛能力。在放牧條件下5～6月齡可達(dá)到23 kg胴體，舍飼條件下，該品種6月齡胴體重可達(dá)26 kg，且脂肪覆蓋均勻，板皮質(zhì)量具佳。澳洲白在我國已多次引入，截至2017年甘肅省已引入發(fā)展數(shù)量達(dá)到2 000只以上[1]。此品種使養(yǎng)殖者能夠在各種養(yǎng)殖條件下用作三元配套的終端父本，可以產(chǎn)出在生長速率、個(gè)體重量、出肉率和出欄周期短等方面理想的商品羔羊。

家畜體尺數(shù)據(jù)直接反映了家畜的體格大小和體軀結(jié)構(gòu)、發(fā)育等情況，與家畜的生理機(jī)能、生產(chǎn)性能、抗病力以及對外界環(huán)境條件的適應(yīng)能力等密切相關(guān)。國內(nèi)外學(xué)者針對有關(guān)澳洲白羊種羊飼養(yǎng)管理技術(shù)和雜交改良地方品種試驗(yàn)等方面的研究較多[2-11]，但對其體質(zhì)量與體尺指標(biāo)相關(guān)性及回歸分析研究尚缺乏報(bào)道。本文通過對以湖羊?yàn)槭荏w采用胚胎移植技術(shù)擴(kuò)繁的澳洲白羊體質(zhì)量與其體尺指標(biāo)進(jìn)行測定，初步掌握澳洲白羊的生長發(fā)育情況，并進(jìn)行了相關(guān)性分析、通徑分析、逐步回歸分析以及對其表型相關(guān)系數(shù)分解，探究其體尺指標(biāo)與體質(zhì)量的關(guān)系，建立體質(zhì)量與其體尺指標(biāo)之間的最優(yōu)回歸方程，以期能為將來的肉羊新品種選育工作提供參考數(shù)據(jù)。

1 材料與方法

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

試驗(yàn)動(dòng)物澳洲白羊?yàn)?017年9月底在甘肅中盛華美羊業(yè)西峰羊場(以湖羊?yàn)槭荏w，從天津奧群公司購買的純種澳洲白A級冷凍胚胎)采用胚胎移植技術(shù)所產(chǎn)，年齡為6月齡左右，健康狀況良好。

1.2 測定項(xiàng)目

隨機(jī)抽取32只(公母羊各16只)健康羊只對其進(jìn)行體質(zhì)量(y)、體尺(xi)測定。體質(zhì)量稱量用電子秤，精度為0.01 kg。體尺指標(biāo)測量參考家畜育種學(xué)[12]，測量工具為測杖、卷尺，測量指標(biāo)包括體高(x1)、體斜長(x2)、胸圍(x3)、胸深(x4)、胸寬(x5)和管圍(x6)。

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

采用Excel 2007對所測數(shù)據(jù)進(jìn)行初步整理，用SPSS 19.0統(tǒng)計(jì)分析軟件進(jìn)行相關(guān)性分析、通徑分析、逐步回歸分析，并對相關(guān)系數(shù)進(jìn)行分解，建立體質(zhì)量與體尺指標(biāo)之間的最優(yōu)回歸方程。試驗(yàn)數(shù)據(jù)均以平均數(shù)±標(biāo)準(zhǔn)差表示。

2 結(jié)果與分析

2.1 體質(zhì)量與體尺指標(biāo)的表型統(tǒng)計(jì)量

表1為澳洲白羊體質(zhì)量與體尺指標(biāo)的表型統(tǒng)計(jì)量。由表1可見，體質(zhì)量(y)的變異系數(shù)最大，為10.0%；其次為胸寬(x5)、胸深(x4)、管圍(x6)、胸圍(x3)、體斜長(x2)，其變異系數(shù)分別為7.0%，6.1%，5.3%，4.8%，4.6%；體高(x1)的變異系數(shù)最小，為4.2%。

表1 澳洲白羊體質(zhì)量與體尺指標(biāo)的表型統(tǒng)計(jì)量

2.2 體質(zhì)量與體尺指標(biāo)的表型相關(guān)性分析

由表2可見，澳洲白羊體質(zhì)量(y)與胸圍(x3)和管圍(x6)相關(guān)性極顯著(P<0.01)，與體高(x1)相關(guān)性顯著(0.010.05)。體質(zhì)量(y)與管圍(x6)的表型相關(guān)系數(shù)最大，為0.680；與胸圍(x3)、體高(x1)、體斜長(x2)、胸深(x4)、胸寬(x5)的表型相關(guān)系數(shù)為0.501，0.393，0.314，0.182，0.154。其中體斜長(x2)與管圍(x6)、胸圍(x3)與胸寬(x5)相關(guān)性極顯著(P<0.01)，其表型相關(guān)系數(shù)分別為0.493，0.450。

表2 澳洲白羊體質(zhì)量與體尺指標(biāo)的表型相關(guān)系數(shù)

2.3 體質(zhì)量與體尺指標(biāo)之間的通徑分析

為了進(jìn)一步揭示體質(zhì)量與體尺指標(biāo)間的相關(guān)關(guān)系，對體質(zhì)量、體尺指標(biāo)進(jìn)行了通徑分析，用體尺指標(biāo)與體質(zhì)量的表型相關(guān)系數(shù)建立以通徑系數(shù)為變量的多元一次方程組，方程組為：

P1+0.291P2+0.272P3+0.079P4+0.042P5-0.050P6=0.393

0.291P1+P2+0.177P3+0.006P4+0.208P5+0.493P6=0.314

0.272P1+0.177P2+P3+0.240P4+0.450P5+0.286P6=0.501

0.079P1+0.006P2+0.240P3+P4+0.110P5+0.313P6=0.182

-0.050P1+0.208P2+0.450P3+0.110P4+P5+0.740P6=0.154

0.345P1+0.493P2+0.286P3+0.313P4+0.740P5+P6=0.680

用矩陣法對多元一次方程組進(jìn)行求解，得到澳洲白羊體尺指標(biāo)對體質(zhì)量的通徑系數(shù)分別為：P1=0.121；P2=-0.081；P3=0.334；P4=-0.099；P5=-0.008；P6=0.614。

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

根據(jù)相關(guān)系數(shù)的組成效應(yīng)，可將澳洲白羊各體尺指標(biāo)(xi)與體質(zhì)量(y)的表型相關(guān)系數(shù)(皮爾森相關(guān)系數(shù))分解為體尺對體質(zhì)量的直接作用和該體尺指標(biāo)通過其他體尺性狀對體質(zhì)量的間接作用兩部分。澳洲白羊體尺對體質(zhì)量相關(guān)系數(shù)的分解見表3。

由表3可知，澳洲白羊的管圍(x6)對體質(zhì)量的直接作用最大，為0.614；胸深(x4)對體質(zhì)量的直接作用最小，為-0.099；體斜長(x2)對體重(y)的綜合間接作用最大，為0.394；其次為體高(x1)、胸深(x4)、胸圍(x3)、胸寬(x5)，分別為：0.271，0.209，0.167，0.161；管圍(x6)的綜合間接作用最小，為0.066。各體尺指標(biāo)對體質(zhì)量的直接作用大小為：管圍(x6)>胸圍(x3)>體高(x1)>胸寬(x5)>體斜長(x2)>胸深(x4)；各體尺指標(biāo)對體質(zhì)量的間接作用大小為：體斜長(x2)>體高(x1)>胸深(x4)>胸圍(x3)>胸寬(x5)>管圍(x6)?？梢钥闯?，體高(x1)對體質(zhì)量既有較強(qiáng)的直接作用，又有較強(qiáng)的間接作用；管圍(x6)、胸圍(x3)主要是通過直接作用對體質(zhì)量產(chǎn)生影響，對體質(zhì)量的間接作用較弱；體斜長(x2)、胸深(x4)主要是通過間接作用對體質(zhì)量產(chǎn)生影響，對體質(zhì)量的直接作用較弱；胸寬(x5)對體質(zhì)量的直接作用、間接作用均較弱。體高(x1)、體斜長(x2)、胸圍(x3)胸深(x4)主要是通過管圍(x6)對體質(zhì)量產(chǎn)生間接作用，胸寬(x5)和管圍(x6)主要是通過胸圍(x3)對體質(zhì)量產(chǎn)生間接作用。

表3 澳洲白羊體尺對體質(zhì)量相關(guān)系數(shù)的分解

2.5 體尺指標(biāo)對體質(zhì)量的逐步回歸分析

2.5.1 體質(zhì)量和各體尺指標(biāo)整體回歸分析模型 從表4～6中可以看出，體高(x1)、體斜長(x2)、胸深(x4)、胸寬(x5)、對體質(zhì)量(y)的影響不顯著，胸圍(x3)、管圍(x6)對體質(zhì)量(y)的影響顯著，因此應(yīng)該剔除自變量體高(x1)、體斜長(x2)、胸深(x4)、胸寬(x5)后，對因變量體質(zhì)量(y)與胸圍(x3)、管圍(x6)進(jìn)行多元逐步回歸分析。

表4 體質(zhì)量和各體尺指標(biāo)整體回歸分析模型

表5 體質(zhì)量和各體尺指標(biāo)整體回歸分析的多元回歸方差分析

表6 體質(zhì)量和各體尺指標(biāo)總體回歸分析的多元回歸系數(shù)參數(shù)估計(jì)值

2.5.2 體質(zhì)量和各體尺指標(biāo)逐步回歸分析模型 逐步回歸分析即按自變量對因變量作用的顯著程度，從大到小逐個(gè)引入回歸方程，即將自變量的t值從大到小逐個(gè)引入回歸方程。由表6可知，各個(gè)自變量的顯著性由大到小依次為管圍(x6)>胸圍(x3)>體高(x1)(不顯著)>胸深(x4)(不顯著)>體斜長(x2)(不顯著)>胸寬(x5)(不顯著)，因此建立了體質(zhì)量(y)和管圍(x6)、胸圍(x3)之間的回歸分析模型2，結(jié)果見表7。由模型2所得到的方程為體質(zhì)量(y)和體尺指標(biāo)之間的最優(yōu)回歸方程為y=-39.698+0.355x3+5.064x6(R=0.752，P<0.05，表7)，標(biāo)準(zhǔn)化后的回歸方程為y=0.334x3+0.585x6(表9)，方差分析的F值為18.815(表8)，達(dá)到極顯著水平(P<0.01)。

表7 澳洲白羊體質(zhì)量和胸圍、管圍指標(biāo)逐步回歸分析模型

表8 體質(zhì)量和胸圍、管圍指標(biāo)逐步回歸分析的多元回歸方差分析

表9 體質(zhì)量和胸圍、管圍指標(biāo)逐步回歸分析的多元回歸系數(shù)參數(shù)估計(jì)值

3 討 論

3.1 體質(zhì)量和體尺指標(biāo)的表型統(tǒng)計(jì)量

變異系數(shù)又稱離散系數(shù)，它反映的是相關(guān)研究變量變化的相對程度指標(biāo)。澳洲白羊的體質(zhì)量與其體尺指標(biāo)的變異系數(shù)由大到小依次為：體質(zhì)量(y)>胸寬(x5)>胸深(x4)>管圍(x6)>胸圍(x3)>體斜長(x2)>體高(x1)，體質(zhì)量的變異系數(shù)為10%，說明澳洲白羊的體質(zhì)量變異幅度較大，體質(zhì)量具有較大的選擇與提高潛力。體高的變異系數(shù)最小，為4.2%，說明體高的選擇與提高潛力很小。

3.2 體質(zhì)量與體尺指標(biāo)的相關(guān)性分析

相關(guān)系數(shù)是用來衡量2個(gè)自變量或1個(gè)自變量與1個(gè)因變量之間相關(guān)程度的指標(biāo)，由于在數(shù)據(jù)分析過程中筆者沒有消除其他研究變量對這兩個(gè)變量的影響，所以表型相關(guān)系數(shù)并不能夠真實(shí)地反映2個(gè)自變量或自變量與因變量之間的相關(guān)關(guān)系，只能反映所研究的2個(gè)自變量或自變量與因變量之間的綜合作用程度。本文通過對澳洲白羊體質(zhì)量與其體尺指標(biāo)進(jìn)行相關(guān)性分析，可以看出澳洲白羊體質(zhì)量(y)與胸圍(x3)和管圍(x6)相關(guān)性極顯著(P<0.01)，與體高(x1)相關(guān)性顯著(0.01

3.3 體質(zhì)量與其體尺指標(biāo)間的通徑分析及其相關(guān)系數(shù)的分解

通徑系數(shù)表示自變量對因變量的直接作用大小，其真實(shí)地反映了2個(gè)變量間的直接相關(guān)程度，它是標(biāo)準(zhǔn)化的偏回歸系數(shù)[18-19]。通過對體質(zhì)量和各體尺指標(biāo)之間的通徑分析以及其相關(guān)系數(shù)的分解可以看出：6月齡澳洲白羊體質(zhì)量(y)與胸圍(x3)和管圍(x6)均表現(xiàn)出極顯著相關(guān)性(P<0.01)。從表3中相關(guān)系數(shù)的分解中可以看出，體高(x1)對體質(zhì)量的直接作用和間接作用均較強(qiáng)；管圍(x6)、胸圍(x3)主要是通過直接作用對體質(zhì)量產(chǎn)生影響，對體質(zhì)量的間接作用較弱；體斜長(x2)、胸深(x4)主要是通過間接作用對體質(zhì)量產(chǎn)生影響，對體質(zhì)量的直接作用較弱；胸寬(x5)對體質(zhì)量的直接作用和間接作用均較弱。通過上述分析說明澳洲白羊的生長發(fā)育是一個(gè)復(fù)雜的有機(jī)整體，各指標(biāo)間有著種種密切而復(fù)雜的聯(lián)系。影響其體質(zhì)量的因素不但多，而且影響途徑亦不同。所以在選種過程中除了加強(qiáng)對體質(zhì)量影響直接作用強(qiáng)的因素如澳洲白羊的管圍、胸圍選擇的同時(shí)，也應(yīng)該關(guān)注其他因素對體質(zhì)量的作用。在選種時(shí)，可以參考本實(shí)驗(yàn)數(shù)據(jù)分析結(jié)果，即本實(shí)驗(yàn)所建立的澳洲白羊體質(zhì)量與其體尺指標(biāo)之間的最優(yōu)回歸方程，進(jìn)行綜合全面考慮，避免盲目選擇，從而能達(dá)到最優(yōu)的選種效果。

3.4 體尺指標(biāo)對體質(zhì)量的逐步回歸分析

通過對澳洲白羊體質(zhì)量與其體尺指標(biāo)之間的逐步回歸分析，得到它們之間的最優(yōu)回歸方程：y=-39.698+5.064x6+0.355x3(R=0.752，P<0.05)，標(biāo)準(zhǔn)化后的回歸方程：y=0.334x3+0.585x6，方差分析的F值為18.815，達(dá)極顯著水平(P<0.01)。體高(x1)、體斜長(x2)、胸深(x4)和胸寬(x5)回歸系數(shù)顯著程度不高(P值均大于0.05)，所以予以排除，說明這4個(gè)體尺指標(biāo)對澳洲白的體質(zhì)量影響不明顯。這與楊燕等[15]對沂蒙黑山羊、李勝開等[16]對隆林山羊以及王欣榮等[17]對甘南草地型藏羊的研究結(jié)果類似。而他們的回歸方程中除了管圍、胸圍入選外，也入選了體高、體斜長、胸深等因素，這主要是由于羊的品種、飼養(yǎng)管理水平、生活環(huán)境等諸多因素不同而造成的。從方程可知，當(dāng)胸圍保持不變時(shí)，管圍每增加1 cm，澳洲白羊的體質(zhì)量增加0.585 kg；當(dāng)管圍保持不變時(shí)，胸圍每增加1 cm，杜泊羊的體質(zhì)量增加0.334 kg?；貧w方程中R2=0.565，說明該回歸方程的回歸平方和占總變異的56.5%，回歸方程顯著存在，所選體尺指標(biāo)對體質(zhì)量有較大的決定作用，方程擬合度較高，具有一定的參考價(jià)值。剩余因子e=1-R2=0.435，說明對澳洲白羊體質(zhì)量產(chǎn)生影響的自變量不止管圍、胸圍，還有一些其他影響因素沒有考慮到，因此下一步需要對影響澳洲白羊體質(zhì)量的因素(包括澳洲白羊的體尺指標(biāo)、營養(yǎng)水平、生存環(huán)境)進(jìn)行全面分析。

4 結(jié) 論

在本試驗(yàn)條件下，6月齡澳洲白羊管圍(x6)和胸圍(x3)是影響其體質(zhì)量的主要因素，體高(x1)、體斜長(x2)、胸圍(x3)、胸深(x4)、主要是通過管圍(x6)對體質(zhì)量產(chǎn)生間接作用，胸寬(x5)和管圍(x6)，主要是通過胸圍(x3)對體質(zhì)量產(chǎn)生間接作用。其體質(zhì)量和體尺指標(biāo)之間的最優(yōu)回歸方程y=-39.698+5.064x6+0.355x3(R=0.752，P<0.05)，標(biāo)準(zhǔn)化后的回歸方程：y=0.334x3+0.585x6(F=18.815，P<0.01)。