Grace Avecilla

Copy number variants (CNVs) are DNA duplications and deletions that cause genetic variation, underlying rapid adaptive evolution. CNVs often confer selective advantages but can also incur fitness costs. Evolution of Saccharomyces cerevisiae in nutrient-limited chemostats recurrently selects for amplifications of nutrient transporter genes. However, their fate upon return to a non-selective environment remains unknown. To investigate CNV fitness and stability upon removing the original selection pressure, we studied 15 CNV lineages (11 segmental and 4 whole-chromosomal amplifications) selected in nitrogen-limited chemostats. CNV stability was monitored using fluorescent reporters during propagation in nutrient-rich batch cultures for 110 to 220 generations. All aneuploid lineages showed rapid CNV loss and reversion to a single-copy genotype, whereas segmental amplifications were remarkably stable; one of the 11 strains reverted. Pairwise fitness competitions in rich media revealed strong fitness defects associated solely with CNVs that reverted; reversion led to increased fitness. Using simulation-based inference to estimate reversion rates and fitness effects, we determined negative selection as the primary driver of CNV loss. Whole-genome sequencing revealed that reversion of aneuploids and a segmental amplification left no evidence of prior CNV existence, rendering revertant genomes indistinguishable from the single-copy ancestor. Detailed characterization of a partial revertant identified chromosomal translocation, suggesting that extant CNVs can undergo structural diversification. Our findings provide novel evidence that most segmental CNVs adapted to nitrogen limitation are stable upon removal of selection, but costly gene amplifications are readily reversible. Together, these highlight the importance of CNVs in both long-term genome evolution and rapid, reversible adaptation to transient selection.

Copy number variants (CNVs) are DNA duplications and deletions that cause genetic variation, underlying rapid adaptive evolution. CNVs often confer selective advantages, but can also incur fitness costs. Evolution of Saccharomyces cerevisiae in nutrient-limited chemostats recurrently selects for amplifications of nutrient transporter genes. However, their fate upon return to a non-selective environment remains unknown. To investigate CNV fitness and stability upon removing the original selection pressure, we studied 15 CNV lineages (11 segmental, 4 whole-chromosomal amplifications) selected in nitrogen-limited chemostats. CNV stability was monitored using fluorescent reporters during propagation in nutrient-rich batch cultures for 110–220 generations. All aneuploid lineages showed rapid CNV loss and reversion to a single-copy genotype, whereas segmental amplifications were remarkably stable– one of 11 strains reverted. Pairwise fitness competitions in rich media revealed strong fitness defects associated solely with CNVs that reverted; reversion led to increased fitness. Using simulation-based inference to estimate reversion rates and fitness effects, we determined negative selection as the primary driver of CNV loss. Whole-genome sequencing revealed that reversion of aneuploids and a segmental amplification left no evidence of prior CNV existence, rendering revertant genomes indistinguishable from the single-copy ancestor. Detailed characterization of a partial revertant identified chromosomal translocation, suggesting that extant CNVs can undergo structural diversification. Our findings provide novel evidence that most segmental CNVs adapted to nitrogen limitation are stable upon removal of selection, but costly gene amplifications are readily reversible. Together, these highlight the importance of CNVs in both long-term genome evolution and rapid, reversible adaptation to transient selection.

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 investigatedSaccharomyces 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 andBMH1that 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.

Copy number variants (CNVs), duplications and deletions of genomic content, 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 have CNVs of variable structures but with multiple copies of the gene GAP1 . Although beneficial in glutamine-limited chemostats, CNVs result in decreased fitness compared with the euploid ancestor in rich media. We used transposon mutagenesis to investigate mutational tolerance and genetic interactions with CNVs. We find that CNVs confer novel mutational tolerance in amplified essential genes and novel genetic interactions. We validated a novel genetic interaction with BMH1 . CNV strains have increased mutational tolerance in genes related to translation, and reduced mutational tolerance in genes related to mitochondrial function. We performed RNAseq and found that transcriptional dosage compensation does not affect the majority of genes amplified by CNVs. Furthermore, we do not find that CNV strains exhibit previously described transcriptional signatures of aneuploidy. Instead, CNV strains exhibit downregulation of genes involved in cellular respiration, nucleoside biosynthetic processes, and small molecule metabolism, and upregulation of genes involved in transposition, nucleic acid metabolic processes, and siderophore transport. 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.

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⁻⁴ 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.

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 two 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 the observed CNV adaptation dynamics to estimate the rate at which beneficial CNVs are introduced through de novo mutation and their fitness effects using simulation-based Bayesian likelihood-free inference approaches. We tested the suitability of two evolutionary models: a standard Wright-Fisher model and a chemostat growth model. We evaluated two 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 yeast as 10 ^−4.7 -10 ⁻⁴ per cell division, and a selection coefficient of 0.04 - 0.1 per generation for GAP1 CNVs in glutamine-limited chemostats. We experimentally validated our estimates using barcode lineage tracking and pairwise fitness assays. Our results are consistent with a high beneficial CNV supply rate that is 10-fold greater than the estimated rates of beneficial single-nucleotide mutations, explaining their outsized importance in rapid adaptive evolution. More generally, our study demonstrates the utility of novel simulation-based likelihood-free inference methods for inferring the rates and effects of evolutionary processes from empirical data.

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 10²–10⁴ 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. Duplications and deletions of genomic sequence, known as copy number variants, are a common source of genetic diversity across all domains of life. Copy number variants play a crucial role in driving evolutionary processes but can also cause genetic disease and cancer. Although copy number variants are important drivers of diversity, adaptation, and disease, the underlying dynamics of their formation and selection are poorly understood. Copy number variants are difficult to detect, especially when present at low frequencies in heterogenous evolving populations. To overcome this challenge, we developed a novel fluorescent reporter that allows us to visualize copy number variants as they emerge and to track them throughout hundreds of generations of laboratory evolution. We show that copy number variants arise early and repeatedly, that they are diverse in size and copy number, and that they are generated at a high rate, leading to competition among cells containing different copy number variants. Molecular characterization of copy number variants indicates that many of them are likely generated by errors during DNA replication. This method is broadly applicable to studying the molecular mechanisms underlying formation of copy number variants, as well as their role in driving evolutionary processes and cancer.