Yeast ARMs (DNA at-risk motifs) can reveal sources of genome instability

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

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.

Site-specific MCM sumoylation prevents genome rearrangements by controlling origin-bound MCM

Timely completion of eukaryotic genome duplication requires coordinated DNA replication initiation at multiple origins. Replication begins with the loading of the Mini-Chromosome Maintenance (MCM) complex, proceeds by the activation of the Cdc45-MCM-GINS (CMG) helicase, and ends with CMG removal after chromosomes are fully replicated. Post-translational modifications on the MCM and associated factors ensure an orderly transit of these steps. Although the mechanisms of CMG activation and removal are partially understood, regulated MCM loading is not, leaving an incomplete understanding of how DNA replication begins. Here we describe a site-specific modification of Mcm3 by the Small Ubiquitin-like MOdifier (SUMO). Mutations that prevent this modification reduce the MCM loaded at replication origins and lower CMG levels, resulting in impaired cell growth, delayed chromosomal replication, and the accumulation of gross chromosomal rearrangements (GCRs). These findings demonstrate the existence of a SUMO-dependent regulation of origin-bound MCM and show that this pathway is needed to prevent genome rearrangements.

Faithful replication of the genome is essential for the survival and health of all living organisms. The eukaryotic genome presents a unique and difficult challenge: its enormous size demands the coordinated action of numerous DNA replication origins to ensure timely completion of genome duplication. Although the mechanisms that control the activation and removal of DNA replisome are partially understood, whether and how cells regulate the loading of the Mini-Chromosome Maintenance (MCM) complex, the precursor of the DNA replisome, at replication origins are not. Because mutations to MCM-loading factors and enzymes that catalyze reversible protein sumoylation cause substantial gross chromosomal rearrangements (GCRs) that characterize the cancer genome, understanding regulated MCM loading is one of the most pressing questions in the field. Here, we identified a site-specific SUMO modification of MCM and found that mutation disabling this modification causes severe growth defect and impaired DNA replication. These defects are attributable to reduced MCM at DNA replication origins, resulting in a lower DNA replisome level and a dramatic accumulation of GCRs. Thus, these findings identify a hitherto unknown regulatory mechanism: Site-specific MCM sumoylation regulates origin-bound MCM, and this prevents genome rearrangements.

UV-exposure, endogenous DNA damage, and DNA replication errors shape the spectra of genome changes in human skin

Human skin is continuously exposed to environmental DNA damage leading to the accumulation of somatic mutations over the lifetime of an individual. Mutagenesis in human skin cells can be also caused by endogenous DNA damage and by DNA replication errors. The contributions of these processes to the somatic mutation load in the skin of healthy humans has so far not been accurately assessed because the low numbers of mutations from current sequencing methodologies preclude the distinction between sequencing errors and true somatic genome changes. In this work, we sequenced genomes of single cell-derived clonal lineages obtained from primary skin cells of a large cohort of healthy individuals across a wide range of ages. We report here the range of mutation load and a comprehensive view of the various somatic genome changes that accumulate in skin cells. We demonstrate that UV-induced base substitutions, insertions and deletions are prominent even in sun-shielded skin. In addition, we detect accumulation of mutations due to spontaneous deamination of methylated cytosines as well as insertions and deletions characteristic of DNA replication errors in these cells. The endogenously induced somatic mutations and indels also demonstrate a linear increase with age, while UV-induced mutation load is age-independent. Finally, we show that DNA replication stalling at common fragile sites are potent sources of gross chromosomal rearrangements in human cells. Thus, somatic mutations in skin of healthy individuals reflect the interplay of environmental and endogenous factors in facilitating genome instability and carcinogenesis.

Skin forms the first barrier against a variety of environmental toxins and DNA damaging agents. Additionally, DNA of skin cells suffer from endogenous damage and errors during replication. Altogether, these lesions cause a variety of genome changes resulting in disease including cancer. However, the accurate measurement of the range and complete spectrum of genome changes in healthy skin was missing due to technical or biological limitations of prior studies. We present here accurate measurements of the various types of somatic genome changes that we found in skin fibroblasts and melanocytes from 21 donors ranging in ages from 25 to 79 years, which allowed to distinguish age related from age independent changes. Our cohort contains both White and African American donors, allowing an estimation of the impacts of skin color on mutagenesis. As a result, we revealed the complete spectrum and determined the range of somatic genome changes and their etiologies in healthy human skin fibroblasts and melanocytes and highlighted molecular mechanisms underlying these changes. Therefore, our study introduces a base line for defining disease levels of genome instability in skin.

Replication Stress and Consequential Instability of the Genome and Epigenome

Cells must faithfully duplicate their DNA in the genome to pass their genetic information to the daughter cells. To maintain genomic stability and integrity, double-strand DNA has to be replicated in a strictly regulated manner, ensuring the accuracy of its copy number, integrity and epigenetic modifications. However, DNA is constantly under the attack of DNA damage, among which oxidative DNA damage is the one that most frequently occurs, and can alter the accuracy of DNA replication, integrity and epigenetic features, resulting in DNA replication stress and subsequent genome and epigenome instability. In this review, we summarize DNA damage-induced replication stress, the formation of DNA secondary structures, peculiar epigenetic modifications and cellular responses to the stress and their impact on the instability of the genome and epigenome mainly in eukaryotic cells.

Pathways and Mechanisms that Prevent Genome Instability in Saccharomyces cerevisiae

Genome rearrangements result in mutations that underlie many human diseases, and ongoing genome instability likely contributes to the development of many cancers. The tools for studying genome instability in mammalian cells are limited, whereas model organisms such as Saccharomyces cerevisiae are more amenable to these studies. Here, we discuss the many genetic assays developed to measure the rate of occurrence of Gross Chromosomal Rearrangements (called GCRs) in S. cerevisiae. These genetic assays have been used to identify many types of GCRs, including translocations, interstitial deletions, and broken chromosomes healed by de novo telomere addition, and have identified genes that act in the suppression and formation of GCRs. Insights from these studies have contributed to the understanding of pathways and mechanisms that suppress genome instability and how these pathways cooperate with each other. Integrated models for the formation and suppression of GCRs are discussed.

Mechanisms of germ line genome instability

During meiosis, numerous DNA double-strand breaks (DSBs) are formed as part of the normal developmental program. This seemingly destructive behavior is necessary for successful meiosis, since repair of the DSBs through homologous recombination (HR) helps to produce physical links between the homologous chromosomes essential for correct chromosome segregation later in meiosis. However, DSB formation at such a massive scale also introduces opportunities to generate gross chromosomal rearrangements. In this review, we explore ways in which meiotic DSBs can result in such genomic alterations.

Yeast Assay Highlights the Intrinsic Genomic Instability of Human PML Intron 6 over Intron 3 and the Role of Replication Fork Proteins

Human acute promyelocytic leukemia (APL) is characterized by a specific balanced translocation t(15;17)(q22;q21) involving the PML and RARA genes. In both de novo and therapy-related APL, the most frequent PML breakpoints are located within intron 6, and less frequently in intron 3; the precise mechanisms by which these breakpoints arise and preferentially in PML intron 6 remain unsolved. To investigate the intrinsic properties of the PML intron sequences in vivo, we designed Saccharomyces cerevisiae strains containing human PML intron 6 or intron 3 sequences inserted in yeast chromosome V and measured gross chromosomal rearrangements (GCR). This approach provided evidence that intron 6 had a superior instability over intron 3 due to an intrinsic property of the sequence and identified the 3’ end of intron 6 as the most susceptible to break. Using yeast strains invalidated for genes that control DNA replication, we show that this differential instability depended at least upon Rrm3, a DNA helicase, and Mrc1, the human claspin homolog. GCR induction by hydrogen peroxide, a general genotoxic agent, was also dependent on genetic context. We conclude that: 1) this yeast system provides an alternative approach to study in detail the properties of human sequences in a genetically controlled situation and 2) the different susceptibility to produce DNA breaks in intron 6 versus intron 3 of the human PML gene is likely due to an intrinsic property of the sequence and is under replication fork genetic control.

Distinct SUMO ligases cooperate with Esc2 and Slx5 to suppress duplication-mediated genome rearrangements

Suppression of duplication-mediated gross chromosomal rearrangements (GCRs) is essential to maintain genome integrity in eukaryotes. Here we report that SUMO ligase Mms21 has a strong role in suppressing GCRs in Saccharomyces cerevisiae, while Siz1 and Siz2 have weaker and partially redundant roles. Understanding the functions of these enzymes has been hampered by a paucity of knowledge of their substrate specificity in vivo. Using a new quantitative SUMO-proteomics technology, we found that Siz1 and Siz2 redundantly control the abundances of most sumoylated substrates, while Mms21 more specifically regulates sumoylation of RNA polymerase-I and the SMC-family proteins. Interestingly, Esc2, a SUMO-like domain-containing protein, specifically promotes the accumulation of sumoylated Mms21-specific substrates and functions with Mms21 to suppress GCRs. On the other hand, the Slx5-Slx8 complex, a SUMO-targeted ubiquitin ligase, suppresses the accumulation of sumoylated Mms21-specific substrates. Thus, distinct SUMO ligases work in concert with Esc2 and Slx5-Slx8 to control substrate specificity and sumoylation homeostasis to prevent GCRs.

High-resolution mapping of spontaneous mitotic recombination hotspots on the 1.1 Mb arm of yeast chromosome IV

Although homologous recombination is an important pathway for the repair of double-stranded DNA breaks in mitotically dividing eukaryotic cells, these events can also have negative consequences, such as loss of heterozygosity (LOH) of deleterious mutations. We mapped about 140 spontaneous reciprocal crossovers on the right arm of the yeast chromosome IV using single-nucleotide-polymorphism (SNP) microarrays. Our mapping and subsequent experiments demonstrate that inverted repeats of Ty retrotransposable elements are mitotic recombination hotspots. We found that the mitotic recombination maps on the two homologs were substantially different and were unrelated to meiotic recombination maps. Additionally, about 70% of the DNA lesions that result in LOH are likely generated during G1 of the cell cycle and repaired during S or G2. We also show that different genetic elements are associated with reciprocal crossover conversion tracts depending on the cell cycle timing of the initiating DSB.

Double-strand breaks (DSBs) are DNA lesions that can be fatal to a cell if left unrepaired. They can be caused by exogenous sources, such as gamma radiation, or endogenous stresses, such as high levels of transcription. Yeast cells primarily repair DSBs that are initiated outside of meiosis by mitotic recombination, which can result in physical exchanges between chromosomes, known as crossovers. We created a mitotic recombination map of one chromosome arm, representing 10% of the genome. This recombination map allows us to determine which regions of the chromosome arm are more susceptible to DNA damage than other regions. We were able to determine that most DSBs that result in detectable genomic changes were initiated prior to DNA replication and that some secondary DNA structures can be recombination hotspots. Recombination can also occur during meiosis, as a method of ensuring proper chromosome segregation. However, previously reported meiotic recombination maps have no correlation with our mitotic recombination map.

Mismatch repair balances leading and lagging strand DNA replication fidelity

The two DNA strands of the nuclear genome are replicated asymmetrically using three DNA polymerases, α, δ, and ε. Current evidence suggests that DNA polymerase ε (Pol ε) is the primary leading strand replicase, whereas Pols α and δ primarily perform lagging strand replication. The fact that these polymerases differ in fidelity and error specificity is interesting in light of the fact that the stability of the nuclear genome depends in part on the ability of mismatch repair (MMR) to correct different mismatches generated in different contexts during replication. Here we provide the first comparison, to our knowledge, of the efficiency of MMR of leading and lagging strand replication errors. We first use the strand-biased ribonucleotide incorporation propensity of a Pol ε mutator variant to confirm that Pol ε is the primary leading strand replicase in Saccharomyces cerevisiae. We then use polymerase-specific error signatures to show that MMR efficiency in vivo strongly depends on the polymerase, the mismatch composition, and the location of the mismatch. An extreme case of variation by location is a T-T mismatch that is refractory to MMR. This mismatch is flanked by an AT-rich triplet repeat sequence that, when interrupted, restores MMR to >95% efficiency. Thus this natural DNA sequence suppresses MMR, placing a nearby base pair at high risk of mutation due to leading strand replication infidelity. We find that, overall, MMR most efficiently corrects the most potentially deleterious errors (indels) and then the most common substitution mismatches. In combination with earlier studies, the results suggest that significant differences exist in the generation and repair of Pol α, δ, and ε replication errors, but in a generally complementary manner that results in high-fidelity replication of both DNA strands of the yeast nuclear genome.

The stability of complex and highly organized nuclear genomes partly depends on the ability of mismatch repair (MMR) to correct a variety of different mismatches generated as the leading and lagging strand templates are copied by three polymerases, each with different fidelity. Here we provide the first comparison, to our knowledge, of the efficiency of MMR of leading and lagging strand replication errors. We first confirm that Pol ε is the primary leading strand replicase, complementing earlier assignment of Pols α and δ as the primary lagging strand replicases. We then show that MMR efficiency in vivo strongly depends on the polymerase that generates the mismatch and on the composition and location of mismatches. In one extreme case, a flanking triplet repeat sequence eliminates MMR altogether. Overall, MMR is most efficient for mismatches generated at the highest rates and having the most deleterious potential, thereby ultimately achieving high-fidelity replication of both DNA strands.

Analysis of microsatellite mutations in buccal cells from a case-control study for lung cancer

Exposure to tobacco carcinogens is the major cause of human lung cancer, but even heavy smokers have only about a 10% life-time risk of developing lung cancer. Currently used screening processes, based largely on age and exposure status, have proven to be of limited clinical utility in predicting cancer risk. More precise methods of assessing an individual's risk of developing lung cancer are needed. Because of their sensitivity to DNA damage, microsatellites are potentially useful for the assessment of somatic mutational load in normal cells. We assessed mutational load using hypermutable microsatellites in buccal cells obtained from lung carcinoma cases and controls to test if such a measure could be used to estimate lung cancer risk. There was no significant association between smoking status and mutation frequency with any of the markers tested. No significant association between case status and mutation frequency was observed. Age was significantly related to mutation frequency in the microsatellite marker D7S1482. These observations indicate that somatic mutational load, as measured using mutation frequency of microsatellites in buccal cells, increases with increasing age but that subjects who develop lung cancer have a similar mutational load as those who remain cancer free. This finding suggests that mutation frequency of microsatellite mutations in buccal cells may not be a promising biomarker for lung cancer risk.

Guanine repeat-containing sequences confer transcription-dependent instability in an orientation-specific manner in yeast

Non-B DNA structures are a major contributor to the genomic instability associated with repetitive sequences. Immunoglobulin switch Mu (Sμ) region sequence is comprised of guanine-rich repeats and has high potential for forming G4 DNA, in which one strand of DNA folds into an array of guanine quartets. Taking advantage of the genetic tractability of Saccharomyces cerevisiae, we developed a recombination assay to investigate mechanisms involved in maintaining stability of G-rich repetitive sequence. By embedding Sμ sequence within recombination substrates under the control of a tetracycline-regulatable promoter, we demonstrate that the rate and orientation of transcription both affect the stability of Sμ sequence. In particular, the greatest instability was observed under high-transcription conditions when the Sμ sequence was oriented with the C-rich strand as the transcription template. The effect of transcription orientation was enhanced in the absence of the Type IB topoisomerase Top1, possibly due to enhanced R-loop formation. Loss of Sgs1 helicase and RNase H activity also increased instability, suggesting they may cooperatively function to reduce the formation of non-B DNA structures in highly transcribed regions. Finally, the Sμ sequence was unstable when transcription elongation was perturbed due to a defective THO complex. In a THO-deficient background, there was further exacerbation of orientation-dependent instability associated with the ectopically expressed, single-strand cytosine deaminase AID. The implications of our findings to understanding instability associated with potential G4 DNA forming sequences are discussed.

Damage-induced localized hypermutability

Genome instability continuously presents perils of cancer, genetic disease and death of a cell or an organism. At the same time, it provides for genome plasticity that is essential for development and evolution. We address here the genome instability confined to a small fraction of DNA adjacent to free DNA ends at uncapped telomeres and double-strand breaks. We found that budding yeast cells can tolerate nearly 20 kilobase regions of subtelomeric single-strand DNA that contain multiple UV-damaged nucleotides. During restoration to the double-strand state, multiple mutations are generated by error-prone translesion synthesis. Genome-wide sequencing demonstrated that multiple regions of damage-induced localized hypermutability can be tolerated, which leads to the simultaneous appearance of multiple mutation clusters in the genomes of UV- irradiated cells. High multiplicity and density of mutations suggest that this novel form of genome instability may play significant roles in generating new alleles for evolutionary selection as well as in the incidence of cancer and genetic disease.

The in vitro fidelity of yeast DNA polymerase δ and polymerase ε holoenzymes during dinucleotide microsatellite DNA synthesis

Elucidating the sources of genetic variation within microsatellite alleles has important implications for understanding the etiology of human diseases. Mismatch repair is a well described pathway for the suppression of microsatellite instability. However, the cellular polymerases responsible for generating microsatellite errors have not been fully described. We address this gap in knowledge by measuring the fidelity of recombinant yeast polymerase δ (Pol δ) and ɛ (Pol ɛ) holoenzymes during synthesis of a [GT/CA] microsatellite. The in vitro HSV-tk forward assay was used to measure DNA polymerase errors generated during gap-filling of complementary GT(10) and CA(10)-containing substrates and ∼90 nucleotides of HSV-tk coding sequence surrounding the microsatellites. The observed mutant frequencies within the microsatellites were 4 to 30-fold higher than the observed mutant frequencies within the coding sequence. More specifically, the rate of Pol δ and Pol ɛ misalignment-based insertion/deletion errors within the microsatellites was ∼1000-fold higher than the rate of insertion/deletion errors within the HSV-tk gene. Although the most common microsatellite error was the deletion of a single repeat unit, ∼ 20% of errors were deletions of two or more units for both polymerases. The differences in fidelity for wild type enzymes and their exonuclease-deficient derivatives were ∼2-fold for unit-based microsatellite insertion/deletion errors. Interestingly, the exonucleases preferentially removed potentially stabilizing interruption errors within the microsatellites. Since Pol δ and Pol ɛ perform not only the bulk of DNA replication in eukaryotic cells but also are implicated in performing DNA synthesis associated with repair and recombination, these results indicate that microsatellite errors may be introduced into the genome during multiple DNA metabolic pathways.

Methods to study mitotic homologous recombination and genome stability

Spontaneous mitotic recombination occurs in response to DNA damage incurred during DNA replication or from lesions that do not block replication but leave recombinogenic substrates such as single-stranded DNA gaps. Other types of damages result in general genome instability such as chromosome loss, chromosome fragmentation, and chromosome rearrangements. The genome is kept intact through recombination, repair, replication, checkpoints, and chromosome organization functions. Therefore when these pathways malfunction, genomic instabilities occur. Here we outline some general strategies to monitor a subset of the genomic instabilities: spontaneous mitotic recombination and chromosome loss, in both haploid and diploid cells. The assays, while not inclusive of all genome instability assays, give a broad assessment of general genome damage or inability to repair damage in various genetic backgrounds.

Metabolism of postsynaptic recombination intermediates

DNA double strand breaks and blocked or collapsed DNA replication forks are potentially genotoxic lesions that can result in deletions, aneuploidy or cell death. Homologous recombination (HR) is an essential process employed during repair of these forms of damage. HR allows for accurate restoration of the damaged DNA through use of a homologous template for repair. Although inroads have been made towards understanding the mechanisms of HR, ambiguity still surrounds aspects of the process. Until recently, relatively little was known concerning metabolism of postsynaptic RAD51 filaments or how synthesis dependent strand annealing intermediates are processed. This review discusses recent findings implicating RTEL1, HELQ and the Caenorhabditis elegans RAD51 paralog RFS-1 in post-strand exchange events during HR.

Post-replication repair suppresses duplication-mediated genome instability

RAD6 is known to suppress duplication-mediated gross chromosomal rearrangements (GCRs) but not single-copy sequence mediated GCRs. Here, we found that the RAD6- and RAD18-dependent post-replication repair (PRR) and the RAD5-, MMS2-, UBC13-dependent error-free PRR branch acted in concert with the replication stress checkpoint to suppress duplication-mediated GCRs formed by homologous recombination (HR). The Rad5 helicase activity, but not its RING finger, was required to prevent duplication-mediated GCRs, although the function of Rad5 remained dependent upon modification of PCNA at Lys164. The SRS2, SGS1, and HCS1 encoded helicases appeared to interact with Rad5, and epistasis analysis suggested that Srs2 and Hcs1 act upstream of Rad5. In contrast, Sgs1 likely functions downstream of Rad5, potentially by resolving DNA structures formed by Rad5. Our analysis is consistent with models in which PRR prevents replication damage from becoming double strand breaks (DSBs) and/or regulates the activity of HR on DSBs.

Genome instability is a hallmark of many cancers and underlies many inherited disorders that cause a predisposition to cancer. The human genome has many different types of duplicated sequences that can lead to genome instability by recombination-mediated pathways. We previously discovered that duplication-mediated chromosomal rearrangements are suppressed by a number of pathways. Some of these pathways were specific to rearrangements between genomic duplications. Here, we have performed a detailed analysis of pathways dependent upon RAD6, and have discovered that the error-free branch of post-replication repair (PRR) either is as an alternative to homologous recombination or prevents the generation of homologous recombination intermediates. Both of these functions could lead to genomic instability in the context of genomes containing substantial amounts of duplications. The extreme sensitivity of our assay to post-replication repair defects reveals substantial complexity in the interaction of PRR defects, suggesting the presence of many alternative PRR pathways. Together, the results emphasize the importance for appropriately balancing different repair pathways to maintain global genomic stability and highlight a number of defects that could underlie genome instabilities in some cancers.

Measurements of spontaneous rates of mutations in the recent past and the near future

The rate of spontaneous mutation in natural populations is a fundamental parameter for many evolutionary phenomena. Because the rate of mutation is generally low, most of what is currently known about mutation has been obtained through indirect, complex and imprecise methodological approaches. However, in the past few years genome-wide sequencing of closely related individuals has made it possible to estimate the rates of mutation directly at the level of the DNA, avoiding most of the problems associated with using indirect methods. Here, we review the methods used in the past with an emphasis on next generation sequencing, which may soon make the accurate measurement of spontaneous mutation rates a matter of routine.

Leaping forks at inverted repeats

Genome rearrangements are often associated with genome instability observed in cancer and other pathological disorders. Different types of repeat elements are common in genomes and are prone to instability. S-phase checkpoints, recombination, and telomere maintenance pathways have been implicated in suppressing chromosome rearrangements, but little is known about the molecular mechanisms and the chromosome intermediates generating such genome-wide instability. In the December 15, 2009, issue of Genes & Development, two studies by Paek and colleagues (2861-2875) and Mizuno and colleagues (pp. 2876-2886), demonstrate that nearby inverted repeats in budding and fission yeasts recombine spontaneously and frequently to form dicentric and acentric chromosomes. The recombination mechanism underlying this phenomenon does not appear to require double-strand break formation, and is likely caused by a replication mechanism involving template switching.

The pol3-t hyperrecombination phenotype and DNA damage-induced recombination in Saccharomyces cerevisiae is RAD50 dependent

The DNA polymerase delta (POL3/CDC2) allele pol3-t of Saccharomyces cerevisiae has previously been shown to be sensitive to methylmethanesulfonate (MMS) and has been proposed to be involved in base excision repair. Our results, however, show that the pol3-t mutation is synergistic for MMS sensitivity with MAG1, a known base excision repair gene, but it is epistatic with rad50Delta, suggesting that POL3 may be involved not only in base excision repair but also in a RAD50 dependent function. We further studied the interaction of pol3-t with rad50Delta by examining their effect on spontaneous, MMS-, UV-, and ionizing radiation-induced intrachromosomal recombination. We found that rad50Delta completely abolishes the elevated spontaneous frequency of intrachromosomal recombination in the pol3-t mutant and significantly decreases UV- and MMS-induced recombination in both POL3 and pol3-t strains. Interestingly, rad50Delta had no effect on gamma-ray-induced recombination in both backgrounds between 0 and 50 Gy. Finally, the deletion of RAD50 had no effect on the elevated frequency of homologous integration conferred by the pol3-t mutation. RAD50 is possibly involved in resolution of replication forks that are stalled by mutagen-induced external DNA damage, or internal DNA damage produced by growing the pol3-t mutant at the restrictive temperature.

Specific pathways prevent duplication-mediated genome rearrangements

We have investigated the ability of different regions of the left arm of Saccharomyces cerevisiae chromosome V to participate in the formation of gross chromosomal rearrangements (GCRs). We found that the 4.2 kb HXT13 DSF1 region sharing divergent homology with chromosomes IV, X, and XIV, similar to mammalian segmental duplications, was “at-risk” for participating in duplication-mediated GCRs generated by homologous recombination. Numerous genes and pathways, including SGS1, TOP3, RMI1, SRS2, RAD6, SLX1, SLX4, SLX5, MSH2, MSH6, RAD10 and the DNA replication stress checkpoint requiring MRC1 and TOF1 were highly specific for suppressing these GCRs compared to GCRs mediated by single copy sequences. These results indicate that the mechanisms for formation and suppression of rearrangements occurring in regions containing “at risk” sequences differ from those occurring in regions of single copy sequence. This explains how extensive genome instability is prevented in eukaryotic cells whose genomes contain numerous divergent repeated sequences.

Size-dependent palindrome-induced intrachromosomal recombination in yeast

Palindromic and quasi-palindromic sequences are important DNA motifs found in various cis-acting genetic elements, but are also known to provoke different types of genetic alterations. The instability of such motifs is clearly size-related and depends on their potential to adopt secondary structures known as hairpins and cruciforms. Here we studied the influence of palindrome size on recombination between two directly repeated copies of the yeast CYC1 gene leading to the loss of the intervening sequence ("pop-out" recombination). We show that palindromes inserted either within one copy or between the two copies of the CYC1 gene become recombinogenic only when they attain a certain critical size and we estimate this critical size to be about 70 bp. With the longest palindrome used in this study (150 bp) we observed a more than 20-fold increase in the pop-out recombination. In the sae2/com1 mutant the palindrome-stimulated recombination was completely abolished. Suppression of palindrome recombinogenicity may be crucial for the maintenance of genetic stability in organisms containing a significant number of large palindromes in their genomes, like humans.

Mutagenic and recombinagenic responses to defective DNA polymerase delta are facilitated by the Rev1 protein in pol3-t mutants of Saccharomyces cerevisiae

Defective DNA replication can result in substantial increases in the level of genome instability. In the yeast Saccharomyces cerevisiae, the pol3-t allele confers a defect in the catalytic subunit of replicative DNA polymerase delta that results in increased rates of mutagenesis, recombination, and chromosome loss, perhaps by increasing the rate of replicative polymerase failure. The translesion polymerases Pol eta, Pol zeta, and Rev1 are part of a suite of factors in yeast that can act at sites of replicative polymerase failure. While mutants defective in the translesion polymerases alone displayed few defects, loss of Rev1 was found to suppress the increased rates of spontaneous mutation, recombination, and chromosome loss observed in pol3-t mutants. These results suggest that Rev1 may be involved in facilitating mutagenic and recombinagenic responses to the failure of Pol delta. Genome stability, therefore, may reflect a dynamic relationship between primary and auxiliary DNA polymerases.

Mutagenicity of tamoxifen DNA adducts in human endometrial cells and in silico prediction of p53 mutation hotspots

Tamoxifen elevates the risk of endometrial tumours in women and alpha-(N(2)-deoxyguanosinyl)-tamoxifen adducts are reportedly present in endometrial tissue of patients undergoing therapy. Given the widespread use of tamoxifen there is considerable interest in elucidating the mechanisms underlying treatment-associated cancer. Using a combined experimental and multivariate statistical approach we have examined the mutagenicity and potential consequences of adduct formation by reactive intermediates in target uterine cells. pSP189 plasmid containing the supF gene was incubated with alpha-acetoxytamoxifen or 4-hydroxytamoxifen quinone methide (4-OHtamQM) to generate dG-N(2)-tamoxifen and dG-N(2)-4-hydroxytamoxifen, respectively. Plasmids were replicated in Ishikawa cells then screened in Escherichia coli. Treatment with both alpha-acetoxytamoxifen and 4-OHtamQM caused a dose-related increase in adduct levels, resulting in a damage-dependent increase in mutation frequency for alpha-acetoxytamoxifen; 4-OHtamQM had no apparent effect. Only alpha-acetoxytamoxifen generated statistically different supF mutation spectra relative to the spontaneous pattern, with most mutations being GC-->TA transversions. Application of the LwPy53 algorithm to the alpha-acetoxytamoxifen spectrum predicted strong GC-->TA hotspots at codons 244 and 273. These signature alterations do not correlate with current reports of the mutations observed in endometrial carcinomas from treated women, suggesting that dG-N(2)-tam adduct formation in the p53 gene is not a prerequisite for endometrial cancer initiation in women.

Chronic oxidative DNA damage due to DNA repair defects causes chromosomal instability in Saccharomyces cerevisiae

Oxidative DNA damage is likely to be involved in the etiology of cancer and is thought to accelerate tumorigenesis via increased mutation rates. However, the majority of malignant cells acquire a specific type of genomic instability characterized by large-scale genomic rearrangements, referred to as chromosomal instability (CIN). The molecular mechanisms underlying CIN are not entirely understood. We utilized Saccharomyces cerevisiae as a model system to delineate the relationship between genotoxic stress and CIN. It was found that elevated levels of chronic, unrepaired oxidative DNA damage caused chromosomal aberrations at remarkably high frequencies under both selective and nonselective growth conditions. In this system, exceeding the cellular capacity to appropriately manage oxidative DNA damage resulted in a "gain-of-CIN" phenotype and led to profound karyotypic instability. These results illustrate a novel mechanism for genome destabilization that is likely to be relevant to human carcinogenesis.

Microsatellite mutations in buccal cells are associated with aging and head and neck carcinoma

Carcinogen exposure from tobacco smoking is the major cause of upper aerodigestive tract cancer, yet heavy smokers only have about a 10% life-time risk of developing one of these cancers. Current technologies allow only limited prediction of cancer risk and there are no approved screening methods applicable to the general population. We developed a method to assess somatic mutational load using small-pool PCR (SP-PCR) and analysed mutations in DNA isolated from cells obtained by mouth rinse. Mutation levels in the hypermutable tetranucleotide marker D7S1482 were analysed in specimens from 25 head and neck squamous carcinoma (HNSCC) cases and 31 controls and tested for associations with age, smoking history and cancer status. We found a significant association between mutation frequency and age (P=0.021, Generalized Linear Model (GLM), N=56), but no influence of smoking history. Cases had higher mutation frequencies than controls when corrected for the effects of age, a difference that was statistically significant in the subgroup of 10 HNSCC patients who were treated with surgery only (P=0.017, GLM, N=41). We also present evidence that cancer status is linked to levels of nonunique, and presumably clonally derived, mutations in D7S1482. Insertion mutations were observed in 833 (79%) of 1058 alleles, of which 457 (43%) could be explained by insertion of a single repeat unit; deletion mutations were found in 225 (21%) of tested alleles. In conclusion, we demonstrate that the sensitive detection of single molecule mutations in clinical specimens is feasible by SP-PCR. Our study confirms an earlier report that microsatellite mutations increase with age and is the first to provide evidence that these mutations may be associated with cancer status in individual subjects.

ATR protects the genome against CGG.CCG-repeat expansion in Fragile X premutation mice

Fragile X mental retardation syndrome is a repeat expansion disease caused by expansion of a CGG.CCG-repeat tract in the 5' UTR of the FMR1 gene. In humans, small expansions occur more frequently on paternal transmission while large expansions are exclusively maternal in origin. It has been suggested that expansion is the result of aberrant DNA replication, repair or recombination. To distinguish amongst these possibilities we crossed mice containing 120 CGG.CCG-repeats in the 5' UTR of the mouse Fmr1 gene to mice with mutations in ATR, a protein important in the cellular response to stalled replication forks and bulky DNA lesions. We show here that ATR heterozygosity results in increased expansion rates of maternally, but not paternally, transmitted alleles. In addition, age-related somatic expansions occurred in mice of both genders that were not seen in ATR wild-type animals. Some ATR-sensitive expansion occurs in postmitotic cells including haploid gametes suggesting that aberrant DNA repair is responsible. Our data suggest that two mechanisms of repeat expansion exist that may explain the small and large expansions seen in humans. In addition, our data provide an explanation for the maternal bias of large expansions in humans and the lower incidence of these expansions in mice.

Measuring the rate of gross chromosomal rearrangements in Saccharomyces cerevisiae: A practical approach to study genomic rearrangements observed in cancer

Gross chromosomal rearrangements (GCRs), including translocations, deletions, amplifications and aneuploidy are frequently observed in various types of human cancers. Despite their clear importance in carcinogenesis, the molecular mechanisms by which GCRs are generated and held in check are poorly understood. By using a GCR assay, which can measure the rate of accumulation of spontaneous GCRs in Saccharomyces cerevisiae, we have found that many proteins involved in DNA replication, DNA repair, DNA recombination, checkpoints, chromosome remodeling, and telomere maintenance, play crucial roles in GCR metabolism. We describe here the theoretical background and practical procedures of this GCR assay. We will explain the breakpoint structure and DNA damage that lead to GCR formation. We will also summarize the pathways that suppress and enhance GCR formation. Finally, we will briefly describe similar assays developed by others and discuss their potential in studying GCR metabolism.

Unusually long palindromes are abundant in mitochondrial control regions of insects and nematodes

Background

Palindromes are known to be involved in a variety of biological processes. In the present investigation we carried out a comprehensive analysis of palindromes in the mitochondrial control regions (CRs) of several animal groups to study their frequency, distribution and architecture to gain insights into the origin of replication of mtDNA.

Methodology/Principal Findings

Many species of Arthropoda, Nematoda, Mollusca and Annelida harbor palindromes and inverted repeats (IRs) in their CRs. Lower animals like cnidarians and higher animal groups like chordates are almost devoid of palindromes and IRs. The study revealed that palindrome occurrence is positively correlated with the AT content of CRs, and that IRs are likely to give rise to longer palindromes.

Conclusions/Significance

The present study attempts to explain possible reasons and gives in silico evidence for absence of palindromes and IRs from CR of vertebrate mtDNA and acquisition and retention of the same in insects. Study of CRs of different animal phyla uncovered unique architecture of this locus, be it high abundance of long palindromes and IRs in CRs of Insecta and Nematoda, or short IRs of 10–20 nucleotides with a spacer region of 12–14 bases in subphylum Chelicerata, or nearly complete of absence of any long palindromes and IRs in Vertebrata, Cnidaria and Echinodermata.

Genetic exchange between homeologous sequences in mammalian chromosomes is averted by local homology requirements for initiation and resolution of recombination

We examined the mechanism by which recombination between imperfectly matched sequences (homeologous recombination) is suppressed in mammalian chromosomes. DNA substrates were constructed, each containing a thymidine kinase (tk) gene disrupted by insertion of an XhoI linker and referred to as a "recipient" gene. Each substrate also contained one of several "donor" tk sequences that could potentially correct the recipient gene via recombination. Each donor sequence either was perfectly homologous to the recipient gene or contained homeologous sequence sharing only 80% identity with the recipient gene. Mouse Ltk(-) fibroblasts were stably transfected with the various substrates and tk(+) segregants produced via intrachromosomal recombination were recovered. We observed exclusion of homeologous sequence from gene conversion tracts when homeologous sequence was positioned adjacent to homologous sequence in the donor but not when homeologous sequence was surrounded by homology in the donor. Our results support a model in which homeologous recombination in mammalian chromosomes is suppressed by a nondestructive dismantling of mismatched heteroduplex DNA (hDNA) intermediates. We suggest that mammalian cells do not dismantle mismatched hDNA by responding to mismatches in hDNA per se but rather rejection of mismatched hDNA appears to be driven by a requirement for localized homology for resolution of recombination.

Tetraplex DNA transitions within the human c-myc promoter detected by multivariate curve resolution of fluorescence resonance energy transfer

The nuclease hypersensitive element (NHE) III(I) of the c-myc promoter regulates the expression of oncogene c-myc and hence is an important anti-cancer target. Paranemic secondary structure formation within the promoter has been implicated in mechanistic regulation models. Here, it is shown that two monomeric tetraplexes form within the c-myc promoter, which coexist in solution. The development and application of a new experimental approach for detection of conformation transitions in nucleic acids [which exploits the sensitivity of fluorescence resonance energy transfer (FRET) for theoretical spectral resolution by multivariate curve resolution-alternating least-squares (MCR-ALS) method] has been used for this study. The pK(a) for tetraplex transitions are centered around 5.9 +/- 0.2 (between two intercalation topologies) and 6.8 +/- 0.1 (tetraplex to random coil). The presence of two tetraplexes has been further confirmed by S1 nuclease digestion. Finally, it is established that MCR-ALS analysis of FRET at different temperatures, pH, and salt concentrations allows resolution of pure species. Results are discussed in the light of recent observations implicating paranemic DNA motifs within the c-myc NHE in regulation of the oncogene. This method has several advantages over other methods vis-à-vis, high sensitivity and linear detection over a wide concentration range and, particularly, potential applications in intracellular probing.

Y chromosome instability in testicular cancer

Approximately 15-25% of male infertility cases carry extensive azoospermic factor (AZF) deletions. Moreover, about 80% of Finnish testicular germ cell tumors (TGCT) and about 23-25% of TGCTs from other geographic regions carry short and interstitial AZF deletions. In infertility cases the AZF deficiency occurs in the germ cells of the proband father giving rise to mosaic sperm populations comprising non-deleted and deleted sperms. Fertilization of an oocyte by a Y deleted sperm will give rise to an AZF-deleted and infertile F1 male. In TGCTs the AZF deletions take place in the initial stages of embryogenesis producing individuals that are a mosaic of Y deleted and non-deleted cell lineages. Carcinoma in situ (CIS) is a premalignant lesion that some believe may develop in gonads of male embryos before the ninth week of age due to transformation of a totipotent primordial germ cell. If the transformed cell carries AZF deletions the resultant CIS will also have Y deletions. CIS will differentiate into seminoma or into embryonal carcinoma and non-seminomas in about 1 x 10(-3) of the young adults carrying premalignant CIS outgrowths; if the CIS lesion has AZF deletions the derived forms of testicular cancer will also exhibit these deletions. AZF deletions play no role in the development of testicular cancers. On the other hand, they are a marker of Y chromosome instability and eventually of a more generalized pattern of genome instability associated with the appearance of TGCT. Genetic factors such as malfunction of metabolizing genes, DNA repairing genes, Y-linked or X-linked genes have been considered as possible causes of AZF deletions in testicular cancer. Yet, the exact identification of the genes involved remains elusive. AZF deletions have also been identified in non-Hodgkin lymphomas and in colorectal cancers, two forms of malignancy that have been found to be associated with TGCTs.

Mutagenesis and the three R's in yeast

Mutagenesis is a prerequisite for evolution and also is an important contributor to human diseases. Most mutations in actively dividing cells originate during DNA replication as errors introduced when copying an undamaged DNA template or during the bypass of DNA lesions. In addition, mutations can be introduced during the repair of DNA double-strand breaks by either homologous recombination or non-homologous end-joining pathways. Finally, although generally considered to be a very high-fidelity process, the excision repair of DNA damage may be an important contributor to mutagenesis in non-dividing cells. In this review, we will discuss the well-known contributions of DNA replication to mutagenesis in Saccharomyces cerevisiae, as well as the less-appreciated contributions of recombination and repair to mutagenesis in this organism.

Palindrome content of the yeast Saccharomyces cerevisiae genome

Palindromic sequences are important DNA motifs involved in the regulation of different cellular processes, but are also a potential source of genetic instability. In order to initiate a systematic study of palindromes at the whole genome level, we developed a computer program that can identify, locate and count palindromes in a given sequence in a strictly defined way. All palindromes, defined as identical inverted repeats without spacer DNA, can be analyzed and sorted according to their size, frequency, GC content or alphabetically. This program was then used to prepare a catalog of all palindromes present in the chromosomal DNA of the yeast Saccharomyces cerevisiae. For each palindrome size, the observed palindrome counts were significantly different from those in the randomly generated equivalents of the yeast genome. However, while the short palindromes (2-12 bp) were under-represented, the palindromes longer than 12 bp were over-represented, AT-rich and preferentially located in the intergenic regions. The 44-bp palindrome found between the genes CDC53 and LYS21 on chromosome IV was the longest palindrome identified and contained only two C-G base pairs. Avoidance of coding regions was also observed for palindromes of 4-12 bp, but was less pronounced. Dinucleotide analysis indicated a strong bias against palindromic dinucleotides that could explain the observed short palindrome avoidance. We discuss some possible mechanisms that may influence the evolutionary dynamics of palindromic sequences in the yeast genome.

When, where and how the bridge breaks: anaphase bridge breakage plays a crucial role in gene amplification and HSR generation

Amplified genes are frequently localized on extrachromosomal double minutes (DMs) or in chromosomal homogeneously staining regions (HSRs). We previously showed that a plasmid bearing a mammalian replication initiation region could efficiently generate DMs and HSRs after transfection into human tumor cell lines. The Breakage-Fusion-Bridge (BFB) cycle model, a classical model that explains how HSRs form, could also be used to explain how the transfected plasmids generate HSRs. The BFB cycle model involves anaphase bridge formation due to the presence of dicentric chromosomes, followed by the breakage of the bridge. In this study, we used our plasmid-based model system to analyze how anaphase bridges break during mitosis. Dual-color fluorescence in situ hybridization analyses revealed that anaphase bridges were most frequently severed in their middle irrespective of their lengths, which suggests that a structurally fragile site exists in the middle of the anaphase bridge. Breakage of the chromosomal bridges occurred prior to nuclear membrane reformation and the completion of cytokinesis, which indicates that mechanical tension rather than cytokinesis is primarily responsible for severing anaphase bridges. Time-lapse observation of living cells revealed that the bridges rapidly shrink after being severed. If HSR length was extended too far, the bridge could no longer be resolved and became tangled depending on the tension. The unbroken bridge appeared to inhibit the completion of cytokinesis. These observations strongly suggest that anaphase bridges are highly elastic and that the length of the spindle axis determines the maximal HSR length.

Abundance, distribution, and mutation rates of homopolymeric nucleotide runs in the genome of Caenorhabditis elegans

Homopolymeric nucleotide runs, also called mononucleotide microsatellites, are a ubiquitous, dominant, and mutagenic feature of eukaryotic genomes. A clear understanding of the forces that shape patterns of homopolymer evolution, however, is lacking. We provide a focused investigation of the abundance, chromosomal distribution, and mutation spectra of the four strand-specific homopolymer types (A, T, G, C) >or=8 bp in the genome of Caenorhabditis elegans. A and T homopolymers vastly outnumber G and C HPs, and the run-length distributions of A and T homopolymers differ significantly from G and C homopolymers. A scanning window analysis of homopolymer chromosomal distribution reveals distinct clusters of homopolymer density in autosome arms that are regions of high recombination in C. elegans. Dramatic biases are detected among closely spaced homopolymers; for instance, we observe 994 A homopolymers immediately followed by a T homopolymer (5' to 3') and only 8 instances of T homopolymers directly followed by an A homopolymer. Empirical homopolymer mutation assays in a set of C. elegans mutation-accumulation lines reveal an approximately 20-fold higher mutation rate for G and C homopolymers compared to A and T homopolymers. Nuclear A and T homopolymers are also found to mutate approximately 100-fold more slowly than mitochondrial A and T homopolymers. This integrative approach yields a total nuclear genome-wide homopolymer mutation rate estimate of approximately 1.6 mutations per genome per generation.

Cis-acting transmission of genomic instability

Genomic instability is a highly pleiotropic phenotype, which may reflect a variety of underlying mechanisms. Destabilization has been shown in some cases to involve mutational alteration or inactivation of trans-acting cellular factors, for example, p53 or mismatch repair functions. However, aspects of instability are not well explained by mutational inactivation of trans-acting factors, and other epigenetic and cis-acting mechanisms have recently been proposed. The trans and cis models result in divergent predictions for the distribution of instability-associated genetic alterations within the genome, and for the inheritance of genomic instability among sibling sub-clones of unstable parents. These predictions have been tested in this study primarily by tracking the karyotypic distribution of chromosomal rearrangements in clones and sub-clones exhibiting radiation-induced genomic instability; inheritance of mutator phenotypes was also analyzed. The results indicate that genomic instability is unevenly transmitted to sibling sub-clones, that chromosomal rearrangements within unstable clones are non-randomly distributed throughout the karyotype, and that the majority of chromosomal rearrangements associated with instability affect trisomic chromosomal segments. Observations of instability in trisomic regions suggests that in addition to promoting further alterations in chromosomal number, aneuploidy can affect the recovery of structural rearrangements. In summary, these findings cannot be fully explained by invoking a homogeneously distributed factor acting in trans, but do provide support for previous suggestions that genomic instability may in part be driven by a cis-acting mechanism.

High frequency of mosaicism among patients with neurofibromatosis type 1 (NF1) with microdeletions caused by somatic recombination of the JJAZ1 gene

Detailed analyses of 20 patients with sporadic neurofibromatosis type 1 (NF1) microdeletions revealed an unexpected high frequency of somatic mosaicism (8/20 [40%]). This proportion of mosaic deletions is much higher than previously anticipated. Of these deletions, 16 were identified by a screen of unselected patients with NF1. None of the eight patients with mosaic deletions exhibited the mental retardation and facial dysmorphism usually associated with NF1 microdeletions. Our study demonstrates the importance of a general screening for NF1 deletions, regardless of a special phenotype, because of a high estimated number of otherwise undetected mosaic NF1 microdeletions. In patients with mosaicism, the proportion of cells with the deletion was 91%-100% in peripheral leukocytes but was much lower (51%-80%) in buccal smears or peripheral skin fibroblasts. Therefore, the analysis of other tissues than blood is recommended, to exclude mosaicism with normal cells in patients with NF1 microdeletions. Furthermore, our study reveals breakpoint heterogeneity. The classic 1.4-Mb deletion was found in 13 patients. These type I deletions encompass 14 genes and have breakpoints in the NF1 low-copy repeats. However, we identified a second major type of NF1 microdeletion, which spans 1.2 Mb and affects 13 genes. This type II deletion was found in 8 (38%) of 21 patients and is mediated by recombination between the JJAZ1 gene and its pseudogene. The JJAZ1 gene, which is completely deleted in patients with type I NF1 microdeletions and is disrupted in deletions of type II, is highly expressed in brain structures associated with learning and memory. Thus, its haploinsufficiency might contribute to mental impairment in patients with constitutional NF1 microdeletions. Conspicuously, seven of the eight mosaic deletions are of type II, whereas only one was a classic type I deletion. Therefore, the JJAZ1 gene is a preferred target of strand exchange during mitotic nonallelic homologous recombination. Although type I NF1 microdeletions occur by interchromosomal recombination during meiosis, our findings imply that type II deletions are mediated by intrachromosomal recombination during mitosis. Thus, NF1 microdeletions acquired during mitotic cell divisions differ from those occurring in meiosis and are caused by different mechanisms.

Twinkle and POLG defects enhance age-dependent accumulation of mutations in the control region of mtDNA

Autosomal dominant and/or recessive progressive external ophthalmoplegia (ad/arPEO) is associated with mtDNA mutagenesis. It can be caused by mutations in three nuclear genes, encoding the adenine nucleotide translocator 1, the mitochondrial helicase Twinkle or DNA polymerase gamma (POLG). How mutations in these genes result in progressive accumulation of multiple mtDNA deletions in post- mitotic tissues is still unclear. A recent hypothesis suggested that mtDNA replication infidelity could promote slipped mispairing, thereby stimulating deletion formation. This hypothesis predicts that mtDNA of ad/arPEO patients will contain frequent mutations throughout; in fact, our analysis of muscle from ad/arPEO patients revealed an age-dependent, enhanced accumulation of point mutations in addition to deletions, but specifically in the mtDNA control region. Both deleted and non-deleted mtDNA molecules showed increased point mutation levels, as did mtDNAs of patients with a single mtDNA deletion, suggesting that point mutations do not cause multiple deletions. Deletion breakpoint analysis showed frequent breakpoints around homopolymeric runs, which could be a signature of replication stalling. Therefore, we propose replication stalling as the principal cause of deletion formation.

Context of deletions and insertions in human coding sequences

We studied the dependence of the rate of short deletions and insertions on their contexts using the data on mutations within coding exons at 19 human loci that cause mendelian diseases. We confirm that periodic sequences consisting of three to five or more nucleotides are mutagenic. Mutability of sequences with strongly biased nucleotide composition is also elevated, even when mutations within homonucleotide runs longer than three nucleotides are ignored. In contrast, no elevated mutation rates have been detected for imperfect direct or inverted repeats. Among known candidate contexts, the indel context GTAAGT and regions with purine-pyrimidine imbalance between the two DNA strands are mutagenic in our sample, and many others are not mutagenic. Data on mutation hot spots suggest two novel contexts that increase the deletion rate. Comprehensive analysis of mutability of all possible contexts of lengths four, six, and eight indicates a substantially elevated deletion rate within YYYTG and similar sequences, which is one of the two contexts revealed by the hot spots. Possible contexts that increase the insertion rate (AT(A/C)(A/C)GCC and TACCRC) and decrease deletion (TATCGC) or insertion (GCGG) rates have also been identified. Two-thirds of deletions remove a repeat, and over 80% of insertions create a repeat, i.e., they are duplications.

Reciprocal crossovers and a positional preference for strand exchange in recombination events resulting in deletion or duplication of chromosome 17p11.2

Smith-Magenis syndrome (SMS) is caused by an approximately 4-Mb heterozygous interstitial deletion on chromosome 17p11.2 in approximately 80%-90% of affected patients. Three large ( approximately 200 kb), complex, and highly homologous ( approximately 98%) low-copy repeats (LCRs) are located inside or flanking the SMS common deletion. These repeats, also known as "SMS-REPs," are termed "distal," "middle," and "proximal." The directly oriented distal and proximal copies act as substrates for nonallelic homologous recombination resulting in both the deletion associated with SMS and the reciprocal duplication: dup(17)(p11.2p11.2). Using restriction enzyme cis-morphism analyses and direct sequencing, we mapped the regions of strand exchange in 16 somatic-cell hybrids that harbor only the recombinant SMS-REP. Our studies showed that the sites of crossovers were distributed throughout the region of homology between the distal and proximal SMS-REPs. However, despite approximately 170 kb of high homology, 50% of the recombinant junctions occurred in a 12.0-kb region within the KER gene clusters. DNA sequencing of this hotspot (positional preference for strand exchange) in seven recombinant SMS-REPs narrowed the crossovers to an approximately 8-kb interval. Four of them occurred in a 1,655-bp region rich in polymorphic nucleotides that could potentially reflect frequent gene conversion. For further evaluation of the strand exchange frequency in patients with SMS, novel junction fragments from the recombinant SMS-REPs were identified. As predicted by the reciprocal-recombination model, junction fragments were also identified from this hotspot region in patients with dup(17)(p11.2p11.2), documenting reciprocity of the positional preference for strand exchange. Several potential cis-acting recombination-promoting sequences were identified within the hotspot. It is interesting that we found 2.1-kb AT-rich inverted repeats flanking the proximal and middle KER gene clusters but not the distal one. The role of any or all of these in stimulating double-strand breaks around this positional recombination hotspot remains to be explored.

Molecular characterization of Brucella abortus chromosome II recombination

Large-scale genomic rearrangements including inversions, deletions, and duplications are significant in bacterial evolution. The recently completed Brucella melitensis 16M and Brucella suis 1330 genomes have facilitated the investigation of such events in the Brucella spp. Suppressive subtractive hybridization (SSH) was employed in identifying genomic differences between B. melitensis 16M and Brucella abortus 2308. Analysis of 45 SSH clones revealed several deletions on chromosomes of B. abortus and B. melitensis that encoded proteins of various metabolic pathways. A 640-kb inversion on chromosome II of B. abortus has been reported previously (S. Michaux Charachon, G. Bourg, E. Jumas Bilak, P. Guigue Talet, A. Allardet Servent, D. O'Callaghan, and M. Ramuz, J. Bacteriol. 179:3244-3249, 1997) and is further described in this study. One end of the inverted region is located on a deleted TATGC site between open reading frames BMEII0292 and BMEII0293. The other end inserted at a GTGTC site of the cyclic-di-GMP phosphodiesterase A (PDEA) gene (BMEII1009), dividing PDEA into two unequal DNA segments of 160 and 977 bp. As a consequence of inversion, the 160-bp segment that encodes the N-terminal region of PDEA was relocated at the opposite end of the inverted chromosomal region. The splitting of the PDEA gene most likely inactivated the function of this enzyme. A recombination mechanism responsible for this inversion is proposed.

Trinucleotide repeat instability: a hairpin curve at the crossroads of replication, recombination, and repair

The trinucleotide repeats that expand to cause human disease form hairpin structures in vitro that are proposed to be the major source of their genetic instability in vivo. If a replication fork is a train speeding along a track of double-stranded DNA, the trinucleotide repeats are a hairpin curve in the track. Experiments have demonstrated that the train can become derailed at the hairpin curve, resulting in significant damage to the track. Repair of the track often results in contractions and expansions of track length. In this review we introduce the in vitro evidence for why CTG/CAG and CCG/CGG repeats are inherently unstable and discuss how experiments in model organisms have implicated the replication, recombination and repair machinery as contributors to trinucleotide repeat instability in vivo.

Theoretical analysis of mutation hotspots and their DNA sequence context specificity

Mutation frequencies vary significantly along nucleotide sequences such that mutations often concentrate at certain positions called hotspots. Mutation hotspots in DNA reflect intrinsic properties of the mutation process, such as sequence specificity, that manifests itself at the level of interaction between mutagens, DNA, and the action of the repair and replication machineries. The hotspots might also reflect structural and functional features of the respective DNA sequences. When mutations in a gene are identified using a particular experimental system, resulting hotspots could reflect the properties of the gene product and the mutant selection scheme. Analysis of the nucleotide sequence context of hotspots can provide information on the molecular mechanisms of mutagenesis. However, the determinants of mutation frequency and specificity are complex, and there are many analytical methods for their study. Here we review computational approaches for analyzing mutation spectra (distribution of mutations along the target genes) that include many mutable (detectable) positions. The following methods are reviewed: derivation of a consensus sequence, application of regression approaches to correlate nucleotide sequence features with mutation frequency, mutation hotspot prediction, analysis of oligonucleotide composition of regions containing mutations, pairwise comparison of mutation spectra, analysis of multiple spectra, and analysis of "context-free" characteristics. The advantages and pitfalls of these methods are discussed and illustrated by examples from the literature. The most reliable analyses were obtained when several methods were combined and information from theoretical analysis and experimental observations was considered simultaneously. Simple, robust approaches should be used with small samples of mutations, whereas combinations of simple and complex approaches may be required for large samples. We discuss several well-documented studies where analysis of mutation spectra has substantially contributed to the current understanding of molecular mechanisms of mutagenesis. The nucleotide sequence context of mutational hotspots is a fingerprint of interactions between DNA and DNA repair, replication, and modification enzymes, and the analysis of hotspot context provides evidence of such interactions.

Double-stranded RNA-mediated gene silencing in fission yeast

Double-stranded RNA (dsRNA) can specifically inhibit gene expression in a variety of organisms by invoking post-transcriptional degradation of homologous mRNA. Here we show that dsRNA-mediated gene regulation also occurs in the fission yeast Schizosaccharomyces pombe. We present evidence that: (i) reporter gene silencing is significantly enhanced when additional non-coding sense RNA is co-expressed with antisense RNA; (ii) expression of a panhandle RNA also silences target gene expression; (iii) expression of dsRNA is associated with siRNAs; (iv) a novel host-encoded factor which enhances antisense RNA gene silencing also enhances panhandle RNA-mediated gene inhibition. Both the exogenously introduced lacZ and c-myc genes are shown to be susceptible to dsRNA- mediated gene silencing in this model. Taken together, these data indicate that RNA-mediated gene silencing can occur through a RNAi-like mechanism in fission yeast.

Cadmium is a mutagen that acts by inhibiting mismatch repair

Most errors that arise during DNA replication can be corrected by DNA polymerase proofreading or by post-replication mismatch repair (MMR). Inactivation of both mutation-avoidance systems results in extremely high mutability that can lead to error catastrophe. High mutability and the likelihood of cancer can be caused by mutations and epigenetic changes that reduce MMR. Hypermutability can also be caused by external factors that directly inhibit MMR. Identifying such factors has important implications for understanding the role of the environment in genome stability. We found that chronic exposure of yeast to environmentally relevant concentrations of cadmium, a known human carcinogen, can result in extreme hypermutability. The mutation specificity along with responses in proofreading-deficient and MMR-deficient mutants indicate that cadmium reduces the capacity for MMR of small misalignments and base-base mismatches. In extracts of human cells, cadmium inhibited at least one step leading to mismatch removal. Together, our data show that a high level of genetic instability can result from environmental impediment of a mutation-avoidance system.

Characterization of the hyperrecombination phenotype of the pol3-t mutation of Saccharomyces cerevisiae

The DNA polymerase delta (Pol3p/Cdc2p) allele pol3-t of Saccharomyces cerevisiae has previously been shown to increase the frequency of deletions between short repeats (several base pairs), between homologous DNA sequences separated by long inverted repeats, and between distant short repeats, increasing the frequency of genomic deletions. We found that the pol3-t mutation increased intrachromosomal recombination events between direct DNA repeats up to 36-fold and interchromosomal recombination 14-fold. The hyperrecombination phenotype of pol3-t was partially dependent on the Rad52p function but much more so on Rad1p. However, in the double-mutant rad1 Delta rad52 Delta, the pol3-t mutation still increased spontaneous intrachromosomal recombination frequencies, suggesting that a Rad1p Rad52p-independent single-strand annealing pathway is involved. UV and gamma-rays were less potent inducers of recombination in the pol3-t mutant, indicating that Pol3p is partly involved in DNA-damage-induced recombination. In contrast, while UV- and gamma-ray-induced intrachromosomal recombination was almost completely abolished in the rad52 or the rad1 rad52 mutant, there was still good induction in those mutants in the pol3-t background, indicating channeling of lesions into the above-mentioned Rad1p Rad52p-independent pathway. Finally, a heterozygous pol3-t/POL3 mutant also showed an increased frequency of deletions and MMS sensitivity at the restrictive temperature, indicating that even a heterozygous polymerase delta mutation might increase the frequency of genetic instability.

Molecular characteristics of spontaneous deletions in the hyperthermophilic archaeon Sulfolobus acidocaldarius

Prokaryotic genomes acquire and eliminate blocks of DNA sequence by lateral gene transfer and spontaneous deletion, respectively. The basic parameters of spontaneous deletion, which are expected to influence the course of genome evolution, have not been determined for any hyperthermophilic archaeon. We therefore screened a number of independent pyrimidine auxotrophs of Sulfolobus acidocaldarius for deletions and sequenced those detected. Deletions accounted for only 0.4% of spontaneous pyrE mutations, corresponding to a frequency of about 10(-8) per cell. Nucleotide sequence analysis of five independent deletions showed no significant association of the endpoints with short direct repeats, despite the fact that several such repeats occur within the pyrE gene and that duplication mutations in pyrE reverted at high frequencies. Endpoints of the spontaneous deletions did not coincide with short inverted repeats or potential stem-loop structures. No consensus sequence common to all the deletions could be identified, although two deletions showed the potential of being stabilized by octanucleotide sequences elsewhere in pyrE, and another pair of deletions shared an octanucleotide at their 3' ends. The unusually low frequency and low sequence dependence of spontaneous deletions in the S. acidocaldarius pyrE gene compared to other genetic systems could not be explained in terms of possible constraints imposed by the 5-fluoroorotate selection.