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

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

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

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

1 材料與方法

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

1.2 實驗方法

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

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

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

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

表2 產(chǎn)蛋率擬合模型

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

2 結(jié)果與分析

2.1 3 種模型產(chǎn)蛋率曲線擬合結(jié)果 3 種數(shù)學(xué)模型擬合得到的各模型周產(chǎn)蛋率參數(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 種產(chǎn)蛋率擬合模型參數(shù)值

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

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

2.3 3 種模型累計產(chǎn)蛋數(shù)曲線擬合結(jié)果 烏蒙鳳雞零世代累計產(chǎn)蛋數(shù)3 種非線性模型擬合參數(shù)估計值和擬合度見表6。把參數(shù)估計值A(chǔ)、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 模型對烏蒙鳳雞累計產(chǎn)蛋數(shù)擬合度分別為0.989、0.996、0.998，以Von Bertalanffy 模型擬合度最高，用其計算出的拐點產(chǎn)蛋周齡為33.950 周，拐點產(chǎn)蛋數(shù)為20.253個，極限產(chǎn)蛋數(shù)為68.354 個。Logistic 模型拐點產(chǎn)蛋周齡為35.236 周，拐點產(chǎn)蛋數(shù)為25.904 個，極限產(chǎn)蛋數(shù)為51.808 個。Gompertz 模型的拐點產(chǎn)蛋周齡為34.345周，拐點產(chǎn)蛋數(shù)為22.191 個，極限產(chǎn)蛋數(shù)為60.316 個。

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

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

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

3 討論

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

4 結(jié)論

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