Inter-individual variation of DNA methylation and its implications for large-scale epigenome mapping

Genomic DNA methylation profiles exhibit substantial variation within the human population, with important functional implications for gene regulation. So far little is known about the characteristics and determinants of DNA methylation variation among healthy individuals. We performed bioinformatic analysis of high-resolution methylation profiles from multiple individuals, uncovering complex patterns of inter-individual variation that are strongly correlated with the local DNA sequence. CpG-rich regions exhibit low and relatively similar levels of DNA methylation in all individuals, but the sequential order of the (few) methylated among the (many) unmethylated CpGs differs randomly across individuals. In contrast, CpG-poor regions exhibit substantially elevated levels of inter-individual variation, but also significant conservation of specific DNA methylation patterns between unrelated individuals. This observation has important implications for experimental analysis of DNA methylation, e.g. in the context of epigenome projects. First, DNA methylation mapping at single-CpG resolution is expected to uncover informative DNA methylation patterns for the CpG-poor bulk of the human genome. Second, for CpG-rich regions it will be sufficient to measure average methylation levels rather than assaying every single CpG. We substantiate these conclusions by an in silico benchmarking study of six widely used methods for DNA methylation mapping. Based on our findings, we propose a cost-optimized two-track strategy for mammalian methylome projects.     

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.     

DNA Methylation Analysis of Chromosome 21 Gene Promoters at Single Base Pair and Single Allele Resolution

Differential DNA methylation is an essential epigenetic signal for gene regulation, development, and disease processes. We mapped DNA methylation patterns of 190 gene promoter regions on chromosome 21 using bisulfite conversion and subclone sequencing in five human cell types. A total of 28,626 subclones were sequenced at high accuracy using (long-read) Sanger sequencing resulting in the measurement of the DNA methylation state of 580427 CpG sites. Our results show that average DNA methylation levels are distributed bimodally with enrichment of highly methylated and unmethylated sequences, both for amplicons and individual subclones, which represent single alleles from individual cells. Within CpG-rich sequences, DNA methylation was found to be anti-correlated with CpG dinucleotide density and GC content, and methylated CpGs are more likely to be flanked by AT-rich sequences. We observed over-representation of CpG sites in distances of 9, 18, and 27 bps in highly methylated amplicons. However, DNA sequence alone is not sufficient to predict an amplicon's DNA methylation status, since 43% of all amplicons are differentially methylated between the cell types studied here. DNA methylation in promoter regions is strongly correlated with the absence of gene expression and low levels of activating epigenetic marks like H3K4 methylation and H3K9 and K14 acetylation. Utilizing the single base pair and single allele resolution of our data, we found that i) amplicons from different parts of a CpG island frequently differ in their DNA methylation level, ii) methylation levels of individual cells in one tissue are very similar, and iii) methylation patterns follow a relaxed site-specific distribution. Furthermore, iv) we identified three cases of allele-specific DNA methylation on chromosome 21. Our data shed new light on the nature of methylation patterns in human cells, the sequence dependence of DNA methylation, and its function as epigenetic signal in gene regulation. Further, we illustrate genotype–epigenotype interactions by showing novel examples of allele-specific methylation.     

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.     

The molecular basis for the lack of immunostimulatory activity of vertebrate DNA

Macrophages and B cells are activated by unmethylated CpG-containing sequences in bacterial DNA. The lack of activity of self DNA has generally been attributed to CpG suppression and methylation, although the role of methylation is in doubt. The frequency of CpG in the mouse genome is 12.5% of Escherichia coli, with unmethylated CpG occurring at approximately 3% the frequency of E. coli. This suppression of CpG alone is insufficient to explain the inactivity of self DNA; vertebrate DNA was inactive at 100 micro g/ml, 3000 times the concentration at which E. coli DNA activity was observed. We sought to resolve why self DNA does not activate macrophages. Known active CpG motifs occurred in the mouse genome at 18% of random occurrence, similar to general CpG suppression. To examine the contribution of methylation, genomic DNAs were PCR amplified. Removal of methylation from the mouse genome revealed activity that was 23-fold lower than E. coli DNA, although there is only a 7-fold lower frequency of known active CpG motifs in the mouse genome. This discrepancy may be explained by G-rich sequences such as GGAGGGG, which potently inhibited activation and are found in greater frequency in the mouse than the E. coli genome. In summary, general CpG suppression, CpG methylation, inhibitory motifs, and saturable DNA uptake combined to explain the inactivity of self DNA. The immunostimulatory activity of DNA is determined by the frequency of unmethylated stimulatory sequences within an individual DNA strand and the ratio of stimulatory to inhibitory sequences.     

CpG_MI: a novel approach for identifying functional CpG islands in mammalian genomes

CpG islands (CGIs) are CpG-rich regions compared to CpG-depleted bulk DNA of mammalian genomes and are generally regarded as the epigenetic regulatory regions in association with unmethylation, promoter activity and histone modifications. Accurate identification of CpG islands with epigenetic regulatory function in bulk genomes is of wide interest. Here, the common features of functional CGIs are identified using an average mutual information method to differentiate functional CGIs from the remaining CGIs. A new approach (CpG mutual information, CpG_MI) was further explored to identify functional CGIs based on the cumulative mutual information of physical distances between two neighboring CpGs. Compared to current approaches, CpG_MI achieved the highest prediction accuracy. This approach also identified new functional CGIs overlapping with gene promoter regions which were missed by other algorithms. Nearly all CGIs identified by CpG_MI overlapped with histone modification marks. CpG_MI could also be used to identify potential functional CGIs in other mammalian genomes, as the CpG dinucleotide contents and cumulative mutual information distributions are almost the same among six mammalian genomes in our analysis. It is a reliable quantitative tool for the identification of functional CGIs from bulk genomes and helps in understanding the relationships between genomic functional elements and epigenomic modifications.     

Methylation status of CpG-rich islands on active and inactive mouse X chromosomes

Single copy probes derived from CpG-rich island clones fromEag I andNot I linking libraries and nine rare-cutter restriction endonucleases were used to investigate the methylation status of CpG-rich islands on the inactive and active X chromosomes (Chr) of the mouse. Thirteen of the 14 probes used detected CpG-rich islands in genomic DNA. The majority of island CpGs detected by rare-cutter restriction endonucleases were methylated on the inactive X Chr and unmethylated on the active X Chr, but some heterogeneity within the cell population used to make genomic DNA was detected. The CpG-rich islands detected by two putative pseudoautosomal probes remained unmethylated on both the active and inactive X Chrs. Otherwise, distance from the X Chr inactivation center did not affect the methylation profile of CpG-rich islands. We conclude that methylation of CpG-rich islands is a general feature of X Chr inactivation.     

De novo DNA methylation at nonrandom founder sites 5' from an unmethylated minimal origin of DNA replication in latent Epstein-Barr virus genomes

Latent episomal genomes of Epstein-Barr virus, a human gammaherpesvirus, represent a suitable model system for studying replication and methylation of chromosomal DNA in mammals. We analyzed the methylation patterns of CpG dinucleotides in the latent origin of DNA replication of Epstein-Barr virus using automated fluorescent genomic sequencing of bisulfite-modified DNA samples. We observed that the minimal origin of DNA replication was unmethylated in 8 well-characterized human cell lines or clones carrying latent Epstein-Barr virus genomes as well as in a prototype virus producer marmoset cell line. This observation suggests that unmethylated DNA domains can function as initiation sites or zones of DNA replication in human cells. Furthermore, 5' from this unmethylated region we observed focal points of de novo DNA methylation in nonrandom positions in the majority of Burkitt's lymphoma cell lines and clones studied while the corresponding CpG dinucleotides in viral genomes carried by lymphoblastoid cell lines and marmoset cells were completely unmethylated. Clustering of highly methylated CpG dinucleotides suggests that de novo methylation of unmethylated double-stranded episomal viral genomes starts at discrete founder sites in vivo. This is the first comparative high-resolution methylation analysis of a latent viral origin of DNA replication in human cells.     

Alu repeated DNAs are differentially methylated in primate germ cells.

A significant fraction of Alu repeats in human sperm DNA, previously found to be unmethylated, is nearly completely methylated in DNA from many somatic tissues. A similar fraction of unmethylated Alus is observed here in sperm DNA from rhesus monkey. However, Alus are almost completely methylated at the restriction sites tested in monkey follicular oocyte DNA. The Alu methylation patterns in mature male and female monkey germ cells are consistent with Alu methylation in human germ cell tumors. Alu sequences are hypomethylated in seminoma DNAs and more methylated in a human ovarian dysgerminoma. These results contrast with methylation patterns reported for germ cell single-copy, CpG island, satellite, and L1 sequences. The function of Alu repeats is not known, but differential methylation of Alu repeats in the male and female germ lines suggests that they may serve as markers for genomic imprinting or in maintaining differences in male and female meiosis.     

A sensitive new method for rapid detection of abnormal methylation patterns in global DNA and within CpG islands.

To assess alterations in DNA methylation density in both global DNA and within CpG islands, we have developed a simple method based on the use of methylation-sensitive restriction endonucleases that leave a 5' guanine overhang after DNA cleavage, with subsequent single nucleotide extension with radiolabeled [(3)H]dCTP. The methylation-sensitive restriction enzymes HpaII and AciI have relatively frequent recognition sequences at CpG sites that occur randomly throughout the genome. BssHII is a methylation sensitive enzyme that similarly leaves a guanine overhang, but the recognition sequence is nonrandom and occurs predominantly at unmethylated CpG sites within CpG islands. The selective use of these enzymes can be used to screen for alterations in genome-wide methylation and CpG island methylation status, respectively. The extent of [(3)H]dCTP incorporation opposite the exposed guanine after restriction enzyme treatment is directly proportional to the number of unmethylated (cleaved) CpG sites. The "cytosine-extension assay" has several advantages over existing methods because (a) radiolabel incorporation is independent of the integrity of the DNA, (b) methylation detection does not require PCR amplification or DNA methylase reactions, and (c) it is applicable to ng quantities of DNA. Using DNA extracted from normal human liver and from human hepatocellular carcinoma, the applicability of the assay is demonstrated by the detection of an increase in genome-wide hypomethylation and CpG island hypermethylation in the tumor DNA.     

Enzymatic properties of recombinant Dnmt3a DNA methyltransferase from mouse: the enzyme modifies DNA in a non-processive manner and also methylates non-CpG [correction of non-CpA] sites

We present the first in vitro study investigating the catalytic properties of a mammalian de novo DNA methyltransferase. Dnmt3a from mouse was cloned and expressed in Escherichia coli. It was shown to be catalytically active in E. coli cells in vivo. The methylation activity of the purified protein was highest at pH 7.0 and 30 mM KCl. Our data show that recombinant Dnmt3a protein is indeed a de novo methyltransferase, as it catalyzes the transfer of methyl groups to unmethylated substrates with similar efficiency as to hemimethylated substrates. With oligonucleotide substrates, the catalytic activity of Dnmt3a is similar to that of Dnmt1: the K(m) values for the unmethylated and hemimethylated oligonucleotide substrates are 2.5 microM, and the k(cat) values are 0.05 h(-1) and 0.07 h(-1), respectively. The enzyme catalyzes the methylation of DNA in a distributive manner, suggesting that Dnmt3a and Dnmt1 may cooperate during de novo methylation of DNA. Further, we investigated the methylation activity of Dnmt3a at non-canonical sites. Even though the enzyme shows maximum activity at CpG sites, with oligonucleotide substrates, a high methylation activity was also found at CpA sites, which are modified only twofold slower than CpG sites. Therefore, the specificity of Dnmt3a is completely different from that of the maintenance methyltransferase Dnmt1, which shows a 40 to 50-fold preference for hemimethylated over unmethylated CpG sites and has almost no methylation activity at non-CpG sites.