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.

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

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

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

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

Small RNAs Are Implicated in Regulation of Gene and Transposable Element Expression in the Protist Trichomonas vaginalis.

Trichomonas vaginalis is the causative agent of trichomoniasis, the most prevalent nonviral sexually transmitted infection worldwide. Repetitive elements, including transposable elements (TEs) and virally derived repeats, comprise more than half of the ∼160-Mb T. vaginalis genome. An intriguing question is how the parasite controls its potentially lethal complement of mobile elements, which can disrupt transcription of protein-coding genes and genome functions. In this study, we generated high-throughput RNA sequencing (RNA-Seq) and small RNA-Seq data sets in triplicate for the T. vaginalis G3 reference strain and characterized the mRNA and small RNA populations and their mapping patterns along all six chromosomes. Mapping the RNA-Seq transcripts to the genome revealed that the majority of genes predicted within repetitive elements are not expressed. Interestingly, we identified a novel species of small RNA that maps bidirectionally along the chromosomes and is correlated with reduced protein-coding gene expression and reduced RNA-Seq coverage in repetitive elements. This novel small RNA family may play a regulatory role in gene and repetitive element expression. Our results identify a possible small RNA pathway mechanism by which the parasite regulates expression of genes and TEs and raise intriguing questions as to the role repeats may play in shaping T. vaginalis genome evolution and the diversity of small RNA pathways in general.IMPORTANCE Trichomoniasis, caused by the protozoan Trichomonas vaginalis, is the most common nonviral sexually transmitted infection in humans. The millions of cases each year have sequelae that may include complications during pregnancy and increased risk of HIV infection. Given its evident success in this niche, it is paradoxical that T. vaginalis harbors in its genome thousands of transposable elements that have the potential to be extremely detrimental to normal genomic function. In many organisms, transposon expression is regulated by the activity of endogenously expressed short (∼21 to 35 nucleotides [nt]) small RNA molecules that effect gene silencing by targeting mRNAs for degradation or by recruiting epigenetic silencing machinery to locations in the genome. Our research has identified small RNA molecules correlated with reduced expression of T. vaginalis genes and transposons. This suggests that a small RNA pathway is a major contributor to gene expression patterns in the parasite and opens up new avenues for investigation into small RNA biogenesis, function, and diversity.

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

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