Printer Friendly

Genomic heritability of bovine growth using a mixed model.


Genetic improvement of bovine growth is one of the most important tasks for a higher productivity in the beef cattle industry, and thus its heritability as the portion of genetic variability to phenotypic variability has been a great concern to cattle geneticists. Heritability for growth traits of beef cattle has been estimated in a variety of breeds and populations. Traditionally, they were estimated based on pedigree information. This heritability estimated by the traditional method is hereafter called "classical heritability".

Nowadays high throughput genotyping has been available. The use of commercial DNA microarray chip is currently a routine in identifying genomewide association signals for bovine complex phenotypes (Van Tassell et al., 2008; Snelling et al., 2010; Kim et al., 2011; Lu et al., 2013). Abundant information on tens of thousands of markers from the chip experiment actually enables us not only to test their genetic associations with phenotypes but also to assess genetic variability among individuals. Incorporating the genetic variability into the genetic association analysis in a mixed model would control population stratification and result in heritability estimate explained by polygenic effects (Yang et al., 2011). The objective of this study was to examine heritability estimates for growth traits of Korean cattle using genomic information of single nucleotide polymorphism (SNP) genotypes obtained with a DNA microarray chip.


Animals and phenotypes

An indigenous beef cattle breed called Korean cattle (also named Hanwoo or Bos taurus coreanae) was used to estimate heritability of body weight. For characteristics and historical details of the breed, see Kim and Lee (2000) and Lee and Pollak (2002). Body weights of 367 steers under a controlled feeding program were collected from Korean Cattle Improvement Center, National Agricultural Cooperative Federation in Korea. The feeding program was described in Kim et al. (2011). Body weights were available at 6, 9, 12, 15, 18, 21, and 24 months. The descriptive statistics for body weights are summarized in Table 1. Their meat samples were obtained at slaughter in spring and fall of 2007 through 2009.


Genomic DNA was extracted from the tissue samples of animals using a commercial kit from Qiagen (QIAGEN Inc., Valencia, CA, USA). The bovine Illumina Infinium assay (BovineSNP50, Illumina, San Diego, CA, USA) with 54,001 SNPs across the whole bovine genome was used for genotyping. Genotypes were determined based on custom clusters developed with more than 7,000 samples using Beadstudio software provided by Illumina. Every signal cluster was manually reviewed and appropriately fitted.

We excluded the samples with call rate <98% from the analysis of estimating heritability. Exclusion criteria for SNPs were call rate <98%, minor allele frequency <0.05, and deviation from Hardy-Weinberg equilibrium. The Hardy-Weinberg disequilibrium was tested using the chi-square statistic with 2 degrees of freedom (p < 0.001).

Statistical analysis

A mixed model was employed to estimate genetic and environmental variance components as following:

y = X[beta] + g + e

where y is a vector of body weights, [beta] is a vector of fixed effects for contemporary group (year and season, 6 levels) effects, X is the corresponding incidence matrix, g is a vector of random polygenic effects under the assumption of N(0, G[[sigma].sup.2.sub.g]) and e is a vector of random environmental effects under the assumption of N(0,|[[sigma].sup.2.sub.e]). G denotes a genetic relationship matrix with elements of pairwise relationship coefficient (VanRaden, 2008; Yang et al., 2011), and the pairwise relationship was estimated by the following formula:

[[psi].sup.-1] [[psi].summation over (i=1)] ([n.sub.ij] - [2p.sub.i]) ([n.sub.ik] - [2p.sub.i])/[2p.sub.i](1-[p.sub.i])

where [psi] is the number of SNPs, [n.sub.ij] ([n.sub.ik]) is the number (0, 1, or 2) of the reference allele for the ith SNP of the jth (kth) individuals, and [p.sub.i] is the frequency of the reference allele of the ith SNP. The [[sigma].sup.2.sub.g] and [[sigma].sup.2.sub.e] are polygenic and environmental variance components, and | is an identity matrix.

The genetic and environmental variance components were estimated using restricted maximum likelihood (REML). Genome-wide Complex Trait Analysis (GCTA, version 2.0) which is a freely available program at was used to estimate the variance components (Yang et al., 2011). Heritability was calculated using the REML estimates as [[sigma].sup.2]/[[sigma].sup.2.sub.g] + [[sigma].sup.2.sub.g] + [[sigma].sup.2.sub.e], and its variance was estimated using the negative inverse of the average information matrix.

Further analysis was conducted to examine chromosomal contribution to heritability of body weight. We first calculated genetic relationship coefficients using genomic information by chromosomes and then to estimate heritability with the genetic relationship matrix for each chromosome.

We also conducted a genetic association analysis to estimate regression coefficients in testing associations of SNPs with body weights. Correlation of the regression coefficients between body weights at different age was estimated.


Heritability was estimated using 364 samples and 39,112 SNPs remained after the quality assurance filtering. The number of SNPs used in the analysis indicated that reasonably informative SNPs were included in the BovineSNP50 BeadChip for Korean cattle as reported before (Kim et al., 2011). Genetic association of body weight with any individual SNP was not identified by the Bonferroni multiple test (p > 1.83 x [10.sup.-7]).

Heritability estimates for body weights of Korean cattle ranged from 0.58 to 0.76 (Figure 1). They showed a trend that the heritability was larger at a smaller age. For example, the heritability estimate was largest for body weight at age of 6 months and smallest for body weight at age of 2 years.

The chromosomal contribution to heritability of body weight was heterogeneous by age in months (Figure 2). The contribution of chromosome 24 to heritability decreased with a larger age whereas the contribution of chromosome 13 increased with a larger age.

Regression coefficients were obtained from testing genetic associations of SNPs with body weights. Their correlation estimates between body weights at different age ranged from 0.66 to 0.98 (Figure 3). They decreased with a large difference in age.



The current study showed heritability estimates of 0.58 to 0.76 for body weights of beef cattle using genomic information of SNP genotypes obtained with the BovineSNP50 BeadChip. They were larger than classical heritability estimated using pedigree information in previous studies. For example, the classical heritability estimates for body weights of beef cattle ranged from 0.32 to 0.52 for Angus cattle (Davis and Simmen, 2006), from 0.18 to 0.49 for Korean cattle (Kim et al., 2000), from 0.19 to 0.25 for Simmental cattle (Lee et al., 1997), and from 0.26 to 0.35 for the US Meat Animal Research Center III composite breed (Rolfe et al., 2011). We suspect that the larger heritability might be produced at least partially by the use of genomic relationships although it might be influenced by the heterogeneous breeds, selection programs, and environments.


Since the SNP markers of the DNA chip used in this study were sparsely located in the bovine genome, heritability estimated with an increasing number of SNP markers would have been probably even larger than those estimated in the current study as suggested by simulation studies (Zhang et al., 2011; Ryoo and Lee, 2014). Most markers significantly associated with body weights were expected to link to functional genetic variants, and moreover both functional variant itself and variants in strong linkage with it might have been excluded from the analysis.

Thus, the phenotypic variability for growth of beef cattle might be explained by genomic information better than by pedigree information. This is opposite from a general tendency for human heritability. That is, heritabilities for human complex traits have been underestimated with genomic information comparing to those with pedigree information (Manolio et al., 2009). It has been called 'missing heritability', and this missing heritability might be attributed to the genetic relationships diluted by noncontributing markers (Ryoo and Lee, 2014) and to the inflated heritability from shared environmental factors (Zuk et al., 2012). In practice, heritability might be underestimated by diluted genetic relationships with genomic information or with pedigree information. Furthermore, the underestimation becomes more serious with wrong information, especially with wrong pedigree or misidentification (Lee and Pollak, 1997).


The current study showed a trend that the heritability estimates for body weight decreased at older ages. This suggests a necessity of an early genetic evaluation for the growth traits of animals because their genetic merits might be more explained at an earlier age, 6 months in the current study. Additional association studies of body weight at an early age with previously identified variants for growth of Korean cattle might be useful for identifying stronger signals. The genomic sequence variants identified for association with growth or growth-related traits of Korean cattle might include intragenic SNPs of the genes encoding dishevelled segment polarity protein 1 for eye muscle area (Kim et al., 2011), mitogen-activated protein kinase kinase 6 (MAP2K6) and uncoupling protein 2 for carcass weight (Ryu et al., 2012), and phosphodiesterase 1B (PDE1B) for Longissimus dorsi muscle area (Shin et al., 2012) and intergenic SNPs of rs110228023 for eye muscle area (Kim et al., 2011) and rs29012331 for yearling weight (Kim et al., 2014). Furthermore, associations of variants identified for meat quality with an early body weight might be useful to determine whether their alleles are beneficial or conflicting to both quantity and quality of meat in breeding strategy. For example, sequence variants of human immunodeficiency virus type I enhancer binding protein 3 (Kim et al., 2011) and MAP2K6 (carcass weight) (Ryu et al., 2012) were associated with marbling score of Korean cattle, and variants of thyroid hormone responsive (La et al., 2013), peroxisome proliferator-activated receptor gamma (Oh et al., 2012b), PDE1B (Shin et al., 2012), and fatty acid synthase (Oh et al., 2012a) were associated with fatty acid composition. We suspect their heterogeneous effects on body weights by age because the correlation of regression coefficients for genetic associations decreased with a large difference in age.

This study showed a novel heritability estimate of bovine growth using genomic information, which was larger than those traditionally estimated with pedigree information. Also, the heritability estimates by age revealed a trend that the heritability increased for body weight at an earlier age. An early genetic evaluation using genomic information might be suggested to improve genetic merits for growth of Korean cattle.


This study was supported by the Basic Science Research Program of the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science, and Technology (Grant No. 2012002096).


Davis, M. E. and R. C. M. Simmen. 2006. Genetic parameter estimates for serum insulin-like growth factor I concentrations, and body weight and weight gains in Angus beef cattle divergently selected for serum insulin-like growth factor I concentration. J. Anim. Sci. 84:2299-2308.

Kim, J. B. and C. Lee. 2000. Historical look at the genetic improvement in Korean cattle. Asian Australas. J. Anim. Sci. 13:1467-1481.

Kim, Y, J. Ryu, J. Woo, J. B. Kim, C. Y Kim, and C. Lee. 2011. Genome-wide association study reveals five nucleotide sequence variants for carcass traits in beef cattle. Anim. Genet. 42:361-365.

Kim, Y, J. Ryu, and C. Lee. 2014. Replicated association of single-nucleotide marker on chromosome 6 with bovine yearling weight using a mixed model analysis. Anim. Genet. 45:151-153.

La, B., D. Oh, Y. Lee, S. Shin, C. Lee, E. Chung, and J. Yeo. 2013. Association of bovine fatty acid composition with novel missense nucleotide polymorphism in the thyroid hormone-responsive (THRSP) gene. Anim. Genet. 44:118.

Lee, C. and E. J. Pollak. 1997. Influence of sire misidentification on sire x year interaction variance and direct-maternal genetic covariance for weaning weight in beef cattle. J. Anim. Sci. 75: 2858-2863.

Lee, C., C. P. Van Tassell, and E. J. Pollak. 1997. Estimation of genetic variance and covariance components for weaning weight in Simmental cattle. J. Anim. Sci. 75:325-330.

Lee, C. and E. J. Pollak. 2002. Genetic antagonism between body weight and milk production in beef cattle. J. Anim. Sci. 80:316-321.

Lu, D., S. Miller, M. Sargolzaei, M. Kelly, G Vander Voort, T. Caldwell, Z. Wang, G Plastow, and S. Moore. 2013. Genome-wide association analyses for growth and feed efficiency traits in beef cattle. J. Anim. Sci. 91:3612-3633.

Manolio, T. A., F. S. Collins, N. J. Cox, D. B. Goldstein, L. A. Hindorff, D. J. Hunter, M. I. McCarthy, E. M. Ramos, L. R. Cardon, A. Chakravarti, J. H. Cho, A. E. Guttmacher, A. Kong, L. Kruglyak, E. Mardis, C. N. Rotimi, M. Slatkin, D. Valle, A. S. Whittemore, M. Boehnke, A. G. Clark, E. E. Eichler, G. Gibson, J. L. Haines, T. F. Mackay, S. A. McCarroll, and P. M. Visscher. 2009. Finding the missing heritability of complex diseases. Nature 461:747-753.

Oh, D., Y Lee, B. La, J. Yeo, E. Chung, Y. Kim, and C. Lee. 2012a. Fatty acid composition of beef is associated with exonic nucleotide variants of the gene encoding FASN. Mol. Biol. Rep. 39:4083-4090.

Oh, D., Y. Lee, C. Lee, E. Chung, and J. Yeo. 2012b. Association of bovine fatty acid composition with missense nucleotide polymorphism in exon7 of peroxisome proliferator-activated receptor gamma gene. Anim. Genet. 43:474.

Rolfe, K. M., W. M. Snelling, M. K. Nielsen, H. C. Freetly, C. L. Ferrell, and T G. Jenkins. 2011. Genetic and phenotypic parameter estimates for feed intake and other traits in growing beef cattle, and opportunities for selection. J. Anim. Sci. 89:3452-3429.

Ryoo, H. and C. Lee. 2014. Underestimation of heritability using a mixed model with a polygenic covariance structure in a genome-wide association study for complex traits. Eur. J. Hum. Genet. 22:851-854.

Ryu, J., Y. Kim, C. Kim, J. Kim, and C. Lee. 2012. Association of bovine carcass phenotypes with genes in an adaptive thermogenesis pathway. Mol. Biol. Rep. 39:1441-1445.

Shin, S., J. Heo, J. Yeo, C. Lee, and E. Chung. 2012. Genetic association of Phosphodiesterase 1B (PDE1B) with carcass traits in Korean cattle. Mol. Biol. Rep. 39:4869-4874.

Snelling, W. M., M. F. Allan, J. W. Keele, L. A. Kuehn, T. McDaneld, T. P. L. Smith, T. S. Sonstegard, R. M. Thallman, and G. L. Bennett. 2010. Genome-wide association study of growth in crossbred beef cattle. J. Anim. Sci. 88:837-848.

VanRaden, P. M. 2008. Efficient methods to compute genomic predictions. J. Dairy Sci. 91:4414-4423.

Van Tassell, C. P., T. P. L. Smith, L. K. Matukumalli, J. F. Taylor, R. D. Schnabel, C. T. Lawley, C. D. Haudenschild, S. S. Moore, W. C. Warren, and T. S. Sonstegard. 2008. SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries. Nat. Methods 5:247-252.

Yang, J., S. H. Lee, M. E. Goddard, and P. M. Visscher. 2011. GCTA: A tool for genome-wide complex trait analysis. The Am. J. Hum. Genet. 88:76-82.

Zuk, O., E. Hechter, S. R. Sunyaev, and E. S. Lander. 2012. The mystery of missing heritability: Genetic interactions create phantom heritability. Proc. Natl. Acad. Sci. 109:1193-1198.

Zhang, Z., X. Ding, J. Liu, Q. Zhang, and D. J. de Koning. 2011. Accuracy of genomic prediction using low-density marker panels. J. Dairy Sci. 94:3642-3650.

Jihye Ryu and Chaeyoung Lee *

Department of Bioinformatics and Life Science, Soongsil University, Seoul 156-743, Korea

* Corresponding Author: Chaeyoung Lee. Tel: +82-2-820-0455,

Fax: +82-2-824-4383, E-mail:

Submitted Apr. 17, 2014; Revised Jun. 24, 2014; Accepted Jul. 8, 2014
Table 1. Basic statistics for body weights (kg) of Korean cattle
analyzed in this study

Age in  Mean    SD     CV     Minimum   Maximum

6       167.2   28.7   17.1   92.0      262.5
9       231.1   28.3   12.3   153.7     314.5
12      326.0   39.8   12.2   221.2     451.6
15      395.1   37.1   9.4    292.2     510.6
18      489.2   44.4   9.1    370.0     617.5
21      567.7   53.0   9.3    422.3     723.0
24      650.1   61.8   9.5    454.3     818.6

SD, standard deviation; CV, coefficient of variation.
COPYRIGHT 2014 Asian - Australasian Association of Animal Production Societies
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Ryu, Jihye; Lee, Chaeyoung
Publication:Asian - Australasian Journal of Animal Sciences
Article Type:Report
Date:Nov 1, 2014
Previous Article:A review of recent developments in buffalo reproduction--a review.
Next Article:Polymorphism, genetic effect and association with egg production traits of chicken matrix metalloproteinases 9 promoter.

Terms of use | Copyright © 2017 Farlex, Inc. | Feedback | For webmasters