Marker-Based Estimation of Heritability in Immortal Populations

Abstract

Heritability is a central parameter in quantitative genetics, from both an evolutionary and a breeding perspective. For plant traits heritability is traditionally estimated by comparing within-and between-genotype variability. This approach estimates broad-sense heritability and does not account for different genetic relatedness. With the availability of high-density markers there is growing interest in marker-based estimates of narrow-sense heritability, using mixed models in which genetic relatedness is estimated from genetic markers. Such estimates have received much attention in human genetics but are rarely reported for plant traits. A major obstacle is that current methodology and software assume a single phenotypic value per genotype, hence requiring genotypic means. An alternative that we propose here is to use mixed models at the individual plant or plot level. Using statistical arguments, simulations, and real data we investigate the feasibility of both approaches and how these affect genomic prediction with the best linear unbiased predictor and genome-wide association studies. Heritability estimates obtained from genotypic means had very large standard errors and were sometimes biologically unrealistic. Mixed models at the individual plant or plot level produced more realistic estimates, and for simulated traits standard errors were up to 13 times smaller. Genomic prediction was also improved by using these mixed models, with up to a 49% increase in accuracy. For genome-wide association studies on simulated traits, the use of individual plant data gave almost no increase in power. The new methodology is applicable to any complex trait where multiple replicates of individual genotypes can be scored. This includes important agronomic crops, as well as bacteria and fungi.KEYWORDS marker-based estimation of heritability; GWAS; genomic prediction; Arabidopsis thaliana; one-vs. two-stage approaches N ARROW-SENSE heritability is an important parameter in quantitative genetics, determining the response to selection and representing the proportion of phenotypic variance that is due to additive genetic effects (Jacquard 1983;Ritland 1996;Visscher et al. 2006Visscher et al. , 2008Holland et al. 2010;Sillanpaa 2011). This definition of heritability goes back to Fisher (1918) and Wright (1920) almost a century ago. In plant species for which replicates of the same genotype are available (inbred lines, doubled haploids, clones), a different form of heritability, broadsense heritability, is traditionally estimated by the intraclass correlation coefficient for genotypic effects, using estimates for within-and between-genotype variance. Broad-sense heritability is also referred to as repeatability and gives the proportion of phenotypic variance explained by heritable (additive) and nonheritable (dominance, epistasis) genetic variance.With the arrival of high-density genotyping there is growing interest in marker-based estimation of narrow-sense heritability (WTCCC 2007;Yang et al. 2010Yang et al. , 2011Vatti...

Genome-Wide Association Mapping Reveals That Specific and Pleiotropic Regulatory Mechanisms Fine-Tune Central Metabolism and Growth in Arabidopsis

“…Broad-sense heritability (H 2 ) for biomass, structural components, metabolites, and several enzyme activities (aINV, nINV, PEPC, AGP, UGP, NRVs, and GDH) was moderate to high (20-74%), with marker-based h 2 (Kruijer et al, 2015) very close to H 2 values. This indicates that additive genetic variation is the main source of phenotypic variation and that epistasis does not play a substantial role in the variation of these traits.…”

“…To identify which factors contribute to trait variation, we calculated both broad-sense heritability (H 2 ) and marker-based heritability (h 2 ) using individual replicates (Kruijer et al, 2015) and averages (i.e., line-based heritability) for each trait and experiment (Table 1; Supplemental Data Set 2 and Supplemental Figure 1C). For estimates of marker-based heritability, we compared six models with different kinship matrices (EMMAX, GAPIT, PK1, PK2, PK3, and PK4).…”

Genetic Variance Partitioning and Genome-Wide Prediction with Allele Dosage Information in Autotetraploid Potato

“…To assess the value of tetraploid allele dosage, predictions were made with "diploidized" marker data (Gdip), in which the three heterozygotes were recoded to be identical. (Kruijer et al 2015).…”

Genome-Wide Association Mapping and Genomic Prediction Elucidate the Genetic Architecture of Morphological Traits in Arabidopsis

“…To assess to what extent the morphological variation is defined by the underlying genetic variation, markerbased heritability (h 2 ) was estimated as described before (Kruijer et al, 2015). The value h 2 ranged from 0.42 for RGRbv to 0.93 for FT, with the majority of traits having a heritability higher than 0.70, indicating that most of the measured variation could be attributed to genetic factors (Supplemental Table S1).…”

“…where y i,j is the phenotypic response of replicate j of genotype i; m is the intercept; G = (G 1 ,.,G n ) is the vector of random genetic effects; and the errors E i,j have independent normal distributions with variance s 2 E , which is the residual variance for a single individual (Kruijer et al, 2015). The vector G has a multivariate N(0, s 2 A K) distribution, and the genetic relatedness matrix K is estimated from standardized SNP-scores.…”