吉林芦花鸡生长曲线拟合及体重与体尺的相关性

柳俭强， 刘 臣， 张芳毓， 武 斌， 曹 阳， 肖 成， 李 欣， 刘艳光，2， 张立春*

(1.吉林省农业科学院 畜牧分院， 吉林 公主岭 136100； 2.延边大学 农学院， 吉林 延吉 130021)

吉林芦花鸡是吉林省优良地方家禽品种之一，属于肉蛋兼用型鸡种，因全身羽毛黑白相间、斑纹宽窄均匀，俗称为“芦花鸡”。主要分布于吉林公主岭市、双辽市、永吉县、乾安县等[1]。因其具有良好的繁殖性能，肉质鲜美、营养价值高等优点，深受消费者欢迎。目前，吉林芦花鸡选育工作已取得一定进展，但鉴于其属地方品种，在生长速度上与培育品种和引进品种还存在一定差距。

动物生长发育因受遗传、环境等因素影响，生长规律呈非线性特点。对于一些生长发育指标需要通过非线性模型进行拟合，选择出理想的生长模型，才能相对真实地体现生长发育规律，从而预测家禽未来的生长走势，进一步提高品种选育的效果。研究家禽生长发育规律目前最为常用的是Logistic、Gompertz和 Bertalanffy 3种非线性曲线模型[2-4]，已广泛应用到实际生产中，但其应用于吉林芦花鸡生长发育规律的研究鲜见报道。为此，笔者运用以上3种非线性生长模型对吉林芦花鸡生长曲线进行拟合，并对体重与体尺指标进行相关性分析，旨在为吉林芦花鸡品种资源的保存与选育提供参考。

1材料与方法

1.1材料

吉林芦花鸡雏鸡公母鸡各400只，由吉林省农业科学院畜牧分院养禽研究室提供，体态健康。

1.2方法

1.2.1饲养管理公母鸡每组设置20个重复，每个重复20只，在相同环境条件下进行规范的饲养管理，并按照规定的免疫程序进行免疫，自由采食、饮水。

1.2.2测定指标每个重复随机挑10只，每组共200只，分别在0周龄、1周龄、3周龄、5周龄、7周龄、9周龄、11周龄、13周龄、15周龄、17周龄和21周龄测定体重，并在120日龄时测定体尺，包括体斜长、胫长、龙骨长、胫骨长、胸深、胸宽、盆骨宽和胫围8项指标。

1.2.3生长曲线拟合根据称重结果，选用Logistic、Gompertz和 Bertalanffy非线性曲线模型对体重生长进行拟合，具体表达式与参数见表1。

1.3统计分析

利用Excel整理数据，用SPSS 24.0建立数据库，拟合3种模型曲线A、B、k的最优估计值，建模后根据拟合度R2进行评价，并得出拐点体重、拐点周龄和最大周增重；对体重和体尺进行相关性分析。结果用平均值±标准差表示。

表1Logistic等3种非线性模型曲线表达式与参数

注：Wt为t周龄时体重，A为极限生长量，B为参数，e为自然对数的底，k为生长速率，t为周龄。

Note:Wtis the weight at t weeks of age,Ais the limit growth amount,Bis a parameter, e is the base of the natural logarithm,kis the growth rate, andtis the age of weeks.

2结果与分析

2.1吉林芦花鸡的生长发育情况

从表2看出，吉林芦花鸡0～21周龄公母鸡的体重分别在40.08～2 953.84 g和39.90～2 110.24 g。从图1可知，公、母鸡在前4周生长较为缓慢，差异不明显；整个生长期公鸡的生长速度明显快于母鸡，17周后增重明显减慢。

2.2生长拟合曲线

从表3可知，吉林芦花鸡公鸡的3种生长曲线模型拟合度(R2)均超过0.99，其中，Gompertz曲线拟合度最高，达到0.995；拐点体重为1 318.957 g，拐点周龄为10.14周，最大周增重为188.61 g。母鸡Gompertz和Bertalanffy生长曲线模型的拟合度均为0.992，拐点体重、周龄和最大周增重差异不大，均可较好模拟母鸡的生长曲线。同一个生长曲线模型的公鸡拐点周龄晚于母鸡、体重高于母鸡，说明公鸡生长周期比母鸡长。

表2吉林芦花鸡0～21周龄的体重

Table 2 Weights of Jilin Luhua chicken from 0 to 21 weeks of age g

周龄 Week age公鸡 Cock母鸡 Hen040.08±3.7639.90±3.35167.76±8.3558.62±8.643214.52±38.00196.56±39.765414.37±88.14382.46±88.517802.63±161.25672.57±134.4391 151.19±217.89937.76±186.86111 538.47±342.711 189.76±237.84131 789.56±176.331 382.08±137.89172 372.80±231.611 839.40±177.98212 953.84±241.992 110.24±201.39

图10～21周龄吉林芦花鸡累积生长曲线(左)与周增重曲线(右)

Fig.1 Cumulative growth curve (left) and weight gain curve (right) of Jilin Luhua chicken aged 0 to 21 weeks

表3 3种生长曲线模型的拟合度和参数估计值

2.3吉林芦花鸡生长曲线与3种生长曲线模型的比较

从图2和表4可知，公鸡初生重Logistic曲线的拟合值与实测值偏差较大，拐点体重与实测值较接近，整个生长周期与另2个曲线相比较也有一定差异；0～4周Gompertz曲线的拟合值与实测值较接近，4～11周与Bertalanffy相比差异不大，11周之后模拟值进一步提高，与实测值更为接近；Bertalanffy在0～4周与实测值差异不大并在4周后逐渐缩小与实测值的差距。母鸡Logistic曲线的拟合值与初生重差异大，不适合预测母鸡的生长发育规律，Gompertz、Bertalanffy曲线的拟合值与整体实测值较接近。

图2 吉林芦花鸡公鸡(左)、母鸡(右)生长曲线与3种拟合曲线

表4 吉林芦花鸡3种生长曲线模型拟合值与实测值

2.4120日龄吉林芦花鸡的体重、体尺及其相关性

2.4.1体重与体尺指标从表5看出，公鸡的体重、体斜长、龙骨长、胸深、胸宽、骨盆宽、胫长、胫围等7个性状值均高于母鸡，分别比母鸡高43.59%、11.53%、14.34%、13.85%、20.04%、31.12%、27.97%和18.55%。

2.4.2 相关性从表6看出，公鸡、母鸡的体重与7个体尺性状间均存在不同程度的相关性。公鸡的体重与胫长、胸深、胫围均呈极显著正相关；与龙骨长、骨盆宽呈显著正相关。母鸡的体重与体斜长、胫长、龙骨长、胸深和胫围均呈极显著正相关；与骨盆宽呈显著正相关。

表5120日龄吉林芦花鸡的体重与体尺指标

表6吉林芦花鸡体重与体尺的相关系数

注: *表示在0.05水平(双侧)上显著相关，**表示在0.01水平(双侧)上显著相关。

Note: * indicates a significant correlation at the 0.05 level (both sides) and ** indicates a significant correlation at the 0.01 level (both sides).

3结论与讨论

3.1吉林芦花鸡公鸡、母鸡的生长发育情况

生长曲线拟合分析实质是把不同时期体重资料提供的信息整合成为少数几个参数，是研究家禽生长发育规律的主要方法之一[5-6]。通过分析不同性别鸡的生长曲线，确定重要拐点及日龄，结合合理的日粮营养，使其发挥出最佳的生产潜力。运用Logistic、Gompertz和Bertalanffy模型均能很好模拟出不同性别吉林芦花鸡的生长曲线，拟合度均超过0.99。其中，Gompertz模型的拟合度达到0.995；公母鸡拐点周龄分别为10.14周和8.98周，晚于ROSS308[7]；与Rhode Island Reds[8]、盐津乌骨鸡[9]、腾冲雪鸡[10]、独龙鸡[11]的生长曲线拟合一致；拐点体重与实测值更接近，模拟效果最好，说明Gompertz模型适合描述禽类早期生长迅速的生长过程[12]。边鸡、武定鸡的拟合以Bertalanffy模型最佳[13-14]；对早期生长速度慢、拐点日龄迟的藏鸡Logistic模型更适合[15]。鸡品种、生长周期和生长速度不同造成生长曲线的最佳模型不一,拟合吉林芦花鸡生长曲线的最佳模型是Gompertz。

研究结果表明，吉林芦花鸡公鸡与母鸡的体重生长曲线不同，前4周公母鸡生长缓慢，无明显差异。与浦东鸡公、母鸡3周龄前体重相一致的结果[6]略有不同，可能由于品种差异导致。吉林芦花鸡从第5周开始公母鸡体重差异变化显著，与拉萨白鸡在42日龄后公母体重差异极显著的结果[4]相一致。吉林芦花鸡母鸡生长速度高峰比公鸡早2周左右，与独龙鸡在7周龄周增重达到最大的结果[11]相一致。因此，吉林芦花鸡在饲养过程中建议公母要分开饲养，并在生长高峰期提高营养和饲养管理水平，使各自发挥出最大生产潜力。

3.2吉林芦花鸡体重与体尺间的相关性

体重与体尺指标是重要的表型性状，研究彼此之间的相关关系，对于提高品种选育和经济效益有着重要的意义。研究结果表明，120日龄吉林芦花鸡的体重与体斜长、胫长等7个体尺指标间存在不同程度的正相关关系，与对略阳乌鸡[16]、Ⅰ系雁荡土鸡[17]的研究结果一致。吉林芦花鸡公母鸡体重与胫长、胸深均呈现极显著正相关，与对四名香鸡[18]的研究结果一致，吉林芦花鸡早期选育可将体尺性状的相关性作为参考指标，其中胫长和胸深是重点，公鸡可增加胫围作为参考指标，母鸡可增加体斜长、龙骨长和胫围作为参考指标。