Differential Regulation of Cryptic Genetic Variation Shapes the Genetic Interactome Underlying Complex Traits

Dissecting the Genetic Regulation of Yeast Growth Plasticity in Response to Environmental Changes

Variable individual responses to environmental changes, such as phenotype plasticity, are heritable, with some genotypes being robust and others plastic. This variation for plasticity contributes to variance in complex traits as genotype-by-environment interactions (G × E). However, the genetic basis of this variability in responses to the same external stimuli is still largely unknown. In an earlier study of a large haploid segregant yeast population, genotype-by-genotype-by-environment interactions were found to make important contributions to the release of genetic variation in growth responses to alterations of the growth medium. Here, we explore the genetic basis for heritable variation of different measures of phenotype plasticity in the same dataset. We found that the central loci in the environmentally dependent epistatic networks were associated with overall measures of plasticity, while the specific measures of plasticity identified a more diverse set of loci. Based on this, a rapid one-dimensional genome-wide association (GWA) approach to overall plasticity is proposed as a strategy to efficiently identify key epistatic loci contributing to the phenotype plasticity. The study thus provided both analytical strategies and a deeper understanding of the complex genetic regulation of phenotype plasticity in yeast growth.

Dynamic genetic architecture of yeast response to environmental perturbation shed light on origin of cryptic genetic variation

Cryptic genetic variation could arise from, for example, Gene-by-Gene (G-by-G) or Gene-by-Environment (G-by-E) interactions. The underlying molecular mechanisms and how they influence allelic effects and the genetic variance of complex traits is largely unclear. Here, we empirically explored the role of environmentally influenced epistasis on the suppression and release of cryptic variation by reanalysing a dataset of 4,390 haploid yeast segregants phenotyped on 20 different media. The focus was on 130 epistatic loci, each contributing to segregant growth in at least one environment and that together explained most (69–100%) of the narrow sense heritability of growth in the individual environments. We revealed that the epistatic growth network reorganised upon environmental changes to alter the estimated marginal (additive) effects of the individual loci, how multi-locus interactions contributed to individual segregant growth and the level of expressed genetic variance in growth. The estimated additive effects varied most across environments for loci that were highly interactive network hubs in some environments but had few or no interactors in other environments, resulting in changes in total genetic variance across environments. This environmentally dependent epistasis was thus an important mechanism for the suppression and release of cryptic variation in this population. Our findings increase the understanding of the complex genetic mechanisms leading to cryptic variation in populations, providing a basis for future studies on the genetic maintenance of trait robustness and development of genetic models for studying and predicting selection responses for quantitative traits in breeding and evolution.

Many biological traits are polygenic, with complex interplay between underlying genes and the surrounding environment. As a result, individuals with the same allele might have distinctive phenotypes due to differences in the polygenic background and/or the environment. Such differences often create additional genetic variation that is highly relevant to quantitative and evolutionary genetics by limiting our ability to accurately predict the phenotypes in medical or agricultural applications and providing opportunities for long term evolution. Previously, yeast growth regulating genes were found to be organised in large interacting networks. Here, we found that these networks were reorganised upon environmental changes, and that this resulted in altered effect sizes of individual genes, and how the whole network contributed to growth and the level of total genetic variance, providing a basis for future studies on the genetic maintenance of trait robustness and development of genetic models for studying and predicting selection responses for quantitative traits.

Extensive impact of low-frequency variants on the phenotypic landscape at population-scale

Genome-wide association studies (GWAS) allow to dissect complex traits and map genetic variants, which often explain relatively little of the heritability. One potential reason is the preponderance of undetected low-frequency variants. To increase their allele frequency and assess their phenotypic impact in a population, we generated a diallel panel of 3025 yeast hybrids, derived from pairwise crosses between natural isolates and examined a large number of traits. Parental versus hybrid regression analysis showed that while most phenotypic variance is explained by additivity, a third is governed by non-additive effects, with complete dominance having a key role. By performing GWAS on the diallel panel, we found that associated variants with low frequency in the initial population are overrepresented and explain a fraction of the phenotypic variance as well as an effect size similar to common variants. Overall, we highlighted the relevance of low-frequency variants on the phenotypic variation.

QTL Mapping on a Background of Variance Heterogeneity

Standard QTL mapping procedures seek to identify genetic loci affecting the phenotypic mean while assuming that all individuals have the same residual variance. But when the residual variance differs systematically between groups, perhaps due to a genetic or environmental factor, such standard procedures can falter: in testing for QTL associations, they attribute too much weight to observations that are noisy and too little to those that are precise, resulting in reduced power and and increased susceptibility to false positives. The negative effects of such "background variance heterogeneity" (BVH) on standard QTL mapping have received little attention until now, although the subject is closely related to work on the detection of variance-controlling genes. Here we use simulation to examine how BVH affects power and false positive rate for detecting QTL affecting the mean (mQTL), the variance (vQTL), or both (mvQTL). We compare linear regression for mQTL and Levene's test for vQTL, with tests more recently developed, including tests based on the double generalized linear model (DGLM), which can model BVH explicitly. We show that, when used in conjunction with a suitable permutation procedure, the DGLM-based tests accurately control false positive rate and are more powerful than the other tests. We also find that some adverse effects of BVH can be mitigated by applying a rank inverse normal transform. We apply our novel approach, which we term "mean-variance QTL mapping", to publicly available data on a mouse backcross and, after accommodating BVH driven by sire, detect a new mQTL for bodyweight.

QTL Mapping on a Background of Variance Heterogeneity

Standard QTL mapping procedures seek to identify genetic loci affecting the phenotypic mean while assuming that all individuals have the same residual variance. But when the residual variance differs systematically between groups, perhaps due to a genetic or environmental factor, such standard procedures can falter: in testing for QTL associations, they attribute too much weight to observations that are noisy and too little to those that are precise, resulting in reduced power and and increased susceptibility to false positives. The negative effects of such "background variance heterogeneity" (BVH) on standard QTL mapping have received little attention until now, although the subject is closely related to work on the detection of variance-controlling genes. Here we use simulation to examine how BVH affects power and false positive rate for detecting QTL affecting the mean (mQTL), the variance (vQTL), or both (mvQTL). We compare linear regression for mQTL and Levene's test for vQTL, with tests more recently developed, including tests based on the double generalized linear model (DGLM), which can model BVH explicitly. We show that, when used in conjunction with a suitable permutation procedure, the DGLM-based tests accurately control false positive rate and are more powerful than the other tests. We also find that some adverse effects of BVH can be mitigated by applying a rank inverse normal transform. We apply our novel approach, which we term "mean-variance QTL mapping", to publicly available data on a mouse backcross and, after accommodating BVH driven by sire, detect a new mQTL for bodyweight.

Resolving the Complex Genetic Basis of Phenotypic Variation and Variability of Cellular Growth

In all organisms, the majority of traits vary continuously between individuals. Explaining the genetic basis of quantitative trait variation requires comprehensively accounting for genetic and nongenetic factors as well as their interactions. The growth of microbial cells can be characterized by a lag duration, an exponential growth phase, and a stationary phase. Parameters that characterize these growth phases can vary among genotypes (phenotypic variation), environmental conditions (phenotypic plasticity), and among isogenic cells in a given environment (phenotypic variability). We used a high-throughput microscopy assay to map genetic loci determining variation in lag duration and exponential growth rate in growth rate-limiting and nonlimiting glucose concentrations, using segregants from a cross of two natural isolates of the budding yeast, Saccharomyces cerevisiae We find that some quantitative trait loci (QTL) are common between traits and environments whereas some are unique, exhibiting gene-by-environment interactions. Furthermore, whereas variation in the central tendency of growth rate or lag duration is explained by many additive loci, differences in phenotypic variability are primarily the result of genetic interactions. We used bulk segregant mapping to increase QTL resolution by performing whole-genome sequencing of complex mixtures of an advanced intercross mapping population grown in selective conditions using glucose-limited chemostats. We find that sequence variation in the high-affinity glucose transporter HXT7 contributes to variation in growth rate and lag duration. Allele replacements of the entire locus, as well as of a single polymorphic amino acid, reveal that the effect of variation in HXT7 depends on genetic, and allelic, background. Amplifications of HXT7 are frequently selected in experimental evolution in glucose-limited environments, but we find that HXT7 amplifications result in antagonistic pleiotropy that is absent in naturally occurring variants of HXT7 Our study highlights the complex nature of the genotype-to-phenotype map within and between environments.