Quantitative Analysis of the Rates for Repeat-Mediated Genome Instability in a Yeast Experimental System

Eukaryotic Pif1 helicase unwinds G-quadruplex and dsDNA using a conserved wedge

G-quadruplexes (G4s) formed by guanine-rich nucleic acids induce genome instability through impeding DNA replication fork progression. G4s are stable DNA structures, the unfolding of which require the functions of DNA helicases. Pif1 helicase binds preferentially to G4 DNA and plays multiple roles in maintaining genome stability, but the mechanism by which Pif1 unfolds G4s is poorly understood. Here we report the co-crystal structure of Saccharomyces cerevisiae Pif1 (ScPif1) bound to a G4 DNA with a 5' single-stranded DNA (ssDNA) segment. Unlike the Thermus oshimai Pif1-G4 structure, in which the 1B and 2B domains confer G4 recognition, ScPif1 recognizes G4 mainly through the wedge region in the 1A domain that contacts the 5' most G-tetrad directly. A conserved Arg residue in the wedge is required for Okazaki fragment processing but not for mitochondrial function or for suppression of gross chromosomal rearrangements. Multiple substitutions at this position have similar effects on resolution of DNA duplexes and G4s, suggesting that ScPif1 may use the same wedge to unwind G4 and dsDNA. Our results reveal the mechanism governing dsDNA unwinding and G4 unfolding by ScPif1 helicase that can potentially be generalized to other eukaryotic Pif1 helicases and beyond.

Two residues in the DNA binding site of Pif1 helicase are essential for nuclear functions but dispensable for mitochondrial respiratory growth

Pif1 helicase functions in both the nucleus and mitochondria. Pif1 tightly couples ATP hydrolysis, single-stranded DNA translocation, and duplex DNA unwinding. We investigated two Pif1 variants (F723A and T464A) that have each lost one site of interaction of the protein with the DNA substrate. Both variants exhibit minor reductions in affinity for DNA and ATP hydrolysis but have impaired DNA unwinding activity. However, these variants translocate on single-stranded DNA faster than the wildtype enzyme and can slide on the DNA substrate in an ATP-independent manner. This suggests they have lost their grip on the DNA, interfering with coupling ATP hydrolysis to translocation and unwinding. Yeast expressing these variants have increased gross chromosomal rearrangements, increased telomere length, and can overcome the lethality of dna2Δ, similar to phenotypes of yeast lacking Pif1. However, unlike pif1Δ mutants, they are viable on glycerol containing media and maintain similar mitochondrial DNA copy numbers as Pif1 wildtype. Overall, our data indicate that a tight grip of the trailing edge of the Pif1 enzyme on the DNA couples ATP hydrolysis to DNA translocation and DNA unwinding. This tight grip appears to be essential for the Pif1 nuclear functions tested but is dispensable for mitochondrial respiratory growth.

Massive contractions of myotonic dystrophy type 2-associated CCTG tetranucleotide repeats occur via double-strand break repair with distinct requirements for DNA helicases

Myotonic dystrophy type 2 (DM2) is a genetic disease caused by expanded CCTG DNA repeats in the first intron of CNBP. The number of CCTG repeats in DM2 patients ranges from 75 to 11,000, yet little is known about the molecular mechanisms responsible for repeat expansions or contractions. We developed an experimental system in Saccharomyces cerevisiae that enables the selection of large-scale contractions of (CCTG)100 within the intron of a reporter gene and subsequent genetic analysis. Contractions exceeded 80 repeat units, causing the final repetitive tract to be well below the threshold for disease. We found that Rad51 and Rad52 are involved in these massive contractions, indicating a mechanism that uses homologous recombination. Srs2 helicase was shown previously to stabilize CTG, CAG, and CGG repeats. Loss of Srs2 did not significantly affect CCTG contraction rates in unperturbed conditions. In contrast, loss of the RecQ helicase Sgs1 resulted in a 6-fold decrease in contraction rate with specific evidence that helicase activity is required for large-scale contractions. Using a genetic assay to evaluate chromosome arm loss, we determined that CCTG and reverse complementary CAGG repeats elevate the rate of chromosomal fragility compared to a short-track control. Overall, our results demonstrate that the genetic control of CCTG repeat contractions is notably distinct among disease-causing microsatellite repeat sequences.

Elevated MSH2 MSH3 expression interferes with DNA metabolism in vivo

The Msh2-Msh3 mismatch repair (MMR) complex in Saccharomyces cerevisiae recognizes and directs repair of insertion/deletion loops (IDLs) up to ∼17 nucleotides. Msh2-Msh3 also recognizes and binds distinct looped and branched DNA structures with varying affinities, thereby contributing to genome stability outside post-replicative MMR through homologous recombination, double-strand break repair (DSBR) and the DNA damage response. In contrast, Msh2-Msh3 promotes genome instability through trinucleotide repeat (TNR) expansions, presumably by binding structures that form from single-stranded (ss) TNR sequences. We previously demonstrated that Msh2-Msh3 binding to 5' ssDNA flap structures interfered with Rad27 (Fen1 in humans)-mediated Okazaki fragment maturation (OFM) in vitro. Here we demonstrate that elevated Msh2-Msh3 levels interfere with DNA replication and base excision repair in vivo. Elevated Msh2-Msh3 also induced a cell cycle arrest that was dependent on RAD9 and ELG1 and led to PCNA modification. These phenotypes also required Msh2-Msh3 ATPase activity and downstream MMR proteins, indicating an active mechanism that is not simply a result of Msh2-Msh3 DNA-binding activity. This study provides new mechanistic details regarding how excess Msh2-Msh3 can disrupt DNA replication and repair and highlights the role of Msh2-Msh3 protein abundance in Msh2-Msh3-mediated genomic instability.

Large-scale expansions of Friedreich's ataxia GAA•TTC repeats in an experimental human system: role of DNA replication and prevention by LNA-DNA oligonucleotides and PNA oligomers

Friedreich's ataxia (FRDA) is caused by expansions of GAA•TTC repeats in the first intron of the human FXN gene that occur during both intergenerational transmissions and in somatic cells. Here we describe an experimental system to analyze large-scale repeat expansions in cultured human cells. It employs a shuttle plasmid that can replicate from the SV40 origin in human cells or be stably maintained in S. cerevisiae utilizing ARS4-CEN6. It also contains a selectable cassette allowing us to detect repeat expansions that accumulated in human cells upon plasmid transformation into yeast. We indeed observed massive expansions of GAA•TTC repeats, making it the first genetically tractable experimental system to study large-scale repeat expansions in human cells. Further, GAA•TTC repeats stall replication fork progression, while the frequency of repeat expansions appears to depend on proteins implicated in replication fork stalling, reversal, and restart. Locked nucleic acid (LNA)-DNA mixmer oligonucleotides and peptide nucleic acid (PNA) oligomers, which interfere with triplex formation at GAA•TTC repeats in vitro, prevented the expansion of these repeats in human cells. We hypothesize, therefore, that triplex formation by GAA•TTC repeats stall replication fork progression, ultimately leading to repeat expansions during replication fork restart.

Efficient, robust, and versatile fluctuation data analysis using MLE MUtation Rate calculator (mlemur)

The fluctuation assay remains an important tool for analyzing the levels of mutagenesis in microbial populations. The mutant counts originating from some average number of mutations are usually assumed to obey the Luria-Delbrück distribution. While several tools for estimating mutation rates are available, they sometimes lack accuracy or versatility under non-standard conditions. In this work, extensions to the Luria-Delbrück protocol to account for phenotypic lag and cellular death with either perfect or partial plating were developed. Hence, the novel MLE MUtation Rate calculator, or mlemur, is the first tool that provides a user-friendly graphical interface allowing the researchers to model their data with consideration for partial plating, differential growth of mutants and non-mutants, phenotypic lag, cellular death, variability of the final number of cells, post-exponential-phase mutations, and the size of the inoculum. Additionally, mlemur allows the users to incorporate most of these special conditions at the same time to obtain highly accurate estimates of mutation rates and P values, confidence intervals for an arbitrary function of data (such as fold), and perform power analysis and sample size determination for the likelihood ratio test. The accuracy of point and interval estimates produced by mlemur against historical and simulated fluctuation experiments are assessed. Both mlemur and the analyses in this work might be of great help when evaluating fluctuation experiments and increase the awareness of the limitations of the widely-used Lea-Coulson formulation of the Luria-Delbrück distribution in the more realistic biological contexts.

The RNA export and RNA decay complexes THO and TRAMP prevent transcription-replication conflicts, DNA breaks, and CAG repeat contractions

Expansion of structure-forming CAG/CTG repetitive sequences is the cause of several neurodegenerative disorders and deletion of repeats is a potential therapeutic strategy. Transcription-associated mechanisms are known to cause CAG repeat instability. In this study, we discovered that Thp2, an RNA export factor and member of the THO (suppressors of transcriptional defects of hpr1Δ by overexpression) complex, and Trf4, a key component of the TRAMP (Trf4/5-Air1/2-Mtr4 polyadenylation) complex involved in nuclear RNA polyadenylation and degradation, are necessary to prevent CAG fragility and repeat contractions in a Saccharomyces cerevisiae model system. Depletion of both Thp2 and Trf4 proteins causes a highly synergistic increase in CAG repeat fragility, indicating a complementary role of the THO and TRAMP complexes in preventing genome instability. Loss of either Thp2 or Trf4 causes an increase in RNA polymerase stalling at the CAG repeats and other genomic loci, as well as genome-wide transcription-replication conflicts (TRCs), implicating TRCs as a cause of CAG fragility and instability in their absence. Analysis of the effect of RNase H1 overexpression on CAG fragility, RNAPII stalling, and TRCs suggests that RNAPII stalling with associated R-loops are the main cause of CAG fragility in the thp2Δ mutants. In contrast, CAG fragility and TRCs in the trf4Δ mutant can be compensated for by RPA overexpression, suggesting that excess unprocessed RNA in TRAMP4 mutants leads to reduced RPA availability and high levels of TRCs. Our results show the importance of RNA surveillance pathways in preventing RNAPII stalling, TRCs, and DNA breaks, and show that RNA export and RNA decay factors work collaboratively to maintain genome stability.

Uncontrolled transposition following RNAi loss causes hypermutation and antifungal drug resistance in clinical isolates of Cryptococcus neoformans

Cryptococcus neoformans infections cause approximately 15% of AIDS-related deaths owing to a combination of limited antifungal therapies and drug resistance. A collection of clinical and environmental C. neoformans isolates were assayed for increased mutation rates via fluctuation analysis, and we identified two hypermutator C. neoformans clinical isolates with increased mutation rates when exposed to the combination of rapamycin and FK506. Sequencing of drug target genes found that Cnl1 transposon insertions conferred the majority of resistance to rapamycin and FK506 and could also independently cause resistance to 5-fluoroorotic acid and the clinically relevant antifungal 5-flucytosine. Whole-genome sequencing revealed both hypermutator genomes harbour a nonsense mutation in the RNA-interference component ZNF3 and hundreds of Cnl1 elements organized into massive subtelomeric arrays on each of the fourteen chromosomes. Quantitative trait locus mapping in 28 progeny derived from a cross between a hypermutator and wild-type identified a locus associated with hypermutation that included znf3. CRISPR editing of the znf3 nonsense mutation abolished hypermutation and restored small-interfering-RNA production. We conclude that hypermutation and drug resistance in these clinical isolates result from RNA-interference loss and accumulation of Cnl1 elements.

Bursts of Genomic Instability Potentiate Phenotypic and Genomic Diversification in Saccharomyces cerevisiae

How microbial cells leverage their phenotypic potential to survive in a changing environment is a complex biological problem, with important implications for pathogenesis and species evolution. Stochastic phenotype switching, a particularly fascinating adaptive approach observed in numerous species across the tree of life, introduces phenotypic diversity into a population through mechanisms which have remained difficult to define. Here we describe our investigations into the mechanistic basis of colony morphology phenotype switching which occurs in populations of a pathogenic isolate of Saccharomyces cerevisiae, YJM311. We observed that clonal populations of YJM311 cells produce variant colonies that display altered morphologies and, using whole genome sequence analysis, discovered that these variant clones harbored an exceptional collection of karyotypes newly altered by de novo structural genomic variations (SVs). Overall, our analyses indicate that copy number alterations, more often than changes in allelic identity, provide the causative basis of this phenotypic variation. Individual variants carried between 1 and 16 de novo copy number variations, most of which were whole chromosomal aneuploidies. Notably, we found that the inherent stability of the diploid YJM311 genome is comparable to that of domesticated laboratory strains, indicating that the collections of SVs harbored by variant clones did not arise by a chronic chromosomal instability (CIN) mechanism. Rather, our data indicate that these variant clones acquired such complex karyotypic configurations simultaneously, during stochastic and transient episodes of punctuated systemic genomic instability (PSGI). Surprisingly, we found that the majority of these highly altered variant karyotypes were propagated with perfect fidelity in long-term passaging experiments, demonstrating that high aneuploidy burdens can often be conducive with prolonged genomic integrity. Together, our results demonstrate that colony morphology switching in YJM311 is driven by a stochastic process in which genome stability and plasticity are integrally coupled to phenotypic heterogeneity. Consequently, this system simultaneously introduces both phenotypic and genomic variation into a population of cells, which can, in turn perpetuate population diversity for many generations thereafter.

Complex mutation profiles in mismatch repair and ribonucleotide reductase mutants reveal novel repair substrate specificity of MutS homolog (MSH) complexes

Determining mutation signatures is standard for understanding the etiology of human tumors and informing cancer treatment. Multiple determinants of DNA replication fidelity prevent mutagenesis that leads to carcinogenesis, including the regulation of free deoxyribonucleoside triphosphate pools by ribonucleotide reductase and repair of replication errors by the mismatch repair system. We identified genetic interactions between rnr1 alleles that skew and/or elevate deoxyribonucleoside triphosphate levels and mismatch repair gene deletions. These defects indicate that the rnr1 alleles lead to increased mutation loads that are normally acted upon by mismatch repair. We then utilized a targeted deep-sequencing approach to determine mutational profiles associated with mismatch repair pathway defects. By combining rnr1 and msh mutations to alter and/or increase deoxyribonucleoside triphosphate levels and alter the mutational load, we uncovered previously unreported specificities of Msh2-Msh3 and Msh2-Msh6. Msh2-Msh3 is uniquely able to direct the repair of G/C single-base deletions in GC runs, while Msh2-Msh6 specifically directs the repair of substitutions that occur at G/C dinucleotides. We also identified broader sequence contexts that influence variant profiles in different genetic backgrounds. Finally, we observed that the mutation profiles in double mutants were not necessarily an additive relationship of mutation profiles in single mutants. Our results have implications for interpreting mutation signatures from human tumors, particularly when mismatch repair is defective.

Partners in crime: Tbf1 and Vid22 promote expansions of long human telomeric repeats at an interstitial chromosome position in yeast

In humans, telomeric repeats (TTAGGG)_n are known to be present at internal chromosomal sites. These interstitial telomeric sequences (ITSs) are an important source of genomic instability, including repeat length polymorphism, but the molecular mechanisms responsible for this instability remain to be understood. Here, we studied the mechanisms responsible for expansions of human telomeric (Htel) repeats that were artificially inserted inside a yeast chromosome. We found that Htel repeats in an interstitial chromosome position are prone to expansions. The propensity of Htel repeats to expand depends on the presence of a complex of two yeast proteins: Tbf1 and Vid22. These two proteins are physically bound to an interstitial Htel repeat, and together they slow replication fork progression through it. We propose that slow progression of the replication fork through the protein complex formed by the Tbf1 and Vid22 partners at the Htel repeat cause DNA strand slippage, ultimately resulting in repeat expansions.

Genome-scale patterns in the loss of heterozygosity incidence in Saccharomyces cerevisiae

Former studies have established that loss of heterozygosity can be a key driver of sequence evolution in unicellular eukaryotes and tissues of metazoans. However, little is known about whether the distribution of loss of heterozygosity events is largely random or forms discernible patterns across genomes. To initiate our experiments, we introduced selectable markers to both arms of all chromosomes of the budding yeast. Subsequent extensive assays, repeated over several genetic backgrounds and environments, provided a wealth of information on the genetic and environmental determinants of loss of heterozygosity. Three findings stand out. First, the number of loss of heterozygosity events per unit time was more than 25 times higher for growing than starving cells. Second, loss of heterozygosity was most frequent when regions of homology around a recombination site were identical, about a half-% sequence divergence was sufficient to reduce its incidence. Finally, the density of loss of heterozygosity events was highly dependent on the genome's physical architecture. It was several-fold higher on short chromosomal arms than on long ones. Comparably large differences were seen within a single arm where regions close to a centromere were visibly less affected than regions close, though usually not strictly adjacent, to a telomere. We suggest that the observed uneven distribution of loss of heterozygosity events could have been caused not only by an uneven density of initial DNA damages. Location-depended differences in the mode of DNA repair, or its effect on fitness, were likely to operate as well.

VID22 counteracts G-quadruplex-induced genome instability

Genome instability is a condition characterized by the accumulation of genetic alterations and is a hallmark of cancer cells. To uncover new genes and cellular pathways affecting endogenous DNA damage and genome integrity, we exploited a Synthetic Genetic Array (SGA)-based screen in yeast. Among the positive genes, we identified VID22, reported to be involved in DNA double-strand break repair. vid22Δ cells exhibit increased levels of endogenous DNA damage, chronic DNA damage response activation and accumulate DNA aberrations in sequences displaying high probabilities of forming G-quadruplexes (G4-DNA). If not resolved, these DNA secondary structures can block the progression of both DNA and RNA polymerases and correlate with chromosome fragile sites. Vid22 binds to and protects DNA at G4-containing regions both in vitro and in vivo. Loss of VID22 causes an increase in gross chromosomal rearrangement (GCR) events dependent on G-quadruplex forming sequences. Moreover, the absence of Vid22 causes defects in the correct maintenance of G4-DNA rich elements, such as telomeres and mtDNA, and hypersensitivity to the G4-stabilizing ligand TMPyP4. We thus propose that Vid22 is directly involved in genome integrity maintenance as a novel regulator of G4 metabolism.

Mapping single-cell-resolution cell phylogeny reveals cell population dynamics during organ development

Mapping the cell phylogeny of a complex multicellular organism relies on somatic mutations accumulated from zygote to adult. Available cell barcoding methods can record about three mutations per barcode, enabling only low-resolution mapping of the cell phylogeny of complex organisms. Here we developed SMALT, a substitution mutation-aided lineage-tracing system that outperforms the available cell barcoding methods in mapping cell phylogeny. We applied SMALT to Drosophila melanogaster and obtained on average more than 20 mutations on a three-kilobase-pair barcoding sequence in early-adult cells. Using the barcoding mutations, we obtained high-quality cell phylogenetic trees, each comprising several thousand internal nodes with 84-93% median bootstrap support. The obtained cell phylogenies enabled a population genetic analysis that estimates the longitudinal dynamics of the number of actively dividing parental cells (Np) in each organ through development. The Np dynamics revealed the trajectory of cell births and provided insight into the balance of symmetric and asymmetric cell division.

Rad9-mediated checkpoint activation is responsible for elevated expansions of GAA repeats in CST-deficient yeast

Large-scale expansion of (GAA)n repeats in the first intron of the FXN gene is responsible for the severe neurodegenerative disease, Friedreich's ataxia in humans. We have previously conducted an unbiased genetic screen for GAA repeat instability in a yeast experimental system. The majority of genes that came from this screen encoded the components of DNA replication machinery, strongly implying that replication irregularities are at the heart of GAA repeat expansions. This screen, however, also produced two unexpected hits: members of the CST complex, CDC13 and TEN1 genes, which are required for telomere maintenance. To understand how the CST complex could affect intra-chromosomal GAA repeats, we studied the well-characterized temperature-sensitive cdc13-1 mutation and its effects on GAA repeat instability in yeast. We found that in-line with the screen results, this mutation leads to ∼10-fold increase in the rate of large-scale expansions of the (GAA)100 repeat at semi-permissive temperature. Unexpectedly, the hyper-expansion phenotype of the cdc13-1 mutant largely depends on activation of the G2/M checkpoint, as deletions of individual genes RAD9, MEC1, RAD53, and EXO1 belonging to this pathway rescued the increased GAA expansions. Furthermore, the hyper-expansion phenotype of the cdc13-1 mutant depended on the subunit of DNA polymerase δ, Pol32. We hypothesize, therefore, that increased repeat expansions in the cdc13-1 mutant happen during post-replicative repair of nicks or small gaps within repetitive tracts during the G2 phase of the cell cycle upon activation of the G2/M checkpoint.

A selection-based next generation sequencing approach to develop robust, genotype-specific mutation profiles in Saccharomyces cerevisiae

Distinct mutation signatures arise from environmental exposures and/or from defects in metabolic pathways that promote genome stability. The presence of a particular mutation signature can therefore predict the underlying mechanism of mutagenesis. These insults to the genome often alter dNTP pools, which itself impacts replication fidelity. Therefore, the impact of altered dNTP pools should be considered when making mechanistic predictions based on mutation signatures. We developed a targeted deep-sequencing approach on the CAN1 gene in Saccharomyces cerevisiae to define information-rich mutational profiles associated with distinct rnr1 backgrounds. Mutations in the activity and selectivity sites of rnr1 lead to elevated and/or unbalanced dNTP levels, which compromises replication fidelity and increases mutation rates. The mutation spectra of rnr1Y285F and rnr1Y285A alleles were characterized previously; our analysis was consistent with this prior work but the sequencing depth achieved in our study allowed a significantly more robust and nuanced computational analysis of the variants observed, generating profiles that integrated information about mutation spectra, position effects, and sequence context. This approach revealed previously unidentified, genotype-specific mutation profiles in the presence of even modest changes in dNTP pools. Furthermore, we identified broader sequence contexts and nucleotide motifs that influenced variant profiles in different rnr1 backgrounds, which allowed specific mechanistic predictions about the impact of altered dNTP pools on replication fidelity.

Relocation of Collapsed Forks to the Nuclear Pore Complex Depends on Sumoylation of DNA Repair Proteins and Permits Rad51 Association

Expanded CAG repeats form stem-loop secondary structures that lead to fork stalling and collapse. Previous work has shown that these collapsed forks relocalize to nuclear pore complexes (NPCs) in late S phase in a manner dependent on replication, the nucleoporin Nup84, and the Slx5 protein, which prevents repeat fragility and instability. Here, we show that binding of the Smc5/6 complex to the collapsed fork triggers Mms21-dependent sumoylation of fork-associated DNA repair proteins, and that RPA, Rad52, and Rad59 are the key sumoylation targets that mediate relocation. The SUMO interacting motifs of Slx5 target collapsed forks to the NPC. Notably, Rad51 foci only co-localize with the repeat after it is anchored to the nuclear periphery and Rad51 exclusion from the early collapsed fork is dependent on RPA sumoylation. This pathway may provide a mechanism to constrain recombination at stalled or collapsed forks until it is required for fork restart.

webSalvador: a Web Tool for the Luria-Delbrük Experiment

Existing Web tools for the Luria-Delbrück fluctuation experiment do not offer many desirable capabilities that are vital to mutation research. webSalvador offers these capabilities via a user interface that allows researchers to access most of the functions in the R package rSalvador without having to learn the R language.

Novel mechanisms of TolC-independent decreased bile-salt susceptibility in Escherichia coli

Bile salts, including sodium deoxycholate (DOC), are secreted into the intestine to aid fat digestion and contribute to antimicrobial protection. Gram-negative pathogens such as Escherichia coli, however, are highly resistant to DOC, using multiple mechanisms of which the multidrug efflux pump AcrAB-TolC is the dominant one. Given that TolC-mediated efflux masks the interaction of DOC with potential targets, we sought to identify those targets by identifying genes whose mutations cause an increase in the MIC to DOC relative to the ∆tolC parental strain, that lacks TolC-associated functional efflux pumps. Using a mutant screen, we isolated twenty independent spontaneous mutants that had a higher MICDOC than the E. coli parental ∆tolC strain. Whole genome sequencing of these mutants mapped most mutations to the ptsI or cyaA gene. Analysis of knock-out mutants and complementation showed that elimination of PtsI, a component of the carbohydrate phosphotransferase system, or one of the two key proteins involved in cAMP synthesis and signaling, adenylate cyclase (CyaA) or cAMP receptor protein (Crp) causes low-level increased resistance of a ∆tolC E. coli strain to DOC.

Phenotypic and molecular evolution across 10,000 generations in laboratory budding yeast populations

Laboratory experimental evolution provides a window into the details of the evolutionary process. To investigate the consequences of long-term adaptation, we evolved 205 Saccharomyces cerevisiae populations (124 haploid and 81 diploid) for ~10,000,000 generations in three environments. We measured the dynamics of fitness changes over time, finding repeatable patterns of declining adaptability. Sequencing revealed that this phenotypic adaptation is coupled with a steady accumulation of mutations, widespread genetic parallelism, and historical contingency. In contrast to long-term evolution in E. coli, we do not observe long-term coexistence or populations with highly elevated mutation rates. We find that evolution in diploid populations involves both fixation of heterozygous mutations and frequent loss-of-heterozygosity events. Together, these results help distinguish aspects of evolutionary dynamics that are likely to be general features of adaptation across many systems from those that are specific to individual organisms and environmental conditions.

Factors enforcing the species boundary between the human pathogens Cryptococcus neoformans and Cryptococcus deneoformans

Hybridization has resulted in the origin and variation in extant species, and hybrids continue to arise despite pre- and post-zygotic barriers that limit their formation and evolutionary success. One important system that maintains species boundaries in prokaryotes and eukaryotes is the mismatch repair pathway, which blocks recombination between divergent DNA sequences. Previous studies illuminated the role of the mismatch repair component Msh2 in blocking genetic recombination between divergent DNA during meiosis. Loss of Msh2 results in increased interspecific genetic recombination in bacterial and yeast models, and increased viability of progeny derived from yeast hybrid crosses. Hybrid isolates of two pathogenic fungal Cryptococcus species, Cryptococcus neoformans and Cryptococcus deneoformans, are isolated regularly from both clinical and environmental sources. In the present study, we sought to determine if loss of Msh2 would relax the species boundary between C. neoformans and C. deneoformans. We found that crosses between these two species in which both parents lack Msh2 produced hybrid progeny with increased viability and high levels of aneuploidy. Whole-genome sequencing revealed few instances of recombination among hybrid progeny and did not identify increased levels of recombination in progeny derived from parents lacking Msh2. Several hybrid progeny produced structures associated with sexual reproduction when incubated alone on nutrient-rich medium in light, a novel phenotype in Cryptococcus. These findings represent a unique, unexpected case where rendering the mismatch repair system defective did not result in increased meiotic recombination across a species boundary. This suggests that alternative pathways or other mismatch repair components limit meiotic recombination between homeologous DNA and enforce species boundaries in the basidiomycete Cryptococcus species.

Several mechanisms enforce species boundaries by either preventing the formation of hybrid zygotes, known as pre-zygotic barriers, or preventing the viability and fecundity of hybrids, known as post-zygotic barriers. Despite these barriers, interspecific hybrids form at an appreciable frequency, such as hybrid isolates of the human fungal pathogenic species, Cryptococcus neoformans and Cryptococcus deneoformans, which are regularly isolated from both clinical and environmental sources. C. neoformans x C. deneoformans hybrids are typically highly aneuploid, sterile, and display phenotypes intermediate to those of either parent, although self-fertile isolates and transgressive phenotypes have been observed. One important mechanism known to enforce species boundaries or lead to incipient speciation is the DNA mismatch repair system, which blocks recombination between divergent DNA sequences during meiosis. The aim of this study was to determine if genetically deleting the DNA mismatch repair component Msh2 would relax the species boundary between C. neoformans and C. deneoformans. Progeny derived from C. neoformans x C. deneoformans crosses in which both parental strains lacked Msh2 had higher viability, and unlike previous studies in Saccharomyces, these Cryptococcus hybrid progeny had higher levels of aneuploidy and no observable increase in meiotic recombination at the whole-genome level.

Characterization of systemic genomic instability in budding yeast

Conventional models of genome evolution are centered around the principle that mutations form independently of each other and build up slowly over time. We characterized the occurrence of bursts of genome-wide loss-of-heterozygosity (LOH) in Saccharomyces cerevisiae, providing support for an additional nonindependent and faster mode of mutation accumulation. We initially characterized a yeast clone isolated for carrying an LOH event at a specific chromosome site, and surprisingly found that it also carried multiple unselected rearrangements elsewhere in its genome. Whole-genome analysis of over 100 additional clones selected for carrying primary LOH tracts revealed that they too contained unselected structural alterations more often than control clones obtained without any selection. We also measured the rates of coincident LOH at two different chromosomes and found that double LOH formed at rates 14- to 150-fold higher than expected if the two underlying single LOH events occurred independently of each other. These results were consistent across different strain backgrounds and in mutants incapable of entering meiosis. Our results indicate that a subset of mitotic cells within a population can experience discrete episodes of systemic genomic instability, when the entire genome becomes vulnerable and multiple chromosomal alterations can form over a narrow time window. They are reminiscent of early reports from the classic yeast genetics literature, as well as recent studies in humans, both in cancer and genomic disorder contexts. The experimental model we describe provides a system to further dissect the fundamental biological processes responsible for punctuated bursts of structural genomic variation.

Punctuated Aneuploidization of the Budding Yeast Genome

Remarkably complex patterns of aneuploidy have been observed in the genomes of many eukaryotic cell types, ranging from brewing yeasts to tumor cells. Such aberrant karyotypes are generally thought to take shape progressively over many generations, but evidence also suggests that genomes may undergo faster modes of evolution. Here, we used diploid Saccharomyces cerevisiae cells to investigate the dynamics with which aneuploidies arise. We found that cells selected for the loss of a single chromosome often acquired additional unselected aneuploidies concomitantly. The degrees to which these genomes were altered fell along a spectrum, ranging from simple events affecting just a single chromosome, to systemic events involving many. The striking complexity of karyotypes arising from systemic events, combined with the high frequency at which we detected them, demonstrates that cells can rapidly achieve highly altered genomic configurations during temporally restricted episodes of genomic instability.

Identifying Small Molecules That Promote Quasipalindrome-Associated Template-Switch Mutations in Escherichia coli

DNA can assemble into non-B form structures that stall replication and cause genomic instability. One such secondary structure results from an inverted DNA repeat that can assemble into hairpin and cruciform structures during DNA replication. Quasipalindromes (QP), imperfect inverted repeats, are sites of mutational hotspots. Quasipalindrome-associated mutations (QPMs) occur through a template-switch mechanism in which the replicative polymerase stalls at a QP site and uses the nascent strand as a template instead of the correct template strand. This mutational event causes the QP to become a perfect or more perfect inverted repeat. Since it is not fully understood how template-switch events are stimulated or repressed, we designed a high-throughput screen to discover drugs that affect these events. QP reporters were engineered in the Escherichia coli lacZ gene to allow us to study template-switch events specifically. We tested 700 compounds from the NIH Clinical Collection through a disk diffusion assay and identified 11 positive hits. One of the hits was azidothymidine (zidovudine, AZT), a thymidine analog and DNA chain terminator. The other ten were found to be fluoroquinolone antibiotics, which induce DNA-protein crosslinks. This work shows that our screen is useful in identifying small molecules that affect quasipalindrome-associated template-switch mutations. We are currently assessing more small molecule libraries and applying this method to study other types of mutations.

Role of error-prone DNA polymerases in spontaneous mutagenesis in Caulobacter crescentus

Spontaneous mutations are important players in evolution. Nevertheless, there is a paucity of information about the mutagenic processes operating in most bacterial species. In this work, we implemented two forward mutational markers for studies in Caulobacter crescentus. We confirmed previous results in which A:T → G:C transitions are the most prevalent type of spontaneous base substitutions in this organism, although there is considerable deviation from this trend in one of the loci analyzed. We also investigated the role of dinB and imuC, encoding error-prone DNA polymerases, in spontaneous mutagenesis in this GC-rich organism. Both dinB and imuC mutant strains show comparable mutation rates to the parental strain. Nevertheless, both strains show differences in the base substitution patterns, and the dinB mutant strain shows a striking reduction in the number of spontaneous -1 deletions and an increase in C:G → T:A transitions in both assays.

Ras-Mediated Activation of NF-κB and DNA Damage Response in Carcinogenesis

Cancer is a multi-step process during which cells acquire mutations that eventually lead to uncontrolled cell growth and division and evasion of programmed cell death. The oncogenes such as Ras and c-Myc may be responsible in all three major stages of cancer i.e., early, intermediate, and late. The NF-κB has been shown to control the expression of genes linked with tumor pathways such as chronic inflammation, tumor cell survival, anti-apoptosis, proliferation, invasion, and angiogenesis. In the last few decades, various biomarker pathways have been identified that play a critical role in carcinogenesis such as Ras, NF-κB and DNA damage.

Contribution of GO System Glycosylases to Mutation Prevention in Caulobacter crescentus

8-oxo-7,8-dihydroguanine, commonly referred to as 8-oxoG, is considered one of the most predominant oxidative lesions formed in DNA. Due to its ability to pair with adenines in its syn configuration, this lesion has a strong mutagenic potential in both eukaryotes and prokaryotes. Escherichia coli cells are endowed with the GO system, which protects them from the mutagenic properties of this lesion when formed both in cellular DNA and the nucleotide pool. MutY and MutM (Fpg) DNA glycosylases are crucial components of the GO system. A strong mutator phenotype of the Escherichia coli mutM mutY double mutant underscores the importance of 8-oxoG repair for genomic stability. Here, we report that in Caulobacter crescentus, a widely studied alpha-proteobacterium with a GC-rich genome, the combined lack of MutM and MutY glycosylases produces a more modest mutator phenotype when compared to E. coli. Genetic analysis indicates that other glycosylases and other repair pathways do not act synergistically with the GO system for spontaneous mutation prevention. We also show that there is not a statistically significant difference in the spontaneous levels 8-oxodGuo in E. coli and C. crescentus, suggesting that other yet to be identified differences in repair or replication probably account for the differential importance of the GO system between these two species. Environ. Mol. Mutagen. 61:246-255, 2020. © 2019 Wiley Periodicals, Inc.

Large-scale contractions of Friedreich's ataxia GAA repeats in yeast occur during DNA replication due to their triplex-forming ability

Expansions of GAA repeats cause a severe hereditary neurodegenerative disease, Friedreich’s ataxia. In this study, we characterized the mechanisms of GAA repeat contractions in a yeast experimental system. These mechanisms might, in the long run, aid development of a therapy for this currently incurable disease. We show that GAA repeats contract during DNA replication, which can explain the high level of somatic instability of this repeat in patient tissues. We also provided evidence that a triple-stranded DNA structure is at the heart of GAA repeat instability. This discovery highlights the role of triplex DNA in genome instability and human disease.

Friedreich’s ataxia (FRDA) is a human hereditary disease caused by the presence of expanded (GAA)_n repeats in the first intron of the FXN gene [V. Campuzano et al., Science 271, 1423–1427 (1996)]. In somatic tissues of FRDA patients, (GAA)_n repeat tracts are highly unstable, with contractions more common than expansions [R. Sharma et al., Hum. Mol. Genet. 11, 2175–2187 (2002)]. Here we describe an experimental system to characterize GAA repeat contractions in yeast and to conduct a genetic analysis of this process. We found that large-scale contraction is a one-step process, resulting in a median loss of ∼60 triplet repeats. Our genetic analysis revealed that contractions occur during DNA replication, rather than by various DNA repair pathways. Repeats contract in the course of lagging-strand synthesis: The processivity subunit of DNA polymerase δ, Pol32, and the catalytic domain of Rev1, a translesion polymerase, act together in the same pathway to counteract contractions. Accumulation of single-stranded DNA (ssDNA) in the lagging-strand template greatly increases the probability that (GAA)_n repeats contract, which in turn promotes repeat instability in rfa1, rad27, and dna2 mutants. Finally, by comparing contraction rates for homopurine-homopyrimidine repeats differing in their mirror symmetry, we found that contractions depend on a repeat’s triplex-forming ability. We propose that accumulation of ssDNA in the lagging-strand template fosters the formation of a triplex between the nascent and fold-back template strands of the repeat. Occasional jumps of DNA polymerase through this triplex hurdle, result in repeat contractions in the nascent lagging strand.

Genetic Assays to Study Repeat Fragility in Saccharomyces cerevisiae

Trinucleotide repeats are common in the human genome and can undergo changes in repeat number and cause length-dependent chromosome fragility. Expanded CAG repeats have been linked to over 14 human diseases and are considered hotspots for breakage and genomic rearrangement. Here we describe two Saccharomyces cerevisiae based assays that evaluate the rate of chromosome breakage that occurs within a repeat tract (fragility), with variations that allow the role of transcription to be evaluated. The first fragility assay utilizes end-loss and subsequent telomere addition as the main mode of repair of a yeast artificial chromosome (YAC). The second fragility assay relies on the fact that a chromosomal break stimulates recombination-mediated repair. A PCR-based assay can be used to evaluate instability of the repeat in the same conditions used to measure repeat fragility. These assays have contributed to understanding the genetic mechanisms that cause chromosome breaks and tract-length changes at unstable trinucleotide repeats.

RNA-DNA hybrids promote the expansion of Friedreich's ataxia (GAA)n repeats via break-induced replication

Expansion of simple DNA repeats is responsible for numerous hereditary diseases in humans. The role of DNA replication, repair and transcription in the expansion process has been well documented. Here we analyzed, in a yeast experimental system, the role of RNA–DNA hybrids in genetic instability of long (GAA)_n repeats, which cause Friedreich’s ataxia. Knocking out both yeast RNase H enzymes, which counteract the formation of RNA–DNA hybrids, increased (GAA)_n repeat expansion and contraction rates when the repetitive sequence was transcribed. Unexpectedly, we observed a similar increase in repeat instability in RNase H-deficient cells when we either changed the direction of transcription-replication collisions, or flipped the repeat sequence such that the (UUC)_n run occurred in the transcript. The increase in repeat expansions in RNase H-deficient strains was dependent on Rad52 and Pol32 proteins, suggesting that break-induced replication (BIR) is responsible for this effect. We conclude that expansions of (GAA)_n repeats are induced by the formation of RNA–DNA hybrids that trigger BIR. Since this stimulation is independent of which strand of the repeat (homopurine or homopyrimidine) is in the RNA transcript, we hypothesize that triplex H-DNA structures stabilized by an RNA–DNA hybrid (H-loops), rather than conventional R-loops, could be responsible.

Guidelines for DNA recombination and repair studies: Cellular assays of DNA repair pathways

Understanding the plasticity of genomes has been greatly aided by assays for recombination, repair and mutagenesis. These assays have been developed in microbial systems that provide the advantages of genetic and molecular reporters that can readily be manipulated. Cellular assays comprise genetic, molecular, and cytological reporters. The assays are powerful tools but each comes with its particular advantages and limitations. Here the most commonly used assays are reviewed, discussed, and presented as the guidelines for future studies.