何愛華，李志剛，張字文，余 龍

（云南邦格農(nóng)業(yè)集團(tuán)有限公司，云南 昆明 650200）

滇南小耳豬繁殖性狀的相關(guān)及通徑分析

何愛華，李志剛，張字文，余 龍

（云南邦格農(nóng)業(yè)集團(tuán)有限公司，云南 昆明 650200）

本研究對120頭滇南小耳豬的產(chǎn)總仔數(shù)（X1），產(chǎn)活仔數(shù)（X2），出生窩重（X3），出生均重（X4），35日齡斷奶仔數(shù)（X5），35日齡斷奶個(gè)體重（X6），35日齡斷奶窩重（Y），7個(gè)繁殖性狀進(jìn)行了相關(guān)分析和通徑分析，為其早期選種提供理論依據(jù)。研究結(jié)果表明，各個(gè)繁殖性狀在表型值上都存在不同程度的相關(guān)，其中總產(chǎn)仔數(shù)（X1）與產(chǎn)活仔數(shù)（X2）相關(guān)系數(shù)最大（0.89402），而35日齡斷奶個(gè)體重（X6）和35日齡斷奶仔豬數(shù)（X5），對35日齡斷奶窩重（Y）有極顯著的影響（P＜0.01），在此基礎(chǔ)上，建立估計(jì)35日齡斷奶窩重（Y）的最優(yōu)回歸方程為Y=4.3576-0.59901X5+0.1375X6。

滇南小耳豬；繁殖性狀；通徑分析；相關(guān)系數(shù)；偏相關(guān)系數(shù)；回歸方程

滇南小耳豬，是我國優(yōu)良的地方品種，國家級畜禽遺傳資源保護(hù)品種，是云南省南方地區(qū)開展特色優(yōu)質(zhì)豬肉產(chǎn)業(yè)化生產(chǎn)的寶貴資源[1]。滇南小耳豬肉質(zhì)鮮美，富含多種營養(yǎng)物質(zhì)，但存在產(chǎn)仔數(shù)少、繁殖能力低等缺點(diǎn)，導(dǎo)致養(yǎng)殖規(guī)模小、養(yǎng)殖數(shù)量少，市場供不應(yīng)求。因此重視育種工作中滇南小耳豬的繁殖性狀的選育，一方面可提高養(yǎng)殖戶的生產(chǎn)效益及經(jīng)濟(jì)水平，另一方面可使滇南小耳豬的養(yǎng)殖規(guī)模不斷增大、遺傳資源得以傳承。本文根據(jù)西雙版納邦格牧業(yè)科技有限公司滇南小耳豬核心種豬場120窩母豬的記錄資料，對繁殖性狀的多個(gè)指標(biāo)進(jìn)行了相關(guān)分析，從而找出其相互間的更深層次的關(guān)系，為滇南小耳豬的選種過程進(jìn)行早期估測、為其繁殖性狀的選育工作提供理論基礎(chǔ)。

1 材料與方法

1.1 材料 所用數(shù)據(jù)材料來自云南省西雙版納州邦格牧業(yè)科技有限公司，選用2014年5月和6月的滇南小耳豬各胎次母豬的產(chǎn)仔哺乳記錄，共120窩。

1.2 飼養(yǎng)管理 母豬的飼養(yǎng)管理遵循現(xiàn)代化、規(guī)?；B(yǎng)殖流程。妊娠期母豬采用半限位欄飼養(yǎng)，以便有足夠的活動(dòng)空間，采用自動(dòng)投喂器同步落料，每天飼喂2次，早上9:00和下午5:00；飼喂量：妊娠前、中期0.8～1.0 kg/d，妊娠后期1.0～1.5 kg/ d，下午添加一次青綠飼料1 kg/（頭·d）（新鮮香蕉枝干、紅薯藤等）。臨產(chǎn)前7 d清洗母豬后進(jìn)入產(chǎn)床單元房（高架床）產(chǎn)仔哺乳，自由采食，保溫設(shè)施齊備，采用全進(jìn)全出模式。

仔豬出生稱重，7日齡誘食，35日齡斷奶；全場實(shí)行統(tǒng)一的飼養(yǎng)管理和免疫接種。

1.3 繁殖性狀 本研究的繁殖性狀包括產(chǎn)總仔數(shù)（X1）、產(chǎn)活仔數(shù)（X2）、出生窩重（X3）、出生個(gè)體重（X4）、斷奶仔數(shù)（X5）、35日齡斷奶個(gè)體重（X6）、斷奶窩重（Y），共7個(gè)性狀。

1.4 統(tǒng)計(jì)分析方法 首先在EXCEL上建立原始數(shù)據(jù)庫，然后用SPSS 19.0軟件進(jìn)行統(tǒng)計(jì)分析。對所有試驗(yàn)數(shù)據(jù)整理計(jì)算出各個(gè)性狀的表型值參數(shù)以及性狀間的簡單相關(guān)系數(shù)、通徑系數(shù)等，并建立最優(yōu)回歸方程。

2 結(jié)果與分析

2.1 繁殖性狀的表型值參數(shù) 滇南小耳豬繁殖性能表型值參數(shù)結(jié)果見表1。

表1 滇南小耳豬母豬繁殖性能的表型參數(shù)Table1 Diannan sma ll-ear pigs on phenotypic param eters for rep roductive trait

由表1可見，產(chǎn)總仔數(shù)（X1）9.14±2.30頭，產(chǎn)活仔數(shù)（X2）8.41±2.03頭，出生個(gè)體重（X4）0.64± 0.08 kg，斷奶仔數(shù)（X5）7.14±1.36頭，35日齡斷奶個(gè)體重（X6）4.59±1.02 kg，可見其整個(gè)哺乳期間日增重低。繁殖性狀的變異系數(shù)較大，產(chǎn)總仔數(shù)（X1）、產(chǎn)活仔數(shù)（X2）、出生窩重（X3）、斷奶個(gè)體重（X6）、斷奶窩重（Y）變異系數(shù)皆在20%以上；其中出生個(gè)體重（X4）的變異系數(shù)最低，斷奶窩重（Y）的變異系數(shù)最大。變異系數(shù)大，選擇差就大，說明通過選育提高產(chǎn)總仔數(shù)、出生窩重、斷奶窩重的潛力較大。2.2 繁殖性狀間的簡單相關(guān)系數(shù) 繁殖性狀間的簡單相關(guān)系數(shù)結(jié)果見表2。

表2 繁殖性狀間的相關(guān)系數(shù)Table2 The Correlation coefficient in reproductive traits

由表2中可以看出，產(chǎn)總仔數(shù)（X1）與產(chǎn)活仔數(shù)（X2）相關(guān)系數(shù)最大（0.89402）；出生窩重（X3）與產(chǎn)總仔數(shù)（X1）、產(chǎn)活仔數(shù)（X2）、出生個(gè)體重（X4）呈正相關(guān)（0.78284、0.87072、0.44718），相關(guān)均極顯著（P＜0.01）。出生個(gè)體重（X4）與產(chǎn)總仔數(shù)（X1）、產(chǎn)活仔數(shù)（X2）、斷奶仔數(shù)（X5）呈負(fù)相關(guān)（-0.03021、-0.03827、-0.00651），相關(guān)性不顯著（P＞0.05）；與出生窩重（X3）呈正相關(guān)（0.44718），相關(guān)極顯著（P＜0.01）；這結(jié)果與呂繹如等[2]、姿義洲等[3]研究滇南小耳豬的結(jié)果一致。

斷奶窩重（Y）與產(chǎn)總仔數(shù)（X1）、產(chǎn)活仔數(shù)（X2）、出生窩重（X3）、斷奶仔數(shù)（X4）、斷奶個(gè)體重（X6）均呈正相關(guān)，相關(guān)均極顯著（P＜0.01）；與斷奶窩重（Y）相關(guān)系數(shù)（0.7635）最大的繁殖性狀為35日齡斷奶個(gè)體重（X6），這與劉慶雨等[4]、陳海燕等[5]的研究相近。

表3 繁殖性能的偏相關(guān)系數(shù)Tab le3 The partial corre lation coefficient in rep roductive traits

2.3 繁殖性狀間的偏相關(guān)系數(shù) 見表3.

從表3中可以看出，出生個(gè)體重（X4）與產(chǎn)活仔數(shù)（X2）的偏相關(guān)系數(shù)最大（0.9986），呈極顯著（P＜0.01），斷奶窩重（Y）和斷奶仔數(shù)（X5）的相關(guān)系數(shù)次之（0.99586），呈極顯著（P＜0.01）。

2.4 繁殖性狀間的通徑分析 根據(jù)通徑分析原理，以產(chǎn)總仔數(shù)（X1）、產(chǎn)活仔數(shù)（X2）、出生窩重（X3）、出生個(gè)體重（X4）、斷奶仔數(shù)（X5）、斷奶個(gè)體重（X6）為自變量，以斷奶窩重（Y）為依變量進(jìn)行通徑分析，列出如下正規(guī)方程組：

解方程組得：P1=0.0386，P2=-0.0027，P3=-0.0277，P4=0.0027，P5=0.6272，P6=0.7610。

相關(guān)分析僅僅反映了變量的表面關(guān)系，通徑分析更準(zhǔn)確的反映了變量的真實(shí)關(guān)系，通過通徑分析更可能的反映出多個(gè)性狀間的相對影響程度，將相關(guān)系數(shù)剖分為自變量對依變量的直接效應(yīng)。從表4中可看出，對斷奶窩重（Y）直接影響最大是35日齡斷奶個(gè)體重（X6），其系數(shù)是0.7610，其次是斷奶仔數(shù)（X5）（0.6272），這與李明等[6]、蘇從文等[7]研究結(jié)果中對斷奶窩重的直接影響最大是斷奶仔數(shù)有所差異。

自變量對依變量除直接影響外，還通過其他自變量對依變量產(chǎn)生間接影響。因此，對相關(guān)系數(shù)進(jìn)行剖析，分析了間接影響，表4中結(jié)果表明，產(chǎn)活仔數(shù)（X2）和出生窩重（X3）等性狀雖然對斷奶窩重（Y）的直接影響最小（-0.0027，-0.0277），但通過斷奶個(gè)體重（X6）的間接影響，最后對斷奶窩重（Y）間接效應(yīng)分別為0.5084和0.5364；而斷奶仔數(shù)（X5）對斷奶窩重（Y）的直接作用和間接作用都很大，分別為0.6272和0.8050。

表4 繁殖性狀的間的通徑分析Table4 The path analysis in rep roductive traits

實(shí)踐中，除了研究原因?qū)Y(jié)果的相對影響程度外，還要研究原因?qū)Y(jié)果的絕對影響程度。絕對影響程度用決定系數(shù)表示。表5中計(jì)算結(jié)果表明，對斷奶窩重（Y）的影響最大的是斷奶個(gè)體重（X6）（0.582926），其次是斷奶仔數(shù)（X5）（0.393365），這與陶勇等[8]研究結(jié)果不一致。總的決定系數(shù)0.9865，誤差系數(shù)0.0135，說明大多數(shù)影響斷奶窩重（Y）的性狀已經(jīng)分析到，而其他因素在本研究中影響較小。

表5 繁殖性狀間的決定系數(shù)剖析Table5 The determ ination coefficient in reproductive traits

2.5 建立最優(yōu)回歸方程 通過逐步回歸分析，逐步剔除了產(chǎn)總仔數(shù)（X1）、產(chǎn)活仔數(shù)（X2）、出生窩重（X3）、出生個(gè)體重（X4）共4個(gè)性狀，剩余2個(gè)性狀斷奶仔數(shù)（X5）、斷奶個(gè)體重（X6）與斷奶窩重（Y）間存在極顯著的線性關(guān)系，相關(guān)指數(shù)R2=0.9673，在此基礎(chǔ)上建立了估計(jì)滇南小耳豬斷奶窩重的最優(yōu)回歸方程為：Y=4.3576-0.59901X5+0.1375X6。若用Y與X5、X6之間的最優(yōu)回歸方程來估測Y，其可靠程度達(dá)96.73%，表明影響斷奶窩重（Y）的主要性狀都已包含在方程中，該方程能夠較真實(shí)地反映這些性狀對斷奶窩重（Y）的影響。

3 討論

通徑系數(shù)能夠有效地表示相關(guān)變量之間原因?qū)Y(jié)果的直接效應(yīng)和間接效應(yīng)，它比相關(guān)和回歸分析更為精確，對多變數(shù)資料的統(tǒng)計(jì)分析更符合實(shí)際，更有利于指導(dǎo)育種工作[9-10]。本研究分析結(jié)果可以看出，滇南小耳豬繁殖性能的變異系數(shù)大，說明選育提高有較大的遺傳潛力。出生個(gè)體重（X4）與產(chǎn)活仔數(shù)（X2）呈負(fù)相關(guān)，說明若過分追求產(chǎn)仔數(shù)的提高就會(huì)導(dǎo)致頭均重的下降，因本研究中滇南小耳豬的出生個(gè)體重0.64 kg且其奶頭數(shù)保持在5對，且很少超過5對，若產(chǎn)仔數(shù)增加、出生個(gè)體重降低，會(huì)影響仔豬的成活率，故須把選育重點(diǎn)放在斷奶窩重（Y）、斷奶仔數(shù)（X5）等性狀上。

對其斷奶窩重（Y）的主要決定因素為斷奶仔數(shù)（X5），其次為斷奶個(gè)體重（X6）。在豬生產(chǎn)學(xué)中，斷奶窩重（Y）體現(xiàn)了母豬的繁殖力[11]，斷奶仔數(shù)（X5）是決定母豬年生產(chǎn)能力的重要指標(biāo)，它影響著年生產(chǎn)量和獲得的經(jīng)濟(jì)效益。因此，加強(qiáng)對斷奶仔數(shù)（X5）和斷奶個(gè)體重（X6）的選育不僅提高母豬的繁殖能力，更在一定程度上提高產(chǎn)量、增加經(jīng)濟(jì)收入。但由于繁殖性狀屬于低遺傳力性狀，受環(huán)境因素的影響比較大。所以，在實(shí)際生產(chǎn)中，除了通過對繁殖性狀選擇提高斷奶窩重（Y）外，重點(diǎn)還要改善飼養(yǎng)管理?xiàng)l件和營養(yǎng)水平，加強(qiáng)仔豬疾病的防治，經(jīng)過一段時(shí)間的選育滇南小耳豬會(huì)在體重、窩重性能上會(huì)有較大的提高，從而為促進(jìn)種豬的生產(chǎn)效益起到積極的作用。

參考文獻(xiàn)

[1]李志剛，何愛華，張字文.滇南小耳豬胴體品質(zhì)及肉質(zhì)的研究[A].第十二屆（2014）中國豬業(yè)發(fā)展大會(huì)論文集[C].2014年12屆中國豬業(yè)發(fā)展大會(huì)：99-102.

[2]呂繹如，黃安樞，王家明.滇南小耳豬繁殖性能初析[J].云南畜牧獸醫(yī)，1982，1：1-12.

[3]姿義洲，張念嚴(yán).滇南小耳豬幾個(gè)性狀遺傳參數(shù)的研究[J].遺傳，1982，4（3）：11-14.

[4]劉慶雨，張樹敏，于永生.松遼黑豬繁殖性狀的相關(guān)及通徑分析[J].吉林畜牧獸醫(yī)，2012，4：28-30.

[5]陳海燕，傅衍，陸建定.利用通徑分析法研究臺系杜洛克的繁殖性狀[J].畜牧生產(chǎn)，2003（4）：14-15.

[6]李明，孫桂榮，劉孟洲，等.丹系長白母豬繁殖性狀通徑分析及最優(yōu)回歸方程的建立[J].中國農(nóng)學(xué)通報(bào)，2006，22（9）：35-38.

[7]蘇從文，宗紹志.大白豬繁殖性狀的相關(guān)及通徑分析[J].畜牧獸醫(yī)雜志，2010，29（3）：14-17.

[8]陶勇，董曉君，張紹祥.蘇太豬繁殖性狀間相關(guān)分析[J].當(dāng)代畜牧，2007，12：24-25.

[9]吳其文，陳建錨，肖天放，等．莆田黑豬繁殖性狀的通徑分析[J].福建畜牧獸醫(yī)，2009，31（2）：4-6.

[10]張毅，向釗，朱丹，等.榮昌豬繁殖性狀的研究[J].中國畜牧雜志，2004，40（5）：52-55.

[11]楊公社.豬生產(chǎn)學(xué)[M].北京：中國農(nóng)業(yè)出版社，2002.

Correlation and Path Analysis for Reproductive Traits ofDiannan Small-ear Pigs

He Aihua,Li Zhigang,Zhang Ziwen,Yu Long
（Yunnan bankgreen group Co.,Ltd.Yunnan Kunming 650200）

The 120 Diannan ear pigs sows were used to carry cor relation and path analysis on reproductive t raits(total number of farrowed Piglets(X1),the number of l iving newborn Piglets(X2), the weight of a nest of newborn Piglets(X3),the average weight of newborn Piglets(X4),the number of weaned Piglets(X5),the average weight of 35-day old weaned piglets(X6),the weight of a nest of 35-day old weaned piglets(Y))in the paper.The resul ts showed that each reproductive t rait had various degrees cor relation with the other,the biggest cor relation coef ficient appeared in the relationship of the total number of farrowed newborn Piglets(X1)and the number of l iving newborn Piglets(X2)(0.89402).The two factors,the number of weaned Piglets(X5)and the average weight of 35-day old Piglets(X6)had remarkable ef fect on the weight of a nest of weaned 35-day old Piglets(Y)(P＜0.01).Based on above,in order to estimate the weight of a nest of weaned 35-day old Piglets(Y),the best regression mode was established as fol lows:Y=4.3576-0.59901X5+ 0.1375X6.

Diannan smal l-ear Pigs;Reproductive trait;Path analysis;Cor relation coef f icient; Partial Cor relation coef ficient;Regression mode

S828.3

1672-9692(2015)07-0021-06

2015-05-14

何愛華（1984-），女，碩士學(xué)位，中級畜牧師，主要從事動(dòng)物遺傳育種與生產(chǎn)研究。

云南省院省校科技合作專項(xiàng)（2013-01-03-B-H-02-0001）