5種模型估算杜洛克豬的剩余采食量及與飼料效率等指標(biāo)的相關(guān)分析

李 亮，樊福好，曹長仁，謝水華，劉建營，郭建超

（廣東省畜牧技術(shù)推廣總站，廣東 廣州 510520）

隨著飼料原料價格波動和環(huán)保日漸被重視，為了更系統(tǒng)化地提高飼料效率，研究人員從遺傳、配合飼料、飼料添加劑、生產(chǎn)管理等多個方面進(jìn)行了探索改善動物的飼料效率并取得大量成就，為養(yǎng)殖產(chǎn)業(yè)抵抗風(fēng)險、穩(wěn)定發(fā)展提供了強(qiáng)大的動力[1]。剩余采食量（Residual feed intake,RFI）在1963年即被定義為維持體代謝過程和生產(chǎn)所需期望飼料采食量與實際采食量的差值，用于估測動物的飼料效率[2]。且有研究證明RFI是獨(dú)立于采食量、平均日增重等表型的獨(dú)立性狀[3]，其變異可以遺傳并呈現(xiàn)中到高的遺傳力[4]；所以利用RFI選種可以降低飼料采食量而不降低增重，是提高動物種群經(jīng)濟(jì)效益的有效選擇性狀，被認(rèn)為可以作為一種更加準(zhǔn)確的評價飼料效率的指標(biāo)，在選種上具備更大的優(yōu)勢，可以有更多應(yīng)用[5]。

前期對農(nóng)業(yè)農(nóng)村部種豬質(zhì)量監(jiān)督檢驗測試中心（廣州）種豬測定數(shù)據(jù)統(tǒng)計發(fā)現(xiàn)不同品種間豬的各項測定指標(biāo)會有差異[6]，各個品種群內(nèi)RFI與各項指標(biāo)的關(guān)系也呈現(xiàn)出了一定的關(guān)系[7]。本研究以2017—2019年集中測定的杜洛克種公豬測定數(shù)據(jù)為基礎(chǔ)，使用統(tǒng)計軟件SPSS，分析不同RFI評估模型以及RFI與背膘、眼肌面積、日增重等標(biāo)準(zhǔn)測定指標(biāo)之間的關(guān)系，以期為后期研究提供更多的基礎(chǔ)數(shù)據(jù)和研究思路，助力育種和生產(chǎn)的發(fā)展。

1 材料與方法

1.1 數(shù)據(jù)采集

在2017—2019年農(nóng)業(yè)農(nóng)村部種豬質(zhì)量監(jiān)督檢驗測試中心（廣州）所進(jìn)行的種豬測定所得數(shù)據(jù)中，選擇杜洛克種公豬共359頭。受測種公豬于體重25 kg左右開始用美國奧斯本和廣東廣興自動喂料系統(tǒng)測定其肥育期的采食量（Feed Intake,FI）、平均日采食量（Average feed intake,ADFI）、增重（Weight gain,WG）、平均日增重（Average daily weight gain,ADG）、料重比（Feed conversion ratio,FCR）、眼肌面積（loin muscle area,LMA）等指標(biāo)，使用荷蘭Pie Medical公司Aquila Vet型B超儀進(jìn)行背膘厚、眼肌面積的測量，根據(jù)《種豬生產(chǎn)性能測定規(guī)程》要求，計算校正至100 kg體重日齡、校正至100 kg體重背膘厚等性狀。

1.2 估算方法

RFI的估算采取表1中的5種估算方法[8]。表1式中，*表示乘號；ADFI為平均日采食量；OnBW為測定初期體重；OffBW為測定末期體重；MWT為測定中期體重；ADG為平均日增重；BF為平均背膘厚；LMA為眼肌面積；b1～b5為對應(yīng)性狀的回歸系數(shù)。

表1 不同模型的RFI估算公式

在個體生產(chǎn)性能記錄的基礎(chǔ)上，利用SPSS建立多重線性回歸模型，得到全群RFI的線性表型回歸系數(shù)，代入公式算出每個個體的RFI，再通過SPSS相關(guān)功能分析不同模型RFI與測定指標(biāo)的相關(guān)性。將RFI升序排列并等分為高、中、低3組，分析測定指標(biāo)的組間差異及分組的影響。選擇SPSS數(shù)據(jù)分析軟件的單因素ANOVA或非參數(shù)檢驗進(jìn)行數(shù)據(jù)分析，應(yīng)用事后LSD或Tamhane進(jìn)行多重比較，結(jié)果以平均數(shù)±標(biāo)準(zhǔn)差表示。使用Pearson或Spearman分析各指標(biāo)的相關(guān)性。

2 分析與討論

2.1 不同估算模型結(jié)果與測定指標(biāo)相關(guān)性

分析不同模型的估算RFI與測定指標(biāo)相關(guān)性，結(jié)果如表2。

表2 不同模型的估算RFI與測定指標(biāo)相關(guān)性

各模型估算的RFI與杜洛克豬群的測定期料重比都有極顯著正相關(guān)（P＜0.01），說明這幾個模型估算的RFI都可以很好地反映杜洛克種豬的飼料效率的差異；RFIa、RFIb、RFId、RFIe 與校正至 100 kg體重背膘厚呈顯著負(fù)相關(guān)（P＜0.05），說明大部分模型的估測都與瘦肉率呈現(xiàn)正相關(guān)，但RFIc與背膘厚呈極顯著正相關(guān)（P＜0.01），有點(diǎn)特別。此外，各RFI又都表現(xiàn)出與眼肌面積負(fù)相關(guān)，但基本都不顯著，也許低RFI的杜洛克豬雖然背膘更薄，同時眼肌面積可能也會更小。但總體而言，低RFI的杜洛克豬瘦肉率可能會更高。也許RFI與背膘厚的負(fù)相關(guān)性與其估算公式中背膘厚的系數(shù)是負(fù)數(shù)有關(guān)，但ADG的系數(shù)也是負(fù)數(shù)，而各模型估算的RFI與ADG均呈現(xiàn)正相關(guān)，所以相關(guān)系數(shù)與估算公式的系數(shù)關(guān)系不大，更重要的是與背膘厚、眼肌面積數(shù)據(jù)的變化趨勢有關(guān)。

校正至100 kg體重日齡與各RFI都呈現(xiàn)負(fù)相關(guān)，且與 RFIc、RFId、RFIe顯著負(fù)相關(guān)（P＜0.05），可能低RFI的校正至100 kg體重日齡會較大，可能比高RFI組長得稍慢，但其與RFIa、RFIb兩個模型的結(jié)果差異不顯著（P＞0.05）；同時5組RFI與校正30～100 kg體重日增重呈正相關(guān)（P＞0.05），表示低RFI組此階段日增重會更小，但同時其校正至100 kg體重日齡又較大。分析其原因，有可能是低RFI的杜洛克豬雖然料重比低、瘦肉率高，但可能日增重較小、生長相對較慢。與長白、大白的相關(guān)數(shù)據(jù)相比，杜洛克的生長速度、增重都屬中等偏低[6]，可能RFI這種特點(diǎn)亦與其品種特性有關(guān)。此時RFI的這種規(guī)律有可能要求育種時要評估飼料效率與生長速度之間的取舍與平衡。

2.2 RFI低、中、高組測定指標(biāo)的差異

對各模型估算的RFI分為低、中、高3組，測定指標(biāo)組間差異對比見表3。

表3 各模型低、中、高RFI組間測定指標(biāo)差異

各個模型估算的RFI不同組間的測定期料重比都有極顯著差異（P＜0.01），平均值從低到高的排序均為：低RFI組＜中RFI組＜高RFI組，與表2傳遞的信息相互呼應(yīng)，很有力地再次揭示了杜洛克豬RFI與飼料效率的密切相關(guān)。

對于校正至100 kg體重日齡和校正30～100 kg體重日增重而言，a、b、e 3個模型的RFI分組都沒有表現(xiàn)出明顯的變化梯度，但c、d兩個模型一致體現(xiàn)出顯著或極顯著的梯度差異，這一點(diǎn)與表1的結(jié)果相呼應(yīng)，可能RFI低的杜洛克豬在100 kg之前日增重確實比較低、生長相對較慢。

對校正至100 kg體重膘厚而言，5組RFI呈現(xiàn)出了兩種不同的趨勢，a、b、d、e 4組低RFI組的數(shù)值均是最大值、中RFI組次之、高RFI組最小，而c組正好相反。c組與其他各組的不同點(diǎn)在于估算模型的自變量不包含BF，至于這是不是它的結(jié)果與眾不同的原因，有待研究。同時，對校正至100 kg體重眼肌面積而言，除e組外，各組均未顯現(xiàn)梯度變化趨勢。結(jié)合日增重指標(biāo)，可以說a、b模型低RFI組瘦肉率顯著低于高RFI組（P＜0.05），d組也有這個趨勢但不顯著（P＞0.05），而 c、e組的結(jié)論剛好相反（P＜0.05），這跟表2顯示的信息基本一致。暫時少數(shù)服從多數(shù)，認(rèn)為低RFI組瘦肉率更高，也許是c組模型本身有問題，有待進(jìn)一步研究。

2.3 各RFI分組與測定指標(biāo)的相關(guān)性

對各模型估算的RFI分為低、中、高3組，測定分組操作與各指標(biāo)的相關(guān)系數(shù)，結(jié)果見表4。

表4 根據(jù)RFI大小分組與各指標(biāo)的相關(guān)系數(shù)

表4顯示，各RFI的分組情況與豬群測定期料重比均有極顯著相關(guān)（P＜0.01），同時與校正至100 kg體重眼肌面積、校正30～100 kg體重日增重之間亦表現(xiàn)出與表2、表3相印證之信息，揭示依大小為RFI分組確有必要，依此作選種依據(jù)應(yīng)符合生產(chǎn)所需。其他指標(biāo)與分組無顯著相關(guān)（P＞0.05），但其均值亦顯示少許趨勢。此外，由于a、b、d模型估算的RFI與瘦肉率都呈現(xiàn)顯著負(fù)相關(guān)（P＜0.05）或類似趨勢，但c、e模型估算的結(jié)果似乎可以得出相反結(jié)論，與以前的研究結(jié)論接近[7]，可能是由于品種的關(guān)系，也可能是由于模型估算公式需要更多研究改善。

3 結(jié)論與展望

5種模型估算的RFI均與飼料效率有極顯著相關(guān)，可有效評估杜洛克種豬飼料效率差異；5種模型的RFI與背膘厚、眼肌面積的關(guān)系表現(xiàn)出兩種相反的趨勢，a、b、d、e 模型中，低 RFI組背膘更厚，眼肌面積趨勢不明顯；c模型中，低RFI組背膘更薄，眼肌面積趨勢不明顯，所以可能低RFI組背膘更厚。同時c、d、e模型低RFI組表現(xiàn)出與生長速度更慢（P＜0.05），c、d模型低RFI組表現(xiàn)出日增重更低的趨勢（P＜0.05）。其呈現(xiàn)的規(guī)律提示可能低RFI豬雖然飼料效率較好，但生長會更慢、瘦肉率會更低，這可能與品種有關(guān)系，但不同的估算公式亦造成的不同結(jié)論，還需要進(jìn)一步試驗研究公式間的異同。