Origin-Dependent Inverted-Repeat Amplification: Tests of a Model for Inverted DNA Amplification

Transcription as source of genetic heterogeneity in budding yeast

Transcription presents challenges to genome stability both directly, by altering genome topology and exposing single-stranded DNA to chemical insults and nucleases, and indirectly by introducing obstacles to the DNA replication machinery. Such obstacles include the RNA polymerase holoenzyme itself, DNA-bound regulatory factors, G-quadruplexes and RNA-DNA hybrid structures known as R-loops. Here, we review the detrimental impacts of transcription on genome stability in budding yeast, as well as the mitigating effects of transcription-coupled nucleotide excision repair and of systems that maintain DNA replication fork processivity and integrity. Interactions between DNA replication and transcription have particular potential to induce mutation and structural variation, but we conclude that such interactions must have only minor effects on DNA replication by the replisome with little if any direct mutagenic outcome. However, transcription can significantly impair the fidelity of replication fork rescue mechanisms, particularly Break Induced Replication, which is used to restart collapsed replication forks when other means fail. This leads to de novo mutations, structural variation and extrachromosomal circular DNA formation that contribute to genetic heterogeneity, but only under particular conditions and in particular genetic contexts, ensuring that the bulk of the genome remains extremely stable despite the seemingly frequent interactions between transcription and DNA replication.

A unifying model that explains the origins of human inverted copy number variants

With the release of the telomere-to-telomere human genome sequence and the availability of both long-read sequencing and optical genome mapping techniques, the identification of copy number variants (CNVs) and other structural variants is providing new insights into human genetic disease. Different mechanisms have been proposed to account for the novel junctions in these complex architectures, including aberrant forms of DNA replication, non-allelic homologous recombination, and various pathways that repair DNA breaks. Here, we have focused on a set of structural variants that include an inverted segment and propose that they share a common initiating event: an inverted triplication with long, unstable palindromic junctions. The secondary rearrangement of these palindromes gives rise to the various forms of inverted structural variants. We postulate that this same mechanism (ODIRA: origin-dependent inverted-repeat amplification) that creates the inverted CNVs in inherited syndromes also generates the palindromes found in cancers.

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.

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.

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.

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.

A unifying model for extrachromosomal circular DNA load in eukaryotic cells

Extrachromosomal circular DNA (eccDNA) with exons and whole genes are common features of eukaryotic cells. Work from especially tumours and the yeast Saccharomyces cerevisiae has revealed that eccDNA can provide large selective advantages and disadvantages. Besides the phenotypic effect due to expression of an eccDNA fragment, eccDNA is different from other mutations in that it is released from 1:1 segregation during cell division. This means that eccDNA can quickly change copy number, pickup secondary mutations and reintegrate into a chromosome to establish substantial genetic variation that could not have evolved via canonical mechanisms. We propose a unifying 5-factor model for conceptualizing the eccDNA load of a eukaryotic cell, emphasizing formation, replication, segregation, selection and elimination. We suggest that the magnitude of these sequential events and their interactions determine the copy number of eccDNA in mitotically dividing cells. We believe that our model will provide a coherent framework for eccDNA research, to understand its biology and the factors that can be manipulated to modulate eccDNA load in eukaryotic cells.

Contingency and selection in mitochondrial genome dynamics

High frequencies of mutant mitochondrial DNA (mtDNA) in human cells lead to cellular defects that are associated with aging and disease. Yet much remains to be understood about the dynamics of the generation of mutant mtDNAs and their relative replicative fitness that informs their fate within cells and tissues. To address this, we utilize long-read single-molecule sequencing to track mutational trajectories of mtDNA in the model organism Saccharomyces cerevisiae. This model has numerous advantages over mammalian systems due to its much larger mtDNA and ease of artificially competing mutant and wild-type mtDNA copies in cells. We show a previously unseen pattern that constrains subsequent excision events in mtDNA fragmentation in yeast. We also provide evidence for the generation of rare and contentious non-periodic mtDNA structures that lead to persistent diversity within individual cells. Finally, we show that measurements of relative fitness of mtDNA fit a phenomenological model that highlights important biophysical parameters governing mtDNA fitness. Altogether, our study provides techniques and insights into the dynamics of large structural changes in genomes that we show are applicable to more complex organisms like humans.

Extrachromosomal circular DNA in cancer: history, current knowledge, and methods

Extrachromosomal circular DNA (eccDNA) is a closed-circle, nuclear, nonplasmid DNA molecule found in all tested eukaryotes. eccDNA plays important roles in cancer pathogenesis, evolution of tumor heterogeneity, and therapeutic resistance. It is known under many names, including very large cancer-specific circular extrachromosomal DNA (ecDNA), which carries oncogenes and is often amplified in cancer cells. Our understanding of eccDNA has historically been limited and fragmented. To provide better a context of new and previous research on eccDNA, in this review we give an overview of the various names given to eccDNA at different times. We describe the different mechanisms for formation of eccDNA and the methods used to study eccDNA thus far. Finally, we explore the potential clinical value of eccDNA.

Extrachromosomal Circular DNA: Category, Biogenesis, Recognition, and Functions

Extrachromosomal circular DNA (eccDNA), existing as double-stranded circular DNA, is derived and free from chromosomes. It is common in eukaryotes but has a strong heterogeneity in count, length, and origin. It has been demonstrated that eccDNA could function in telomere and rDNA maintenance, aging, drug resistance, tumorigenesis, and phenotypic variations of plants and animals. Here we review the current knowledge about eccDNA in category, biogenesis, recognition, and functions. We also provide perspectives on the potential implications of eccDNA in life science.

Recombination-mediated genome rearrangements

Homologous recombination (HR) is a universal DNA double-strand break (DSB) repair pathway that uses an intact DNA molecule as a template. Signature HR reactions are homology search and DNA strand invasion catalyzed by the prototypical RecA-ssDNA filament (Rad51 and Dmc1 in eukaryotes), which produces heteroduplex DNA-containing joint molecules (JMs). These reactions uniquely infringe on the DNA strands association established at replication, on the basis of substantial sequence similarity. For that reason, and despite the high fidelity of its templated nature, DSB repair by HR authorizes the alteration of genome structure, guided by repetitive DNA elements. The resulting structural variations (SVs) can involve vast genomic regions, potentially affecting multiple coding sequences and regulatory elements at once, with possible pathological consequences. Here, we discuss recent advances in our understanding of genetic and molecular vulnerabilities of HR leading to SVs, and of the various fidelity-enforcing factors acting across scales on the balancing act of this complex pathway. An emphasis is put on extra-chomosomal DNAs, both product of, and substrate for HR-mediated chromosomal rearrangements.

Anatomy and evolution of a DNA replication origin

DNA amplification occurs at the DNA puff II/9A locus in the fungus fly Sciara coprophila. As a foundation to study the molecular mechanism for the initiating events of II/9A DNA re-replication, we have sequenced 14 kb spanning a DNase hypersensitive site (DHS) upstream of the 1 kb amplification origin and through transcription units II/9-1 and II/9-2 downstream of the origin. These elements are annotated as well as the ORC binding site at the origin and the transition point (TP) between continuous and discontinuous DNA syntheses that marks the origin of bidirectional replication at the nucleotide level. A 9 bp motif found at the TP is repeated near the other end of the 1 kb ORI and may identify a putative second TP. The steroid hormone ecdysone induces DNA amplification as well as transcription and puffing at locus II/9A. Within the 14 kb, several matches to the ecdysone response element (EcRE) consensus sequence were identified, including some in the amplification origin region. EcRE O-P is at a central axis of a remarkable symmetry, equidistant to the TPs that are themselves equidistant to EcRE O-1 and EcRE O-2. DNA sequence alterations have occurred throughout the II/9A region in a newly discovered polymorphism (#2). Polymorphism #2 is not specific to developmental stage, sex, or tissue, and it does not impair DNA amplification. The DHS, both 9 bp TP sequences, and EcREs O-1, O-P, and O-2 are conserved between the polymorphism #1 and #2 sequences, suggesting their functional importance and retention during evolutionary selection. Moreover, a 72 bp sequence in the Sciara DHS at DNA puff II/9A is conserved in DNA puff C-3 of Rhynchosciara americana. Comparisons are discussed between the Sciara II/9A amplicon and the chorion locus amplicon on the third chromosome of Drosophila.

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.

Extrachromosomal DNA: Redefining the pathogenesis of glioblastoma

Glioblastoma is an incurable most prevalent primary malignant brain tumor in adults. Surgery followed by radiotherapy with concomitant chemotherapy is the standard of care in patients with glioblastoma. Although, prognosis remains poor with a median survival in the range of 12-15 months. Over the decades of research has identified the gene mutation, angiogenesis, cell signaling for the development novel therapeutics. However, recent understanding on extrachromosomal DNA (ecDNA) put extra-layer of complexity in glioblastoma pathogenesis. These ecDNAs are present in significantly higher copy number in the nucleus of the cancer cells and contains several oncogenes which are instrumental for intra-tumoral genetic heterogeneity, accelerated tumor evolution and therapy resistance. In this review, we will discuss the current understanding on biogenesis, disease progression and potential therapeutic implications of ecDNAs in glioblastoma.

Inverted duplicate DNA sequences increase translocation rates through sequencing nanopores resulting in reduced base calling accuracy

Inverted duplicated DNA sequences are a common feature of structural variants (SVs) and copy number variants (CNVs). Analysis of CNVs containing inverted duplicated DNA sequences using nanopore sequencing identified recurrent aberrant behavior characterized by low confidence, incorrect and missed base calls. Inverted duplicate DNA sequences in both yeast and human samples were observed to have systematic elevation in the electrical current detected at the nanopore, increased translocation rates and decreased sampling rates. The coincidence of inverted duplicated DNA sequences with dramatically reduced sequencing accuracy and an increased translocation rate suggests that secondary DNA structures may interfere with the dynamics of transit of the DNA through the nanopore.

Expandable and reversible copy number amplification drives rapid adaptation to antifungal drugs

Previously, we identified long repeat sequences that are frequently associated with genome rearrangements, including copy number variation (CNV), in many diverse isolates of the human fungal pathogen Candida albicans (Todd et al., 2019). Here, we describe the rapid acquisition of novel, high copy number CNVs during adaptation to azole antifungal drugs. Single-cell karyotype analysis indicates that these CNVs appear to arise via a dicentric chromosome intermediate and breakage-fusion-bridge cycles that are repaired using multiple distinct long inverted repeat sequences. Subsequent removal of the antifungal drug can lead to a dramatic loss of the CNV and reversion to the progenitor genotype and drug susceptibility phenotype. These findings support a novel mechanism for the rapid acquisition of antifungal drug resistance and provide genomic evidence for the heterogeneity frequently observed in clinical settings.

The adaptive potential of circular DNA accumulation in ageing cells

Carefully maintained and precisely inherited chromosomal DNA provides long-term genetic stability, but eukaryotic cells facing environmental challenges can benefit from the accumulation of less stable DNA species. Circular DNA molecules lacking centromeres segregate randomly or asymmetrically during cell division, following non-Mendelian inheritance patterns that result in high copy number instability and massive heterogeneity across populations. Such circular DNA species, variously known as extrachromosomal circular DNA (eccDNA), microDNA, double minutes or extrachromosomal DNA (ecDNA), are becoming recognised as a major source of the genetic variation exploited by cancer cells and pathogenic eukaryotes to acquire drug resistance. In budding yeast, circular DNA molecules derived from the ribosomal DNA (ERCs) have been long known to accumulate with age, but it is now clear that aged yeast also accumulate other high-copy protein-coding circular DNAs acquired through both random and environmentally-stimulated recombination processes. Here, we argue that accumulation of circular DNA provides a reservoir of heterogeneous genetic material that can allow rapid adaptation of aged cells to environmental insults, but avoids the negative fitness impacts on normal growth of unsolicited gene amplification in the young population.

Transcription-induced formation of extrachromosomal DNA during yeast ageing

Extrachromosomal circular DNA (eccDNA) facilitates adaptive evolution by allowing rapid and extensive gene copy number variation and is implicated in the pathology of cancer and ageing. Here, we demonstrate that yeast aged under environmental copper accumulate high levels of eccDNA containing the copper-resistance gene CUP1. Transcription of the tandemly repeated CUP1 gene causes CUP1 eccDNA accumulation, which occurs in the absence of phenotypic selection. We have developed a sensitive and quantitative eccDNA sequencing pipeline that reveals CUP1 eccDNA accumulation on copper exposure to be exquisitely site specific, with no other detectable changes across the eccDNA complement. eccDNA forms de novo from the CUP1 locus through processing of DNA double-strand breaks (DSBs) by Sae2, Mre11 and Mus81, and genome-wide analyses show that other protein coding eccDNA species in aged yeast share a similar biogenesis pathway. Although abundant, we find that CUP1 eccDNA does not replicate efficiently, and high-copy numbers in aged cells arise through frequent formation events combined with asymmetric DNA segregation. The transcriptional stimulation of CUP1 eccDNA formation shows that age-linked genetic change varies with transcription pattern, resulting in gene copy number profiles tailored by environment.

Transcription can cause the de novo formation of protein-coding extrachromosomal DNA that accumulates in ageing yeast cells; these extrachromosomal circular DNA molecules form frequently by a DNA double strand break repair mechanism.

Fitness benefits of loss of heterozygosity in Saccharomyces hybrids

With two genomes in the same organism, interspecific hybrids have unique fitness opportunities and costs. In both plants and yeasts, wild, pathogenic, and domesticated hybrids may eliminate portions of one parental genome, a phenomenon known as loss of heterozygosity (LOH). Laboratory evolution of hybrid yeast recapitulates these results, with LOH occurring in just a few hundred generations of propagation. In this study, we systematically looked for alleles that are beneficial when lost in order to determine how prevalent this mode of adaptation may be and to determine candidate loci that might underlie the benefits of larger-scale chromosome rearrangements. These aims were accomplished by mating Saccharomyces uvarum with the S. cerevisiae deletion collection to create hybrids such that each nonessential S. cerevisiae allele is deleted. Competitive fitness assays of these pooled, barcoded, hemizygous strains, and accompanying controls, revealed a large number of loci for which LOH is beneficial. We found that the fitness effects of hemizygosity are dependent on the species context, the selective environment, and the species origin of the deleted allele. Further, we found that hybrids have a wider distribution of fitness consequences versus matched S. cerevisiae hemizygous diploids. Our results suggest that LOH can be a successful strategy for adaptation of hybrids to new environments, and we identify candidate loci that drive the chromosomal rearrangements observed in evolution of yeast hybrids.

Phenotypic and Genotypic Consequences of CRISPR/Cas9 Editing of the Replication Origins in the rDNA of Saccharomyces cerevisiae

The complex structure and repetitive nature of eukaryotic ribosomal DNA (rDNA) is a challenge for genome assembly, thus the consequences of sequence variation in rDNA remain unexplored. However, renewed interest in the role that rDNA variation may play in diverse cellular functions, aside from ribosome production, highlights the need for a method that would permit genetic manipulation of the rDNA. Here, we describe a clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-based strategy to edit the rDNA locus in the budding yeast Saccharomyces cerevisiae, developed independently but similar to one developed by others. Using this approach, we modified the endogenous rDNA origin of replication in each repeat by deleting or replacing its consensus sequence. We characterized the transformants that have successfully modified their rDNA locus and propose a mechanism for how CRISPR/Cas9-mediated editing of the rDNA occurs. In addition, we carried out extended growth and life span experiments to investigate the long-term consequences that altering the rDNA origin of replication have on cellular health. We find that long-term growth of the edited clones results in faster-growing suppressors that have acquired segmental aneusomy of the rDNA-containing region of chromosome XII or aneuploidy of chromosomes XII, II, or IV. Furthermore, we find that all edited isolates suffer a reduced life span, irrespective of their levels of extrachromosomal rDNA circles. Our work demonstrates that it is possible to quickly, efficiently, and homogeneously edit the rDNA origin via CRISPR/Cas9.

Genome plasticity in Candida albicans is driven by long repeat sequences

Genome rearrangements resulting in copy number variation (CNV) and loss of heterozygosity (LOH) are frequently observed during the somatic evolution of cancer and promote rapid adaptation of fungi to novel environments. In the human fungal pathogen Candida albicans, CNV and LOH confer increased virulence and antifungal drug resistance, yet the mechanisms driving these rearrangements are not completely understood. Here, we unveil an extensive array of long repeat sequences (65-6499 bp) that are associated with CNV, LOH, and chromosomal inversions. Many of these long repeat sequences are uncharacterized and encompass one or more coding sequences that are actively transcribed. Repeats associated with genome rearrangements are predominantly inverted and separated by up to ~1.6 Mb, an extraordinary distance for homology-based DNA repair/recombination in yeast. These repeat sequences are a significant source of genome plasticity across diverse strain backgrounds including clinical, environmental, and experimentally evolved isolates, and represent previously uncharacterized variation in the reference genome.

An evolving view of copy number variants

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

De novo origination of MIRNAs through generation of short inverted repeats in target genes

MIRNA (MIR) gene origin and early evolutionary processes, such as hairpin precursor sequence origination, promoter activity acquirement and the sequence of these two processes, are fundamental and fascinating subjects. Three models, including inverted gene duplication, spontaneous evolution and transposon transposition, have been proposed for de novo origination of hairpin precursor sequence. However, these models still open to discussion. In addition, de novo origination of MIR gene promoters has not been well investigated. Here, I systematically investigated the origin of evolutionarily young polyphenol oxidase gene (PPO)-targeting MIRs, including MIR1444, MIR058 and MIR12112, and a genomic region termed AasPPO-as-hp, which contained a hairpin-forming sequence. I found that MIR058 precursors and the hairpin-forming sequence of AasPPO-as-hp originated in an ancient PPO gene through forming short inverted repeats. Palindromic-like sequences and imperfect inverted repeats in the ancient PPO gene contributed to initiate the generation of short inverted repeats probably by causing errors during DNA duplication. Analysis of MIR058 and AasPPO-as-hp promoters showed that they originated in the 3'-flanking region of the ancient PPO gene. Promoter activities were gained by insertion of a CAAT-box and multiple-copper-response element (CuRE)-containing miniature inverted-repeat transposable element (MITE) in the upstream of AT-rich TATA-box-like sequence. Gain of promoter activities occurred before hairpin-forming sequence origination. Sequence comparison of MIR1444, MIR058 and MIR12112 promoters showed frequent birth and death of CuREs, indicating copper could be vital for the origination and evolution of PPO-targeting MIRs. Based on the evidence obtained, a novel model for plant MIR origination and evolution is proposed.

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.

Discoveries of Extrachromosomal Circles of DNA in Normal and Tumor Cells

While the vast majority of cellular DNA in eukaryotes is contained in long linear strands in chromosomes, we have long recognized some exceptions like mitochondrial DNA, plasmids in yeasts, and double minutes (DMs) in cancer cells where the DNA is present in extrachromosomal circles. In addition, specialized extrachromosomal circles of DNA (eccDNA) have been noted to arise from repetitive genomic sequences like telomeric DNA or rDNA. Recently eccDNA arising from unique (nonrepetitive) DNA have been discovered in normal and malignant cells, raising interesting questions about their biogenesis, function and clinical utility. Here, we review recent results and future directions of inquiry on these new forms of eccDNA.

Loss of Heterozygosity Drives Adaptation in Hybrid Yeast

Hybridization is often considered maladaptive, but sometimes hybrids can invade new ecological niches and adapt to novel or stressful environments better than their parents. The genomic changes that occur following hybridization that facilitate genome resolution and/or adaptation are not well understood. Here, we examine hybrid genome evolution using experimental evolution of de novo interspecific hybrid yeast Saccharomyces cerevisiae × Saccharomyces uvarum and their parentals. We evolved these strains in nutrient-limited conditions for hundreds of generations and sequenced the resulting cultures identifying numerous point mutations, copy number changes, and loss of heterozygosity (LOH) events, including species-biased amplification of nutrient transporters. We focused on a particularly interesting example, in which we saw repeated LOH at the high-affinity phosphate transporter gene PHO84 in both intra- and interspecific hybrids. Using allele replacement methods, we tested the fitness of different alleles in hybrid and S. cerevisiae strain backgrounds and found that the LOH is indeed the result of selection on one allele over the other in both S. cerevisiae and the hybrids. This is an example where hybrid genome resolution is driven by positive selection on existing heterozygosity and demonstrates that even infrequent outcrossing may have lasting impacts on adaptation.

Differential paralog divergence modulates genome evolution across yeast species

Evolutionary outcomes depend not only on the selective forces acting upon a species, but also on the genetic background. However, large timescales and uncertain historical selection pressures can make it difficult to discern such important background differences between species. Experimental evolution is one tool to compare evolutionary potential of known genotypes in a controlled environment. Here we utilized a highly reproducible evolutionary adaptation in Saccharomyces cerevisiae to investigate whether experimental evolution of other yeast species would select for similar adaptive mutations. We evolved populations of S. cerevisiae, S. paradoxus, S. mikatae, S. uvarum, and interspecific hybrids between S. uvarum and S. cerevisiae for ~200–500 generations in sulfate-limited continuous culture. Wild-type S. cerevisiae cultures invariably amplify the high affinity sulfate transporter gene, SUL1. However, while amplification of the SUL1 locus was detected in S. paradoxus and S. mikatae populations, S. uvarum cultures instead selected for amplification of the paralog, SUL2. We measured the relative fitness of strains bearing deletions and amplifications of both SUL genes from different species, confirming that, converse to S. cerevisiae, S. uvarum SUL2 contributes more to fitness in sulfate limitation than S. uvarum SUL1. By measuring the fitness and gene expression of chimeric promoter-ORF constructs, we were able to delineate the cause of this differential fitness effect primarily to the promoter of S. uvarum SUL1. Our data show evidence of differential sub-functionalization among the sulfate transporters across Saccharomyces species through recent changes in noncoding sequence. Furthermore, these results show a clear example of how such background differences due to paralog divergence can drive changes in genome evolution.

Both comparative genomics and experimental evolution are powerful tools that can be used to make inferences about evolutionary processes. Together, these approaches provide the opportunity to observe evolutionary adaptation over millions of years where selective history is largely unknown, and over short timescales under controlled selective pressures in the laboratory. We have used comparative experimental evolution to observe the evolutionary fate of an adaptive mutation, and determined to what degree the outcome is conditional on the genetic background. We evolved several populations of different yeast species for over 200 generations in sulfate-limited conditions to determine how the differences in genomic context can alter evolutionary routes when challenged with a nutrient limitation selection pressure. We find that the gene encoding a high affinity sulfur transporter becomes amplified in most species of Saccharomyces, except in S. uvarum, in which the amplification of the paralogous sulfate transporter gene SUL2 is recovered. We attribute this change in amplification preference to mutations in the non-coding region of SUL1, likely due to reduced expression of this gene in S. uvarum. We conclude that the adaptive mutations selected for in each organism depend on the genomic context, even when faced with the same environmental condition.

Ontogeny of Unstable Chromosomes Generated by Telomere Error in Budding Yeast

DNA replication errors at certain sites in the genome initiate chromosome instability that ultimately leads to stable genomic rearrangements. Where instability begins is often unclear. And, early instability may form unstable chromosome intermediates whose transient nature also hinders mechanistic understanding. We report here a budding yeast model that reveals the genetic ontogeny of genome rearrangements, from initial replication error to unstable chromosome formation to their resolution. Remarkably, the initial error often arises in or near the telomere, and frequently forms unstable chromosomes. Early unstable chromosomes may then resolve to an internal "collection site" where a dicentric forms and resolves to an isochromosome (other outcomes are possible at each step). The initial telomere-proximal unstable chromosome is increased in mutants in telomerase subunits, Tel1, and even Rad9, with no known telomere-specific function. Defects in Tel1 and in Rrm3, a checkpoint protein kinase with a role in telomere maintenance and a DNA helicase, respectively, synergize dramatically to generate unstable chromosomes, further illustrating the consequence of replication error in the telomere. Collectively, our results suggest telomeric replication errors may be a common cause of seemingly unrelated genomic rearrangements located hundreds of kilobases away.