Effect of DNA cytosine methylation upon deamination-induced mutagenesis in a natural target sequence in duplex DNA

The clock-like accumulation of germline and somatic mutations can arise from the interplay of DNA damage and repair

The rates at which mutations accumulate across human cell types vary. To identify causes of this variation, mutations are often decomposed into a combination of the single-base substitution (SBS) "signatures" observed in germline, soma, and tumors, with the idea that each signature corresponds to one or a small number of underlying mutagenic processes. Two such signatures turn out to be ubiquitous across cell types: SBS signature 1, which consists primarily of transitions at methylated CpG sites thought to be caused by spontaneous deamination, and the more diffuse SBS signature 5, which is of unknown etiology. In cancers, the number of mutations attributed to these 2 signatures accumulates linearly with age of diagnosis, and thus the signatures have been termed "clock-like." To better understand this clock-like behavior, we develop a mathematical model that includes DNA replication errors, unrepaired damage, and damage repaired incorrectly. We show that mutational signatures can exhibit clock-like behavior because cell divisions occur at a constant rate and/or because damage rates remain constant over time, and that these distinct sources can be teased apart by comparing cell lineages that divide at different rates. With this goal in mind, we analyze the rate of accumulation of mutations in multiple cell types, including soma as well as male and female germline. We find no detectable increase in SBS signature 1 mutations in neurons and only a very weak increase in mutations assigned to the female germline, but a significant increase with time in rapidly dividing cells, suggesting that SBS signature 1 is driven by rounds of DNA replication occurring at a relatively fixed rate. In contrast, SBS signature 5 increases with time in all cell types, including postmitotic ones, indicating that it accumulates independently of cell divisions; this observation points to errors in DNA repair as the key underlying mechanism. Thus, the two "clock-like" signatures observed across cell types likely have distinct origins, one set by rates of cell division, the other by damage rates.

Antimicrobial resistance in the wild: Insights from epigenetics

Spreading of bacterial and fungal strains that are resistant to antimicrobials poses a serious threat to the well-being of humans, animals, and plants. Antimicrobial resistance has been mainly investigated in clinical settings. However, throughout their evolutionary history microorganisms in the wild have encountered antimicrobial substances, forcing them to evolve strategies to combat antimicrobial action. It is well known that many of these strategies are based on genetic mechanisms, but these do not fully explain important aspects of the antimicrobial response such as the rapid development of resistance, reversible phenotypes, and hetero-resistance. Consequently, attention has turned toward epigenetic pathways that may offer additional insights into antimicrobial mechanisms. The aim of this review is to explore the epigenetic mechanisms that confer antimicrobial resistance, focusing on those that might be relevant for resistance in the wild. First, we examine the presence of antimicrobials in natural settings. Then we describe the documented epigenetic mechanisms in bacteria and fungi associated with antimicrobial resistance and discuss innovative epigenetic editing techniques to establish causality in this context. Finally, we discuss the relevance of these epigenetic mechanisms on the evolutionary dynamics of antimicrobial resistance in the wild, emphasizing the critical role of priming in the adaptation process. We underscore the necessity of incorporating non-genetic mechanisms into our understanding of antimicrobial resistance evolution. These mechanisms offer invaluable insights into the dynamics of antimicrobial adaptation within natural ecosystems.

IDH mutation, glioma immunogenicity, and therapeutic challenge of primary mismatch repair deficient IDH-mutant astrocytoma PMMRDIA: a systematic review

In 2021, Suwala et al. described Primary Mismatch Repair Deficient IDH-mutant Astrocytoma (PMMRDIA) as a distinct group of gliomas. In unsupervised clustering, PMMRDIA forms distinct cluster, separate from other IDH-mutant gliomas, including IDH-mutant gliomas with secondary mismatch repair (MMR) deficiency. In the published cohort, three patients received treatment with an immune checkpoint blocker (ICB), yet none exhibited a response, which aligns with existing knowledge about the decreased immunogenicity of IDH-mutant gliomas in comparison to IDH-wildtype. In the case of PMMRDIA, the inherent resistance to the standard-of-care temozolomide caused by MMR deficiency is an additional challenge. It is known that a gain-of-function mutation of IDH1/2 genes produces the oncometabolite R-2-hydroxyglutarate (R-2-HG), which increases DNA and histone methylation contributing to the characteristic glioma-associated CpG island methylator phenotype (G-CIMP). While other factors could be involved in remodeling the tumor microenvironment (TME) of IDH-mutant gliomas, this systematic review emphasizes the role of R-2-HG and the subsequent G-CIMP in immune suppression. This highlights a potential actionable pathway to enhance the response of ICB, which might be relevant for addressing the unmet therapeutic challenge of PMMRDIA.

Disentangling sources of clock-like mutations in germline and soma

The rates of mutations vary across cell types. To identify causes of this variation, mutations are often decomposed into a combination of the single base substitution (SBS) "signatures" observed in germline, soma and tumors, with the idea that each signature corresponds to one or a small number of underlying mutagenic processes. Two such signatures turn out to be ubiquitous across cell types: SBS signature 1, which consists primarily of transitions at methylated CpG sites caused by spontaneous deamination, and the more diffuse SBS signature 5, which is of unknown etiology. In cancers, the number of mutations attributed to these two signatures accumulates linearly with age of diagnosis, and thus the signatures have been termed "clock-like." To better understand this clock-like behavior, we develop a mathematical model that includes DNA replication errors, unrepaired damage, and damage repaired incorrectly. We show that mutational signatures can exhibit clock-like behavior because cell divisions occur at a constant rate and/or because damage rates remain constant over time, and that these distinct sources can be teased apart by comparing cell lineages that divide at different rates. With this goal in mind, we analyze the rate of accumulation of mutations in multiple cell types, including soma as well as male and female germline. We find no detectable increase in SBS signature 1 mutations in neurons and only a very weak increase in mutations assigned to the female germline, but a significant increase with time in rapidly-dividing cells, suggesting that SBS signature 1 is driven by rounds of DNA replication occurring at a relatively fixed rate. In contrast, SBS signature 5 increases with time in all cell types, including post-mitotic ones, indicating that it accumulates independently of cell divisions; this observation points to errors in DNA repair as the key underlying mechanism. Thus, the two "clock-like" signatures observed across cell types likely have distinct origins, one set by rates of cell division, the other by damage rates.

A combinatorial system to examine the enzymatic repair of multiply damaged DNA substrates

DNA damage drives genetic mutations that underlie the development of cancer in humans. Multiple pathways have been described in mammalian cells which can repair this damage. However, most work to date has focused upon single lesions in DNA. We present here a combinatorial system which allows assembly of duplexes containing single or multiple types of damage by ligating together six oligonucleotides containing damaged or modified bases. The combinatorial system has dual fluorescent labels allowing examination of both strands simultaneously, in order to study interactions or competition between different DNA repair pathways. Using this system, we demonstrate how repair of oxidative damage in one DNA strand can convert a mispaired T:G deamination intermediate into a T:A mutation. We also demonstrate that slow repair of a T:G mispair, relative to a U:G mispair, by the human methyl-binding domain 4 DNA glycosylase provides a competitive advantage to competing repair pathways, and could explain why CpG dinucleotides are hotspots for C to T mutations in human tumors. Data is also presented that suggests repair of closely spaced lesions in opposing strands can be repaired by a combination of short and long-patch base excision repair and simultaneous repair of multiply damage sites can potentially lead to lethal double strand breaks.

The influence of DNA methylation on the sequence specificity of UVB- and UVC-induced DNA damage

Ultraviolet light (UV) is one of the most common DNA damaging agents in the human environment. This paper examined the influence of DNA methylation on the level of UVB- and UVC-induced DNA damage. A purified DNA sequence containing CpG dinucleotides was methylated with a CpG methylase. We employed the linear amplification technique and the end-labelling approach followed by capillary electrophoresis with laser-induced fluorescence to investigate the sequence specificity of UV-induced DNA damage. The linear amplification technique mainly detects cyclobutane pyrimidine dimer (CPD) adducts, while the end-labelling approach mainly detects 6-4 photoproduct (6-4PP) lesions. The levels of CPD and 6-4PP adducts detected in methylated/unmethylated labelled sequences were analysed. The comparison showed that 5-methyl-cytosine significantly reduced the level of both CPD and 6-4PP adducts after UVB (308 nm) and UVC (254 nm) irradiation compared with the non-methylated counterpart.

Role of DNA Methylation in Mediating Genetic Risk of Psychiatric Disorders

Psychiatric disorders are common, complex, and heritable conditions estimated to be the leading cause of disability worldwide. The last decade of research in genomics of psychiatry, performed by multinational, and multicenter collaborative efforts on hundreds of thousands of mental disorder cases and controls, provided invaluable insight into the genetic risk variants of these conditions. With increasing cohort sizes, more risk variants are predicted to be identified in the near future, but there appears to be a knowledge gap in understanding how these variants contribute to the pathophysiology of psychiatric disorders. Majority of the identified common risk single-nucleotide polymorphisms (SNPs) are non-coding but are enriched in regulatory regions of the genome. It is therefore of great interest to study the impact of identified psychiatric disorders' risk SNPs on DNA methylation, the best studied epigenetic modification, playing a pivotal role in the regulation of transcriptomic processes, brain development, and functioning. This work outlines the mechanisms through which risk SNPs can impact DNA methylation levels and provides a summary of current evidence on the role of DNA methylation in mediating the genetic risk of psychiatric disorders.

Quantifying Influences on Intragenomic Mutation Rate

We report work to quantify the impact on the probability of human genome polymorphism both of recombination and of sequence context at different scales. We use population-based analyses of data on human genetic variants obtained from the public Ensembl database. For recombination, we calculate the variance due to recombination and the probability that a recombination event causes a mutation. We employ novel statistical procedures to take account of the spatial auto-correlation of recombination and mutation rates along the genome. Our results support the view that genomic diversity in recombination hotspots arises largely from a direct effect of recombination on mutation rather than predominantly from the effect of selective sweeps. We also use the statistic of variance due to context to compare the effect on the probability of polymorphism of contexts of various sizes. We find that when the 12 point mutations are considered separately, variance due to context increases significantly as we move from 3-mer to 5-mer and from 5-mer to 7-mer contexts. However, when all mutations are considered in aggregate, these differences are outweighed by the effect of interaction between the central base and its immediate neighbors. This interaction is itself dominated by the transition mutations, including, but not limited to, the CpG effect. We also demonstrate strand-asymmetry of contextual influence in intronic regions, which is hypothesized to be a result of transcription coupled DNA repair. We consider the extent to which the measures we have used can be used to meaningfully compare the relative magnitudes of the impact of recombination and context on mutation.

Cytosine Methylation Affects the Mutability of Neighboring Nucleotides in Germline and Soma

Methylated cytosines deaminate at higher rates than unmethylated cytosines, and the lesions they produce are repaired less efficiently. As a result, methylated cytosines are mutational hotspots. Here, combining rare polymorphism and base-resolution methylation data in humans, Arabidopsis thaliana, and rice (Oryza sativa), we present evidence that methylation state affects mutation dynamics not only at the focal cytosine but also at neighboring nucleotides. In humans, contrary to prior suggestions, we find that nucleotides in the close vicinity (±3 bp) of methylated cytosines mutate less frequently. Reduced mutability around methylated CpGs is also observed in cancer genomes, considering single nucleotide variants alongside tissue-of-origin-matched methylation data. In contrast, methylation is associated with increased neighborhood mutation risk in A. thaliana and rice. The difference in neighborhood mutation risk is less pronounced further away from the focal CpG and modulated by regional GC content. Our results are consistent with a model where altered risk at neighboring bases is linked to lesion formation at the focal CpG and subsequent long-patch repair. Our findings indicate that cytosine methylation has a broader mutational footprint than is commonly assumed.

The role of epigenetics in lysosomal storage disorders: Uncharted territory

The study of the contribution of epigenetic mechanisms, including DNA methylation, histone modifications, and microRNAs, to human disease has enhanced our understanding of different cellular processes and diseased states, as well as the effect of environmental factors on phenotypic outcomes. Epigenetic studies may be particularly relevant in evaluating the clinical heterogeneity observed in monogenic disorders. The lysosomal storage disorders are Mendelian disorders characterized by a wide spectrum of associated phenotypes, ranging from neonatal presentations to symptoms that develop in late adulthood. Some lack a tight genotype/phenotype correlation. While epigenetics may explain some of the discordant phenotypes encountered in patients with the same lysosomal storage disorder, especially among patients sharing the same genotype, to date, few studies have focused on these mechanisms. We review three common epigenetic mechanisms, DNA methylation, histone modifications, and microRNAs, and highlight their applications to phenotypic variation and therapeutics. Three specific lysosomal storage diseases, Gaucher disease, Fabry disease, and Niemann-Pick type C disease are presented as prototypical disorders with vast clinical heterogeneity that may be impacted by epigenetics. Our goal is to motivate researchers to consider epigenetics as a mechanism to explain the complexities of biological functions and pathologies of these rare disorders.

Source of CpG Depletion in the HIV-1 Genome

The dinucleotide CpG is highly underrepresented in the genome of human immunodeficiency virus type 1 (HIV-1). To identify the source of CpG depletion in the HIV-1 genome, we investigated two biological mechanisms: (1) CpG methylation-induced transcriptional silencing and (2) CpG recognition by Toll-like receptors (TLRs). We hypothesized that HIV-1 has been under selective evolutionary pressure by these mechanisms leading to the reduction of CpG in its genome. A CpG depleted genome would enable HIV-1 to avoid methylation-induced transcriptional silencing and/or to avoid recognition by TLRs that identify foreign CpG sequences. We investigated these two hypotheses by determining the sequence context dependency of CpG depletion and comparing it with that of CpG methylation and TLR recognition. We found that in both human and HIV-1 genomes the CpG motifs flanked by T/A were depleted most and those flanked by C/G were depleted least. Similarly, our analyses of human methylome data revealed that the CpG motifs flanked by T/A were methylated most and those flanked by C/G were methylated least. Given that a similar CpG depletion pattern was observed for the human genome within which CpGs are not likely to be recognized by TLRs, we argue that the main source of CpG depletion in HIV-1 is likely host-induced methylation. Analyses of CpG motifs in over 100 viruses revealed that this unique CpG representation pattern is specific to the human and simian immunodeficiency viruses.

Mechanisms of base substitution mutagenesis in cancer genomes

Cancer genome sequence data provide an invaluable resource for inferring the key mechanisms by which mutations arise in cancer cells, favoring their survival, proliferation and invasiveness. Here we examine recent advances in understanding the molecular mechanisms responsible for the predominant type of genetic alteration found in cancer cells, somatic single base substitutions (SBSs). Cytosine methylation, demethylation and deamination, charge transfer reactions in DNA, DNA replication timing, chromatin status and altered DNA proofreading activities are all now known to contribute to the mechanisms leading to base substitution mutagenesis. We review current hypotheses as to the major processes that give rise to SBSs and evaluate their relative relevance in the light of knowledge acquired from cancer genome sequencing projects and the study of base modifications, DNA repair and lesion bypass. Although gene expression data on APOBEC3B enzymes provide support for a role in cancer mutagenesis through U:G mismatch intermediates, the enzyme preference for single-stranded DNA may limit its activity genome-wide. For SBSs at both CG:CG and YC:GR sites, we outline evidence for a prominent role of damage by charge transfer reactions that follow interactions of the DNA with reactive oxygen species (ROS) and other endogenous or exogenous electron-abstracting molecules.

Cytosine hypomethylation at CHG and CHH sites in the pleiotropic mutants of Mendelian inheritance in Catharanthus roseus

The 5S and 18S rDNA sequences of Catharanthus roseus cv 'Nirmal' (wild type) and its leafless inflorescence (lli), evergreen dwarf (egd) and irregular leaf lamina (ill) single mutants and lli egd, lli ill and egd ill double mutants were characterized. The lli, egd and ill mutants of Mendelian inheritance bore the names after their most conspicuous morphological feature(s). They had been chemically induced and isolated for their salt tolerance. The double mutants were isolated as morphological segregants from crosses between single mutants. The morphological features of the two parents accompanied salt tolerance in the double mutants. All the six mutants were hypomethylated at repeat sequences, upregulated and downregulated for many genes and carried pleiotropic alterations for several traits. Here the 5S and 18S rDNAs of C. roseus were found to be relatively low in cytosine content. Cytosines were preponderantly in CG context (53%) and almost all of them were methylated (97%). The cytosines in CHH and CHG (where H = A, T or C) contexts were largely demethylated (92%) in mutants. The demethylation was attributable to reduced expression of RDR2 and DRM2 led RNA dependant DNA methylation and CMT3 led maintenance methylation pathways. Mutants had gained some cytosines by substitution of C at T sites. These perhaps arose on account of errors in DNA replication, mediated by widespread cytosine demethylation at CHG and CHH sites. It was concluded that the regulation of cytosine ethylation mechanisms was disturbed in the mutants. ILL, EGD and LLI genes were identified as the positive regulators of other genes mediating the RdDM and CMT3 pathways, for establishment and maintenance of cytosine methylation in C. roseus.

Cytosine-to-uracil deamination by SssI DNA methyltransferase

The prokaryotic DNA(cytosine-5)methyltransferase M.SssI shares the specificity of eukaryotic DNA methyltransferases (CG) and is an important model and experimental tool in the study of eukaryotic DNA methylation. Previously, M.SssI was shown to be able to catalyze deamination of the target cytosine to uracil if the methyl donor S-adenosyl-methionine (SAM) was missing from the reaction. To test whether this side-activity of the enzyme can be used to distinguish between unmethylated and C5-methylated cytosines in CG dinucleotides, we re-investigated, using a sensitive genetic reversion assay, the cytosine deaminase activity of M.SssI. Confirming previous results we showed that M.SssI can deaminate cytosine to uracil in a slow reaction in the absence of SAM and that the rate of this reaction can be increased by the SAM analogue 5'-amino-5'-deoxyadenosine. We could not detect M.SssI-catalyzed deamination of C5-methylcytosine ((m5)C). We found conditions where the rate of M.SssI mediated C-to-U deamination was at least 100-fold higher than the rate of (m5)C-to-T conversion. Although this difference in reactivities suggests that the enzyme could be used to identify C5-methylated cytosines in the epigenetically important CG dinucleotides, the rate of M.SssI mediated cytosine deamination is too low to become an enzymatic alternative to the bisulfite reaction. Amino acid replacements in the presumed SAM binding pocket of M.SssI (F17S and G19D) resulted in greatly reduced methyltransferase activity. The G19D variant showed cytosine deaminase activity in E. coli, at physiological SAM concentrations. Interestingly, the C-to-U deaminase activity was also detectable in an E. coli ung (+) host proficient in uracil excision repair.

Transposable elements as genetic regulatory substrates in early development

The abundance and ancient origins of transposable elements (TEs) in eukaryotic genomes has spawned research into the potential symbiotic relationship between these elements and their hosts. In this review, we introduce the diversity of TEs, discuss how distinct classes are uniquely regulated in development, and describe how they appear to have been coopted for the purposes of gene regulation and the orchestration of a number of processes during early embryonic development. Although young, active TEs play an important role in somatic tissues and evolution, we focus mostly on the contributions of the older, fixed elements in mammalian genomes. We also discuss major challenges inherent in the study of TEs and contemplate future experimental approaches to further investigate how they coordinate developmental processes.

The curious chemical biology of cytosine: deamination, methylation, and oxidation as modulators of genomic potential

A multitude of functions have evolved around cytosine within DNA, endowing the base with physiological significance beyond simple information storage. This versatility arises from enzymes that chemically modify cytosine to expand the potential of the genome. Some modifications alter coding sequences, such as deamination of cytosine by AID/APOBEC enzymes to generate immunologic or virologic diversity. Other modifications are critical to epigenetic control, altering gene expression or cellular identity. Of these, cytosine methylation is well understood, in contrast to recently discovered modifications, such as oxidation by TET enzymes to 5-hydroxymethylcytosine. Further complexity results from cytosine demethylation, an enigmatic process that impacts cellular pluripotency. Recent insights help us to propose an integrated DNA demethylation model, accounting for contributions from cytosine oxidation, deamination, and base excision repair. Taken together, this rich medley of alterations renders cytosine a genomic "wild card", whose context-dependent functions make the base far more than a static letter in the code of life.

On the sequence-directed nature of human gene mutation: the role of genomic architecture and the local DNA sequence environment in mediating gene mutations underlying human inherited disease

Different types of human gene mutation may vary in size, from structural variants (SVs) to single base-pair substitutions, but what they all have in common is that their nature, size and location are often determined either by specific characteristics of the local DNA sequence environment or by higher order features of the genomic architecture. The human genome is now recognized to contain "pervasive architectural flaws" in that certain DNA sequences are inherently mutation prone by virtue of their base composition, sequence repetitivity and/or epigenetic modification. Here, we explore how the nature, location and frequency of different types of mutation causing inherited disease are shaped in large part, and often in remarkably predictable ways, by the local DNA sequence environment. The mutability of a given gene or genomic region may also be influenced indirectly by a variety of noncanonical (non-B) secondary structures whose formation is facilitated by the underlying DNA sequence. Since these non-B DNA structures can interfere with subsequent DNA replication and repair and may serve to increase mutation frequencies in generalized fashion (i.e., both in the context of subtle mutations and SVs), they have the potential to serve as a unifying concept in studies of mutational mechanisms underlying human inherited disease.

Characterisation of DNA methylation status using spectroscopy (mid-IR versus Raman) with multivariate analysis

Methylation status plays important roles in the regulation of gene expression and significantly influences the dynamics, bending and flexibility of DNA. The aim of this study was to determine whether attenuated total reflection Fourier-transform infrared (ATR-FTIR) or Raman spectroscopy with subsequent multivariate analysis could determine methylation patterning in oligonucleotides variously containing 5-methylcytosine, cytosine and guanine bases. Applied to Low-E reflective glass slides, 10 independent spectral acquisitions were acquired per oligonucleotide sample. Resultant spectra were baseline-corrected and vector normalised over the 1750 cm(-1) -760 cm(-1) (for ATR-FTIR spectroscopy) or the 1750 cm(-1) -600 cm(-1) (for Raman spectroscopy) regions. Data were then analysed using principal component analysis (PCA) coupled with linear discriminant analysis (LDA). Exploiting this approach, biomolecular signatures enabling sensitive and specific discrimination of methylation patterning were derived. For DNA sequence and methylation analysis, this approach has the potential to be an important tool, especially when material is scarce.

The parental non-equivalence of imprinting control regions during mammalian development and evolution

In mammals, imprinted gene expression results from the sex-specific methylation of imprinted control regions (ICRs) in the parental germlines. Imprinting is linked to therian reproduction, that is, the placenta and imprinting emerged at roughly the same time and potentially co-evolved. We assessed the transcriptome-wide and ontology effect of maternally versus paternally methylated ICRs at the developmental stage of setting of the chorioallantoic placenta in the mouse (8.5dpc), using two models of imprinting deficiency including completely imprint-free embryos. Paternal and maternal imprints have a similar quantitative impact on the embryonic transcriptome. However, transcriptional effects of maternal ICRs are qualitatively focused on the fetal-maternal interface, while paternal ICRs weakly affect non-convergent biological processes, with little consequence for viability at 8.5dpc. Moreover, genes regulated by maternal ICRs indirectly influence genes regulated by paternal ICRs, while the reverse is not observed. The functional dominance of maternal imprints over early embryonic development is potentially linked to selection pressures favoring methylation-dependent control of maternal over paternal ICRs. We previously hypothesized that the different methylation histories of ICRs in the maternal versus the paternal germlines may have put paternal ICRs under higher mutational pressure to lose CpGs by deamination. Using comparative genomics of 17 extant mammalian species, we show here that, while ICRs in general have been constrained to maintain more CpGs than non-imprinted sequences, the rate of CpG loss at paternal ICRs has indeed been higher than at maternal ICRs during evolution. In fact, maternal ICRs, which have the characteristics of CpG-rich promoters, have gained CpGs compared to non-imprinted CpG-rich promoters. Thus, the numerical and, during early embryonic development, functional dominance of maternal ICRs can be explained as the consequence of two orthogonal evolutionary forces: pressure to tightly regulate genes affecting the fetal-maternal interface and pressure to avoid the mutagenic environment of the paternal germline.

Suicidal function of DNA methylation in age-related genome disintegration

This article is dedicated to the 60th anniversary of 5-methylcytosine discovery in DNA. Cytosine methylation can affect genetic and epigenetic processes, works as a part of the genome-defense system and has mutagenic activity; however, the biological functions of this enzymatic modification are not well understood. This review will put forward the hypothesis that the host-defense role of DNA methylation in silencing and mutational destroying of retroviruses and other intragenomic parasites was extended during evolution to most host genes that have to be inactivated in differentiated somatic cells, where it acquired a new function in age-related self-destruction of the genome. The proposed model considers DNA methylation as the generator of 5mC>T transitions that induce 40-70% of all spontaneous somatic mutations of the multiple classes at CpG and CpNpG sites and flanking nucleotides in the p53, FIX, hprt, gpt human genes and some transgenes. The accumulation of 5mC-dependent mutations explains: global changes in the structure of the vertebrate genome throughout evolution; the loss of most 5mC from the DNA of various species over their lifespan and the Hayflick limit of normal cells; the polymorphism of methylation sites, including asymmetric mCpNpN sites; cyclical changes of methylation and demethylation in genes. The suicidal function of methylation may be a special genetic mechanism for increasing DNA damage and the programmed genome disintegration responsible for cell apoptosis and organism aging and death.

Density Functional Study of the Influence of C5 Cytosine Substitution in Base Pairs with Guanine

The present study employs density-functional electronic structure methods to investigate the effect of chemical modification at the C5 position of cytosine. A series of experimentally motivated chemical modifications are considered, including alkyl, halogen, aromatic, fused ring, and strong σ and π withdrawing functional groups. The effect of these modifications on cytosine geometry, electronic structure, proton affinities, gas phase basicities, cytosine-guanine base-pair hydrogen bond network and corresponding nucleophilicity at guanine are examined. Ultimately, these results play a part in dissecting the effect of endogenous cytosine methylation on the reactivity of neighboring guanine toward carcinogens and DNA alkylating agents.

Influence of cytosine methylation on ultraviolet-induced cyclobutane pyrimidine dimer formation in genomic DNA

The ultraviolet (UV) component of sunlight is the main cause of skin cancer. More than 50% of all non-melanoma skin cancers and >90% of squamous cell carcinomas in the US carry a sunlight-induced mutation in the p53 tumor suppressor gene. These mutations have a strong tendency to occur at methylated cytosines. Ligation-mediated PCR (LMPCR) was used to compare at nucleotide resolution DNA photoproduct formation at dipyrimidine sites either containing or lacking a methylated cytosine. For this purpose, we exploited the fact that the X chromosome is methylated in females only on the inactive X chromosome, and that the FMR1 (fragile-X mental retardation 1) gene is methylated only in fragile-X syndrome male patients. Purified genomic DNA was irradiated with UVC (254nm), UVB (290-320nm) or monochromatic UVB (302 and 313nm) to determine the effect of different wavelengths on cyclobutane pyrimidine dimer (CPD) formation along the X-linked PGK1 (phosphoglycerate kinase 1) and FMR1 genes. We show that constitutive methylation of cytosine increases the frequency of UVB-induced CPD formation by 1.7-fold, confirming that methylation per se is influencing the probability of damage formation. This was true for both UVB sources used, either broadband or monochromatic, but not for UVC. Our data prove unequivocally that following UVB exposure methylated cytosines are significantly more susceptible to CPD formation compared with unmethylated cytosines.

Substitution of adenine for guanine in the quadruplex-forming human telomere DNA sequence G(3)(T(2)AG(3))(3)

We have studied the formation and structural properties of quadruplexes of the human telomeric DNA sequence G(3)(T(2)AG(3))(3) and related sequences in which each guanine base was replaced by an adenine base. None of these single base substitutions hindered the formation of antiparallel quadruplexes, as shown by circular dichroism, gel electrophoresis, and UV thermal stability measurements in NaCl solutions. Effect of substitution did differ, however, depending on the position of the substituted base. The A-for-G substitution in the middle quartet of the antiparallel basket scaffold led to the most distorted and least stable structures and these sequences preferred to form bimolecular quadruplexes. Unlike G(3)(T(2)AG(3))(3), no structural transitions were observed for the A-containing analogs of G(3)(T(2)AG(3))(3) when sodium ions were replaced by potassium ions. The basic quadruplex topology remained the same for all sequences studied in both salts. As in vivo misincorporation of A for a G in the telomeric sequence is possible and potassium is a physiological salt, these findings may have biological relevance.

Significance of CpG methylation for solar UV-induced mutagenesis and carcinogenesis in skin

Mutations detected in the p53 gene in human nonmelanoma skin cancers show a highly UV-specific mutation pattern, a dominance of C --> T base substitutions at dipyrimidine sites plus frequent CC --> TT tandem substitutions, indicating a major involvement of solar UV in the skin carcinogenesis. These mutations also have another important characteristic of frequent occurrences at CpG dinucleotide sites, some of which actually show prominent hotspots in the p53 gene. Although mammalian solar UV-induced mutation spectra were studied intensively in the aprt gene using rodent cultured cells and the UV-specific mutation pattern was confirmed, the second characteristic of the p53 mutations in human skin cancers had not been reproduced. However, studies with transgenic mouse systems developed thereafter for mutation research, which harbor methyl CpG-abundant transgenes as mutation markers, yielded complete reproductions of the situation of the human skin cancer mutations in terms of both the UV-specific pattern and the frequent occurrence at CpG sites. In this review, we evaluate the significance of the CpG methylation for solar UV mutagenesis in the mammalian genome, which would lead to skin carcinogenesis. We propose that the UV-specific mutations at methylated CpG sites, C --> T transitions at methyl CpG-associated dipyrimidine sites, are a solar UV-specific mutation signature, and have estimated the wavelength range effective for the solar-UV-specific mutation as 310-340 nm. We also recommend the use of methyl CpG-enriched sequences as mutational targets for studies on solar-UV genotoxicity for human, rather than conventional mammalian mutational marker genes such as the aprt and hprt genes.

Cytosine methylation by the SuaI restriction-modification system: implications for genetic fidelity in a hyperthermophilic archaeon

5-methylcytosine in chromosomal DNA represents a potential source of frequent spontaneous mutation for hyperthermophiles. To determine the relevance of this threat for the archaeon Sulfolobus acidocaldarius, the mode of GGCC methylation by its restriction-modification system, SuaI, was investigated. Distinct isoschizomers of the SuaI endonuclease were used to probe the methylation state of GGCC in native S. acidocaldarius DNA. In addition, the methylation sensitivity of the SuaI endonuclease was determined with synthetic oligonucleotide substrates and modified natural DNAs. The results show that the SuaI system uses N(4) methylation to block cleavage of its recognition site, thereby avoiding the creation of G. T mismatches by spontaneous deamination at extremely high temperature.

An overview of the analysis of DNA methylation in mammalian genomes

DNA methylation at position C5 of the pyrimidine ring of cytosine in mammalian genomes has received a great deal of research interest due to its importance in many biological phenomena. It is associated with events such as epigenetic gene silencing and the maintenance of genome integrity. Aberrant DNA methylation, particularly that of chromosomal regions called CpG islands, is an important step in carcinogenesis. In order to elucidate methylation profiling of complex genomes, various methods have been developed. Many of these methods are based on the differential reactivity of cytosine and 5-methylcytosine to various chemicals. The combined use of these chemical reactions and other preexisting methods has enabled the discrimination of cytosine and 5-methylcytosine in complex genomes. The use of proteins that preferentially bind to methylated DNA has also successfully been used to discriminate between methylated and unmethylated sites. The chemical and structural dissection of the in vivo processes of enzymatic methylation and the binding of methyl-CpG binding proteins provides evidence for the complex mechanisms that nature has acquired. In this review we summarize the methods available for the discrimination between cytosine and 5-methylcytosine in complex genomes.

Spontaneous mutation of the lacI transgene in rodents: absence of species, strain, and insertion-site influence

Comparison of spontaneous mutation spectra derived from different transgenic constructs can provide valuable insights for interpreting the mechanisms of spontaneous mutation. In this study, spontaneous mutation frequencies and spectra of the lacI transgene are compared in the liver of C57BL/6, B6C3F1, and BC-1 mice and F344 rats. Before correction for clonal expansion, the mutant frequency varied from 2.6 +/- 0.45 to 5.0 +/- 2.4 x 10(-5). Correction for potential clonal expansion reduced the range in mutation frequency to between 2.3 +/- 0.45 and 3.5 +/- 2.0 x 10(-5). There is thus no statistical difference in spontaneous mutation frequency between the different strains and species. G:C --> A:T transitions and to a lesser extent, G:C --> T:A transversions dominate the mutational spectra in all of these animals. In three strains of mice, G:C --> A:T transitions account for 50.7-53.3% of mutation, 81.7-83.8% of which involve CpG sites, whereas G:C --> T:A transversions account for 17.8-32.9% of mutations with 43.2-50.0% found at CpG sites. In rats, G:C --> A:T transitions account for 38.0% of the spectra, 70.0% of which involve CpG sites, whereas G:C --> T:A transversions account for 23.0% of the spectra, 70.0% of which involve CpG sites. The distribution of other classes of mutations is also very similar. We conclude that, despite reports about species and strain differences in induced mutation, spontaneous mutations in the lacI transgene appear to be similar, regardless of genomic location, rodent strain, or species. In addition to insights into spontaneous mutation, this study also provides essential data for comparison with and interpretation of induced mutations.

Through a glass, darkly: reflections of mutation from lacI transgenic mice

The study of mutational frequency (Mf) and specificity in aging Big Blue lacI transgenic mice provides a unique opportunity to determine mutation rates (MR) in vivo in different tissues. We found that MR are not static, but rather, vary with the age or developmental stage of the tissue. Although Mf increase more rapidly early in life, MR are actually lower in younger animals than in older animals. For example, we estimate that the changes in Mf are 4.9x10(-8) and 1.1 x 10(-8) mutations/base pair/month in the livers of younger mice (<1. 5 months old) and older mice (> or =1.5 months old), respectively (a 4-fold decrease), and that the MR are 3.9 x 10(-9) and 1.3 x 10(-7) mutations/base pair/cell division, respectively ( approximately 30-fold increase). These data also permit an estimate of the MR of GC --> AT transitions occurring at 5'-CpG-3' (CpG) dinucleotide sequences. Subsequently, the contribution of these transitions to age-related demethylation of genomic DNA can be evaluated. Finally, to better understand the origin of observed Mf, we consider the contribution of various factors, including DNA damage and repair, by constructing a descriptive mutational model. We then apply this model to estimate the efficiency of repair of deaminated 5-methylcytosine nucleosides occurring at CpG dinucleotide sequences, as well as the influence of the Msh2(-/-) DNA repair defect on overall DNA repair efficiency in Big Blue mice. We conclude that even slight changes in DNA repair efficiency could lead to significant increases in mutation frequencies, potentially contributing significantly to human pathogenesis, including cancer.

Distribution of spontaneous CpG-associated G:C --> A:T mutations in the lacZ gene of Muta mice: effects of CpG methylation, the sequence context of CpG sites, and severity of mutations on the activity of the lacZ gene product

In our previous study using transgenic Muta mice, G:C --> A:T transitions at 5'-CG-3' (CpG) sites, which are the most common mammalian spontaneous mutation, were detected in 197 of 330 spontaneous lacZ mutants. These transitions were recovered at only 27 of the 357 mutable G:C pairs within CpG sites where the transition could produce a missense or termination codon in the lacZ gene. To address the underlying mechanism for the uneven distribution of mutated CpG sites, the CpG methylation status of the Muta lacZ gene was analyzed by a bisulfite method. All the CpG sites examined in the coding region were evenly methylated at a high level, and no site-specific methylation was evident. Analysis of the sequence context around the mutated CpG sites, however, revealed that 21 of these 27 sites contained a CpG flanked by a pyrimidine on the 5' side, and that 187 of the 197 mutants resulted from substitutions at these sites. Moreover, we found five hotspots among those sites, the location of which was intimately related to the enzymatic activity of the gene product: one site produced a nonsense codon; three sites, one of which corresponded to the nucleophile at the active site, resided in the substrate-binding pocket; and the other site was located in a region conserved in the beta-galactosidase family. These results strongly suggest that recovery of lacZ mutations at each site largely depend on the adjacent sequence context and the extent to which the mutation damages the enzymatic activity of the gene product.

5-methylcytosine at HpaII sites in p53 is not hypermutable after UVC irradiation

Using a yeast based p53 functional assay we previously demonstrated that the UVC-induced p53 mutation spectrum appears to be indistinguishable from the one observed in Non Melanoma Skin Cancer (NMSC). However, position 742 (codon 248, CpG site) represented the major hot spot in NMSC but was not found mutated in the yeast system. In order to determine whether UVC-induced mutagenic events may be facilitated at methylated cytosine (5mC), a yeast expression vector harbouring a human wild-type p53 cDNA (pLS76) was methylated in vitro by HpaII methylase. Methylation induced 98% protection to HpaII endonuclease. Unmethylated and methylated pLS76 vectors were then UVC irradiated (lambda(max): 254 nm) and transfected into a yeast strain containing the ADE2 gene regulated by a p53-responsive promoter. The results revealed that: (i) 5mC at HpaII sites did not cause any difference in the UVC-induced survival and/or mutagenicity; (ii) none of the 20 mutants derived from methylated pLS76 showed p53 mutations targeted at HpaII sites; (iii) the UVC-induced p53 mutation spectra derived from methylated and unmethylated pLS76 were indistinguishable not only when classes of mutations and hot spots were concerned, but also when compared through a rigorous statistical test to estimate their relatedness (P = 0.85); (iv) the presence of 5mC did not increase the formation of photo-lesions at codon 248, as determined by using a stop polymerase assay. Although based on a limited number of mutants, these results suggest that the mere presence of 5mC at position 742 does not cause a dramatic increase of its mutability after UVC irradiation. We propose that position 742 is a hot spot in NMSC either because of mutagenic events at 5mC caused by other UV components of solarlight and/or because not all the NMSC are directly correlated with UV mutagenesis but may have a "spontaneous" origin.

Reactivity of guanine at m5CpG steps in DNA: evidence for electronic effects transmitted through the base pairs

BACKGROUND

Mitomycin C (MC), a DNA cross-linking and alkylating agent, targets guanines in the m5CpG sequence with 2-3-fold preference over guanines in unmethylated CpG. Benzo[a]pyrenediolepoxide (BPDE) and several other aromatic carcinogens form guanine adducts with an identical selectivity for m5CpG, and in certain cancers G to T transversion mutation 'hotspots' in the p53 tumor suppressor gene are more frequent at this sequence than at guanines in other sequences. MC appears suitable to probe the general mechanism of this selectivity.

RESULTS

A 162-bp DNA fragment containing C, m5C or f5C (5-fluoro cytosine) at all cytosine positions was cross-linked by MC at guanines in CpG steps. The extent of cross-linking increased in the order f5C < C < m5C. Monoalkylation or cross-linking of duplex 12-mer oligonucleotides containing a single CpG, f5CpG or m5CpG step gave yields of adducts that increased in the same order. The rates showed a correlation with the Hammett sigma constant of the methyl and fluoro substituents of the cytosine. Only the base-pair cytosine substituent influenced reactivity of guanine.

CONCLUSIONS

The 2-amino group of guanine in the m5CpG sequence of DNA has a greater nucleophilic reactivity with mitomycin than CpG. Evidence is presented for a novel mechanism: transmission of the electron-donating effect of the 5-methyl substituent of the cytosine to guanine through H-bonding of the m5C.G base pair. The results explain the enhanced reaction of BPDE at m5CpG in DNA and the origin of G-T mutational hotspots in the p53 gene in cancer.

Quantitative whole-genome analysis of DNA-protein interactions by in vivo methylase protection in E. coli

A global methylation-based technique was used to identify, display, and quantitate the in vivo occupancy of numerous protein-binding sites within the Escherichia coli genome. The protein occupancy profiles of these sites showed variation across different growth conditions and genetic backgrounds. Of the 25 sites identified in this study, 24 occurred within 5' noncoding regions. Protein occupancy at 13 of these sites was supported by independent biochemical and genetic evidence. Most of the remaining 12 sites fell upstream of genes with no previously known function. A multivariate statistical analysis was utilized to group such uncharacterized genes with well-characterized ones, providing insights into their function based on a common pattern of transcriptional regulation.

HpaII methyltransferase is mutagenic in Escherichia coli

A genetic reversion assay to study C-to-T mutations within CG sites in DNA is described. It was used to demonstrate that the presence of HpaII methyltransferase (MTase) in Escherichia coli causes a substantial increase in C-to-T mutations at CG sites. This is similar to the known mutagenic effects of E. coli MTase Dcm within its own recognition sequence. With this genetic system, a homolog of an E. coli DNA repair gene in Haemophilus parainfluenzae was tested for antimutagenic activity. Unexpectedly, the homolog was found to have little effect on the reversion frequency. The system was also used to show that HpaII and SssI MTases can convert cytosine to uracil in vitro. These studies define 5-methylcytosine as an intrinsic mutagen and further elaborate the mutagenic potential of cytosine MTases.