基于系谱和基因组信息估计荷斯坦青年母牛体重性状遗传参数

doi:10.11843/j.issn.0366-6964.2022.11.005

摘要/Abstract

摘要： 旨在设计利用不同信息来源的模型估计荷斯坦后备牛不同月龄体重性状的遗传参数。本研究于2014—2020年测定并收集了7 122头荷斯坦牛32 338条0~12月龄体重数据，分别利用系谱信息（linear mixed model with pedigree relationship matrix，LM_A）和系谱-基因组信息构建亲缘关系矩阵（linear mixed model with genotype-pedigree joint relationship matrix，LM_H），基于母体效应动物模型估计初生重，基于是否考虑初生重作为协变量的单性状动物模型估计2~12月龄各月龄体重遗传力，并利用双性状动物模型估计初生重与其它月龄体重的遗传相关。结果显示，对于初生重，根据赤池信息量准则（Akaike information criterion，AIC），LM_H方法的拟合程度显著优于LM_A方法，但两种方法估计的遗传参数相差不大：直接遗传力分别为0.30和0.32，母体遗传力分别为0.08和0.09，个体直接遗传效应和母体遗传效应遗传相关系数分别为-0.65和-0.64；对于2~12月龄体重，LM_A和LM_H两种方法估计的校正初生重后的各月龄体重遗传力分别为0.15~0.55和0.28~0.49，未校正初生重的各月龄体重遗传力分别为0.16~0.54和0.28~0.51。初生重与2、5月龄体重之间为高遗传相关（相关系数>0.6）。5月龄后，各月龄体重与初生重的遗传相关系数随着时间间隔的增加而减小。相较于LM_A，LM_H方法更稳定，AIC值较小（即拟合优度较大），遗传参数标准误较小。综上，采用LM_H方法估计目标性状可获得更准确、更稳定的遗传参数。本研究为建立中国荷斯坦牛生长性状基因组选择体系提供了理论依据。

关键词: 荷斯坦牛, 体重, 遗传参数, 系谱, 基因组信息

Abstract: The aim of this study was to estimate genetic parameters of body weight (BW) at different ages in Holstein heifers using models considering different source of information. A data set including 32 338 BW records of 0-12 months old on 7 122 female Holstein cattle measured from 2019 to 2020 was collected. Variance components and parameters were estimated using linear mixed model with pedigree relationship matrix (LM_A) and genotype-pedigree joint relationship matrix (LM_H). For calf birth weight (CBW), animal models with maternal genetic effects were applied. The single-trait animal models with or without a CBW as a covariate were used for the estimation of heritability of monthly body weight from month 2 to 12. Genetic correlations between CAW and monthly BW were estimated using bivariate animal models. The results showed that models with LM_H had better goodness of fit based on Akaike information criterion (AIC), although similar parameters for CBW were obtained through LM_A and LM_H methods. For CBW, direct heritabilities estimated by LM_A and LM_H were 0.30 and 0.32, respectively, and maternal heritabilities were 0.08 and 0.09, respectively. The correlation coefficients between individual direct genetic effect and maternal genetic effect was antagonistic for LM_A (-0.65) and LM_H (-0.64). Monthly BW traits had moderate to high heritabilities. Estimates from LM_A and LM_H with model including CBW as a covariate ranged from 0.15 to 0.55 and 0.28 to 0.49, respectively, and heritabilities from LM_A and LM_H with model without CBW ranged from 0.16 to 0.54 and 0.28 to 0.51, respectively. BW at age of month 2 and 5 were highly genetically correlated with CBW (correlation coefficient was higher than 0.6). After 5 months old, genetic correlations between CBW and BW at different months were decreased with increasing time spans. Compared with LM_A, the LM_H method was more stable, with a smaller AIC value (that is, a larger goodness of fit) and a smaller standard error of the genetic parameters. In conclusion, using the LM_H method to estimate target traits can obtain more accurate and stable genetic parameters. This study provides a theoretical basis for the establishment of a genomic selection system for growth traits in Chinese Holstein cattle.

Key words: Holstein cattle, body weight, genetic parameters, pedigree, genomic information

中图分类号:

S823.2

常瑶, 苏国生, 李艳华, 李想, 麻柱, 王雅春. 基于系谱和基因组信息估计荷斯坦青年母牛体重性状遗传参数[J]. 畜牧兽医学报, 2022, 53(11): 3759-3768.

CHANG Yao, SU Guosheng, LI Yanhua, LI Xiang, MA Zhu, WANG Yachun. Estimating Genetic Parameters for Body Weights using Pedigree and Genotype-pedigree based Approaches in Holstein Heifers[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(11): 3759-3768.

参考文献 35

[1]	YIN T, KÖNIG S.Genomic predictions of growth curves in Holstein dairy cattle based on parameter estimates from nonlinear models combined with different kernel functions[J].J Dairy Sci, 2020, 103(8):7222-7237.
[2]	BYRNE T J, SANTOS B F S, AMER P R, et al.New breeding objectives and selection indices for the Australian dairy industry[J].J Dairy Sci, 2016, 99(10):8146-8167.
[3]	YIN T, KONIG S.Genetic parameters for body weight from birth to calving and associations between weights with test-day, health, and female fertility traits[J].J Dairy Sci, 2018, 101(3):2158-2170.
[4]	PRYCE J E, GONZALEZ-RECIO O, NIEUWHOF G, et al.Hot topic:definition and implementation of a breeding value for feed efficiency in dairy cows[J].J Dairy Sci, 2015, 98(10):7340-7350.
[5]	SPURLOCK D M, DEKKERS J C M, FERNANDO R, et al.Genetic parameters for energy balance, feed efficiency, and related traits in Holstein cattle[J].J Dairy Sci, 2012, 95(9):5393-5402.
[6]	HEINRICHS A J, HARGROVE G L.Standards of weight and height for Holstein heifers[J].J Dairy Sci, 1987, 70(3):653-660.
[7]	HEINRICHS A J, LOSINGER W C.Growth of Holstein dairy heifers in the United States[J].J Anim Sci, 1998, 76(5):1254-1260.
[8]	ABDALLAH J M, MCDANIEL B T.Genetic parameters and trends of milk, fat, days Open, and body weight after calving in North Carolina experimental herds[J].J Dairy Sci, 2000, 83(6):1364-1370.
[9]	VAN ELZAKKER P J M, VAN ARENDONK J A M.Feed intake, body weight and milk production:genetic analysis of different measurements in lactating dairy heifers[J].Livest Prod Sci, 1993, 37(1):37-51.
[10]	DEMEKE S, NESER F W C, SCHOEMAN S J.Variance components and genetic parameters for early growth traits in a mixed population of purebred Bos indicus and crossbred cattle[J].Livest Prod Sci, 2003, 84(1):11-21.
[11]	COFFEY M P, HICKEY J, BROTHERSTONE S.Genetic aspects of growth of Holstein-Friesian dairy cows from birth to maturity[J].J Dairy Sci, 2006, 89(1):322-329.
[12]	KOENEN E P C, GROEN A F.Genetic analysis of growth patterns of black and white dairy heifers[J].J Dairy Sci, 1996, 79(3):495-501.
[13]	BROTHERSTONE S, COFFEY M P, BANOS G.Genetic parameters of growth in dairy cattle and associations between growth and health traits[J].J Dairy Sci, 2007, 90(1):444-450.
[14]	BANOS G, WIGGANS G R, POWELL R L.Impact of paternity errors in cow identification on genetic evaluations and international comparisons[J].J Dairy Sci, 2001, 84(11):2523-2529.
[15]	VISSCHER P M, WOOLLIAMS J A, SMITH D, et al.Estimation of pedigree errors in the UK dairy population using microsatellite markers and the impact on selection[J].J Dairy Sci, 2002, 85(9):2368-2375.
[16]	周磊, 刘林, 初芹, 等.奶牛亲子鉴定应用的标记和方法研究进展[J].中国奶牛, 2011(2):26-29.ZHOU L, LIU L, CHU Q, et al.Progress on label and method of paternity test in dairy cattle[J].China Dairy Cattle, 2011(2):26-29.(in Chinese)
[17]	WELLER J I, FELDMESSER E, GOLIK M, et al.Factors affecting incorrect paternity assignment in the Israeli Holstein population[J].J Dairy Sci, 2004, 87(8):2627-2640.
[18]	杨湛澄, 闫青霞, 黄河天, 等.中国荷斯坦公牛系谱完整性研究[J].中国奶牛, 2017(2):22-25.YANG Z C, YAN Q X, HUANG H T, et al.Evalutation of pedigree completeness in Chinese Holstein bulls[J].China Dairy Cattle, 2017(2):22-25.(in Chinese)
[19]	MISZTAL I, AGUILAR I, LEGARRA A, et al.Choice of parameters for single-step genomic evaluation for type[J].J Dairy Sci, 2010, 93(S1):166.
[20]	MOTULSKY H J, BROWN R E.Detecting outliers when fitting data with nonlinear regression-a new method based on robust nonlinear regression and the false discovery rate[J].BMC Bioinformatics, 2006, 7:123.
[21]	BROWNING B L, BROWNING S R.A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals[J].Am J Hum Genet, 2009, 84(2):210-223.
[22]	MISZTAL I, TSURUTA S, LOURENCO D, et al.Manual for BLUPF90 family of programs[EB/OL].(2014-05-12).http://nce.ads.uga.edu/wiki/lib/exe/fetch.php?media=blupf90_all2.pdf.
[23]	BURNHAM K P, ANDERSON D R.Multimodel inference:understanding AIC and BIC in model selection[J].Sociol Method Res, 2004, 33(2):261-304.
[24]	于志, 高艳霞, 曹玉凤, 等.中国荷斯坦犊牛和育成牛生长发育规律的研究[J].中国畜牧兽医, 2014, 41(6):121-125.YU Z, GAO Y X, CAO Y F, et al.Study on growth and development pattern of Chinese Holstein calf and heifer[J].China Animal Husbandry & Veterinary Medicine, 2014, 41(6):121-125.(in Chinese)
[25]	JOHANSON J M, BERGER P J, TSURUTA S, et al.A Bayesian threshold-linear model evaluation of perinatal mortality, dystocia, birth weight, and gestation length in a Holstein herd[J].J Dairy Sci, 2011, 94(1):450-460.
[26]	EVERETT R W, MAGEE W T.Maternal ability and genetic ability of birth weight and gestation length[J].J Dairy Sci, 1965, 48(7):957-961.
[27]	CLÉMENT V, BIBÉ B, VERRIER É, et al.Simulation analysis to test the influence of model adequacy and data structure on the estimation of genetic parameters for traits with direct and maternal effects[J].Genet Sel Evol, 2001, 33(4):369-395.
[28]	HEYDARPOUR M, SCHAEFFER L R, YAZDI M H.Influence of population structure on estimates of direct and maternal parameters[J].J Anim Breed Genet, 2008, 125(2):89-99.
[29]	BALDI F, ALBUQUERQUE L G, ALENCAR M M.Random regression models on Legendre polynomials to estimate genetic parameters for weights from birth to adult age in Canchim cattle[J].J Anim Breed Genet, 2010, 127(4):289-299.
[30]	周身娉, 付帝生, 杨明, 等.仔猪初生重对育肥期生长性状及采食性状的影响[J].中国畜牧杂志, 2019, 55(1):41-45.ZHOU S P, FU D S, YANG M, et al.Effect of Duroc birth weight on their growth and feeding traits in fattening stage[J].Chinese Journal of Animal Science, 2019, 55(1):41-45.(in Chinese)
[31]	WATANABE R N, DO NASCIMENTO G B, BERNARDES P A, et al.Inclusion of genomic information in estimation of genetic parameters for body weights and visual scores in Nelore cattle[J].Rev Bras Zootec, 2021, 50:e20200077.
[32]	LOBERG A, DVRR J W, FIKSE W F, et al.Estimates of genetic variance and variance of predicted genetic merits using pedigree or genomic relationship matrices in six Brown Swiss cattle populations for different traits[J].J Anim Breed Genet, 2015, 132(5):376-385.
[33]	SCHROOTEN I C.Genomic selection in dairy cattle[J].Wur Wageningen Ur, 2011, 95(4):2097-2107.
[34]	VEERKAMP R F, MULDER H A, THOMPSON R, et al.Genomic and pedigree-based genetic parameters for scarcely recorded traits when some animals are genotyped[J].J Dairy Sci, 2011, 94(8):4189-4197.
[35]	GROEN A F, VOS H.Genetic parameters for body weight and growth in Dutch Black and White replacement stock[J].Livest Prod Sci, 1995, 41(3):201-206.