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.     

DNA methylation variation in cloned mice

Summary: Mammalian cloning has been accomplished in several mammalian species by nuclear transfer. However, the production rate of cloned animals is quite low, and many cloned offspring die or show abnormal symptoms. A possible cause of the low success rate of cloning and abnormal symptoms in many cloned animals is the incomplete reestablishment of DNA methylation after nuclear transfer. We first analyzed tissue‐specific methylation patterns in the placenta, skin, and kidney of normal B6D2F1 mice. There were seven spots/CpG islands (0.5% of the total CpG islands detected) methylated differently in the three different tissues examined. In the placenta and skin of two cloned fetuses, a total of four CpG islands were aberrantly methylated or unmethylated. Interestingly, three of these four loci corresponded to the tissue‐specific loci in the normal control fetuses. The extent of aberrant methylation of genomic DNA varied between the cloned animals. In cloned animals, aberrant methylation occurred mainly at tissue‐specific methylated loci. Individual cloned animals have different methylation aberrations. In other words, cloned animals are by no means perfect copies of the original animals as far as the methylation status of genomic DNA is concerned. genesis 30:45–50, 2001. © 2001 Wiley‐Liss, Inc.     

CpG mutation rates in the human genome are highly dependent on local GC content

CpG dinucleotides mutate at a high rate because cytosine is vulnerable to deamination, cytosines in CpG dinucleotides are often methylated, and deamination of 5-methylcytosine (5mC) produces thymidine. Previous experiments have shown that DNA melting is the rate-limiting step in cytosine deamination. Here we show, through the analysis of human single-nucleotide polymorphisms (SNPs), that the mutation rate produced by 5mC deamination is highly dependent on local GC content. In fact, linear regression analysis showed that the log(10) of the 5mC mutation rates (inferred from SNP frequencies) had slopes of -3 when graphed with respect to the GC content of neighboring sequences. This is the ideal slope that would be expected if the correlation between CpG underrepresentation and GC content had been solely caused by DNA melting. Moreover, this same result was obtained regardless of the SNP locations (all SNPs versus only SNPs in noncoding intergenic regions, excluding CpG islands) and regardless of the lengths over which GC content was calculated (SNP sequences with a modal length of 564 bp versus genomic contigs with a modal length of 163 kb). Several alternative interpretations are discussed.     

Synthesis of universal unmethylated control DNA by nested whole genome amplification with phi29 DNA polymerase

Optimization of highly sensitive methods to detect methylation of CpG islands in gene promoter regions requires adequate methylated and unmethylated control DNA. Whereas universal methylated control DNA is available, universal unmethylated control (UUC) DNA has not been made because demethylase is not available to remove methyl groups from all methylated cytosines. On the basis that DNA synthesized by DNA polymerase does not contain methylated cytosines, we developed a method to create UUC DNA by nested whole genome amplification (WGA) with phi29 DNA polymerase. Contamination of the template genomic DNA in UUC was only 3.1 x 10(-7), below the detection limit of sensitive methods used for methylation studies such as methylation-specific PCR. Assessment of microsatellite markers demonstrated that even nested phi29 WGA achieves highly accurate and homogeneous amplification with very low amounts of genomic DNA as an initial template. The UUC DNA created by nested phi29 WGA is practically very useful for methylation analysis.     

Genome-wide analysis of DNA methylation status of CpG islands in embryoid bodies, teratomas, and fetuses

Differentiation of embryonic stem (ES) cells into embryoid bodies (EBs) provides an in vitro system for the study of early lineage determination during mammalian development. We have previously reported that there are 247 CpG islands that potentially have tissue-dependent and differentially methylated regions (T-DMRs). This provided evidence that the formation of DNA methylation patterns at CpG islands is a crucial epigenetic event underlying mammalian development. Here we present an analysis by the restriction landmark genomic scanning (RLGS) using NotI as a landmark enzyme of the genome-wide methylation status of CpG islands of ES cells and EBs and of teratomas produced from ES cells. These results are considered in relation to the methylation status of CpG islands of genomic DNA from normal fetus (10.5 dpc) and adult tissues. We have prepared a DNA methylation panel that consists of 259 T-DMRs and includes novel T-DMRs that are distinctly methylated or unmethylated in the teratomas. The DNA methylation pattern was complex and differed for the ES cells, EBs, and teratomas, providing evidence that differentiation of cells involves both de novo DNA methylation as well as demethylation. Comparison of the numbers of T-DMRs, that were differentially methylated or unmethylated among the cells and tissue types studied, revealed that the teratomas were the most epigenetically different from ES cells. Thus, analysis of the DNA methylation profiles prepared in this study provides new insights into the differentiation of ES cells and development of fetus, EB, teratoma, and somatic tissues.     

Linking the epigenome to the genome: correlation of different features to DNA methylation of CpG islands

DNA methylation of CpG islands plays a crucial role in the regulation of gene expression. More than half of all human promoters contain CpG islands with a tissue-specific methylation pattern in differentiated cells. Still today, the whole process of how DNA methyltransferases determine which region should be methylated is not completely revealed. There are many hypotheses of which genomic features are correlated to the epigenome that have not yet been evaluated. Furthermore, many explorative approaches of measuring DNA methylation are limited to a subset of the genome and thus, cannot be employed, e.g., for genome-wide biomarker prediction methods. In this study, we evaluated the correlation of genetic, epigenetic and hypothesis-driven features to DNA methylation of CpG islands. To this end, various binary classifiers were trained and evaluated by cross-validation on a dataset comprising DNA methylation data for 190 CpG islands in HEPG2, HEK293, fibroblasts and leukocytes. We achieved an accuracy of up to 91% with an MCC of 0.8 using ten-fold cross-validation and ten repetitions. With these models, we extended the existing dataset to the whole genome and thus, predicted the methylation landscape for the given cell types. The method used for these predictions is also validated on another external whole-genome dataset. Our results reveal features correlated to DNA methylation and confirm or disprove various hypotheses of DNA methylation related features. This study confirms correlations between DNA methylation and histone modifications, DNA structure, DNA sequence, genomic attributes and CpG island properties. Furthermore, the method has been validated on a genome-wide dataset from the ENCODE consortium. The developed software, as well as the predicted datasets and a web-service to compare methylation states of CpG islands are available at http://www.cogsys.cs.uni-tuebingen.de/software/dna-methylation/.     

Analysis of repetitive element DNA methylation by MethyLight

Repetitive elements represent a large portion of the human genome and contain much of the CpG methylation found in normal human postnatal somatic tissues. Loss of DNA methylation in these sequences might account for most of the global hypomethylation that characterizes a large percentage of human cancers that have been studied. There is widespread interest in correlating the genomic 5-methylcytosine content with clinical outcome, dietary history, lifestyle, etc. However, a high-throughput, accurate and easily accessible technique that can be applied even to paraffin-embedded tissue DNA is not yet available. Here, we report the development of quantitative MethyLight assays to determine the levels of methylated and unmethylated repeats, namely, Alu and LINE-1 sequences and the centromeric satellite alpha (Satα) and juxtacentromeric satellite 2 (Sat2) DNA sequences. Methylation levels of Alu, Sat2 and LINE-1 repeats were significantly associated with global DNA methylation, as measured by high performance liquid chromatography, and the combined measurements of Alu and Sat2 methylation were highly correlative with global DNA methylation measurements. These MethyLight assays rely only on real-time PCR and provide surrogate markers for global DNA methylation analysis. We also describe a novel design strategy for the development of methylation-independent MethyLight control reactions based on Alu sequences depleted of CpG dinucleotides by evolutionary deamination on one strand. We show that one such Alu-based reaction provides a greatly improved detection of DNA for normalization in MethyLight applications and is less susceptible to normalization errors caused by cancer-associated aneuploidy and copy number changes.     

Densely methylated sequences that are preferentially localized at telomere-proximal regions of human chromosomes

We have constructed a library of densely methylated DNA sequences from human blood DNA by selecting fragments with a high affinity for a methyl-CpG binding domain (MBD) column. PCR analysis of the library confirmed the presence of known densely methylated CpG island sequences. Analysis of random clones, however, showed that the library was dominated by sequences whose G+C content and CpG frequency were intermediate between those of bulk genomic DNA and bona fide CpG islands. When human chromosomes were probed with the library by fluorescent in situ hybridisation (FISH), the predominant sites of labelling were at terminal regions of many chromosomes, approximately corresponding to T-bands. Analysis of the methylation status of random clones indicated that all were heavily methylated at CpGs in blood DNA, but many were under-methylated in sperm DNA. Lack of methylation in germ cells may reduce CpG depletion at some sub-terminal sequences and result in a high density of methyl-CpG when these regions become methylated in somatic cells.     

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.     

Effect of CpG Island Methylation on MicroRNA Expression in the k-562 Cell Line

To test the hypothesis that methylation of a CpG island is associated with regulation of microRNA expression, we investigated CpG islands in the upstream sequences of microRNA precursors (pre-miRNAs) through bioinformatic analysis and determined whether the CpG islands were methylated by methylation-specific PCR in the k-562 cell line. We used 5-azacytidine for DNA demethylation, and changes in microRNA expression were detected by microarray assay, RT-PCR, and real-time PCR after 5-azacytidine induction. We showed that the CpG islands in the upstream regions of 18 pre-miRNAs were methylated, including miR-663, miR-369, miR-615, and miR-410, and promoter activity was detected in the upstream region of pre-miR-663. We found that a decrease in methylation of a CpG island could up-regulate the expression of miR-663, suggesting that miR-663 could be regulated by DNA methylation. Expression levels of miR-369, miR-615, and miR-410 were not regulated by DNA methylation in this cell line.     

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     

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.     

Methylation and sequence analysis around EagI sites: identification of 28 new CpG islands in XQ24-XQ28.

Thirty-two probes for CpG islands of the distal long arm of the human X chromosome have been identified. From a genomic library of DNA of the hamster-human cell hybrid X3000.1 digested with the rare cutter restriction enzyme EagI, 53 different human clones have been isolated and characterized by methylation and sequence analysis. The characteristic pattern of DNA methylation of CpG islands at the 5' end of genes of the X chromosome has been used to distinguish between EagI sites in CpG islands versus isolated EagI sites. The sequence analysis has confirmed and completed the characterization showing that sequences at the 5' end of known genes were among the clones defined CpG islands and that the non-CpG islands clones were mostly repetitive sequences with a non-methylated or variably methylated EagI site. Thus, since clones corresponding to repetitive sequences can be easily identified by sequencing, such libraries are a very good source of CpG islands. The methylation analysis of 28 different new probes allows to state that demethylation of CpG islands of the active X and methylation of those on the inactive X chromosome are the general rule. Moreover, the finding, in all instances, of methylation differences between male and female DNA is in very strong support of the notion that most genes of the distal long arm of the X chromosome are subject to X inactivation.