烏蒙鳳雞零世代產蛋曲線擬合研究

王天松，王 梅，張 勇*，郭徵力，莫先艇，駱科印，楊 紅，曾姍姍，左 麗

（1.貴州大學動物科學學院，高原山地動物遺傳育種與繁殖教育部重點實驗室，貴州貴陽 550025；2.六盤水市農業(yè)農村局，貴州六盤水 553001；3.畢節(jié)市畜牧站，貴州畢節(jié) 551700）

烏蒙鳳雞原產于貴州省六盤水市，其中心產區(qū)海拔高度介于1 400～1 900 m，產區(qū)特殊的喀斯特地形地貌減少了雞種間的基因交流，烏蒙鳳雞基因較為純合，以獨特的鳳頭性狀而保存下來。烏蒙鳳雞是以肉用性能為主的地方雞種，具有體型適中，雌雄雞鳳頭明顯、耐粗飼、抗逆性強、覓食能力強及突出的高原特色等特點[1]。目前，烏蒙鳳雞雌雞群體羽色多為麻羽和黃羽，少數(shù)為白羽和黑羽，具有明顯的冠羽特征，初產日齡在181～210 d。烏蒙鳳雞經保種擴繁和推廣，已申報獲得國家農產品地理標志登記保護及農產品地理標志登記證書。

產蛋量是雞最重要的數(shù)量性狀，受品種、年齡、光照、營養(yǎng)水平等諸多因素影響，然而產蛋量變化通常遵循一定規(guī)律，服從一定曲線分布[2]。如果將日產蛋率及累計產蛋數(shù)標記在坐標軸上，將各點連接起來，則可構成產蛋曲線，產蛋曲線數(shù)學模型可從動態(tài)角度展示產蛋率隨時間變化的規(guī)律[3]。雞產蛋曲線不僅反映其產蛋率高低，還可以反映雞不同時期的飼養(yǎng)管理水平，可以指導生產，使雞在標準體重下發(fā)揮最佳產蛋水平[4]。烏蒙鳳雞作為新發(fā)現(xiàn)的遺傳資源，到目前為止，相關研究工作報道較少，因此本實驗開展烏蒙鳳雞的產蛋曲線擬合研究，探索烏蒙鳳雞產蛋規(guī)律，為進一步加大其選育力度、提高產蛋率、促進其保種擴繁奠定基礎。

1 材料與方法

1.1 實驗材料 實驗在六盤水金鳳棲自然牧業(yè)有限公司烏蒙鳳雞保種擴繁場進行。選取零世代烏蒙鳳雞母雞360 羽，相同條件下單籠飼養(yǎng)，母雞開產后分別記錄28～40 周齡產蛋數(shù)。

1.2 實驗方法

1.2.1 實驗期日糧營養(yǎng)水平 實驗雞不同生長階段的主要營養(yǎng)需要見表1，根據產蛋期采食量、產蛋率、蛋重、蛋重周增長及體重變化等調整飼料配方。實驗期間，定期觀察雞群狀況，做好死亡數(shù)、產蛋率、存欄數(shù)和產蛋量等相關記錄。

表1 烏蒙鳳雞營養(yǎng)需要量

1.2.2 產蛋性能統(tǒng)計 記錄實驗雞的開產日齡（日產蛋率達5%時）、產蛋量，匯總周產蛋量和存欄雞數(shù)，統(tǒng)計各周齡產蛋率及累計產蛋數(shù)，計算公式如下：周產蛋率（%）=各周產蛋量/統(tǒng)計期內總飼養(yǎng)母禽只數(shù)×100%；周累計產蛋數(shù)（個）=周累計產蛋個數(shù)/存欄母雞數(shù)。

1.3 數(shù)據處理 采用Excel 2010 進行數(shù)據初處理，用SPSS 22.0 軟件計算產蛋率和累計產蛋數(shù)模型擬合參數(shù)值。伍德模型、楊寧模型、分室模型[5-7]對產蛋率曲線擬合公式見表2。Logistic、Gompertz 和Von Bertalanffy模型[8-10]對累計產蛋數(shù)曲線擬合公式見表3。

表2 產蛋率擬合模型

表3 累計產蛋數(shù)擬合模型

2 結果與分析

2.1 3 種模型產蛋率曲線擬合結果 3 種數(shù)學模型擬合得到的各模型周產蛋率參數(shù)估計值見表4。把參數(shù)估計值a、b、c、d 分別代入3 個模型公式得到：伍德模型y（t）=（8.758E-44）×t40.489×exp（-1.159t）；楊寧模型y（t）=221.095×exp（-0.035t）/[1+exp（-1.665×（t-29.569））]；分室模型y（t）=392.429×[1-exp（-0.423（t-27.898））]×exp（-0.049t）。由表4 可知，楊寧模型擬合度R2最高，為0.940，伍德模型和分室模型R2分別為0.589、0.849。

表4 烏蒙鳳雞3 種產蛋率擬合模型參數(shù)值

2.2 產蛋率觀測值與擬合值比較 各周齡產蛋率觀測值和模型擬合值比較結果見表5。烏蒙鳳雞28 周齡產蛋率為10%，28～32 周齡產蛋率持續(xù)上升，32 周產蛋率達76.11%。產蛋高峰期為31～40 周齡，產蛋率在52%以上。從3 個模型擬合值來看，28 周產蛋率伍德模型擬合值為27.89%，大于實際觀測值，35 周齡達到產蛋高峰，產蛋率為70.24%，晚于實際產蛋高峰，與觀測值相差5.87%。28 周產蛋率楊寧模型擬合值為5.70%，小于實際觀測值，但與觀測值最相近，32 周齡達到最高產蛋率71.29%，與觀測值相差4.82%。28 周產蛋率分室模型擬合值為4.18%，產蛋高峰在33 周齡，產蛋率為68.29%，與觀測值相差7.82%。綜上，楊寧模型更適合烏蒙鳳雞產蛋率曲線的擬合。

表5 烏蒙鳳雞周產蛋率實際觀測值與3 種模型擬合值

2.3 3 種模型累計產蛋數(shù)曲線擬合結果 烏蒙鳳雞零世代累計產蛋數(shù)3 種非線性模型擬合參數(shù)估計值和擬合度見表6。把參數(shù)估計值A、B、k 分別代入3 個模型公式得到：Logistic：y（t）=51.808/[1+22.944×exp（-0.433×t）]；Gompertz：y（t）=60.316×exp[-4.508×exp（-0.234×t）]；Bertalanffy：y（t）=68.354×[1-0.899×exp（-0.166×t）]3。由表6 可知，Logistic、Gompertz 和Von Bertalanffy 模型對烏蒙鳳雞累計產蛋數(shù)擬合度分別為0.989、0.996、0.998，以Von Bertalanffy 模型擬合度最高，用其計算出的拐點產蛋周齡為33.950 周，拐點產蛋數(shù)為20.253個，極限產蛋數(shù)為68.354 個。Logistic 模型拐點產蛋周齡為35.236 周，拐點產蛋數(shù)為25.904 個，極限產蛋數(shù)為51.808 個。Gompertz 模型的拐點產蛋周齡為34.345周，拐點產蛋數(shù)為22.191 個，極限產蛋數(shù)為60.316 個。

表6 烏蒙鳳雞3 種累計產蛋數(shù)擬合模型參數(shù)值

2.4 累計產蛋數(shù)觀測值與擬合值比較 由表7 可知，28～30 周齡，Logistic 擬合值與實際觀測值差值最大，其次是Gompertz 模型，Von Bertalanffy 模型最接近觀測值。3 種模型在28～30 周齡和35～39 周齡的擬合值大于觀測值，31～34 周齡和40 周齡的擬合值小于觀測值。綜上，Von Bertalanffy 模型擬合累計產蛋數(shù)效果最佳，優(yōu)于Logistic 模型和Gompertz 模型，為最佳擬合模型。

表7 烏蒙鳳雞累計產蛋數(shù)模型擬合值與實際觀測值

3 討論

影響雞產蛋率的因素包括遺傳因素、環(huán)境衛(wèi)生、飼養(yǎng)管理等方面[11]，產蛋曲線能更好地反映雞群的產蛋規(guī)律性。烏蒙鳳雞作為新近發(fā)現(xiàn)的種質資源，其產蛋性狀的可塑性可提高烏蒙鳳雞的選育效果。許多畜牧科技工作者應用楊寧模型、分室模型和伍德模型對京海黃雞[12]、二郎山山地雞[13]、淅川烏骨雞[14]和豫粉1 號蛋雞祖代D 系種雞[15]等家禽品種（配套系）開展了產蛋曲線擬合研究及分析，這對提高各雞種產蛋性能具有重要意義。本研究以零世代烏蒙鳳雞為研究對象，采用不同模型進行產蛋率和累計產蛋數(shù)的曲線擬合，結果發(fā)現(xiàn)楊寧模型對烏蒙鳳雞產蛋率曲線擬合效果最好，這與李鴻飛等[16]研究“豫粉1 號”商品代雞（19～63 周齡）和陸雪林等[17]研究美國七彩山雞27～55 周齡和黑化山雞30～55 周齡的產蛋曲線結果存在差異。原因一方面可能在于不同品種或品系適用于不同模型，另一方面在于數(shù)據處理和生長期長短存在差異而導致擬合結果不同。李鴻飛等[18]對“豫粉1 號”配套系商品代蛋雞（21～63 周齡）和王志成等[19]對淮南麻黃雞和文昌雞雜交后代母雞（23～58 周齡）進行的曲線擬合結果與本研究相一致。在實際生產中，為有效提高雞群產蛋性能，需要對雞群的累計產蛋數(shù)以及產蛋拐點充分了解。在雞的產蛋曲線擬合研究方面，有學者已對京海黃雞[20]、蘆花雞H 系雞[21]、彭縣黃雞[22]和舊院黑雞[23]進行了累計產蛋數(shù)曲線擬合分析，擬合結果與本研究的擬合結果一致，以Von Bertalanffy 模型的擬合效果最佳。綜上，在實際生產中可用楊寧模型和Von Bertalanffy 模型提高烏蒙鳳雞產蛋性能、選育進展和保種擴繁效果。隨著數(shù)學模型不斷優(yōu)化，產蛋曲線描述更具準確性，其預測結果也將越來越貼近真實情況，為提高烏蒙鳳雞的產蛋水平提供更好的方法和手段。

4 結論

本實驗運用伍德模型、分室模型、楊寧模型對烏蒙鳳雞零世代產蛋率曲線進行擬合，結合Logistic、Gompertz 和Bertalanffy 3 種模型對其累計產蛋數(shù)進行擬合分析，結果表明，楊寧模型對烏蒙鳳雞產蛋率曲線擬合效果最好，Von Bertalanffy 模型擬合累計產蛋數(shù)曲線效果最佳。