吉林矮小芦花鸡早期生长曲线拟合及相关分析

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

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

吉林地方品种以芦花鸡较为常见，通过多世代选育，培育出个体小、基础代谢低、饲料报酬高、抗病抗应激能力强的肉蛋兼用型品系吉林矮小芦花鸡，具有良好的市场前景。

家禽生长发育规律呈非线性，通过非线性模型进行拟合分析，可较为动态了解其生长发育趋势与规律，指导家禽的选育和饲养管理，使其生产潜力发挥到最大[1]。目前最为常用的是Logistic、Gompertz 和 Bertalanffy 3种非线性曲线模型，已经广泛用于家禽生长发育规律研究中[2-3]。本试验通过体重和胫长指标，并结合3种曲线模型对吉林矮小芦花鸡生长发育规律进行拟合分析，旨在为其早期的选育提供科学依据。

1 材料与方法

1.1 试验材料

选择吉林省农业科学院畜牧分院试验鸡场的吉林芦花鸡在人工选择的基础上，通过纯繁获得吉林矮小芦花鸡(dw)和吉林芦花鸡(Dw)各200只，公母各100只，共400只。

1.2 饲养管理

确定分组后将每只试验鸡穿翅号，自由采食和饮水，日常免疫接种与管理。饲养过程分为1～5周的育雏期和6～13周的育成期。

1.3 测定指标

1.3.1 体重

用电子天平空腹称量试验鸡的初生重以及0、1、3、5、7、9、11、13周龄的体重，并做好记录。

1.3.2 胫长和胫骨长

用游标卡尺测量试验鸡7、9、11、13周龄的胫长和胫骨长(在测量胫骨长时，若长度超出游标卡尺最大范围，可用皮尺替换)，并做好记录。

1.4 生长曲线拟合

选择Logistic、Gompertz 和 Bertalanffy 3种非线性曲线模型对0～13周生长曲线进行拟合，具体表达式与参数见表1。

表1 3种非线性模型曲线

1.5 数据统计与分析

利用Excel进行初步数据整理，结合SPSS 24.0单因素方差分析和Duncan’s多重比较，数值用“平均数±标准误”表示，拟合3种非线性曲线，分别计算A、B、K的最优估计值，建模后再计算出拐点体重、周龄以及最大周增重。然后对体重与胫长和胫骨长进行相关性分析。

2 结果与分析

2.1 体重实际生长情况

由表2和图1可知，dw母鸡初生重与dw公鸡、Dw母鸡差异不显著(P>0.05)，与Dw公鸡差异显著(P<0.05)；dw母鸡从3周开始生长速度明显加快，但与另外3组相比较体重保持最低，差异显著(P<0.05)，在第7周时达到高峰，随后生长速度放缓，与Dw母鸡相比，从5周开始，生长趋势大致保持一致。dw公鸡5～9周增长迅速，并在9周时达到高峰，比Dw公鸡高峰期提前了2周左右，随后放缓在11周后又略有上升。在整个早期生长阶段，dw公鸡始终低于Dw公鸡，差异显著(P<0.05)。早期同一生长阶段体重大致为：Dw公鸡>Dw母鸡≈dw公鸡>dw母鸡。

表2 0～13周龄体重测定值比较 g

图1 0～13周龄吉林矮小芦花鸡周增重

2.2 生长曲线拟合分析

由表3可知，所有鸡种生长曲线模型拟合度均超过了0.99，其中dw母鸡Gompertz和Bertalanffy模型拟合度达到0.995，Gompertz模型拐点体重为607.924 g，拐点周龄为7.345，最大周增重达到113.682 g；dw公鸡Logistic和Gompertz模型拟合度相同；Dw母鸡Gompertz模型拟合度最高，达到0.995，拐点体重为808.599 g，拐点周龄7.342，最大周增重为164.146 g；Dw公鸡Logistic和Gompertz模型拟合度均为0.995，拐点体重、周龄、最大周增重差异不明显，均可作为适合模型，另外，Dw、dw公鸡的拐点周龄均晚于Dw、dw母鸡。

表3 3种生长曲线模型的拟合参数

2.3 生长曲线与3种拟合曲线的比较

如表4、表5所示，通过比较dw和Dw鸡实际测定值和3种生长曲线模型的估计值发现，Logistic模型的初生重估计值远高于实际测定值，Bertalanffy模型的初生重估计值远低于实际测定值，而Gompertz模型的初生重估计值与实际测定值较为接近，整个早期生长阶段的模拟值与实际测定值吻合度较高，因此Gompertz模型拟合效果优于Logistic和Bertalanffy模型。

表4 dw鸡实际测定值与3种拟合曲线的比较

表5 Dw鸡实际测定值与3种拟合曲线的比较

2.4 胫长和胫骨长的发育情况

各组鸡胫长和胫骨长发育情况见表6、表7。同一阶段dw母鸡胫长均显著低于Dw母鸡(P<0.05)，随着周龄增加，dw母鸡胫长分别占Dw母鸡的82.7%、82.4%、75.5%和72.2%，差距逐渐变大，dw母鸡在7～9周胫长增速明显，9周之后增长速度放缓，Dw母鸡在7～11周增速逐渐加快，11周以后增速放缓；同一阶段dw母鸡胫骨长均显著低于Dw母鸡(P<0.05)，在观测期中，dw母鸡胫骨长分别占Dw母鸡的85.2%、80.0%、81.1%和68.4%，差距逐渐拉大，dw母鸡胫骨长在7～11周平稳增加，11周以后趋于平稳，而Dw母鸡胫骨长在观测期一直逐渐增加；dw母鸡胫长与胫骨长生长趋势几乎一致，Dw母鸡胫长与胫骨长生长趋势大致相同。

表6 7～13周龄胫长比较 mm

表7 7～13周龄胫骨长比较 mm

同一阶段dw与Dw公鸡胫长相比较差异显著(P<0.05)，dw公鸡胫长先增加在9周后逐渐稳定，dw公鸡胫长分别占Dw公鸡的86.3%、84.9%、73.9%和73.7%，而Dw公鸡胫长呈现7～11周逐渐增加，11周后增速放缓，胫长趋于稳定，dw公鸡胫骨长分别占Dw公鸡的85.4%、81.9%、74.9%和73.0%；同一阶段dw与Dw公鸡胫骨长相比较差异显著(P<0.05)，Dw公鸡胫骨长逐渐增加，dw公鸡胫骨的增长速度呈先增加后减少的趋势，dw公鸡胫长与胫骨长生长趋势无明显差别，Dw公鸡胫长与胫骨长生长趋势大致相同。

2.5 体重与胫长、胫骨长的相关性分析

体重与胫长和胫骨长的结果见表8和表9。不同周龄各组合的胫长与体重呈现出极显著正相关(P<0.01)，其中dw与Dw母鸡分别在7周龄和9周龄相关系数最大，为0.934和0.891；dw和Dw公鸡均在9周龄相关系数最大；dw鸡胫长与体重的相关程度高于正常的Dw鸡。不同周龄各组合的胫骨长与体重除dw母鸡13周龄呈现出显著正相关(P<0.05)外，其余均呈现出极显著正相关(P<0.01)，其中dw和Dw母鸡胫骨长与体重的相关程度分别在7、9周龄时最大，而dw和Dw公鸡则表现在9和11周龄。

表8 不同周龄胫长与体重相关性分析

表9 不同周龄胫骨长与体重相关性分析

胫长与胫骨长的相关性结果如表10所示，7～9周龄各组合胫长与胫骨长的相关性达到极显著水平(P<0.01)；11周龄dw母鸡和公鸡胫长与胫骨长的相关性达到极显著水平(P<0.01)，Dw母鸡和Dw公鸡胫长与胫骨长的相关性达到显著水平(P<0.05)；13周龄各组合胫长与胫骨长的相关性达到显著水平(P<0.05)。

表10 不同周龄胫长与胫骨长相关性分析

3 讨论

3.1 生长曲线拟合特点

动物生长发育具有相对稳定的基本特征,这是用数学模型描述动物体生长过程的基础[4]。非线性生长模型是根据生物体生物量增加，描述其随时间规律性的变化的过程[5]。动物生长发育在遗传、环境及营养等多种因素的作用下，不同品种间差异大，因此不同品种要用不同的生长模型去拟合。目前常用的3种非线性生长模型是Logistic、Gompertz和Bertalanffy模型，Logistic模型适用于早期生长速度慢的品种，Gompertz则适用于早期生长速度快的品种，Bertalanffy适用于生长较为迟缓的动物[6]。本研究利用这3个非线性曲线模型对0～13周的吉林矮小芦花鸡和正常芦花鸡进行了拟合分析。结果表明Gompertz模型拟合度(R2)最高达到0.995，拟合值与实际测定值最为接近，效果最佳，这一结果与Mohammed[7]、李青青等[8]、杜德英等[9]、付睿琦等[10]对罗斯308、独龙鸡、文昌鸡和北京油鸡的研究结果相一致。Manjula等[11]和包艳波等[12]则认为Logistic模型更适用于韩国本地鸡(KNC)和大骨鸡的生长曲线拟合。姜自琴等[13]认为，Bertallanffy 模型在拟合度和预测武定鸡体重方面效果最好。吉林矮小芦花鸡和正常吉林芦花鸡中，公鸡的生长周期始终要比母鸡长，拐点周龄相对较晚,所以在饲养时要注意分群，并在此之前提供充足的养分，来满足其生长发育的需求。因此，Gompertz模型对吉林矮小芦花鸡的拟合效果最佳，拐点周龄及体重与实测值最为接近。

3.2 胫长、胫骨长和体重的相关性分析

骨骼作为家禽机体的重要组成，在其生长发育过程中有着重要的意义，对于生产潜力挖掘产生了重要影响[14]。胫长和胫骨长则是衡量家禽早期生长发育的重要指标，与体重结合进行相关性分析，可为吉林矮小型芦花鸡的选育提供有价值的参考。本研究结果表明，dw和Dw公母鸡在各自周龄胫长、胫骨长与体重呈极显著正相关，这与陈锐等[15]认为略阳乌鸡胫长与体质量呈极显著正相关的结果相一致。叶伟庆等[16]的研究结果表明，信宜怀乡鸡体质量与胫长极显著正相关。成海荣等[17]提出，胫骨长度和重量是反映肉鸡骨骼发育的重要指标。本研究发现，dw和Dw公母鸡胫骨相对生长在9～11周龄时达到最高，这与宁康健等[18]发现青脚麻公母鸡胫骨的相对生长均在28日龄达到最高结果不一致。这可能是由于品种、营养、环境等因素造成的胫骨发育的差异。吉林矮小芦花鸡各周龄体重与胫长和胫骨长呈极显著正相关，胫长和胫骨长可作为反映生长发育的主要指标。