Variation in efficiency of DNA mismatch repair at different sites in the yeast genome

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.

Mlh1 heterozygosity and promoter methylation associates with microsatellite instability in mouse sperm

DNA mismatch repair (MMR) proteins play an important role in maintaining genome stability, both in somatic and in germline cells. Loss of MLH1, a central MMR protein, leads to infertility and to microsatellite instability (MSI) in spermatocytes, however, the effect of Mlh1 heterozygosity on germline genome stability remains unexplored. To test the effect of Mlh1 heterozygosity on MSI in mature sperm, we combined mouse genetics with single-molecule PCR that detects allelic changes at unstable microsatellites. We discovered 4.5% and 5.9% MSI in sperm of 4- and 12-month-old Mlh1^+/− mice, respectively, and that Mlh1 promoter methylation in Mlh1^+/− sperm correlated with higher MSI. No such elevated MSI was seen in non-proliferating somatic cells. Additionally, we show contrasting dynamics of deletions versus insertions at unstable microsatellites (mononucleotide repeats) in sperm.

DNA mechanics and its biological impact

Almost all nucleoprotein interactions and DNA manipulation events involve mechanical deformations of DNA. Extraordinary progresses in single-molecule, structural, and computational methods have characterized the average mechanical properties of DNA, such as bendability and torsional rigidity, in high resolution. Further, the advent of sequencing technology has permitted measuring, in high-throughput, how such mechanical properties vary with sequence and epigenetic modifications along genomes. We review these recent technological advancements, and discuss how they have contributed to the emerging idea that variations in the mechanical properties of DNA play a fundamental role in regulating, genome-wide, diverse processes involved in chromatin organization.

Advances in understanding the evolution of fungal genome architecture

Diversity within the fungal kingdom is evident from the wide range of morphologies fungi display as well as the various ecological roles and industrial purposes they serve. Technological advances, particularly in long-read sequencing, coupled with the increasing efficiency and decreasing costs across sequencing platforms have enabled robust characterization of fungal genomes. These sequencing efforts continue to reveal the rampant diversity in fungi at the genome level. Here, we discuss studies that have furthered our understanding of fungal genetic diversity and genomic evolution. These studies revealed the presence of both small-scale and large-scale genomic changes. In fungi, research has recently focused on many small-scale changes, such as how hypermutation and allelic transmission impact genome evolution as well as how and why a few specific genomic regions are more susceptible to rapid evolution than others. High-throughput sequencing of a diverse set of fungal genomes has also illuminated the frequency, mechanisms, and impacts of large-scale changes, which include chromosome structural variation and changes in chromosome number, such as aneuploidy, polyploidy, and the presence of supernumerary chromosomes. The studies discussed herein have provided great insight into how the architecture of the fungal genome varies within species and across the kingdom and how modern fungi may have evolved from the last common fungal ancestor and might also pave the way for understanding how genomic diversity has evolved in all domains of life.

Effects of OsMSH6 Mutations on Microsatellite Stability and Homeologous Recombination in Rice

DNA mismatch repair (MMR) system is important for maintaining DNA replication fidelity and genome stability by repairing erroneous deletions, insertions and mis-incorporation of bases. With the aim of deciphering the role of the MMR system in genome stability and recombination in rice, we investigated the function of OsMSH6 gene, an import component of the MMR system. To achieve this goal, homeologous recombination and endogenous microsatellite stability were evaluated by using rice mutants carrying a Tos17 insertion into the OsMSH6 gene. Totally 60 microsatellites were analyzed and 15 distributed on chromosome 3, 6, 8, and 10 showed instability in three OsMSH6 mutants, D6011, NF7784 and NF9010, compared with the wild type MSH6WT (the control). The disruption of OsMSH6 gene is associated with modest increases in homeologous recombination, ranging from 2.0% to 32.5% on chromosome 1, 3, 9, and 10 in the BCF2 populations of the mutant ND6011 and NF9010. Our results suggest that the OsMSH6 plays an important role in ensuring genome stability and genetic recombination, providing the first evidence for the MSH6 gene in maintaining microsatellite stability and restricting homeologous recombination in plants.

Coordinated protein and DNA conformational changes govern mismatch repair initiation by MutS

MutS homologs identify base-pairing errors made in DNA during replication and initiate their repair. In the presence of adenosine triphosphate, MutS induces DNA bending upon mismatch recognition and subsequently undergoes conformational transitions that promote its interaction with MutL to signal repair. In the absence of MutL, these transitions lead to formation of a MutS mobile clamp that can move along the DNA. Previous single-molecule FRET (smFRET) studies characterized the dynamics of MutS DNA-binding domains during these transitions. Here, we use protein–DNA and DNA–DNA smFRET to monitor DNA conformational changes, and we use kinetic analyses to correlate DNA and protein conformational changes to one another and to the steps on the pathway to mobile clamp formation. The results reveal multiple sequential structural changes in both MutS and DNA, and they suggest that DNA dynamics play a critical role in the formation of the MutS mobile clamp. Taking these findings together with data from our previous studies, we propose a unified model of coordinated MutS and DNA conformational changes wherein initiation of mismatch repair is governed by a balance of DNA bending/unbending energetics and MutS conformational changes coupled to its nucleotide binding properties.

Scales and mechanisms of somatic mutation rate variation across the human genome

Cancer genome sequencing has revealed that somatic mutation rates vary substantially across the human genome and at scales from megabase-sized domains to individual nucleotides. Here we review recent work that has both revealed the major mutation biases that operate across the genome and the molecular mechanisms that cause them. The default mutation rate landscape in mammalian genomes results in active genes having low mutation rates because of a combination of factors that increase DNA repair: early DNA replication, transcription, active chromatin modifications and accessible chromatin. Therefore, either an increase in the global mutation rate or a redistribution of mutations from inactive to active DNA can increase the rate at which consequential mutations are acquired in active genes. Several environmental carcinogens and intrinsic mechanisms operating in tumor cells likely cause cancer by this second mechanism: by specifically increasing the mutation rate in active regions of the genome.

Cooperation between non-essential DNA polymerases contributes to genome stability in Saccharomyces cerevisiae

DNA polymerases influence genome stability through their involvement in DNA replication, response to DNA damage, and DNA repair processes. Saccharomyces cerevisiae possess four non-essential DNA polymerases, Pol λ, Pol η, Pol ζ, and Rev1, which have varying roles in genome stability. In order to assess the contribution of the non-essential DNA polymerases in genome stability, we analyzed the pol4Δ rev1Δ rev3Δ rad30Δ quadruple mutant in microhomology mediated repair, due to recent studies linking some of these DNA polymerases to this repair pathway. Our results suggest that the length and quality of microhomology influence both the overall efficiency of repair and the involvement of DNA polymerases. Furthermore, the non-essential DNA polymerases demonstrate overlapping and redundant functions when repairing double-strand breaks using short microhomologies containing mismatches. Then, we examined genome-wide mutation accumulation in the pol4Δ rev1Δ rev3Δ rad30Δ quadruple mutant compared to wild type cells. We found a significant decrease in the overall rate of mutation accumulation in the quadruple mutant cells compared to wildtype, but an increase in frameshift mutations and a shift towards transversion base-substitution with a preference for G:C to T:A or C:G. Thus, the non-essential DNA polymerases have an impact on the nature of the mutational spectrum. The sequence and functional homology shared between human and S. cerevisiae non-essential DNA polymerases suggest these DNA polymerases may have a similar role in human cells.

GC content elevates mutation and recombination rates in the yeast Saccharomyces cerevisiae

The chromosomes of many eukaryotes have regions of high GC content interspersed with regions of low GC content. In the yeast Saccharomyces cerevisiae, high-GC regions are often associated with high levels of meiotic recombination. In this study, we constructed URA3 genes that differ substantially in their base composition [URA3-AT (31% GC), URA3-WT (43% GC), and URA3-GC (63% GC)] but encode proteins with the same amino acid sequence. The strain with URA3-GC had an approximately sevenfold elevated rate of ura3 mutations compared with the strains with URA3-WT or URA3-AT About half of these mutations were single-base substitutions and were dependent on the error-prone DNA polymerase ζ. About 30% were deletions or duplications between short (5-10 base) direct repeats resulting from DNA polymerase slippage. The URA3-GC gene also had elevated rates of meiotic and mitotic recombination relative to the URA3-AT or URA3-WT genes. Thus, base composition has a substantial effect on the basic parameters of genome stability and evolution.

Molecular phylogeography of East Asian Boea clarkeana (Gesneriaceae) in relation to habitat restriction

Subfamily Cyrtandroideae (Gesneriaceae) comprises a broadly distributed group of rocky-slope herbs, with China being the center of its distributional range. The normal growth of many species within the family is particularly dependent on special habitats. Due to the paucity of molecular studies, very little is known regarding East Asian herb phylogeographic pattern. Here, we investigate the molecular phylogeography of Boea clarkeana Hemsl., a unique resurrection herb endemic to China, focusing on geographically restrictive effects of habitat distribution on evolutionary history. Variation in three chloroplast DNA (cpDNA) intergenic spacers (psbA-trnH, rps12-rpl20, and trnL-trnF), the ribosomal internal transcribed spacer (ITS) and simple sequence repeats in expressed sequence tags (EST-SSRs) was investigated across 18 populations to assess genetic diversity, genetic structure and historical dynamics. Genetic diversity was low within populations (cpDNA, hS = 0.03, πS×10(3) = 0.17; ITS, hS = 0.16, πS×10(3) = 0.43) but high for species (cpDNA, hT = 0.82, πT×10(3) = 3.12; ITS, hT = 0.88, πT×10(3) = 6.39); 76 alleles were detected in this highly inbred species (FIS = 0.22), with a significantly low average of 1.34 alleles per locus. No cpDNA or ITS haplotypes were shared between regions. Based on cpDNA results, the Mt. Huangshan-Tianmu and Mt. Qinling-Daba haplotypes are ancestral; these two regions represent potential refugia. Although no evidence of significant retreat-migration phenomena during glacial cycles was detected, interglacial range expansion from northern Mt. Qinling-Daba was identified (121,457 yr BP). Rapid agricultural growth caused bottlenecks in many populations, especially on Mt. Huang-Tianmu. Habitat restriction and fragmentation, weak seed and pollen dispersal abilities, and long-term isolation caused by human-induced or environmental changes are considered the main causes of extinction of several populations and low genetic diversity within populations and regions. These analyses clarify the effects of habitat restriction on B. clarkeana, representing an evolutionary reference for similar gesneriads, and enrich our understanding of the molecular phylogeography of East Asian rocky-slope herbs.

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.

Explosive mutation accumulation triggered by heterozygous human Pol ε proofreading-deficiency is driven by suppression of mismatch repair

Tumors defective for DNA polymerase (Pol) ε proofreading have the highest tumor mutation burden identified. A major unanswered question is whether loss of Pol ε proofreading by itself is sufficient to drive this mutagenesis, or whether additional factors are necessary. To address this, we used a combination of next generation sequencing and in vitro biochemistry on human cell lines engineered to have defects in Pol ε proofreading and mismatch repair. Absent mismatch repair, monoallelic Pol ε proofreading deficiency caused a rapid increase in a unique mutation signature, similar to that observed in tumors from patients with biallelic mismatch repair deficiency and heterozygous Pol ε mutations. Restoring mismatch repair was sufficient to suppress the explosive mutation accumulation. These results strongly suggest that concomitant suppression of mismatch repair, a hallmark of colorectal and other aggressive cancers, is a critical force for driving the explosive mutagenesis seen in tumors expressing exonuclease-deficient Pol ε.

Unique molecular mechanisms for maintenance and alteration of genetic information in the budding yeast Saccharomyces cerevisiae

The high-fidelity transmission of genetic information is crucial for the survival of organisms, the cells of which have the ability to protect DNA against endogenous and environmental agents, including reactive oxygen species (ROS), ionizing radiation, and various chemical compounds. The basis of protection mechanisms has been evolutionarily conserved from yeast to humans; however, each organism often has a specialized mode of regulation that uses different sets of machineries, particularly in lower eukaryotes. The divergence of molecular mechanisms among related organisms has provided insights into the evolution of cellular machineries to a higher architecture. Uncommon characteristics of machineries may also contribute to the development of new applications such as drugs with novel mechanisms of action. In contrast to the cellular properties for maintaining genetic information, living organisms, particularly microbes, inevitably undergo genetic alterations in order to adapt to environmental conditions. The maintenance and alteration of genetic information may be inextricably linked to each other. In this review, we describe recent findings on the unconventional molecular mechanisms of DNA damage response and DNA double-strand break (DSB) repair in the budding yeast Saccharomyces cerevisiae. We also introduce our previous research on genetic and phenotypic instabilities observed in a clonal population of clinically-derived S. cerevisiae. The molecular mechanisms of this case were associated with mutations to generate tyrosine-inserting tRNA-Tyr ochre suppressors and the position effects of mutation frequencies among eight tRNA-Tyr loci dispersed in the genome. Phenotypic variations among different strain backgrounds have also been observed by another type of nonsense suppressor, the aberrant form of the translation termination factor. Nonsense suppressors are considered to be responsible for the genome-wide translational readthrough of termination codons, including natural nonsense codons. The nonsense suppressor-mediated acquisition of phenotypic variations may be advantageous for adaptation to environmental conditions and survival during evolution.

Involvement of DNA mismatch repair in the maintenance of heterochromatic DNA stability in Saccharomyces cerevisiae

Heterochromatin contains a significant part of nuclear DNA. Little is known about the mechanisms that govern heterochromatic DNA stability. We show here that in the yeast Saccharomyces cerevisiae (i) DNA mismatch repair (MMR) is required for the maintenance of heterochromatic DNA stability, (ii) MutLα (Mlh1-Pms1 heterodimer), MutSα (Msh2-Msh6 heterodimer), MutSβ (Msh2-Msh3 heterodimer), and Exo1 are involved in MMR at heterochromatin, (iii) Exo1-independent MMR at heterochromatin frequently leads to the formation of Pol ζ-dependent mutations, (iv) MMR cooperates with the proofreading activity of Pol ε and the histone acetyltransferase Rtt109 in the maintenance of heterochromatic DNA stability, (v) repair of base-base mismatches at heterochromatin is less efficient than repair of base-base mismatches at euchromatin, and (vi) the efficiency of repair of 1-nt insertion/deletion loops at heterochromatin is similar to the efficiency of repair of 1-nt insertion/deletion loops at euchromatin.

Eukaryotic mismatch repair is an important intracellular process that defends DNA against mutations. Inactivation of mismatch repair in human cells strongly increases the risk of cancer initiation and development. Although significant progress has been made in understanding mismatch repair at euchromatin, mismatch repair at heterochromatin is not well understood. Baker’s yeast is a key model organism to study mismatch repair. We determined that in baker’s yeast (1) mismatch repair protects heterochromatic DNA from mutations, (2) the MutLα, MutSα, MutSβ, and Exo1 proteins play important roles in mismatch repair at heterochromatin, (3) Exo1-independent mismatch repair at heterochromatin is an error-prone process; (4) mismatch repair cooperates with two other intracellular processes to protect the stability of heterochromatic DNA; and (5) the efficiency of repair of base-base mismatches at heterochromatin is lower than the efficiency of repair of base-base mismatches at euchromatin, but the efficiency of 1-nt insertion/deletion loop repair at heterochromatin is similar to the efficiency of 1-nt insertion/deletion loop repair at euchromatin.

Understanding how mismatch repair proteins participate in the repair/anti-recombination decision

Mismatch repair (MMR) systems correct DNA mismatches that result from DNA polymerase misincorporation errors. Mismatches also appear in heteroduplex DNA intermediates formed during recombination between nearly identical sequences, and can be corrected by MMR or removed through an unwinding mechanism, known as anti-recombination or heteroduplex rejection. We review studies, primarily in baker's yeast, which support how specific factors can regulate the MMR/anti-recombination decision. Based on recent advances, we present models for how DNA structure, relative amounts of key repair proteins, the timely localization of repair proteins to DNA substrates and epigenetic marks can modulate this critical decision.

The Rate and Spectrum of Spontaneous Mutations in Mycobacterium smegmatis, a Bacterium Naturally Devoid of the Postreplicative Mismatch Repair Pathway

Mycobacterium smegmatis is a bacterium that is naturally devoid of known postreplicative DNA mismatch repair (MMR) homologs, mutS and mutL, providing an opportunity to investigate how the mutation rate and spectrum has evolved in the absence of a highly conserved primary repair pathway. Mutation accumulation experiments of M. smegmatis yielded a base-substitution mutation rate of 5.27 × 10⁻¹⁰ per site per generation, or 0.0036 per genome per generation, which is surprisingly similar to the mutation rate in MMR-functional unicellular organisms. Transitions were found more frequently than transversions, with the A:T→G:C transition rate significantly higher than the G:C→A:T transition rate, opposite to what is observed in most studied bacteria. We also found that the transition-mutation rate of M. smegmatis is significantly lower than that of other naturally MMR-devoid or MMR-knockout organisms. Two possible candidates that could be responsible for maintaining high DNA fidelity in this MMR-deficient organism are the ancestral-like DNA polymerase DnaE1, which contains a highly efficient DNA proofreading histidinol phosphatase (PHP) domain, and/or the existence of a uracil-DNA glycosylase B (UdgB) homolog that might protect the GC-rich M. smegmatis genome against DNA damage arising from oxidation or deamination. Our results suggest that M. smegmatis has a noncanonical Dam (DNA adenine methylase) methylation system, with target motifs differing from those previously reported. The mutation features of M. smegmatis provide further evidence that genomes harbor alternative routes for improving replication fidelity, even in the absence of major repair pathways.

Regulation of mismatch repair by histone code and posttranslational modifications in eukaryotic cells

DNA mismatch repair (MMR) protects genome integrity by correcting DNA replication-associated mispairs, modulating DNA damage-induced cell cycle checkpoints and regulating homeologous recombination. Loss of MMR function leads to cancer development. This review describes progress in understanding how MMR is carried out in the context of chromatin and how chromatin organization/compaction, epigenetic mechanisms and posttranslational modifications of MMR proteins influence and regulate MMR in eukaryotic cells.

The Rate and Molecular Spectrum of Spontaneous Mutations in the GC-Rich Multichromosome Genome of Burkholderia cenocepacia

Spontaneous mutations are ultimately essential for evolutionary change and are also the root cause of many diseases. However, until recently, both biological and technical barriers have prevented detailed analyses of mutation profiles, constraining our understanding of the mutation process to a few model organisms and leaving major gaps in our understanding of the role of genome content and structure on mutation. Here, we present a genome-wide view of the molecular mutation spectrum in Burkholderia cenocepacia, a clinically relevant pathogen with high %GC content and multiple chromosomes. We find that B. cenocepacia has low genome-wide mutation rates with insertion-deletion mutations biased toward deletions, consistent with the idea that deletion pressure reduces prokaryotic genome sizes. Unlike prior studies of other organisms, mutations in B. cenocepacia are not AT biased, which suggests that at least some genomes with high %GC content experience unusual base-substitution mutation pressure. Importantly, we also observe variation in both the rates and spectra of mutations among chromosomes and elevated G:C > T:A transversions in late-replicating regions. Thus, although some patterns of mutation appear to be highly conserved across cellular life, others vary between species and even between chromosomes of the same species, potentially influencing the evolution of nucleotide composition and genome architecture.

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.

New insights into the mechanism of DNA mismatch repair

The genome of all organisms is constantly being challenged by endogenous and exogenous sources of DNA damage. Errors like base:base mismatches or small insertions and deletions, primarily introduced by DNA polymerases during DNA replication are repaired by an evolutionary conserved DNA mismatch repair (MMR) system. The MMR system, together with the DNA replication machinery, promote repair by an excision and resynthesis mechanism during or after DNA replication, increasing replication fidelity by up-to-three orders of magnitude. Consequently, inactivation of MMR genes results in elevated mutation rates that can lead to increased cancer susceptibility in humans. In this review, we summarize our current understanding of MMR with a focus on the different MMR protein complexes, their function and structure. We also discuss how recent findings have provided new insights in the spatio-temporal regulation and mechanism of MMR.

Asymmetric Context-Dependent Mutation Patterns Revealed through Mutation-Accumulation Experiments

Despite the general assumption that site-specific mutation rates are independent of the local sequence context, a growing body of evidence suggests otherwise. To further examine context-dependent patterns of mutation, we amassed 5,645 spontaneous mutations in wild- type (WT) and mismatch-repair deficient (MMR(-)) mutation-accumulation (MA) lines of the gram-positive model organism Bacillus subtilis. We then analyzed>7,500 spontaneous base-substitution mutations across B. subtilis, Escherichia coli, and Mesoplasma florum WT and MMR(-) MA lines, finding a context-dependent mutation pattern that is asymmetric around the origin of replication. Different neighboring nucleotides can alter site-specific mutation rates by as much as 75-fold, with sites neighboring G:C base pairs or dimers involving alternating pyrimidine-purine and purine-pyrimidine nucleotides having significantly elevated mutation rates. The influence of context-dependent mutation on genome architecture is strongest in M. florum, consistent with the reduced efficiency of selection in organisms with low effective population size. If not properly accounted for, the disparities arising from patterns of context-dependent mutation can significantly influence interpretations of positive and purifying selection.

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.

Precise estimates of mutation rate and spectrum in yeast

Mutation is the ultimate source of genetic variation. The most direct and unbiased method of studying spontaneous mutations is via mutation accumulation (MA) lines. Until recently, MA experiments were limited by the cost of sequencing and thus provided us with small numbers of mutational events and therefore imprecise estimates of rates and patterns of mutation. We used whole-genome sequencing to identify nearly 1,000 spontaneous mutation events accumulated over ∼311,000 generations in 145 diploid MA lines of the budding yeast Saccharomyces cerevisiae. MA experiments are usually assumed to have negligible levels of selection, but even mild selection will remove strongly deleterious events. We take advantage of such patterns of selection and show that mutation classes such as indels and aneuploidies (especially monosomies) are proportionately much more likely to contribute mutations of large effect. We also provide conservative estimates of indel, aneuploidy, environment-dependent dominant lethal, and recessive lethal mutation rates. To our knowledge, for the first time in yeast MA data, we identified a sufficiently large number of single-nucleotide mutations to measure context-dependent mutation rates and were able to (i) confirm strong AT bias of mutation in yeast driven by high rate of mutations from C/G to T/A and (ii) detect a higher rate of mutation at C/G nucleotides in two specific contexts consistent with cytosine methylation in S. cerevisiae.

New insights and challenges in mismatch repair: getting over the chromatin hurdle

DNA mismatch repair (MMR) maintains genome stability primarily by repairing DNA replication-associated mispairs. Because loss of MMR function increases the mutation frequency genome-wide, defects in this pathway predispose affected individuals to cancer. The genes encoding essential eukaryotic MMR activities have been identified, as the recombinant proteins repair 'naked' heteroduplex DNA in vitro. However, the reconstituted system is inactive on nucleosome-containing heteroduplex DNA, and it is not understood how MMR occurs in vivo. Recent studies suggest that chromatin organization, nucleosome assembly/disassembly factors and histone modifications regulate MMR in eukaryotic cells, but the complexity and importance of the interaction between MMR and chromatin remodeling has only recently begun to be appreciated. This article reviews recent progress in understanding the mechanism of eukaryotic MMR in the context of chromatin structure and dynamics, considers the implications of these recent findings and discusses unresolved questions and challenges in understanding eukaryotic MMR.

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.

Extensive microsatellite variation in rice induced by introgression from wild rice (Zizania latifolia Griseb.)

BACKGROUND

It is widely accepted that interspecific hybridization may induce genomic instability in the resultant hybrids. However, few studies have been performed on the genomic analysis of homoploid hybrids and introgression lines. We have reported previously that by introgressive hybridization, a set of introgression lines between rice (Oryza sativa L.) and wild rice (Zizania latifolia Griseb.) was successfully generated, and which have led to the release of several cultivars.

METHODOLOGY

Using 96 microsatellite markers located in the nuclear and organelle genomes of rice, we investigated microsatellite stability in three typical introgression lines. Expression of a set of mismatch repair (MMR) genes and microsatellite-containing genes was also analyzed.

RESULTS/CONCLUSIONS

Compared with the recipient rice cultivar (Matsumae), 55 of the 96 microsatellite loci revealed variation in one or more of the introgression lines, and 58.2% of the altered alleles were shared by at least two lines, indicating that most of the alterations had occurred in the early stages of introgression before their further differentiation. 73.9% of the non-shared variations were detected only in one introgression line, i.e. RZ2. Sequence alignment showed that the variations included substitutions and indels that occurred both within the repeat tracts and in the flanking regions. Interestingly, expression of a set of MMR genes altered dramatically in the introgression lines relative to their rice parent, suggesting participation of the MMR system in the generation of microsatellite variants. Some of the altered microsatellite loci are concordant with changed expression of the genes harboring them, suggesting their possible cis-regulatory roles in controlling gene expression. Because these genes bear meaningful homology to known-functional proteins, we conclude that the introgression-induced extensive variation of microsatellites may have contributed to the novel phenotypes in the introgression lines.

Pms2 suppresses large expansions of the (GAA·TTC)n sequence in neuronal tissues

Expanded trinucleotide repeat sequences are the cause of several inherited neurodegenerative diseases. Disease pathogenesis is correlated with several features of somatic instability of these sequences, including further large expansions in postmitotic tissues. The presence of somatic expansions in postmitotic tissues is consistent with DNA repair being a major determinant of somatic instability. Indeed, proteins in the mismatch repair (MMR) pathway are required for instability of the expanded (CAG·CTG)(n) sequence, likely via recognition of intrastrand hairpins by MutSβ. It is not clear if or how MMR would affect instability of disease-causing expanded trinucleotide repeat sequences that adopt secondary structures other than hairpins, such as the triplex/R-loop forming (GAA·TTC)(n) sequence that causes Friedreich ataxia. We analyzed somatic instability in transgenic mice that carry an expanded (GAA·TTC)(n) sequence in the context of the human FXN locus and lack the individual MMR proteins Msh2, Msh6 or Pms2. The absence of Msh2 or Msh6 resulted in a dramatic reduction in somatic mutations, indicating that mammalian MMR promotes instability of the (GAA·TTC)(n) sequence via MutSα. The absence of Pms2 resulted in increased accumulation of large expansions in the nervous system (cerebellum, cerebrum, and dorsal root ganglia) but not in non-neuronal tissues (heart and kidney), without affecting the prevalence of contractions. Pms2 suppressed large expansions specifically in tissues showing MutSα-dependent somatic instability, suggesting that they may act on the same lesion or structure associated with the expanded (GAA·TTC)(n) sequence. We conclude that Pms2 specifically suppresses large expansions of a pathogenic trinucleotide repeat sequence in neuronal tissues, possibly acting independently of the canonical MMR pathway.

Identification of a novel deletion mutant strain in Saccharomyces cerevisiae that results in a microsatellite instability phenotype

The DNA mismatch repair (MMR) pathway corrects specific types of DNA replication errors that affect microsatellites and thus is critical for maintaining genomic integrity. The genes of the MMR pathway are highly conserved across different organisms. Likewise, defective MMR function universally results in microsatellite instability (MSI) which is a hallmark of certain types of cancer associated with the Mendelian disorder hereditary nonpolyposis colorectal cancer. (Lynch syndrome). To identify previously unrecognized deleted genes or loci that can lead to MSI, we developed a functional genomics screen utilizing a plasmid containing a microsatellite sequence that is a host spot for MSI mutations and the comprehensive homozygous diploid deletion mutant resource for Saccharomyces cerevisiae. This pool represents a collection of non-essential homozygous yeast diploid (2N) mutants in which there are deletions for over four thousand yeast open reading frames (ORFs). From our screen, we identified a deletion mutant strain of the PAU24 gene that leads to MSI. In a series of validation experiments, we determined that this PAU24 mutant strain had an increased MSI-specific mutation rate in comparison to the original background wildtype strain, other deletion mutants and comparable to a MMR mutant involving the MLH1 gene. Likewise, in yeast strains with a deletion of PAU24, we identified specific de novo indel mutations that occurred within the targeted microsatellite used for this screen.

Escherichia coli frameshift mutation rate depends on the chromosomal context but not on the GATC content near the mutation site

Different studies have suggested that mutation rate varies at different positions in the genome. In this work we analyzed if the chromosomal context and/or the presence of GATC sites can affect the frameshift mutation rate in the Escherichia coli genome. We show that in a mismatch repair deficient background, a condition where the mutation rate reflects the fidelity of the DNA polymerization process, the frameshift mutation rate could vary up to four times among different chromosomal contexts. Furthermore, the mismatch repair efficiency could vary up to eight times when compared at different chromosomal locations, indicating that detection and/or repair of frameshift events also depends on the chromosomal context. Also, GATC sequences have been proved to be essential for the correct functioning of the E. coli mismatch repair system. Using bacteriophage heteroduplexes molecules it has been shown that GATC influence the mismatch repair efficiency in a distance- and number-dependent manner, being almost nonfunctional when GATC sequences are located at 1 kb or more from the mutation site. Interestingly, we found that in E. coli genomic DNA the mismatch repair system can efficiently function even if the nearest GATC sequence is located more than 2 kb away from the mutation site. The results presented in this work show that even though frameshift mutations can be efficiently generated and/or repaired anywhere in the genome, these processes can be modulated by the chromosomal context that surrounds the mutation site.

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

Evolutionary theory assumes that mutations occur randomly in the genome; however, studies performed in a variety of organisms indicate the existence of context-dependent mutation biases. Sources of mutagenesis variation across large genomic contexts (e.g., hundreds of bases) have not been identified. Here, we use high-coverage whole-genome sequencing of a conditional mismatch repair mutant line of diploid yeast to identify mutations that accumulated after 160 generations of growth. The vast majority of the mutations accumulated as insertion/deletions (in/dels) in homopolymeric [poly(dA:dT)] and repetitive DNA tracts. Surprisingly, the likelihood of an in/del mutation in a given poly(dA:dT) tract is increased by the presence of nearby poly(dA:dT) tracts in up to a 1,000 bp region centered on the given tract. Our work suggests that specific mutation hot spots can contribute disproportionately to the genetic variation that is introduced into populations and provides long-range genomic sequence context that contributes to mutagenesis.

Overview for the histone codes for DNA repair

DNA damage occurs continuously as a result of various factors-intracellular metabolism, replication, and exposure to genotoxic agents, such as ionizing radiation and chemotherapy. If left unrepaired, this damage could result in changes or mutations within the cell genomic material. There are a number of different pathways that the cell can utilize to repair these DNA breaks. However, it is of utmost interest to know how the DNA damage is signaled to the various DNA pathways. As DNA damage occurs within the chromatin, we postulate that modifications of histones are important for signaling the position of DNA damage, recruiting the DNA repair proteins to the site of damage, and creating an open structure such that the repair proteins can access the site of damage. We discuss the modifications that occur on the histones and the manner in which they relate to the type of damage that has occurred as well as the DNA repair pathways that are activated.

Mutator suppression and escape from replication error-induced extinction in yeast

Cells rely on a network of conserved pathways to govern DNA replication fidelity. Loss of polymerase proofreading or mismatch repair elevates spontaneous mutation and facilitates cellular adaptation. However, double mutants are inviable, suggesting that extreme mutation rates exceed an error threshold. Here we combine alleles that affect DNA polymerase δ (Pol δ) proofreading and mismatch repair to define the maximal error rate in haploid yeast and to characterize genetic suppressors of mutator phenotypes. We show that populations tolerate mutation rates 1,000-fold above wild-type levels but collapse when the rate exceeds 10⁻³ inactivating mutations per gene per cell division. Variants that escape this error-induced extinction (eex) rapidly emerge from mutator clones. One-third of the escape mutants result from second-site changes in Pol δ that suppress the proofreading-deficient phenotype, while two-thirds are extragenic. The structural locations of the Pol δ changes suggest multiple antimutator mechanisms. Our studies reveal the transient nature of eukaryotic mutators and show that mutator phenotypes are readily suppressed by genetic adaptation. This has implications for the role of mutator phenotypes in cancer.

Organisms strike a balance between genetic continuity and change. Most cells are well adapted to their niches and therefore invest heavily in mechanisms that maintain accurate DNA replication. When cell populations are confronted with changing environmental conditions, “mutator” clones with high mutation rates emerge and readily adapt to the new conditions by rapidly acquiring beneficial mutations. However, deleterious mutations also accumulate, raising the question: what level of mutational burden can cell populations sustain before collapsing? Here we experimentally determine the maximal mutation rate in haploid yeast. We observe that yeast can withstand a 1,000-fold increase in mutation rate without losing colony forming capacity. Yet no strains survive a 10,000-fold increase in mutation rate. Escape mutants with an “anti-mutator” phenotype frequently emerge from cell populations undergoing this error-induced extinction. The diversity of antimutator changes suggests that strong mutator phenotypes in nature may be inherently transient, ensuring that rapid adaptation is followed by genetic attenuation which preserves the beneficial, adaptive mutations. These observations are relevant to microbial populations during infection as well as the somatic evolution of cancer cells.

Mutation rates across budding yeast chromosome VI are correlated with replication timing

Previous experimental studies suggest that the mutation rate is nonuniform across the yeast genome. To characterize this variation across the genome more precisely, we measured the mutation rate of the URA3 gene integrated at 43 different locations tiled across Chromosome VI. We show that mutation rate varies 6-fold across a single chromosome, that this variation is correlated with replication timing, and we propose a model to explain this variation that relies on the temporal separation of two processes for replicating past damaged DNA: error-free DNA damage tolerance and translesion synthesis. This model is supported by the observation that eliminating translesion synthesis decreases this variation.

Evidence of association between nucleosome occupancy and the evolution of transcription factor binding sites in yeast

Background

Divergence of transcription factor binding sites is considered to be an important source of regulatory evolution. The associations between transcription factor binding sites and phenotypic diversity have been investigated in many model organisms. However, the understanding of other factors that contribute to it is still limited. Recent studies have elucidated the effect of chromatin structure on molecular evolution of genomic DNA. Though the profound impact of nucleosome positions on gene regulation has been reported, their influence on transcriptional evolution is still less explored. With the availability of genome-wide nucleosome map in yeast species, it is thus desirable to investigate their impact on transcription factor binding site evolution. Here, we present a comprehensive analysis of the role of nucleosome positioning in the evolution of transcription factor binding sites.

Results

We compared the transcription factor binding site frequency in nucleosome occupied regions and nucleosome depleted regions in promoters of old (orthologs among Saccharomycetaceae) and young (Saccharomyces specific) genes; and in duplicate gene pairs. We demonstrated that nucleosome occupied regions accommodate greater binding site variations than nucleosome depleted regions in young genes and in duplicate genes. This finding was confirmed by measuring the difference in evolutionary rates of binding sites in sensu stricto yeasts at nucleosome occupied regions and nucleosome depleted regions. The binding sites at nucleosome occupied regions exhibited a consistently higher evolution rate than those at nucleosome depleted regions, corroborating the difference in the selection constraints at the two regions. Finally, through site-directed mutagenesis experiment, we found that binding site gain or loss events at nucleosome depleted regions may cause more expression differences than those in nucleosome occupied regions.

Conclusions

Our study indicates the existence of different selection constraint on binding sites at nucleosome occupied regions than at the nucleosome depleted regions. We found that the binding sites have a different rate of evolution at nucleosome occupied and depleted regions. Finally, using transcription factor binding site-directed mutagenesis experiment, we confirmed the difference in the impact of binding site changes on expression at these regions. Thus, our work demonstrates the importance of composite analysis of chromatin and transcriptional evolution.

Context-dependent codon partition models provide significant increases in model fit in atpB and rbcL protein-coding genes

BACKGROUND

Accurate modelling of substitution processes in protein-coding sequences is often hampered by the computational burdens associated with full codon models. Lately, codon partition models have been proposed as a viable alternative, mimicking the substitution behaviour of codon models at a low computational cost. Such codon partition models however impose independent evolution of the different codon positions, which is overly restrictive from a biological point of view. Given that empirical research has provided indications of context-dependent substitution patterns at four-fold degenerate sites, we take those indications into account in this paper.

RESULTS

We present so-called context-dependent codon partition models to assess previous empirical claims that the evolution of four-fold degenerate sites is strongly dependent on the composition of its two flanking bases. To this end, we have estimated and compared various existing independent models, codon models, codon partition models and context-dependent codon partition models for the atpB and rbcL genes of the chloroplast genome, which are frequently used in plant systematics. Such context-dependent codon partition models employ a full dependency scheme for four-fold degenerate sites, whilst maintaining the independence assumption for the first and second codon positions.

CONCLUSIONS

We show that, both in the atpB and rbcL alignments of a collection of land plants, these context-dependent codon partition models significantly improve model fit over existing codon partition models. Using Bayes factors based on thermodynamic integration, we show that in both datasets the same context-dependent codon partition model yields the largest increase in model fit compared to an independent evolutionary model. Context-dependent codon partition models hence perform closer to codon models, which remain the best performing models at a drastically increased computational cost, compared to codon partition models, but remain computationally interesting alternatives to codon models. Finally, we observe that the substitution patterns in both datasets are drastically different, leading to the conclusion that combined analysis of these two genes using a single model may not be advisable from a context-dependent point of view.

Factors affecting germline mutations in a hypervariable microsatellite: a comparative analysis of six species of swallows (Aves: Hirundinidae)

Microsatellites mutate frequently by replication slippage. Empirical evidence shows that the probability of such slippage mutations may increase with the length of the repeat region as well as exposure to environmental mutagens, but the mutation rate can also differ between the male and female germline. It has been hypothesized that more intense sexual selection or sperm competition can also lead to elevated mutation rates, but the empirical evidence is inconclusive. Here, we analyzed the occurrence of germline slippage mutations in the hypervariable pentanucleotide microsatellite locus HrU10 across six species of swallow (Aves: Hirundinidae). These species exhibit marked differences in the length range of the microsatellite, as well as differences in the intensity of sperm competition. We found a strong effect of microsatellite length on the probability of mutation, but no residual effect of species or their level of sperm competition when the length effect was accounted for. Neither could we detect any difference in mutation rate between tree swallows (Tachycineta bicolor) breeding in Hamilton Harbour, Ontario, an industrial site with previous documentation of elevated mutation rates for minisatellite DNA, and a rural reference population. However, our cross-species analysis revealed two significant patterns of sex differences in HrU10 germline mutations: (1) mutations in longer alleles occurred typically in the male germline, those in shorter alleles in the female germline, and (2) male germline mutations were more often expansions than contractions, whereas no directional bias was evident in the female germline. These results indicate some fundamental differences in male and female gametogenesis affecting the probability of slippage mutations. Our study also reflects the value of a comparative, multi-species approach for locus-specific mutation analyses, through which a wider range of influential factors can be assessed than in single-species studies.

Break-induced replication is highly inaccurate

DNA replication initiated by one-ended homologous recombination at a double-strand break is highly inaccurate, as it greatly stimulates frameshift mutations over the entire path of the replication fork.

DNA must be synthesized for purposes of genome duplication and DNA repair. While the former is a highly accurate process, short-patch synthesis associated with repair of DNA damage is often error-prone. Break-induced replication (BIR) is a unique cellular process that mimics normal DNA replication in its processivity, rate, and capacity to duplicate hundreds of kilobases, but is initiated at double-strand breaks (DSBs) rather than at replication origins. Here we employed a series of frameshift reporters to measure mutagenesis associated with BIR in Saccharomyces cerevisiae. We demonstrate that BIR DNA synthesis is intrinsically inaccurate over the entire path of the replication fork, as the rate of frameshift mutagenesis during BIR is up to 2,800-fold higher than during normal replication. Importantly, this high rate of mutagenesis was observed not only close to the DSB where BIR is less stable, but also far from the DSB where the BIR replication fork is fast and stabilized. We established that polymerase proofreading and mismatch repair correct BIR errors. Also, dNTP levels were elevated during BIR, and this contributed to BIR-related mutagenesis. We propose that a high level of DNA polymerase errors that is not fully compensated by error-correction mechanisms is largely responsible for mutagenesis during BIR, with Pol δ generating many of the mutagenic errors. We further postulate that activation of BIR in eukaryotic cells may significantly contribute to accumulation of mutations that fuel cancer and evolution.

Accurate transmission of genetic information requires the precise replication of parental DNA. Mutations (which can be beneficial or deleterious) arise from errors that remain uncorrected. DNA replication occurs during S-phase of the cell cycle and is extremely accurate due to highly selective DNA polymerases coupled with effective error-correction mechanisms. In contrast, DNA synthesis associated with short-patch DNA repair is often error-prone. Break-induced replication (BIR) presents an interesting case of large-scale DNA duplication that occurs in the context of DNA repair. In this study we employed a yeast-based system to investigate the level of mutagenesis associated with BIR compared to mutagenesis during normal DNA replication. We report that frameshifts, which are the most deleterious kind of point mutation, are much more frequent during BIR than during normal DNA replication. Surprisingly, we observed that the majority of mutations associated with BIR were created by polymerases responsible for normal DNA replication, which are assumed to be highly precise. Overall, we propose that BIR is a novel source of mutagenesis that may contribute to disease genesis and evolution.

Simple sequence repeat variation in the Daphnia pulex genome

Background

Simple sequence repeats (SSRs) are highly variable features of all genomes. Their rapid evolution makes them useful for tracing the evolutionary history of populations and investigating patterns of selection and mutation across gnomes. The recently sequenced Daphnia pulex genome provides us with a valuable data set to study the mode and tempo of SSR evolution, without the inherent biases that accompany marker selection.

Results

Here we catalogue SSR loci in the Daphnia pulex genome with repeated motif sizes of 1-100 nucleotides with a minimum of 3 perfect repeats. We then used whole genome shotgun reads to determine the average heterozygosity of each SSR type and the relationship that it has to repeat number, motif size, motif sequence, and distribution of SSR loci. We find that SSR heterozygosity is motif specific, and positively correlated with repeat number as well as motif size. For non-repeat unit polymorphisms, we identify a motif-dependent end-nucleotide polymorphism bias that may contribute to the patterns of abundance for specific homopolymers, dimers, and trimers. Our observations confirm the high frequency of multiple unit variation (multistep) at large microsatellite loci, and further show that the occurrence of multiple unit variation is dependent on both repeat number and motif size. Using the Daphnia pulex genetic map, we show a positive correlation between dimer and trimer frequency and recombination.

Conclusions

This genome-wide analysis of SSR variation in Daphnia pulex indicates that several aspects of SSR variation are motif dependent and suggests that a combination of unit length variation and end repeat biased base substitution contribute to the unique spectrum of SSR repeat loci.

Genome-wide model for the normal eukaryotic DNA replication fork

To investigate DNA replication enzymology across the nuclear genome of budding yeast, deep sequencing was used to establish the pattern of uncorrected replication errors generated by an asymmetric mutator variant of DNA polymerase δ (Pol δ). Sequencing of 16 genomes identified 1,206-bp substitutions generated over 33 generations by L612M Pol δ in a mismatch repair defective strain. Alignment of sequences flanking these substitutions identified "hotspot" motifs for Pol δ replication errors. The substitutions were distributed evenly across all 16 chromosomes. The vast majority were transitions that occurred with a strand bias that varied in a predictable manner relative to known functional origins of replication. This strand bias strongly supports the idea that Pol δ is primarily a lagging strand polymerase during replication across the entire nuclear genome.

Detection of heterozygous mutations in the genome of mismatch repair defective diploid yeast using a Bayesian approach

DNA replication errors that escape polymerase proofreading and mismatch repair (MMR) can lead to base substitution and frameshift mutations. Such mutations can disrupt gene function, reduce fitness, and promote diseases such as cancer and are also the raw material of molecular evolution. To analyze with limited bias genomic features associated with DNA polymerase errors, we performed a genome-wide analysis of mutations that accumulate in MMR-deficient diploid lines of Saccharomyces cerevisiae. These lines were derived from a common ancestor and were grown for 160 generations, with bottlenecks reducing the population to one cell every 20 generations. We sequenced to between 8- and 20-fold coverage one wild-type and three mutator lines using Illumina Solexa 36-bp reads. Using an experimentally aware Bayesian genotype caller developed to pool experimental data across sequencing runs for all strains, we detected 28 heterozygous single-nucleotide polymorphisms (SNPs) and 48 single-nt insertion/deletions (indels) from the data set. This method was evaluated on simulated data sets and found to have a very low false-positive rate (∼6 × 10(-5)) and a false-negative rate of 0.08 within the unique mapping regions of the genome that contained at least sevenfold coverage. The heterozygous mutations identified by the Bayesian genotype caller were confirmed by Sanger sequencing. All of the mutations were unique to a given line, except for a single-nt deletion mutation which occurred independently in two lines. All 48 indels, composed of 46 deletions and two insertions, occurred in homopolymer (HP) tracts [i.e., 47 poly(A) or (T) tracts, 1 poly(G) or (C) tract] between 5 and 13 bp long. Our findings are of interest because HP tracts are present at high levels in the yeast genome (>77,400 for 5- to 20-nt HP tracts), and frameshift mutations in these regions are likely to disrupt gene function. In addition, they demonstrate that the mutation pattern seen previously in mismatch repair defective strains using a limited number of reporters holds true for the entire genome.

Defective mismatch repair, microsatellite mutation bias, and variability in clinical cancer phenotypes

Microsatellite instability is associated with 10% to 15% of colorectal, endometrial, ovarian, and gastric cancers, and has long been used as a diagnostic tool for hereditary nonpolyposis colorectal carcinoma-related cancers. Tumor-specific length alterations within microsatellites are generally accepted to be a consequence of strand slippage events during DNA replication, which are uncorrected due to a defective postreplication mismatch repair (MMR) system. Mutations arising within microsatellites associated with critical target genes are believed to play a causative role in the evolution of MMR-defective tumors. In this review, we summarize current evidence of mutational biases within microsatellites arising as a consequence of intrinsic DNA sequence effects as well as variation in MMR efficiency. Microsatellite mutational biases are generally not considered during clinical testing; however, we suggest that such biases may be clinically significant as a factor contributing to phenotypic variation among microsatellite instability-positive tumors.

Genomic mutation rates: what high-throughput methods can tell us

High-throughput DNA analyses are increasingly being used to detect rare mutations in moderately sized genomes. These methods have yielded genome mutation rates that are markedly higher than those obtained using pre-genomic strategies. Recent work in a variety of organisms has shown that mutation rate is strongly affected by sequence context and genome position. These observations suggest that high-throughput DNA analyses will ultimately allow researchers to identify trans-acting factors and cis sequences that underlie mutation rate variation. Such work should provide insights on how mutation rate variability can impact genome organization and disease progression.

Evidence that nucleosomes inhibit mismatch repair in eukaryotic cells

The influence of chromatin structure on DNA metabolic processes, including DNA replication and repair, has been a matter of intensive studies in recent years. Although the human mismatch repair (MMR) reaction has been reconstituted using purified proteins, the influence of chromatin structure on human MMR is unknown. This study examines the interaction between human MutSalpha and a mismatch located within a nucleosome or between two nucleosomes. The results show that, whereas MutSalpha specifically recognizes both types of nucleosomal heteroduplexes, the protein bound the mismatch within a nucleosome with much lower efficiency than a naked heteroduplex or a heterology free of histone proteins but between two nucleosomes. Additionally, MutSalpha displays reduced ATPase- and ADP-binding activity when interacting with nucleosomal heteroduplexes. Interestingly, nucleosomes block ATP-induced MutSalpha sliding along the DNA helix when the mismatch is in between two nucleosomes. These findings suggest that nucleosomes may inhibit MMR in eukaryotic cells. The implications of these findings for our understanding of eukaryotic MMR are discussed.

The polymerase eta translesion synthesis DNA polymerase acts independently of the mismatch repair system to limit mutagenesis caused by 7,8-dihydro-8-oxoguanine in yeast

Reactive oxygen species are ubiquitous mutagens that have been linked to both disease and aging. The most studied oxidative lesion is 7,8-dihydro-8-oxoguanine (GO), which is often miscoded during DNA replication, resulting specifically in GC --> TA transversions. In yeast, the mismatch repair (MMR) system repairs GO.A mismatches generated during DNA replication, and the polymerase eta (Poleta) translesion synthesis DNA polymerase additionally promotes error-free bypass of GO lesions. It has been suggested that Poleta limits GO-associated mutagenesis exclusively through its participation in the filling of MMR-generated gaps that contain GO lesions. In the experiments reported here, the SUP4-o forward-mutation assay was used to monitor GC --> TA mutation rates in strains defective in MMR (Msh2 or Msh6) and/or in Poleta activity. The results clearly demonstrate that Poleta can function independently of the MMR system to prevent GO-associated mutations, presumably through preferential insertion of cytosine opposite replication-blocking GO lesions. Furthermore, the Poleta-dependent bypass of GO lesions is more efficient on the lagging strand of replication and requires an interaction with proliferating cell nuclear antigen. These studies establish a new paradigm for the prevention of GO-associated mutagenesis in eukaryotes.

From bacteria to plants: a compendium of mismatch repair assays

Mismatch repair (MMR) system maintains genome integrity by correcting mispaired or unpaired bases which have escaped the proofreading activity of DNA polymerases. The basic features of the pathway have been highly conserved throughout evolution, although the nature and number of the proteins involved in the mechanism vary from prokaryotes to eukaryotes and even between humans and plants. Cells deficient in MMR genes have been observed to display a mutator phenotype characterized by an increased rate in spontaneous mutation, instability of microsatellite sequences and illegitimate recombination between diverged DNA sequences. Studies of the mutator phenotype have demonstrated a critical role for the MMR system in mutation avoidance and genetic stability. Here, we briefly review our current knowledge of the MMR mechanism and then focus on the in vivo biochemical and genetic assays used to investigate the function of the MMR proteins in processing DNA mismatches generated during replication and mitotic recombination in Escherichia coli, Saccharomyces cerevisiae, Homo sapiens and Arabidopsis thaliana. An overview of the biochemical assays developed to study mismatch correction in vitro is also provided.

Every microsatellite is different: Intrinsic DNA features dictate mutagenesis of common microsatellites present in the human genome

Microsatellite sequences are ubiquitous in the human genome and are important regulators of genome function. Here, we examine the mutational mechanisms governing the stability of highly abundant mono-, di-, and tetranucleotide microsatellites. Microsatellite mutation rate estimates from pedigree analyses and experimental models range from a low of approximately 10(-6) to a high of approximately 10(-2) mutations per locus per generation. The vast majority of observed mutational variation can be attributed to features intrinsic to the allele itself, including motif size, length, and sequence composition. A greater than linear relationship between motif length and mutagenesis has been observed in several model systems. Motif sequence differences contribute up to 10-fold to the variation observed in human cell mutation rates. The major mechanism of microsatellite mutagenesis is strand slippage during DNA synthesis. DNA polymerases produce errors within microsatellites at a frequency that is 10- to 100-fold higher than the frequency of frameshifts in coding sequences. Motif sequence significantly affects both polymerase error rate and specificity, resulting in strand biases within complementary microsatellites. Importantly, polymerase errors within microsatellites include base substitutions, deletions, and complex mutations, all of which produced interrupted alleles from pure microsatellites. Postreplication mismatch repair efficiency is affected by microsatellite motif size and sequence, also contributing to the observed variation in microsatellite mutagenesis. Inhibition of DNA synthesis within common microsatellites is highly sequence-dependent, and is positively correlated with the production of errors. DNA secondary structure within common microsatellites can account for some DNA polymerase pause sites, and may be an important factor influencing mutational specificity.

Analysis of genomic instability in adult-onset celiac disease patients by microsatellite instability and loss of heterozygosis

BACKGROUND AND AIMS

Malignant complications of celiac disease (CD) include carcinomas and lymphomas. The genetic basis behind cancer development in CD is not known, but acquisition of genetic abnormalities and genomic instability has been involved. The aim of this study was to explore molecular characteristics of genomic instability in CD patients by analyzing microsatellite instability (MSI) and loss of heterozygosis (LOH) with carefully selected microsatellites.

METHODS

We genotyped small bowel biopsies and peripheral blood samples from 20 untreated CD patients using five microsatellites related to MMR genes (panel A), and five repeats associated with tumor suppressor genes, chromosome instability, inflammation, and cancer (panel B).

RESULTS

Genomic instability was found in seven out of 20 (35%) cases at: D5S107, D18S58, GSTP, TP53 or DCC, being TP53 the most frequently affected (five out of seven cases; 71%). Microsatellite alterations were significantly found using panel B markers (P=0.04). No cases with high frequency of MSI and replication error phenotype were detected. Only one case displayed MSI-L alone. Three patients exhibited LOH and three other cases showed LOH with low level of MSI, being classified as having chromosome instability phenotype.

CONCLUSION

Two novel observations were found in this study: first, the finding that non-neoplastic cells from a group of untreated CD patients present genomic instability at nucleotide level; and second, the advantage to use carefully selected microsatellites to identify celiac patients with molecular instability. Our data support the existence of chromosome instability phenotype in CD, suggesting that stable and unstable patients are genomically distinct subtypes that may follow a different evolution.

Comparative genomics and molecular dynamics of DNA repeats in eukaryotes

Repeated elements can be widely abundant in eukaryotic genomes, composing more than 50% of the human genome, for example. It is possible to classify repeated sequences into two large families, "tandem repeats" and "dispersed repeats." Each of these two families can be itself divided into subfamilies. Dispersed repeats contain transposons, tRNA genes, and gene paralogues, whereas tandem repeats contain gene tandems, ribosomal DNA repeat arrays, and satellite DNA, itself subdivided into satellites, minisatellites, and microsatellites. Remarkably, the molecular mechanisms that create and propagate dispersed and tandem repeats are specific to each class and usually do not overlap. In the present review, we have chosen in the first section to describe the nature and distribution of dispersed and tandem repeats in eukaryotic genomes in the light of complete (or nearly complete) available genome sequences. In the second part, we focus on the molecular mechanisms responsible for the fast evolution of two specific classes of tandem repeats: minisatellites and microsatellites. Given that a growing number of human neurological disorders involve the expansion of a particular class of microsatellites, called trinucleotide repeats, a large part of the recent experimental work on microsatellites has focused on these particular repeats, and thus we also review the current knowledge in this area. Finally, we propose a unified definition for mini- and microsatellites that takes into account their biological properties and try to point out new directions that should be explored in a near future on our road to understanding the genetics of repeated sequences.