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.     

Gene-associated CpG islands in plants as revealed by analyses of genomic sequences

We screened plant genome sequences, primarily from rice and Arabidopsis thaliana, for CpG islands, and identified DNA segments rich in CpG dinucleotides within these sequences. These CpG-rich clusters appeared in the analysed sequences as discrete peaks and occurred at the frequencies of one per 4.7 kb in rice and one per 4.0 kb in A. thaliana. In rice and A. thaliana, most of the CpG-rich clusters were associated with genes, which suggests that these clusters are useful landmarks in genome sequences for identifying genes in plants with small genomes. In contrast, in plants with larger genomes, only a few of the clusters were associated with genes. These plant CpG-rich clusters satisfied the criteria used for identifying human CpG islands, which suggests that these CpG clusters may be regarded as plant CpG islands. The position of each island relative to the 5'-end of its associated gene varied considerably. Genes in the analysed sequences were grouped into five classes according to the position of the CpG islands within their associated genes. A large proportion of the genes belonged to one of two classes, in which a CpG island occurred near the 5'-end of the gene or covered the whole gene region. The position of a plant CpG island within its associated gene appeared to be related to the extent of tissue-specific expression of the gene; the CpG islands of most of the widely expressed rice genes occurred near the 5'-end of the genes.     

Reconstructing the ancestral germ line methylation state of young repeats

One of the key objectives of comparative genomics is the characterization of the forces that shape genomes over the course of evolution. In the last decades, evidence has been accumulated that for vertebrate genomes also epigenetic modifications have to be considered in this context. Especially, the elevated mutation frequency of 5-methylcytosine (5mC) is assumed to facilitate the depletion of CpG dinucleotides in species that exhibit global DNA methylation. For instance, the underrepresentation of CpG dinucleotides in many mammalian genomes is attributed to this effect, which is only neutralized in so-called CpG islands (CGIs) that are preferentially unmethylated and thus partially protected from rapid CpG decay. For primate-specific CpG-rich transposable elements from the ALU family, it is unclear whether their elevated CpG frequency is caused by their small age or by the absence of DNA methylation. In consequence, these elements are often misclassified in CGI annotations. We present a method for the estimation of germ line methylation from pairwise ancestral-descendant alignments. The approach is validated in a simulation study and tested on DNA repeats from the AluSx family. We conclude that a predicted unmethylated state in the germ line is highly correlated with epigenetic activity of the respective genomic region. Thus, CpG-rich repeats can be facilitated as in silico probes for the epigenetic potential of their genomic neighborhood.     

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.     

Tissue-specific methylation of a CpG island in transgenic mice.

Clustering of CpG dinucleotides in CpG-rich islands is a characteristic feature of mammalian genomes. Such CpG islands are frequently associated with genes and usually hypomethylated, regardless of the gene activity. This is the case for the CpG island of the murine Thy-1 gene. A transgenic line containing multiple copies of a truncated, concatemeric CpG island from the Thy-1.1 allele (Thy-1.2 background) showed that a stable fraction (approx. 0.20) became fully methylated in somatic tissues of homozygous mice with respect to testable restriction sites, while the remaining copies were methylation-free, i.e., this methylation appears to be an 'all-or-none' phenomenon. DNA from extraembryonic tissues (placenta and yolk sac) and epididymal sperm showed, however, an even higher degree of methylation in two distinct patterns. In the extraembryonic tissue, partial methylation of each copy was seen, whereas in sperm a high degree of 'all-or-none' methylation (greater than 0.35) was observed.     

High-throughput assay of DNA methylation based on methylation-specific primer and SAGE

Mapping of genomic DNA methylation is a dispensable part of functional genome. We have developed a novel method based on methylation-specific primer and serial analysis of gene expression, called MSP-SAGE, with potential of high-throughput quantification of genomic DNA methylation. We used a 6-mer methylation-specific primer to extend the methylated CpG sequences other than non-methylated CpG sequences. The 17 bp tags contained methylated CpG sequence, which were obtained from extended methylation sequence by digestion of restriction endonuclease, and then the tags were concatenated and cloned for sequencing. We can identify the locations of methylation according to the sequences of tags and quantify the methylation status from the frequency of the tags. MSP-SAGE has a good linearity in a broad methylation range from 5% to 100% with good accuracy and high precision. The proof-of-principle study shows that MSP-SAGE is a reliable high-throughput assay for quantification of DNA methylation.     

Bio-CAP: a versatile and highly sensitive technique to purify and characterise regions of non-methylated DNA

Across vertebrate genomes methylation of cytosine residues within the context of CpG dinucleotides is a pervasive epigenetic mark that can impact gene expression and has been implicated in various developmental and disease-associated processes. Several biochemical approaches exist to profile DNA methylation, but recently an alternative approach based on profiling non-methylated CpGs was developed. This technique, called CxxC affinity purification (CAP), uses a ZF-CxxC (CxxC) domain to specifically capture DNA containing clusters of non-methylated CpGs. Here we describe a new CAP approach, called biotinylated CAP (Bio-CAP), which eliminates the requirement for specialized equipment while dramatically improving and simplifying the CxxC-based DNA affinity purification. Importantly, this approach isolates non-methylated DNA in a manner that is directly proportional to the density of non-methylated CpGs, and discriminates non-methylated CpGs from both methylated and hydroxymethylated CpGs. Unlike conventional CAP, Bio-CAP can be applied to nanogram quantities of genomic DNA and in a magnetic format is amenable to efficient parallel processing of samples. Furthermore, Bio-CAP can be applied to genome-wide profiling of non-methylated DNA with relatively small amounts of input material. Therefore, Bio-CAP is a simple and streamlined approach for characterizing regions of the non-methylated DNA, whether at specific target regions or genome wide.     

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.     

Extensive methylation of a part of the CpG island located 3.0-4.5 kbp upstream to the chicken alpha-globin gene cluster may contribute to silencing the globin genes in non-erythroid cells

Here, we show that in the chicken genome, the domain of alpha-globin genes is preceded by a CpG island of which the downstream part ( approximately 0.65 kbp) is heavily methylated in lymphoid cells; it is either non-methylated or undermethylated in erythroid cells. Recombinant plasmids were constructed with the corresponding DNA fragment (called "uCpG") placed upstream to a reporter CAT gene expressed from the promoter of the alpha(D) chicken globin gene. Selective methylation of CpG dinucleotides within the uCpG fragment suppressed fivefold the expression of the CAT gene, when neither this gene itself nor the alpha(D) promoter were methylated. Methylation of CpG dinucleotides within the alpha(D) gene promoter did not modify the suppression effect exerted by methylated uCpG. We interpret these results within the frame of the hypothesis postulating, that methylation of the upstream CpG island of the chicken alpha-globin gene domain may play an essential role in silencing the alpha-globin genes in non-erythroid cells.     

Heterologous CpG island becomes extensively methylated in the genome of transgenic mice

It is demonstrated that a heterologous (chicken) CpG island containing five Sp1 canonical recognition sequences becomes highly methylated in the genome of transgenic mice bearing one or several copies of the transgene. Similar levels of methylation of the chicken CpG island were observed in different tissues of transgenic mice except the brain where the level of methylation of this chicken CpG-rich fragment was significantly lower than in other tissues. Analysis of susceptibility of the "transgenic" CpG island to Hpa II and Msp I restriction nucleases revealed an unusual methylation pattern interfering with the action of both of these enzymes. A conclusion has been drawn that heterologous CpG island per se does not contain all necessary signals permitting to maintain its own non-methylated status in the genome of transgenic animals.     

Epigenetic Natural Variation in Arabidopsis thaliana

Cytosine methylation of repetitive sequences is widespread in plant genomes, occurring in both symmetric (CpG and CpNpG) as well as asymmetric sequence contexts. We used the methylation-dependent restriction enzyme McrBC to profile methylated DNA using tiling microarrays of Arabidopsis Chromosome 4 in two distinct ecotypes, Columbia and Landsberg erecta. We also used comparative genome hybridization to profile copy number polymorphisms. Repeated sequences and transposable elements (TEs), especially long terminal repeat retrotransposons, are densely methylated, but one third of genes also have low but detectable methylation in their transcribed regions. While TEs are almost always methylated, genic methylation is highly polymorphic, with half of all methylated genes being methylated in only one of the two ecotypes. A survey of loci in 96 Arabidopsis accessions revealed a similar degree of methylation polymorphism. Within-gene methylation is heritable, but is lost at a high frequency in segregating F₂ families. Promoter methylation is rare, and gene expression is not generally affected by differences in DNA methylation. Small interfering RNA are preferentially associated with methylated TEs, but not with methylated genes, indicating that most genic methylation is not guided by small interfering RNA. This may account for the instability of gene methylation, if occasional failure of maintenance methylation cannot be restored by other means.     

A nuclear protein that binds preferentially to methylated DNA in vitro may play a role in the inaccessibility of methylated CpGs in mammalian nuclei

The effects of DNA methylation on gene expression and chromatin structure suggest the existence of a mechanism in the nucleus capable of distinguishing methylated and non-methylated sequences. We report the finding of a nuclear protein in several vertebrate tissues and cell lines that binds preferentially to methylated DNA in vitro. Its lack of sequence-specific requirements makes it potentially capable of binding to any methylated sequence in mammalian nuclei. An in vivo counterpart of these results is that methylated CpGs are inaccessible to nucleases within nuclei. In contrast, non-methylated CpG sites, located mainly at CpG islands, and restriction sites not containing this dinucleotide, are relatively accessible. The possibility that DNA methylation acts through binding to specific proteins that could alter chromatin structure is discussed.     

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.