DNA methylation of distal regulatory sites characterizes dysregulation of cancer genes
© Aran et al.; licensee BioMed Central Ltd. 2013
Received: 23 August 2012
Accepted: 12 March 2013
Published: 12 March 2013
Abnormal epigenetic marking is well documented in gene promoters of cancer cells, but the study of distal regulatory siteshas lagged behind.We performed a systematic analysis of DNA methylation sites connected with gene expression profilesacross normal and cancerous human genomes.
Utilizing methylation and expression data in 58 cell types, we developed a model for methylation-expression relationships in gene promoters and extrapolated it to the genome. We mapped numerous sites at which DNA methylation was associated with expression of distal genes. These sites bind transcription factors in a methylation-dependent manner, and carry the chromatin marks of a particular class of transcriptional enhancers. In contrast to the traditional model of one enhancer site per cell type, we found that single enhancer sites may define gradients of expression levels across many different cell types. Strikingly, the identified sites were drastically altered in cancers: hypomethylated enhancer sites associated with upregulation of cancer-related genes and hypermethylated sites with downregulation. Moreover, the association between enhancer methylation and gene deregulation in cancerwas significantly stronger than the association of promoter methylationwith gene deregulation.
Methylation of distal regulatory sites is closely related to gene expression levels across the genome. Single enhancers may modulate ranges of cell-specific transcription levels, from constantlyopen promoters. In contrast to the remote relationships between promoter methylation and gene dysregulation in cancer, altered methylation of enhancer sites is closely related to gene expression profiles of transformed cells.
KeywordsCancer DNA methylation distal control elements enhancers epigenomics gene-enhancer pairing gene regulation machine-learning
DNA methylation is a key determinant of regulatory chromatin complexes transmitted through cell divisions. Relationships between DNA methylation and gene expression levels were first recognized at gene promoters and CpG islands , but were recently observed across the genome [2–6]. In gene promoters, DNA methylation mediates silencing by altering the binding of transcription modulators to their DNA targets , but the cause and function of expression-associated methylation away from gene promoters has been elusive. The relationships between aberrant DNA methylation and the altered expression profiles of cancer cells are also not well understood: the predominant pattern in gene promoters is denovo methylation of polycomb-repressed genes [8–11]. Since these genes are already inactive in the normal tissue and generally remain inactive in the cancer [12, 13], it is hard to establish direct contribution to the cancer process . Away from gene promoters, perturbed methylation has been associated with reduced genomic stability [15, 16] and global silencing of large chromatin domains [17, 18], but the effect on transcription of particular genes remains unknown.
Transcriptional enhancers support tissue-specific expression profiles through physical interactions with gene promoters. Based on analyses of chromatin structures, hundreds ofthousands of transcriptional enhancers were predicted in mammalian genomes (reviewed in ). These sites bind chromatin-modulating factors, interact with distal promoters through DNA loops, and demonstrate a unique pattern of DNA methylation . In some particular examples, enhancer sites turned out to be less methylated in cells expressing the controlled genes than in non-expressing cells [21–24]. Taking a more systematic approach, Wienchet al. recently revealed activity-dependent methylation in a group of distal enhancers in the mouse , and Bock et al. demonstrated correlations between enhancer methylation and expression of particular developmental genes during mouse tissue differentiation.
We conducted a methodical analysis of distal DNA methylation sites associated with gene expression in normal and cancerous human cells. The results suggest that enhancer methylation corresponds closely with expression profiles of cancer genes in transformed cells.
Results and discussion
To explore relationships between DNA methylation and gene expression levels across the genome, we took the following strategy. First, we developed a model for the relationships between methylation in gene promoters and gene expression using a machine-learning algorithm. Then, we applied the derived model to methylation sites away from gene promoters and characterized the chromatin structure and binding profile of the discovered sites. Finally, we explored how the identified sites were altered in cancer.
We analyzed available DNA methylation and gene expression data for 58 human cell types (Table S1 in Additional file 1). The methylation data were produced by two different assays: reduced representation bisulfite sequencing (RRBS) andInfinium HumanMethylation450 BeadChip(Illuminainc., San Diego, CA, USA). Following verification of high agreement between the assays (Figures S1A and S2A, Bin Additional file 1), we combined the two datasets. We then analyzed variation in methylation levels across the panel of cell types (see Materials and methods for detailed description).Sites with methylation levels that did not change across the cell types were eliminated, and the remaining 670, 906 variable methylation sites (VMSs) were used in the study (Figure S1Bin Additional file 1). VMSs in gene promoterswere negatively correlated with expression whereas VMSs in gene bodies were positively correlated with expression (Figure S2 in Additional file 1), as previously reported . The connection with expression levels was most evident at -500 to +2, 000 base pairs relative to transcription start sites (TSSs) (Figure S2A-Cin Additional file 1). Therefore, we used the CpG methylation sites in this range for model development.
Promoter-based model of regulatory methylation sites
Relating genes to distal regulatory sitesusing DNA methylation
We then sorted out the sites within 5 kb of promoters or alternative promoters of the associated genes. The remaining 1,911 pairs included 486 unique genes corresponding with 1,041 unique distalmethylation regions (we arbitrarily clustered CpGs with <3 kb between them into a single region). On average, we obtained 2.14distal regions and 3.93 methylation sites per associated gene. The following analyses were done on these 1,911 high-scoring CpG-gene pairs.
To confirm that the high-scoring methylation sites wereindeed regulatory sites, we analyzed their chromatin signatures. High-scoring sites were enriched by histone modifications that have been previously associated with transcription activation, including histone 3 lysine 4 methylations, lysine 27 acetylation, lysine 9 acetylation and histone variant H2A.Z, but not with CTCF binding (Figure 2C). Further analysis of chromatin states defined by the ChromHMM algorithm (based on various chromatin factors ) revealed significant enrichment of the high-scoring sites within promoters and within a particular class (ChromHMM state 4) of strong enhancers (P = 2.4e-15; Figure 2D). In agreement with this, VMSs were hypomethylated when included in enhancer chromatin, compared to the methylation of the same sites in cells where they were included in non-regulatory chromatin (Figure S3 in Additional file 1). In addition, the highest frequency of hypomethylated sites was observed in state 4 enhancers (Figure S3Cin Additional file 1). We also analyzed the protein-binding capacities of the identified distal sites. The mapped methylation sites bind a larger number of transcription factors than expected. Moreover, unmethylated sites bind more factors than methylated ones (Figure 2E).Finally, we found that the high-scoring methylation sites are evolutionarily conserved (Figure 2F). Together, these findings suggest that distal expression-related methylation sites populate a particular class of transcriptional enhancers, which bind transcription factors in a methylation-depended manner.
Validation of enhancer-promoter pairing
To validate distal enhancer-promoter interactions, we compared our predictions with a recent comprehensive study of long-range DNA interactions, measured bythe chromosome conformation capture carbon copy (5C) technique in three cell types . We identified enhancer sites (n = 318)assessed by both methylation score and 5Cprobes, and determined whether the methylation-based gene-enhancer pairs also reveal physical DNA interactions. Between-assay agreements were markedly increased with methylation scores:at scores ≥0.85, 5C interactions were found in 53% of the cases, compared with 29% of lower score interactions (P = 0.021). A representative example of a long-range enhancer-promoter interaction captured by the both assay is shown in Figure S4 in Additional file 1. We also analyzed the frequency of DNA interactions in 21 enhancer-TSS pairs predicted by our assay, versus the interactions of the same enhancers with the other 264 TSSs in the regions (10 kb to 1 Mb from the enhancer site). The frequency of 5C interactions was significantly higher with the TSS predicted by our assay than with the other TSSs (308.5 ±738.2 junction sequence reads in the high-scoring pairs, versus 90.7 ±277.7 reads in the other pairs, P = 0.009). Thus, thelarge-scale looping analysis affirms the fidelity of enhancer-promoter pairings based on DNA methylation signals. Given that our study examined many more cell types than the three analyzed in the 5C study, and thus could identify additional enhancers not active in these three cell types, this level of between-assay agreement is probably an underestimation of the actual degree of valid pairs.
Since certain epigenetic profiles may not be fully preserved in cell lines , we also confirmed the ability to map enhancer-promoter interactions in freshlyobtained (uncultured) tissue samples. Analysis of seven tissue biopsies from two donors ([GSE:30654]) revealed significant intersection (P <0.01) with the enhancer-gene pairs observed in the original set of 58 cell types (see Materials and methods for detailed description). Thus, a significant number of the gene-enhancer interactions predicted in cell cultures were replicated in uncultured tissue samples. This analysis also indicates a significant degree of overlapping between gene-enhancer pairing in diverse tissue collections.
Enhancer methylation defines gradients of cell-type transcription levels from permissive promoters
Altered enhancer methylation predicts changes in the expression profiles of cancer genes
To examine whether these patterns are exclusive to mammary cells, we repeated the analysis in lung epithelia. The analysis of normal human bronchial epithelial cells versus epithelial lung carcinoma (A549) cells revealed similar relationships between enhancer and promoter methylation and gene expression (Figure 4A, middle row). We have further analyzed a set of 18 normal cell types of various tissue origins versus 18 cancers of various tissue origins (Table S2 in Additional file 1). Once again, the methylation of enhancers was drastically altered and was tightly connected with alteration of gene expression, whereas promoters showed smaller alterations and weaker correlation with expression (R = -0.51 versus R = -0.24); P <1e-20; Figure 4A, lower row). Thus, alteration of enhancer methylation is common to many cancer types, and is associated with substantial changes in gene expression.
Since our enhancers were specifically enriched in the chromatin class defined as ChromHMM state 4 (Figure 2C), and since enhancers in this state showed the highest methylation differences compared to non-regulatory chromatin (Figure S3 in Additional file 1), we hypothesized that state-4 enhancers would show the closest association with changes in gene expression in cancers. Indeed, this group of enhancers showed an even tighter correlation with gene expression levels compared with the entire group of enhancers (R = -0.58 versus R = -0.51; Figure 4B).
We also sought to analyze what characteristics the altered genes might share in addition to association with enhancer methylation. For this, we first analyzed whether the same genes tended to be upregulated or downregulated across cancer types, and then asked whether they belonged to defined gene-ontology (GO) categories. As shown in Figure 4, a majority of the genes were hypermethylated and downregulated, and a smaller fraction of the genes showed the opposite pattern. We found no significant overlap between the genes that were hypermethylated and downregulated in mammary, lung or in the sets of various cell types. These genes also showed no significant GO clustering. Thus, the hypermethylated genes might be cancer-type-specific (for example, cell-type-specific genes that are downregulated in the derived cancer). By contrast, the hypomethylated,upregulated genes significantly repeated between cancer types (Figure 4C, left), and 58 of them appeared in all cancer types examined (P <e-20, compared with the null expectation of 2.3 overlapping genes; Table S3 in Additional file 1). GO analysis showed that the hypomethylated and upregulated genes were frequently involved in cell growing-related functions (Figure 4C right; an extended list of significant GO clusters is given in Table S4 in Additional file 1). Thus, the hypomethylated enhancers associate with genes that are active in many cancer types and support cell proliferation.
An interesting question raised from the above is whether enhancer methylation may also distinguish between cancer and normally-dividing cells. To examine this, we compared Epstein-Barr virus-immortalizedlymphoblastoid(GM12878)cells, which are rapidly growing in culture but otherwise have normal karyotype and no particular characteristics of cancer cells, with lymphoblastoid (Jukat) cells derived from an acute leukemia.We identified 74 genes that were hypomethylated and upregulated in the leukemia compared to the normal lymphoblastoids. As expected, these were not proliferative genes - the proliferative genes were already active in the dividing normal cells, so they were not expected to show differences compared to the cancer. Instead, this analysis captured genes that were involved in cancer-related processes other than cellproliferation, such as the PRAME gene, which is believed to inhibit myeloid differentiation in certain myeloid leukemias (Table S5 in Additional file 1). Of the 74 upregulated genes, 31 were also hypomethylated and upregulated in the sets of normal versus cancer cell types (Figure 4), a significantly higher number than expected by chance (P = 7e-14). Thus, enhancer hypomethylation may associate with upregulation of genes involved in a variety of cancer-related pathways.
Enhancer methylation maybe specifically altered in cancer
With these baseline analyses in hand, we then asked whether any particular group of methylation sites behaved as expected from the global trends, or stand alone. Figure 5F,G shows the average alteration of a given genomic site in cancer, as a function of its methylation level in the normal cells. For example, if a given site was 10% methylated in the normal cells, its average methylation shift in the cancer, based on the analyses shown in Figure 5A,B, was about plus 30% methylation. We found that promoters (green lines in Figure 5E,G) were behaving more or less as expected from the global trend of the genome (black lines). In other words, the alteration of promoter methylation in cancer was similar to the alteration of non-promoter sites with similar normal methylation levels. Looking more carefully into the promoter classes, it appeared that active promoters tended to maintain low methylation against the global trend, whereas poised (polycomb-repressed) promoters gained more methylation than expected from the global trend. Yet, the overall alteration of promoter methylation in cancer was relatively close to that expected from the global trend. The shift in the methylationof the high-scoring enhancers, however, was clearly beyond what could be expected from the global trend (Figure 5F,G). Thus, the alteration of the high-scoring enhancer sites in cancer is probably not simply the outcome of the global loss of methylation control in cancer cells.Rather, these sites may be under the effect of a more specific mechanism, such as targeted methylation or demethylation and/or cell selection.
The mapping of associationsbetween distal regulatory sites and the genes they control is a challenging task, which only recently began to be confronted on the genome-wide scale. Attempts to predict gene-enhancer pairs were based on the profiling of chromatin states or transcription factor binding [27, 34], or of long-range DNA looping [28, 29].Here, we show that enhancers can be also associated with genes using DNA methylation. In contrast to the above mapping approaches, methylation data are readily available and are highly quantitative, and thus may enhance mapping of gene-enhancer pairing.
We found that distal expression-related methylation sites are abundant in the human genome, co-localizing with enhancer chromatin marks, and are more predictive of expression levels then promoter methylations. While not all distal regulatory sites in the genome must exhibited promoter-like methylation, we showed that a large number of enhancer sites demonstrate reverse correlation between methylation and expression, as in gene promoters.
We have further shown that hypomethylation state is directly related to enhancer activity across cell types (Figure S3 in Additional file 1). The observation that hypomethylated enhancers bind more transcription factors than methylated ones (Figure 2E) suggests a possible mechanism underlying the connection between DNA methylation and enhancer activity. Consistent with this possibility, high-scoring enhancers are particularly enriched within a defined chromatin state (Figures 2D), which is particularly hypomethylated compared to non-regulatory chromatin (Figure S3Cin Additional file 1). Whether this chromatin state holds particularly active enhancers, or perhaps a unique class of methylation-related enhancers, remains to be elucidated.
The range of cell types we analyzed in this study was determined by the availability of methylation and expression data. In addition, the RRBS and Infinium HumanMethylation450 BeadChip methylation data we used provide limited genomic coverage and are strongly biased towards promoters and certain other portions of the genome, while enhancers are not efficiently targeted. Because of this, it is likely that more complete methylomic coverage will expose many additional enhancer-gene pairs. Whole-genome bisulfite sequencing approaches have recently become popular and whole methylome analyses of human tissues are rapidly accumulating. Utilizing our approach, these additional data should allow the production of more comprehensive maps of enhancer-gene pairing across tissues, cell types and conditions.
Unmethylated promoters are permissive for, but do not necessarily determine, transcription initiation. We showed that enhancer methylation associates with cell-type-specific expression levels, even when the promoter is constantly unmethylated (Figure 3). Moreover, enhancer methylation characterizes small (and larger) expression differences. Thus, enhancers are not just on-off switches of cell-type transcription levels, as previously suggested, but may also mediate ranges of expression levels across multiple cell types. In contrast to the traditional model of one enhancer site per cell type, we suggest that a gradient of methylation states at a single enhancer site may direct distinct expression levels in many different cell types (Figure 3C).
In occasional examples, enhancer methylation level has been suggested to be associated with the control of cancer-related genes[35–37]. However, to our knowledge this is the first report on a global association between perturbed enhancer methylation and aberrant expression of cancer genes. We have shown that hypomethylated enhancers associated with the upregulation of many cancer genes controlling various cellular functions (Figure 4C), some of them involved in cell proliferation and some in other cancer-related processes. Moreover, many of these hypomethylated genes were found in most cancer types examined, suggesting a pan-cancer mechanism. However, the larger group of hypermethylated enhancers seemed to target cancer-type-specific genes. Given the limited genomic coverage of this study, many additional cancer-related enhancers are expected in the genome.
To date, almost all studies of cancer-related methylation have focused on gene promoters and CpG islands. Among these, the predominant event in cancers is hypermethylation of polycomb-repressed promoters [9–11]. This hypermethylation does not directly affect expression levels, as the associated genes are inactive in the normal tissue and generally remain inactive in the cancer(although it may limit the potential for re-activation of silenced genes in the cancer). Here, we established a very different occurrence in the other large group of regulatory sites - the transcriptional enhancers. These sites are drastically altered in cancers, to both hypo- and hypermethylation, and are closely related to substantial modifications in the expression levels of cancer-related genes (Figure 4). Moreover, their aberrant methylation in cancers might derive from targeted methylation or demethylation or from selection of the altered cells (Figure 5). Whether targeted or selected, aberrant enhancer methylation may be involved in important events during cancer development.
We have uncovered a class of distal methylation sites that closely describe cell-type transcription levels. These sites reside in a particular subclass of transcriptional enhancers and are associated with cell-type-specific enhancer activity, possibly through communication with the binding of transcription factors. Methylation levels of these enhancers associate with gradual expression differences across cell types, even when the linked promoters are consistently unmethylated across the cell types. The radical changes in methylation of these sites in cancer is beyond that expected from the general profile of the cancer methylome, and may reflect specific targeting of the methylation and demethylation machinery to these sites, and/or functional contribution to tumor development. Further analyses of these sites may provide crucial information about paradigms of gene expression control in normal and cancerous cells.
Materials and methods
DNA methylation data for 58 cells types, for which expression data were also available. Methylation data were the average of two RRBS experiments (52 of the cell types, as detailed in Table S1 in Additional file 1) or Infinium HumanMethylation450 BeadChip data (36 of the RRBS cell types plus additional six cell types for which only BeadChip data were available).
Expression data for the above cell types produced by the Affymetrix Human Exon array. Gene-level expression data for 17,862 genes was extracted from the raw data following Robust Multi-chip Average normalization.
Histone chromatin immunoprecipitation (ChIP)-sequencing peaks and ChromHMM annotation data available for six of the above cell types, including K562, human mammary epithelial cells, human skeletal muscle myoblasts, normal human lung fibroblasts, GM12878 and HepG2 cells.
ChIP-sequencing peaks data for 122 DNA binding factors in various cell types.
Human whole-genome methylation data were from Lister et al. . Genomic locations are according to human genome version hg19 of the human genome. Gene sizes, start and end sites, promoters and alternative promoters are according to the RefSeq database (National Center for Biotechnology Information, Bethesda, MD, USA). Methylation and expression data for somatic tissues of two individuals are from ([GSE:30654] patients #1 and 2).
Data filtering and categorization
VMSs were CpGs with standard deviation >5 after omitting outliers utilizing the Thomson's Tau method (alpha = 0.1%). A 'sample' (gene-CpG pair) was the two-dimensional matrix describing the methylation levels of a given CpG site as a function of the expression levels of a given gene across the cell types. For the RRBS data,only methylation sites that were sequenced at least 10times in at least 60% (80% for the Support Vector Machine (SVM) algorithm) of the cell types were included.
The following fourfeatures were extracted for each CpG-gene sample: positive Pearson coefficient (0 if negative);negative Pearson coefficient (0 if positive); positive Spearman coefficient (0 if negative); and negative Spearman coefficient (0 if positive). Based on these features, we learned a model distinguishing false from true samples, applying a machine-learning algorithm. The SVMmap algorithm is a SVM algorithm for predicting rankings. It performs supervised learning using binary labeled training examples, with the goal of optimizing mean average precision (MAP) . We used SVM-MAP for ranking the CpG-gene samples in 3,434 gene queries (RRBS data in 52 cell types) containing 15,205 true samples and a 50-fold excess of randomly selected false samples, where the goal wasfitting a model that ranks the true samples above the false samples. The regularization parameter was optimized by dividing the gene queries into training and test sets (4:1 ratio) and performing 10-fold cross-validation. Query P -values (the probability of obtaining the top ranked true sample by chance) were calculated for the test set queries as follows:
The resulting model for the methylation-expression relationship consisted of the following learned features: positive Pearson coefficient = -0.0058, negative Pearson coefficient = 0.0104, positive Spearman coefficient = -0.0038, and negative Spearman coefficient = 0.0088.
Applying the learned model to the datasets
Frequency map of high-scoring sites
The distribution of methylation-expression relationships of the high-scoring sites over the cell types is a contour representation of the model learned to distinguish true from random pairs. The frequency map shown in Figure 2B is a representation of the corresponding scatter plot (SVM-MAP for 25 × 25 grids of high-scoring pairs), with smoothing of each point around a sphere (radius = 3).
Enrichment of chromatin marks around the high-scoring sites
ChIP-seq data ('modification score') were downloaded from the ENCODE website, averaged across the analyzed cell types, normalized to a 0 to 1 scale, and averaged across the high-scoring methylation sites. Loess smoothing with a span of 10% was applied to the presented data.
Enrichment of chromatin states around the high-scoring sites
ChromHMM states were identified for the CpG sites across the cell types. (A state was called when a given site was found in a given state in at least one of the cell types. Sites may be related to more than one state). The actual number of sites in a given chromatin category and a given score level ('real'), was divided by the number obtained from shuffling the expression data between genes 10 times ('random'). P -values were calculated based on the normal distribution of the shuffled data.
Binding of transcription factors to the high-scoring sites
ChIP-sequencing peaks data for 122 DNA binding factors in fourcell types(GM12878, HepG2, HeLaS3, K562 were downloaded. For every high-scoring site we counted the number of factors that bind, and averaged across the cell types. Random expectation and real versus random P-values were obtained as above. The P -value for the difference between methylated and unmethylated sites was calculated using the Wilcoxon rank-sum test for difference between averages.
Conserved sequences around the high-scoring sites
PhastConsdata (phastCons46wayPrimates and phastCons46wayPlacental) were downloaded from the University of California, Santa Cruz genome browser and plotted around the high-scoring sites.
Agreement with gene-enhancer pairing by long-range chromatin interactions
From 5C data produced at the University of Massachusetts in threecell lines (Gm12878, HeLA-S3, K562) we sorted out the enhancer sitesthat were probed by the 5C, for which we also have at least one CpG site available to our analysis at±500 base pairs from the probed enhancer. We then located the TSS that obtained the highest score in relation to the enhancer CpG. This yielded 318 interactions, at distances of 10 kb to 1 Mb between the CpG side and the TSS side. We then counted the number of cases in which the maximum number of interaction sequence reads (average of two repeats over three cell types) washigher than the average interaction read number (98.23 reads). P -values were obtained by comparing the binomial probability of the agreement in high versus low scoringsites. We also compared the number of 5C interaction reads inhigh-scoring (score≥0.85) enhancer-TSS pairs, versusthe number interactions between these enhancers and other TSSs at 10 kbto 1 Mb from the enhancer sites. The P -value of the difference was calculated by the Wilcoxon Rank-Sum test.
Replication in uncultured samples
Methylation and expression data for bladder, lymph node, ureter, lung, stomach, skeletal muscles and adipose tissue biopsies of twoindividuals were downloaded from a recently-published dataset [GSE:30654]. Out of the enhancer-promoter interactions that were predicted in the original set of the 58 cell types, 876 were also assessable in the new dataset (that is, methylation data were available for the relevant CpG sites, and expression data were available for the interacting genes). We re-scored these gene-CpG pairs, using our model, in the new dataset. If the interactions that were predicted in the 58 cell types were not relevant to interactions in the uncultured tissue samples, than we expected that 20.5 ±4.5out of the 876 would (by chance) obtain high-scoring (score>0.9) interactionsin individual #1, and 19.9 ±3.5in individual #2. Instead, we observed 29 in individual #1 (P = 0.0145) and 28 (P = 0.0084) in individual #2. The P-value (normal distribution, two-sided) indicated in the result and discussion section is for the average of both individuals (expected = 20.21±3.206, observed = 28.5 interactions, P = 0.0098).
Alteration in cancer
Differences betweennormal and cancer cells were measured by Robust Multi-chip Average expression units and by methylation percentages. The distributions of methylation and expression differences between normal and cancer samples (Figure 4) were smoothed (color density representation) using a kernel density estimate (transformation function = x^3). Linear regressions were performed (red lines and Pearson coefficients in Figure 4). The P-value of the frequencies of signals in the upper-left and lower-right quarters, in enhancers versus promoters, were calculated using normal distribution of two-proportion z test. P -valuesof the overlap between genes in the cancer types were calculated based on the expectation from random intersections between the groups.
Gene ontology analysis
Analysis of GO terms of thehypomethylatedupregulated genes was done with the GOrilla tool , using the list of high-scoring genes as a background.
chromosome conformation capture carbon copy technique
mean average precision
reduced representation bisulfite sequencing
Support Vector Machine
transcription start site
variable methylation site.
We are grateful to Shai Shalev-Shwartz for his helpful advice regarding the development of the machine-learning algorithm. Thanks are also due to Yehudit Bergman, Gadi Glaser, Amir Eden, Benjamin Glaser, Reuven Laskov and Howard Cedar for helpful discussions and comments. This work was supported by the Adams Fellowship Program of the Israel Academy of Sciences and Humanities (to SS), by the I-CORE Program of the Planning and Budgeting Committee and the Israel Science Foundation (Grant No. 41/11 to AH), and by the Israel Science Foundation (Grant No. 887/08to AH).
- Klose RJ, Bird AP: Genomic DNA methylation: the mark and its mediators. Trends Biochem Sci. 2006, 31 (2): 89-97. 10.1016/j.tibs.2005.12.008.PubMedView ArticleGoogle Scholar
- Hellman A, Chess A: Gene body-specific methylation on the active × chromosome. Science. 2007, 315 (5815): 1141-1143. 10.1126/science.1136352.PubMedView ArticleGoogle Scholar
- Suzuki MM, Bird A: DNA methylation landscapes: provocative insights from epigenomics. Nat Rev Genet. 2008, 9 (6): 465-476. 10.1038/nrg2341.PubMedView ArticleGoogle Scholar
- Zemach A, McDaniel IE, Silva P, Zilberman D: Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science. 2010, 328 (5980): 916-919. 10.1126/science.1186366.PubMedView ArticleGoogle Scholar
- Feng S, Cokus SJ, Zhang X, Chen PY, Bostick M, Goll MG, Hetzel J, Jain J, Strauss SH, Halpern ME, Ukomadu C, Sadler KC, Pradhan S, Pellegrini M, Jacobsen SE: Conservation and divergence of methylation patterning in plants and animals. Proc Natl Acad Sci USA. 2010, 107 (19): 8689-8694. 10.1073/pnas.1002720107.PubMedPubMed CentralView ArticleGoogle Scholar
- Aran D, Toperoff G, Rosenberg M, Hellman A: Replication timing-related and gene body-specific methylation of active human genes. Hum Mol Genet. 2011, 20 (4): 670-680. 10.1093/hmg/ddq513.PubMedView ArticleGoogle Scholar
- Cedar H, Bergman Y: Programming of DNA methylation patterns. Annu Rev Biochem. 2012, 81: 97-117. 10.1146/annurev-biochem-052610-091920.PubMedView ArticleGoogle Scholar
- Gal-Yam EN, Egger G, Iniguez L, Holster H, Einarsson S, Zhang X, Lin JC, Liang G, Jones PA, Tanay A: Frequent switching of Polycomb repressive marks and DNA hypermethylation in the PC3 prostate cancer cell line. Proc Natl Acad Sci USA. 2008, 105 (35): 12979-12984. 10.1073/pnas.0806437105.PubMedPubMed CentralView ArticleGoogle Scholar
- Schlesinger Y, Straussman R, Keshet I, Farkash S, Hecht M, Zimmerman J, Eden E, Yakhini Z, Ben-Shushan E, Reubinoff BE, Bergman Y, Simon I, Cedar H: Polycomb-mediated methylation on Lys27 of histone H3 pre-marks genes for de novo methylation in cancer. Nat Genet. 2007, 39 (2): 232-236. 10.1038/ng1950.PubMedView ArticleGoogle Scholar
- Ohm JE, McGarvey KM, Yu X, Cheng L, Schuebel KE, Cope L, Mohammad HP, Chen W, Daniel VC, Yu W, Berman DM, Jenuwein T, Pruitt K, Sharkis SJ, Watkins DN, Herman JG, Baylin SB: A stem cell-like chromatin pattern may predispose tumor suppressor genes to DNA hypermethylation and heritable silencing. Nat Genet. 2007, 39 (2): 237-242. 10.1038/ng1972.PubMedPubMed CentralView ArticleGoogle Scholar
- Widschwendter M, Fiegl H, Egle D, Mueller-Holzner E, Spizzo G, Marth C, Weisenberger DJ, Campan M, Young J, Jacobs I, Laird PW: Epigenetic stem cell signature in cancer. Nat Genet. 2007, 39 (2): 157-158. 10.1038/ng1941.PubMedView ArticleGoogle Scholar
- Easwaran H, Johnstone SE, Van Neste L, Ohm J, Mosbruger T, Wang Q, Aryee MJ, Joyce P, Ahuja N, Weisenberger D, Collisson E, Zhu J, Yegnasubramanian S, Matsui W, Baylin SB: A DNA hypermethylation module for the stem/progenitor cell signature of cancer. Genome Res. 2012, 22 (5): 837-849. 10.1101/gr.131169.111.PubMedPubMed CentralView ArticleGoogle Scholar
- Sproul D, Nestor C, Culley J, Dickson JH, Dixon JM, Harrison DJ, Meehan RR, Sims AH, Ramsahoye BH: Transcriptionally repressed genes become aberrantly methylated and distinguish tumors of different lineages in breast cancer. Proc Natl Acad Sci USA. 2011, 108 (11): 4364-4369. 10.1073/pnas.1013224108.PubMedPubMed CentralView ArticleGoogle Scholar
- De Carvalho DD, Sharma S, You JS, Su SF, Taberlay PC, Kelly TK, Yang X, Liang G, Jones PA: DNA methylation screening identifies driver epigenetic events of cancer cell survival. Cancer Cell. 2012, 21 (5): 655-667. 10.1016/j.ccr.2012.03.045.PubMedPubMed CentralView ArticleGoogle Scholar
- Lengauer C, Kinzler KW, Vogelstein B: DNA methylation and genetic instability in colorectal cancer cells. Proc Natl Acad Sci USA. 1997, 94 (6): 2545-2550. 10.1073/pnas.94.6.2545.PubMedPubMed CentralView ArticleGoogle Scholar
- Eden A, Gaudet F, Waghmare A, Jaenisch R: Chromosomal instability and tumors promoted by DNA hypomethylation. Science. 2003, 300 (5618): 455-10.1126/science.1083557.PubMedView ArticleGoogle Scholar
- Hon GC, Hawkins RD, Caballero OL, Lo C, Lister R, Pelizzola M, Valsesia A, Ye Z, Kuan S, Edsall LE, Camargo AA, Stevenson BJ, Ecker JR, Bafna V, Strausberg RL, Simpson AJ, Ren B: Global DNA hypomethylation coupled to repressive chromatin domain formation and gene silencing in breast cancer. Genome Res. 2012, 22 (2): 246-258. 10.1101/gr.125872.111.PubMedPubMed CentralView ArticleGoogle Scholar
- Raddatz G, Gao Q, Bender S, Jaenisch R, Lyko F: Dnmt3a protects active chromosome domains against cancer-associated hypomethylation. PLoS Genet. 2012, 8 (12): e1003146-10.1371/journal.pgen.1003146.PubMedPubMed CentralView ArticleGoogle Scholar
- Bulger M, Groudine M: Functional and mechanistic diversity of distal transcription enhancers. Cell. 2011, 144 (3): 327-339. 10.1016/j.cell.2011.01.024.PubMedPubMed CentralView ArticleGoogle Scholar
- Stadler MB, Murr R, Burger L, Ivanek R, Lienert F, Scholer A, Wirbelauer C, Oakeley EJ, Gaidatzis D, Tiwari VK, Schubeler D: DNA-binding factors shape the mouse methylome at distal regulatory regions. Nature. 2011, 480 (7378): 490-495.PubMedGoogle Scholar
- Xu J, Watts JA, Pope SD, Gadue P, Kamps M, Plath K, Zaret KS, Smale ST: Transcriptional competence and the active marking of tissue-specific enhancers by defined transcription factors in embryonic and induced pluripotent stem cells. Genes Dev. 2009, 23 (24): 2824-2838. 10.1101/gad.1861209.PubMedPubMed CentralView ArticleGoogle Scholar
- Sandovici I, Smith NH, Nitert MD, Ackers-Johnson M, Uribe-Lewis S, Ito Y, Jones RH, Marquez VE, Cairns W, Tadayyon M, O'Neill LP, Murrell A, Ling C, Constancia M, Ozanne SE: Maternal diet and aging alter the epigenetic control of a promoter-enhancer interaction at the Hnf4a gene in rat pancreatic islets. Proc Natl Acad Sci USA. 2011, 108 (13): 5449-5454. 10.1073/pnas.1019007108.PubMedPubMed CentralView ArticleGoogle Scholar
- Carone BR, Fauquier L, Habib N, Shea JM, Hart CE, Li R, Bock C, Li C, Gu H, Zamore PD, Meissner A, Weng Z, Hofmann HA, Friedman N, Rando OJ: Paternally induced transgenerational environmental reprogramming of metabolic gene expression in mammals. Cell. 2010, 143 (7): 1084-1096. 10.1016/j.cell.2010.12.008.PubMedPubMed CentralView ArticleGoogle Scholar
- Taberlay PC, Kelly TK, Liu CC, You JS, De Carvalho DD, Miranda TB, Zhou XJ, Liang G, Jones PA: Polycomb-repressed genes have permissive enhancers that initiate reprogramming. Cell. 2011, 147 (6): 1283-1294. 10.1016/j.cell.2011.10.040.PubMedPubMed CentralView ArticleGoogle Scholar
- Wiench M, John S, Baek S, Johnson TA, Sung MH, Escobar T, Simmons CA, Pearce KH, Biddie SC, Sabo PJ, Thurman RE, Stamatoyannopoulos JA, Hager GL: DNA methylation status predicts cell type-specific enhancer activity. EMBO J. 2011, 30 (15): 3028-3039. 10.1038/emboj.2011.210.PubMedPubMed CentralView ArticleGoogle Scholar
- Bock C, Beerman I, Lien WH, Smith ZD, Gu H, Boyle P, Gnirke A, Fuchs E, Rossi DJ, Meissner A: DNA methylation dynamics during in vivo differentiation of blood and skin stem cells. Mol Cell. 2012, 47: 633-647. 10.1016/j.molcel.2012.06.019.PubMedPubMed CentralView ArticleGoogle Scholar
- Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, Zhang X, Wang L, Issner R, Coyne M, Ku M, Durham T, Kellis M, Bernstein BE: Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011, 473 (7345): 43-49. 10.1038/nature09906.PubMedPubMed CentralView ArticleGoogle Scholar
- Li G, Ruan X, Auerbach RK, Sandhu KS, Zheng M, Wang P, Poh HM, Goh Y, Lim J, Zhang J, Sim HS, Peh SQ, Mulawadi FH, Ong CT, Orlov YL, Hong S, Zhang Z, Landt S, Raha D, Euskirchen G, Wei CL, Ge W, Wang H, Davis C, Fisher-Aylor KI, Mortazavi A, Gerstein M, Gingeras T, Wold B, Sun Y, et al: Extensive promoter-centered chromatin interactions provide a topological basis for transcription regulation. Cell. 2012, 148 (12): 84-98.PubMedPubMed CentralView ArticleGoogle Scholar
- Sanyal A, Lajoie BR, Jain G, Dekker J: The long-range interaction landscape of gene promoters. Nature. 2012, 489 (7414): 109-113. 10.1038/nature11279.PubMedPubMed CentralView ArticleGoogle Scholar
- Nazor KL, Altun G, Lynch C, Tran H, Harness JV, Slavin I, Garitaonandia I, Muller FJ, Wang YC, Boscolo FS, Fakunle E, Dumevska B, Lee S, Park HS, Olee T, D'Lima DD, Semechkin R, Parast MM, Galat V, Laslett AL, Schmidt U, Keirstead HS, Loring JF, Laurent LC: Recurrent variations in DNA methylation in human pluripotent stem cells and their differentiated derivatives. Cell Stem Cell. 2012, 10 (5): 620-634. 10.1016/j.stem.2012.02.013.PubMedPubMed CentralView ArticleGoogle Scholar
- Visel A, Rubin EM, Pennacchio LA: Genomic views of distant-acting enhancers. Nature. 2009, 461 (7261): 199-205. 10.1038/nature08451.PubMedPubMed CentralView ArticleGoogle Scholar
- Akhtar-Zaidi B, Cowper-Sal-lari R, Corradin O, Saiakhova A, Bartels CF, Balasubramanian D, Myeroff L, Lutterbaugh J, Jarrar A, Kalady MF, Willis J, Moore JH, Tesar PJ, Laframboise T, Markowitz S, Lupien M, Scacheri PC: Epigenomic enhancer profiling defines a signature of colon cancer. Science. 2012, 336 (6082): 736-739. 10.1126/science.1217277.PubMedPubMed CentralView ArticleGoogle Scholar
- Oehler VG, Guthrie KA, Cummings CL, Sabo K, Wood BL, Gooley T, Yang T, Epping MT, Shou Y, Pogosova-Agadjanyan E, Ladne P, Stirewalt DL, Abkowitz JL, Radich JP: The preferentially expressed antigen in melanoma (PRAME) inhibits myeloid differentiation in normal hematopoietic and leukemic progenitor cells. Blood. 2009, 114 (15): 3299-3308. 10.1182/blood-2008-07-170282.PubMedPubMed CentralView ArticleGoogle Scholar
- Rodelsperger C, Guo G, Kolanczyk M, Pletschacher A, Kohler S, Bauer S, Schulz MH, Robinson PN: Integrative analysis of genomic, functional and protein interaction data predicts long-range enhancer-target gene interactions. Nucleic Acids Res. 2011, 39 (7): 2492-2502. 10.1093/nar/gkq1081.PubMedPubMed CentralView ArticleGoogle Scholar
- Tekpli X, Zienolddiny S, Skaug V, Stangeland L, Haugen A, Mollerup S: DNA methylation of the CYP1A1 enhancer is associated with smoking-induced genetic alterations in human lung. Int J Cancer. 2012, 131 (7): 1509-1516. 10.1002/ijc.27421.PubMedView ArticleGoogle Scholar
- Steinmann K, Richter AM, Dammann RH: Epigenetic silencing of erythropoietin in human cancers. Genes Cancer. 2011, 2 (1): 65-73. 10.1177/1947601911405043.PubMedPubMed CentralView ArticleGoogle Scholar
- Angrisano T, Sacchetti S, Natale F, Cerrato A, Pero R, Keller S, Peluso S, Perillo B, Avvedimento VE, Fusco A, Bruni CB, Lembo F, Santoro M, Chiariotti L: Chromatin and DNA methylation dynamics during retinoic acid-induced RET gene transcriptional activation in neuroblastoma cells. Nucleic Acids Res. 2011, 39 (6): 1993-2006. 10.1093/nar/gkq864.PubMedPubMed CentralView ArticleGoogle Scholar
- A user's guide to the encyclopedia of DNA elements (ENCODE). PLoS Biol. 2011, 9 (4): e1001046-10.1371/journal.pbio.1001046.
- Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon G, Tonti-Filippini J, Nery JR, Lee L, Ye Z, Ngo QM, Edsall L, Antosiewicz-Bourget J, Stewart R, Ruotti V, Millar AH, Thomson JA, Ren B, Ecker JR: Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 2009, 462 (7271): 315-322. 10.1038/nature08514.PubMedPubMed CentralView ArticleGoogle Scholar
- Yue Y, Finley T, Radlinski F, Joachims T: A support vector method for optimizing average precision. Proceedings of the 30th annual international ACM SIGIR conference on Research and development in information retrieval: July 23-27 2007; Amsterdam, The Netherlands. 2007, New York: ACM, 271-278.Google Scholar
- Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z: GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics. 2009, 10: 48-10.1186/1471-2105-10-48.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.