An Evolutionary Perspective on Epistasis and the Missing Heritability

<div><p>The relative importance between additive and non-additive genetic variance has been widely argued in quantitative genetics. By approaching this question from an evolutionary perspective we show that, while additive variance can be maintained under selection at a low level for some patterns of epistasis, the majority of the genetic variance that will persist is actually non-additive. We propose that one reason that the problem of the “missing heritability” arises is because the additive genetic variation that is estimated to be contributing to the variance of a trait will most likely be an artefact of the non-additive variance that can be maintained over evolutionary time. In addition, it can be shown that even a small reduction in linkage disequilibrium between causal variants and observed SNPs rapidly erodes estimates of epistatic variance, leading to an inflation in the perceived importance of additive effects. We demonstrate that the perception of independent additive effects comprising the majority of the genetic architecture of complex traits is biased upwards and that the search for causal variants in complex traits under selection is potentially underpowered by parameterising for additive effects alone. Given dense SNP panels the detection of causal variants through genome-wide association studies may be improved by searching for epistatic effects explicitly.</p> </div

Multiple-trait QTL mapping for body and organ weights in a cross between NMRI8 and DB/2 mice

This publication is with permission of the rights owner freely accessible due to an Alliance licence and a national licence (funded by the DFG, German Research Foundation) respectively.Multiple-trait analyses have been shown to improve the detection of quantitative trait loci (QTLs) with multiple effects. Here we applied a multiple-trait approach on obesity- and growth-related traits that were surveyed in 275 F2 mice generated from an intercross between the high body weight selected line NMRI8 and DBA/2 as lean control. The parental lines differed 2·5-fold in body weight at the age of 6 weeks. Within the F2 population, the correlations between body weight and weights of abdominal fat weight, muscle, liver and kidney at the age of 6 weeks were about 0·8. A least squares multiple-trait QTL analysis was performed on these data to understand more precisely the cause of the genetic correlation between body weight, body composition traits and weights of inner organs. Regions on Chr 1, 2, 7 and 14 for body weights at different early ages and regions on Chr 1, 2, 4, 7, 14, 17 and 19 for organ weights at 6 weeks were found to have significant multiple effects at the genome-wide level.Peer Reviewe

Investigating pedigree- and SNP-associated components of heritability in a wild population of Soay sheep

Estimates of narrow sense heritability derived from genomic data that contain related individuals may be biased due to the within-family effects such as dominance, epistasis and common environmental factors. However, for many wild populations, removal of related individuals from the data would result in small sample sizes. In 2013, Zaitlen et al. proposed a method to estimate heritability in populations that include close relatives by simultaneously fitting an identity-by-state (IBS) genomic relatedness matrix (GRM) and an identity-by-descent (IBD) GRM. The IBD GRM is identical to the IBS GRM, except relatedness estimates below a specified threshold are set to 0. We applied this method to a sample of 8557 wild Soay sheep from St. Kilda, with genotypic information for 419,281 single nucleotide polymorphisms. We aimed to see how this method would partition heritability into population-level (IBS) and family-associated (IBD) variance for a range of genetic architectures, and so we focused on a mixture of polygenic and monogenic traits. We also implemented a variant of the model in which the IBD GRM was replaced by a GRM constructed from SNPs with low minor allele frequency to examine whether any additive genetic variance is captured by rare alleles. Whilst the inclusion of the IBD GRM did not significantly improve the fit of the model for the monogenic traits, it improved the fit for some of the polygenic traits, suggesting that dominance, epistasis and/or common environment not already captured by the non-genetic random effects fitted in our models may influence these traits.</p

Investigating pedigree- and SNP-associated components of heritability in a wild population of Soay sheep

Estimates of narrow sense heritability derived from genomic data that contain related individuals may be biased due to the within-family effects such as dominance, epistasis and common environmental factors. However, for many wild populations, removal of related individuals from the data would result in small sample sizes. In 2013, Zaitlen et al. proposed a method to estimate heritability in populations that include close relatives by simultaneously fitting an identity-by-state (IBS) genomic relatedness matrix (GRM) and an identity-by-descent (IBD) GRM. The IBD GRM is identical to the IBS GRM, except relatedness estimates below a specified threshold are set to 0. We applied this method to a sample of 8557 wild Soay sheep from St. Kilda, with genotypic information for 419,281 single nucleotide polymorphisms. We aimed to see how this method would partition heritability into population-level (IBS) and family-associated (IBD) variance for a range of genetic architectures, and so we focused on a mixture of polygenic and monogenic traits. We also implemented a variant of the model in which the IBD GRM was replaced by a GRM constructed from SNPs with low minor allele frequency to examine whether any additive genetic variance is captured by rare alleles. Whilst the inclusion of the IBD GRM did not significantly improve the fit of the model for the monogenic traits, it improved the fit for some of the polygenic traits, suggesting that dominance, epistasis and/or common environment not already captured by the non-genetic random effects fitted in our models may influence these traits.</p

Haplotyping via minimum recombinant paradigm

<p>Abstract</p> <p>Background</p> <p>Haplotypes can increase the power of gene detection over genotypes and are essential to estimate linkage disequilibrium.</p> <p>Methods</p> <p>Haplotyping was based on the minimum recombinant paradigm, whereby a phase is obtained only if it uniquely minimises the number of recombinants within a full sib family. Performance of this method was tested across three different data sets, consisting of genotypes and pedigree.</p> <p>Results</p> <p>The percentage of phased alleles ranged from ~80% to ~95%, and the percentage of correct phases reached ~99% in all cases. A measure of uncertainty was obtained via simulations. A partial haplotyping algorithm consisting of four deterministic rules was almost as effective as a full one consisting of six deterministic rules, and took up to 5 times less time to compute.</p> <p>Conclusion</p> <p>Haplotyping via the minimum recombinant paradigm is consistently reliable and computationally efficient. A single simulation is enough to produce a population-wide uncertainty estimate associated with a set of all reconstructed haplotypes.</p

The impact of SNP density on quantitative genetic analyses of body size traits in a wild population of Soay sheep

Understanding the genetic architecture underpinning quantitative traits in wild populations is pivotal to understanding the processes behind trait evolution. The ‘animal model’ is a popular method for estimating quantitative genetic parameters such as heritability and genetic correlation and involves fitting an estimate of relatedness between individuals in the study population. Genotypes at genome‐wide markers can be used to estimate relatedness; however, relatedness estimates vary with marker density, potentially affecting results. Increasing density of markers is also expected to increase the power to detect quantitative trait loci (QTL). In order to understand how the density of genetic markers affects the results of quantitative genetic analyses, we estimated heritability and performed genome‐wide association studies (GWAS) on five body size traits in an unmanaged population of Soay sheep using two different SNP densities: a dataset of 37,037 genotyped SNPs and an imputed dataset of 417,373 SNPs. Heritability estimates did not differ between the two SNP densities, but the high‐density imputed SNP dataset revealed four new SNP‐trait associations that were not found with the lower density dataset, as well as confirming all previously‐found QTL. We also demonstrated that fitting fixed and random effects in the same step as performing GWAS is a more powerful approach than pre‐correcting for covariates in a separate model

Comparative analysis of quantitative trait loci for body weight, growth rate and growth curve parameters from 3 to 72 weeks of age in female chickens of a broiler-layer cross

Background: Comparisons of quantitative trait loci (QTL) for growth and parameters of growth curves assist in understanding the genetics and ultimately the physiology of growth. Records of body weight at 3, 6, 12, 24, 48 and 72 weeks of age and growth rate between successive age intervals of about 500 F female chickens of the Roslin broiler-layer cross were available for analysis. These data were analysed to detect and compare QTL for body weight, growth rate and parameters of the Gompertz growth function.Results: Over 50 QTL were identified for body weight at specific ages and most were also detected in the nearest preceding and/or subsequent growth stage. The sum of the significant and suggestive additive effects for bodyweight at specific ages accounted for 23-43% of the phenotypic variation. A single QTL for body weight on chromosome 4 at 48 weeks of age had the largest additive effect (550.4 ± 68.0 g, 11.5% of the phenotypic variation) and a QTL at a similar position accounted 14.5% of the phenotypic variation at 12 weeks of age. Age specific QTL for growth rate were detected suggesting that there are specific genes that affect developmental processes during the different stages of growth. Relatively few QTL influencing Gompertz growth curve parameters were detected and overlapped with loci affecting growth rate. Dominance effects were generally not significant but from 12 weeks of age they exceeded the additive effect in a few cases. No evidence for epistatic QTL pairs was found.Conclusions: The results confirm the location for body weight and body weight gain during growth that were identified in previous studies and were consistent with QTL for the parameters of the Gompertz growth function. Chromosome 4 explained a relatively large proportion of the observed growth variation across the different ages, and also harboured most of the detected QTL for Gompertz parameters, confirming its importance in controlling growth. Very few QTL were detected for body weight or gain at 48 and 72 weeks of age, probably reflecting the effect of differences in reproduction and random environmental effects