Mutation hot spots in yeast caused by long-range clustering of homopolymeric sequences

MutLα suppresses error-prone DNA mismatch repair and preferentially protects noncoding DNA from mutations

The DNA mismatch repair (MMR) system promotes genome stability and protects humans from certain types of cancer. Its primary function is the correction of DNA polymerase errors. MutLα is an important eukaryotic MMR factor. We have examined the contributions of MutLα to maintaining genome stability. We show here that loss of MutLα in yeast increases the genome-wide mutation rate by ∼130-fold and generates a genome-wide mutation spectrum that consists of small indels and base substitutions. We also show that loss of yeast MutLα leads to error-prone MMR that produces T > C base substitutions in 5'-ATA-3' sequences. In agreement with this finding, our examination of human whole-genome DNA sequencing data has revealed that loss of MutLα in induced pluripotent stem cells triggers error-prone MMR that leads to the formation of T > C mutations in 5'-NTN-3' sequences. Our further analysis has shown that MutLα-independent MMR plays a role in suppressing base substitutions in N3 homopolymeric runs. In addition, we describe that MutLα preferentially protects noncoding DNA from mutations. Our study defines the contributions of MutLα-dependent and independent mechanisms to genome-wide MMR.

Cellular Responses to Widespread DNA Replication Stress

Replicative DNA polymerases are blocked by nearly all types of DNA damage. The resulting DNA replication stress threatens genome stability. DNA replication stress is also caused by depletion of nucleotide pools, DNA polymerase inhibitors, and DNA sequences or structures that are difficult to replicate. Replication stress triggers complex cellular responses that include cell cycle arrest, replication fork collapse to one-ended DNA double-strand breaks, induction of DNA repair, and programmed cell death after excessive damage. Replication stress caused by specific structures (e.g., G-rich sequences that form G-quadruplexes) is localized but occurs during the S phase of every cell division. This review focuses on cellular responses to widespread stress such as that caused by random DNA damage, DNA polymerase inhibition/nucleotide pool depletion, and R-loops. Another form of global replication stress is seen in cancer cells and is termed oncogenic stress, reflecting dysregulated replication origin firing and/or replication fork progression. Replication stress responses are often dysregulated in cancer cells, and this too contributes to ongoing genome instability that can drive cancer progression. Nucleases play critical roles in replication stress responses, including MUS81, EEPD1, Metnase, CtIP, MRE11, EXO1, DNA2-BLM, SLX1-SLX4, XPF-ERCC1-SLX4, Artemis, XPG, FEN1, and TATDN2. Several of these nucleases cleave branched DNA structures at stressed replication forks to promote repair and restart of these forks. We recently defined roles for EEPD1 in restarting stressed replication forks after oxidative DNA damage, and for TATDN2 in mitigating replication stress caused by R-loop accumulation in BRCA1-defective cells. We also discuss how insights into biological responses to genome-wide replication stress can inform novel cancer treatment strategies that exploit synthetic lethal relationships among replication stress response factors.

Genome-wide contributions of the MutSα- and MutSβ-dependent DNA mismatch repair pathways to the maintenance of genetic stability in S. cerevisiae

The DNA mismatch repair (MMR) system is a major DNA repair system that suppresses inherited and sporadic cancers in humans. In eukaryotes the MutSα-dependent and MutSβ-dependent MMR pathways correct DNA polymerase errors. Here, we investigated these two pathways on a whole-genome level in S. cerevisiae. We found that inactivation of MutSα-dependent MMR by deletion of the MSH6 gene increases the genome-wide mutation rate by ∼17-fold, and loss of MutSβ-dependent MMR via deletion of MSH3 elevates the genome-wide mutation rate by ∼4-fold. We also found that MutSα-dependent MMR does not show a preference for protecting coding or noncoding DNA from mutations, whereas MutSβ-dependent MMR preferentially protects noncoding DNA from mutations. The most frequent mutations in the msh6Δ strain are C>T transitions, whereas 1-6-bp deletions are the most common genetic alterations in the msh3Δ strain. Strikingly, MutSα-dependent MMR is more important than MutSβ-dependent MMR for protection from 1-bp insertions, while MutSβ-dependent MMR has a more critical role in the defense against 1-bp deletions and 2-6-bp indels. We also determined that a mutational signature of yeast MSH6 loss is similar to mutational signatures of human MMR deficiency. Furthermore, our analysis showed that compared to other 5'-NCN-3' trinucleotides, 5'-GCA-3' trinucleotides are at the highest risk of accumulating C>T transitions at the central position in the msh6Δ cells and that the presence of a G/A base at the -1 position is important for the efficient MutSα-dependent suppression of C>T transitions. Our results highlight key differences between the roles of the MutSα-dependent and MutSβ-dependent MMR pathways.

Heterogeneous Mutation Rates and Spectra in Yeast Hybrids

Mutation rates and spectra vary between species and among populations. Hybridization can contribute to this variation, but its role remains poorly understood. Estimating mutation rates requires controlled conditions where the effect of natural selection can be minimized. One way to achieve this is through mutation accumulation experiments coupled with genome sequencing. Here, we investigate 400 mutation accumulation lines initiated from 11 genotypes spanning intralineage, interlineage, and interspecific crosses of the yeasts Saccharomyces paradoxus and S. cerevisiae and propagated for 770 generations. We find significant differences in mutation rates and spectra among crosses, which are not related to the level of divergence of parental strains but are specific to some genotype combinations. Differences in number of generations and departures from neutrality play a minor role, whereas polyploidy and loss of heterozygosity impact mutation rates in some of the hybrid crosses in an opposite way.

Stability across the Whole Nuclear Genome in the Presence and Absence of DNA Mismatch Repair

We describe the contribution of DNA mismatch repair (MMR) to the stability of the eukaryotic nuclear genome as determined by whole-genome sequencing. To date, wild-type nuclear genome mutation rates are known for over 40 eukaryotic species, while measurements in mismatch repair-defective organisms are fewer in number and are concentrated on Saccharomyces cerevisiae and human tumors. Well-studied organisms include Drosophila melanogaster and Mus musculus, while less genetically tractable species include great apes and long-lived trees. A variety of techniques have been developed to gather mutation rates, either per generation or per cell division. Generational rates are described through whole-organism mutation accumulation experiments and through offspring–parent sequencing, or they have been identified by descent. Rates per somatic cell division have been estimated from cell line mutation accumulation experiments, from systemic variant allele frequencies, and from widely spaced samples with known cell divisions per unit of tissue growth. The latter methods are also used to estimate generational mutation rates for large organisms that lack dedicated germlines, such as trees and hyphal fungi. Mechanistic studies involving genetic manipulation of MMR genes prior to mutation rate determination are thus far confined to yeast, Arabidopsis thaliana, Caenorhabditis elegans, and one chicken cell line. A great deal of work in wild-type organisms has begun to establish a sound baseline, but far more work is needed to uncover the variety of MMR across eukaryotes. Nonetheless, the few MMR studies reported to date indicate that MMR contributes 100-fold or more to genome stability, and they have uncovered insights that would have been impossible to obtain using reporter gene assays.

Pervasive function and evidence for selection across standing genetic variation in S. cerevisiae

Quantitative genetics aims to map genotype to phenotype, often with the goal of understanding how organisms evolved. However, it remains unclear whether the genetic variants identified are exemplary of evolution. Here we analyzed progeny of two wild Saccharomyces cerevisiae isolates to identify 195 loci underlying complex metabolic traits, resolving 107 to single polymorphisms with diverse molecular mechanisms. More than 20% of causal variants exhibited patterns of emergence inconsistent with neutrality. Moreover, contrary to drift-centric expectation, variation in diverse wild yeast isolates broadly exhibited this property: over 30% of shared natural variants exhibited phylogenetic signatures suggesting that they are not neutral. This pattern is likely attributable to both homoplasy and balancing selection on ancestral polymorphism. Variants that emerged repeatedly were more likely to have done so in isolates from the same ecological niche. Our results underscore the power of super-resolution mapping of ecologically relevant traits in understanding adaptation and evolution.

Muver, a computational framework for accurately calling accumulated mutations

BACKGROUND

Identification of mutations from next-generation sequencing data typically requires a balance between sensitivity and accuracy. This is particularly true of DNA insertions and deletions (indels), that can impart significant phenotypic consequences on cells but are harder to call than substitution mutations from whole genome mutation accumulation experiments. To overcome these difficulties, we present muver, a computational framework that integrates established bioinformatics tools with novel analytical methods to generate mutation calls with the extremely low false positive rates and high sensitivity required for accurate mutation rate determination and comparison.

RESULTS

Muver uses statistical comparison of ancestral and descendant allelic frequencies to identify variant loci and assigns genotypes with models that include per-sample assessments of sequencing errors by mutation type and repeat context. Muver identifies maximally parsimonious mutation pathways that connect these genotypes, differentiating potential allelic conversion events and delineating ambiguities in mutation location, type, and size. Benchmarking with a human gold standard father-son pair demonstrates muver's sensitivity and low false positive rates. In DNA mismatch repair (MMR) deficient Saccharomyces cerevisiae, muver detects multi-base deletions in homopolymers longer than the replicative polymerase footprint at rates greater than predicted for sequential single-base deletions, implying a novel multi-repeat-unit slippage mechanism.

CONCLUSIONS

Benchmarking results demonstrate the high accuracy and sensitivity achieved with muver, particularly for indels, relative to available tools. Applied to an MMR-deficient Saccharomyces cerevisiae system, muver mutation calls facilitate mechanistic insights into DNA replication fidelity.

Phenotype variability and allelic heterogeneity in KMT2B-Associated disease

BACKGROUND

Mutations in Lysine-Specific Histone Methyltransferase 2B gene (KMT2B) have been reported to be associated with complex early-onset dystonia. Almost all reported KMT2B mutations occurred de novo in the paternal germline or in the early development of the patient. We describe clinico-genetic features on four Japanese patients with novel de novo mutations and demonstrate the phenotypic spectrum of KMT2B mutations.

METHODS

We performed genetic studies, including trio-based whole exome sequencing (WES), in a cohort of Japanese patients with a seemingly sporadic early-onset generalized combined dystonia. Potential effects by the identified nucleotide variations were evaluated biologically. Genotype-phenotype correlations were also investigated.

RESULTS

Four patients had de novo heterozygous mutations in KMT2B, c.309delG, c.1656dupC, c.3325_3326insC, and c.5636delG. Biological analysis of KMT2B mRNA levels showed a reduced expression of mutant transcript frame. All patients presented with motor milestone delay, microcephaly, mild psychomotor impairment, childhood-onset generalized dystonia and superimposed choreoathetosis or myoclonus. One patient cannot stand due to axial hypotonia associated with cerebellar dysfunction. Three patients had bilateral globus pallidal deep brain stimulation (DBS) with excellent or partial response.

CONCLUSIONS

We further demonstrate the allelic heterogeneity and phenotypic variations of KMT2B-associated disease. Haploinsufficiency is one of molecular pathomechanisms underlying the disease. Cardinal clinical features include combined dystonia accompanying mild psychomotor disability. Cerebellum would be affected in KMT2B-associated disease.

Comparative Genomics of Two Sequential Candida glabrata Clinical Isolates

Candida glabrata is an important fungal pathogen which develops rapid antifungal resistance in treated patients. It is known that azole treatments lead to antifungal resistance in this fungal species and that multidrug efflux transporters are involved in this process. Specific mutations in the transcriptional regulator PDR1 result in upregulation of the transporters. In addition, we showed that the PDR1 mutations can contribute to enhance virulence in animal models. In this study, we were interested to compare genomes of two specific C. glabrata-related isolates, one of which was azole susceptible (DSY562) while the other was azole resistant (DSY565). DSY565 contained a PDR1 mutation (L280F) and was isolated after a time-lapse of 50 d of azole therapy. We expected that genome comparisons between both isolates could reveal additional mutations reflecting host adaptation or even additional resistance mechanisms. The PacBio technology used here yielded 14 major contigs (sizes 0.18-1.6 Mb) and mitochondrial genomes from both DSY562 and DSY565 isolates that were highly similar to each other. Comparisons of the clinical genomes with the published CBS138 genome indicated important genome rearrangements, but not between the clinical strains. Among the unique features, several retrotransposons were identified in the genomes of the investigated clinical isolates. DSY562 and DSY565 each contained a large set of adhesin-like genes (101 and 107, respectively), which exceed by far the number of reported adhesins (63) in the CBS138 genome. Comparison between DSY562 and DSY565 yielded 17 nonsynonymous SNPs (among which the was the expected PDR1 mutation) as well as small size indels in coding regions (11) but mainly in adhesin-like genes. The genomes contained a DNA mismatch repair allele of MSH2 known to be involved in the so-called hyper-mutator phenotype of this yeast species and the number of accumulated mutations between both clinical isolates is consistent with the presence of a MSH2 defect. In conclusion, this study is the first to compare genomes of C. glabrata sequential clinical isolates using the PacBio technology as an approach. The genomes of these isolates taken in the same patient at two different time points exhibited limited variations, even if submitted to the host pressure.

Extremely Rare Polymorphisms in Saccharomyces cerevisiae Allow Inference of the Mutational Spectrum

The characterization of mutational spectra is usually carried out in one of three ways-by direct observation through mutation accumulation (MA) experiments, through parent-offspring sequencing, or by indirect inference from sequence data. Direct observations of spontaneous mutations with MA experiments are limited, given (i) the rarity of spontaneous mutations, (ii) applicability only to laboratory model species with short generation times, and (iii) the possibility that mutational spectra under lab conditions might be different from those observed in nature. Trio sequencing is an elegant solution, but it is not applicable in all organisms. Indirect inference, usually from divergence data, faces no such technical limitations, but rely upon critical assumptions regarding the strength of natural selection that are likely to be violated. Ideally, new mutational events would be directly observed before the biased filter of selection, and without the technical limitations common to lab experiments. One approach is to identify very young mutations from population sequencing data. Here we do so by leveraging two characteristics common to all new mutations-new mutations are necessarily rare in the population, and absent in the genomes of immediate relatives. From 132 clinical yeast strains, we were able to identify 1,425 putatively new mutations and show that they exhibit extremely low signatures of selection, as well as display a mutational spectrum that is similar to that identified by a large scale MA experiment. We verify that population sequencing data are a potential wealth of information for inferring mutational spectra, and should be considered for analysis where MA experiments are infeasible or especially tedious.

Sister chromatid telomere fusions, but not NHEJ-mediated inter-chromosomal telomere fusions, occur independently of DNA ligases 3 and 4

Telomeres shorten with each cell division and can ultimately become substrates for nonhomologous end-joining repair, leading to large-scale genomic rearrangements of the kind frequently observed in human cancers. We have characterized more than 1400 telomere fusion events at the single-molecule level, using a combination of high-throughput sequence analysis together with experimentally induced telomeric double-stranded DNA breaks. We show that a single chromosomal dysfunctional telomere can fuse with diverse nontelomeric genomic loci, even in the presence of an otherwise stable genome, and that fusion predominates in coding regions. Fusion frequency was markedly increased in the absence of TP53 checkpoint control and significantly modulated by the cellular capacity for classical, versus alternative, nonhomologous end joining (NHEJ). We observed a striking reduction in inter-chromosomal fusion events in cells lacking DNA ligase 4, in contrast to a remarkably consistent profile of intra-chromosomal fusion in the context of multiple genetic knockouts, including DNA ligase 3 and 4 double-knockouts. We reveal distinct mutational signatures associated with classical NHEJ-mediated inter-chromosomal, as opposed to alternative NHEJ-mediated intra-chromosomal, telomere fusions and evidence for an unanticipated sufficiency of DNA ligase 1 for these intra-chromosomal events. Our findings have implications for mechanisms driving cancer genome evolution.

The Effect of Local Sequence Context on Mutational Bias of Genes Encoded on the Leading and Lagging Strands

All organisms must replicate their genetic information accurately to ensure its faithful transmission. DNA polymerase errors provide an important source of genetic variation that can drive evolution. Understanding the origins of genetic variation will inform our understanding of evolution and the development of genetic diseases. A number of factors have been proposed to influence mutagenesis [1-10]. Here, we used mutation accumulation lines, whole-genome sequencing, and whole-transcriptome analysis to study the locations and rate at which mutations arise in bacteria with as little selection bias as possible [11, 12]. Our analysis of greater than 7,000 replication errors in over 180 sequenced lines that underwent a total of more than 370,000 generations has provided new insights into how DNA polymerase errors sculpt genetic variation and drive evolution. Homopolymer run enrichment outside of genes causes insertions and deletions in these regions. Genes encoded in the lagging strand are transcribed such that RNA polymerase and DNA polymerase collide head-on. Head-on genes have been proposed to mutate at a higher rate than genes transcribed codirectionally with DNA polymerase progression due to conflicts between transcription and DNA replication [6, 10]. We did not detect associations between the number of base pair substitutions in genes and their orientation or expression. Strikingly, any higher mutation rate for head-on genes can be explained by differing sequence composition between the leading and lagging strands and the error bias for DNA polymerase in specific sequence contexts. Therefore, we find local sequence context is the major determinant of mutagenesis in bacteria.

The Eukaryotic Mismatch Recognition Complexes Track with the Replisome during DNA Synthesis

During replication, mismatch repair proteins recognize and repair mispaired bases that escape the proofreading activity of DNA polymerase. In this work, we tested the model that the eukaryotic mismatch recognition complex tracks with the advancing replisome. Using yeast, we examined the dynamics during replication of the leading strand polymerase Polε using Pol2 and the eukaryotic mismatch recognition complex using Msh2, the invariant protein involved in mismatch recognition. Specifically, we synchronized cells and processed samples using chromatin immunoprecipitation combined with custom DNA tiling arrays (ChIP-chip). The Polε signal was not detectable in G1, but was observed at active origins and replicating DNA throughout S-phase. The Polε signal provided the resolution to track origin firing timing and efficiencies as well as replisome progression rates. By detecting Polε and Msh2 dynamics within the same strain, we established that the mismatch recognition complex binds origins and spreads to adjacent regions with the replisome. In mismatch repair defective PCNA mutants, we observed that Msh2 binds to regions of replicating DNA, but the distribution and dynamics are altered, suggesting that PCNA is not the sole determinant for the mismatch recognition complex association with replicating regions, but may influence the dynamics of movement. Using biochemical and genomic methods, we provide evidence that both MutS complexes are in the vicinity of the replisome to efficiently repair the entire spectrum of mutations during replication. Our data supports the model that the proximity of MutSα/β to the replisome for the efficient repair of the newly synthesized strand before chromatin reassembles.

During replication, errors that escape the replication machinery are identified and repaired by DNA mismatch repair proteins. A mismatch in the helix is recognized by MutS homologs and subsequent events include excision of the error-containing strand followed by re-synthesis. A critical step in this process is directing repair to the newly synthesized strand. Current data suggest that transient discontinuities in the DNA backbone, known as nicks, generated during replication serve as the strand discrimination signals. Additionally, proteins that package DNA have the capacity to block mismatch recognition and are known to rapidly assemble behind the replication fork. Thus, there must be a short window of opportunity for the mismatch recognition complexes to scan for mismatches and access the strand discrimination signals. To address these issues, we tested the model that the mismatch recognition complexes track with the replisome. We employed high resolution genomic methods to determine that during replication, the mismatch recognition complexes bind origins of replication and advances with the replisome. The findings support the hypothesis that the mismatch recognition proteins track with the DNA replication machinery to accurately survey and repair the newly synthesized strands while the DNA is unpackaged and strand specificity signals are accessible.

A Genetic Incompatibility Accelerates Adaptation in Yeast

During mismatch repair (MMR) MSH proteins bind to mismatches that form as the result of DNA replication errors and recruit MLH factors such as Mlh1-Pms1 to initiate excision and repair steps. Previously, we identified a negative epistatic interaction involving naturally occurring polymorphisms in the MLH1 and PMS1 genes of baker’s yeast. Here we hypothesize that a mutagenic state resulting from this negative epistatic interaction increases the likelihood of obtaining beneficial mutations that can promote adaptation to stress conditions. We tested this by stressing yeast strains bearing mutagenic (incompatible) and non-mutagenic (compatible) mismatch repair genotypes. Our data show that incompatible populations adapted more rapidly and without an apparent fitness cost to high salt stress. The fitness advantage of incompatible populations was rapid but disappeared over time. The fitness gains in both compatible and incompatible strains were due primarily to mutations in PMR1 that appeared earlier in incompatible evolving populations. These data demonstrate a rapid and reversible role (by mating) for genetic incompatibilities in accelerating adaptation in eukaryotes. They also provide an approach to link experimental studies to observational population genomics.

In nature, bacterial populations with high mutation rates can adapt faster to new environments by acquiring beneficial mutations. However, such populations also accumulate harmful mutations that reduce their fitness. We show that the model eukaryote baker’s yeast can use a similar mutator strategy to adapt to new environments. The mutator state that we observed resulted from an incompatibility involving two genes, MLH1 and PMS1, that work together to remove DNA replication errors through a spellchecking mismatch repair mechanism. This incompatibility can occur through mating between baker’s yeast from different genetic backgrounds, yielding mutator offspring containing an MLH1-PMS1 combination not present in either parent. Interestingly, these offspring adapted more rapidly to stress, compared to the parental strains, and did so without an overall loss in fitness. DNA sequencing analyses of baker’s yeast strains from across the globe support the presence of incompatible hybrid yeast strains in nature. These observations provide a powerful model to understand how the segregation of defects in DNA mismatch repair can serve as an effective strategy to enable eukaryotes to adapt to changing environments.

Dissecting genetic and environmental mutation signatures with model organisms

Deep sequencing has impacted on cancer research by enabling routine sequencing of genomes and exomes to identify genetic changes associated with carcinogenesis. Researchers can now use the frequency, type, and context of all mutations in tumor genomes to extract mutation signatures that reflect the driving mutational processes. Identifying mutation signatures, however, may not immediately suggest a mechanism. Consequently, several recent studies have employed deep sequencing of model organisms exposed to discrete genetic or environmental perturbations. These studies exploit the simpler genomes and availability of powerful genetic tools in model organisms to analyze mutation signatures under controlled conditions, forging mechanistic links between mutational processes and signatures. We discuss the power of this approach and suggest that many such studies may be on the horizon.

Differences in genome-wide repeat sequence instability conferred by proofreading and mismatch repair defects

Mutation rates are used to calibrate molecular clocks and to link genetic variants with human disease. However, mutation rates are not uniform across each eukaryotic genome. Rates for insertion/deletion (indel) mutations have been found to vary widely when examined in vitro and at specific loci in vivo. Here, we report the genome-wide rates of formation and repair of indels made during replication of yeast nuclear DNA. Using over 6000 indels accumulated in four mismatch repair (MMR) defective strains, and statistical corrections for false negatives, we find that indel rates increase by 100 000-fold with increasing homonucleotide run length, representing the greatest effect on replication fidelity of any known genomic parameter. Nonetheless, long genomic homopolymer runs are overrepresented relative to random chance, implying positive selection. Proofreading defects in the replicative polymerases selectively increase indel rates in short repetitive tracts, likely reflecting the distance over which Pols δ and ϵ interact with duplex DNA upstream of the polymerase active site. In contrast, MMR defects hugely increase indel mutagenesis in long repetitive sequences. Because repetitive sequences are not uniformly distributed among genomic functional elements, the quantitatively different consequences on genome-wide repeat sequence instability conferred by defects in proofreading and MMR have important biological implications.

Eukaryotic Mismatch Repair in Relation to DNA Replication

Three processes act in series to accurately replicate the eukaryotic nuclear genome. The major replicative DNA polymerases strongly prevent mismatch formation, occasional mismatches that do form are proofread during replication, and rare mismatches that escape proofreading are corrected by mismatch repair (MMR). This review focuses on MMR in light of increasing knowledge about nuclear DNA replication enzymology and the rate and specificity with which mismatches are generated during leading- and lagging-strand replication. We consider differences in MMR efficiency in relation to mismatch recognition, signaling to direct MMR to the nascent strand, mismatch removal, and the timing of MMR. These studies are refining our understanding of relationships between generating and repairing replication errors to achieve accurate replication of both DNA strands of the nuclear genome.

Heterogeneous polymerase fidelity and mismatch repair bias genome variation and composition

Mutational heterogeneity must be taken into account when reconstructing evolutionary histories, calibrating molecular clocks, and predicting links between genes and disease. Selective pressures and various DNA transactions have been invoked to explain the heterogeneous distribution of genetic variation between species, within populations, and in tissue-specific tumors. To examine relationships between such heterogeneity and variations in leading- and lagging-strand replication fidelity and mismatch repair, we accumulated 40,000 spontaneous mutations in eight diploid yeast strains in the absence of selective pressure. We found that replicase error rates vary by fork direction, coding state, nucleosome proximity, and sequence context. Further, error rates and DNA mismatch repair efficiency both vary by mismatch type, responsible polymerase, replication time, and replication origin proximity. Mutation patterns implicate replication infidelity as one driver of variation in somatic and germline evolution, suggest mechanisms of mutual modulation of genome stability and composition, and predict future observations in specific cancers.

Genome destabilizing mutator alleles drive specific mutational trajectories in Saccharomyces cerevisiae

In addition to environmental factors and intrinsic variations in base substitution rates, specific genome-destabilizing mutations can shape the mutational trajectory of genomes. How specific alleles influence the nature and position of accumulated mutations in a genomic context is largely unknown. Understanding the impact of genome-destabilizing alleles is particularly relevant to cancer genomes where biased mutational signatures are identifiable. We first created a more complete picture of cellular pathways that impact mutation rate using a primary screen to identify essential Saccharomyces cerevisiae gene mutations that cause mutator phenotypes. Drawing primarily on new alleles identified in this resource, we measure the impact of diverse mutator alleles on mutation patterns directly by whole-genome sequencing of 68 mutation-accumulation strains derived from wild-type and 11 parental mutator genotypes. The accumulated mutations differ across mutator strains, displaying base-substitution biases, allele-specific mutation hotspots, and break-associated mutation clustering. For example, in mutants of POLα and the Cdc13–Stn1–Ten1 complex, we find a distinct subtelomeric bias for mutations that we show is independent of the target sequence. Together our data suggest that specific genome-instability mutations are sufficient to drive discrete mutational signatures, some of which share properties with mutation patterns seen in tumors. Thus, in a population of cells, genome-instability mutations could influence clonal evolution by establishing discrete mutational trajectories for genomes.

Extensive variation in the density and distribution of DNA polymorphism in sorghum genomes

Sorghum genotypes currently used for grain production in the United States were developed from African landraces that were imported starting in the mid-to-late 19^th century. Farmers and plant breeders selected genotypes for grain production with reduced plant height, early flowering, increased grain yield, adaptation to drought, and improved resistance to lodging, diseases and pests. DNA polymorphisms that distinguish three historically important grain sorghum genotypes, BTx623, BTx642 and Tx7000, were characterized by genome sequencing, genotyping by sequencing, genetic mapping, and pedigree-based haplotype analysis. The distribution and density of DNA polymorphisms in the sequenced genomes varied widely, in part because the lines were derived through breeding and selection from diverse Kafir, Durra, and Caudatum race accessions. Genomic DNA spanning dw1 (SBI-09) and dw3 (SBI-07) had identical haplotypes due to selection for reduced height. Lower SNP density in genes located in pericentromeric regions compared with genes located in euchromatic regions is consistent with background selection in these regions of low recombination. SNP density was higher in euchromatic DNA and varied >100-fold in contiguous intervals that spanned up to 300 Kbp. The localized variation in DNA polymorphism density occurred throughout euchromatic regions where recombination is elevated, however, polymorphism density was not correlated with gene density or DNA methylation. Overall, sorghum chromosomes contain distal euchromatic regions characterized by extensive, localized variation in DNA polymorphism density, and large pericentromeric regions of low gene density, diversity, and recombination.

Genome duplication and mutations in ACE2 cause multicellular, fast-sedimenting phenotypes in evolved Saccharomyces cerevisiae

Laboratory evolution of the yeast Saccharomyces cerevisiae in bioreactor batch cultures yielded variants that grow as multicellular, fast-sedimenting clusters. Knowledge of the molecular basis of this phenomenon may contribute to the understanding of natural evolution of multicellularity and to manipulating cell sedimentation in laboratory and industrial applications of S. cerevisiae. Multicellular, fast-sedimenting lineages obtained from a haploid S. cerevisiae strain in two independent evolution experiments were analyzed by whole genome resequencing. The two evolved cell lines showed different frameshift mutations in a stretch of eight adenosines in ACE2, which encodes a transcriptional regulator involved in cell cycle control and mother-daughter cell separation. Introduction of the two ace2 mutant alleles into the haploid parental strain led to slow-sedimenting cell clusters that consisted of just a few cells, thus representing only a partial reconstruction of the evolved phenotype. In addition to single-nucleotide mutations, a whole-genome duplication event had occurred in both evolved multicellular strains. Construction of a diploid reference strain with two mutant ace2 alleles led to complete reconstruction of the multicellular-fast sedimenting phenotype. This study shows that whole-genome duplication and a frameshift mutation in ACE2 are sufficient to generate a fast-sedimenting, multicellular phenotype in S. cerevisiae. The nature of the ace2 mutations and their occurrence in two independent evolution experiments encompassing fewer than 500 generations of selective growth suggest that switching between unicellular and multicellular phenotypes may be relevant for competitiveness of S. cerevisiae in natural environments.

Mutation rates, spectra, and genome-wide distribution of spontaneous mutations in mismatch repair deficient yeast

DNA mismatch repair is a highly conserved DNA repair pathway. In humans, germline mutations in hMSH2 or hMLH1, key components of mismatch repair, have been associated with Lynch syndrome, a leading cause of inherited cancer mortality. Current estimates of the mutation rate and the mutational spectra in mismatch repair defective cells are primarily limited to a small number of individual reporter loci. Here we use the yeast Saccharomyces cerevisiae to generate a genome-wide view of the rates, spectra, and distribution of mutation in the absence of mismatch repair. We performed mutation accumulation assays and next generation sequencing on 19 strains, including 16 msh2 missense variants implicated in Lynch cancer syndrome. The mutation rate for DNA mismatch repair null strains was approximately 1 mutation per genome per generation, 225-fold greater than the wild-type rate. The mutations were distributed randomly throughout the genome, independent of replication timing. The mutation spectra included insertions/deletions at homopolymeric runs (87.7%) and at larger microsatellites (5.9%), as well as transitions (4.5%) and transversions (1.9%). Additionally, repeat regions with proximal repeats are more likely to be mutated. A bias toward deletions at homopolymers and insertions at (AT)_n microsatellites suggests a different mechanism for mismatch generation at these sites. Interestingly, 5% of the single base pair substitutions might represent double-slippage events that occurred at the junction of immediately adjacent repeats, resulting in a shift in the repeat boundary. These data suggest a closer scrutiny of tumor suppressors with homopolymeric runs with proximal repeats as the potential drivers of oncogenesis in mismatch repair defective cells.

A versatile microsatellite instability reporter system in human cells

Here, we report the investigation of microsatellite instability (MSI) in human cells with a newly developed reporter system based on fluorescence. We composed a vector into which microsatellites of different lengths and nucleotide composition can be introduced between a functional copy of the fluorescent protein mCherry and an out-of-frame copy of EGFP; in vivo frameshifting will lead to EGFP expression, which can be quantified by fluorescence activated cell sorting (FACS). Via targeted recombineering, single copy reporters were introduced in HEK293 and MCF-7 cells. We found predominantly -1 and +1 base pair frameshifts, the levels of which are kept in tune by mismatch repair. We show that tract length and composition greatly influences MSI. In contrast, a tracts' potential to form a G-quadruplex structure, its strand orientation or its transcriptional status is not affecting MSI. We further validated the functionality of the reporter system for screening microsatellite mutagenicity of compounds and for identifying modifiers of MSI: using a retroviral miRNA expression library, we identified miR-21, which targets MSH2, as a miRNA that induces MSI when overexpressed. Our data also provide proof of principle for the strategy of combining fluorescent reporters with next-generation sequencing technology to identify genetic factors in specific pathways.

Two mechanisms produce mutation hotspots at DNA breaks in Escherichia coli

Mutation hotspots and showers occur across phylogeny and profoundly influence genome evolution, yet the mechanisms that produce hotspots remain obscure. We report that DNA double-strand breaks (DSBs) provoke mutation hotspots via stress-induced mutation in Escherichia coli. With tet reporters placed 2 kb to 2 Mb (half the genome) away from an I-SceI site, RpoS/DinB-dependent mutations occur maximally within the first 2 kb and decrease logarithmically to ∼60 kb. A weak mutation tail extends to 1 Mb. Hotspotting occurs independently of I-site/tet-reporter-pair position in the genome, upstream and downstream in the replication path. RecD, which allows RecBCD DSB-exonuclease activity, is required for strong local but not long-distance hotspotting, indicating that double-strand resection and gap-filling synthesis underlie local hotspotting, and newly illuminating DSB resection in vivo. Hotspotting near DSBs opens the possibility that specific genomic regions could be targeted for mutagenesis, and could also promote concerted evolution (coincident mutations) within genes/gene clusters, an important issue in the evolution of protein functions.