牟腾慧,龙广丽,刘 洋,林家栋,张福平
(贵州大学高原山地动物遗传育种与繁殖教育部重点实验室,贵州大学贵州省动物遗传育种与繁殖重点实验室,贵州大学动物科学学院,贵州贵阳 550025)
威宁鸡是贵州省西部威宁地区的一个中型蛋肉兼用地方品种,以耐粗饲、抗寒、肉质鲜嫩、蛋黄品质优良而著称。威宁鸡体躯结构匀称、体质壮实、骨骼粗壮,其公鸡腿高粗壮,母鸡腿长短适中、胫部间有距毛,喙、胫爪为黑色[1]。
家禽生长曲线拟合分析是研究家禽生长发育规律的重要方法[2]。近年来,国内对地方品种鸡如琅琊鸡、双莲鸡、边鸡等生长发育规律均有研究[3-7]。本文采用家禽常见的3 种非线性生长曲线模型(Logistic 模型、Von Вertalanffy 模 型和Gompertz 模型)[8-9]对威 宁 鸡0~16 周龄生长发育规律进行拟合分析,通过生长曲线的参数值来反映威宁鸡公母鸡在生长率、最大体重方面的差异性,对威宁鸡的遗传特点和规律进行深入了解,并根据其生长曲线拟合分析结果为公母鸡制定合理饲养方案,提高威宁鸡生产效率,同时也为威宁鸡的资源保护利用开发和产业发展提供一定的科学依据和数据参考。
1.1 实验动物 本研究中威宁鸡由贵州大学科研鸡场提供。选择刚出壳的威宁鸡雏鸡180 只(公母各半,60只备用)进行个体编号,网上平养,自然通风。实验前对育雏室、育成室进行消毒,实验期间育雏室鸡舍温度30~33℃、24 h 光照,育成室进行自然光照,采用常规免疫程序进行免疫。
实验动物按照实验日龄分为2 个阶段,即0~6 周龄和6~16 周龄。0~6 周龄营养水平为:代谢能11.66 MJ/kg,粗蛋白19.9%,钙1%,有效磷0.37%,蛋氨酸0.9%,赖氨酸0.44%;6~16 周龄营养水平为:代谢能12.05 MJ/kg,粗蛋白17.49%,钙0.88%,有效磷0.33%,蛋氨酸0.9%,赖氨酸0.34%。
1.2 体重测定 雏鸡出雏时称测初生重,之后于2、4、6、8、10、12、14、16 周龄,随机选择120 只实验鸡(公母各半)禁食12 h 后进行单独称重,并计算其平均体重。
1.3 统计方法
1.3.1 拟合曲线模型 采用3 种常用非线性生长模型Logistic 模 型、Von Вertalanffy 模 型 和Gompertz 模 型对初生重及各周龄体重进行拟合与分析,表1 为3 种生长曲线模型及特征参数。
1.3.2 统计分析 本文运用Excel建立数据库,采用SPSS 20.0软件中的非线性回归分析拟合Logistic、Gompertz 和Вertalanffy 3 种模型的生长曲线,拟合计算出A、В、k 最优估计值,建立生长模型,分别计算出拐点体重、拐点周龄,根据拟合度(R2)评价生长模型,并利用Excel 作生长曲线图。
表1 3种生长曲线模型及特征参数
应用SPSS22.0 软件中的Data Editor 界面输入周龄和体重数据,依次选择Analtze、Regression、Nonlinera Regression 程序,在Dependent 文本框中选择“Weight”为变量,并在Model Expression分别输入Logistic(Gompertz或 Вertalanffy)公式进行生长曲线拟合算出模型参数A、В、k 的最优估计值,建立生长模型,计算拐点体重、拐点周龄,根据R2评价生长模型。
2.1 威宁鸡实际生长情况 由表2 可知,威宁鸡公、母鸡的平均增重分别为1 600.71、1 206.39 g,平均周增重分别为100.04、75.40 g。其中0~10 周是威宁鸡增重上升期,10 周龄、8 周龄、6 周龄增重达到高峰,相较前2 周公鸡分别增重264.50、252.08、228.80 g,母鸡分别增重202.29、195.73、176.14 g。12~16 周体重增重开始呈现下降,16 周龄、14 周龄、12 周龄相较前两周,公鸡分别增重212.19、216.21、202.53 g,母鸡分别增重169.92、144.16、133.75 g。由上述结果可知,威宁鸡生长速度在8~10 周龄达到最快,其次是6~8 周龄,12 周龄后生长趋于缓慢。
2.2 生长曲线拟合分析
2.2.1 3种生长曲线模型的参数估计和拟合度 由表3可知,威宁鸡公鸡、母鸡的3 种模型R2均在0.99 以上,均能很好地模拟威宁鸡的生长曲线,公鸡R2分别达到0.996、0.999 和0.999,母 鸡R2分 别 为0.995、0.999和0.999。由Logistic 模型拟合的公、母鸡的拐点体重最大为884.792、657.693 g,其次为Gompertz 模型,威宁鸡公、母鸡的拐点体重分别为800.903、583.856 g,Von Вertalanffy 模型拐点体重最小,威宁鸡公、母鸡的拐点体重分别为778.856、555.275 g。由Logistic 模型拟合的公、母鸡拐点周龄最大,分别为9.251、8.861 周,其次仍是Gompertz 模型,分别为8.564、7.991 周,最小的为Von Вertalanffy 模型,分别为8.365、7.648 周。
表2 威宁鸡公母鸡实际生长情况
表3 3种非线性生长曲线模型拟合结果
2.2.2 3种模型生长曲线拟合值与实测值的比较 由图1、2 可知,3 种生长曲线模型均能很好地模拟威宁公鸡、母鸡的生长曲线,并且其模型拟合的生长曲线与公鸡、母鸡生长发育实测值曲线基本吻合。由表4 可知,3 种模型预测的0 周龄公鸡、母鸡体重均与实际体重有一定的差异,其中Logistic 模型预测的0 周龄体重与实际相差最大,公鸡、母鸡预测值分别为83.396、69.923 g。Gompertz 模型预测的0 周龄公鸡、母鸡体重分别为39.715、35.575 g,与实测值差异最小。综合来看,Gompertz 模型拟合效果最好。
表4 威宁鸡周龄体重实测值与3 种模型预测值比较
本实验中,威宁鸡公母鸡前4 周体重差异不大,生长曲线基本相同,4 周后公鸡、母鸡体重出现极大分化,在6~10 周龄生长速度均加快,公鸡体重增重明显高于母鸡直至实验结束,这可能是由于营养需要不一所致。其中公母鸡均在10 周龄增重达到最大值,分别为132.25、101.15 g。在整个饲养过程中,应根据其公母鸡后期营养需求对饲养方案进行调整,特别是在威宁鸡6 周龄之后。其次,在6 周龄后由于威宁鸡公母鸡饲养方案不同,可将公母鸡分群饲养,这样既可以提升公母鸡生产效率,也能避免饲养过程中饲料浪费。
威宁鸡在0~16 周的曲线拟合结果以Gompertz 模型和Von Вertalanffy 模型的模拟效果更好,R2均达0.999,虽两者预测值与实测值都较为接近,但0~16 周龄预测值中初生重Gompertz 模型预测值与实测值更为接近,因此Gompertz 模型能够很好地拟合威宁鸡的生长曲线。不同品种的鸡种由于其自身品种、生长饲养环境等因素差异,可根据各生长曲线模型特点,选择最符合其生长发育规律的模型。如朱志明等[10]、曹爱青等[11]、吴艳等[5]研究认为藏鸡、朗德鹅、双莲鸡生长曲线为Gompertz模型拟合最好,姜自琴等[12]、李俊营等[13]认为Von Вertalanffy 模型能更好地拟合武定鸡、淮南麻黄鸡的生长曲线,而阳光远等[14]、李洪林等[15]则认为Logistic 模型能更好地拟合白宜黑鸡、贵州黄鸡慢羽系的生长曲线。
本实验发现Gompertz 模型拟合效果最好,与威宁鸡的实际生长测量最为接近,其R2在0.999。通过Gompertz 模型得到公鸡拐点体重为800.903 g,拐点周龄为8.564 周,母鸡拐点体重为583.856 g,拐点周龄为7.991 周。