用随机回归模型估计宁夏地区荷斯坦牛头胎测定日产奶量遗传参数

【目的】通过拟合泌乳曲线对影响宁夏地区荷斯坦牛测定日产奶量的诸多因素进行分析,并采用随机回归测定日模型估计遗传参数,为后续估计乳成分和体细胞数遗传参数及个体育种值估计奠定基础,同时也为制定适合该群体的优化育种方案提供参数。【方法】收集了2009—2013年间,宁夏地区24个牛场38 592头荷斯坦牛原始DHI数据,共包含550 078条测定日产奶量记录。质控标准为:产犊月龄为22—36 mo、泌乳天数为5—305 d及日产奶量为5.9—53.3 kg。质控后,共得到14 320头母牛的127 478条测定日记录(占原始记录的23.14%)。将有表型的母牛至少向上追溯三代(父母、祖父母、外祖父母),得到遗传评估所用的系谱文件,共包含24 272头奶牛。使用Excel 2013计算各测定日的平均产奶量,利用SAS 9.1软件的非线性回归过程(NLIN)结合Wood模型拟合群体水平的泌乳曲线。采用随机回归测定日动物模型,通过DMU 5.2软件对头胎产奶量进行遗传参数估计。随机回归测定日模型中包括了一般固定效应、固定回归项及随机回归项三类效应。研究中使用场-测定日作为固定效应、产犊年季效应作为固定回归项,并将加性遗传效应和永久环境效应作为随机回归项。测定日子模型采用4阶勒项德多项式(Legendre polynomials)。根据宁夏地区气候特点,将产犊季节划分为春(3月11日至5月20日)、夏(5月21日至8月25日)、秋(8月26日至10月15日)和冬(10月16日至3月10日)四季。【结果】宁夏地区荷斯坦牛头胎测定日产奶量群体均值为29.66 kg,在泌乳90 d时达到泌乳高峰,高峰产奶量为31.84 kg。通过不同产犊年份、产犊季节和牧场效应组合下的泌乳曲线拟合效果分析,以上因素通过影响泌乳曲线拟合效果进而影响头胎产奶量。在使用测定日模型估计头胎遗传参数时均应通过子模型拟合考虑其泌乳曲线的差异,在产奶量遗传评估模型中予以考虑。产奶量性状5—305d的遗传力在0.08—0.29之间,305d产奶量的总遗传力为0.16。【结论】通过使用wood模型拟合宁夏地区荷斯坦牛头胎泌乳曲线,优化了适合于该群体的测定日模型。利用修改后的测定日模型,估计了宁夏地区荷斯坦牛测定日产奶量的遗传参数。305 d遗传力(0.16)略低于其他群体的同类文研究结果。在测定日模型评估生产性能记录时,应注意对残差方差的处理并关注构建分子亲缘关系矩阵所需系谱的完整性,以保证遗传评估的可靠性。本研究为日后使用多性状随机回归测定日模型估计泌乳性状遗传参数研究奠定了基础。

江苏省荷斯坦牛日产奶量Wood模型DHI分析

利用Wood泌乳曲线模型对江苏省12个奶牛场2017~2019年荷斯坦牛398449条DHI日产奶量记录数据作曲线拟合,分析牧场规模、胎次、产犊季节、产犊间隔和305 d产奶量等因素对泌乳曲线相关参数影响。测定日产奶量泌乳曲线为标准泌乳曲线,高峰产奶量为39.08 kg,高峰产奶日为88.20 d,泌乳持久力为2.85,平均拟合度为0.9920;规模在1000头以下牧场荷斯坦牛泌乳曲线拟合度最大(0.9939),规模为1000~2000头牧场拟合度最小(0.9916);1~5胎奶牛泌乳曲线拟合度随胎次增加逐渐降低;冬季产犊拟合度最高(0.9923),夏秋两季产犊拟合度最低(0.9919);产犊间隔为366~400 d时拟合度最大(0.9923),300~365 d时拟合度最小(0.9920);305 d产奶量在11001~13000 kg拟合度最高(0.9955),305 d产奶量在3000~5000 kg拟合度最低(0.9864)。Wood模型适合江苏省荷斯坦牛日产奶量变化曲线拟合分析,牧场规模、胎次、产犊季节、产犊间隔、305 d产奶量等因素对Wood泌乳曲线模型参数均有影响。