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.     

Patterns of evolution and host gene mimicry in influenza and other RNA viruses

It is well known that the dinucleotide CpG is under-represented in the genomic DNA of many vertebrates. This is commonly thought to be due to the methylation of cytosine residues in this dinucleotide and the corresponding high rate of deamination of 5-methycytosine, which lowers the frequency of this dinucleotide in DNA. Surprisingly, many single-stranded RNA viruses that replicate in these vertebrate hosts also have a very low presence of CpG dinucleotides in their genomes. Viruses are obligate intracellular parasites and the evolution of a virus is inexorably linked to the nature and fate of its host. One therefore expects that virus and host genomes should have common features. In this work, we compare evolutionary patterns in the genomes of ssRNA viruses and their hosts. In particular, we have analyzed dinucleotide patterns and found that the same patterns are pervasively over- or under-represented in many RNA viruses and their hosts suggesting that many RNA viruses evolve by mimicking some of the features of their host's genes (DNA) and likely also their corresponding mRNAs. When a virus crosses a species barrier into a different host, the pressure to replicate, survive and adapt, leaves a footprint in dinucleotide frequencies. For instance, since human genes seem to be under higher pressure to eliminate CpG dinucleotide motifs than avian genes, this pressure might be reflected in the genomes of human viruses (DNA and RNA viruses) when compared to those of the same viruses replicating in avian hosts. To test this idea we have analyzed the evolution of the influenza virus since 1918. We find that the influenza A virus, which originated from an avian reservoir and has been replicating in humans over many generations, evolves in a direction strongly selected to reduce the frequency of CpG dinucleotides in its genome. Consistent with this observation, we find that the influenza B virus, which has spent much more time in the human population, has adapted to its human host and exhibits an extremely low CpG dinucleotide content. We believe that these observations directly show that the evolution of RNA viral genomes can be shaped by pressures observed in the host genome. As a possible explanation, we suggest that the strong selection pressures acting on these RNA viruses are most likely related to the innate immune response and to nucleotide motifs in the host DNA and RNAs.     

Why is CpG suppressed in the genomes of virtually all small eukaryotic viruses but not in those of large eukaryotic viruses?

Dinucleotide over- and underrepresentation is evaluated in all available completely sequenced DNA or RNA viral genomes, ranging in size from 3 to 250 kb (available RNA viruses fall into the small-virus category). The dinucleotide CpG is statistically underrepresented (suppressed) in all but four of the small viruses (more than 75 with lengths of < 30 kb) but has normal relative abundances in most large viruses (> or = 30 kb). Most retrotransposons in eukaryotic species also show low CpG relative abundances. Interpretations, especially in some cases of DNA viruses or viruses with a DNA intermediate, might relate to methylation effects and modes of viral integration and excision. Other possible contributing factors relate to dinucleotide stacking energies, special mutation mechanisms, and evolutionary events.     

A human CpG island randomly inserted into a plant genome is protected from methylation

In vertebrate genomes the dinucleotide CpG is heavily methylated, except in CpG islands, which are normally unmethylated. It is not clear why the CpG islands are such poor substrates for DNA methyltransferase. Plant genomes display methylation, but otherwise the genomes of plants and animals represent two very divergent evolutionary lines. To gain a further understanding of the resistance of CpG islands to methylation, we introduced a human CpG island from the proteasome-like subunit I gene into the genome of the plant Arabidopsis thaliana. Our results show that prevention of methylation is an intrinsic property of CpG islands, recognized even if a human CpG island is transferred to a plant genome. Two different parts of the human CpG island – the promoter region/ first exon and exon2–4 – both displayed resistance against methylation, but the promoter/ exon1 construct seemed to be most resistant. In contrast, certain sites in a plant CpG-rich region used as a control transgene were always methylated. The frequency of silencing of the adjacent nptII (Km^R) gene in the human CpG constructs was lower than observed for the plant CpG-rich region. These results have implications for understanding DNA methylation, and for construction of vectors that will reduce transgene silencing.     

Evolutionary Basis of Codon Usage and Nucleotide Composition Bias in Vertebrate DNA Viruses

Understanding the extent and causes of biases in codon usage and nucleotide composition is essential to the study of viral evolution, particularly the interplay between viruses and host cells or immune responses. To understand the common features and differences among viruses we analyzed the genomic characteristics of a representative collection of all sequenced vertebrate-infecting DNA viruses. This revealed that patterns of codon usage bias are strongly correlated with overall genomic GC content, suggesting that genome-wide mutational pressure, rather than natural selection for specific coding triplets, is the main determinant of codon usage. Further, we observed a striking difference in CpG content between DNA viruses with large and small genomes. While the majority of large genome viruses show the expected frequency of CpG, most small genome viruses had CpG contents far below expected values. The exceptions to this generalization, the large gammaherpesviruses and iridoviruses and the small dependoviruses, have sufficiently different life-cycle characteristics that they may help reveal some of the factors shaping the evolution of CpG usage in viruses. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Nicolas Galtier]     

DNA Sequence Explains Seemingly Disordered Methylation Levels in Partially Methylated Domains of Mammalian Genomes

For the most part metazoan genomes are highly methylated and harbor only small regions with low or absent methylation. In contrast, partially methylated domains (PMDs), recently discovered in a variety of cell lines and tissues, do not fit this paradigm as they show partial methylation for large portions (20%–40%) of the genome. While in PMDs methylation levels are reduced on average, we found that at single CpG resolution, they show extensive variability along the genome outside of CpG islands and DNase I hypersensitive sites (DHS). Methylation levels range from 0% to 100% in a roughly uniform fashion with only little similarity between neighboring CpGs. A comparison of various PMD-containing methylomes showed that these seemingly disordered states of methylation are strongly conserved across cell types for virtually every PMD. Comparative sequence analysis suggests that DNA sequence is a major determinant of these methylation states. This is further substantiated by a purely sequence based model which can predict 31% (R²) of the variation in methylation. The model revealed CpG density as the main driving feature promoting methylation, opposite to what has been shown for CpG islands, followed by various dinucleotides immediately flanking the CpG and a minor contribution from sequence preferences reflecting nucleosome positioning. Taken together we provide a reinterpretation for the nucleotide-specific methylation levels observed in PMDs, demonstrate their conservation across tissues and suggest that they are mainly determined by specific DNA sequence features.     

Different measures of “genome-wide” DNA methylation exhibit unique properties in placental and somatic tissues

DNA methylation of CpGs located in two types of repetitive elements—LINE1 (L1) and Alu—is used to assess “global” changes in DNA methylation in studies of human disease and environmental exposure. L1 and Alu contribute close to 30% of all base pairs in the human genome and transposition of repetitive elements is repressed through DNA methylation. Few studies have investigated whether repetitive element DNA methylation is associated with DNA methylation at other genomic regions, or the biological and technical factors that influence potential associations. Here, we assess L1 and Alu DNA methylation by Pyrosequencing of consensus sequences and using subsets of probes included in the Illumina Infinium HumanMethylation27 BeadChip array. We show that evolutionary age and assay method affect the assessment of repetitive element DNA methylation. Additionally, we compare Pyrosequencing results for repetitive elements to average DNA methylation of CpG islands, as assessed by array probes classified into strong, weak and non-islands. We demonstrate that each of these dispersed sequences exhibits different patterns of tissue-specific DNA methylation. Correlation of DNA methylation suggests an association between L1 and weak CpG island DNA methylation in some of the tissues examined. We caution, however, that L1, Alu and CpG island DNA methylation are distinct measures of dispersed DNA methylation and one should not be used in lieu of another. Analysis of DNA methylation data is complex and assays may be influenced by environment and pathology in different or complementary ways.     

Parvovirus b19 DNA CpG dinucleotide methylation and epigenetic regulation of viral expression

CpG DNA methylation is one of the main epigenetic modifications playing a role in the control of gene expression. For DNA viruses whose genome has the ability to integrate in the host genome or to maintain as a latent episome, a correlation has been found between the extent of DNA methylation and viral quiescence. No information is available for Parvovirus B19, a human pathogenic virus, which is capable of both lytic and persistent infections. Within Parvovirus B19 genome, the inverted terminal regions display all the characteristic signatures of a genomic CpG island; therefore we hypothesised a role of CpG dinucleotide methylation in the regulation of viral genome expression.The analysis of CpG dinucleotide methylation of Parvovirus B19 DNA was carried out by an aptly designed quantitative real-time PCR assay on bisulfite-modified DNA. The effects of CpG methylation on the regulation of viral genome expression were first investigated by transfection of either unmethylated or in vitro methylated viral DNA in a model cell line, showing that methylation of viral DNA was correlated to lower expression levels of the viral genome. Then, in the course of in vitro infections in different cellular environments, it was observed that absence of viral expression and genome replication were both correlated to increasing levels of CpG methylation of viral DNA. Finally, the presence of CpG methylation was documented in viral DNA present in bioptic samples, indicating the occurrence and a possible role of this epigenetic modification in the course of natural infections.The presence of an epigenetic level of regulation of viral genome expression, possibly correlated to the silencing of the viral genome and contributing to the maintenance of the virus in tissues, can be relevant to the balance and outcome of the different types of infection associated to Parvovirus B19.     

Causal analysis of CpG suppression in the Mycoplasma genome

Some bacterial genomes are known to have low CpG dinucleotide frequencies. While their causes are not clearly understood, the frequency of CpG is suppressed significantly in the genome of Mycoplasma genitalium, but not in that of Mycoplasma pneumoniae. We compared orthologous gene pairs of the two closely related species to analyze CpG substitution patterns between these two genomes. We also divided genome sequences into three regions: protein-coding, noncoding, and RNA-coding, and obtained the CpG frequencies for each region for each organism. It was found that the observed/expected ratio of CpG dinucleotides is low in both the protein-coding and noncoding regions; while that ratio is in the normal range in the RNA-coding region. Our results indicate that CpG suppression of the Mycoplasma genome is not caused by (1) biased usage amino acid; (2) biased usage of synonymous codon; or (3) methylation effects by the CpG methyltransferase in the genomes of their hosts. Instead, we consider it likely that a certain global pressure, such as genome-wide pressure for the advantages of DNA stability or replication, has the effect of decreasing CpG over the entire genome, which, in turn, resulted in the biased codon usage.     

Comparative Analysis of CpG Islands among HBV Genotypes

DNA methylation is being increasingly recognized to play a role in regulation of hepatitis B virus (HBV) gene expression. The aim of this study was to compare the CpG island distribution among different HBV genotypes. We analyzed 176 full-length HBV genomic sequences obtained from the GenBank database, belonging to genotypes A through J, to identify the CpG islands in the HBV genomes. Our results showed that while 79 out of 176 sequences contained three conventional CpG islands (I–III) as previously described, 83 HBV sequences harbored only two of the three known islands. Novel CpG islands were identified in the remaining 14 HBV isolates and named as CpG island IV, V, and VI. Among the eight known HBV genotypes and two putative genotypes, while HBV genomes containing three CpG islands were predominant in genotypes A, B, D, E, and I; genotypes C, F, G, and H tended to contain only two CpG islands (II and III). In conclusion, the CpG islands, which are potential targets for DNA methylation mediated by the host functions, differ among HBV genotypes, and these genotype-specific differences in CpG island distribution could provide new insights into the understanding of epigenetic regulation of HBV gene expression and hepatitis B disease outcome.     

Segmental distribution of genes harboring a CpG island-like region on rice chromosomes

In plant genomes, there exist discrete regions rich in CpG dinucleotides, namely CpG clusters. In rice, most of these CpG clusters are associated with genes. Rice genes are grouped into one of the five classes according to the position of an associated CpG cluster. Among them, class 1 genes, which harbor a CpG cluster at the 5′-terminus, share similarities with human genes having CpG islands. In the present study, by analyzing plant genome sequence data, primarily from rice, we investigated the chromosomal distribution of genes of each class, mainly class 1 genes. Class 1 genes were not uniformly distributed across the rice genome, but were clustered into discrete chromosomal segments. EST-based analysis of the distribution of expressed genes indicates that this segmental distribution of class 1 genes caused a preferential distribution of expressed genes within class 1 gene-rich segments. We then compared the methylation status of genes of each class to examine the possibility that differential DNA methylation, if any, is relevant to the observed differential expression level of genes inside and outside the class 1 segments. The difference in the methylation level between these genes was revealed to be fairly small, which does not support the above-mentioned possibility. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Prenatal antidepressant exposure associated with CYP2E1 DNA methylation change in neonates

Some but not all neonates are affected by prenatal exposure to serotonin reuptake inhibitor antidepressants (SRI) and maternal mood disturbances. Distinguishing the impact of these 2 exposures is challenging and raises critical questions about whether pharmacological, genetic, or epigenetic factors can explain the spectrum of reported outcomes. Using unbiased DNA methylation array measurements followed by a detailed candidate gene approach, we examined whether prenatal SRI exposure was associated with neonatal DNA methylation changes and whether such changes were associated with differences in birth outcomes. Prenatal SRI exposure was first associated with increased DNA methylation status primarily at CYP2E1(β_Non-exposed = 0.06, β_SRI-exposed = 0.30, FDR = 0); however, this finding could not be distinguished from the potential impact of prenatal maternal depressed mood. Then, using pyrosequencing of CYP2E1 regulatory regions in an expanded cohort, higher DNA methylation status—both the mean across 16 CpG sites (P < 0.01) and at each specific CpG site (P < 0.05)—was associated with exposure to lower 3rd trimester maternal depressed mood symptoms only in the SRI-exposed neonates, indicating a maternal mood x SRI exposure interaction. In addition, higher DNA methylation levels at CpG2 (P = 0.04), CpG9 (P = 0.04) and CpG10 (P = 0.02), in the interrogated CYP2E1 region, were associated with increased birth weight independently of prenatal maternal mood, SRI drug exposure, or gestational age at birth. Prenatal SRI antidepressant exposure and maternal depressed mood were associated with altered neonatal CYP2E1 DNA methylation status, which, in turn, appeared to be associated with birth weight.     

The majority of methylated deoxycytidines in human DNA are not in the CpG dinucleotide

The only natural postsynthetic modification known to occur in mammalian DNA is the methylation in the 5 position of deoxycytidines. Of the four 5'-CpN-3' dinucleotides (ie. CpG, CpC, CpA, and CpT), the dinucleotide which contains the highest proportion of deoxycytidines methylated is CpG, with 40 to 80% methylation in different mammalian genomes. It has also been shown that CpA, CpT, and CpC are methylated as well but to a much lower extent. Here we report the result of a full nearest neighbour analysis (together with quantitation of methylation levels in the 4 CpN dinucleotides) for DNA from human spleen. Using the values we have calculated the overall frequencies for all the methylated dinucleotides in the human genome. Because of the relative underrepresentation (by 7 to 10 fold) of the CpG dinucleotide, only 45.5% of total mC was present in mCpG, with 54.5% in mCpA, mCpT plus mCpC. These calculations have implications for studies into the function and significance of DNA methylation in mammalian cells.