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.     

Unmethylated CpG islands associated with genes in higher plant DNA

The genomes of many higher plant species are the most highly methylated among eukaryotes. We report here that in spite of their heavy methylation, genomic DNAs from four plant species contain a fraction that is very rich in non-methylated sites. The fraction was characterized in maize where it represents about 2.5% of the total nuclear genome. In order to establish the genomic origin of the fraction, three maize genes containing clustered CpG were tested for methylation and were found to be non-methylated in the CpG-rich regions. By contrast, tested CpGs were methylated in a gene whose sequence showed no clustering of CpG. These observations suggest that the CpG-rich fraction of plants is at least partially derived from non-methylated regions that are associated with genes. A similar phenomenon has been described in vertebrate genomes. We discuss the evolution of CpG islands in both groups of organisms, and their possible uses in mapping and gene isolation in plants.     

The loss of CpG dinucleotides from DNA. I. Methylated and non-methylated genome compartments in eukaryotes with different levels of 5-methylcytosine in DNA

The methylation of cytosine residues in CpG significantly increases the frequency of m5CpG----TpG transitions in DNA and CpG dinucleotides are eliminated from the genome (CpG-suppression). In the millions of years of vertebrates evolution about 3 mol% of 5-methylcytosine have disappeared from their genome, i.e., 2-3-fold more than the amount persisting in the DNA of the now extant species. A computer analysis has been carried out of neighboring b.p. frequencies in more than 2500 sequenced genes of different species in the EMBL bank with an overall extension of over 3000 kb. It has been found that CpG methylated sites exhibit a highly irregular distribution pattern in the genome of eucaryotes. The majority of the vertebrate sequences (92%) bears the impress of a significant lack of CpG and an excess of TpG+CpA; therefore they may be referred to the genome methylated compartment. A group of genes has been discovered (about 8%) where CpG must have never been subjected to methylation. In invertebrates, such a nonmethylated compartment makes up 59% of the genome and in eubacteria--85%. A brief list of genes, belonging to the methylated and the non-methylated compartments of the invertebrate and yeast genome, is given. It has been established that the mean value of CpG-suppression in genes is directly proportional to the methylation level of total DNA in different species.     

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.     

Methylation sites in angiosperm genes

Summary The extent to which CpG dinucleotides were depleted in a large set of angiosperm genes was, on average, very similar to the extent of CpG depletion in total angiosperm genomic DNA and far less than the extent of CpG depletion in vertebrate genes. Gene sequences from Arabidopsis thaliana, a dicotyledonous species with relatively low levels of total 5-methylcytosine, were just as CpG depleted as the angiosperm genes in general. Furthermore, levels of TpG and CpA, the potential deamination mutation products of methylated CpG, were elevated in A. thaliana genes, supporting a high rate of deamination mutation as the cause of the CpG deficiency. Using a method that takes into account the dinucleotide frequencies within each sequence of interest, we calculated the expected frequencies of CpNpG trinucleotides, which are also highly methylated in angiosperm genomes. CpNpG trinucleotides were not extensively enriched or depleted in the angiosperm genes. Two hypotheses could account for our results. Differential depletion of CpG and CpNpG within angiosperm genes and differential depletion of CpG in angiosperm and vertebrate genes could arise from different efficiencies of mismatch repair or from different levels of cytosine methylation in the cell lineages that contribute to germ cells.Offprint requests to: M. Gardiner-Garden     

Structural and functional features of the 5-methylcytosine distribution in the eukaryotic genome

DNA methylation is an integral part of the mechanism of a remodeling and modification of the chromatin structure. The global complex net of chromatin modification and remodeling reactions is still to be determined, and studies of the mechanisms controlling the epigenetic processes of histone modification and DNA methylation are in their infancy. Cytosine methylation occurs predominantly in CpG sequences of the eukaryotic genome, and it also takes place at symmetric CpHpG and nonsymmetric CpHpH sites (where H is A, T, or C). The modification efficiency of the three types of DNA methylation sites depends on their genomic localization. Different regions of the eukaryotic genome are remarkable for their methylation features: CpG-islands, CpG-island shores, differentially methylated regions of imprinted genes, and regions of nonalternative site-specific modification. The three canonical sites (CpG, CpHpG, and CpHpH) differ in DNA methylation efficiency depending on their nucleotide context. An epigenetic code of DNA methylation can be assumed with context differences playing a specific functional role. The review summarizes the main up-to-date data on the structural and functional features of site-specific cytosine methylation in eukaryotic genomes. Pathogenesis-related alterations in the methylation pattern of the eukaryotic genome are considered.     

A species-generalized probabilistic model-based definition of CpG islands

The DNA of most vertebrates is depleted in CpG dinucleotides, the target for DNA methylation. The remaining CpGs tend to cluster in regions referred to as CpG islands (CGI). CGI have been useful as marking functionally relevant epigenetic loci for genome studies. For example, CGI are enriched in the promoters of vertebrate genes and thought to play an important role in regulation. Currently, CGI are defined algorithmically as an observed-to-expected ratio (O/E) of CpG greater than 0.6, G+C content greater than 0.5, and usually but not necessarily greater than a certain length. Here we find that the current definition leaves out important CpG clusters associated with epigenetic marks, relevant to development and disease, and does not apply at all to nonvertabrate genomes. We propose an alternative Hidden Markov model-based approach that solves these problems. We fit our model to genomes from 30 species, and the results support a new epigenomic view toward the development of DNA methylation in species diversity and evolution. The O/E of CpG in islands and nonislands segregated closely phylogenetically and showed substantial loss in both groups in animals of greater complexity, while maintaining a nearly constant difference in CpG O/E between islands and nonisland compartments. Lists of CGI for some species are available at http://www.rafalab.org.     

The loss of CpG dinucleotides from DNA. III. Methylation and evolution of histone genes

From nucleotide sequences of more than 70 histones genes in 15 species of eucaryotes the probable frequency was determined for CpG----TpG + CpA substitutions, occurring as a result of deamination of 5-methylcytosine residues in DNA. It was found that histone genes differ in the character of CpG methylation with respect to the species studied and may be divided into three groups differing in the value of CpG suppression. In one of them, M-, CpG dinucleotides must have not been methylated throughout the existence of these genes; in another, M+, nearly every other CpG has undergone transition. In the third group, M +/-, no more than 20% of CpG have steadily undergone methylation (and mutation). The CpG deficiency in M+ and M +/- histone genes is in general proportional to the level of methylation of total DNA in different species. It has been noted that the genes of different core histones in the same organism are characterized, as a rule, by the same type of CpG methylation and belong to the same group. Genes H1 and H5 show a higher level of CpG suppression and thus have a higher degree of methylation than the genes of core histones from the same organism. The most conserved among the histone genes, those for H3 and H4 in particular, must have not been methylated in the majority of the species studied. The distribution of methylated and non-methylated spacers and coding sequences of histone genes of man, mouse, hen and yeast reveals a mosaic pattern. It has been found that 5'-flanked regions in most cases are methylated more than respective genes, while the G + C content in them is significantly lower, compared with the coding gene sequences. The absence of methylation in the 5'-regulatory regions does not appear to be mandatory for histone genes. It has been established that the genes of the same histones may differ in the level of methylation even in more or less closely related species. Group M- comprises genes of core histones of man, hen, sea urchin, Drosophila, Neurospora and wheat; group M +/- includes analogous genes of mouse, Xenopus, trout and sea urchins. The results obtained testify against the possible universal involvement of methylation in the regulation of histone gene expression.     

Presence and role of cytosine methylation in DNA viruses of animals

Nucleotide composition varies greatly among DNA viruses of animals, yet the evolutionary pressures and biological mechanisms driving these patterns are unclear. One of the most striking discrepancies lies in the frequency of CpG (the dinucleotide CG, linked by a phosphate group), which is underrepresented in most small DNA viruses (those with genomes below 10 kb) but not in larger DNA viruses. Cytosine methylation might be partially responsible, but research on this topic has focused on a few virus groups. For several viruses that integrate their genome into the host genome, the methylation status during this stage has been studied extensively, and the relationship between methylation and viral-induced tumor formation has been examined carefully. However, for actively replicating viruses—particularly small DNA viruses—the methylation status of CpG motifs is rarely known and the effects on the viral life cycle are obscure. In vertebrate host genomes, most cytosines at CpG sites are methylated, which in vertebrates acts to regulate gene expression and facilitates the recognition of unmethylated, potentially pathogen-associated DNA. Here we briefly introduce cytosine methylation before reviewing what is currently known about CpG methylation in DNA viruses.     

Evolutionary conservation of DNA methylation in CpG sites within ultraconserved noncoding elements

Ultraconserved noncoding elements (UCNEs) constitute less than 1 Mb of vertebrate genomes and are impervious to accumulating mutations. About 4000 UCNEs exist in vertebrate genomes, each at least 200 nucleotides in length, sharing greater than 95% sequence identity between human and chicken. Despite extreme sequence conservation over 400 million years of vertebrate evolution, we show both ordered interspecies and within-species interindividual variation in DNA methylation in these regions. Here, we surveyed UCNEs with high CpG density in 56 species finding half to be intermediately methylated and the remaining near 0% or 100%. Intermediately methylated UCNEs displayed a greater range of methylation between mouse tissues. In a human population, most UCNEs showed greater variation than the LINE1 transposon, a frequently used epigenetic biomarker. Global methylation was found to be inversely correlated to hydroxymethylation across 60 vertebrates. Within UCNEs, DNA methylation is flexible, conserved between related species, and relaxed from the underlying sequence selection pressure, while remaining heritable through speciation.     

Allele-Specific Transcriptome and Methylome Analysis Reveals Stable Inheritance and Cis-Regulation of DNA Methylation in Nasonia

Gene expression divergence between closely related species could be attributed to both cis- and trans- DNA sequence changes during evolution, but it is unclear how the evolutionary dynamics of epigenetic marks are regulated. In eutherian mammals, biparental DNA methylation marks are erased and reset during gametogenesis, resulting in paternal or maternal imprints, which lead to genomic imprinting. Whether DNA methylation reprogramming exists in insects is not known. Wasps of the genus Nasonia are non-social parasitoids that are emerging as a model for studies of epigenetic processes in insects. In this study, we quantified allele-specific expression and methylation genome-wide in Nasonia vitripennis and Nasonia giraulti and their reciprocal F1 hybrids. No parent-of-origin effect in allelic expression was found for >8,000 covered genes, suggesting a lack of genomic imprinting in adult Nasonia. As we expected, both significant cis- and trans- effects are responsible for the expression divergence between N. vitripennis and N. giraulti. Surprisingly, all 178 differentially methylated genes are also differentially methylated between the two alleles in F1 hybrid offspring, recapitulating the parental methylation status with nearly 100% fidelity, indicating the presence of strong cis-elements driving the target of gene body methylation. In addition, we discovered that total and allele-specific expression are positively correlated with allele-specific methylation in a subset of the differentially methylated genes. The 100% cis-regulation in F1 hybrids suggests the methylation machinery is conserved and DNA methylation is targeted by cis features in Nasonia. The lack of genomic imprinting and parent-of-origin differentially methylated regions in Nasonia, together with the stable inheritance of methylation status between generations, suggests either a cis-regulatory motif for methylation at the DNA level or highly stable inheritance of an epigenetic signal in Nasonia.