乌蒙凤鸡零世代产蛋曲线拟合研究

王天松，王 梅，张 勇*，郭徵力，莫先艇，骆科印，杨 红，曾姗姗，左 丽

（1.贵州大学动物科学学院，高原山地动物遗传育种与繁殖教育部重点实验室，贵州贵阳 550025；2.六盘水市农业农村局，贵州六盘水 553001；3.毕节市畜牧站，贵州毕节 551700）

乌蒙凤鸡原产于贵州省六盘水市，其中心产区海拔高度介于1 400～1 900 m，产区特殊的喀斯特地形地貌减少了鸡种间的基因交流，乌蒙凤鸡基因较为纯合，以独特的凤头性状而保存下来。乌蒙凤鸡是以肉用性能为主的地方鸡种，具有体型适中，雌雄鸡凤头明显、耐粗饲、抗逆性强、觅食能力强及突出的高原特色等特点[1]。目前，乌蒙凤鸡雌鸡群体羽色多为麻羽和黄羽，少数为白羽和黑羽，具有明显的冠羽特征，初产日龄在181～210 d。乌蒙凤鸡经保种扩繁和推广，已申报获得国家农产品地理标志登记保护及农产品地理标志登记证书。

产蛋量是鸡最重要的数量性状，受品种、年龄、光照、营养水平等诸多因素影响，然而产蛋量变化通常遵循一定规律，服从一定曲线分布[2]。如果将日产蛋率及累计产蛋数标记在坐标轴上，将各点连接起来，则可构成产蛋曲线，产蛋曲线数学模型可从动态角度展示产蛋率随时间变化的规律[3]。鸡产蛋曲线不仅反映其产蛋率高低，还可以反映鸡不同时期的饲养管理水平，可以指导生产，使鸡在标准体重下发挥最佳产蛋水平[4]。乌蒙凤鸡作为新发现的遗传资源，到目前为止，相关研究工作报道较少，因此本实验开展乌蒙凤鸡的产蛋曲线拟合研究，探索乌蒙凤鸡产蛋规律，为进一步加大其选育力度、提高产蛋率、促进其保种扩繁奠定基础。

1 材料与方法

1.1 实验材料 实验在六盘水金凤栖自然牧业有限公司乌蒙凤鸡保种扩繁场进行。选取零世代乌蒙凤鸡母鸡360 羽，相同条件下单笼饲养，母鸡开产后分别记录28～40 周龄产蛋数。

1.2 实验方法

1.2.1 实验期日粮营养水平 实验鸡不同生长阶段的主要营养需要见表1，根据产蛋期采食量、产蛋率、蛋重、蛋重周增长及体重变化等调整饲料配方。实验期间，定期观察鸡群状况，做好死亡数、产蛋率、存栏数和产蛋量等相关记录。

表1 乌蒙凤鸡营养需要量

1.2.2 产蛋性能统计 记录实验鸡的开产日龄（日产蛋率达5%时）、产蛋量，汇总周产蛋量和存栏鸡数，统计各周龄产蛋率及累计产蛋数，计算公式如下：周产蛋率（%）=各周产蛋量/统计期内总饲养母禽只数×100%；周累计产蛋数（个）=周累计产蛋个数/存栏母鸡数。

1.3 数据处理 采用Excel 2010 进行数据初处理，用SPSS 22.0 软件计算产蛋率和累计产蛋数模型拟合参数值。伍德模型、杨宁模型、分室模型[5-7]对产蛋率曲线拟合公式见表2。Logistic、Gompertz 和Von Bertalanffy模型[8-10]对累计产蛋数曲线拟合公式见表3。

表2 产蛋率拟合模型

表3 累计产蛋数拟合模型

2 结果与分析

2.1 3 种模型产蛋率曲线拟合结果 3 种数学模型拟合得到的各模型周产蛋率参数估计值见表4。把参数估计值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 种产蛋率拟合模型参数值

2.2 产蛋率观测值与拟合值比较 各周龄产蛋率观测值和模型拟合值比较结果见表5。乌蒙凤鸡28 周龄产蛋率为10%，28～32 周龄产蛋率持续上升，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 种模型累计产蛋数曲线拟合结果 乌蒙凤鸡零世代累计产蛋数3 种非线性模型拟合参数估计值和拟合度见表6。把参数估计值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 模型对乌蒙凤鸡累计产蛋数拟合度分别为0.989、0.996、0.998，以Von Bertalanffy 模型拟合度最高，用其计算出的拐点产蛋周龄为33.950 周，拐点产蛋数为20.253个，极限产蛋数为68.354 个。Logistic 模型拐点产蛋周龄为35.236 周，拐点产蛋数为25.904 个，极限产蛋数为51.808 个。Gompertz 模型的拐点产蛋周龄为34.345周，拐点产蛋数为22.191 个，极限产蛋数为60.316 个。

表6 乌蒙凤鸡3 种累计产蛋数拟合模型参数值

2.4 累计产蛋数观测值与拟合值比较 由表7 可知，28～30 周龄，Logistic 拟合值与实际观测值差值最大，其次是Gompertz 模型，Von Bertalanffy 模型最接近观测值。3 种模型在28～30 周龄和35～39 周龄的拟合值大于观测值，31～34 周龄和40 周龄的拟合值小于观测值。综上，Von Bertalanffy 模型拟合累计产蛋数效果最佳，优于Logistic 模型和Gompertz 模型，为最佳拟合模型。

表7 乌蒙凤鸡累计产蛋数模型拟合值与实际观测值

3 讨论

影响鸡产蛋率的因素包括遗传因素、环境卫生、饲养管理等方面[11]，产蛋曲线能更好地反映鸡群的产蛋规律性。乌蒙凤鸡作为新近发现的种质资源，其产蛋性状的可塑性可提高乌蒙凤鸡的选育效果。许多畜牧科技工作者应用杨宁模型、分室模型和伍德模型对京海黄鸡[12]、二郎山山地鸡[13]、淅川乌骨鸡[14]和豫粉1 号蛋鸡祖代D 系种鸡[15]等家禽品种（配套系）开展了产蛋曲线拟合研究及分析，这对提高各鸡种产蛋性能具有重要意义。本研究以零世代乌蒙凤鸡为研究对象，采用不同模型进行产蛋率和累计产蛋数的曲线拟合，结果发现杨宁模型对乌蒙凤鸡产蛋率曲线拟合效果最好，这与李鸿飞等[16]研究“豫粉1 号”商品代鸡（19～63 周龄）和陆雪林等[17]研究美国七彩山鸡27～55 周龄和黑化山鸡30～55 周龄的产蛋曲线结果存在差异。原因一方面可能在于不同品种或品系适用于不同模型，另一方面在于数据处理和生长期长短存在差异而导致拟合结果不同。李鸿飞等[18]对“豫粉1 号”配套系商品代蛋鸡（21～63 周龄）和王志成等[19]对淮南麻黄鸡和文昌鸡杂交后代母鸡（23～58 周龄）进行的曲线拟合结果与本研究相一致。在实际生产中，为有效提高鸡群产蛋性能，需要对鸡群的累计产蛋数以及产蛋拐点充分了解。在鸡的产蛋曲线拟合研究方面，有学者已对京海黄鸡[20]、芦花鸡H 系鸡[21]、彭县黄鸡[22]和旧院黑鸡[23]进行了累计产蛋数曲线拟合分析，拟合结果与本研究的拟合结果一致，以Von Bertalanffy 模型的拟合效果最佳。综上，在实际生产中可用杨宁模型和Von Bertalanffy 模型提高乌蒙凤鸡产蛋性能、选育进展和保种扩繁效果。随着数学模型不断优化，产蛋曲线描述更具准确性，其预测结果也将越来越贴近真实情况，为提高乌蒙凤鸡的产蛋水平提供更好的方法和手段。

4 结论

本实验运用伍德模型、分室模型、杨宁模型对乌蒙凤鸡零世代产蛋率曲线进行拟合，结合Logistic、Gompertz 和Bertalanffy 3 种模型对其累计产蛋数进行拟合分析，结果表明，杨宁模型对乌蒙凤鸡产蛋率曲线拟合效果最好，Von Bertalanffy 模型拟合累计产蛋数曲线效果最佳。