The Genomic Architecture of Interactions Between Natural Genetic Polymorphisms and Environments in Yeast Growth

The loci of environmental adaptation in a model eukaryote

While the underlying genetic changes have been uncovered in some cases of adaptive evolution, the lack of a systematic study prevents a general understanding of the genomic basis of adaptation. For example, it is unclear whether protein-coding or noncoding mutations are more important to adaptive evolution and whether adaptations to different environments are brought by genetic changes distributed in diverse genes and biological processes or concentrated in a core set. We here perform laboratory evolution of 3360 Saccharomyces cerevisiae populations in 252 environments of varying levels of stress. We find the yeast adaptations to be primarily fueled by large-effect coding mutations overrepresented in a relatively small gene set, despite prevalent antagonistic pleiotropy across environments. Populations generally adapt faster in more stressful environments, partly because of greater benefits of the same mutations in more stressful environments. These and other findings from this model eukaryote help unravel the genomic principles of environmental adaptation.

What Has Genomics Taught An Evolutionary Biologist?

Genomics, an interdisciplinary field of biology on the structure, function, and evolution of genomes, has revolutionized many subdisciplines of life sciences, including my field of evolutionary biology, by supplying huge data, bringing high-throughput technologies, and offering a new approach to biology. In this review, I describe what I have learned from genomics and highlight the fundamental knowledge and mechanistic insights gained. I focus on three broad topics that are central to evolutionary biology and beyond-variation, interaction, and selection-and use primarily my own research and study subjects as examples. In the next decade or two, I expect that the most important contributions of genomics to evolutionary biology will be to provide genome sequences of nearly all known species on Earth, facilitate high-throughput phenotyping of natural variants and systematically constructed mutants for mapping genotype-phenotype-fitness landscapes, and assist the determination of causality in evolutionary processes using experimental evolution.

Interpretation of Manhattan Plots and Other Outputs of Genome-Wide Association Studies

With increasing marker density, estimation of recombination rate between a marker and a causal mutation using linkage analysis becomes less important. Instead, linkage disequilibrium (LD) becomes the major indicator for gene mapping through genome-wide association studies (GWAS). In addition to the linkage between the marker and the causal mutation, many other factors may contribute to the LD, including population structure and cryptic relationships among individuals. As statistical methods and software evolve to improve statistical power and computing speed in GWAS, the corresponding outputs must also evolve to facilitate the interpretation of input data, the analytical process, and final association results. In this chapter, our descriptions focus on (1) considerations in creating a Manhattan plot displaying the strength of LD and locations of markers across a genome; (2) criteria for genome-wide significance threshold and the different appearance of Manhattan plots in single-locus and multiple-locus models; (3) exploration of population structure and kinship among individuals; (4) quantile-quantile (QQ) plot; (5) LD decay across the genome and LD between the associated markers and their neighbors; (6) exploration of individual and marker information on Manhattan and QQ plots via interactive visualization using HTML. The ultimate objective of this chapter is to help users to connect input data to GWAS outputs to balance power and false positives, and connect GWAS outputs to the selection of candidate genes using LD extent.

Different genetic basis for alcohol dehydrogenase activity and plasticity in a novel alcohol environment for Drosophila melanogaster

Phenotypic plasticity is known to enhance population persistence, facilitate adaptive evolution and initiate novel phenotypes in novel environments. How plasticity can contribute or hinder adaptation to different environments hinges on its genetic architecture. Even though plasticity in many traits is genetically controlled, whether and how plasticity's genetic architecture might change in novel environments is still unclear. Because much of gene expression can be environmentally influenced, each environment may trigger different sets of genes that influence a trait. Using a quantitative trait loci (QTL) approach, we investigated the genetic basis of plasticity in a classic functional trait, alcohol dehydrogenase (ADH) activity in D. melanogaster, across both historical and novel alcohol environments. Previous research in D. melanogaster has also demonstrated that ADH activity is plastic in response to alcohol concentration in substrates used by both adult flies and larvae. We found that across all environments tested, ADH activity was largely influenced by a single QTL encompassing the Adh-coding gene and its known regulatory locus, delta-1. After controlling for the allelic variation of the Adh and delta-1 loci, we found additional but different minor QTLs in the 0 and 14% alcohol environments. In contrast, we discovered no major QTL for plasticity itself, including the Adh locus, regardless of the environmental gradients. This suggests that plasticity in ADH activity is likely influenced by many loci with small effects, and that the Adh locus is not environmentally sensitive to dietary alcohol.

Antagonistic pleiotropy conceals molecular adaptations in changing environments

The importance of positive selection in molecular evolution is debated. Evolution experiments under invariant laboratory conditions typically show a higher rate of nonsynonymous nucleotide changes than the rate of synonymous changes, demonstrating prevalent molecular adaptations. Natural evolution inferred from genomic comparisons, however, almost always exhibits the opposite pattern even among closely related conspecifics, which is indicative of a paucity of positive selection. Here we hypothesize that this apparent contradiction is at least in part attributable to ubiquitous and frequent environmental changes in nature, causing nonsynonymous mutations that are beneficial at one time to become deleterious soon after because of antagonistic pleiotropy and hindering their fixations relative to synonymous mutations despite continued population adaptations. To test this hypothesis, we performed yeast evolution experiments in changing and corresponding constant environments, followed by genome sequencing of the evolving populations. We observed a lower nonsynonymous to synonymous rate ratio in antagonistic changing environments than in the corresponding constant environments, and the population dynamics of mutations supports our hypothesis. These findings and the accompanying population genetic simulations suggest that molecular adaptation is consistently underestimated in nature due to the antagonistic fitness effects of mutations in changing environments.

Quantitative Trait Nucleotides Impacting the Technological Performances of Industrial Saccharomyces cerevisiae Strains

The budding yeast Saccharomyces cerevisiae is certainly the prime industrial microorganism and is related to many biotechnological applications including food fermentations, biofuel production, green chemistry, and drug production. A noteworthy characteristic of this species is the existence of subgroups well adapted to specific processes with some individuals showing optimal technological traits. In the last 20 years, many studies have established a link between quantitative traits and single-nucleotide polymorphisms found in hundreds of genes. These natural variations constitute a pool of QTNs (quantitative trait nucleotides) that modulate yeast traits of economic interest for industry. By selecting a subset of genes functionally validated, a total of 284 QTNs were inventoried. Their distribution across pan and core genome and their frequency within the 1,011 Saccharomyces cerevisiae genomes were analyzed. We found that 150 of the 284 QTNs have a frequency lower than 5%, meaning that these variants would be undetectable by genome-wide association studies (GWAS). This analysis also suggests that most of the functional variants are private to a subpopulation, possibly due to their adaptive role to specific industrial environment. In this review, we provide a literature survey of their phenotypic impact and discuss the opportunities and the limits of their use for industrial strain selection.

Patterns and Mechanisms of Diminishing Returns from Beneficial Mutations

Diminishing returns epistasis causes the benefit of the same advantageous mutation smaller in fitter genotypes and is frequently observed in experimental evolution. However, its occurrence in other contexts, environment dependence, and mechanistic basis are unclear. Here, we address these questions using 1,005 sequenced segregants generated from a yeast cross. Under each of 47 examined environments, 66-92% of tested polymorphisms exhibit diminishing returns epistasis. Surprisingly, improving environment quality also reduces the benefits of advantageous mutations even when fitness is controlled for, indicating the necessity to revise the global epistasis hypothesis. We propose that diminishing returns originates from the modular organization of life where the contribution of each functional module to fitness is determined jointly by the genotype and environment and has an upper limit, and demonstrate that our model predictions match empirical observations. These findings broaden the concept of diminishing returns epistasis, reveal its generality and potential cause, and have important evolutionary implications.

CCR5-∆32 is deleterious in the homozygous state in humans

We use the genotyping and death register information of 409,693 individuals of British ancestry to investigate fitness effects of the CCR5-∆32 mutation. We estimate a 21% increase in the all-cause mortality rate in individuals who are homozygous for the ∆32 allele. A deleterious effect of the ∆32/∆32 mutation is also independently supported by a significant deviation from the Hardy-Weinberg equilibrium (HWE) due to a deficiency of ∆32/∆32 individuals at the time of recruitment.

Neutral Theory and Phenotypic Evolution

Although the neutral theory of molecular evolution was proposed to explain DNA and protein sequence evolution, in principle it could also explain phenotypic evolution. Nevertheless, overall, phenotypes should be less likely than genotypes to evolve neutrally. I propose that, when phenotypic traits are stratified according to a hierarchy of biological organization, the fraction of evolutionary changes in phenotype that are adaptive rises with the phenotypic level considered. Consistently, molecular traits are frequently found to evolve neutrally whereas a large, random set of organismal traits were recently reported to vary largely adaptively. Many more studies of unbiased samples of phenotypic traits are needed to test the general validity of this hypothesis.

Environment-dependent pleiotropic effects of mutations on the maximum growth rate r and carrying capacity K of population growth

Maximum growth rate per individual (r) and carrying capacity (K) are key life-history traits that together characterize the density-dependent population growth and therefore are crucial parameters of many ecological and evolutionary theories such as r/K selection. Although r and K are generally thought to correlate inversely, both r/K tradeoffs and trade-ups have been observed. Nonetheless, neither the conditions under which each of these relationships occur nor the causes of these relationships are fully understood. Here, we address these questions using yeast as a model system. We estimated r and K using the growth curves of over 7,000 yeast recombinants in nine environments and found that the r-K correlation among genotypes changes from 0.53 to -0.52 with the rise of environment quality, measured by the mean r of all genotypes in the environment. We respectively mapped quantitative trait loci (QTLs) for r and K in each environment. Many QTLs simultaneously influence r and K, but the directions of their effects are environment dependent such that QTLs tend to show concordant effects on the two traits in poor environments but antagonistic effects in rich environments. We propose that these contrasting trends are generated by the relative impacts of two factors-the tradeoff between the speed and efficiency of ATP production and the energetic cost of cell maintenance relative to reproduction-and demonstrate an agreement between model predictions and empirical observations. These results reveal and explain the complex environment dependency of the r-K relationship, which bears on many ecological and evolutionary phenomena and has biomedical implications.

Dissection of the molecular bases of genotype x environment interactions: a study of phenotypic plasticity of Saccharomyces cerevisiae in grape juices

Background

The ability of a genotype to produce different phenotypes according to its surrounding environment is known as phenotypic plasticity. Within different individuals of the same species, phenotypic plasticity can vary greatly. This contrasting response is caused by gene-by-environment interactions (GxE). Understanding GxE interactions is particularly important in agronomy, since selected breeds and varieties may have divergent phenotypes according to their growing environment. Industrial microbes such as Saccharomyces cerevisiae are also faced with a large range of fermentation conditions that affect their technological properties. Finding the molecular determinism of such variations is a critical task for better understanding the genetic bases of phenotypic plasticity and can also be helpful in order to improve breeding methods.

Results

In this study we implemented a QTL mapping program using two independent cross (~ 100 progeny) in order to investigate the molecular basis of yeast phenotypic response in a wine fermentation context. Thanks to whole genome sequencing approaches, both crosses were genotyped, providing saturated genetic maps of thousands of markers. Linkage analyses allowed the detection of 78 QTLs including 21 with significant interaction with the environmental conditions. Molecular dissection of a major QTL demonstrated that the sulfite pump Ssu1p has a pleiotropic effect and impacts the phenotypic plasticity of several traits.

Conclusions

The detection of QTLs and their interactions with environment emphasizes the complexity of yeast industrial traits. The validation of the interaction of SSU1 allelic variants with the nature of the fermented juice increases knowledge about the impact of the sulfite pump during fermentation. All together these results pave the way for exploiting and deciphering the genetic determinism of phenotypic plasticity.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5145-4) contains supplementary material, which is available to authorized users.

Multi-environment fitness landscapes of a tRNA gene

A fitness landscape (FL) describes the genotype-fitness relationship in a given environment. To explain and predict evolution, it is imperative to measure the FL in multiple environments because the natural environment changes frequently. Using a high-throughput method that combines precise gene replacement with next-generation sequencing, we determine the in vivo FL of a yeast tRNA gene comprising over 23,000 genotypes in four environments. Although genotype-by-environment interaction (G×E) is abundantly detected, its pattern is so simple that we can transform an existing FL to that in a new environment with fitness measures of only a few genotypes in the new environment. Under each environment, we observe prevalent, negatively biased epistasis between mutations (G×G). Epistasis-by-environment interaction (G×G×E) is also prevalent, but trends in epistasis difference between environments are predictable. Our study thus reveals simple rules underlying seemingly complex FLs, opening the door to understanding and predicting FLs in general.

Gene-gene and gene-environment interactions in complex traits in yeast

One of the fundamental questions in biology is how the genotype regulates the phenotype. An increasing number of studies indicate that, in most cases, the effect of a genetic locus on the phenotype is context-dependent, i.e. it is influenced by the genetic background and the environment in which the phenotype is measured. Still, the majority of the studies, in both model organisms and humans, that map the genetic regulation of phenotypic variation in complex traits primarily identify additive loci with independent effects. This does not reflect an absence of the contribution of genetic interactions to phenotypic variation, but instead is a consequence of the technical limitations in mapping gene-gene interactions (GGI) and gene-environment interactions (GEI). Yeast, with its detailed molecular understanding, diverse population genomics and ease of genetic manipulation, is a unique and powerful resource to study the contributions of GGI and GEI in the regulation of phenotypic variation. Here we review studies in yeast that have identified GGI and GEI that regulate phenotypic variation, and discuss the contribution of these findings in explaining missing heritability of complex traits, and how observations from these GGI and GEI studies enhance our understanding of the mechanisms underlying genetic robustness and adaptability that shape the architecture of the genotype-phenotype map.