蜀興1 號肉兔生長性狀校正公式的制定

李叢艷，謝曉紅，雷 岷，梅秀麗，鄭 潔，郭志強，李鈺瑩，楊 銳，鄺良德

（四川省畜牧科學研究院，動物遺傳育種四川省重點實驗室，四川成都 610066）

生長性狀是畜禽最重要的經濟性狀之一，它在育種方案制定、遺傳評估、選種選配等多個重要環節中均起到重要作用。根據農業農村部《畜禽新品種配套系審定和畜禽遺傳資源鑒定技術規范試行 修訂稿》的規定，要求測定肉兔固定日齡（84 日齡）的體重，而在實際生產中，人們對肉兔的飼養時間關注度較低，出欄體重才是衡量其經濟價值的直接標準。因此，將肉兔生長性狀從測定固定日齡體重轉變為測定達固定體重日齡更具現實意義。生長性狀校正公式對于動物遺傳評估具有重要意義，特別是在生豬上，已經有了深入地研究和較為成熟的理論體系。目前，達100 kg 體重日齡的通用校正公式廣泛應用于大白、長白和杜洛克豬的性能測定[1]，瘦肉型豬[2]或地方豬種[3]達50 kg 體重日齡校正公式也有報道。然而在家兔上，類似的生長性狀校正公式尚未見報道。西南地區（特別是四川）是我國肉兔消費第一大區，消費量占全國70% 以上，人們普遍喜好出欄體重2 000 g 左右的鮮活兔。蜀興1 號是四川省畜牧科學研究院針對這一獨特的消費習慣育成的肉兔配套系，該配套系2 000 g 體重上市，屠宰率高、耗料少。針對該配套系制定生長性狀校正公式，能提高目標性狀選擇的準確性，指導生產。本研究以蜀興1 號肉兔為研究對象，探索制定肉兔達2 000 g 體重日齡的生長性狀校正公式，為蜀興1 號肉兔商品兔的適宜出欄時間提供依據，也為今后其他兔種目標性狀校正公式的制定提供借鑒。

1 材料與方法

1.1 實驗材料 實驗選擇四川省畜牧科學研究院種兔場4 周齡斷奶蜀興1 號商品代肉兔352 只（公母各半），測定其28～84 日齡的個體重，每周稱測1 次。所有實驗兔由同一人員飼養，2 只1 籠（1 公1 母），人工投料，自動飲水，42 日齡前適當限飼，此后自由采食，于上午07:00—08:00 和下午17:00—18:00 按時飼喂。飼料營養水平為消化能10.50 MJ/kg（理論值）、粗蛋白15.82%、粗纖維15.39%、鈣1.02%、總磷0.71%（實測值）。

1.2 研究方法

1.2.1 體重的性別差異研究 根據測定結果，利用ANOVA 程序進行性別間差異顯著性分析，以確定是否分性別進行生長曲線擬合。

1.2.2 動物生長模型的最優選擇 選擇Logistic、Gompertz 以及Von Bertallanffy 模型，利用數學軟件SPSS 17.0 進行擬合，采用F 檢驗法進行顯著性檢驗，通過比較RSS、AdjustedR2、RMSE 可以判定擬合結果的優劣。RSS 越小，RMSE 也越小，AdjustedR2則會越大，表示擬合結果更好。

表1 待定的生長模型

1.2.3 生長性狀校正公式的制定 根據生長發育測定結果對生長模型進行最優選擇，再進行校正公式的制定。其中，校正公式：目標性狀校正值=目標性狀實測值*C；C 為校正系數，C=性狀的目標預測值/ 性狀的實際預測值。

2 結果

2.1 體重的性別差異 方差分析結果表明（表2），不同日齡的蜀興1 號肉兔商品兔公母間生長發育沒有顯著差異，因此，在對生長曲線擬合時可以忽略性別效應，直接進行群體擬合。

表2 蜀興1 號肉兔不同性別間體重差異

2.2 動物生長模型擬合 由表3 可知，3 個模型擬合結果均較好，其AdjustedR2均在0.98 以上，三者中Gompertz 模型的RSS 和RMSE最小，AdjustedR2最大，擬合效果最優，表明蜀興1 號肉兔生長曲線的最佳擬合模型為Gompertz 模型（圖1），方程表達式為其拐點時間為52.65 d，拐點體重為1 525.20 g，最大日增重為41.18 g。因此，進一步選擇Gompertz 模型進行生長性狀校正公式的制定。

表3 動物生長模型方程擬合結果

圖1 蜀興1 號肉兔生長曲線

2.3 生長性狀校正公式制定 根據目標市場需求，將2 000 g 體重作為校正公式的基準。由于目標性狀校正值=目標性狀實測值*C，則t2000=t實*C，校正系數。由于則t=所以，代入數值整理得校正公式：

其中：t校正為校正為2 000 g 體重時的日齡，t實為實際測量時的日齡，y實為實際測量時的體重，單位為g。上述推導過程發現，如果實測體重大于等于a 值或小于等于a ?eb時，則校正公式會無意義，因此校正公式的限制條件為：65.82 g＜y實＜4 145.925 g。

根據校正公式，將各實測日齡體重校正到達2 000 g體重日齡（表4），結果表明：除28 日齡外，其余測定階段的校正結果均在64～65 d 之間，校正結果較為一致。由于63 日齡體重離目標體重最接近，校正結果與實際結果最接近為64.52±3.94 d，因此將其他測定階段校正結果與63 日齡校正結果進行相關分析（表4），結果發現相關性均達極顯著水平，相關系數在63 日齡前隨測定日齡的增加而增加，在63 日齡后隨著測定日齡的增加而降低，其中56～77 日齡（體重范圍為1 674.33～2 464.67g）的相關系數均大于0.96，校正效果較好。以上結果表明：蜀興1 號肉兔適宜的出欄時間為65 d。

表4 不同測定階段校正結果及與63 日齡校正結果的相關系數

3 討 論

3.1 生長模型的選擇 動物生長發育的基本特征具有相對穩定性，可以用數學模型來描述其過程，其實質是將一系列日齡-體重數據壓縮為數學模型中的幾個參數，有效消除部分試驗誤差的影響。本研究選用的Logistic、Gompertz 和Von Bertallanffy 非線性模型被廣泛應用于豬[4]、雞[5]、牛[6]、羊[7]、兔[8-9]、驢[10]等畜種的生長曲線擬合上，在動態反應畜禽生長過程和生長規律、監測畜禽健康狀況、加快畜禽品種選育過程中發揮了重要作用。本研究發現，3 個模型的擬合度均在0.98 以上，能較好地模擬蜀興1 號肉兔體重增長的趨勢和特點，綜合擬合度和多種參數可知，用Gompertz 模型來擬合蜀興1 號肉兔的體重生長發育過程效果更優，這一結果與袁經緯[5]、張樹山[11]、張年[12]等對藏雞、巴馬小型豬、湖北黑頭羊的擬合結果一致，表明蜀興1號肉兔早期生長發育快，符合一般畜禽的生長發育規律。

在家兔的相關研究中，最佳生長模型的選擇不完全相同，多個學者在同時采用Logistic、Gompertz 和Von Bertallanffy 模型進行擬合分析時，發現3 種模型的擬合效果均較好，其中楊翠軍等[13]發現限飼條件下0～180 日齡彩色獺兔的最佳擬合模型為Von Bertalanffy模型，丁鵬等[14]發現加利福尼亞兔×新西蘭大白兔雜交后的橫交（F2）群體35～84 日齡的最佳擬合模型為Logistic 模型，牛曉艷等[15]發現海貍色和白色獺兔1～5月齡的最佳擬合曲線均為Gompertz 和Von Bertalanffy模型，二者的擬合度相同，李叢艷等[9]發現不同營養條件下35～84 日齡齊卡新西蘭兔的最佳擬合曲線均為Gompertz 模型。以上結果的差異可能是由于品種、年齡、營養水平等因素造成的，應該根據家兔的品種特點及生長階段等選擇符合其生長發育規律的模型。

3.2 生長性狀校正公式制定 人們在進行目標性狀測定時，往往不能按照設定的理想時間或體重進行準確測量，特別是進行群體測定時，由于不能實時監測兔只的體重，且其體重受采食和排泄狀態等因素影響，因此很難直接判斷其達目標體重的實際日齡，這就需要利用生長性狀校正公式對實際測定值進行校正以達到同一標準，從而得到較為準確的估計值，便于不同個體、群體或不同場間生長速度的客觀比較分析。生長性狀校正公式是體重和日齡的函數，受生長曲線的影響，不同品種、營養、環境等因素不同，生長性狀校正公式也不同。在瘦肉型豬育種上，生長性狀校正公式分別有加拿大、美國、丹麥、瑞典等校正公式[16]，我國瘦肉型豬育種先后采用了加拿大校正公式[17]和美國校正公式[18]，并開始制定符合國內瘦肉型豬和地方豬種生長規律的生長性狀校正公式。陳斌[19]利用Logistc 和Gompertz 方程對大白、長白和杜洛克豬進行了生長曲線擬合，并利用擬合準確度更高的Lotistic 方程制定了后備豬達100 kg 體重日齡的校正系數，檢驗后發現校正的系數準確度較高，效果較好。劉望宏[16]利用Logistic 模型、Gompertz 模型、Quadratic（二次型）模型、Polynomial（多項式）模型、Saturation 模型及線性方程擬合了大白、長白和杜洛克豬的日齡-體重數據，發現3 個品種的最適研究模型為Saturation 方程W=t/（a+bt），并通過推導t=aw/（1-bw）制定了達50 kg 體重日齡的校正公式。楊岸奇[3]利用Lotistic 模型、Gompertz 模型、Saturation 模型以及Brody 模型擬合了沙子嶺豬的生長曲線，發現前3 個模型的擬合結果均較好，并根據這3 個模型制定了沙子嶺豬達50 kg 體重日齡校正公式，其中Saturation 模型的通用性更佳。張哲等[1]、鄭浩等[20]分別基于美國校正公式和美國及加拿大校正公式制定了中國瘦肉型豬達100 kg 體重日齡的校正公式和里岔黑豬、巴里雜交豬的校正系數，發現其校正效果均較好。以上結果表明，以生長曲線為基礎進行特定品種生長性狀校正公式的制定是可行的。

本研究借鑒陳斌[19]、劉望宏[16]、楊岸奇等[3]的研究思路，采用畜禽常用的Logistc、Gompertz 和Von Bertallanffy 模型擬合了28～84 日齡生長曲線，以擬合效果最好的Gompertz 模型為基礎，以市場需求的2 000 g體重為目標性狀，制定了蜀興1 號肉兔的生長性狀較正公式。由于蜀興1 號肉兔84 日齡前的性別效應不顯著，因此，該校正公式不分性別。終測體重與目標體重越接近，校正日齡的誤差就越小。本研究中測定的63 日齡體重1 947.94 g 與目標體重2 000 g 最接近，通過該測定結果校正的蜀興1 號肉兔達2 000 g 的平均日齡為64.52 d。在實際生產中，終測體重與目標體重越接近，操作起來就越困難，在選擇終測體重范圍時，既要考慮校正結果的可靠性，又要兼顧可操作性。本研究通過對各測定階段的校正結果與63 日齡體重校正結果進行相關分析，發現其相關性均達極顯著水平，且其相關系數大小在63 日齡前隨測定體重的增加而增加，在63 日齡后隨著測定體重的增加而降低，表明測定體重越接近2 000 g，校正效果越好，同時發現，終測體重高于2 000 g時其與63 日齡體重校正結果的相關系數較大，而低于2 000 g 時相關系數較小，這可能是由于前者是校正已經達到的體重，是既定事實，而后者校正的是未來的體重，具有不確定性的原因。其中56～77 日齡體重（即體重范圍為1 674.33～2 464.67 g）校正系數與63 日齡體重校正系數的相關系數均大于0.96，因此，實測體重1 674.33～2 464.67 g 時校正達2 000 g 體重日齡較為適宜。生長性狀校正公式受許多因素的影響，如品種、性別、營養水平、增重等，因此，本研究制定的生長性狀校正公式可以作為參考，下一步可進一步擴大測定群體，以獲得更準確的校正公式。

4 結 論

本研究利用動物生長模型擬合了蜀興1 號肉兔的生長曲線，并以擬合效果最佳的Gompertz 模型為基礎制定了蜀興1 號肉兔達2 000 g 體重日齡的生長性狀校正公式，估計得到其達2 000 g 體重日齡為64.52 d，適宜的實測體重范圍為1 674.33～2 464.67 g。結果表明：蜀興1 號肉兔適宜的出欄時間為65 d。