快慢羽太行鸡生长规律和Logistic生长曲线拟合分析

张乐超，张秀玲，刘春杨，刘小辉，李丽莎，魏彦辉，李兰会*，李祥龙

（1.河北农业大学动物科技学院，河北保定 071 000；2. 河北科技师范学院动物科技学院，河北秦皇岛 066004）

快慢羽太行鸡生长规律和Logistic生长曲线拟合分析

张乐超1，张秀玲1，刘春杨1，刘小辉1，李丽莎1，魏彦辉1，李兰会1*，李祥龙2

（1.河北农业大学动物科技学院，河北保定 071 000；2. 河北科技师范学院动物科技学院，河北秦皇岛 066004）

为了解快慢羽太行鸡的生长发育规律，对其体重和体尺指标及其相关性以及Logistic拟合生长曲线进行分析。结果表明:慢羽公鸡生长最快，其次为快羽公鸡，快羽和慢羽母鸡相似最低，其Logistic生长曲线分别为Y=1.760/（1+13.445e-0.034t）、Y=1.519/（1+14.066e-0.036t）、Y=1.377/（1+11.503e-0.031t）和Y=1.314/（1+10.286e-0.031t）。当快羽母鸡生长到77～79日龄体重达到0.698 kg左右，慢羽母鸡生长到74～76日龄体重达到0.657 kg左右，快羽公鸡生长到72～74日龄体重达到0.760 kg，慢羽公鸡生长到75～77日龄体重达到0.880 kg，此时有最大的生长速率。

快慢羽；太行鸡；体重；相关性；Logistic生长曲线

太行鸡主产于河北境内，分布于太行山脉沿线的保定、石家庄、邢台、邯郸等地区，是长期自然选育适应当地环境特点形成的蛋肉兼用型地方品种。因其腿纤细、体型瘦小如柴，而且常栖息于柴禾之上，被当地百姓俗称为“柴鸡”[1]。当地群众长期粗放饲养形成独特的品种特征--体型矮小、结构匀称，羽毛紧凑，骨骼纤细，耐粗饲、抗逆性强，遗传性能稳定、抗病力强、有良好适应性，并且所产蛋、肉品质良好[2]。

按照生长速度划分，鸡的羽毛可以分为快、慢羽（亦称早羽和迟羽）2种，标准为翼羽和尾羽长出的快慢和早迟。控制这一性状的基因有位于常染色体上的T 基因（T、tS、t）和位于Z染色体上控制羽速的K 基因（Kn、Ks、K、k）。羽速基因在早期可以影响雏鸡主覆翼羽、尾羽的生长，也会对一些组织器官产生影响，进而影响鸡整体的生长发育，并且其对鸡的生产性能、抗病力、繁殖性能都有影响。

体重和体尺是非常重要的生长发育指标，是鸡生长发育的基本特征，通常直接或间接影响生产性能[3]，河田鸡[4]、大骨鸡[5]、坝上长尾鸡[6]等品种的体重与体尺之间均存在着一定相关性。本研究对快慢羽太行公母鸡的体重和体尺性状变化以及体重体尺间的相关性进行分析，以期为太行鸡选种选育提供依据。

生长曲线的分析和拟合是研究畜禽生长发育规律的主要方法之一。不同个体生长曲线的特征被视为重要的综合经济指标，用于家禽品种的选育，可为选育适合不同要求的家禽品种提供依据[7]。安义瓦灰鸡[8]、坝上长尾鸡[9]等品种鸡的体重拟合分析都有相关研究，但对于快慢羽太行公母鸡体重拟合分析还未见报道，快慢羽基因是一对控制羽毛生长速度的伴性遗传基因，对鸡早期生长速度和体重、羽毛生长发育、产蛋性状、抗病力等生产性状都有一定的影响。本研究采用Logistic非线性数学模型对太行鸡体重累积生长曲线进行拟合与分析，对比快慢羽太行公母鸡生长规律差异，以期用于指导太行鸡的生产实践，为提高鸡群的生产性能、建立更科学的饲养管理体系提供依据。

1 材料与方法

1.1 实验材料 随机选取112只太行鸡体重体尺数据进行统计分析，其中34只快羽母鸡，26只慢羽母鸡，32只快羽公鸡，20只慢羽公鸡。

1.2 实验方法 按常规饲养管理方式，执行常规免疫程序和规范的饲养管理，实行自由采食、自由饮水。早上空腹称体重，按快慢羽公母分别计算平均体重。

根据初步模拟分析，选用Logistic数学模型拟合体质量的累积生长过程。模型表达式及相关指标计算式见表1。参数A为极限体重（成熟体重），k为瞬时相对生长率，B为常数尺度，w为体重。比较生长差异，主要看初始体重、拐点体重、拐点日龄这3个参数。

表1 曲线模型

1.3 统计分析 采用Excel 2007整理数据，利用SPSS19.0的One-Way ANOVA分析不同品系鸡体重和体尺指标的差异性，体重和体尺指标间的相关性采用SPSS19.0统计软件中的双变量相关分析。采用SPSS19.0统计软件回归分析程序中的非线性回归分析子程序，得到相应各个参数（A，B，k）的估计值与拟合度（R2值）[7]，相关参数（拐点周龄、拐点体质量、初始体质量等）计算与图形的绘制用Excel软件完成[10]，结果表示为平均值±标准差。

2 结 果

2.1 太行鸡体重体尺性状变化 由表2可知，30日龄慢羽公鸡显著重于快羽母鸡（P＜0.05），51日龄开始快羽和慢羽公鸡分别显著（P＜0.05）和极显著（P＜0.01）重于母鸡，73～139日龄慢羽公鸡极显著重于快羽公鸡（P＜0.01）。0～139日龄快慢羽母鸡间体重没有显著差异（P＞0.05）。

不同日龄太行鸡羽速类型和性别的体尺变化特征（表3）分述如下。体长变化：107～139日龄慢羽公鸡显著高于快羽公鸡（P＜0.05），快慢羽母鸡间差异不显著（P＞0.05）。龙骨长变化：快慢羽公鸡间龙骨长差异不显著（P＞0.05）；73和139日龄慢羽母鸡显著低于快羽母鸡（P＜0.05）。胫围变化：51～107日龄慢羽公鸡显著高于快羽公鸡（P＜0.05），139日龄快慢羽公鸡间差异不显著（P＞0.05）；51日龄慢羽母鸡显著高于快羽母鸡（P＜0.05），73～139日龄快慢羽母鸡间差异不显著（P＞0.05）。胫长变化：107～139日龄慢羽公鸡显著高于快羽公鸡（P＜0.05），快慢羽母鸡间差异不显著（P＞0.05）。胸宽变化：107～139日龄慢羽公鸡显著高于快羽公鸡（P＜0.05），快慢羽母鸡间差异不显著（P＞0.05）。胸深变化： 107日龄慢羽公鸡显著高于快羽公鸡（P＜0.05），139日龄慢羽公鸡与快羽公鸡差异不显著（P＞0.05）；快慢羽母鸡间差异不显著（P＞0.05）。

2.2 太行鸡体重体尺指标相关性 由表4可知，体重和体尺指标均呈极显著正相关（P＜0.01）。体重与体长的相关性最大为0.968，与胫围的相关性最小为0.880；龙骨与体长的相关性最大为0.967，与胫围的相关性最小为0.871；胫围与胫长的相关性最大为0.917，与胸宽的相关性最小为0.849。胫长与胸深的相关性最大为0.930，与胸宽的相关性最小为0.868。胸宽与体长的相关性最大为0.943。胸深与体长的相关性最大为0.957，与胫围的相关性最小为0.880。体长与胫围相关性最小为0.882。

表2 不同日龄快慢羽太行公母鸡体重变化 kg

表3 快慢羽太行公母鸡体尺性状变化

表4 太行鸡体重和体尺指标相关性

2.3 快慢羽太行公母鸡累积体重的Logistic模型比较。

2.3.1 快慢羽太行公母鸡Logistic模型参数 由表5可知，慢羽公鸡的极限体重最大可达1.760 kg，慢羽母鸡极限体重最小为1.314 kg，快羽母鸡和快羽公鸡的极限体重分别为1.377 kg和1.519 kg。快羽公鸡的拐点日龄最早为73.438 d，快羽母鸡最晚为78.794 d。慢羽公鸡拐点体重最大为0.880 kg，慢羽母鸡最小为0.657 kg。最大日增重，由大到小依次为慢羽公鸡、快羽公鸡、快羽母鸡和慢羽母鸡。

2.3.2 快慢羽太行公母鸡累积体重拟合生长曲线快羽母鸡、慢羽母鸡、快羽公鸡和慢羽公鸡的体重模型分别为Y=1.377/（1+11.503e-0.031t）、Y=1.314/（ 1+10.286e-0.031t）、Y=1.519/（1+1 4.0 6 6 e-0.0 3 6 t） 和 Y=1.7 6 0/（1+13.445e-0.034t），拟合生长曲线见图1，横坐标表示日龄，纵坐标表示体重（ kg）。

结合图1和表3可知，当快羽母鸡77～79日龄时体重达到0.698 kg左右，慢羽母鸡74～76日龄时体重达到0.657 kg左右，快羽公鸡72～74日龄时体重达到0.760 kg，慢羽公鸡75～77日龄时体重达到0.880 kg，此时有最大的生长速率，骨骼和肌肉的发育最为迅速，体重的累积生长最快。

图1 快慢羽太行公母鸡Logistic生长拟合曲线

3 讨 论

3.1 太行鸡体重和体尺 在育成初期，羽速对鸡的生长发育有一定的影响。随着鸡日龄的增加，这种影响逐渐转变为性别因素占据主导。在育成期，公鸡发育更加迅速，母鸡开始为产蛋期做准备，公母鸡体型差异变大。通过太行鸡生长指标分析发现：慢羽公鸡生长最快，其次为快羽公鸡，快羽和慢羽母鸡生长最慢，并且二者没有差异。慢羽公鸡从51日龄开始，体长、龙骨长、胫围、胫长、胸宽和胸深均超出快羽和慢羽母鸡，107日龄的体长、胫围和胸深超出快羽公鸡，139日龄的胫长和胸宽超出快羽公鸡。相关性分析表明，太行鸡体重和体尺指标两两间均有极显著正相关性，其中体长、龙骨和胸深与体重的相关性较大。

3.2 生长规律的差异 目前，许多学者对禽类的生长规律做了细致的研究，通过各种数学模型预测各品种不同生长阶段的体重[11-12]。李莉等[13]对桂香鸡生长规律的研究表明，在预测各阶段重要生长拐点时，以Logistic模型较为准确。生长曲线被看作是重要的综合经济指标，用于家禽品种的选育，可为选育适合不同要求的家禽品种提供依据。本实验发现以Logistic模型预测太行鸡的拐点周龄在10～11周，拐点体重在650～880 g。韩占兵等[14]研究发现，丝毛乌骨鸡的拐点周龄为8.75周，拐点体重为510.8 g，拐点周龄小于太行鸡，拐点体重也小于太行鸡。付睿琦等[15]研究北京油鸡发现，公鸡拐点周龄和拐点体重为10.8周和986.3 g，母鸡拐点周龄和拐点体重为7.8周和639.4 g，公鸡的拐点日龄和母鸡的拐点体重与太行鸡差异不大。说明河北太行鸡生长发育特点与北京油鸡相似，推测其营养需求相似。

表5 Logistic生长模型分析的参数估计和各指标计算值

早期的拐点体重和拐点日龄具有重要的市场影响，有助于提高该阶段的饲养管理水平，节约饲养成本，增加经济效益。本实验研究表明，在30日龄前快慢羽太行公母鸡体重没有发生显著差异，在30～51日龄慢羽公鸡与快慢羽母鸡体重差异到达显著水平，51日龄后快慢羽公鸡与快慢羽母鸡达到极显著水平，但快羽和慢母鸡间体重差异不大。Logistic模型曲线拟合结果表明，快羽母鸡、慢羽母鸡、快羽公鸡和慢羽公鸡生长拐点分别为78.8、75.2、73.4 d和76.4 d，此时分别达到了生长高峰期，因而，在饲养过程中应特别注意快羽母鸡79 d、慢羽母鸡76 d、快羽公鸡74 d、慢羽公鸡77 d前的饲料营养的全价性和平衡性，为太行鸡的良好生长发育打下坚实的基础，使其达到体重标准。快羽母鸡、慢羽母鸡、快羽公鸡和慢羽公鸡的拐点体重分别为0.689、0.657、0.760 kg和0.880 kg，可为太行鸡的饲养和出栏日龄提供参考。

4 结 论

慢羽公鸡的生长速度最快，其次为快羽公鸡，快慢羽母鸡生长速度最慢。根据Logistic拟合生长曲线，快慢羽母鸡78日龄、慢羽公鸡76日龄、快羽公鸡73日龄为体重增加速度由快到慢的转折点。因此，在饲养过程中在各品系的最快生长速度时期应特别注意饲料营养的平衡性，使达到体重标准，为太行鸡的良好生长发育和发挥高水平产蛋性能打下坚实的基础。

[1] 刘树常, 张永镇, 申明, 等. 河北省畜牧志[M]. 北京:农业出版社, 1993, (7): 137-139.

[2] 张秀玲, 李兰会, 李祥龙, 等. 河北太行鸡遗传资源调研报告[J]. 山东农业大学学报(自然科学版), 2015, 46(5): 725-729.

[3] 秦巧梅. 鲁西斗鸡体尺体重及屠宰性状间的典型性相关分析[J]. 中国畜牧兽医, 2012, 39(6): 153-156.

[4] 朱志明, 杨武, 缪中纬, 等. 河田鸡体重,体尺及屠宰性状测定及其相关性分析[J]. 中国农学通报, 2012, 28(35): 25-28.

[5] 包艳波, 徐日福, 王志贤, 等. 大骨鸡体重和体斜长生长曲线拟合及其相关分析[J]. 中国畜牧杂志, 2013, 49(3): 5-8.

[6] 刘小辉, 李祥龙, 逯春香, 等. 基坝上长尾鸡体重体尺性状相关及回归分析[J]. 甘肃畜牧兽医, 2015, 45(2): 42-45.

[7] 戴国俊, 王金玉, 杨建生, 等. 应用统计软件SPSS 拟合生长曲线方程[J]. 畜牧与兽医, 2006, 38(9): 28-30.

[8] 谢明贵, 谢金防, 刘林秀, 等. 安义瓦灰鸡早期生长规律研究[J]. 中国家禽, 2012, 34( 7): 16-22.

[9] 刘小辉, 李祥龙, 李 杰, 等. 坝上长尾鸡育雏育成期生长曲线拟合研究[J]. 浙江农业学报, 2015, 27(5): 746-750.

[10] 唐五湘. Excel 在统计中的应用[M]. 北京:电子工业出版社, 2002: 71-189.

[11] Aggrrey S E. Comparison of three nonlinear and spline regression models for describing chickengrowth curve[J]. Poult Sci, 2002,81(12):1782-1788.

[12] Mignon-grasteau S, pilest M, Varona L, et al. Genetic analysis ofgrowth curve parameters for male and female chickens resulting from selection on shape ofgrowth curve[J]. J Anim Sci, 2000, 78(10):2512-2524.

[13] 李莉, 许晟玮, 梁明振, 等. 桂香鸡生长规律的研究[J].动物营养学报, 2014, 26(7): 1787-1795.

[14] 韩占兵, 张立恒, 刘健, 等. 丝毛乌骨鸡生长规律和生长曲线拟合研究[J]. 畜牧与兽医, 2016, 48(1): 53-46.

[15] 付睿琦, 赵桂平, 刘冉冉, 等. 北京油鸡生长曲线拟合与分析[J]. 中国家禽, 2014, 36(3): 6-9.

Abstract:To study the growth and development character of the early and late feathering Taihang chickens, body weight and size and their correlation and Logistic growth curve fitting model were analyzed. Results showed that the late feathering male Taihang chicken had the fastest growth rate, and the early feathering male Taihang chickens growed faster, while the early and late feathering female Taihang chickens had the slowest growth rate. The Logistic growth ftting curve were Y=1.760/(1+13.445e-0.034 t), Y=1.519/(1+14.066e-0.036 t), Y=1.377/(1+11.503e-0.031 t) and Y=1.314/(1+10.286e-0.031 t) respectively. The weight and the age of the four various chickens that had the highest growth rate were predicted to be 0.698 kg and 77-79 day,0.657 kg and 74-76 day,0.760kg and 72-74 day,0.880 kg and 75-77 day respectively.

Growth Character and Logistic Growth Curve Fitting Model of the Early and Late Feathering Taihang Chickens

ZHANG Le-chao1, ZHANG Xiu-ling1,LⅠU Chun-yang1, LⅠU Xiao-hui1, LⅠ Li-sha1, WEⅠ Yan-hui1, LⅠ Lan-hui1*LⅠ Xiang-long2
（1.College of Animal Science and Technology Agricultural University of Hebei，Hebei Baoding 071000，China；2.College of Animal Science and Technology, Hebei Normal University of Science&Technology，Hebei Qinhuangdao 066004，China）

Early and late feathering; Taihang chickens; Body weight; Correlation; Logistic growth curve

S831.5

A

10.19556/j.0258-7033.2017-01-017

2016-06-21 ；

2016-08-30

河北省现代农业技术体系蛋鸡产业创新团队（HBCT2013090206），河北省自然基金项目（C2015204176）

张乐超（1991-），男，河北省石家庄人，在读硕士，研究方向为动物遗传育种与繁殖，E-mail:1596653383@qq.com

*通讯作者: 李兰会（1972-），河北枣强人，副教授，硕士生导师，研究方向为动物遗传育种，E-mail:lanhuili13@163.com