- Open Access
Tissue of origin determines cancer-associated CpG island promoter hypermethylation patterns
Genome Biologyvolume 13, Article number: R84 (2012)
Aberrant CpG island promoter DNA hypermethylation is frequently observed in cancer and is believed to contribute to tumor progression by silencing the expression of tumor suppressor genes. Previously, we observed that promoter hypermethylation in breast cancer reflects cell lineage rather than tumor progression and occurs at genes that are already repressed in a lineage-specific manner. To investigate the generality of our observation we analyzed the methylation profiles of 1,154 cancers from 7 different tissue types.
We find that 1,009 genes are prone to hypermethylation in these 7 types of cancer. Nearly half of these genes varied in their susceptibility to hypermethylation between different cancer types. We show that the expression status of hypermethylation prone genes in the originator tissue determines their propensity to become hypermethylated in cancer; specifically, genes that are normally repressed in a tissue are prone to hypermethylation in cancers derived from that tissue. We also show that the promoter regions of hypermethylation-prone genes are depleted of repetitive elements and that DNA sequence around the same promoters is evolutionarily conserved. We propose that these two characteristics reflect tissue-specific gene promoter architecture regulating the expression of these hypermethylation prone genes in normal tissues.
As aberrantly hypermethylated genes are already repressed in pre-cancerous tissue, we suggest that their hypermethylation does not directly contribute to cancer development via silencing. Instead aberrant hypermethylation reflects developmental history and the perturbation of epigenetic mechanisms maintaining these repressed promoters in a hypomethylated state in normal cells.
Aberrant DNA hypermethylation of CpG island (CGI) promoters (promoter hypermethylation) occurs in many cancers. This epigenetic reprogramming is associated with the absence of transcription and can occur at a number of known tumor suppressor genes, suggesting that it contributes to tumor progression by silencing the expression of affected genes . Although this model has been hugely influential, the significance of hypermethylation at CGIs in cancer has long been debated and questioned [2–4]. Also, despite intense study, the mechanisms directing promoter hypermethylation in cancer remain elusive and it is unclear whether the same mechanism operates in different cancer types. In colorectal cancer, a CGI hypermethylator phenotype (termed CIMP) has been described where hundreds of CGIs become coordinately hypermethylated during tumor progression [5, 6]. Similar methylator phenotypes have been reported to occur in cancers originating from other tissues [7–9]. In these cases, it is particularly unclear whether hypermethylation is the primary event responsible for the silencing of target genes, however based on the propensity of large numbers of genes to become re-activated by exposure to DNA de-methylating drugs, it has been suggested that this might be the case .
Hypermethylation also plays a role in the regulation of some genes during normal development, particularly at imprinted loci and at CGI promoters on the inactive X-chromosome (Xi) in female mammalian cells [11, 12]. During X-inactivation CGI hypermethylation occurs after gene silencing has already taken place [13, 14] and the initial silencing event does not require DNA methyltransferases [15, 16]. Absence of the maintenance methyltransferase, Dnmt1, in mice can lead to reactivation of the Xi later in development suggesting that in this case CGI hypermethylation acts as a stabilizing factor that maintains silencing . Where the temporal dynamics of gene inactivation have been studied for autosomal genes, hypermethylation occurs subsequent to repression by other mechanisms .
We have recently shown that genes whose promoters are hypermethylated in breast cancer cell lines and tumors are already repressed in the putative lineage of origin and that when methylation is removed in cancer cell lines, either pharmacologically or genetically, most hypermethylated genes do not become re-activated . This implies that the majority of cancer-associated CGI hypermethylation does not contribute to tumor progression under the classic model because it occurs at genes that are already switched off. Others have shown that hypermethylation of APC frequently occurs in gastric cancer, but at a promoter that is not utilized in normal gastric tissue  and that RUNX3, whose tumor suppressor gene status is largely based on the fact that it is frequently methylated in gastric cancer, is never expressed in the gastrointestinal epithelial cells that give rise to these tumors .
Here, we explore the generality of our observations in breast cancer by analyzing data derived from 1,154 tumors arising in 7 different human tissues. We show that variability in promoter CGI hypermethylation patterns between tumors is explained by variability in gene expression patterns between normal tissues and it is genes that are repressed in the pre-cancerous tissue that become preferentially hypermethylated in tumors. Our study represents the first comprehensive analysis of promoter CGI hypermethylation in different human cancers and we propose that the hypermethylation of repressed CGI promoters is a common feature of most cancers.
Tissue of origin determines promoter hypermethylation patterns in cancers
We have previously shown that cell lineage determines promoter hypermethylation patterns in breast cancer . To examine the generality of these observations in cancers arising in other tissues, we collected methylation profiling data from 1,149 tumors of 7 different cancer types: breast (Gene Expression Omnibus, , GEO:GSE31979), colorectal (GEO:GSE25062), prostate (GEO:GSE26126), lung (The Cancer Genome Atlas, TCGA) and ovarian tumors (TCGA), along with acute-myeloid leukemias (AMLs, TCGA) and glioblastomas (TCGA) [5, 8, 23–25]. These datasets were all generated using Illumina Infinium HumanMethylation27 BeadChip methylation arrays, facilitating their cross comparison. We used these data to define sets of genes that were frequently aberrantly hypermethylated in each of the seven cancer types (See Additional file 1, unmethylated in the corresponding normal tissue and methylated in >20% of cancer samples, see methods for details). Our analyses were limited to genes possessing CGI promoters because the hypermethylation of non-CGI promoters is not always associated with transcriptional repression [26, 27]. The number of frequently hypermethylated genes varied between cancer types with the greatest number found in colorectal and lung tumors (382 and 396 genes, respectively) and the least found in ovarian tumors (100 genes) (See Additional file 2, Figure S1A). To assess the reproducibility of these lists, we derived a second set of genes frequently aberrantly hypermethylated in breast tumors from a meta-analysis of three studies [7, 18, 28]. Of these 316 genes, 81.5% (256) were found in our original list, a highly significant overlap (P < 2 × 10-16, Fisher's exact test), demonstrating the reproducibility of our methodology. In total, 1,009 genes were prone to hypermethylation by this analysis in at least one type of cancer, including a number reported to be frequently hypermethylated in cancer (for example, APC, DAPK1, ESR1, GSTP1, SFRP genes and HOX genes) [29–31]. None of the 1,009 gene sets were common to all cancer types and roughly half (503 genes) were unique to a single cancer type.
The overall levels of DNA methylation at these 1,009 hypermethylation-prone genes varied dramatically within cancer types but were highest in colorectal tumors and lowest in ovarian tumors (Figure 1a andAdditional file 2 Figure S1B). Examination of the methylation profiles of the 1,009 genes in the different cancer samples revealed that 220 of the genes were consistently methylated in cancers of different tissues (in at least 5% of samples for each tissue, Figure 1a). However, 446 of the genes had variable methylation profiles and were hypermethylated in some cancer types but not in others (Figure 1a, tick marks). For example, 86 of the 1,009 hypermethylation prone genes were never methylated in breast tumors but were methylated in at least one other cancer type. To systematically analyze sources of variation in the methylation profiles of the 1,149 samples, we performed principal component analysis (PCA) on the methylation data for the set of 1,009 hypermethylation prone genes . The first principal component accounted for around 66% of the variance in the data and was significantly correlated with the median methylation level of the 1,009 hypermethylation prone genes (Figure 1b, R = 0.90, P < 2 × 10-16). The next three components of the data accounted for 10.4% of the variance in the data and clearly separated out the samples into the seven different tissue types (Figure 1c). These analyses indicate that a substantial number of genes are prone to hypermethylation in multiple cancer types but that the susceptibility of many other genes to hypermethylation in cancer is determined by tissue-type specific factors.
Genes prone to hypermethylation in cancer are not constitutively expressed
Having defined genes that were prone to hypermethylation in cancer, we next examined which factors affected their propensity to become hypermethylated. As a control, we derived a second set of genes that were resistant to hypermethylation in cancer (those that were never methylated in any of the 1,149 cancer samples tested; 2,123 genes). The hypermethylation-prone and -resistant gene sets were associated with different Gene Ontology (GO) terms (Figure 2a). In particular, resistant genes were enriched in housekeeping terms such as 'Mitotic Cell Cycle', and 'RNA Processing and Macromolecule Catabolic Process' whereas prone genes were enriched in developmental terms such as 'System Development' and 'Organ Development'. We have previously shown that genes hypermethylated in breast cancer cell lines are expressed in a tissue-specific fashion in normal tissues and these functional terms might suggest that genes hypermethylated in diverse primary cancers also have tissue-specific expression patterns .
We used a method based on information theory to directly quantify the degree of tissue-specificity in a given gene's expression pattern across nine normal tissues that were profiled by high-throughput mRNA sequencing (RNA-seq, Sequence Read Archive, SRA:SRA008403) [33–35], with a higher score equating to a more tissue-specific pattern of expression. Hypermethylation prone genes were significantly more tissue-specific than hypermethylation resistant genes (Figure 2b). We observed similar results when we defined the specificity of expression from a panel of 36 tissues profiled on microarrays (See Additional file 2, Figure S2A, GEO:GSE2361)  or varied the thresholds used to define hypermethylation prone genes (See Additional file 2, Figure S2B). Furthermore, genes frequently hypermethylated in each of the seven different cancers were also found to have tissue-specific expression patterns (See Additional file 2, Figure S2C) as were genes found to be hypermethylated in colorectal tumors by alternative methylation profiling techniques (methyl-binding domain pull-down and sequencing, MBD-seq, or whole genome bisulfite sequencing [37–39], Additional file 2, Figure S2D, SRA:SRA029584 and [40, 41]). Therefore, genes prone to hypermethylation in cancer are robustly associated with tissue-specific expression patterns in normal tissues. One possibility is that hypermethylation selectively accumulates at tissue specific genes because the disruption of many housekeeping genes might be cell-lethal. However, we found that a set of CGI promoter genes reported as recurrently mutated in breast tumors showed no preference towards either tissue specific or housekeeping expression patterns in normal tissues implying that the disruption of housekeeping genes is not necessarily lethal, at least to breast tumor cells (See Additional file 2, Figure S2E). Our analyses show that genes that are prone to hypermethylation in cancer are distinguished from those resistant to hypermethylation by their regulated expression pattern in normal tissues.
Aberrantly hypermethylated genes have conserved promoter regions
Based on genes hypermethylated in multiple cancer cell lines, one study has suggested that the transcriptional start sites (TSSs) of genes prone to hypermethylation are depleted of repetitive elements . We investigated whether this was also true of our set of hypermethylation prone genes derived from primary cancers. In our analyses, all three major classes of repetitive elements (LINEs, SINEs and long terminal repeats (LTRs)) were depleted from the TSSs of CGI promoters and to a lesser extent non-CGI promoters (See Additional file 2, Figure S3A). However, genes prone to hypermethylation in cancer had a significantly greater depletion of repetitive elements than hypermethylation resistant genes (Figure 3a). The greater depletion from the promoters of hypermethylation prone genes could be caused by an unknown activity of repetitive elements in protecting CGIs from aberrant hypermethylation as has been previously suggested . However, this model is inconsistent with both the hypermethylation of repetitive elements in normal tissues and their hypomethylation in cancer . Based on our observation that hypermethylation prone genes have tissue-specific expression patterns (Figure 2b), we considered an alternative scenario. The expression pattern of tissue-specific genes is often regulated by elements that lie distant to their promoter . The insertion of a transposable element close to a tissue-specific gene might be detrimental to its regulation because it could directly disrupt one of these regulatory elements or interrupt their interaction with the gene promoter. The depletion of repetitive elements seen at hypermethylation prone genes could, therefore, reflect an evolutionary need to preserve the proper developmental regulation of these genes.
Many of the bioinformatic techniques used to discover functional elements in the human genome use comparisons of the genomes of multiple species to infer their presence through evolutionary conservation . Therefore, a testable consequence of our hypothesis regarding the presence of regulatory elements in the vicinity of hypermethylation prone promoters is that we should detect a greater degree of evolutionary constraint or conservation around these promoters. We quantified the level of evolutionary conservation around transcription start sites using two different measures: one based on the rate of nucleotide substitutions between species  and the other based on the measurement of the rate of insertions and deletions between species . The profiles of these scores mirrored that of repetitive elements and the greatest conservation was seen directly over the TSS (Figure 3b). Conservation was greater downstream of the TSS relative to the upstream region, probably due to the presence of exonic sequences. However, hypermethylation-prone genes had significantly higher levels of conservation as measured by both scores, at the TSS and extending into the upstream and downstream regions (Figure 3b). Similar results were observed for hypermethylation prone genes defined from either MBD-seq or whole-genome bisulfite sequencing profiling of colorectal tumors [37–39] (See Additional file 2, Figure S3B) suggesting that this property was not an artifact of data generated from Illumina Infinium arrays.
Recently, long intergenic non-coding RNAs (lincRNAs) have been proposed to play a cis regulatory role at some tissue specific genes . Their presence is therefore a further surrogate of regulatory complexity at nearby genes, so we asked whether lincRNAs were enriched at hypermethylation prone genes. As predicted, we found that hypermethylation prone genes were significantly enriched in neighboring lincRNAs defined in a recent comprehensive analysis of human tissues when compared to hypermethylation resistant genes (Figure 3d) . Thus, hypermethylation prone genes are normally expressed in a tissue-specific manner and the vicinity of their promoters is depleted of repeats and is evolutionarily conserved compared to hypermethylation resistant genes. We propose that these characteristics result from an evolutionary need to preserve regulatory elements required for the proper regulation of genes prone to hypermethylation in cancer during normal development.
Variation in hypermethylation patterns in tumors is determined by gene expression patterns in the tissue of origin
Although repeat occupancy and conservation differ between hypermethylation-prone and -resistant genes, these factors displayed overlapping distributions for the two gene sets (See Additional file 2, Figure S3C and D). For example, some hypermethylation prone genes completely lacked SINE elements in the vicinity of their TSSs but other hypermethylation prone genes were found with more SINE elements than the average hypermethylation resistant gene (See Additional file 2, Figure S3C). Also, repeat occupancy and evolutionary conservation are invariant between different tissues and so do not explain the variable susceptibility of some genes to hypermethylation between cancers of different tissues (Figure 1a). Therefore, there must be other determinants of a gene's susceptibility to hypermethylation in a particular cancer.
To uncover such determinants, we considered genes with variable methylation between tumors (VM genes, 446 hypermethylation prone genes defined as being never hypermethylated in at least one cancer type, see Additional file 3 and Figure 1a, tick marks). For comparison, we also defined a set of 220 consistently methylated (CM) genes that are methylated in all 7 cancer types (≥5% of samples of each tumor type, see Additional file 4). Both VM and CM genes were expressed in a more tissue specific fashion, depleted in repetitive elements and evolutionarily conserved compared to hypermethylation resistant genes (Figure 4a andAdditional file 2, Figure S4A and B). However, the expression of CM genes in normal tissues was significantly more tissue-specific than VM genes (Figure 4a). This suggests an inverse relationship between a gene's breadth of expression in normal tissues and the number of cancers in which it becomes hypermethylated; that is, genes that are expressed in fewer tissues become hypermethylated in more tumor types. In support of this relationship, we observed a significant correlation between a gene's specificity of expression in normal tissues and the number of tumors in which it was frequently hypermethylated (See Additional file 2, Figure S4C).
We have previously demonstrated that a gene's expression status in normal cells is linked to its susceptibility to hypermethylation in breast cancer by showing that genes repressed in a lineage-specific fashion in the normal breast are prone to hypermethylation in different subtypes of breast cancer cell lines and tumors . We, therefore, examined whether gene expression patterns in normal tissues might explain the differential susceptibility to hypermethylation for VM genes in cancer. Examination of the list of VM genes along with their susceptibility suggested this might be the case. For example, PAX6 is prone to hypermethylation in cancer but not in glioblastomas (See Additional file 3). The gene is vital for the normal development of the brain and its expression persists into adulthood . Similarly, GFI1 is prone to hypermethylation in cancer but not in AML and is vital for normal hematopoiesis; mice and humans lacking functional GFI1 are neutropenic suggesting that GFI1 functions in myleopoiesis and is expressed in the cells from which AMLs originate [51, 52].
We tested if normal expression patterns determined hypermethylation susceptibility by considering VM genes with differential susceptibility in individual cancer types. VM genes that were frequently hypermethylated in breast tumors (67 genes) were repressed in the cells of origin of most breast tumors, luminal epithelial cells , as compared to normal breast stromal cells (Figure 4b, GEO:GSE16997). Conversely, VM genes that were never hypermethylated in breast tumors were active in luminal epithelial cells (Figure 4b, 86 genes). Similarly, VM genes resistant to hypermethylation in colorectal tumors, glioblastomas and lung tumors were significantly more active in the corresponding normal tissue than VM genes prone to hypermethylation in the same tumor type (Figure 4c, SRA:SRA008403), and genes that were hypermethylated in colorectal tumors, as defined by MDB-seq or whole-genome bisulfite sequencing, were also significantly less active than those that did not become hypermethylated (See Additional file 2, Figure S4D). Furthermore, expression status in normal tissues was predictive of aberrant hypermethylation in cancer as genes which were repressed in normal colon compared to normal liver were significantly more likely to be hypermethylated in colorectal tumors than genes that are active in normal colon but repressed in the normal liver (Figure 4d, one-sided Wilcoxon rank sum test P = 1.6 × 10-7, GEO:GSE13471). We have previously shown a similar preference for genes specifically repressed in luminal epithelial cells to be hypermethylated in breast tumors . Together these analyses suggest that variability in promoter hypermethylation patterns between cancer types results from the variability in gene expression patterns in normal tissues, and that genes that become hypermethylated in cancer are repressed in the pre-cancerous tissue of origin.
The aberrant hypermethylation of CGI promoters is of interest because it correlates with gene silencing and can occur at tumor suppressor genes . Here, we present the first comprehensive analysis of CGI promoter hypermethylation in multiple cancer types and show that the genes that are hypermethylated are already repressed in the normal tissues that give rise to these tumors (Figure 5a). The potential for a gene to act as a tumor suppressor is dependent on its activity in a particular cellular context. Our study, therefore, demonstrates that the major contribution of general CGI promoter hypermethylation to cancer cannot be the silencing of tumor suppressor genes because it affects genes that are already repressed in pre-cancerous tissue.
During normal development, DNA methylation accumulates at loci that are already repressed and may facilitate stable transcriptional repression rather than directly cause silencing . We have previously demonstrated that breast cancer cell lines and tumors of different lineages preferentially hypermethylate genes that are already silent in their equivalent normal cells . Our current study extends this to tumors arising in different tissues and suggests that the hypermethylation of repressed genes represents a universal principle across all cancers. Other studies of individual genes also support this conclusion, for example, RUNX3 is frequently hypermethylated in gastric cancers but is never expressed in normal gastric epithelia . Therefore, CGI promoter hypermethylation in cancer shares features with processes that occur in normal cells and does not necessarily represent a de novo aberrant mechanism.
Our analyses show that this model applies to the majority of hypermethylated genes found in tumors, but it has been proposed that within each tumor a few 'driver' genes are directly repressed by hypermethylation . Under this scenario, the hypermethylation of repressed genes could be a 'passenger' event and is a surrogate of epigenetic dysregulation. An analogous model is proposed for genetic mutations in cancer [56, 57]. Known tumor suppressor genes are hypermethylated in the tumors we analyzed but methylation of these genes generally occurs much more rarely than the hypermethylation of repressed genes, suggesting that a driver/passenger model may in fact apply (See Additional file 2, Table S1). For example, the hypermethylation of BRCA1 only occurs in 12% of ovarian cancers and 2% of breast cancers. We find that APC is hypermethylated more frequently (for example, in 33% of colorectal cancers) but it has multiple TSSs and a promoter that is repressed in normal gastric tissue has been shown to be the site of hypermethylation in gastric cancers . Therefore, the significance of frequent APC hypermethylation depends on whether it occurs at the major promoter in these tissues. It is unclear whether the hypermethylation of these potential driver genes occurs as a by-product of the process that results in the hypermethylation of repressed genes or by an alternative mechanism (for example, the direct selection of epimutations). It is known that MLH1 is frequently hypermethylated in colorectal tumors that possess a CIMP phenotype . However, if methylator phenotypes do generally contribute to the repression of driver genes, we would expect tumors with higher levels of promoter CGI hypermethylation to demonstrate more aggressive clinical behavior because they would be statistically more likely to have inactivated more tumor suppressor genes. Tumors with methylator phenotypes in colorectal cancer, breast cancer and glioblastoma correlate with better clinical prognosis [7, 8, 58].
Our results confirm a previous observation that the promoters of genes prone to aberrant hypermethylation in cancer are depleted of repetitive elements . However, we suggest that this occurs due to an evolutionary need to preserve the regulation of these genes in normal development rather than as a direct protective effect of repeats as was suggested . Our conclusion that hypermethylation resistant genes are primarily housekeeping genes is supported by another study that associated the presence of motifs for general transcription factors with resistance to hypermethylation in cancer . The features we associate with hypermethylation prone genes overlap with those of genes regulated by Polycomb repressive complexes. For example, the prototypical gene of this class, HOX genes, are found in clusters that are devoid of repetitive elements and are regulated by lincRNAs [60, 61]. Previous studies have linked promoter hypermethylation in cancer to Polycomb and the overlap between Polycomb-marked genes in embryonic stem (ES) cells and genes hypermethylated in cancer has been noted . However, the profile of Polycomb marks in a single cell type is constant and does not account for the variability in hypermethylated genes between cancers of different tissues. Furthermore, although sets of hypermethylated genes are statistically enriched in these Polycomb-marked genes, only 30% to 60% of hypermethylated genes carry these marks in ES cells [5, 18]. Polycomb-marked genes in ES cells carry bivalent histone modifications and are differentially activated or repressed in alternative cell lineages as differentiation proceeds [62, 63]. It is therefore possible that these genes are prone to methylation because they can be repressed in a tissue-specific fashion rather than because of their association with Polycomb in ES cells.
We find that the aberrant hypermethylation of repressed genes occurs in all cancer types analyzed implying that a common mechanism might be responsible for promoter hypermethylation in all cancers. Although the exact mechanism remains unknown, our results mean that proposed mechanisms must account for the specificity of hypermethylation for repressed genes (Figure 5b). For example, if aberrant hypermethylation results from the loss of an activity protecting CGIs from hypermethylation in normal cells  then the specificity of hypermethylation for repressed genes implies that different factors are responsible for maintaining hypomethylation at repressed and active CGIs or that hypomethylation is maintained at active CGIs via multiple redundant mechanisms that are not all present at the CGI promoters of repressed genes. TET (ten-eleven translocation) hydroxylase enzymes may be capable of mediating this protective activity through their proposed role in DNA demethylation  and inhibition of their enzymatic activity in cancer correlates with the hypermethylation of CGIs . Aberrant hypermethylation could also result from the recruitment of DNA methyltransferases (DNMTs) by transcription factors [66, 67]. However, transcription factors also activate genes and it remains to be demonstrated how these interactions might result in the specific hypermethylation of repressed genes. Over-expression of DNMT3B promotes tumorigenesis in a mouse model of colorectal cancer and is associated with the hypermethylation of specific genes . Higher DNMT3B levels have also been associated with the CIMP phenotype in human colorectal tumors [69, 70]. One of these studies also determined the stage in tumorigenesis at which different genes became hypermethylated showing that the repressed gene RUNX3 was the earliest CGI promoter to show significant change [20, 69], suggesting that differences in the expression of DNMTs could be linked to the hypermethylation of repressed genes.
Here we have shown that differences exist in the aberrant hypermethylation profiles of cancers arising in different tissue contexts. However, our results also make it clear that there is heterogeneity in the methylation profiles within particular types of cancer (Figure 1a). It is unclear how this heterogeneity arises but some mutations may play a direct role in its generation, for example, those that inhibit TET enzyme activity [8, 65]. Colorectal cancer has previously been split into at least three groups based on methylation profiles: non-CIMP tumors, CIMP-high tumors associated with BRAF mutations and CIMP-low mutations associated with KRAS mutations [5, 71]. Interestingly, a recent study suggested that CIMP-low tumors hypermethylate a subset of the genes hypermethylated in CIMP-high tumors rather than distinct sets of genes . Our own results might also suggest that variation between cancers in a given tissue can manifest itself as variable levels of methylation at methylation prone genes rather than the hypermethylation of alternative gene sets (see Figure 1a). We have previously shown that differences in the hypermethylation profiles of breast cancer subtypes of putatively different cells of origin can arise because of differences in gene expression in normal cell populations . Taken together, these results suggest that the hypermethylation of genes that are repressed in the normal cells of origin can account for the heterogeneity of tumor methylation profiles and variation in aberrant hypermethylation arises due to variations in the cells of origin or other factors, such as mutations, that influence the strength of the repressed gene methylator phenotype.
Recently, hydroxymethylated cytosine (hmC) has been re-discovered as a DNA modification present at significant levels in mammalian cells . The Illumina arrays that were used to generate most of the datasets we have analyzed are unable to distinguish methylated cytosine (mC) from hmC  and the results we present may relate to hmC rather than mC marked promoters in cancer. However, we have confirmed that these results equally apply in additional datasets derived by MBD pull-down, which is specific for 5mC (See Additional file 2, Figures S2D, S3B and S4D). In addition, hmC appears to be generally depleted in cancer [74–76]. It is likely, therefore, that repressed genes are prone to hypermethylation rather than hyperhydroxymethylation.
In summary, our results argue that the bulk of aberrant promoter hypermethylation in cancer occurs predominantly at genes that are repressed in pre-cancerous tissue and therefore does not directly contribute to tumor progression by silencing tumor suppressor genes. This epigenetic alteration is common to all the cancer types we have analyzed implying that a common mechanism is responsible for promoter hypermethylation at repressed genes in all cancers. Future research in this field should, therefore, focus on confirming whether aberrant hypermethylation does directly suppress rare driver genes and if the mechanism responsible for driver gene suppression is the same as that acting at repressed genes. Finally, we would suggest that researchers must exercise caution in assigning a tumor suppressor status to a gene based on its propensity to become hypermethylated in cancer.
Materials and methods
All statistical analyses were performed using the R statistical software (version 2.12.1) . Additional packages used are mentioned under the appropriate section.
Gene expression and methylation data used in this study were taken from previously published studies. The sources of the data are indicated in Table 1 and the number of samples in each dataset in Table 2.
In order to apply a consistent annotation to the data used in this study, all data were re-annotated to Ensembl 54 gene IDs (NCBI36). CpG probes from the Illumina Infinium arrays were mapped to the closest Ensembl gene based on TSS location using custom Perl and R scripts. CpGs that ambiguously mapped to more than one gene ID were removed from the analysis. CGI locations were taken from those biologically defined in a recent study . Similarly, expression data were mapped as previously described for Illumina expression arrays  or using publically available re-annotations for Affymetrix expression arrays . RNAseq data were mapped to Ensembl gene IDs as described below. lincRNAs were mapped to neighboring Ensembl gene IDs as described below.
Processing of methylation data
For data originating from Infinium methylation arrays, beta values were used as a measure of the methylation level at a given CpG probe (derived from the intensity of the methylated, I meth , and unmethylated, I unmeth , allele probes: I meth /(I meth + I unmeth ) ). We have previously shown that these are a reliable estimate of the level of methylation at a locus . These data were then filtered to remove unreliable values based on the detection P-value from the Infinium arrays (threshold 0.01). Methylation data originating from other techniques (MDB-seq or whole-genome bisulfite sequencing) were either downloaded as processed data provided by the authors [37, 38] or processed from raw sequencing files . We first downloaded raw sequencing data from the SRA . We then aligned these reads to the genome using Bowtie (version 0.12.7)  and the BEDtools (version 2.12.0) coverageBED tool to quantify the number of reads at each CGI . The read counts of CGIs were then normalized for CGI length and the total number of reads per sample to obtain a reads per kb per million mapped reads (RPKM) value for each CGI and the mean value taken from replicates of individual samples.
Definition of hypermethylation-prone and -resistant genes
Hypermethylation-prone and -resistant genes were defined from Illumina infinium array data using beta value cutoffs (roughly equating to percent methylation divided by 100). Previously, we have shown that probes with beta values <0.3 represent unmethylated areas of the genome  and we therefore defined unmethylated probes on this basis. In cell lines, we have previously shown that probes with beta >0.7 represent genomic loci that are fully methylated . However, in a preliminary analysis, we found that in the breast tumor samples used here, probes that had beta values >0.7 were also all methylated in normal breast tissue (data not shown). Probes that were aberrantly hypermethylated in these tumors had lower beta values because of the mix of cancerous and normal tissue in the samples analyzed. In this study, we therefore set a beta value threshold of >0.3 to define methylated probes. We only considered probes that were located within a CGI and within 200bp of a TSS that were unmethylated in all available normal samples from that tissue when defining gene sets (the 'all' genes control set for each tumor type). Frequently hypermethylated genes for each cancer were defined as genes satisfying these criteria that were methylated in at least 20% of tumor samples. Similarly, hypermethylation resistant genes satisfied these criteria but were not found to be methylated in any of the tumors. Genes present in both lists were then excluded from the analysis as being of ambiguous status to control for the presence of multiple probes at some genes.
To ensure that the method of gene selection did not bias our results, we also carried out analyses in which parameters were varied (See Additional file 2, Figure S2B and data not shown). We considered two major variations: we varied the threshold used to define aberrantly methylated genes and we varied the threshold required to call genes frequently aberrantly hypermethylated. In the first case, aberrantly hypermethylated genes were defined as those for which no probes had beta >0.3 in normal tissue and for which their mean beta value was >0.5 in at least 20% of cancers of a given type. In the second case, we varied the percent of samples required for a gene to be defined as frequently hypermethylated from 10% to 50%.
Two of the datasets used did not contain normal samples to define probes' normal tissue methylation status. In these cases, we made use of other datasets. For glioblastoma we used the fetal and adult brain samples from Sproul et al. . For AML we used the whole blood, neutrophils, B-cells, CD4 and CD8 T-cells, natural killer cells and CD34+ hematopoeitic stem cells samples from Calvanese et al. . We then defined the 1,009 hypermethylation prone genes from those being frequently hypermethylated in at least 1 of the 7 cancers, and the 2,123 hypermethylation resistant genes as those that were never methylated in any of the tumors analyzed. Consistently and variably prone genes (CM and VM, respectively) were defined as hypermethylation prone genes that were methylated in at least 5% of tumors of each type or never methylated in at least one tumor type, respectively.
To validate the reproducibility of our method of defining hypermethylation prone genes, we compared our list of genes frequently hypermethylated in breast tumors to a second list defined by the cross-comparison of three independent studies [7, 18, 28]. Data from these studies were either downloaded from GEO or from the author's website (GEO: GSE26349 and GSE26990) . Frequently hypermethylated genes were defined as above but only genes that were frequently hypermethylated in all three datasets were included in the analysis.
To define frequently hypermethylated and resistant genes from the Illingworth et al. MBD-seq data, we first generated lists of CGIs that were unmethylated in all of the normal colon samples . We then defined those CGIs that had higher levels of methylation in at least two of the tumor samples when compared to their matched normal samples as frequently hypermethylated CGIs. Resistant CGIs were defined as those that did not show higher levels of methylation in any of the tumors compared to their matched normal tissues. CGIs were assigned to genes if their transcription start site was present in the CGI. Genes present in both frequent and resistant lists were also removed because their status was ambiguous.
To define genes which were hypermethylated in colorectal tumors from the Xu et al. MDB-seq data , we used one-sided Wilcoxon rank sum tests to find CGIs with significantly more reads in tumor samples than normal samples (P < 0.05). CGIs were assigned to genes if their TSS was present in the CGI. Using this methodology, we were unable to define a set of hypermethylation resistant genes.
We defined genes prone to and resistant to hypermethylation from the Berman et al. whole- genome bisulfite sequencing data  as genes with CGI TSSs which were located in regions defined as methylation-prone or -resistant in that study. These regions were downloaded from the author's website. Genes which were defined as both methylation-prone and -resistant were excluded as being of ambiguous status.
Processing of expression data
To process RNA-seq data, raw sequence data for nine human tissues  were downloaded from GEO and converted to FASTQ format using the SRA Toolkit (version 2.1.7). Several technical replicates were available for each tissue. However, we randomly chose a single replicate in each case for simplicity and because different tissues had different numbers of replicates in this dataset. We examined each sample for per-base and per-read quality and over-represented kmers using the FastQC software (version 0.9.4) . Reads were mapped simultaneously to the human genome (NCBI version 36/hg18) and a library containing the sequences of all possible exon splice junctions (Ensembl 54 exons) created using RSEQtools . Reads were mapped using Bowtie (version 0.12.7)  allowing for a maximum of two mismatched bases and reporting the single best alignment for each read. RPKM values were computed for each ENSEMBL gene using reads mapping to exons and junctions of its longest transcript.
Processed Illumina gene expression data were downloaded from GEO and summarized to individual Ensembl IDs by taking the mean value of all probes mapped to that gene . Reported probe detection values were used to remove genes from the analysis for which all probes had a value >0.05 in all samples. Affymetrix expression data were summarized using the RMA algorithm from the Bioconducter affy package and an updated annotation . Detection calls were also generated using the MAS5 algorithm in the Bioconducter affy package.
To define CGI genes that were differentially expressed between normal colon and normal liver, we assayed for differential expression in processed expression microarray data using t-tests and Benjamini-Hochberg correction for multiple testing (assuming unequal variance and with a false discovery rate (FDR) of 5%). We then refined these lists to only those that had CGI promoters and had probes on the Infinium array within 200bp of their TSS that were unmethylated in normal colon (1,456 colon active and 356 colon repressed genes) before comparing methylation levels in colorectal tumors between these groups.
Analysis of GO-terms
To analyze functional terms, Ensembl Biomart was used to map gene identifiers to GO biological process terms (Ensembl 54). Enrichment of specific terms in each gene list was then assessed using Fisher's exact test as compared to all genes present on the Infinium array. Terms that were associated with less than 10 genes on the Infinium arrays were excluded from the analysis. Data were presented as change in the percent of genes in each set as compared to the control.
Definition of tissue specificity of gene expression
The specificity of a gene's expression pattern in normal tissues was measured using a method based on information theory . A low score indicates that a gene is uniformly expressed and a high score indicates that it is expressed specifically in one tissue. For plotting, we calculated the specificity for all genes and then split them into five equally sized groups of increasing specificity. We first removed any genes from the analysis that were potentially unexpressed in all assayed samples. For microarray expression data this was done using Affymetrix MAS5 detection calls by defining genes that were 'absent' in all samples as being unexpressed in all tissues. For RNA-seq, genes unexpressed in all tissues were defined as those that had RPKM values of 0 in all tissues. The specificity of individual gene sets was examined by plotting their distribution across the five specificity groups or by plotting the scores themselves. To test significance, the specificity scores of gene sets were compared using Wilcoxon rank sum tests.
Definition of genes mutated in breast cancer
Genes reported as mutated in breast cancer were defined using the Catalogue of Somatic Mutations in Cancer (COSMIC, ) database . Reports of genes mutated or not mutated in breast cancer samples were downloaded from COSMIC biomart (version 52B) and used to generate two lists of genes: those mutated in at least two samples and those not mutated in any. Lists were further limited to only those genes with CGI promoters (341 mutated and 10,117 non-mutated genes). The control set for this analysis was all CGI promoter genes reported as analyzed in breast cancer by COSMIC (11,022 genes).
Analysis of repetitive elements at promoters
To define the density of repetitive elements around TSSs, repetitive element positions were downloaded from the Repeat Masker track of the University of California, Santa Cruz (UCSC) genome browser (hg18) . Custom R scripts were then used to determine whether a repeat of a given class was present in a particular genomic interval. Repeats were defined as being present if they overlapped this interval. Genes were analyzed with respect to their TSS in non-overlapping 1 kb windows upstream and downstream of the TSS (with respect to the direction of transcription). To compare gene sets we plotted the frequency of repeats found at each window within that set. Differences between sets were tested using Fisher's exact tests based upon the density of repeats within a window ± 2 kb from TSSs.
Analysis of evolutionary conservation at promoters
We defined the level of conservation around gene promoters using two different measurements. The first was based upon the measurement of base substitutions between 17 vertebrate species . To define this score, the 'aggregate' tool from the Galaxy suite of bioinformatic tools was used to generate mean Phastcons (conservation) scores in 500bp windows surrounding each TSS (using the Phastcons 17-vertebrate alignments from genome build hg18). The second score was defined using data on sequences that showed a significant depletion of short insertions and deletions in comparisons of multiple species . We downloaded the locations of these sequences from the UCSC browser (hg18)  and then used the coverageBed tool from the BEDtools suite of bioinformatics tools  to calculate the percent of a given genomic interval occupied by these insertion and deletion purified sequences. Genes were analyzed with respect to their TSS in non-overlapping 500bp windows upstream and downstream of the TSS (with respect to the direction of transcription) by deriving mean scores for genes in the set. The significance of differences between gene sets was tested using the scores calculated for the window -2 to +2 kb from the TSS and Wilcoxon rank sum tests.
Analysis of genes neighboring lincRNAs
We used a recent survey of lincRNAs in the human genome to define genes that had a neighboring lincRNA . The nearest genic neighbor of each lincRNA was defined from that study's supplementary data and the given Refseq IDs mapped to Ensembl gene IDs. Gene sets were compared by examining the proportion of genes that had a lincRNA as their closest neighbor and significance was tested using Fisher's exact tests.
acute myeloid leukemia
Catalogue of Somatic Mutations in Cancer
- ES cell:
embryonic stem cell
Gene Expression Omnibus
long intergenic non-coding RNA
long interspersed element
long terminal repeat
principal component analysis
reads per kb per million mapped reads
short interspersed element
sequence read archive
The Cancer Genome Atlas
transcription start site
Jones PA, Baylin SB: The fundamental role of epigenetic events in cancer. Nat Rev Genet. 2002, 3: 415-428.
Bestor TH: Unanswered questions about the role of promoter methylation in carcinogenesis. Ann N Y Acad Sci. 2003, 983: 22-27. 10.1111/j.1749-6632.2003.tb05959.x.
Deaton AM, Bird A: CpG islands and the regulation of transcription. Genes Dev. 2011, 25: 1010-1022. 10.1101/gad.2037511.
Fearon ER: BRCA1 and E-cadherin promoter hypermethylation and gene inactivation in cancer-association or mechanism?. J Natl Cancer Inst. 2000, 92: 515-517. 10.1093/jnci/92.7.515.
Hinoue T, Weisenberger DJ, Lange CP, Shen H, Byun HM, Van Den Berg D, Malik S, Pan F, Noushmehr H, van Dijk CM, Tollenaar RA, Laird PW: Genome-scale analysis of aberrant DNA methylation in colorectal cancer. Genome Res. 2011, 22: 271-282.
Toyota M, Ahuja N, Ohe-Toyota M, Herman JG, Baylin SB, Issa JP: CpG island methylator phenotype in colorectal cancer. Proc Natl Acad Sci USA. 1999, 96: 8681-8686. 10.1073/pnas.96.15.8681.
Fang F, Turcan S, Rimner A, Kaufman A, Giri D, Morris LG, Shen R, Seshan V, Mo Q, Heguy A, Baylin SB, Ahuja N, Viale A, Massague J, Norton L, Vahdat LT, Moynahan ME, Chan TA: Breast cancer methylomes establish an epigenomic foundation for metastasis. Sci Transl Med. 2011, 3: 75ra25-10.1126/scitranslmed.3001875.
Noushmehr H, Weisenberger DJ, Diefes K, Phillips HS, Pujara K, Berman BP, Pan F, Pelloski CE, Sulman EP, Bhat KP, Verhaak RG, Hoadley KA, Hayes DN, Perou CM, Schmidt HK, Ding L, Wilson RK, Van Den Berg D, Shen H, Bengtsson H, Neuvial P, Cope LM, Buckley J, Herman JG, Baylin SB, Laird PW, Aldape K: Identification of a CpG island methylator phenotype that defines a distinct subgroup of glioma. Cancer Cell. 2010, 17: 510-522. 10.1016/j.ccr.2010.03.017.
Teodoridis JM, Hardie C, Brown R: CpG island methylator phenotype (CIMP) in cancer: causes and implications. Cancer Lett. 2008, 268: 177-186. 10.1016/j.canlet.2008.03.022.
Schuebel KE, Chen W, Cope L, Glockner SC, Suzuki H, Yi JM, Chan TA, Van Neste L, Van Criekinge W, van den Bosch S, van Engeland M, Ting AH, Jair K, Yu W, Toyota M, Imai K, Ahuja N, Herman JG, Baylin SB: Comparing the DNA hypermethylome with gene mutations in human colorectal cancer. PLoS Genet. 2007, 3: 1709-1723.
Ferguson-Smith AC: Genomic imprinting: the emergence of an epigenetic paradigm. Nat Rev Genet. 2011, 12: 565-575.
Wutz A: Gene silencing in X-chromosome inactivation: advances in understanding facultative heterochromatin formation. Nat Rev Genet. 2011, 12: 542-553. 10.1038/nrg3035.
Lock LF, Takagi N, Martin GR: Methylation of the Hprt gene on the inactive × occurs after chromosome inactivation. Cell. 1987, 48: 39-46. 10.1016/0092-8674(87)90353-9.
Okamoto I, Heard E: Lessons from comparative analysis of X-chromosome inactivation in mammals. Chromosome Res. 2009, 17: 659-669. 10.1007/s10577-009-9057-7.
Sado T, Fenner MH, Tan SS, Tam P, Shioda T, Li E: X inactivation in the mouse embryo deficient for Dnmt1: distinct effect of hypomethylation on imprinted and random × inactivation. Dev Biol. 2000, 225: 294-303. 10.1006/dbio.2000.9823.
Sado T, Okano M, Li E, Sasaki H: De novo DNA methylation is dispensable for the initiation and propagation of × chromosome inactivation. Development. 2004, 131: 975-982. 10.1242/dev.00995.
Feldman N, Gerson A, Fang J, Li E, Zhang Y, Shinkai Y, Cedar H, Bergman Y: G9a-mediated irreversible epigenetic inactivation of Oct-3/4 during early embryogenesis. Nat Cell Biol. 2006, 8: 188-194. 10.1038/ncb1353.
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: 4364-4369. 10.1073/pnas.1013224108.
Hosoya K, Yamashita S, Ando T, Nakajima T, Itoh F, Ushijima T: Adenomatous polyposis coli 1A is likely to be methylated as a passenger in human gastric carcinogenesis. Cancer Lett. 2009, 285: 182-189. 10.1016/j.canlet.2009.05.016.
Levanon D, Bernstein Y, Negreanu V, Bone KR, Pozner A, Eilam R, Lotem J, Brenner O, Groner Y: Absence of Runx3 expression in normal gastrointestinal epithelium calls into question its tumour suppressor function. EMBO Mol Med. 2011, 3: 593-604. 10.1002/emmm.201100168.
Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo/]
The Cancer Genome Atlas. [http://cancergenome.nih.gov/]
Fackler MJ, Umbricht CB, Williams D, Argani P, Cruz LA, Merino VF, Teo WW, Zhang Z, Huang P, Visvananthan K, Marks J, Ethier S, Gray JW, Wolff AC, Cope LM, Sukumar S: Genome-wide methylation analysis identifies genes specific to breast cancer hormone receptor status and risk of recurrence. Cancer Res. 2011, 71: 6195-6207. 10.1158/0008-5472.CAN-11-1630.
Kobayashi Y, Absher DM, Gulzar ZG, Young SR, McKenney JK, Peehl DM, Brooks JD, Myers RM, Sherlock G: DNA methylation profiling reveals novel biomarkers and important roles for DNA methyltransferases in prostate cancer. Genome Res. 2011, 21: 1017-1027. 10.1101/gr.119487.110.
Cancer Genome Atlas Research Network: Integrated genomic analyses of ovarian carcinoma. Nature. 2011, 474: 609-615. 10.1038/nature10166.
Boyes J, Bird A: Repression of genes by DNA methylation depends on CpG density and promoter strength: evidence for involvement of a methyl-CpG binding protein. EMBO J. 1992, 11: 327-333.
Weber M, Hellmann I, Stadler MB, Ramos L, Paabo S, Rebhan M, Schubeler D: Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome. Nat Genet. 2007, 39: 457-466. 10.1038/ng1990.
Van der Auwera I, Yu W, Suo L, Van Neste L, van Dam P, Van Marck EA, Pauwels P, Vermeulen PB, Dirix LY, Van Laere SJ: Array-based DNA methylation profiling for breast cancer subtype discrimination. PLoS One. 2010, 5: e12616-10.1371/journal.pone.0012616.
Herman JG, Umar A, Polyak K, Graff JR, Ahuja N, Issa JP, Markowitz S, Willson JK, Hamilton SR, Kinzler KW, Kane MF, Kolodner RD, Vogelstein B, Kunkel TA, Baylin SB: Incidence and functional consequences of hMLH1 promoter hypermethylation in colorectal carcinoma. Proc Natl Acad Sci USA. 1998, 95: 6870-6875. 10.1073/pnas.95.12.6870.
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: 237-242. 10.1038/ng1972.
Shah N, Sukumar S: The Hox genes and their roles in oncogenesis. Nat Rev Cancer. 2010, 10: 361-371. 10.1038/nrc2826.
Ringner M: What is principal component analysis?. Nat Biotechnol. 2008, 26: 303-304. 10.1038/nbt0308-303.
Martinez O, Reyes-Valdes MH: Defining diversity, specialization, and gene specificity in transcriptomes through information theory. Proc Natl Acad Sci USA. 2008, 105: 9709-9714. 10.1073/pnas.0803479105.
Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456: 470-476. 10.1038/nature07509.
Sequence Read Archive. [http://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi]
Ge X, Yamamoto S, Tsutsumi S, Midorikawa Y, Ihara S, Wang SM, Aburatani H: Interpreting expression profiles of cancers by genome-wide survey of breadth of expression in normal tissues. Genomics. 2005, 86: 127-141. 10.1016/j.ygeno.2005.04.008.
Berman BP, Weisenberger DJ, Aman JF, Hinoue T, Ramjan Z, Liu Y, Noushmehr H, Lange CP, van Dijk CM, Tollenaar RA, Van Den Berg D, Laird PW: Regions of focal DNA hypermethylation and long-range hypomethylation in colorectal cancer coincide with nuclear lamina-associated domains. Nat Genet. 2011, 44: 40-46. 10.1038/ng.969.
Illingworth RS, Gruenewald-Schneider U, Webb S, Kerr AR, James KD, Turner DJ, Smith C, Harrison DJ, Andrews R, Bird AP: Orphan CpG islands identify numerous conserved promoters in the mammalian genome. PLoS Genet. 2010, 6: e1001134-10.1371/journal.pgen.1001134.
Xu Y, Hu B, Choi AJ, Gopalan B, Lee BH, Kalady MF, Church JM, Ting AH: Unique DNA methylome profiles in CpG island methylator phenotype colon cancers. Genome Res. 2012, 22: 283-291. 10.1101/gr.122788.111.
Illingworth et al.: 2010, Publication Supplemenatary Dataset S1. [http://www.plosgenetics.org/article/fetchSingleRepresentation.action?uri=info:doi/10.1371/journal.pgen.1001134.s001]
USC Epigenome Center Supplemental Material for Berman et al. 2011, [http://epigenome.usc.edu/publicationdata/berman20101101/]
Estecio MR, Gallegos J, Vallot C, Castoro RJ, Chung W, Maegawa S, Oki Y, Kondo Y, Jelinek J, Shen L, Hartung H, Aplan PD, Czerniak BA, Liang S, Issa JP: Genome architecture marked by retrotransposons modulates predisposition to DNA methylation in cancer. Genome Res. 2010, 20: 1369-1382. 10.1101/gr.107318.110.
Wild L, Flanagan JM: Genome-wide hypomethylation in cancer may be a passive consequence of transformation. Biochim Biophys Acta. 2010, 1806: 50-57.
Kleinjan DA, van Heyningen V: Long-range control of gene expression: emerging mechanisms and disruption in disease. Am J Hum Genet. 2005, 76: 8-32. 10.1086/426833.
Ponting CP, Hardison RC: What fraction of the human genome is functional?. Genome Res. 2011, 21: 1769-1776. 10.1101/gr.116814.110.
Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, Clawson H, Spieth J, Hillier LW, Richards S, Weinstock GM, Wilson RK, Gibbs RA, Kent WJ, Miller W, Haussler D: Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res. 2005, 15: 1034-1050. 10.1101/gr.3715005.
Lunter G, Ponting CP, Hein J: Genome-wide identification of human functional DNA using a neutral indel model. PLoS Comput Biol. 2006, 2: e5-10.1371/journal.pcbi.0020005.
Wang KC, Chang HY: Molecular mechanisms of long noncoding RNAs. Mol Cell. 2011, 43: 904-914. 10.1016/j.molcel.2011.08.018.
Cabili MN, Trapnell C, Goff L, Koziol M, Tazon-Vega B, Regev A, Rinn JL: Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. Genes Dev. 2011, 25: 1915-1927. 10.1101/gad.17446611.
Hanson I, Van Heyningen V: Pax6: more than meets the eye. Trends Genet. 1995, 11: 268-272. 10.1016/S0168-9525(00)89073-3.
Karsunky H, Zeng H, Schmidt T, Zevnik B, Kluge R, Schmid KW, Duhrsen U, Moroy T: Inflammatory reactions and severe neutropenia in mice lacking the transcriptional repressor Gfi1. Nat Genet. 2002, 30: 295-300. 10.1038/ng831.
Person RE, Li FQ, Duan Z, Benson KF, Wechsler J, Papadaki HA, Eliopoulos G, Kaufman C, Bertolone SJ, Nakamoto B, Papayannopoulou T, Grimes HL, Horwitz M: Mutations in proto-oncogene GFI1 cause human neutropenia and target ELA2. Nat Genet. 2003, 34: 308-312. 10.1038/ng1170.
Gusterson B: Do 'basal-like' breast cancers really exist?. Nat Rev Cancer. 2009, 9: 128-134. 10.1038/nrc2571.
Baylin SB, Jones PA: A decade of exploring the cancer epigenome - biological and translational implications. Nat Rev Cancer. 2011, 11: 726-734. 10.1038/nrc3130.
Bird A: DNA methylation patterns and epigenetic memory. Genes Dev. 2002, 16: 6-21. 10.1101/gad.947102.
Haber DA, Settleman J: Cancer: drivers and passengers. Nature. 2007, 446: 145-146. 10.1038/446145a.
Sieber OM, Heinimann K, Tomlinson IP: Genomic instability--the engine of tumorigenesis?. Nat Rev Cancer. 2003, 3: 701-708. 10.1038/nrc1170.
Ogino S, Nosho K, Kirkner GJ, Kawasaki T, Meyerhardt JA, Loda M, Giovannucci EL, Fuchs CS: CpG island methylator phenotype, microsatellite instability, BRAF mutation and clinical outcome in colon cancer. Gut. 2009, 58: 90-96. 10.1136/gut.2008.155473.
Gebhard C, Benner C, Ehrich M, Schwarzfischer L, Schilling E, Klug M, Dietmaier W, Thiede C, Holler E, Andreesen R, Rehli M: General transcription factor binding at CpG islands in normal cells correlates with resistance to de novo DNA methylation in cancer cells. Cancer Res. 2010, 70: 1398-1407. 10.1158/0008-5472.CAN-09-3406.
Simons C, Pheasant M, Makunin IV, Mattick JS: Transposon-free regions in mammalian genomes. Genome Res. 2006, 16: 164-172.
Hekimoglu B, Ringrose L: Non-coding RNAs in polycomb/trithorax regulation. RNA Biol. 2009, 6: 129-137. 10.4161/rna.6.2.8178.
Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, Fry B, Meissner A, Wernig M, Plath K, Jaenisch R, Wagschal A, Feil R, Schreiber SL, Lander ES: A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006, 125: 315-326. 10.1016/j.cell.2006.02.041.
De Gobbi M, Garrick D, Lynch M, Vernimmen D, Hughes JR, Goardon N, Luc S, Lower KM, Sloane-Stanley JA, Pina C, Soneji S, Renella R, Enver T, Taylor S, Jacobsen SE, Vyas P, Gibbons RJ, Higgs DR: Generation of bivalent chromatin domains during cell fate decisions. Epigenetics Chromatin. 2011, 4: 9-10.1186/1756-8935-4-9.
Cimmino L, Abdel-Wahab O, Levine RL, Aifantis I: TET family proteins and their role in stem cell differentiation and transformation. Cell Stem Cell. 2011, 9: 193-204. 10.1016/j.stem.2011.08.007.
Figueroa ME, Abdel-Wahab O, Lu C, Ward PS, Patel J, Shih A, Li Y, Bhagwat N, Vasanthakumar A, Fernandez HF, Tallman MS, Sun Z, Wolniak K, Peeters JK, Liu W, Choe SE, Fantin VR, Paietta E, Lowenberg B, Licht JD, Godley LA, Delwel R, Valk PJ, Thompson CB, Levine RL, Melnick A: Leukemic IDH1 and IDH2 mutations result in a hypermethylation phenotype, disrupt TET2 function, and impair hematopoietic differentiation. Cancer Cell. 2010, 18: 553-567. 10.1016/j.ccr.2010.11.015.
Brenner C, Deplus R, Didelot C, Loriot A, Vire E, De Smet C, Gutierrez A, Danovi D, Bernard D, Boon T, Pelicci PG, Amati B, Kouzarides T, de Launoit Y, Di Croce L, Fuks F: Myc represses transcription through recruitment of DNA methyltransferase corepressor. EMBO J. 2005, 24: 336-346. 10.1038/sj.emboj.7600509.
Hervouet E, Vallette FM, Cartron PF: Dnmt3/transcription factor interactions as crucial players in targeted DNA methylation. Epigenetics. 2009, 4: 487-499. 10.4161/epi.4.7.9883.
Linhart HG, Lin H, Yamada Y, Moran E, Steine EJ, Gokhale S, Lo G, Cantu E, Ehrich M, He T, Meissner A, Jaenisch R: Dnmt3b promotes tumorigenesis in vivo by gene-specific de novo methylation and transcriptional silencing. Genes Dev. 2007, 21: 3110-3122. 10.1101/gad.1594007.
Ibrahim AE, Arends MJ, Silva AL, Wyllie AH, Greger L, Ito Y, Vowler SL, Huang TH, Tavare S, Murrell A, Brenton JD: Sequential DNA methylation changes are associated with DNMT3B overexpression in colorectal neoplastic progression. Gut. 2011, 60: 499-508. 10.1136/gut.2010.223602.
Nosho K, Shima K, Irahara N, Kure S, Baba Y, Kirkner GJ, Chen L, Gokhale S, Hazra A, Spiegelman D, Giovannucci EL, Jaenisch R, Fuchs CS, Ogino S: DNMT3B expression might contribute to CpG island methylator phenotype in colorectal cancer. Clin Cancer Res. 2009, 15: 3663-3671. 10.1158/1078-0432.CCR-08-2383.
Ogino S, Kawasaki T, Kirkner GJ, Loda M, Fuchs CS: CpG island methylator phenotype-low (CIMP-low) in colorectal cancer: possible associations with male sex and KRAS mutations. J Mol Diagn. 2006, 8: 582-588. 10.2353/jmoldx.2006.060082.
Kriaucionis S, Heintz N: The nuclear DNA base 5-hydroxymethylcytosine is present in Purkinje neurons and the brain. Science. 2009, 324: 929-930. 10.1126/science.1169786.
Nestor C, Ruzov A, Meehan R, Dunican D: Enzymatic approaches and bisulfite sequencing cannot distinguish between 5-methylcytosine and 5-hydroxymethylcytosine in DNA. Biotechniques. 2010, 48: 317-319. 10.2144/000113403.
Haffner MC, Chaux A, Meeker AK, Esopi DM, Gerber J, Pellakuru LG, Toubaji A, Argani P, Iacobuzio-Donahue C, Nelson WG, Netto GJ, De Marzo AM, Yegnasubramanian S: Global 5-hydroxymethylcytosine content is significantly reduced in tissue stem/progenitor cell compartments and in human cancers. Oncotarget. 2011, 2: 627-637.
Jin SG, Jiang Y, Qiu R, Rauch TA, Wang Y, Schackert G, Krex D, Lu Q, Pfeifer GP: 5-Hydroxymethylcytosine is strongly depleted in human cancers but its levels do not correlate with IDH1 mutations. Cancer Res. 2011, 71: 7360-7365. 10.1158/0008-5472.CAN-11-2023.
Nestor CE, Ottaviano R, Reddington J, Sproul D, Reinhardt D, Dunican D, Katz E, Dixon JM, Harrison DJ, Meehan R: Tissue-type is a major modifier of the 5-hydroxymethylcytosine content of human genes. Genome Res. 2011, 22: 467-477.
The R Project for Statistical Computing. [http://www.r-project.org/]
Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, Bunney WE, Myers RM, Speed TP, Akil H, Watson SJ, Meng F: Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucleic Acids Res. 2005, 33: e175-10.1093/nar/gni179.
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.
Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26: 841-842. 10.1093/bioinformatics/btq033.
Calvanese V, Fernandez AF, Urdinguio RG, Suarez-Alvarez B, Mangas C, Perez-Garcia V, Bueno C, Montes R, Ramos-Mejia V, Martinez-Camblor P, Ferrero C, Assenov Y, Bock C, Menendez P, Carrera AC, Lopez-Larrea C, Fraga MF: A promoter DNA demethylation landscape of human hematopoietic differentiation. Nucleic Acids Res. 2012, 40: 116-131. 10.1093/nar/gkr685.
Translational Cancer Research Group Antwerp - Epigenetics. [http://www.tcrg.be/en/page6/page13/epigenetics.html]
Habegger L, Sboner A, Gianoulis TA, Rozowsky J, Agarwal A, Snyder M, Gerstein M: RSEQtools: a modular framework to analyze RNA-Seq data using compact, anonymized data summaries. Bioinformatics. 2011, 27: 281-283. 10.1093/bioinformatics/btq643.
Lim E, Vaillant F, Wu D, Forrest NC, Pal B, Hart AH, Asselin-Labat ML, Gyorki DE, Ward T, Partanen A, Feleppa F, Huschtscha LI, Thorne HJ, Fox SB, Yan M, French JD, Brown MA, Smyth GK, Visvader JE, Lindeman GJ: Aberrant luminal progenitors as the candidate target population for basal tumor development in BRCA1 mutation carriers. Nat Med. 2009, 15: 907-913. 10.1038/nm.2000.
Catalogue of Somatic Mutations in Cancer - COSMIC. [http://www.sanger.ac.uk/perl/genetics/CGP/cosmic/]
Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, Beare D, Jia M, Shepherd R, Leung K, Menzies A, Teague JW, Campbell PJ, Stratton MR, Futreal PA: COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer. Nucleic Acids Res. 2011, 39: D945-950. 10.1093/nar/gkq929.
UCSC Genome Browser. [http://genome.ucsc.edu/]
Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, Ji H, Potash JB, Sabunciyan S, Feinberg AP: The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009, 41: 178-186. 10.1038/ng.298.
This work was funded by Breakthrough Breast Cancer and the Medical Research Council. We thank Professor WA Bickmore, Dr JG Prendergast, J Reddington (MRC Human Genetics Unit, University of Edinburgh), Professor CP Ponting and Dr C Nellaker (MRC Functional Genomics Unit, Oxford) for discussion of results and the manuscript.
The authors declare that they have no competing interests.
DS and RRK performed the research and analyzed data. DS, CEN, JMD, AHS, DJH, BHR and RRM helped design the research and review the manuscript. DS and RRM wrote the paper. All authors read and approved the final version of the manuscript for publication.