Molecular specificity, convergence and constraint shape adaptive evolution in nutrient-poor environments

The Whi2-Psr1-Psr2 complex selectively regulates TORC1 and autophagy under low leucine conditions but not nitrogen depletion

Amino acids and ammonia serve as sources of nitrogen for cell growth and were previously thought to have similar effects on yeast. Consistent with this idea, depletion of either of these two nitrogen sources inhibits the target of rapamycin complex 1 (TORC1), leading to induction of macroautophagy/autophagy and inhibition of cell growth. In this study, we show that Whi2 and the haloacid dehalogenase (HAD)-type phosphatases Psr1 and Psr2 distinguish between these two nitrogen sources in Saccharomyces cerevisiae, as the Whi2-Psr1-Psr2 complex inhibits TORC1 in response to low leucine but not in the absence of nitrogen. In contrast, a parallel pathway controlled by Npr2 and Npr3, components of the Seh1-associated complex inhibiting TORC1 (SEACIT), suppress TORC1 under both low leucine- and nitrogen-depletion conditions. Co-immunoprecipitations with mutants of Whi2, Psr1, Psr2 and fragments of Tor1 support the model that Whi2 recruits Psr1 and Psr2 to TORC1. In accordance, the interaction between Whi2 and Tor1 appears to increase under low leucine but decreases under nitrogen-depletion conditions. Although the targets of Psr1 and Psr2 phosphatases are not known, mutation of their active sites abolishes their inhibitory effects on TORC1. Consistent with the conservation of HAD phosphatases across species, human HAD phosphatases CTDSP1 (CTD small phosphatase 1), CTDSP2, and CTDSPL can functionally replace Psr1 and Psr2 in yeast, restoring TORC1 inhibition and autophagy activation in response to low leucine conditions.

Multilevel Gene Expression Changes in Lineages Containing Adaptive Copy Number Variants

Copy number variants (CNVs) are an important class of genetic variation that can mediate rapid adaptive evolution. Whereas, CNVs can increase the relative fitness of the organism, they can also incur a cost due to the associated increased gene expression and repetitive DNA. We previously evolved populations of Saccharomyces cerevisiae over hundreds of generations in glutamine-limited (Gln-) chemostats and observed the recurrent evolution of CNVs at the GAP1 locus. To understand the role that gene expression plays in adaptation, both in relation to the adaptation of the organism to the selective condition and as a consequence of the CNV, we measured the transcriptome, translatome, and proteome of 4 strains of evolved yeast, each with a unique CNV, and their ancestor in Gln- chemostats. We find CNV-amplified genes correlate with higher mRNA abundance; however, this effect is reduced at the level of the proteome, consistent with post-transcriptional dosage compensation. By normalizing each level of gene expression by the abundance of the preceding step we were able to identify widespread differences in the efficiency of each level of gene expression. Genes with significantly different translational efficiency were enriched for potential regulatory mechanisms including either upstream open reading frames, RNA-binding sites for Ssd1, or both. Genes with lower protein expression efficiency were enriched for genes encoding proteins in protein complexes. Taken together, our study reveals widespread changes in gene expression at multiple regulatory levels in lineages containing adaptive CNVs highlighting the diverse ways in which genome evolution shapes gene expression.

Template switching during DNA replication is a prevalent source of adaptive gene amplification

Copy number variants (CNVs) are an important source of genetic variation underlying rapid adaptation and genome evolution. Whereas point mutation rates vary with genomic location and local DNA features, the role of genome architecture in the formation and evolutionary dynamics of CNVs is poorly understood. Previously, we found the GAP1 gene in Saccharomyces cerevisiae undergoes frequent amplification and selection in glutamine-limitation. The gene is flanked by two long terminal repeats (LTRs) and proximate to an origin of DNA replication (autonomously replicating sequence, ARS), which likely promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we evolved engineered strains lacking either the adjacent LTRs, ARS, or all elements in glutamine-limited chemostats. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. Removal of local DNA elements significantly impacts the fitness effect of GAP1 CNVs and the rate of adaptation. In 177 CNV lineages, across all four strains, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, distal ones mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination can mediate gene amplification following de novo retrotransposon events. Our study reveals that template switching during DNA replication is a prevalent source of adaptive CNVs.

A high-resolution two-step evolution experiment in yeast reveals a shift from pleiotropic to modular adaptation

Evolution by natural selection is expected to be a slow and gradual process. In particular, the mutations that drive evolution are predicted to be small and modular, incrementally improving a small number of traits. However, adaptive mutations identified early in microbial evolution experiments, cancer, and other systems often provide substantial fitness gains and pleiotropically improve multiple traits at once. We asked whether such pleiotropically adaptive mutations are common throughout adaptation or are instead a rare feature of early steps in evolution that tend to target key signaling pathways. To do so, we conducted barcoded second-step evolution experiments initiated from 5 first-step mutations identified from a prior yeast evolution experiment. We then isolated hundreds of second-step mutations from these evolution experiments, measured their fitness and performance in several growth phases, and conducted whole genome sequencing of the second-step clones. Here, we found that while the vast majority of mutants isolated from the first-step of evolution in this condition show patterns of pleiotropic adaptation-improving both performance in fermentation and respiration growth phases-second-step mutations show a shift towards modular adaptation, mostly improving respiration performance and only rarely improving fermentation performance. We also identified a shift in the molecular basis of adaptation from genes in cellular signaling pathways towards genes involved in respiration and mitochondrial function. Our results suggest that the genes in cellular signaling pathways may be more likely to provide large, adaptively pleiotropic benefits to the organism due to their ability to coherently affect many phenotypes at once. As such, these genes may serve as the source of pleiotropic adaptation in the early stages of evolution, and once these become exhausted, organisms then adapt more gradually, acquiring smaller, more modular mutations.

Genotype-by-environment interactions shape ubiquitin-proteasome system activity

In genotype-by-environment interactions (GxE), the effect of a genetic variant on a trait depends on the environment. GxE influences numerous organismal traits across eukaryotic life. However, we have a limited understanding of how GxE shapes the molecular processes that give rise to organismal traits. Here, we characterized how GxE shapes protein degradation, an essential molecular process that influences numerous aspects of cellular and organismal physiology. Using the yeast Saccharomyces cerevisiae, we characterized GxE in the activity of the ubiquitin-proteasome system (UPS), the primary protein degradation system in eukaryotes. By mapping genetic influences on the degradation of six substrates that engage multiple distinct UPS pathways across eight diverse environments, we discovered extensive GxE in the genetics of UPS activity. Hundreds of locus effects on UPS activity varied depending on the substrate, the environment, or both. Most of these cases corresponded to loci that were present in one environment but not another ("presence / absence" GxE), while a smaller number of loci had opposing effects in different environments ("sign change" GxE). The number of loci exhibiting GxE, their genomic location, and the type of GxE (presence / absence or sign change) varied across UPS substrates. Loci exhibiting GxE were clustered at genomic regions that contain core UPS genes and especially at regions containing variation that affects the expression of thousands of genes, suggesting indirect contributions to UPS activity. Our results reveal highly complex interactions at the level of substrates and environments in the genetics of protein degradation.

Massively parallel experimental interrogation of natural variants in ancient signaling pathways reveals both purifying selection and local adaptation

The nature of standing genetic variation remains a central debate in population genetics, with differing perspectives on whether common variants are mostly neutral or have functional effects. We address this question by directly mapping the fitness effects of over 9,000 natural variants in the Ras/PKA and TOR/Sch9 pathways-key regulators of cell proliferation in eukaryotes-across four conditions in Saccharomyces cerevisiae. While many variants are neutral in our assay, on the order of 3,500 exhibited significant fitness effects. These non-neutral variants tend to be missense and affect conserved, more densely packed, and less solvent-exposed protein regions. They are also typically younger, occur at lower frequencies, and more often found in heterozygous states, suggesting they are subject to purifying selection. A substantial fraction of non-neutral variants showing strong fitness effects in our experiments, however, is present at high frequencies in the population. These variants show signs of local adaptation as they tend to be found specifically in domesticated strains adapted to human-made environments. Our findings support the view that while common variants are often neutral, a significant proportion have adaptive functional consequences and are driven into the population by local positive selection. This study highlights the potential to explore the functional effects of natural genetic variation on a genome scale with quantitative fitness measurements in the laboratory, bridging the gap between population genetics and functional genomics to understand evolutionary dynamics in the wild.

Template switching during DNA replication is a prevalent source of adaptive gene amplification

Copy number variants (CNVs)-gains and losses of genomic sequences-are an important source of genetic variation underlying rapid adaptation and genome evolution. However, despite their central role in evolution little is known about the factors that contribute to the structure, size, formation rate, and fitness effects of adaptive CNVs. Local genomic sequences are likely to be an important determinant of these properties. Whereas it is known that point mutation rates vary with genomic location and local DNA sequence features, the role of genome architecture in the formation, selection, and the resulting evolutionary dynamics of CNVs is poorly understood. Previously, we have found that the GAP1 gene in Saccharomyces cerevisiae undergoes frequent and repeated amplification and selection under long-term experimental evolution in glutamine-limiting conditions. The GAP1 gene has a unique genomic architecture consisting of two flanking long terminal repeats (LTRs) and a proximate origin of DNA replication (autonomously replicating sequence, ARS), which are likely to promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we performed experimental evolution in glutamine-limited chemostats using engineered strains lacking either the adjacent LTRs, ARS, or all elements. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. We find that although GAP1 CNVs repeatedly form and sweep to high frequency in strains with modified genome architecture, removal of local DNA elements significantly impacts the rate and fitness effect of CNVs and the rate of adaptation. We performed genome sequence analysis to define the molecular mechanisms of CNV formation for 177 CNV lineages. We find that across all four strain backgrounds, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, a distal ARS can mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination mechanisms still mediate gene amplification following de novo insertion of retrotransposon elements at the locus. Our study demonstrates the remarkable plasticity of the genome and reveals that template switching during DNA replication is a frequent source of adaptive CNVs.

Multilevel gene expression changes in lineages containing adaptive copy number variants

Copy-number variants (CNVs) are an important class of recurrent variants that mediate adaptive evolution. While CNVs can increase the relative fitness of the organism, they can also incur a cost. We previously evolved populations of Saccharomyces cerevisiae over hundreds of generations in glutamine-limited (Gln-) chemostats and observed the recurrent evolution of CNVs at the GAP1 locus. To understand the role that expression plays in adaptation, both in relation to the adaptation of the organism to the selective condition, and as a consequence of the CNV, we measured the transcriptome, translatome, and proteome of 4 strains of evolved yeast, each with a unique CNV, and their ancestor in Gln- conditions. We find CNV-amplified genes correlate with higher RNA abundance; however, this effect is reduced at the level of the proteome, consistent with post-transcriptional dosage compensation. By normalizing each level of expression by the abundance of the preceding step we were able to identify widespread divergence in the efficiency of each step in the gene in the efficiency of each step in gene expression. Genes with significantly different translational efficiency were enriched for potential regulatory mechanisms including either upstream open reading frames, RNA binding sites for SSD1, or both. Genes with lower protein expression efficiency were enriched for genes encoding proteins in protein complexes. Taken together, our study reveals widespread changes in gene expression at multiple regulatory levels in lineages containing adaptive CNVs highlighting the diverse ways in which adaptive evolution shapes gene expression.

Copy number variation introduced by a massive mobile element facilitates global thermal adaptation in a fungal wheat pathogen

Copy number variation (CNV) can drive rapid evolution in changing environments. In microbial pathogens, such adaptation is a key factor underpinning epidemics and colonization of new niches. However, the genomic determinants of such adaptation remain poorly understood. Here, we systematically investigate CNVs in a large genome sequencing dataset spanning a worldwide collection of 1104 genomes from the major wheat pathogen Zymoseptoria tritici. We found overall strong purifying selection acting on most CNVs. Genomic defense mechanisms likely accelerated gene loss over episodes of continental colonization. Local adaptation along climatic gradients was likely facilitated by CNVs affecting secondary metabolite production and gene loss in general. One of the strongest loci for climatic adaptation is a highly conserved gene of the NAD-dependent Sirtuin family. The Sirtuin CNV locus localizes to an ~68-kb Starship mobile element unique to the species carrying genes highly expressed during plant infection. The element has likely lost the ability to transpose, demonstrating how the ongoing domestication of cargo-carrying selfish elements can contribute to selectable variation within populations. Our work highlights how standing variation in gene copy numbers at the global scale can be a major factor driving climatic and metabolic adaptation in microbial species.

Aneuploidy Can Be an Evolutionary Diversion on the Path to Adaptation

Aneuploidy is common in eukaryotes, often leading to decreased fitness. However, evidence from fungi and human tumur cells suggests that specific aneuploidies can be beneficial under stressful conditions and facilitate adaptation. In a previous evolutionary experiment with yeast, populations evolving under heat stress became aneuploid, only to later revert to euploidy after beneficial mutations accumulated. It was therefore suggested that aneuploidy is a "stepping stone" on the path to adaptation. Here, we test this hypothesis. We use Bayesian inference to fit an evolutionary model with both aneuploidy and mutation to the experimental results. We then predict the genotype frequency dynamics during the experiment, demonstrating that most of the evolved euploid population likely did not descend from aneuploid cells, but rather from the euploid wild-type population. Our model shows how the beneficial mutation supply-the product of population size and beneficial mutation rate-determines the evolutionary dynamics: with low supply, much of the evolved population descends from aneuploid cells; but with high supply, beneficial mutations are generated fast enough to outcompete aneuploidy due to its inherent fitness cost. Our results suggest that despite its potential fitness benefits under stress, aneuploidy can be an evolutionary "diversion" rather than a "stepping stone": it can delay, rather than facilitate, the adaptation of the population, and cells that become aneuploid may leave less descendants compared to cells that remain diploid.

Evolution of haploid and diploid populations reveals common, strong, and variable pleiotropic effects in non-home environments

Adaptation is driven by the selection for beneficial mutations that provide a fitness advantage in the specific environment in which a population is evolving. However, environments are rarely constant or predictable. When an organism well adapted to one environment finds itself in another, pleiotropic effects of mutations that made it well adapted to its former environment will affect its success. To better understand such pleiotropic effects, we evolved both haploid and diploid barcoded budding yeast populations in multiple environments, isolated adaptive clones, and then determined the fitness effects of adaptive mutations in 'non-home' environments in which they were not selected. We find that pleiotropy is common, with most adaptive evolved lineages showing fitness effects in non-home environments. Consistent with other studies, we find that these pleiotropic effects are unpredictable: they are beneficial in some environments and deleterious in others. However, we do find that lineages with adaptive mutations in the same genes tend to show similar pleiotropic effects. We also find that ploidy influences the observed adaptive mutational spectra in a condition-specific fashion. In some conditions, haploids and diploids are selected with adaptive mutations in identical genes, while in others they accumulate mutations in almost completely disjoint sets of genes.

Copy number variation alters local and global mutational tolerance

Copy number variants (CNVs), duplications and deletions of genomic sequences, contribute to evolutionary adaptation but can also confer deleterious effects and cause disease. Whereas the effects of amplifying individual genes or whole chromosomes (i.e., aneuploidy) have been studied extensively, much less is known about the genetic and functional effects of CNVs of differing sizes and structures. Here, we investigated Saccharomyces cerevisiae (yeast) strains that acquired adaptive CNVs of variable structures and copy numbers following experimental evolution in glutamine-limited chemostats. Although beneficial in the selective environment, CNVs result in decreased fitness compared with the euploid ancestor in rich media. We used transposon mutagenesis to investigate mutational tolerance and genome-wide genetic interactions in CNV strains. We find that CNVs increase mutational target size, confer increased mutational tolerance in amplified essential genes, and result in novel genetic interactions with unlinked genes. We validated a novel genetic interaction between different CNVs and BMH1 that was common to multiple strains. We also analyzed global gene expression and found that transcriptional dosage compensation does not affect most genes amplified by CNVs, although gene-specific transcriptional dosage compensation does occur for ∼12% of amplified genes. Furthermore, we find that CNV strains do not show previously described transcriptional signatures of aneuploidy. Our study reveals the extent to which local and global mutational tolerance is modified by CNVs with implications for genome evolution and CNV-associated diseases, such as cancer.

Pleiotropy drives evolutionary repair of the responsiveness of polarized cell growth to environmental cues

The ability of cells to translate different extracellular cues into different intracellular responses is vital for their survival in unpredictable environments. In Saccharomyces cerevisiae, cell polarity is modulated in response to environmental signals which allows cells to adopt varying morphologies in different external conditions. The responsiveness of cell polarity to extracellular cues depends on the integration of the molecular network that regulates polarity establishment with networks that signal environmental changes. The coupling of molecular networks often leads to pleiotropic interactions that can make it difficult to determine whether the ability to respond to external signals emerges as an evolutionary response to environmental challenges or as a result of pleiotropic interactions between traits. Here, we study how the propensity of the polarity network of S. cerevisiae to evolve toward a state that is responsive to extracellular cues depends on the complexity of the environment. We show that the deletion of two genes, BEM3 and NRP1, disrupts the ability of the polarity network to respond to cues that signal the onset of the diauxic shift. By combining experimental evolution with whole-genome sequencing, we find that the restoration of the responsiveness to these cues correlates with mutations in genes involved in the sphingolipid synthesis pathway and that these mutations frequently settle in evolving populations irrespective of the complexity of the selective environment. We conclude that pleiotropic interactions make a significant contribution to the evolution of networks that are responsive to extracellular cues.

Best Practices in Microbial Experimental Evolution: Using Reporters and Long-Read Sequencing to Identify Copy Number Variation in Experimental Evolution

Copy number variants (CNVs), comprising gene amplifications and deletions, are a pervasive class of heritable variation. CNVs play a key role in rapid adaptation in both natural, and experimental, evolution. However, despite the advent of new DNA sequencing technologies, detection and quantification of CNVs in heterogeneous populations has remained challenging. Here, we summarize recent advances in the use of CNV reporters that provide a facile means of quantifying de novo CNVs at a specific locus in the genome, and nanopore sequencing, for resolving the often complex structures of CNVs. We provide guidance for the engineering and analysis of CNV reporters and practical guidelines for single-cell analysis of CNVs using flow cytometry. We summarize recent advances in nanopore sequencing, discuss the utility of this technology, and provide guidance for the bioinformatic analysis of these data to define the molecular structure of CNVs. The combination of reporter systems for tracking and isolating CNV lineages and long-read DNA sequencing for characterizing CNV structures enables unprecedented resolution of the mechanisms by which CNVs are generated and their evolutionary dynamics.

An improved algorithm for inferring mutational parameters from bar-seq evolution experiments

BACKGROUND

Genetic barcoding provides a high-throughput way to simultaneously track the frequencies of large numbers of competing and evolving microbial lineages. However making inferences about the nature of the evolution that is taking place remains a difficult task.

RESULTS

Here we describe an algorithm for the inference of fitness effects and establishment times of beneficial mutations from barcode sequencing data, which builds upon a Bayesian inference method by enforcing self-consistency between the population mean fitness and the individual effects of mutations within lineages. By testing our inference method on a simulation of 40,000 barcoded lineages evolving in serial batch culture, we find that this new method outperforms its predecessor, identifying more adaptive mutations and more accurately inferring their mutational parameters.

CONCLUSION

Our new algorithm is particularly suited to inference of mutational parameters when read depth is low. We have made Python code for our serial dilution evolution simulations, as well as both the old and new inference methods, available on GitHub ( https://github.com/FangfeiLi05/FitMut2 ), in the hope that it can find broader use by the microbial evolution community.

Paths to adaptation under fluctuating nitrogen starvation: The spectrum of adaptive mutations in Saccharomyces cerevisiae is shaped by retrotransposons and microhomology-mediated recombination

There are many mechanisms that give rise to genomic change: while point mutations are often emphasized in genomic analyses, evolution acts upon many other types of genetic changes that can result in less subtle perturbations. Changes in chromosome structure, DNA copy number, and novel transposon insertions all create large genomic changes, which can have correspondingly large impacts on phenotypes and fitness. In this study we investigate the spectrum of adaptive mutations that arise in a population under consistently fluctuating nitrogen conditions. We specifically contrast these adaptive alleles and the mutational mechanisms that create them, with mechanisms of adaptation under batch glucose limitation and constant selection in low, non-fluctuating nitrogen conditions to address if and how selection dynamics influence the molecular mechanisms of evolutionary adaptation. We observe that retrotransposon activity accounts for a substantial number of adaptive events, along with microhomology-mediated mechanisms of insertion, deletion, and gene conversion. In addition to loss of function alleles, which are often exploited in genetic screens, we identify putative gain of function alleles and alleles acting through as-of-yet unclear mechanisms. Taken together, our findings emphasize that how selection (fluctuating vs. non-fluctuating) is applied also shapes adaptation, just as the selective pressure (nitrogen vs. glucose) does itself. Fluctuating environments can activate different mutational mechanisms, shaping adaptive events accordingly. Experimental evolution, which allows a wider array of adaptive events to be assessed, is thus a complementary approach to both classical genetic screens and natural variation studies to characterize the genotype-to-phenotype-to-fitness map.

Synonymous and Nonsynonymous Substitutions in Dictyostelium discoideum Ammonium Transporter amtA Are Necessary for Functional Complementation in Saccharomyces cerevisiae

Ammonium transporters are present in all three domains of life. They have undergone extensive horizontal gene transfer (HGT), gene duplication, and functional diversification and therefore offer an excellent paradigm to study protein evolution. We attempted to complement a mep1Δmep2Δmep3Δ strain of Saccharomyces cerevisiae (triple-deletion strain), which otherwise cannot grow on ammonium as a sole nitrogen source at concentrations of <3 mM, with amtA of Dictyostelium discoideum, an orthologue of S. cerevisiae MEP2. We observed that amtA did not complement the triple-deletion strain of S. cerevisiae for growth on low-ammonium medium. We isolated two mutant derivatives of amtA (amtA M1 and amtA M2) from a PCR-generated mutant plasmid library that complemented the triple-deletion strain of S. cerevisiae. amtA M1 bears three nonsynonymous and two synonymous substitutions, which are necessary for its functionality. amtA M2 bears two nonsynonymous substitutions and one synonymous substitution, all of which are necessary for functionality. Interestingly, AmtA M1 transports ammonium but does not confer methylamine toxicity, while AmtA M2 transports ammonium and confers methylamine toxicity, demonstrating functional diversification. Preliminary biochemical analyses indicated that the mutants differ in their conformations as well as their mechanisms of ammonium transport. These intriguing results clearly point out that protein evolution cannot be fathomed by studying nonsynonymous and synonymous substitutions in isolation. The above-described observations have significant implications for various facets of biological processes and are discussed in detail. IMPORTANCE Functional diversification following gene duplication is one of the major driving forces of protein evolution. While the role of nonsynonymous substitutions in the functional diversification of proteins is well recognized, knowledge of the role of synonymous substitutions in protein evolution is in its infancy. Using functional complementation, we isolated two functional alleles of the D. discoideum ammonium transporter gene (amtA), which otherwise does not function in S. cerevisiae as an ammonium transporters. One of them is an ammonium transporter, while the other is an ammonium transporter that also confers methylammonium (ammonium analogue) toxicity, suggesting functional diversification. Surprisingly, both alleles require a combination of synonymous and nonsynonymous substitutions for their functionality. These results bring out a hitherto-unknown pathway of protein evolution and pave the way for not only understanding protein evolution but also interpreting single nucleotide polymorphisms (SNPs).

Long-Term Adaptation to Galactose as a Sole Carbon Source Selects for Mutations Outside the Canonical GAL Pathway

Galactose is a secondary fermentable sugar that requires specific regulatory and structural genes for its assimilation, which are under catabolite repression by glucose. When glucose is absent, the catabolic repression is attenuated, and the structural GAL genes are fully activated. In Saccharomyces cerevisiae, the GAL pathway is under selection in environments where galactose is present. However, it is unclear the adaptive strategies in response to long-term propagation in galactose as a sole carbon source in laboratory evolution experiments. Here, we performed a 4,000-generation evolution experiment using 48 diploid Saccharomyces cerevisiae populations to study adaptation in galactose. We show that fitness gains were greater in the galactose-evolved population than in identically evolved populations with glucose as a sole carbon source. Whole-genome sequencing of 96 evolved clones revealed recurrent de novo single nucleotide mutations in candidate targets of selection, copy number variations, and ploidy changes. We find that most mutations that improve fitness in galactose lie outside of the canonical GAL pathway. Reconstruction of specific evolved alleles in candidate target of selection, SEC23 and IRA1, showed a significant increase in fitness in galactose compared to glucose. In addition, most of our evolved populations (28/46; 61%) fixed aneuploidies on Chromosome VIII, suggesting a parallel adaptive amplification. Finally, we show greater loss of extrachromosomal elements in our glucose-evolved lineages compared with previous glucose evolution. Broadly, these data further our understanding of the evolutionary pressures that drive adaptation to less-preferred carbon sources.

Enforcement of Postzygotic Species Boundaries in the Fungal Kingdom

Understanding the molecular basis of speciation is a primary goal in evolutionary biology. The formation of the postzygotic reproductive isolation that causes hybrid dysfunction, thereby reducing gene flow between diverging populations, is crucial for speciation. Using various advanced approaches, including chromosome replacement, hybrid introgression and transcriptomics, population genomics, and experimental evolution, scientists have revealed multiple mechanisms involved in postzygotic barriers in the fungal kingdom. These results illuminate both unique and general features of fungal speciation. Our review summarizes experiments on fungi exploring how Dobzhansky-Muller incompatibility, killer meiotic drive, chromosome rearrangements, and antirecombination contribute to postzygotic reproductive isolation. We also discuss possible evolutionary forces underlying different reproductive isolation mechanisms and the potential roles of the evolutionary arms race under the Red Queen hypothesis and epigenetic divergence in speciation.

Complex Genomic Rearrangements following Selection in a Glutamine-Limited Medium over Hundreds of Generations

Large-scale genomic changes, including copy number variations (CNVs), are frequently observed in long-term evolution experiments (LTEEs). We have previously reported the detection of recurrent CNVs in Saccharomyces cerevisiae populations adapting to glutamine-limited conditions over hundreds of generations. Here, we present the whole-genome sequencing (WGS) assemblies of 7 LTEE strains and their ancestor.

Preliminary study on the molecular features of mutation in multiple primary oral cancer by whole exome sequencing

Multiple primary cancers (MPCs) refer to cancers that occur simultaneously or metachronously in the same individual. The incidence of MPC has increased recently, as the survival time of malignant tumor patients has been greatly prolonged. It is difficult to differentiate MPC from primary cancers (PCs) in the same anatomical region from the clinical manifestation alone. However, their biological behaviors appear to be distinct. In this study, we show that the prognosis of multiple primary oral cancers (MP-OCs) is worse than primary oral cancers (P-OCs). To better understand the molecular mechanisms of MP-OC, we used whole exome sequencing (WES) to analyze samples from 9 patients with MP-OC and 21 patients with P-OC. We found more somatic mutations in MP-OC than in P-OC. MP-OC had more complicated mutation signatures, which were associated with age-related and Apolipoprotein B mRNA Editing Catalytic Polypeptide-like (APOBEC) activity-related signatures. Tumor mutational burden (TMB) and mutant-allele tumor heterogeneity (MATH) of MP-OC trended higher compared to P-OC. KEGG and GO analysis showed the differential pathways of MP-OC versus P-OC. In addition, MP-OC took amplification, not loss, as the main pattern of copy number variation (CNV), while P-OC took both. Lastly, we did not find significantly different mutant germline genes, but MSH-6 mutation may be a potential MP-OC driver. In short, our preliminary results show that MP-OC and P-OC have different molecular characteristics.

Neural networks enable efficient and accurate simulation-based inference of evolutionary parameters from adaptation dynamics

The rate of adaptive evolution depends on the rate at which beneficial mutations are introduced into a population and the fitness effects of those mutations. The rate of beneficial mutations and their expected fitness effects is often difficult to empirically quantify. As these 2 parameters determine the pace of evolutionary change in a population, the dynamics of adaptive evolution may enable inference of their values. Copy number variants (CNVs) are a pervasive source of heritable variation that can facilitate rapid adaptive evolution. Previously, we developed a locus-specific fluorescent CNV reporter to quantify CNV dynamics in evolving populations maintained in nutrient-limiting conditions using chemostats. Here, we use CNV adaptation dynamics to estimate the rate at which beneficial CNVs are introduced through de novo mutation and their fitness effects using simulation-based likelihood-free inference approaches. We tested the suitability of 2 evolutionary models: a standard Wright-Fisher model and a chemostat model. We evaluated 2 likelihood-free inference algorithms: the well-established Approximate Bayesian Computation with Sequential Monte Carlo (ABC-SMC) algorithm, and the recently developed Neural Posterior Estimation (NPE) algorithm, which applies an artificial neural network to directly estimate the posterior distribution. By systematically evaluating the suitability of different inference methods and models, we show that NPE has several advantages over ABC-SMC and that a Wright-Fisher evolutionary model suffices in most cases. Using our validated inference framework, we estimate the CNV formation rate at the GAP1 locus in the yeast Saccharomyces cerevisiae to be 10-4.7 to 10-4 CNVs per cell division and a fitness coefficient of 0.04 to 0.1 per generation for GAP1 CNVs in glutamine-limited chemostats. We experimentally validated our inference-based estimates using 2 distinct experimental methods-barcode lineage tracking and pairwise fitness assays-which provide independent confirmation of the accuracy of our approach. Our results are consistent with a beneficial CNV supply rate that is 10-fold greater than the estimated rates of beneficial single-nucleotide mutations, explaining the outsized importance of CNVs in rapid adaptive evolution. More generally, our study demonstrates the utility of novel neural network-based likelihood-free inference methods for inferring the rates and effects of evolutionary processes from empirical data with possible applications ranging from tumor to viral evolution.

Emergence of Phenotypically Distinct Subpopulations Is a Factor in Adaptation of Recombinant Saccharomyces cerevisiae under Glucose-Limited Conditions

Cells cultured in a nutrient-limited environment can undergo adaptation, which confers improved fitness under long-term energy limitation. We have shown previously how a recombinant Saccharomyces cerevisiae strain, producing a heterologous insulin product, under glucose-limited conditions adapts over time at the average population level. Here, we investigated this adaptation at the single-cell level by application of fluorescence-activated cell sorting (FACS) and showed that the following three apparent phenotypes underlie the adaptive response observed at the bulk level: (i) cells that drastically reduced insulin production (23%), (ii) cells with reduced enzymatic capacity in central carbon metabolism (46%), and (iii) cells that exhibited pseudohyphal growth (31%). We speculate that the phenotypic heterogeneity is a result of different mechanisms to increase fitness. Cells with reduced insulin productivity have increased fitness by reducing the burden of the heterologous insulin production, and the populations with reduced enzymatic capacity of the central carbon metabolism and pseudohyphal growth have increased fitness toward the glucose-limited conditions. The results highlight the importance of considering population heterogeneity when studying adaptation and evolution. IMPORTANCE The yeast Saccharomyces cerevisiae is an attractive microbial host for industrial production and is used widely for manufacturing, e.g., pharmaceuticals. Chemostat cultivation mode is an efficient cultivation strategy for industrial production processes as it ensures a constant, well-controlled cultivation environment. Nevertheless, both the production of a heterologous product and the constant cultivation environment in the chemostat impose a selective pressure on the production organism, which may result in adaptation and loss of productivity. The exact mechanisms behind the observed adaptation and loss of performance are often unidentified. We used a recombinant S. cerevisiae strain producing heterologous insulin and investigated the adaptation occurring during chemostat growth at the single-cell level. We showed that three apparent phenotypes underlie the adaptive response observed at the bulk level in the chemostat. These findings highlight the importance of considering population heterogeneity when studying adaptation in industrial bioprocesses.

Fluctuating Environments Maintain Genetic Diversity through Neutral Fitness Effects and Balancing Selection

Genetic variation is the raw material upon which selection acts. The majority of environmental conditions change over time and therefore may result in variable selective effects. How temporally fluctuating environments impact the distribution of fitness effects and in turn population diversity is an unresolved question in evolutionary biology. Here, we employed continuous culturing using chemostats to establish environments that switch periodically between different nutrient limitations and compared the dynamics of selection to static conditions. We used the pooled Saccharomyces cerevisiae haploid gene deletion collection as a synthetic model for populations comprising thousands of unique genotypes. Using barcode sequencing, we find that static environments are uniquely characterized by a small number of high-fitness genotypes that rapidly dominate the population leading to dramatic decreases in genetic diversity. By contrast, fluctuating environments are enriched in genotypes with neutral fitness effects and an absence of extreme fitness genotypes contributing to the maintenance of genetic diversity. We also identified a unique class of genotypes whose frequencies oscillate sinusoidally with a period matching the environmental fluctuation. Oscillatory behavior corresponds to large differences in short-term fitness that are not observed across long timescales pointing to the importance of balancing selection in maintaining genetic diversity in fluctuating environments. Our results are consistent with a high degree of environmental specificity in the distribution of fitness effects and the combined effects of reduced and balancing selection in maintaining genetic diversity in the presence of variable selection.

Changes in the distribution of fitness effects and adaptive mutational spectra following a single first step towards adaptation

Historical contingency and diminishing returns epistasis have been typically studied for relatively divergent genotypes and/or over long evolutionary timescales. Here, we use Saccharomyces cerevisiae to study the extent of diminishing returns and the changes in the adaptive mutational spectra following a single first adaptive mutational step. We further evolve three clones that arose under identical conditions from a common ancestor. We follow their evolutionary dynamics by lineage tracking and determine adaptive outcomes using fitness assays and whole genome sequencing. We find that diminishing returns manifests as smaller fitness gains during the 2^nd step of adaptation compared to the 1^st step, mainly due to a compressed distribution of fitness effects. We also find that the beneficial mutational spectra for the 2^nd adaptive step are contingent on the 1^st step, as we see both shared and diverging adaptive strategies. Finally, we find that adaptive loss-of-function mutations, such as nonsense and frameshift mutations, are less common in the second step of adaptation than in the first step.

Analyses of both natural and experimental evolution suggest that adaptation depends on the evolutionary past and adaptive potential decreases over time. Here, by tracking yeast adaptation with DNA barcoding, the authors show that such evolutionary phenomena can be observed even after a single adaptive step.

The genetic basis of differential autodiploidization in evolving yeast populations

Spontaneous whole-genome duplication, or autodiploidization, is a common route to adaptation in experimental evolution of haploid budding yeast populations. The rate at which autodiploids fix in these populations appears to vary across strain backgrounds, but the genetic basis of these differences remains poorly characterized. Here, we show that the frequency of autodiploidization differs dramatically between two closely related laboratory strains of Saccharomyces cerevisiae, BY4741 and W303. To investigate the genetic basis of this difference, we crossed these strains to generate hundreds of unique F1 segregants and tested the tendency of each segregant to autodiplodize across hundreds of generations of laboratory evolution. We find that variants in the SSD1 gene are the primary genetic determinant of differences in autodiploidization. We then used multiple laboratory and wild strains of S. cerevisiae to show that clonal populations of strains with a functional copy of SSD1 autodiploidize more frequently in evolution experiments, while knocking out this gene or replacing it with the W303 allele reduces autodiploidization propensity across all genetic backgrounds tested. These results suggest a potential strategy for modifying rates of spontaneous whole-genome duplications in laboratory evolution experiments in haploid budding yeast. They may also have relevance to other settings in which eukaryotic genome stability plays an important role, such as biomanufacturing and the treatment of pathogenic fungal diseases and cancers.

Consequences of Cryopreservation in Diverse Natural Isolates of Saccharomyces cerevisiae

Experimental evolution allows the observation of change over time as laboratory populations evolve in response to novel, controlled environments. Microbial evolution experiments take advantage of cryopreservation to archive experimental populations in glycerol media, creating a frozen, living "fossil" record. Prior research with Escherichia coli has shown that cryopreservation conditions can affect cell viability and that allele frequencies across the genome can change in response to a freeze-thaw event. We expand on these observations by characterizing fitness and genomic consequences of multiple freeze-thaw cycles in diploid yeast populations. Our study system is a highly recombinant Saccharomyces cerevisiae population (SGRP-4X) that harbors standing genetic variation that cryopreservation may threaten. We also investigate the four parental isogenic strains crossed to create the SGRP-4X. We measure cell viability over five consecutive freeze-thaw cycles; whereas we find that viability increases over time in the evolved recombinant populations, we observe no such viability improvements in the parental strains. We also collect genome-wide sequence data from experimental populations initially, after one freeze-thaw, and after five freeze-thaw cycles. In the recombinant evolved populations, we find a region of significant allele frequency change on chromosome 15 containing the ALR1 gene. In the parental strains, we find little evidence for new mutations. We conclude that cryopreserving yeast populations with standing genetic variation may have both phenotypic and genomic consequences, though the same cryopreservation practices may have only small impacts on populations with little or no initial variation.

Adaptive evolution of sulfite tolerance in Brettanomyces bruxellensis

Brettanomyces bruxellensis is considered one of the most problematic microbes associated with wine production. Sulfur dioxide is commonly used to inhibit the growth of B. bruxellensis and limit the potential wine spoilage. Brettanomyces bruxellensis wine isolates can grow at higher concentrations of this preservative than isolates from other sources. Thus, it has been suggested that the use of sulfite may have selected for B. bruxellensis strains better adapted to survive in the winemaking environment. We utilized laboratory adaptive evolution to determine the potential for this to occur. Three B. bruxellensis strains, representative of known genetic variation within the species, were subjected to increasing sublethal sulfur dioxide concentrations. Individual clones isolated from evolved populations displayed enhanced sulfite tolerance, ranging from 1.6 to 2.5 times higher than the corresponding parental strains. Whole-genome sequencing of sulfite-tolerant clones derived from two of the parental strains revealed structural variations affecting 270 genes. The region containing the sulfite efflux pump encoding gene, SSU1, showed clear copy number variants in all sequenced clones. Regardless of parental strain genetic background, SSU1 copy number changes were reproducibly associated with one SSU1 haplotype. This work clearly demonstrates adaptive evolution of B. bruxellensis when exposed to sublethal sulfites and suggests that, similar to Saccharomyces cerevisiae wine yeast, the mechanism responsible involves the gene SSU1.

Genetic interaction network has a very limited impact on the evolutionary trajectories in continuous culture-grown populations of yeast

Background

The impact of genetic interaction networks on evolution is a fundamental issue. Previous studies have demonstrated that the topology of the network is determined by the properties of the cellular machinery. Functionally related genes frequently interact with one another, and they establish modules, e.g., modules of protein complexes and biochemical pathways. In this study, we experimentally tested the hypothesis that compensatory evolutionary modifications, such as mutations and transcriptional changes, occur frequently in genes from perturbed modules of interacting genes.

Results

Using Saccharomyces cerevisiae haploid deletion mutants as a model, we investigated two modules lacking COG7 or NUP133, which are evolutionarily conserved genes with many interactions. We performed laboratory evolution experiments with these strains in two genetic backgrounds (with or without additional deletion of MSH2), subjecting them to continuous culture in a non-limiting minimal medium. Next, the evolved yeast populations were characterized through whole-genome sequencing and transcriptome analyses. No obvious compensatory changes resulting from inactivation of genes already included in modules were identified. The supposedly compensatory inactivation of genes in the evolved strains was only rarely observed to be in accordance with the established fitness effect of the genetic interaction network. In fact, a substantial majority of the gene inactivations were predicted to be neutral in the experimental conditions used to determine the interaction network. Similarly, transcriptome changes during continuous culture mostly signified adaptation to growth conditions rather than compensation of the absence of the COG7, NUP133 or MSH2 genes. However, we noticed that for genes whose inactivation was deleterious an upregulation of transcription was more common than downregulation.

Conclusions

Our findings demonstrate that the genetic interactions and the modular structure of the network described by others have very limited effects on the evolutionary trajectory following gene deletion of module elements in our experimental conditions and has no significant impact on short-term compensatory evolution. However, we observed likely compensatory evolution in functionally related (albeit non-interacting) genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01830-9.

The Role of Structural Variation in Adaptation and Evolution of Yeast and Other Fungi

Mutations in DNA can be limited to one or a few nucleotides, or encompass larger deletions, insertions, duplications, inversions and translocations that span long stretches of DNA or even full chromosomes. These so-called structural variations (SVs) can alter the gene copy number, modify open reading frames, change regulatory sequences or chromatin structure and thus result in major phenotypic changes. As some of the best-known examples of SV are linked to severe genetic disorders, this type of mutation has traditionally been regarded as negative and of little importance for adaptive evolution. However, the advent of genomic technologies uncovered the ubiquity of SVs even in healthy organisms. Moreover, experimental evolution studies suggest that SV is an important driver of evolution and adaptation to new environments. Here, we provide an overview of the causes and consequences of SV and their role in adaptation, with specific emphasis on fungi since these have proven to be excellent models to study SV.

Thermo-adaptive evolution to generate improved Saccharomyces cerevisiae strains for cocoa pulp fermentations

Cocoa pulp fermentation is a consequence of the succession of indigenous yeasts, lactic acid bacteria and acetic acid bacteria that not only produce a diversity of metabolites, but also cause the production of flavour precursors. However, as such spontaneous fermentations are less reproducible and contribute to produce variability, interest in a microbial starter culture is growing that could be used to inoculate cocoa pulp fermentations. This study aimed to generate robust S. cerevisiae strains by thermo-adaptive evolution that could be used in cocoa fermentation. We evolved a cocoa strain in a sugary defined medium at high temperature to improve both fermentation and growth capacity. Moreover, adaptive evolution at high temperature (40 °C) also enabled us to unveil the molecular basis underlying the improved phenotype by analysing the whole genome sequence of the evolved strain. Adaptation to high-temperature conditions occurred at different genomic levels, and promoted aneuploidies, segmental duplication, and SNVs in the evolved strain. The lipid profile analysis of the evolved strain also evidenced changes in the membrane composition that contribute to maintain an appropriate cell membrane state at high temperature. Our work demonstrates that experimental evolution is an effective approach to generate better-adapted yeast strains at high temperature for industrial processes.

The population genetics of ploidy change in unicellular fungi

Changes in ploidy are a significant type of genetic variation, describing the number of chromosome sets per cell. Ploidy evolves in natural populations, clinical populations, and lab experiments, particularly in fungi. Despite a long history of theoretical work on this topic, predicting how ploidy will evolve has proven difficult, as it is often unclear why one ploidy state outperforms another. Here, we review what is known about contemporary ploidy evolution in diverse fungal species through the lens of population genetics. As with typical genetic variants, ploidy evolution depends on the rate that new ploidy states arise by mutation, natural selection on alternative ploidy states, and random genetic drift. However, ploidy variation also has unique impacts on evolution, with the potential to alter chromosomal stability, the rate and patterns of point mutation, and the nature of selection on all loci in the genome. We discuss how ploidy evolution depends on these general and unique factors and highlight areas where additional experimental evidence is required to comprehensively explain the ploidy transitions observed in the field and the lab.

Acetate and glycerol are not uniquely suited for the evolution of cross-feeding in E. coli

The evolution of cross-feeding among individuals of the same species can help generate genetic and phenotypic diversity even in completely homogeneous environments. Cross-feeding Escherichia coli strains, where one strain feeds on a carbon source excreted by another strain, rapidly emerge during experimental evolution in a chemically minimal environment containing glucose as the sole carbon source. Genome-scale metabolic modeling predicts that cross-feeding of 58 carbon sources can emerge in the same environment, but only cross-feeding of acetate and glycerol has been experimentally observed. Here we use metabolic modeling to ask whether acetate and glycerol cross-feeding are especially likely to evolve, perhaps because they require less metabolic change, and thus perhaps also less genetic change than other cross-feeding interactions. However, this is not the case. The minimally required metabolic changes required for acetate and glycerol cross feeding affect dozens of chemical reactions, multiple biochemical pathways, as well as multiple operons or regulons. The complexity of these changes is consistent with experimental observations, where cross-feeding strains harbor multiple mutations. The required metabolic changes are also no less complex than those observed for multiple other of the 56 cross feeding interactions we study. We discuss possible reasons why only two cross-feeding interactions have been discovered during experimental evolution and argue that multiple new cross-feeding interactions may await discovery.

The evolution of cross-feeding interactions, where one organism thrives by consuming the excretions of others, can create diversity even in simple and homogeneous environments. In past work, we had predicted that 58 cross-feeding interactions could evolve in populations of E. coli grown in glucose minimal media, yet only two have been experimentally observed, those involving acetate and glycerol. We hypothesized that multiple mutations might be required for the evolution of computationally predicted but not experimentally observed cross-feeding interactions. To answer this question, we developed a method that searches for the minimal number of metabolic changes required for individuals to change their metabolic state (from an ancestral glucose-consuming state to an evolved state that produces or consumes other metabolite). We observed that the metabolic changes required for the evolution of acetate and glycerol cross-feeding are no less complex than those required for the evolution of the other predicted cross-feeding interactions, which suggests that multiple cross-feeding interactions may still await discovery.

Slow Adaptive Response of Budding Yeast Cells to Stable Conditions of Continuous Culture Can Occur without Genome Modifications

Continuous cultures assure the invariability of environmental conditions and the metabolic state of cultured microorganisms, whereas batch-cultured cells undergo constant changes in nutrients availability. For that reason, continuous culture is sometimes employed in the whole transcriptome, whole proteome, or whole metabolome studies. However, the typical method for establishing uniform growth of a cell population, i.e., by limited chemostat, results in the enrichment of the cell population gene pool with mutations adaptive for starvation conditions. These adaptive changes can skew the results of large-scale studies. It is commonly assumed that these adaptations reflect changes in the genome, and this assumption has been confirmed experimentally in rare cases. Here we show that in a population of budding yeast cells grown for over 200 generations in continuous culture in non-limiting minimal medium and therefore not subject to selection pressure, remodeling of transcriptome occurs, but not as a result of the accumulation of adaptive mutations. The observed changes indicate a shift in the metabolic balance towards catabolism, a decrease in ribosome biogenesis, a decrease in general stress alertness, reorganization of the cell wall, and transactions occurring at the cell periphery. These adaptive changes signify the acquisition of a new lifestyle in a stable nonstressful environment. The absence of underlying adaptive mutations suggests these changes may be regulated by another mechanism.

Evolutionary stalling and a limit on the power of natural selection to improve a cellular module

Cells consist of molecular modules which perform vital biological functions. Cellular modules are key units of adaptive evolution because organismal fitness depends on their performance. Theory shows that in rapidly evolving populations, such as those of many microbes, adaptation is driven primarily by common beneficial mutations with large effects, while other mutations behave as if they are effectively neutral. As a consequence, if a module can be improved only by rare and/or weak beneficial mutations, its adaptive evolution would stall. However, such evolutionary stalling has not been empirically demonstrated, and it is unclear to what extent stalling may limit the power of natural selection to improve modules. Here we empirically characterize how natural selection improves the translation machinery (TM), an essential cellular module. We experimentally evolved populations of Escherichia coli with genetically perturbed TMs for 1,000 generations. Populations with severe TM defects initially adapted via mutations in the TM, but TM adaptation stalled within about 300 generations. We estimate that the genetic load in our populations incurred by residual TM defects ranges from 0.5 to 19%. Finally, we found evidence that both epistasis and the depletion of the pool of beneficial mutations contributed to evolutionary stalling. Our results suggest that cellular modules may not be fully optimized by natural selection despite the availability of adaptive mutations.

Ethanol exposure increases mutation rate through error-prone polymerases

Ethanol is a ubiquitous environmental stressor that is toxic to all lifeforms. Here, we use the model eukaryote Saccharomyces cerevisiae to show that exposure to sublethal ethanol concentrations causes DNA replication stress and an increased mutation rate. Specifically, we find that ethanol slows down replication and affects localization of Mrc1, a conserved protein that helps stabilize the replisome. In addition, ethanol exposure also results in the recruitment of error-prone DNA polymerases to the replication fork. Interestingly, preventing this recruitment through mutagenesis of the PCNA/Pol30 polymerase clamp or deleting specific error-prone polymerases abolishes the mutagenic effect of ethanol. Taken together, this suggests that the mutagenic effect depends on a complex mechanism, where dysfunctional replication forks lead to recruitment of error-prone polymerases. Apart from providing a general mechanistic framework for the mutagenic effect of ethanol, our findings may also provide a route to better understand and prevent ethanol-associated carcinogenesis in higher eukaryotes.

Whereas the toxic effects of ethanol are well-documented, the underlying mechanism is obscure. This study uses the eukaryotic model S. cerevisiae to reveal how exposure to sublethal ethanol concentrations causes DNA replication stress and an increased mutation rate.

Lager-brewing yeasts in the era of modern genetics

The yeast Saccharomyces pastorianus is responsible for the annual worldwide production of almost 200 billion liters of lager-type beer. S. pastorianus is a hybrid of Saccharomyces cerevisiae and Saccharomyces eubayanus that has been studied for well over a century. Scientific interest in S. pastorianus intensified upon the discovery, in 2011, of its S. eubayanus ancestor. Moreover, advances in whole-genome sequencing and genome editing now enable deeper exploration of the complex hybrid and aneuploid genome architectures of S. pastorianus strains. These developments not only provide novel insights into the emergence and domestication of S. pastorianus but also generate new opportunities for its industrial application. This review paper combines historical, technical and socioeconomic perspectives to analyze the evolutionary origin and genetics of S. pastorianus. In addition, it provides an overview of available methods for industrial strain improvement and an outlook on future industrial application of lager-brewing yeasts. Particular attention is given to the ongoing debate on whether current S. pastorianus originates from a single or multiple hybridization events and to the potential role of genome editing in developing industrial brewing yeast strains.

Differential Gene Expression and Allele Frequency Changes Favour Adaptation of a Heterogeneous Yeast Population to Nitrogen-Limited Fermentations

Alcoholic fermentation is fundamentally an adaptation process, in which the yeast Saccharomyces cerevisiae outperforms its competitors and takes over the fermentation process itself. Although wine yeast strains appear to be adapted to the stressful conditions of alcoholic fermentation, nitrogen limitations in grape must cause stuck or slow fermentations, generating significant economic losses for the wine industry. One way to discover the genetic bases that promote yeast adaptation to nitrogen-deficient environments are selection experiments, where a yeast population undergoes selection under conditions of nitrogen restriction for a number of generations, to then identify by sequencing the molecular characteristics that promote this adaptation. In this work, we carried out selection experiments in bioreactors imitating wine fermentation under nitrogen-limited fermentation conditions (SM60), using the heterogeneous SGRP-4X yeast population, to then sequence the transcriptome and the genome of the population at different time points of the selection process. The transcriptomic results showed an overexpression of genes from the NA strain (North American/YPS128), a wild, non-domesticated isolate. In addition, genome sequencing and allele frequency results allowed several QTLs to be mapped for adaptation to nitrogen-limited fermentation. Finally, we validated the ECM38 allele of NA strain as responsible for higher growth efficiency under nitrogen-limited conditions. Taken together, our results revealed a complex pattern of molecular signatures favouring adaptation of the yeast population to nitrogen-limited fermentations, including differential gene expression, allele frequency changes and loss of the mitochondrial genome. Finally, the results suggest that wild alleles from a non-domesticated isolate (NA) may have a relevant role in the adaptation to the assayed fermentation conditions, with the consequent potential of these alleles for the genetic improvement of wine yeast strains.

Aneuploidy in yeast: Segregation error or adaptation mechanism?

Aneuploidy is the loss or gain of chromosomes within a genome. It is often detrimental and has been associated with cell death and genetic disorders. However, aneuploidy can also be beneficial and provide a quick solution through changes in gene dosage when cells face environmental stress. Here, we review the prevalence of aneuploidy in Saccharomyces, Candida, and Cryptococcus yeasts (and their hybrid offspring) and analyse associations with chromosome size and specific stressors. We discuss how aneuploidy, a segregation error, may in fact provide a natural route for the diversification of microbes and enable important evolutionary innovations given the right ecological circumstances, such as the colonisation of new environments or the transition from commensal to pathogenic lifestyle. We also draw attention to a largely unstudied cross link between hybridisation and aneuploidy. Hybrid meiosis, involving two divergent genomes, can lead to drastically increased rates of aneuploidy in the offspring due to antirecombination and chromosomal missegregation. Because hybridisation and aneuploidy have both been shown to increase with environmental stress, we believe it important and timely to start exploring the evolutionary significance of their co-occurrence.

Genetic variants of TORC1 signaling pathway affect nitrogen consumption in Saccharomyces cerevisiae during alcoholic fermentation

In the alcoholic fermentation process, Saccharomyces cerevisiae strains present differences in their nitrogen consumption profiles, these phenotypic outcomes have complex genetic and molecular architectures. In this sense, variations in nitrogen signaling pathways regulated by TORC1 represent one of the main sources of phenotypic diversity in nitrogen consumption. This emphasizes the possible roles that allelic variants from the TORC1 pathway have in the nitrogen consumption differences observed in yeast during the alcoholic fermentation. Here, we studied the allelic diversity in the TORC1 pathway across four yeast strains and determined how these polymorphisms directly impact nitrogen consumption during alcoholic fermentation. Using a reciprocal hemizygosity approach combined with phenotyping under fermentative conditions, we found that allelic variants of GTR1, TOR2, SIT4, SAP185, EAP1, NPR1 and SCH9 underlie differences in the ammonium and amino acids consumption phenotypes. Among these, GTR1 alleles from the Wine/European and West African genetic backgrounds showed the greatest effects on ammonium and amino acid consumption, respectively. Furthermore, we identified allelic variants of SAP185, TOR2, SCH9 and NPR1 from an oak isolate that increased the amino acid consumption preference over ammonium; representing putative candidates coming from a non-domesticated strain that could be used for genetic improvement programs. In conclusion, our results demonstrated that a large number of allelic variants within the TORC1 pathway significantly impacts on regulatory mechanisms of nitrogen assimilation during alcoholic fermentation.

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.

Adaptive Evolution of Escherichia coli to Ciprofloxacin in Controlled Stress Environments: Contrasting Patterns of Resistance in Spatially Varying versus Uniformly Mixed Concentration Conditions

A microfluidic gradient chamber (MGC) and a homogeneous batch culturing system were used to evaluate whether spatial concentration gradients of the antibiotic ciprofloxacin allow development of greater antibiotic resistance in Escherichia coli strain 307 (E. coli 307) compared to exclusively temporal concentration gradients, as indicated in an earlier study. A linear spatial gradient of ciprofloxacin and Luria-Bertani broth (LB) medium was established and maintained by diffusion over 5 days across a well array in the MGC, with relative concentrations along the gradient of 1.7-7.7× the original minimum inhibitory concentration (MIC_original). The E. coli biomass increased in wells with lower ciprofloxacin concentrations, and only a low level of resistance to ciprofloxacin was detected in the recovered cells (∼2× MIC_original). Homogeneous batch culture experiments were performed with the same temporal exposure history to ciprofloxacin concentration, the same and higher initial cell densities, and the same and higher nutrient (i.e., LB) concentrations as in the MGC. In all batch experiments, E. coli 307 developed higher ciprofloxacin resistance after exposure, ranging from 4 to 24× MIC_original in all replicates. Hence, these results suggest that the presence of spatial gradients appears to reduce the driving force for E. coli 307 adaptation to ciprofloxacin, which suggests that results from batch experiments may over predict the development of antibiotic resistance in natural environments.

An evolving view of copy number variants

Copy number variants (CNVs) are regions of the genome that vary in integer copy number. CNVs, which comprise both amplifications and deletions of DNA sequence, have been identified across all domains of life, from bacteria and archaea to plants and animals. CNVs are an important source of genetic diversity, and can drive rapid adaptive evolution and progression of heritable and somatic human diseases, such as cancer. However, despite their evolutionary importance and clinical relevance, CNVs remain understudied compared to single-nucleotide variants (SNVs). This is a consequence of the inherent difficulties in detecting CNVs at low-to-intermediate frequencies in heterogeneous populations of cells. Here, we discuss molecular methods used to detect CNVs, the limitations associated with using these techniques, and the application of new and emerging technologies that present solutions to these challenges. The goal of this short review and perspective is to highlight aspects of CNV biology that are understudied and define avenues for further research that address specific gaps in our knowledge of these complex alleles. We describe our recently developed method for CNV detection in which a fluorescent gene functions as a single-cell CNV reporter and present key findings from our evolution experiments in Saccharomyces cerevisiae. Using a CNV reporter, we found that CNVs are generated at a high rate and undergo selection with predictable dynamics across independently evolving replicate populations. Many CNVs appear to be generated through DNA replication-based processes that are mediated by the presence of short, interrupted, inverted-repeat sequences. Our results have important implications for the role of CNVs in evolutionary processes and the molecular mechanisms that underlie CNV formation. We discuss the possible extension of our method to other applications, including tracking the dynamics of CNVs in models of human tumors.

Evolutionary Transition of GAL Regulatory Circuit from Generalist to Specialist Function in Ascomycetes

The Gal4 transcription factor (TF) controls gene expression by binding the DNA sequence motif CGG(N₁₁)CCG. Well studied versions regulate metabolism of glucose in Candida albicans and galactose in Saccharomyces cerevisiae. Gal4 is also found within Aspergillus species and shows a wide range of potential binding targets. Members of the CTG clade that reassigned CUG codons from leucine to serine lack the Gal80 binding domain of Gal4, and they use the TF to regulate only glycolytic genes. In this clade, the galactose catabolic pathway (also known as the Leloir pathway) genes are regulated by Rtg1/Rtg3. In the WGD species, the complete Gal4/Gal80 module is limited to regulation of the Leloir pathway, while glycolysis is controlled by Gcr1/Gcr2. This shows a switch of Gal4 from a generalist to a specialist within the ascomycetes, and the split of glucose and galactose metabolism into distinct regulatory circuits.

Rewired cellular signaling coordinates sugar and hypoxic responses for anaerobic xylose fermentation in yeast

Microbes can be metabolically engineered to produce biofuels and biochemicals, but rerouting metabolic flux toward products is a major hurdle without a systems-level understanding of how cellular flux is controlled. To understand flux rerouting, we investigated a panel of Saccharomyces cerevisiae strains with progressive improvements in anaerobic fermentation of xylose, a sugar abundant in sustainable plant biomass used for biofuel production. We combined comparative transcriptomics, proteomics, and phosphoproteomics with network analysis to understand the physiology of improved anaerobic xylose fermentation. Our results show that upstream regulatory changes produce a suite of physiological effects that collectively impact the phenotype. Evolved strains show an unusual co-activation of Protein Kinase A (PKA) and Snf1, thus combining responses seen during feast on glucose and famine on non-preferred sugars. Surprisingly, these regulatory changes were required to mount the hypoxic response when cells were grown on xylose, revealing a previously unknown connection between sugar source and anaerobic response. Network analysis identified several downstream transcription factors that play a significant, but on their own minor, role in anaerobic xylose fermentation, consistent with the combinatorial effects of small-impact changes. We also discovered that different routes of PKA activation produce distinct phenotypes: deletion of the RAS/PKA inhibitor IRA2 promotes xylose growth and metabolism, whereas deletion of PKA inhibitor BCY1 decouples growth from metabolism to enable robust fermentation without division. Comparing phosphoproteomic changes across ira2Δ and bcy1Δ strains implicated regulatory changes linked to xylose-dependent growth versus metabolism. Together, our results present a picture of the metabolic logic behind anaerobic xylose flux and suggest that widespread cellular remodeling, rather than individual metabolic changes, is an important goal for metabolic engineering.

Whi2: a new player in amino acid sensing

A critical function of human, yeast, and bacterial cells is the ability to sense and respond to available nutrients such as glucose and amino acids. Cells must also detect declining nutrient levels to adequately prepare for starvation conditions by inhibiting cell growth and activating autophagy. The evolutionarily conserved protein complex TORC1 regulates these cellular responses to nutrients, and in particular to amino acid availability. Recently, we found that yeast Whi2 (Saccharomyces cerevisiae) and a human counterpart, KCTD11, that shares a conserved BTB structural domain, are required to suppress TORC1 activity under low amino acid conditions. Using yeast, the mechanisms were more readily dissected. Unexpectedly, Whi2 suppresses TORC1 activity independently of the well-known SEACIT-GTR pathway, analogous to the GATOR1-RAG pathway in mammals. Instead, Whi2 requires the plasma membrane-associated phosphatases Psr1 and Psr2, which were known to bind Whi2, although their role was unknown. Yeast WHI2 was previously reported to be involved in regulating several fundamental cellular processes including cell cycle arrest, general stress responses, the Ras-cAMP-PKA pathway, autophagy, and mitophagy, and to be frequently mutated in the yeast knockout collections and in genome evolution studies. Most of these observations are likely explained by the ability of Whi2 to inhibit TORC1. Thus, understanding the function of yeast Whi2 will provide deeper insights into the disease-related KCTD family proteins and the pathogenesis of plant and human fungal infections.

Experimental evolution: its principles and applications in developing stress-tolerant yeasts

Stress tolerance and resistance in industrial yeast strains are important attributes for cost-effective bioprocessing. The source of stress-tolerant yeasts ranges from extremophilic environments to laboratory engineered strains. However, industrial stress-tolerant yeasts are very rare in nature as the natural environment forces them to evolve traits that optimize survival and reproduction and not the ability to withstand harsh habitat-irrelevant industrial conditions. Experimental evolution is a frequent method used to uncover the mechanisms of evolution and microbial adaption towards environmental stresses. It optimizes biological systems by means of adaptation to environmental stresses and thus has immense power of development of robust stress-tolerant yeasts. This mini-review briefly outlines the basics and implications of evolution experiments and their applications to industrial biotechnology. This work is meant to serve as an introduction to those new to the field of experimental evolution, and as a guide to biologists working in the field of yeast stress response. Future perspectives of experimental evolution for potential biotechnological applications have also been elucidated.

Single-cell copy number variant detection reveals the dynamics and diversity of adaptation

Copy number variants (CNVs) are a pervasive source of genetic variation and evolutionary potential, but the dynamics and diversity of CNVs within evolving populations remain unclear. Long-term evolution experiments in chemostats provide an ideal system for studying the molecular processes underlying CNV formation and the temporal dynamics with which they are generated, selected, and maintained. Here, we developed a fluorescent CNV reporter to detect de novo gene amplifications and deletions in individual cells. We used the CNV reporter in Saccharomyces cerevisiae to study CNV formation at the GAP1 locus, which encodes the general amino acid permease, in different nutrient-limited chemostat conditions. We find that under strong selection, GAP1 CNVs are repeatedly generated and selected during the early stages of adaptive evolution, resulting in predictable dynamics. Molecular characterization of CNV-containing lineages shows that the CNV reporter detects different classes of CNVs, including aneuploidies, nonreciprocal translocations, tandem duplications, and complex CNVs. Despite GAP1's proximity to repeat sequences that facilitate intrachromosomal recombination, breakpoint analysis revealed that short inverted repeat sequences mediate formation of at least 50% of GAP1 CNVs. Inverted repeat sequences are also found at breakpoints at the DUR3 locus, where CNVs are selected in urea-limited chemostats. Analysis of 28 CNV breakpoints indicates that inverted repeats are typically 8 nucleotides in length and separated by 40 bases. The features of these CNVs are consistent with origin-dependent inverted-repeat amplification (ODIRA), suggesting that replication-based mechanisms of CNV formation may be a common source of gene amplification. We combined the CNV reporter with barcode lineage tracking and found that 102-104 independent CNV-containing lineages initially compete within populations, resulting in extreme clonal interference. However, only a small number (18-21) of CNV lineages ever constitute more than 1% of the CNV subpopulation, and as selection progresses, the diversity of CNV lineages declines. Our study introduces a novel means of studying CNVs in heterogeneous cell populations and provides insight into their dynamics, diversity, and formation mechanisms in the context of adaptive evolution.

Divergent Roles for cAMP-PKA Signaling in the Regulation of Filamentous Growth in Saccharomyces cerevisiae and Saccharomyces bayanus

The cyclic AMP - Protein Kinase A (cAMP-PKA) pathway is an evolutionarily conserved eukaryotic signaling network that is essential for growth and development. In the fungi, cAMP-PKA signaling plays a critical role in regulating cellular physiology and morphological switches in response to nutrient availability. We undertook a comparative investigation of the role that cAMP-PKA signaling plays in the regulation of filamentous growth in two closely related budding yeast species, Saccharomyces cerevisiae and Saccharomyces bayanus Using chemical and genetic perturbations of this pathway and its downstream targets we discovered divergent roles for cAMP-PKA signaling in the regulation of filamentous growth. While cAMP-PKA signaling is required for the filamentous growth response in both species, increasing or decreasing the activity of this pathway leads to drastically different phenotypic outcomes. In S. cerevisiae, cAMP-PKA inhibition ameliorates the filamentous growth response while hyper-activation of the pathway leads to increased filamentous growth; the same perturbations in S. bayanus result in the obverse. Divergence in the regulation of filamentous growth between S. cerevisiae and S. bayanus extends to downstream targets of PKA, including several kinases, transcription factors, and effector proteins. Our findings highlight the potential for significant evolutionary divergence in gene network function, even when the constituent parts of such networks are well conserved.

Experimental Design, Population Dynamics, and Diversity in Microbial Experimental Evolution

In experimental evolution, laboratory-controlled conditions select for the adaptation of species, which can be monitored in real time. Despite the current popularity of such experiments, nature's most pervasive biological force was long believed to be observable only on time scales that transcend a researcher's life-span, and studying evolution by natural selection was therefore carried out solely by comparative means. Eventually, microorganisms' propensity for fast evolutionary changes proved us wrong, displaying strong evolutionary adaptations over a limited time, nowadays massively exploited in laboratory evolution experiments. Here, we formulate a guide to experimental evolution with microorganisms, explaining experimental design and discussing evolutionary dynamics and outcomes and how it is used to assess ecoevolutionary theories, improve industrially important traits, and untangle complex phenotypes. Specifically, we give a comprehensive overview of the setups used in experimental evolution. Additionally, we address population dynamics and genetic or phenotypic diversity during evolution experiments and expand upon contributing factors, such as epistasis and the consequences of (a)sexual reproduction. Dynamics and outcomes of evolution are most profoundly affected by the spatiotemporal nature of the selective environment, where changing environments might lead to generalists and structured environments could foster diversity, aided by, for example, clonal interference and negative frequency-dependent selection. We conclude with future perspectives, with an emphasis on possibilities offered by fast-paced technological progress. This work is meant to serve as an introduction to those new to the field of experimental evolution, as a guide to the budding experimentalist, and as a reference work to the seasoned expert.