DNA motifs associated with aberrant CpG island methylation

Epigenetic silencing involving the aberrant methylation of promoter region CpG islands is widely recognized as a tumor suppressor silencing mechanism in cancer. However, the molecular pathways underlying aberrant DNA methylation remain elusive. Recently we showed that, on a genome-wide level, CpG island loci differ in their intrinsic susceptibility to aberrant methylation and that this susceptibility can be predicted based on underlying sequence context. These data suggest that there are sequence/structural features that contribute to the protection from or susceptibility to aberrant methylation. Here we use motif elicitation coupled with classification techniques to identify DNA sequence motifs that selectively define methylation-prone or methylation-resistant CpG islands. Motifs common to 28 methylation-prone or 47 methylation-resistant CpG island-containing genomic fragments were determined using the MEME and MAST algorithms (). The five most discriminatory motifs derived from methylation-prone sequences were found to be associated with CpG islands in general and were nonrandomly distributed throughout the genome. In contrast, the eight most discriminatory motifs derived from the methylation-resistant CpG islands were randomly distributed throughout the genome. Interestingly, this latter group tended to associate with Alu and other repetitive sequences. Used together, the frequency of occurrence of these motifs successfully discriminated methylation-prone and methylation-resistant CpG island groups with an accuracy of 87% after 10-fold cross-validation. The motifs identified here are candidate methylation-targeting or methylation-protection DNA sequences.     

Histone methylation marks play important roles in predicting the methylation status of CpG islands

The methylation status of CpG islands is highly correlated with gene expression. Current methods for computational prediction of DNA methylation only utilize DNA sequence features. In this study, besides 35 DNA sequence features, we added four histone methylation marks to predict the methylation status of CpG islands, and improved the accuracy to 89.94%. Also we applied our model to predict the methylation pattern of all the CpG islands in the human genome, and the results are consistent with the previous reports. Our results imply the important roles of histone methylation marks in affecting the methylation status of CpG islands. H3K4me enriched in the methylation-resistant CpG islands could disrupt the contacts between nucleosomes, unravel chromatin and make DNA sequences accessible. And the established open environment may be a prerequisite for or a consequence of the function implementation of zinc finger proteins that could protect CpG islands from DNA methylation.     

一种快速高效的全基因组甲基化CpG岛富集方法的建立

高甲基化的CpG岛所致基因表观遗传学转录失活已经成为肿瘤表观基因组学研究的重要内容。现在已有很多检测CpG岛甲基化的方法,但由于各自的局限,还没有建立一种能快速在全基因组水平上进行甲基化CpG岛的富集方法。本研究利用甲基化结合蛋白MBD2b具有特异性结合甲基化DNA的特性, 建立了一种基于DNA免疫共沉淀技术的全基因组甲基化CpG岛的富集方法。在大肠杆菌中表达重组的GST-MBD2b蛋白,通过Glutathione Sepharose 4B对重组蛋白进行纯化,制备成亲和层析柱,利用在不同的盐离子强度下甲基化DNA和非甲基化DNA的结合能力不同,对甲基化DNA进行富集。用甲基化酶SssI处理过的DNA片段与非甲基化DNA片段进行富集效率的检测,发现0.5M KCl的浓度是甲基化DNA片段和非甲基化DNA片段得以分开的临界条件。样品的富集效率用Real Time PCR进行检测。结果表明,这种方法能够实现对全基因组甲基化DNA的有效富集且最高的富集倍数可达到100多倍。富集到的甲基化DNA可以进行后续的定量PCR, DNA测序和全基因组芯片的分析等工作,为大规模分析全基因组CpG 岛甲基化的改变奠定了基础。     

High density DNA methylation array with single CpG site resolution

We have developed a new generation of genome-wide DNA methylation BeadChip which allows high-throughput methylation profiling of the human genome. The new high density BeadChip can assay over 480K CpG sites and analyze twelve samples in parallel. The innovative content includes coverage of 99% of RefSeq genes with multiple probes per gene, 96% of CpG islands from the UCSC database, CpG island shores and additional content selected from whole-genome bisulfite sequencing data and input from DNA methylation experts. The well-characterized Infinium® Assay is used for analysis of CpG methylation using bisulfite-converted genomic DNA. We applied this technology to analyze DNA methylation in normal and tumor DNA samples and compared results with whole-genome bisulfite sequencing (WGBS) data obtained for the same samples. Highly comparable DNA methylation profiles were generated by the array and sequencing methods (average R² of 0.95). The ability to determine genome-wide methylation patterns will rapidly advance methylation research.     

Germline methylation patterns inferred from local nucleotide frequency of repetitive sequences in the human genome

Given the genomic abundance and susceptibility to DNA methylation, interspersed repetitive sequences in the human genome can be exploited as valuable resources in genome-wide methylation studies. To learn about the relationships between DNA methylation and repeat sequences, we performed a global measurement of CpG dinucleotide frequencies for interspersed repetitive sequences and inferred germline methylation patterns in the human genome. Although extensive CpG depletion was observed for most repeat sequences, those in the proximity to CpG islands have been relatively removed from germline methylation being the potential source of germline activation. We also investigated the CpG depletion patterns of Alu pairs to see whether they might play an active role in germline methylation. Two kinds of Alu pairs, direct or inverted pairs classified according to the orientation, showed contrast CpG depletion patterns with respect to separating distance of Alus, i.e., as two Alu elements are more closely spaced in a pair, a higher extent of CpG depletion was observed in inverted orientation and vice versa for directly repetitive Alu pairs. This suggests that specific organization of repetitive sequences, such as inverted Alu pairs, might play a role in triggering DNA methylation consistent with a homology-dependent methylation hypothesis.     

Systematic assessment of reduced representation bisulfite sequencing to human blood samples: A promising method for large-sample-scale epigenomic studies

Complementary to the time- and cost-intensive direct bisulfite sequencing, we applied reduced representation bisulfite sequencing (RRBS) to the human peripheral blood mononuclear cells (PBMC) from YH, the Asian individual whose genome and epigenome has been deciphered in the YH project and systematically assessed the genomic coverage, coverage depth and reproducibility of this technology as well as the concordance of DNA methylation levels measured by RRBS and direct bisulfite sequencing for the detected CpG sites. Our result suggests that RRBS can cover more than half of CpG islands and promoter regions with a good coverage depth and the proportion of the CpG sites covered by the biological replicates reaches 80-90%, indicating good reproducibility. Given a smaller data quantity, RRBS enjoys much better coverage depth than direct bisulfite sequencing and the concordance of DNA methylation levels between the two methods is high. It can be concluded that RRBS is a time and cost-effective sequencing method for unbiased DNA methylation profiling of CpG islands and promoter regions in a genome-wide scale and it is the method of choice to assay certain genomic regions for multiple samples in a rapid way.     

Genome-Wide Locations of Potential Epimutations Associated with Environmentally Induced Epigenetic Transgenerational Inheritance of Disease Using a Sequential Machine Learning Prediction Approach

Environmentally induced epigenetic transgenerational inheritance of disease and phenotypic variation involves germline transmitted epimutations. The primary epimutations identified involve altered differential DNA methylation regions (DMRs). Different environmental toxicants have been shown to promote exposure (i.e., toxicant) specific signatures of germline epimutations. Analysis of genomic features associated with these epimutations identified low-density CpG regions (<3 CpG / 100bp) termed CpG deserts and a number of unique DNA sequence motifs. The rat genome was annotated for these and additional relevant features. The objective of the current study was to use a machine learning computational approach to predict all potential epimutations in the genome. A number of previously identified sperm epimutations were used as training sets. A novel machine learning approach using a sequential combination of Active Learning and Imbalance Class Learner analysis was developed. The transgenerational sperm epimutation analysis identified approximately 50K individual sites with a 1 kb mean size and 3,233 regions that had a minimum of three adjacent sites with a mean size of 3.5 kb. A select number of the most relevant genomic features were identified with the low density CpG deserts being a critical genomic feature of the features selected. A similar independent analysis with transgenerational somatic cell epimutation training sets identified a smaller number of 1,503 regions of genome-wide predicted sites and differences in genomic feature contributions. The predicted genome-wide germline (sperm) epimutations were found to be distinct from the predicted somatic cell epimutations. Validation of the genome-wide germline predicted sites used two recently identified transgenerational sperm epimutation signature sets from the pesticides dichlorodiphenyltrichloroethane (DDT) and methoxychlor (MXC) exposure lineage F3 generation. Analysis of this positive validation data set showed a 100% prediction accuracy for all the DDT-MXC sperm epimutations. Observations further elucidate the genomic features associated with transgenerational germline epimutations and identify a genome-wide set of potential epimutations that can be used to facilitate identification of epigenetic diagnostics for ancestral environmental exposures and disease susceptibility.     

Copy number variations of 11 macronuclear chromosomes and their gene expression in Oxytricha trifallax

Genome-wide methylation studies frequently lack adequate controls to estimate proportions of background reads in the resulting datasets. To generate appropriate control pools, we developed technique termed nMETR (non-methylated tag recovery) based on digestion of genomic DNA with methylation-sensitive restriction enzyme, ligation of adapter oligonucleotide and PCR amplification of non-methylated sites associated with genomic repetitive elements. The protocol takes only two working days to generate amplicons for deep sequencing. We applied nMETR for human DNA using BspFNI enzyme and retrotransposon Alu-specific primers. 454-sequencing enabled identification of 1113 nMETR tag sites, of them ~65% were parts of CpG islands. Representation of reads inversely correlated with methylation levels, thus confirming nMETR fidelity. We created software that eliminates background reads and enables to map and annotate individual tags on human genome. nMETR tags may serve as the controls for large-scale epigenetic studies and for identifying unmethylated transposable elements located close to genomic CpG islands.     

人类基因组DNA甲基化数据分析的研究现状

作为人类基因组最为典型的表观遗传现象,DNA甲基化在多种关键生理活动中扮演重要角色.系统分析基因组尺度的DNA甲基化概况意义重大.从Cp G岛等基本定义出发,阐述了高通量DNA甲基化的检测技术以及针对芯片技术与下一代测序技术的低水平数据处理方法;重点对比了基于机器学习理论对Cp G位点及Cp G岛甲基化水平的预测算法,以及所利用的特征对预测效果的影响与发展趋势;并对DNA差异甲基化在组织特异性、癌症等多种疾病中的计算分析进行了全面的综述.     

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.     

Regionally Specific and Genome-Wide Analyses Conclusively Demonstrate the Absence of CpG Methylation in Human Mitochondrial DNA

Although CpG methylation clearly distributes genome-wide in vertebrate nuclear DNA, the state of methylation in the vertebrate mitochondrial genome has been unclear. Several recent reports using immunoprecipitation, mass spectrometry, and enzyme-linked immunosorbent assay methods concluded that human mitochondrial DNA (mtDNA) has much more than the 2 to 5% CpG methylation previously estimated. However, these methods do not provide information as to the sites or frequency of methylation at each CpG site. Here, we have used the more definitive bisulfite genomic sequencing method to examine CpG methylation in HCT116 human cells and primary human cells to independently answer these two questions. We found no evidence of CpG methylation at a biologically significant level in these regions of the human mitochondrial genome. Furthermore, unbiased next-generation sequencing of sodium bisulfite treated total DNA from HCT116 cells and analysis of genome-wide sodium bisulfite sequencing data sets from several other DNA sources confirmed this absence of CpG methylation in mtDNA. Based on our findings using regionally specific and genome-wide approaches with multiple human cell sources, we can definitively conclude that CpG methylation is absent in mtDNA. It is highly unlikely that CpG methylation plays any role in direct control of mitochondrial function.     

Improved reduced representation bisulfite sequencing for epigenomic profiling of clinical samples

Background

DNA methylation plays crucial roles in epigenetic gene regulation in normal development and disease pathogenesis. Efficient and accurate quantification of DNA methylation at single base resolution can greatly advance the knowledge of disease mechanisms and be used to identify potential biomarkers. We developed an improved pipeline based on reduced representation bisulfite sequencing (RRBS) for cost-effective genome-wide quantification of DNA methylation at single base resolution. A selection of two restriction enzymes (Taq^αI and MspI) enables a more unbiased coverage of genomic regions of different CpG densities. We further developed a highly automated software package to analyze bisulfite sequencing results from the Solexa GAIIx system.

Results

With two sequencing lanes, we were able to quantify ~1.8 million individual CpG sites at a minimum sequencing depth of 10. Overall, about 76.7% of CpG islands, 54.9% of CpG island shores and 52.2% of core promoters in the human genome were covered with at least 3 CpG sites per region.

Conclusions

With this new pipeline, it is now possible to perform whole-genome DNA methylation analysis at single base resolution for a large number of samples for understanding how DNA methylation and its changes are involved in development, differentiation, and disease pathogenesis.     

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.     

CpG island methylation in human lymphocytes is highly correlated with DNA sequence, repeats, and predicted DNA structure

CpG island methylation plays an important role in epigenetic gene control during mammalian development and is frequently altered in disease situations such as cancer. The majority of CpG islands is normally unmethylated, but a sizeable fraction is prone to become methylated in various cell types and pathological situations. The goal of this study is to show that a computational epigenetics approach can discriminate between CpG islands that are prone to methylation from those that remain unmethylated. We develop a bioinformatics scoring and prediction method on the basis of a set of 1,184 DNA attributes, which refer to sequence, repeats, predicted structure, CpG islands, genes, predicted binding sites, conservation, and single nucleotide polymorphisms. These attributes are scored on 132 CpG islands across the entire human Chromosome 21, whose methylation status was previously established for normal human lymphocytes. Our results show that three groups of DNA attributes, namely certain sequence patterns, specific DNA repeats, and a particular DNA structure, are each highly correlated with CpG island methylation (correlation coefficients of 0.64, 0.66, and 0.49, respectively). We predicted, and subsequently experimentally examined 12 CpG islands from human Chromosome 21 with unknown methylation patterns and found more than 90% of our predictions to be correct. In addition, we applied our prediction method to analyzing Human Epigenome Project methylation data on human Chromosome 6 and again observed high prediction accuracy. In summary, our results suggest that DNA composition of CpG islands (sequence, repeats, and structure) plays a significant role in predisposing CpG islands for DNA methylation. This finding may have a strong impact on our understanding of changes in CpG island methylation in development and disease.     

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.