CpG Dinucleotide Frequencies Reveal the Role of Host Methylation Capabilities in Parvovirus Evolution

Parvoviruses are rapidly evolving viruses that infect a wide range of hosts, including vertebrates and invertebrates. Extensive methylation of the parvovirus genome has been recently demonstrated. A global pattern of methylation of CpG dinucleotides is seen in vertebrate genomes, compared to “fractional” methylation patterns in invertebrate genomes. It remains unknown if the loss of CpG dinucleotides occurs in all viruses of a given DNA virus family that infect host species spanning across vertebrates and invertebrates. We investigated the link between the extent of CpG dinucleotide depletion among autonomous parvoviruses and the evolutionary lineage of the infected host. We demonstrate major differences in the relative abundance of CpG dinucleotides among autonomous parvoviruses which share similar genome organization and common ancestry, depending on the infected host species. Parvoviruses infecting vertebrate hosts had significantly lower relative abundance of CpG dinucleotides than parvoviruses infecting invertebrate hosts. The strong correlation of CpG dinucleotide depletion with the gain in TpG/CpA dinucleotides and the loss of TpA dinucleotides among parvoviruses suggests a major role for CpG methylation in the evolution of parvoviruses. Our data present evidence that links the relative abundance of CpG dinucleotides in parvoviruses to the methylation capabilities of the infected host. In sum, our findings support a novel perspective of host-driven evolution among autonomous parvoviruses.     

Genome-scale relationships between cytosine methylation and dinucleotide abundances in animals

In mammalian genomes CpGs occur at one-fifth their expected frequency. This is accepted as resulting from cytosine methylation and deamination of 5-methylcytosine leading to TpG and CpA dinucleotides. The corollary that a CpG deficit should correlate with TpG excess has not hitherto been systematically tested at a genomic level. I analyzed genome sequences (human, chimpanzee, mouse, pufferfish, zebrafish, sea squirt, fruitfly, mosquito, and nematode) to do this and generally to assess the hypothesis that CpG deficit, TpG excess, and other data are accountable in terms of 5-methylcytosine mutation. In all methylated genomes local CpG deficit decreases with higher G + C content. Local TpG surplus, while positively associated with G + C level in mammalian genomes but negatively associated with G + C in nonmammalian methylated genomes, is always explicable in terms of the CpG trend under the methylation model. Covariance of dinucleotide abundances with G + C demonstrates that correlation analyses should control for G + C. Doing this reveals a strong negative correlation between local CpG and TpG abundances in methylated genomes, in accord with the methylation hypothesis. CpG deficit also correlates with CpT excess in mammals, which may reflect enhanced cytosine mutation in the context 5'-YCG-3'. Analyses with repeat-masked sequences show that the results are not attributable to repetitive elements.     

DNA methylation and Mycoplasma genomes

DNA methylation is one of the many hypotheses proposed to explain the observed deficiency in CpG dinucleotides in a variety of genomes covering a wide taxonomic distribution. Recent studies challenged the methylation hypothesis on empirical grounds. First, it cannot explain why the Mycoplasma genitalium genome exhibits strong CpG deficiency without DNA methylation. Second, it cannot explain the great variation in CpG deficiency between M. genitalium and M. pneumoniae that also does not have CpG-specific methyltransferase genes. I analyzed the genomic sequences of these Mycoplasma species together with the recently sequenced genomes of M. pulmonis, Ureaplasma urealyticum, and Staphylococcus aureus, and found the results fully compatible with the methylation hypothesis. In particular, I present compelling empirical evidence to support the following scenario. The common ancestor of the three Mycoplasma species has CpG-specific methyltransferases, and has evolved strong CpG deficiency as a result of the specific DNA methylation. Subsequently, this ancestral genome diverged into M. pulmonis and the common ancestor of M. pneumoniae and M. genitalium. M. pulmonis has retained methyltransferases and exhibits the strongest CpG deficiency. The common ancestor lost the methyltransferase gene and then diverged into M. genitalium and M. pneumoniae. M. genitalium and M. pneumoniae, after losing methylation activities, began to regain CpG dinucleotides through random mutation. M. genitalium evolved more slowly than M. pneumoniae, gained relatively fewer CpG dinucleotides, and is more CpG-deficient.     

Distribution of CWG and CCWGG in the Human Genome

Expression of the bacterial CG methyltransferase M?HhaI in mammalian cells appears to generate significant biological effects, while biological effects of the expression of the non-CG methyltransferase M?EcoRII in human cells have not been detected. The association of cytosine methylation with the CG site in mammals is also associated with clustering of CG sites near 5´ control regions (CG-islands) of human genes. Moreover spontaneous deamination of 5-methylcytosine at these sites is thought to lead to the well known deficiency of CG sites in genomes where endogenous CG methyltransferases are expressed. Since these associations are generally taken to imply a biological function for the CG dinucleotide that is associated with its selective methylation by endogenous DNA methylation systems, we have asked whether or not CWG or CCWGG sites are clustered in regions flanking human genes and whether or not an overall deficiency of CWG or CCWGG occurs in the human genome. Using build 36.1, of the human genome, we inspected the regions flanking the 28,501 well known gene loci in the human genome. Our analysis confirmed the expected clustering of CG sites near the 5´ region of known genes and open reading frames. In contrast to the CG site, neither the CWG site nor the CCWGG site recognized by the bacterial methyltransferase M?EcoRII were clustered in any particular region near known genes and open reading frames. Moreover, neither the CCWGG nor the CWG site was depleted in the human genome, again in sharp contrast to the known genomic deficiency of CpG sites. Our findings suggest that in contrast to CG site recognition, human cytosine methyltransferases recognize CWG and CCWGG only at very low frequency if at all.     

Comparative analyses of DNA methylation and sequence evolution using Nasonia genomes

The functional and evolutionary significance of DNA methylation in insect genomes remains to be resolved. Nasonia is well situated for comparative analyses of DNA methylation and genome evolution, since the genomes of a moderately distant outgroup species as well as closely related sibling species are available. Using direct sequencing of bisulfite-converted DNA, we uncovered a substantial level of DNA methylation in 17 of 18 Nasonia vitripennis genes and a strong correlation between methylation level and CpG depletion. Notably, in the sex-determining locus transformer, the exon that is alternatively spliced between the sexes is heavily methylated in both males and females, whereas other exons are only sparsely methylated. Orthologous genes of the honeybee and Nasonia show highly similar relative levels of CpG depletion, despite ~190 My divergence. Densely and sparsely methylated genes in these species also exhibit similar functional enrichments. We found that the degree of CpG depletion is negatively correlated with substitution rates between closely related Nasonia species for synonymous, nonsynonymous, and intron sites. This suggests that mutation rates increase with decreasing levels of germ line methylation. Thus, DNA methylation is prevalent in the Nasonia genome, may participate in regulatory processes such as sex determination and alternative splicing, and is correlated with several aspects of genome and sequence evolution.     

DNA structure constraint is probably a fundamental factor inducing CpG deficiency in bacteria

MOTIVATION: It has been speculated that CpG dinucleotide deficiency in genomes is a consequence of DNA methylation. However, this hypothesis does not adequately explain CpG deficiency in bacteria. The hypothesis based on DNA structure constraint as an alternative explanation was therefore examined. RESULTS: By comparing real bacterial genomes and Markov artificial genomes in the second order, we found that the core structure of a restricted pattern, the TTCGAA pattern, was under represented in low GC content bacterial genomes regardless of CpG dinucleotide level. This is in contrast to the AACGTT pattern, indicating that the counterselection is context-dependent. Further study discovered nine underrepresented patterns that were supposed to be capable of inducing DNA structure constraint. In summary, most of them are in TTCGNA and TTCGAN patterns in both DNA strands. An explanation is also proposed for the strong correlation between GC content and CpG deficiency. The result of random sequence simulation showed that the occurrences of these patterns were correlated with GC content, as well as the percentage of CpG dinucleotides being trapped in these patterns. Finally, we suggest that the degree of counter-selection against these restricted patterns could be influenced by global GC content of a genome.     

Comparative Genomic Study Reveals a Transition from TA Richness in Invertebrates to GC Richness in Vertebrates at CpG Flanking Sites： An Indication for Context-Dependent Mutagenicity of Methylated CpG Sites

Vertebrate genomes are characterized with CpG deficiency, particularly for GCpoor regions. The GC content-related CpG deficiency is probably caused by context-dependent deamination of methylated CpG sites. This hypothesis was examined in this study by comparing nucleotide frequencies at CpG flanking positions among invertebrate and vertebrate genomes. The finding is a transition of nucleotide preference of 5＇ T to 5＇ A at the invertebrate-vertebrate boundary, indicating that a large number of CpG sites with 5＇ Ts were depleted because of global DNA methylation developed in vertebrates. At genome level, we investigated CpG observed/expected （obs/exp） values in 500 bp fragments, and found that higher CpG obs/exp value is shown in GC-poor regions of invertebrate genomes （except sea urchin） but in GC-rich sequences of vertebrate genomes. We next compared GC content at CpG flanking positions with genomic average, showing that the GC content is lower than the average in invertebrate genomes, but higher than that in vertebrate genomes. These results indicate that although 5＇ T and 5＇ A are different in inducing deamination of methylated CpG sites, GC content is even more important in affecting the deamination rate. In all the tests, the results of sea urchin are similar to vertebrates perhaps due to its fractional DNA methylation. CpG deficiency is therefore suggested to be mainly a result of high mutation rates of methylated CpG sites in GC-poor regions.     

Specificity of MLL1 and TET3 CXXC domains towards naturally occurring cytosine modifications

CXXC domains have traditionally been considered as CpG specific DNA binding domains that are repelled by cytosine modifications. This view has recently been challenged by the demonstration that CXXC domain of TET3 has relaxed sequence specificity and binds with the highest affinity to symmetric DNA duplex containing 5caCpG. Here, we present a comparative analysis of the MLL1-CXXC and TET3-CXXC sequence specificity and tolerance to cytosine modifications (5-methyl, 5-hydroxymethyl, 5-formyl, 5-carboxyl) in CpG and non-CpG context. For the first time, we take into consideration possible interference from cytosine bases elsewhere in the sequence. We show that despite similar overall structure, MLL1-CXXC has greater sequence and modification specificity than TET3-CXXC. MLL1-CXXC is specific only for CpG and does not tolerate any cytosine modifications. In contrast, TET3-CXXC does not require the CpG context of cytosine bases. Methyl-, formyl- and carboxyl-modifications are tolerated by TET3-CXXC, but only preceding G. Based on our and other data we propose a parsimonious model of MLL1-CXXC and TET3-CXXC DNA binding. This model explains why the binding of modified DNA duplexes by TET3-CXXC requires in some cases a register shift and is therefore context-dependent.     

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.     

Thymine DNA glycosylase is essential for active DNA demethylation by linked deamination-base excision repair

DNA methylation is a major epigenetic mechanism for gene silencing. Whereas methyltransferases mediate cytosine methylation, it is less clear how unmethylated regions in mammalian genomes are protected from de novo methylation and whether an active demethylating activity is involved. Here, we show that either knockout or catalytic inactivation of the DNA repair enzyme thymine DNA glycosylase (TDG) leads to embryonic lethality in mice. TDG is necessary for recruiting p300 to retinoic acid (RA)-regulated promoters, protection of CpG islands from hypermethylation, and active demethylation of tissue-specific developmentally and hormonally regulated promoters and enhancers. TDG interacts with the deaminase AID and the damage response protein GADD45a. These findings highlight a dual role for TDG in promoting proper epigenetic states during development and suggest a two-step mechanism for DNA demethylation in mammals, whereby 5-methylcytosine and 5-hydroxymethylcytosine are first deaminated by AID to thymine and 5-hydroxymethyluracil, respectively, followed by TDG-mediated thymine and 5-hydroxymethyluracil excision repair.     

DNA Methylation: Bisulphite Modification and Analysis

Epigenetics describes the heritable changes in gene function that occur independently to the DNA sequence. The molecular basis of epigenetic gene regulation is complex, but essentially involves modifications to the DNA itself or the proteins with which DNA associates. The predominant epigenetic modification of DNA in mammalian genomes is methylation of cytosine nucleotides (5-MeC). DNA methylation provides instruction to gene expression machinery as to where and when the gene should be expressed. The primary target sequence for DNA methylation in mammals is 5''-CpG-3'' dinucleotides (Figure 1). CpG dinucleotides are not uniformly distributed throughout the genome, but are concentrated in regions of repetitive genomic sequences and CpG "islands" commonly associated with gene promoters (Figure 1). DNA methylation patterns are established early in development, modulated during tissue specific differentiation and disrupted in many disease states including cancer. To understand the biological role of DNA methylation and its role in human disease, precise, efficient and reproducible methods are required to detect and quantify individual 5-MeCs.This protocol for bisulphite conversion is the "gold standard" for DNA methylation analysis and facilitates identification and quantification of DNA methylation at single nucleotide resolution. The chemistry of cytosine deamination by sodium bisulphite involves three steps (Figure 2). (1) Sulphonation: The addition of bisulphite to the 5-6 double bond of cytosine (2) Hydrolic Deamination: hydrolytic deamination of the resulting cytosine-bisulphite derivative to give a uracil-bisulphite derivative (3) Alkali Desulphonation: Removal of the sulphonate group by an alkali treatment, to give uracil. Bisulphite preferentially deaminates cytosine to uracil in single stranded DNA, whereas 5-MeC, is refractory to bisulphite-mediated deamination. Upon PCR amplification, uracil is amplified as thymine while 5-MeC residues remain as cytosines, allowing methylated CpGs to be distinguished from unmethylated CpGs by presence of a cytosine "C" versus thymine "T" residue during sequencing.DNA modification by bisulphite conversion is a well-established protocol that can be exploited for many methods of DNA methylation analysis. Since the detection of 5-MeC by bisulphite conversion was first demonstrated by Frommer et al.¹ and Clark et al.², methods based around bisulphite conversion of genomic DNA account for the majority of new data on DNA methylation. Different methods of post PCR analysis may be utilized, depending on the degree of specificity and resolution of methylation required. Cloning and sequencing is still the most readily available method that can give single nucleotide resolution for methylation across the DNA molecule.     

C5 cytosine methylation at CpG sites enhances sequence selectivity of mitomycin C-DNA bonding

We have established that UvrABC nuclease is equally efficient in cutting mitomycin C (MC)-DNA monoadducts formed at different sequences and that the degree of UvrABC cutting represents the extent of drug-DNA bonding. Using this method we determined the effect of C5 cytosine methylation on the DNA monoalkylation by MC and the related analogues N-methyl-7-methoxyaziridinomitosene (MS-NMA) and 10-decarbamoylmitomycin C (DC-MC). We have found that C5 cytosine methylation at CpG sites greatly enhances MC and MS-NMA DNA adduct formation at those sites while reducing adduct formation at non-CpG sequences. In contrast, although DC-MC DNA bonding at CpG sites is greatly enhanced by CpG methylation, its bonding at non-CpG sequences is not appreciably affected. These cumulative results suggest that C5 cytosine methylation at CpG sites enhances sequence selectivity of drug-DNA bonding. We propose that the methylation pattern and status (hypo- or hypermethylation) of genomic DNA may determine the cells' susceptibility to MC and its analogues, and these effects may, in turn, play a crucial role in the antitumor activities of the drugs.     

Dnmt2 is the most evolutionary conserved and enigmatic cytosine DNA methyltransferase in eukaryotes

Dnmt2 is the most strongly conserved cytosine DNA methyltransferase in eukaryotes. It has been found in all organisms possessing methyltransferases of the Dnmt1 and Dnmt3 families, whereas in many others Dnmt2 is the sole cytosine DNA methyltransferase. The Dnmt2 molecule contains all conserved motifs of cytosine DNA methyltransferases. It forms 3D complexes with DNA very similar to those of bacterial DNA methyltransferases and performs cytosine methylation by a catalytic mechanism common to all cytosine DNA methyltransferases. Catalytic activity of the purified Dnmt2 with DNA substrates is very low and could hardly be detected in direct biochemical assays. Dnmt2 is the sole cytosine DNA methyltransferase in Drosophila and other dipteran insects. Its overexpression as a transgene leads to DNA hypermethylation in all sequence contexts and to an extended life span. On the contrary, a null-mutation of the Dnmt2 gene leads to a diminished life span, though no evident anomalies in development are observed. Dnmt2 is also the sole cytosine DNA methyltransferase in several protists. Similar to Drosophila these protists have a very low level of DNA methylation. Some limited genome compartments, such as transposable sequences, are probably the methylation targets in these organisms. Dnmt2 does not participate in genome methylation in mammals, but seems to be an RNA methyltransferase modifying the 38th cytosine residue in anticodon loop of certain tRNAs. This modification enhances stability of tRNAs, especially in stressful conditions. Dnmt2 is the only enzyme known to perform RNA methylation by a catalytic mechanism characteristic of DNA methyltransferases. The Dnmt2 activity has been shown in mice to be necessary for paramutation establishment, though the precise mechanisms of its participation in this form of epigenetic heredity are unknown. It seems likely, that either of the two Dnmt2 activities could become a predominant one during the evolution of different species. The high level of the Dnmt2 evolutionary conservation proves its activity to have a significant adaptive value in natural environment.     

Characterization of the level,target sites and inheritance of cytosine methylation in tomato nuclear DNA

The tomato nuclear genome was determined to have a G+C content of 37% which is among the lowest reported for any plant species. Non-coding regions have a G+C content even lower (32% average) whereas coding regions are considerably richer in G+C (46%).5-methyl cytosine was the only modified base detected and on average 23% of the cytosine residues are methylated. Immature tissues and protoplasts have significantly lower levels of cytosine methylation (average 20%) than mature tissues (average 25%). Mature pollen has an intermediate level of methylation (22%). Seeds gave the highest value (27%), suggesting de novo methylation after pollination and during seed development.Based on isoschizomer studies we estimate 55% of the CpG target sites (detected by Msp I/Hpa II) and 85% of the CpNpG target sites (detected by Bst NI/Eco RI)are methylated. Unmethylated target sites (both CpG and CpNpG) are not randomly distributed throughout the genome, but frequently occur in clusters. These clusters resemble CpG islands recently reported in maize and tobacco.The low G+C content and high levels of cytosine methylation in tomato may be due to previous transitions of 5mCT. This is supported by the fact that G+C levels are lowest in non-coding portions of the genome in which selection is relaxed and thus transitions are more likely to be tolerated. This hypothesis is also supported by the general deficiency of methylation target sites in the tomato genome, especially in non-coding regions.Using methylation isoschizomers and RFLP analysis we have also determined that polymorphism between plants, for cytosine methylation at allelic sites, is common in tomato. Comparing DNA from two tomato species, 20% of the polymorphisms detected by Bst NI/Eco RII could be attributed to differential methylation at the CpNpG target sites. With Msp I/Hpa II, 50% of the polymorphisms were attributable to methylation (CpG and CpNpG sites). Moreover, these polymorphisms were demonstrated to be inherited in a mendelian fashion and to co-segregate with the methylation target site and thus do not represent variation for transacting factors that might be involved in methylation of DNA. The potential role of heritable methylation polymorphism in evolution of gene regulation and in RFLP studies is discussed.     

Inheritance of an epigenetic mark: the CpG DNA methyltransferase 1 is required for de novo establishment of a complex pattern of non-CpG methylation

Site-specific methylation of cytosines is a key epigenetic mark of vertebrate DNA. While a majority of the methylated residues are in the symmetrical (meC)pG:Gp(meC) configuration, a smaller, but significant fraction is found in the CpA, CpT and CpC asymmetric (non-CpG) dinucleotides. CpG methylation is reproducibly maintained by the activity of the DNA methyltransferase 1 (Dnmt1) on the newly replicated hemimethylated substrates (meC)pG:GpC. On the other hand, establishment and hereditary maintenance of non-CpG methylation patterns have not been analyzed in detail. We previously reported the occurrence of site- and allele-specific methylation at both CpG and non-CpG sites. Here we characterize a hereditary complex of non-CpG methylation, with the transgenerational maintenance of three distinct profiles in a constant ratio, associated with extensive CpG methylation. These observations raised the question of the signal leading to the maintenance of the pattern of asymmetric methylation. The complete non-CpG pattern was reinstated at each generation in spite of the fact that the majority of the sperm genomes contained either none or only one methylated non-CpG site. This observation led us to the hypothesis that the stable CpG patterns might act as blueprints for the maintenance of non-CpG DNA methylation. As predicted, non-CpG DNA methylation profiles were abrogated in a mutant lacking Dnmt1, the enzymes responsible for CpG methylation, but not in mutants defective for either Dnmt3a or Dnmt2.     

YccW is the m5C methyltransferase specific for 23S rRNA nucleotide 1962

Methylation at the 5-position of cytosine [m⁵C (5-methylcytidine)] occurs at three RNA nucleotides in Escherichia coli. All these modifications are at highly conserved nucleotides in the rRNAs, and each is catalyzed by its own m⁵C methyltransferase enzyme. Two of the enzymes, RsmB and RsmF, are already known and methylate 16S rRNA at nucleotides C967 and C1407, respectively. Here, we report the identity of the third E. coli m⁵C methyltransferase. Analysis of rRNAs by matrix-assisted laser desorption/ionization mass spectrometry showed that inactivation of the yccW gene leads to loss of m⁵C methylation at nucleotide 1962 in E. coli 23S rRNA. This methylation is restored by complementing the knockout strain with a plasmid-encoded copy of the yccW gene. Purified recombinant YccW protein retains its specificity for C1962 in vitro and methylates naked 23S rRNA isolated from the yccW knockout strain. However, YccW does not methylate assembled 50S subunits, and this is somewhat surprising as the published crystal structures show nucleotide C1962 to be fully accessible at the subunit interface. YccW-directed methylation at nucleotide C1962 is conserved in bacteria, and loss of this methylation in E. coli marginally reduces its growth rate. YccW had previously eluded identification because it displays only limited sequence similarity to the m⁵C methyltransferases RsmB and RsmF and is in fact more similar to known m⁵U (5-methyluridine) RNA methyltransferases. In keeping with the previously proposed nomenclature system for bacterial rRNA methyltransferases, yccW is now designated as the rRNA large subunit methyltransferase gene rlmI.