The repertoire and dynamics of evolutionary adaptations to controlled nutrient-limited environments in yeast

Transcriptomic and proteomic ramifications of segmental amplification in Escherichia coli

Gene amplification can drive adaptation by rapidly increasing the cellular dosage of critical gene products. Segmental amplifications often encompass large genomic regions surrounding the gene(s) under selection for higher dosage. Overexpression of coamplified neighboring genes imposes a substantial metabolic burden. While compensatory mutations can decrease inappropriate overexpression of coamplified genes, it takes time for such mutations to arise. The extent to which intrinsic regulatory mechanisms modulate expression of coamplified genes in the immediate aftermath of segmental amplification is largely unknown. To address the collateral effects of segmental amplification, we evolved replicate cultures of an Escherichia coli mutant under conditions that select for higher dosage of an inefficient enzyme whose weak activity limits growth rate. Segmental amplifications encompassing the gene encoding the weak-link enzyme arose in all populations. Amplified regions ranged in size (33 to 125 kb) and copy number (2 to ≥14 copies). We performed RNA-seq and label-free proteomics to quantify expression of amplified genes present at 2, 6, and 14 copies. mRNA expression generally scales with gene copy number, but protein expression scales less well with both gene copy number and mRNA expression. We characterize the molecular mechanisms underlying discrepancies between gene copy number and expression for several cases. We also show that segmental amplifications can have system-wide consequences by indirectly altering expression of nonamplified genes. Our findings indicate that the fitness benefit derived from segmental amplification depends on the combined effects of amplicon size, gene content, and copy number as well as collateral effects on nonamplified genes.

Resource presentation dictates genetic and phenotypic adaptation in yeast

BACKGROUND

Environments shape adaptive trajectories of populations, often leading to adaptive parallelism in identical, and divergence in different environments. However, how does the likelihood of these possibilities change with minute changes in the environment remain unclear.

RESULTS

In this study, we evolved Saccharomyces cerevisiae in environments which differed only in the manner in which the sugar source is presented to the population. In one set of populations, carbon was presented as a mixture of glucose-galactose, and in the other, as melibiose, a glucose-galactose disaccharide. Since the two environments differed in how the two monosaccharides are packaged, we call these environments 'synonymous'. Our results show that even subtle environmental differences can lead to differing phenotypic responses between the two sets of evolved populations. However, despite different adaptive responses, pleiotropic effects of adaptation are largely predictable. We also show that distinct genomic targets of adaptation between the two sets of evolved populations are functionally convergent.

CONCLUSION

This study highlights how subtle environmental differences dictate phenotypic and genetic adaptation of populations. Additionally, these results also suggest the predictive potential of ancestor's fitness in understanding pleiotropic responses. Our work underscores the importance of studying more such environments to understand the generality of adaptive responses in populations.

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.

Evolution in microbial microcosms is highly parallel, regardless of the presence of interacting species

Evolution often follows similar trajectories in replicate populations, suggesting that it may be predictable. However, populations are naturally embedded in multispecies communities, and the extent to which evolution is contingent on the specific species interacting with the focal population is still largely unexplored. Here, we study adaptations in strains of 11 different species, experimentally evolved both in isolation and in various pairwise co-cultures. Although partner-specific effects are detectable, evolution was mostly shared between strains evolved with different partners; similar changes occurred in strains' growth abilities, in community properties, and in about half of the repeatedly mutated genes. This pattern persisted even in species pre-adapted to the abiotic conditions. These findings indicate that evolution may not always depend strongly on the biotic environment, making predictions regarding coevolutionary dynamics less challenging than previously thought. A record of this paper's transparent peer review process is included in the supplemental information.

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.

Highly parallelized laboratory evolution of wine yeasts for enhanced metabolic phenotypes

Adaptive Laboratory Evolution (ALE) of microorganisms can improve the efficiency of sustainable industrial processes important to the global economy. However, stochasticity and genetic background effects often lead to suboptimal outcomes during laboratory evolution. Here we report an ALE platform to circumvent these shortcomings through parallelized clonal evolution at an unprecedented scale. Using this platform, we evolved 104 yeast populations in parallel from many strains for eight desired wine fermentation-related traits. Expansions of both ALE replicates and lineage numbers broadened the evolutionary search spectrum leading to improved wine yeasts unencumbered by unwanted side effects. At the genomic level, evolutionary gains in metabolic characteristics often coincided with distinct chromosome amplifications and the emergence of side-effect syndromes that were characteristic of each selection niche. Several high-performing ALE strains exhibited desired wine fermentation kinetics when tested in larger liquid cultures, supporting their suitability for application. More broadly, our high-throughput ALE platform opens opportunities for rapid optimization of microbes which otherwise could take many years to accomplish.

Genetic conflicts in budding yeast: The 2μ plasmid as a model selfish element

Antagonistic coevolution, arising from genetic conflict, can drive rapid evolution and biological innovation. Conflict can arise both between organisms and within genomes. This review focuses on budding yeasts as a model system for exploring intra- and inter-genomic genetic conflict, highlighting in particular the 2-micron (2μ) plasmid as a model selfish element. The 2μ is found widely in laboratory strains and industrial isolates of Saccharomyces cerevisiae and has long been known to cause host fitness defects. Nevertheless, the plasmid is frequently ignored in the context of genetic, fitness, and evolution studies. Here, I make a case for further exploring the evolutionary impact of the 2μ plasmid as well as other selfish elements of budding yeasts, discuss recent advances, and, finally, future directions for the field.

Chromosome-Level Assembly Reveals a Fifteen-Chromosome Aneuploid Genome and Environmental Adaptation Strategy of Chinese Traditional Medical Fungus Wolfiporia hoelen

The sclerotia of Wolfiporia hoelen are one of the most important traditional Chinese medicines and foods commonly used in China, Japan, Korea, and other Asian countries. To provide a high-quality reference genome and deepen our understanding of the genome of W. hoelen to elucidate various biological phenomena. In this study, we assembled three genomes of W. hoelen using a combination of Nanopore and Illumina sequencing strategies. The fifteen-chromosome genome L7 of W. hoelen was assembled with two-sided telomere and rDNA sequences for the first time. The chromosome count was subsequently confirmed through collinearity analysis, correcting the previous belief that W. hoelen had only fourteen chromosomes. Moreover, the aneuploid genome was discovered in W. hoelen for the first time through sequencing depth analysis of different chromosomes, and only some strains of W. hoelen exhibit aneuploid genomes. According to the genome analysis of homokaryotic offspring and protoplast-isolated strains, a potential variation in chromosome allocation patterns was revealed. Moreover, the gene function enrichment analysis of genes on reduplicated chromosomes demonstrated that aneuploidy in the genome may be the result of environmental adaptation for W. hoelen. The discovery of an aneuploid genome also provides new ideas for genetic improvement of W. hoelen.

Halophytes and heavy metals: A multi-omics approach to understand the role of gene and genome duplication in the abiotic stress tolerance of Cakile maritima

PREMISE

The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole-genome triplication with closely related salt-sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy.

METHODS

Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium.

RESULTS

Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance.

CONCLUSIONS

These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response.

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.

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.

Effect of yeast chromosome II aneuploidy on malate production in sake brewing

Chromosome aneuploidy is a common phenomenon in industrial yeast. Aneuploidy is considered one of the strategies to enhance the industrial properties of Saccharomyces cerevisiae strains. However, the effects of chromosomal aneuploidy on the brewing properties of sake have not been extensively studied. In this study, sake brewing was performed using a series of genome-wide segmental duplicated laboratory S. cerevisiae strains, and the effects of each segmentally duplicated region on sake brewing were investigated. We found that the duplication of specific chromosomal regions affected the production of organic acids and aromatic compounds in sake brewing. As organic acids significantly influence the taste of sake, we focused on the segmental duplication of chromosome II that alters malate levels. Sake yeast Kyokai No. 901 strains with segmental chromosome II duplication were constructed using a polymerase chain reaction-mediated chromosomal duplication method, and sake was brewed using the resultant aneuploid sake yeast strains. The results showed the possibility of developing sake yeast strains exhibiting low malate production without affecting ethanol production capacity. Our study revealed that aneuploidy in yeast alters the brewing properties; in particular, the aneuploidy of chromosome II alters malate production in sake brewing. In conclusion, aneuploidization can be a novel and useful tool to breed sake yeast strains with improved traits, possessing industrial significance.

Template switching between the leading and lagging strands at replication forks generates inverted copy number variants through hairpin-capped extrachromosomal DNA

Inherited and germ-line de novo copy number variants (CNVs) are increasingly found to be correlated with human developmental and cancerous phenotypes. Several models for template switching during replication have been proposed to explain the generation of these gross chromosomal rearrangements. We proposed a model of template switching (ODIRA-origin dependent inverted repeat amplification) in which simultaneous ligation of the leading and lagging strands at diverging replication forks could generate segmental inverted triplications through an extrachromosomal inverted circular intermediate. Here, we created a genetic assay using split-ura3 cassettes to trap the proposed inverted intermediate. However, instead of recovering circular inverted intermediates, we found inverted linear chromosomal fragments ending in native telomeres-suggesting that a template switch had occurred at the centromere-proximal fork of a replication bubble. As telomeric inverted hairpin fragments can also be created through double strand breaks we tested whether replication errors or repair of double stranded DNA breaks were the most likely initiating event. The results from CRISPR/Cas9 cleavage experiments and growth in the replication inhibitor hydroxyurea indicate that it is a replication error, not a double stranded break that creates the inverted junctions. Since inverted amplicons of the SUL1 gene occur during long-term growth in sulfate-limited chemostats, we sequenced evolved populations to look for evidence of linear intermediates formed by an error in replication. All of the data are compatible with a two-step version of the ODIRA model in which sequential template switching at short inverted repeats between the leading and lagging strands at a replication fork, followed by integration via homologous recombination, generates inverted interstitial triplications.

Recurrent copy number variations in the human fungal pathogen Candida parapsilosis

Candida parapsilosis is an opportunistic fungal pathogen with increasing incidence in hospital settings worldwide; however, we lack a comprehensive understanding of the mechanisms promoting its virulence and drug resistance. Bergin et al. systematically quantify the frequency and effect of copy number variation (CNV) across 170 diverse clinical and environmental isolates of C. parapsilosis (Bergin SA, Zhao F, Ryan AP, Müller CA, Nieduszynski CA, Zhai B, Rolling T, Hohl TM, Morio F, Scully J, Wolfe KH, Butler G, 2022, mBio, https://doi.org/10.1128/mbio.01777-22). Using a combination of both short- and long-read whole genome sequencing techniques, they determine the structure and copy number of two CNVs that arose recurrently throughout the evolution of these isolates. Each CNV predominantly amplifies one coding sequence (ARR3 or RTA3); however, the amplitude and recombination breakpoints are variable across the isolates. Amplification of RTA3 correlates with drug resistance and deletion causes drug susceptibility. This study highlights the need for further research into the mechanisms and dynamics of CNV formation and the impact of these CNVs on virulence and drug resistance across diverse fungal pathogens.

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.

The Dynamic Fungal Genome: Polyploidy, Aneuploidy and Copy Number Variation in Response to Stress

Fungal species have dynamic genomes and often exhibit genomic plasticity in response to stress. This genome plasticity often comes with phenotypic consequences that affect fitness and resistance to stress. Fungal pathogens exhibit genome plasticity in both clinical and agricultural settings and often during adaptation to antifungal drugs, posing significant challenges to human health. Therefore, it is important to understand the rates, mechanisms, and impact of large genomic changes. This review addresses the prevalence of polyploidy, aneuploidy, and copy number variation across diverse fungal species, with special attention to prominent fungal pathogens and model species. We also explore the relationship between environmental stress and rates of genomic changes and highlight the mechanisms underlying genotypic and phenotypic changes. A comprehensive understanding of these dynamic fungal genomes is needed to identify novel solutions for the increase in antifungal drug resistance.

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.

Condition-dependent fitness effects of large synthetic chromosome amplifications

Whole-chromosome aneuploidy and large segmental amplifications can have devastating effects in multicellular organisms, from developmental disorders and miscarriage to cancer. Aneuploidy in single-celled organisms such as yeast also results in proliferative defects and reduced viability. Yet, paradoxically, CNVs are routinely observed in laboratory evolution experiments with microbes grown in stressful conditions. The defects associated with aneuploidy are often attributed to the imbalance of many differentially expressed genes on the affected chromosomes, with many genes each contributing incremental effects. An alternate hypothesis is that a small number of individual genes are large effect 'drivers' of these fitness changes when present in an altered copy number. To test these two views, we have employed a collection of strains bearing large chromosomal amplifications that we previously assayed in nutrient-limited chemostat competitions. In this study, we focus on conditions known to be poorly tolerated by aneuploid yeast-high temperature, treatment with the Hsp90 inhibitor radicicol, and growth in extended stationary phase. To identify potential genes with a large impact on fitness, we fit a piecewise constant model to fitness data across chromosome arms, filtering breakpoints in this model by magnitude to focus on regions with a large impact on fitness in each condition. While fitness generally decreased as the length of the amplification increased, we were able to identify 91 candidate regions that disproportionately impacted fitness when amplified. Consistent with our previous work with this strain collection, nearly all candidate regions were condition specific, with only five regions impacting fitness in multiple conditions.

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.

Hotspot of de novo telomere addition stabilizes linear amplicons in yeast grown in sulfate-limiting conditions

Evolution is driven by the accumulation of competing mutations that influence survival. A broad form of genetic variation is the amplification or deletion of DNA (≥50 bp) referred to as copy number variation (CNV). In humans, CNV may be inconsequential, contribute to minor phenotypic differences, or cause conditions such as birth defects, neurodevelopmental disorders, and cancers. To identify mechanisms that drive CNV, we monitored the experimental evolution of Saccharomyces cerevisiae populations grown under sulfate-limiting conditions. Cells with increased copy number of the gene SUL1, which encodes a primary sulfate transporter, exhibit a fitness advantage. Previously, we reported interstitial inverted triplications of SUL1 as the dominant rearrangement in a haploid population. Here, in a diploid population, we find instead that small linear fragments containing SUL1 form and are sustained over several generations. Many of the linear fragments are stabilized by de novo telomere addition within a telomere-like sequence near SUL1 (within the SNF5 gene). Using an assay that monitors telomerase action following an induced chromosome break, we show that this region acts as a hotspot of de novo telomere addition and that required sequences map to a region of <250 base pairs. Consistent with previous work showing that association of the telomere-binding protein Cdc13 with internal sequences stimulates telomerase recruitment, mutation of a four-nucleotide motif predicted to associate with Cdc13 abolishes de novo telomere addition. Our study suggests that internal telomere-like sequences that stimulate de novo telomere addition can contribute to adaptation by promoting genomic plasticity.

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.

A Second Role for the Second Messenger Cyclic-di-GMP in E. coli: Arresting Cell Growth by Altering Metabolic Flow

c-di-GMP primarily controls motile to sessile transitions in bacteria. Diguanylate cyclases (DGCs) catalyze the synthesis of c-di-GMP from two GTP molecules. Typically, bacteria encode multiple DGCs that are activated by specific environmental signals. Their catalytic activity is modulated by c-di-GMP binding to autoinhibitory sites (I-sites). YfiN is a conserved inner membrane DGC that lacks these sites. Instead, YfiN activity is directly repressed by periplasmic YfiR, which is inactivated by redox stress. In Escherichia coli, an additional envelope stress causes YfiN to relocate to the mid-cell to inhibit cell division by interacting with the division machinery. Here, we report a third activity for YfiN in E. coli, where cell growth is inhibited without YfiN relocating to the division site. This action of YfiN is only observed when the bacteria are cultured on gluconeogenic carbon sources, and is dependent on absence of the autoinhibitory sites. Restoration of I-site function relieves the growth-arrest phenotype, and disabling this function in a heterologous DGC causes acquisition of this phenotype. Arrested cells are tolerant to a wide range of antibiotics. We show that the likely cause of growth arrest is depletion of cellular GTP from run-away synthesis of c-di-GMP, explaining the dependence of growth arrest on gluconeogenic carbon sources that exhaust more GTP during production of glucose. This is the first report of c-di-GMP-mediated growth arrest by altering metabolic flow. IMPORTANCE The c-di-GMP signaling network in bacteria not only controls a variety of cellular processes such as motility, biofilms, cell development, and virulence, but does so by a dizzying array of mechanisms. The DGC YfiN singularly represents the versatility of this network in that it not only inhibits motility and promotes biofilms, but also arrests growth in Escherichia coli by relocating to the mid-cell and blocking cell division. The work described here reveals that YfiN arrests growth by yet another mechanism in E. coli, changing metabolic flow. This function of YfiN, or of DGCs without autoinhibitory I-sites, may contribute to antibiotic tolerant persisters in relevant niches such as the gut where gluconeogenic sugars are found.

Resequencing of a Pekin duck breeding population provides insights into the genomic response to short-term artificial selection

BACKGROUND

Short-term, intense artificial selection drives fast phenotypic changes in domestic animals and leaves imprints on their genomes. However, the genetic basis of this selection response is poorly understood. To better address this, we employed the Pekin duck Z2 pure line, in which the breast muscle weight was increased nearly 3-fold after 10 generations of breeding. We denovo assembled a high-quality reference genome of a female Pekin duck of this line (GCA_003850225.1) and identified 8.60 million genetic variants in 119 individuals among 10 generations of the breeding population.

RESULTS

We identified 53 selected regions between the first and tenth generations, and 93.8% of the identified variations were enriched in regulatory and noncoding regions. Integrating the selection signatures and genome-wide association approach, we found that 2 regions covering 0.36 Mb containing UTP25 and FBRSL1 were most likely to contribute to breast muscle weight improvement. The major allele frequencies of these 2 loci increased gradually with each generation following the same trend. Additionally, we found that a copy number variation region containing the entire EXOC4 gene could explain 1.9% of the variance in breast muscle weight, indicating that the nervous system may play a role in economic trait improvement.

CONCLUSIONS

Our study not only provides insights into genomic dynamics under intense artificial selection but also provides resources for genomics-enabled improvements in duck breeding.

Copy Number Variation and Allele Ratio Analysis in Candida albicans Using Whole Genome Sequencing Data

Whole genome sequencing of human fungal pathogens has revolutionized the speed and accuracy in which sequence variants that cause antifungal resistance can be identified. Genome rearrangements resulting in copy number variation (CNV) are a significant source of acquired antifungal drug resistance across diverse fungal species. Some CNVs are transient in nature, while other CNVs are stable and well tolerated even in the absence of antifungal drugs. By visualizing whole genome sequencing read depth as a function of genomic location, CNVs and CNV breakpoints (genomic positions where the copy number changes occur relative to the rest of the genome) are rapidly identified. A similar analysis can be used to visualize allele ratio changes that occur across the genomes of heterozygous fungal species, both in the presence and absence of CNVs. This protocol walks through the bioinformatic analysis of CNVs and allele ratios utilizing free, open-source visualization tools. We provide code to use with an example dataset (matched antifungal drug-sensitive and drug-resistant Candida albicans isolates) and notes on how to expand this protocol to other fungal genomes.

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.

Systematic Analysis of Copy Number Variations in the Pathogenic Yeast Candida parapsilosis Identifies a Gene Amplification in RTA3 That is Associated with Drug Resistance

We analyzed the genomes of 170 C. parapsilosis isolates and identified multiple copy number variations (CNVs). We identified two genes, RTA3 (CPAR2_104610) and ARR3 (CPAR2_601050), each of which was the target of multiple independent amplification events. Phylogenetic analysis shows that most of these amplifications originated only once. For ARR3, which encodes a putative arsenate transporter, 8 distinct CNVs were identified, ranging in size from 2.3 kb to 10.5 kb with 3 to 23 copies. For RTA3, 16 distinct CNVs were identified, ranging in size from 0.3 kb to 4.5 kb with 2 to ~50 copies. One unusual amplification resulted in a DUP-TRP/INV-DUP structure similar to some human CNVs. RTA3 encodes a putative phosphatidylcholine (PC) floppase which is known to regulate the inward translocation of PC in Candida albicans. We found that an increased copy number of RTA3 correlated with resistance to miltefosine, an alkylphosphocholine drug that affects PC metabolism. Additionally, we conducted an adaptive laboratory evolution experiment in which two C. parapsilosis isolates were cultured in increasing concentrations of miltefosine. Two genes, CPAR2_303950 and CPAR2_102700, coding for putative PC flippases homologous to S. cerevisiae DNF1 gained homozygous protein-disrupting mutations in the evolved strains. Overall, our results show that C. parapsilosis can gain resistance to miltefosine, a drug that has recently been granted orphan drug designation approval by the United States Food and Drug Administration for the treatment of invasive candidiasis, through both CNVs or loss-of-function alleles in one of the flippase genes. IMPORTANCE Copy number variations (CNVs) are an important source of genomic diversity that have been associated with drug resistance. We identify two unusual CNVs in the human fungal pathogen Candida parapsilosis. Both target a single gene (RTA3 or ARR3), and they have occurred multiple times in multiple isolates. The copy number of RTA3, a putative floppase that controls the inward translocation of lipids in the cell membrane, correlates with resistance to miltefosine, a derivative of phosphatidylcholine (PC) that was originally developed as an anticancer drug. In 2021, miltefosine was designated an orphan drug by the United States Food and Drug Administration for the treatment of invasive candidiasis. Importantly, we find that resistance to miltefosine is also caused by mutations in flippases, which control the outward movement of lipids, and that many C. parapsilosis isolates are prone to easily acquiring an increased resistance to miltefosine.

The landscape of transcriptional and translational changes over 22 years of bacterial adaptation

Organisms can adapt to an environment by taking multiple mutational paths. This redundancy at the genetic level, where many mutations have similar phenotypic and fitness effects, can make untangling the molecular mechanisms of complex adaptations difficult. Here, we use the Escherichia coli long-term evolution experiment (LTEE) as a model to address this challenge. To understand how different genomic changes could lead to parallel fitness gains, we characterize the landscape of transcriptional and translational changes across 12 replicate populations evolving in parallel for 50,000 generations. By quantifying absolute changes in mRNA abundances, we show that not only do all evolved lines have more mRNAs but that this increase in mRNA abundance scales with cell size. We also find that despite few shared mutations at the genetic level, clones from replicate populations in the LTEE are remarkably similar in their gene expression patterns at both the transcriptional and translational levels. Furthermore, we show that the majority of the expression changes are due to changes at the transcriptional level with very few translational changes. Finally, we show how mutations in transcriptional regulators lead to consistent and parallel changes in the expression levels of downstream genes. These results deepen our understanding of the molecular mechanisms underlying complex adaptations and provide insights into the repeatability of evolution.

Four layer multi-omics reveals molecular responses to aneuploidy in Leishmania

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond.

Mesostats—A multiplexed, low-cost, do-it-yourself continuous culturing system for experimental evolution of mesocosms

Microbial experimental evolution allows studying evolutionary dynamics in action and testing theory predictions in the lab. Experimental evolution in chemostats (i.e. continuous flow through cultures) has recently gained increased interest as it allows tighter control of selective pressures compared to static batch cultures, with a growing number of efforts to develop systems that are easier and cheaper to construct. This protocol describes the design and construction of a multiplexed chemostat array (dubbed “mesostats”) designed for cultivation of algae in 16 concurrent populations, specifically intended for studying adaptation to herbicides. We also present control data from several experiments run on the system to show replicability, data illustrating the effects of common issues like leaks, contamination and clumps, and outline possible modifications and adaptations of the system for future research.

Reliance of Host-Encoded Regulators of Retromobility on Ty1 Promoter Activity or Architecture

The Ty1 retrotransposon family is maintained in a functional but dormant state by its host, Saccharomyces cerevisiae. Several hundred RHF and RTT genes encoding co-factors and restrictors of Ty1 retromobility, respectively, have been identified. Well-characterized examples include MED3 and MED15, encoding subunits of the Mediator transcriptional co-activator complex; control of retromobility by Med3 and Med15 requires the Ty1 promoter in the U3 region of the long terminal repeat. To characterize the U3-dependence of other Ty1 regulators, we screened a library of 188 known rhf and rtt mutants for altered retromobility of Ty1his3AI expressed from the strong, TATA-less TEF1 promoter or the weak, TATA-containing U3 promoter. Two classes of genes, each including both RHFs and RTTs, were identified. The first class comprising 82 genes that regulated Ty1his3AI retromobility independently of U3 is enriched for RHF genes that restrict the G1 phase of the cell cycle and those involved in transcriptional elongation and mRNA catabolism. The second class of 51 genes regulated retromobility of Ty1his3AI driven only from the U3 promoter. Nineteen U3-dependent regulators (U3DRs) also controlled retromobility of Ty1his3AI driven by the weak, TATA-less PSP2 promoter, suggesting reliance on the low activity of U3. Thirty-one U3DRs failed to modulate P_PSP2-Ty1his3AI retromobility, suggesting dependence on the architecture of U3. To further investigate the U3-dependency of Ty1 regulators, we developed a novel fluorescence-based assay to monitor expression of p22-Gag, a restriction factor expressed from the internal Ty1i promoter. Many U3DRs had minimal effects on levels of Ty1 RNA, Ty1i RNA or p22-Gag. These findings uncover a role for the Ty1 promoter in integrating signals from diverse host factors to modulate Ty1 RNA biogenesis or fate.

Microbial experimental evolution in a massively multiplexed and high-throughput era

Experimental evolution with microbial model systems has transformed our understanding of the basic rules underlying ecology and evolution. Experiments leveraging evolution as a central feature put evolutionary theories to the test, and modern sequencing and engineering tools then characterized the molecular basis of adaptation. As theory and experimentations refined our understanding of evolution, a need to increase throughput and experimental complexity has emerged. Here, we summarize recent technologies that have made high-throughput experiments practical and highlight studies that have capitalized on these tools, defining an exciting new era in microbial experimental evolution. Multiple research directions previously limited by experimental scale are now accessible for study and we believe applying evolutionary lessons from in vitro studies onto these applied settings has the potential for major innovations and discoveries across ecology and medicine.

Genome plasticity in Candida albicans: A cutting-edge strategy for evolution, adaptation, and survival

Candida albicans is the most implicated fungal species that grows as a commensal or opportunistic pathogen in the human host. It is associated with many life-threatening infections, especially in immunocompromised persons. The genome of Candida albicans is very flexible and can withstand a wide assortment of variations in a continuously changing environment. Thus, genome plasticity is central to its adaptation and has long been of considerable interest. C. albicans has a diploid heterozygous genome that is highly dynamic and can display variation from small to large scale chromosomal rearrangement and aneuploidy, which have implications in drug resistance, virulence, and pathogenicity. This review presents an up-to-date overview of recent genomic studies involving C. albicans. It discusses the accumulating evidence that shows how mitotic recombination events, ploidy dynamics, aneuploidy, and loss of heterozygosity (LOH) influence evolution, adaptation, and survival in C. albicans. Understanding the factors that affect the genome is crucial for a proper understanding of species and rapid development and adjustment of therapeutic strategies to mitigate their spread.

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.

Whole-cell modeling in yeast predicts compartment-specific proteome constraints that drive metabolic strategies

When conditions change, unicellular organisms rewire their metabolism to sustain cell maintenance and cellular growth. Such rewiring may be understood as resource re-allocation under cellular constraints. Eukaryal cells contain metabolically active organelles such as mitochondria, competing for cytosolic space and resources, and the nature of the relevant cellular constraints remain to be determined for such cells. Here, we present a comprehensive metabolic model of the yeast cell, based on its full metabolic reaction network extended with protein synthesis and degradation reactions. The model predicts metabolic fluxes and corresponding protein expression by constraining compartment-specific protein pools and maximising growth rate. Comparing model predictions with quantitative experimental data suggests that under glucose limitation, a mitochondrial constraint limits growth at the onset of ethanol formation—known as the Crabtree effect. Under sugar excess, however, a constraint on total cytosolic volume dictates overflow metabolism. Our comprehensive model thus identifies condition-dependent and compartment-specific constraints that can explain metabolic strategies and protein expression profiles from growth rate optimisation, providing a framework to understand metabolic adaptation in eukaryal cells.

Metabolically active organelles compete for cytosolic space and resources during metabolism rewiring. Here, the authors develop a computational model of yeast metabolism and resource allocation to predict condition- and compartment-specific proteome constraints that govern metabolic strategies.

Reactivation of transposable elements following hybridization in fission yeast

Hybridization is thought to reactivate transposable elements (TEs) that were efficiently suppressed in the genomes of the parental hosts. Here, we provide evidence for this “genomic shock hypothesis” in the fission yeast Schizosaccharomyces pombe. In this species, two divergent lineages (Sp and Sk) have experienced recent, likely human-induced, hybridization. We used long-read sequencing data to assemble genomes of 37 samples derived from 31 S. pombe strains spanning a wide range of ancestral admixture proportions. A comprehensive TE inventory revealed exclusive presence of long terminal repeat (LTR) retrotransposons. Sequence analysis of active full-length elements, as well as solo LTRs, revealed a complex history of homologous recombination. Population genetic analyses of syntenic sequences placed insertion of many solo LTRs before the split of the Sp and Sk lineages. Most full-length elements were inserted more recently, after hybridization. With the exception of a single full-length element with signs of positive selection, both solo LTRs and, in particular, full-length elements carry signatures of purifying selection indicating effective removal by the host. Consistent with reactivation upon hybridization, the number of full-length LTR retrotransposons, varying extensively from zero to 87 among strains, significantly increases with the degree of genomic admixture. This study gives a detailed account of global TE diversity in S. pombe, documents complex recombination histories within TE elements, and provides evidence for the “genomic shock hypothesis.”

DNA Repair in Haploid Context

DNA repair is a well-covered topic as alteration of genetic integrity underlies many pathological conditions and important transgenerational consequences. Surprisingly, the ploidy status is rarely considered although the presence of homologous chromosomes dramatically impacts the repair capacities of cells. This is especially important for the haploid gametes as they must transfer genetic information to the offspring. An understanding of the different mechanisms monitoring genetic integrity in this context is, therefore, essential as differences in repair pathways exist that differentiate the gamete’s role in transgenerational inheritance. Hence, the oocyte must have the most reliable repair capacity while sperm, produced in large numbers and from many differentiation steps, are expected to carry de novo variations. This review describes the main DNA repair pathways with a special emphasis on ploidy. Differences between Saccharomyces cerevisiae and Schizosaccharomyces pombe are especially useful to this aim as they can maintain a diploid and haploid life cycle respectively.

Ploidy and recombination proficiency shape the evolutionary adaptation to constitutive DNA replication stress

In haploid budding yeast, evolutionary adaptation to constitutive DNA replication stress alters three genome maintenance modules: DNA replication, the DNA damage checkpoint, and sister chromatid cohesion. We asked how these trajectories depend on genomic features by comparing the adaptation in three strains: haploids, diploids, and recombination deficient haploids. In all three, adaptation happens within 1000 generations at rates that are correlated with the initial fitness defect of the ancestors. Mutations in individual genes are selected at different frequencies in populations with different genomic features, but the benefits these mutations confer are similar in the three strains, and combinations of these mutations reproduce the fitness gains of evolved populations. Despite the differences in the selected mutations, adaptation targets the same three functional modules in strains with different genomic features, revealing a common evolutionary response to constitutive DNA replication stress.

Molecular signatures of resource competition: Clonal interference favors ecological diversification and can lead to incipient speciation

Microbial ecosystems harbor an astonishing diversity that can persist for long times. To understand how such diversity is structured and maintained, ecological and evolutionary processes need to be integrated at similar timescales. Here, we study a model of resource competition that allows for evolution via de novo mutation, and focus on rapidly adapting asexual populations with large mutational inputs, as typical of many bacteria species. We characterize the adaptation and diversification of an initially maladapted population and show how the eco-evolutionary dynamics are shaped by the interaction between simultaneously emerging lineages - clonal interference. We find that in large populations, more intense clonal interference can foster diversification under sympatry, increasing the probability that phenotypically and genetically distinct clusters coexist. In smaller populations, the accumulation of deleterious and compensatory mutations can push further the diversification process and kick-start speciation. Our findings have implications beyond microbial populations, providing novel insights about the interplay between ecology and evolution in clonal populations.

Temperature affects the repeatability of evolution in the microbial eukaryote Tetrahymena thermophila

Evolutionary biologists have long sought to understand what factors affect the repeatability of adaptive outcomes. To better understand the role of temperature in determining the repeatability of adaptive trajectories, we evolved populations of different genotypes of the ciliate Tetrahymena thermophila at low and high temperatures and followed changes in growth rate over 6,500 generations. As expected, growth rate increased with a decelerating rate for all populations; however, there were differences in the patterns of evolution at the two temperatures. The growth rates of the different genotypes tended to converge as evolution proceeded at both temperatures, but this convergence was quicker and more pronounced at the higher temperature. Additionally, over the first 4,000 generations we found greater repeatability of evolution, in terms of change in growth rate, among replicates of the same genotype at the higher temperature. Finally, we found limited evidence of trade-offs in fitness between temperatures, and an asymmetry in the correlated responses, whereby evolution in a high temperature increases growth rate at the lower temperature significantly more than the reverse. These results demonstrate the importance of temperature in determining the repeatability of evolutionary trajectories for the eukaryotic microbe Tetrahymena thermophila and may provide clues to how temperature affects evolution more generally.

Adaptive response to wine selective pressures shapes the genome of a Saccharomyces interspecies hybrid

During industrial processes, yeasts are exposed to harsh conditions, which eventually lead to adaptation of the strains. In the laboratory, it is possible to use experimental evolution to link the evolutionary biology response to these adaptation pressures for the industrial improvement of a specific yeast strain. In this work, we aimed to study the adaptation of a wine industrial yeast in stress conditions of the high ethanol concentrations present in stopped fermentations and secondary fermentations in the processes of champagne production. We used a commercial Saccharomyces cerevisiae × S. uvarum hybrid and assessed its adaptation in a modified synthetic must (M-SM) containing high ethanol, which also contained metabisulfite, a preservative that is used during wine fermentation as it converts to sulfite. After the adaptation process under these selected stressful environmental conditions, the tolerance of the adapted strain (H14A7-etoh) to sulfite and ethanol was investigated, revealing that the adapted hybrid is more resistant to sulfite compared to the original H14A7 strain, whereas ethanol tolerance improvement was slight. However, a trade-off in the adapted hybrid was found, as it had a lower capacity to ferment glucose and fructose in comparison with H14A7. Hybrid genomes are almost always unstable, and different signals of adaptation on H14A7-etoh genome were detected. Each subgenome present in the adapted strain had adapted differently. Chromosome aneuploidies were present in S. cerevisiae chromosome III and in S. uvarum chromosome VII–XVI, which had been duplicated. Moreover, S. uvarum chromosome I was not present in H14A7-etoh and a loss of heterozygosity (LOH) event arose on S. cerevisiae chromosome I. RNA-sequencing analysis showed differential gene expression between H14A7-etoh and H14A7, which can be easily correlated with the signals of adaptation that were found on the H14A7-etoh genome. Finally, we report alterations in the lipid composition of the membrane, consistent with conserved tolerance mechanisms.

Natural variation in the consequences of gene overexpression and its implications for evolutionary trajectories

Copy number variation through gene or chromosome amplification provides a route for rapid phenotypic variation and supports the long-term evolution of gene functions. Although the evolutionary importance of copy-number variation is known, little is understood about how genetic background influences its tolerance. Here, we measured fitness costs of over 4000 overexpressed genes in 15 Saccharomyces cerevisiae strains representing different lineages, to explore natural variation in tolerating gene overexpression (OE). Strain-specific effects dominated the fitness costs of gene OE. We report global differences in the consequences of gene OE, independent of the amplified gene, as well as gene-specific effects that were dependent on the genetic background. Natural variation in the response to gene OE could be explained by several models, including strain-specific physiological differences, resource limitations, and regulatory sensitivities. This work provides new insight on how genetic background influences tolerance to gene amplification and the evolutionary trajectories accessible to different backgrounds.

Multiplexing mutation rate assessment: determining pathogenicity of Msh2 variants in Saccharomyces cerevisiae

Despite the fundamental importance of mutation rate as a driving force in evolution and disease risk, common methods to assay mutation rate are time-consuming and tedious. Established methods such as fluctuation tests and mutation accumulation experiments are low-throughput and often require significant optimization to ensure accuracy. We established a new method to determine the mutation rate of many strains simultaneously by tracking mutation events in a chemostat continuous culture device and applying deep sequencing to link mutations to alleles of a DNA-repair gene. We applied this method to assay the mutation rate of hundreds of Saccharomyces cerevisiae strains carrying mutations in the gene encoding Msh2, a DNA repair enzyme in the mismatch repair pathway. Loss-of-function mutations in MSH2 are associated with hereditary nonpolyposis colorectal cancer, an inherited disorder that increases risk for many different cancers. However, the vast majority of MSH2 variants found in human populations have insufficient evidence to be classified as either pathogenic or benign. We first benchmarked our method against Luria-Delbrück fluctuation tests using a collection of published MSH2 missense variants. Our pooled screen successfully identified previously characterized nonfunctional alleles as high mutators. We then created an additional 185 human missense variants in the yeast ortholog, including both characterized and uncharacterized alleles curated from ClinVar and other clinical testing data. In a set of alleles of known pathogenicity, our assay recapitulated ClinVar's classification; we then estimated pathogenicity for 157 variants classified as uncertain or conflicting reports of significance. This method is capable of studying the mutation rate of many microbial species and can be applied to problems ranging from the generation of high-fidelity polymerases to measuring the frequency of antibiotic resistance emergence.

Proteomic analysis of serial isolates of Trichosporon asahii identifies host-specific adaptations using the TMT/MRM approach

The opportunistic fungal pathogen Trichosporon asahii (T. asahii) is an important causal agent of mortality in immunocompromised patients and associated with frequent relapses, even with sufficient antifungal treatment. Investigating the proteomes of initial and recurrent isolates may help to identify within-host adaptive changes. In this study, using tandem mass tag (TMT)-labeling combined with liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) technology, we analyzed the proteomes of two T. asahii strains that were isolated 15 years apart from the same patient who suffered initial and recurrent episodes of systemic disseminated trichosporonosis. A total of 597 differentially expressed proteins were identified. Functional analysis showed that the increased proteins were primarily concentrated on peptide/protein/energy/drug metabolism and translation. Most of the results were determined to be consistent with the findings of phenotypic assays, such as tests for drug susceptibility, temperature growth, biofilm formation, melanization and paromomycin assays. Moreover, we performed multiple reaction monitoring (MRM) mass spectrometry to verify 27 candidate proteins, and the results of this experiment were also highly consistent with the results of the TMT analysis. Therefore, to the best of our knowledge, these data provide the first molecular evidence of how the T. asahii proteome changes related to host-specific adaptation during human infection. SIGNIFICANCE: Systemic infection with Trichosporon asahii (T. asahii) has recently been recognized as an important causal agent of mortality in immunocompromised patients. Although triazole treatment usually works efficiently in the early phase of infection, many patients relapse. Hence, comparative analyses of the proteomics of initial and recurrent isolates may reveal evidence of adaptive changes within the host. Our study demonstrates that the recurrent strain has undergone proteomic changes using tandem mass tag (TMT)-labeling combined with liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS). Moreover, the results of phenotypic assays, including drug susceptibility, temperature growth, biofilm formation, melanization and paromomycin assays, were highly consistent with the proteomic changes, and multiple reaction monitoring (MRM) verification also showed similar trends to the TMT results. In summary, our study is the first to investigate the adaptation of T. asahii under pressure from antifungal chemotherapy and host immune responses.

Adaptive laboratory evolution of β-caryophyllene producing Saccharomyces cerevisiae

Background

β-Caryophyllene is a plant terpenoid with therapeutic and biofuel properties. Production of terpenoids through microbial cells is a potentially sustainable alternative for production. Adaptive laboratory evolution is a complementary technique to metabolic engineering for strain improvement, if the product-of-interest is coupled with growth. Here we use a combination of pathway engineering and adaptive laboratory evolution to improve the production of β-caryophyllene, an extracellular product, by leveraging the antioxidant potential of the compound.

Results

Using oxidative stress as selective pressure, we developed an adaptive laboratory evolution that worked to evolve an engineered β-caryophyllene producing yeast strain for improved production within a few generations. This strategy resulted in fourfold increase in production in isolated mutants. Further increasing the flux to β-caryophyllene in the best evolved mutant achieved a titer of 104.7 ± 6.2 mg/L product. Genomic analysis revealed a gain-of-function mutation in the a-factor exporter STE6 was identified to be involved in significantly increased production, likely as a result of increased product export.

Conclusion

An optimized selection strategy based on oxidative stress was developed to improve the production of the extracellular product β-caryophyllene in an engineered yeast strain. Application of the selection strategy in adaptive laboratory evolution resulted in mutants with significantly increased production and identification of novel responsible mutations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01598-z.

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.

Transposable Element Mobilization in Interspecific Yeast Hybrids

Barbara McClintock first hypothesized that interspecific hybridization could provide a “genomic shock” that leads to the mobilization of transposable elements (TEs). This hypothesis is based on the idea that regulation of TE movement is potentially disrupted in hybrids. However, the handful of studies testing this hypothesis have yielded mixed results. Here, we set out to identify if hybridization can increase transposition rate and facilitate colonization of TEs in Saccharomyces cerevisiae × Saccharomyces uvarum interspecific yeast hybrids. Saccharomyces cerevisiae have a small number of active long terminal repeat retrotransposons (Ty elements), whereas their distant relative S. uvarum have lost the Ty elements active in S. cerevisiae. Although the regulation system of Ty elements is known in S. cerevisiae, it is unclear how Ty elements are regulated in other Saccharomyces species, and what mechanisms contributed to the loss of most classes of Ty elements in S. uvarum. Therefore, we first assessed whether TEs could insert in the S. uvarum sub-genome of a S. cerevisiae × S. uvarum hybrid. We induced transposition to occur in these hybrids and developed a sequencing technique to show that Ty elements insert readily and nonrandomly in the S. uvarum genome. We then used an in vivo reporter construct to directly measure transposition rate in hybrids, demonstrating that hybridization itself does not alter rate of mobilization. However, we surprisingly show that species-specific mitochondrial inheritance can change transposition rate by an order of magnitude. Overall, our results provide evidence that hybridization can potentially facilitate the introduction of TEs across species boundaries and alter transposition via mitochondrial transmission, but that this does not lead to unrestrained proliferation of TEs suggested by the genomic shock theory.

Expression attenuation as a mechanism of robustness against gene duplication

Gene duplication is ubiquitous and a major driver of phenotypic diversity across the tree of life, but its immediate consequences are not fully understood. Deleterious effects would decrease the probability of retention of duplicates and prevent their contribution to long-term evolution. One possible detrimental effect of duplication is the perturbation of the stoichiometry of protein complexes. Here, we measured the fitness effects of the duplication of 899 essential genes in the budding yeast using high-resolution competition assays. At least 10% of genes caused a fitness disadvantage when duplicated. Intriguingly, the duplication of most protein complex subunits had small to nondetectable effects on fitness, with few exceptions. We selected four complexes with subunits that had an impact on fitness when duplicated and measured the impact of individual gene duplications on their protein-protein interactions. We found that very few duplications affect both fitness and interactions. Furthermore, large complexes such as the 26S proteasome are protected from gene duplication by attenuation of protein abundance. Regulatory mechanisms that maintain the stoichiometric balance of protein complexes may protect from the immediate effects of gene duplication. Our results show that a better understanding of protein regulation and assembly in complexes is required for the refinement of current models of gene duplication.

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.