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.     

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     

The distribution of the dinucleotide CpG and cytosine methylation in the vitellogenin gene family

Summary Sequence data from regions of five vertebrate vitellogenin genes were used to examine the frequency, distribution, and mutability of the dinucleotide CpG, the preferred modification site for eukaryotic DNA methyltransferases. The observed level of the CpG dinucleotide in all five genes was markedly lower than that expected from the known mononucleotide frequencies. CpG suppression was greater in introns than in exons. CpG-containing codons were found to be avoided in the vitellogenin genes, but not completely despite the redundancy of the genetic code. Frequency and distribution patterns of this dinucleotide varied dramatically among these otherwise closely related genes. Dense clusters of CpG dinucleotides tended to appear in regions of either functional or structural interest (e.g., in the transposon-like Vi-element ofXenopus) and these clusters contained 5-methylcytosine (5 mC). 5 mC is known to undergo deamination to form thymidine, but the extent to which this transition occurs in the heavily methylated genomes of vertebrates and its contribution to CpG suppression are still unclear. Sequence comparison of the methylated vitellogenin gene regions identified CT and GA substitutions that were found to occur at relatively high frequencies. The predicted products of CpG deamination, TpG and CpA, were elevated. These findings are consistent with the view that CpG distribution and methylation are interdependent and that deamination of 5 mC plays an important role in promoting evolutionary change at the nucleotide sequence level.     

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.     

The covariation between TpA deficiency, CpG deficiency, and G+C content of human isochores is due to a mathematical artifact

CpG and TpA dinucleotides are underrepresented in the human genome. The CpG deficiency is due to the high mutation rate from C to T in methylated CpG's. The TpA suppression was thought to reflect a counterselection against TpA's destabilizing effect in RNA. Unexpectedly, the TpA and CpG deficiencies vary according to the G+C contents of sequences. It has been proposed that the variation in CpG suppression was correlated with a particular chromatin organization in G+C-rich isochores. Here, we present an improved model of dinucleotide evolution accounting for the overlap between successive dinucleotides. We show that an increased mutation rate from CpG to TpG or CpA induces both an apparent TpA deficiency and a correlation between CpG and TpA deficiencies and G+C content. Moreover, this model shows that the ratio of observed over expected CpG frequency underestimates the real CpG deficiency in G+C-rich sequences. The predictions of our model fit well with observed frequencies in human genomic data. This study suggests that previously published selectionist interpretations of patterns of dinucleotide frequencies should be taken with caution. Moreover, we propose new criteria to identify unmethylated CpG islands taking into account this bias in the measure of CpG depletion.     

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.     

CpG deficiency, dinucleotide distributions and nucleosome positioning

The dinucleotide CpG is deficient in (A + T)-rich regions of vertebrate DNA in both coding and non-coding sequences and there is a corresponding increase above expectation in the occurrence of TpG and CpA. By contrast in (G + C)-rich regions no deficiency of CpG is found. Such (G + C)-rich sequences, containing the expected number of CpG dinucleotides, alternate along the genome with (A + T)-rich sequences which have a lower than expected CpG content. The G + C content of vertebrate DNA can oscillate with a period of 150-200 bp and this may be a factor in positioning nucleosomes. The role of mutagenesis in loss of CpG and increase of A + T, particularly in non-coding regions, is discussed.     

Active mammalian replication origins are associated with a high-density cluster of mCpG dinucleotides.

ori-beta is a well-characterized origin of bidirectional replication (OBR) located approximately 17 kb downstream of the dihydrofolate reductase gene in hamster cell chromosomes. The approximately 2-kb region of ori-beta that exhibits greatest replication initiation activity also contains 12 potential methylation sites in the form of CpG dinucleotides. To ascertain whether DNA methylation might play a role at mammalian replication origins, the methylation status of these sites was examined with bisulfite to chemically distinguish cytosine (C) from 5-methylcytosine (mC). All of the CpGs were methylated, and nine of them were located within 356 bp flanking the minimal OBR, creating a high-density cluster of mCpGs that was approximately 10 times greater than average for human DNA. However, the previously reported densely methylated island in which all cytosines were methylated regardless of their dinucleotide composition was not detected and appeared to be an experimental artifact. A second OBR, located at the 5' end of the RPS14 gene, exhibited a strikingly similar methylation pattern, and the organization of CpG dinucleotides at other mammalian origins revealed the potential for high-density CpG methylation. Moreover, analysis of bromodeoxyuridine-labeled nascent DNA confirmed that active replication origins were methylated. These results suggest that a high-density cluster of mCpG dinucleotides may play a role in either the establishment or the regulation of mammalian replication origins.     

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.     

Factors to preserve CpG-rich sequences in methylated CpG islands

Background

Mammalian CpG islands (CGIs) normally escape DNA methylation in all adult tissues and developmental stages. However, in our previous study we unexpectedly identified many methylated CGIs in human peripheral blood leukocytes. Methylated CpG dinucleotides convert to TpG dinucleotides through deaminization of their cytosine bases more frequently than hypomethylated CpG dinucleotides. Therefore, we wondered how methylated CGIs in germline or non-germline cells maintain their CpG-rich sequences. It is known that events such as germline hypomethylation, CpG selection, biased gene conversion (BGC), and frequent CpG fixation can contribute to the maintenance of CpG-rich sequences in methylated CGIs in germline or non-germline cells. However, it has not been investigated which of the processes maintain CpG-rich sequences of methylated CGIs in each genomic position.

Results

In this study, we comprehensively examined the contribution of the processes described above to the maintenance of CpG-rich sequences in methylated CGIs in germline and non-germline cells which were classified by genomic positions. Approximately 60–80% of CGIs with high methylation in H1 cell line (H1-HM) in all the genomic positions showed a low average CpG → TpG/CpA substitution rate. In contrast, fewer than half the numbers of CGIs with H1-HM in all the genomic positions showed a low average CpG → TpG/CpA substitution rate and low levels of methylation in sperm cells (SPM-LM). Furthermore, a small fraction of CGIs with a low average CpG → TpG/CpA substitution rate and high levels of methylation in sperm cells (SPM-HM) showed CpG selection.On the other hand, independent of the positions in genes, most CGIs with SPM-HM showed a slightly higher average TpG/CpA → CpG substitution rate compared with those with SPM-LM.

Conclusions

Relatively high numbers (approximately 60–80%) of CGIs with H1-HM in all the genomic positions preserve their CpG-rich sequences by a low CpG → TpG/CpA substitution rate caused mainly by their SPM-LM, and for those with SPM-HM partly by CpG selection and TpG/CpA → CpG fixation. BGC has little contribution to the maintenance of CpG-rich sequences of CGIs with SPM-HM which were classified by genomic positions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1286-x) contains supplementary material, which is available to authorized users.     

Methylation of CpG dinucleotides in the lacI gene of the Big Blue® transgenic mouse

Cytosine residues at CpG dinucleotides can be methylated by endogenous methyltransferases in mammalian cells. The resulting 5-methylcytosine base may undergo spontaneous deamination to form thymine causing G/C to A/T transition mutations. Methylated CpGs also can form preferential targets for environmental mutagens and carcinogens. The Big Blue® transgenic mouse has been used to investigate tissue and organ specificity of mutations and to deduce mutational mechanisms in a mammal in vivo. The transgenic mouse contains approximately 40 concatenated lambda-like shuttle vectors, each of which contains one copy of an Escherichia coli lacI gene as a mutational target. lacI mutations in lambda transgenic mice are characterized by a high frequency of spontaneous mutations targeted to CpG dinucleotides suggesting an important contribution from methylation-mediated events. To study the methylation status of CpGs in the lacI gene, we have mapped the distribution of 5-methylcytosines along the DNA-binding domain and flanking sequences of the lacI gene of transgenic mice. We analyzed genomic DNA from various tissues including thymus, liver, testis, and DNA derived from two thymic lymphomas. The mouse genomic DNAs and methylated and unmethylated control DNAs were chemically cleaved, then the positions of 5-methylcytosines were mapped by ligation-mediated PCR which can be used to distinguish methylated from unmethylated cytosines. Our data show that most CpG dinucleotides in the DNA binding domain of the lacI gene are methylated to a high extent (>98%) in all tissues tested; only a few sites are partially (70–90%) methylated. We conclude that tissue-specific methylation is unlikely to contribute significantly to tissue-specific mutational patterns, and that the occurrence of common mutation sites at specific CpGs in the lacI gene is not related to selective methylation of only these sequences. The data confirm previous suggestions that the high frequency of CpG mutations in lacI transgenes is related to the presence of 5-methylcytosine bases.     

The loss of CpC dinucleotides from DNA. II. Methylated and non-methylated genes of vertebrates

The frequencies of neighboring b.p. in more than 1100 genes of vertebrates in the EMBL bank (1000 kb) have been analysed. It has been found that the majority of such genes exhibit a lack of CpG duplexes and an excess of TpG+CpA. The loss of CpG may indicate that the major part of these sites in the genome is methylated and has been subjected to the pressure of CpG----TpG+CpA mutations. The methylated genes grouped into compartment M+ are represented by a fraction of repeated sequences and by genes of the most rapidly diverging families of proteins (globins, immunoglobulins, structural proteins, etc.). The genes of this compartment are characterized by a correlation between the G+C content and the value of CpG-suppression. A group of genes has been detected in which the CpG mutation process has gone so far that nearly all of these dinucleotides have disappeared from DNA. Judging by the value of CpG-suppression, these genes, grouped in the Mo+ compartment, used to be strongly methylated before. However, in the now extant vertebrates they have fully depleted their CpG reserve and for this reason lost the methylation capacity. Transitions in methylated CpG may be one of the sources of spontaneous mutagenesis resulting in the enhanced genetic instability of the cell. A gene compartment has been detected with an intermediate level of CpG deficiency; this compartment has been designated as M+. In these genes only a few of the available CpGs have been steadily methylated (and subjected to mutation). It has been found that the genome of vertebrates contains a specific CpG-rich fraction which exhibits no CpG-suppression, irrespective of the overall content of G+C. Probably, CpG sites have persisted unmethylated throughout the existence of these genes. We suggest them to constitute a M- compartment. This compartment comprises the genes of tRNA and rRNA (5S, 5.8S, 18S, 28S) and small nuclear RNAs U2-U6, as well as the genes of core histones, some enzymes, viruses and 5'-flanking sequences of certain protein-coding genes. In the genome of vertebrates, the genes of the evolutionary most conserved proteins and RNAs have not undergone methylation. A list of genes, belonging to different compartments of the vertebrate genome, is given. Compartment Mo+ constitutes 19%, M(+)--35%, M(+/-)--28% and M(-)--8% of all the vertebrate genes studied. Possible mechanisms, protecting the functionally most significant genes of vertebrates from methylation, and discussed.     

Possible implications of CpG avoidance in the flatworm Schistosoma mansoni

We report the analysis of the biases of CpG, TpG, and CpA of all the DNA sequences data from the Trematode Schistosoma mansoni. Our results show CpG avoidance whereas TpG and CpA frequencies are over the expected values. These characteristics are similar to the biases displayed by methylated genomes, but in platyhelminths 5mC has not been detected by biochemical methods. The possible implications of this CpG shortage are discussed.     

Homology-dependent inactivation of LTR retrotransposons in <Emphasis Type="Italic">Aspergillus fumigatus</Emphasis> and <Emphasis Type="Italic">A. nidulans</Emphasis> genomes

The repeat-induced point mutation mechanism (RIP) is the most intriguing among the known mechanisms of homology-dependent gene inactivation (silencing) because of its ability to produce irreversible mutations in repetitive DNA sequences. Discovered for the first time in Neurospora crassa, RIP is characterized by C:G to T:A transitions in duplicated sequences. The mechanisms and range of occurrence of RIP are still poorly understood. Mobile elements, including retrotransposons, are a common target for the processes that lead to homology-dependent silencing because of their ability to propagate themselves. Comparative analysis of LTR retrotransposons was performed throughout the genomes of two ascomycetes, Aspergillus fumigatus and A. nidulans. “De-RIP” retroelements were reconstructed on the basis of several copies. CpG, CpA, and TpG sites, which are potential targets for mutagenesis, were found at a much lower frequency in mobile elements than in structural genes. The dinucleotide targets of the two species are affected by RIP at different frequencies: mutagenesis occurs at both CpG and CpA sites in A. fumigatus and is confined to CpG dinucleotides in A. nidulans. This work provides a theoretical background for planning the experimental investigation of RIP inactivation in aspergilli.     

Sequence context at human single nucleotide polymorphisms: overrepresentation of CpG dinucleotide at polymorphic sites and suppression of variation in CpG islands

Human polymorphisms originate as mutations, and the influence of context on mutagenesis should be reflected in the distribution of sequences surrounding single nucleotide polymorphisms (SNPs). We have performed a computational survey of nearly two million human SNPs to determine if sequence-dependent hotspots for polymorphism exist in the human genome. Here we show that sequences containing CpG dinucleotides, which occur at low frequencies in the human genome, are 6.7-fold more abundant at polymorphic sites than expected. In contrast, polymorphisms in CpG sequences located within CpG islands, important regulatory regions that modulate gene expression, are 6.8-fold less prevalent than expected. The distribution of polymorphic alleles at CpGs in CpG islands is also significantly different from that in non-island regions. These data strongly support a role for 5-methylcytosine deamination in the generation of human variation, and suggest that variation at CpGs in islands is suppressed.     

Frequency of abnormal human haemoglobins caused by C----T transitions in CpG dinucleotides

A large part of human genetic disease apparently arises from deamination of cytosines in methylated CpG dinucleotides. Their mutation rate is known to be high when C is present as 5-methylcytosine, but is believed to be normal when it is unmethylated. The beta-globin gene contains five, the gamma-globin gene two, and each of the alpha-globin genes contain 35 CpGs. The CpGs in the beta- and gamma-globin genes are methylated, while those in the alpha-globin genes are undermethylated. One would therefore have expected the CpGs to be a frequent source of mutations in the beta- and gamma-globin genes, but not in the alpha-globin genes. In fact, the evidence points to CpGs being a frequent source of mutations in both the alpha- and beta-globin genes. This suggests either that the mutation rates of both methylated and unmethylated CpGs are abnormally high, which conflicts with published evidence, or that there is a finite chance of some CpGs in the alpha-globin genes of certain individuals being methylated and therefore subject to mutation.     

Methylation, mutation and cancer.

The fifth base in human DNA, 5-methylcytosine, is inherently mutagenic. This has led to marked changes in the distribution of the CpG methyl acceptor site and an 80% depletion in its frequency of occurrence in vertebrate DNA. The coding regions of many genes contain CpGs which are methylated in sperm and serve as hot spots for mutation in human genetic diseases. Fully 30-40% of all human germline point mutations are thought to be methylation induced even though the CpG dinucleotide is under-represented and efficient cellular repair systems exist. Importantly, tumor suppressor genes such as p53 also contain methylated CpGs and these serve as hot spots for mutations in some, but not all, human cancers. Comparison of the spectrum of mutations present in this gene in different human cancers allows for predictions to be made on the molecular mechanisms of tumorigenesis.     

The loss of dinucleotides CpG from DNA. IV. Methylation and divergence of genes and pseudogenes of small nuclear RNA

The frequency of neighboring base pairs in nucleotide sequences of over 80 genes and pseudogenes of low molecular weight RNAs U1-U8, 4.5S and 7S in different eukaryotes was determined. The probable frequency of CpG----TpG + CpA substitutions, caused as a result of the deamination of 5-methylcytosine residues in DNA, was determined. It was found that the genes of small RNAs do not reveal a single level of CpG methylation for all the species studied. In most cases CpG in the genes of U1, 4.5S and 7S RNA are methylated, whereas in the genes of U2-U6-RNA these sites must have never been subjected to methylation. Nearly all the investigated pseudogenes of different small RNAs are strongly methylated due to a considerable lack of CpG. It was established that CpG----TpG + CpA transitions may amount to as much third of all the mutations accumulated in the genes of the same RNAs in different species. Such transitions in pseudogenes may account for 40% of all the nucleotide substitutions. This disproportionately high level of mutations in CpG dinucleotides (3-5-fold higher than in other DNA dupletes) must be the direct result of the methylation of these sites. Consequently, CpG methylation causes a dramatic acceleration of the divergence rate of DNA sequences. It has been concluded that protection of most vital genes against methylation is one of the essential conditions for sustaining the high level of stability of the macromolecular structure and for the reliability of macromolecular functioning in a cell.