- Research article
- Open Access
- Open Peer Review
Cancer associated epigenetic transitions identified by genome-wide histone methylation binding profiles in human colorectal cancer samples and paired normal mucosa
BMC Cancervolume 11, Article number: 450 (2011)
Despite their well-established functional roles, histone modifications have received less attention than DNA methylation in the cancer field. In order to evaluate their importance in colorectal cancer (CRC), we generated the first genome-wide histone modification profiles in paired normal colon mucosa and tumor samples.
Chromatin immunoprecipitation and microarray hybridization (ChIP-chip) was used to identify promoters enriched for histone H3 trimethylated on lysine 4 (H3K4me3) and lysine 27 (H3K27me3) in paired normal colon mucosa and tumor samples from two CRC patients and for the CRC cell line HT29.
By comparing histone modification patterns in normal mucosa and tumors, we found that alterations predicted to have major functional consequences were quite rare. Furthermore, when normal or tumor tissue samples were compared to HT29, high similarities were observed for H3K4me3. However, the differences found for H3K27me3, which is important in determining cellular identity, indicates that cell lines do not represent optimal tissue models. Finally, using public expression data, we uncovered previously unknown changes in CRC expression patterns. Genes positive for H3K4me3 in normal and/or tumor samples, which are typically already active in normal mucosa, became hyperactivated in tumors, while genes with H3K27me3 in normal and/or tumor samples and which are expressed at low levels in normal mucosa, became hypersilenced in tumors.
Genome wide histone modification profiles can be used to find epigenetic aberrations in genes associated with cancer. This strategy gives further insights into the epigenetic contribution to the oncogenic process and may identify new biomarkers.
Cancer has been traditionally considered a genetic and cytogenetic disease, but recent years have brought epigenetics to the forefront of cancer research [1, 2]. Altered DNA methylation is nowadays considered a hallmark of neoplasia, including two different phenomena in cancer cells: global hypomethylation and CpG-promoter hypermethylation of tumor suppressor genes [1, 2]. The importance of epigenetic alterations in cancer is further highlighted by their use in diagnosis and by the development of new therapeutic strategies aiming at correcting them .
Although histone modifications play major roles in processes such as transcription, replication and DNA repair, their oncogenic importance is not yet well established. However, several lines of evidence suggest that alterations in histone modifications are crucial in cancer development and progression. Global changes in histone H4 modifications seem to be universal markers of malignant transformation , while other histone marks predict the prognosis of prostate cancer . Two of the most relevant histone modifications, both in general and from a cancer perspective, are histone H3 lysine 4 trimethylation (H3K4me3), which is found in promoters of active genes and histone H3 lysine 27 trimethylation (H3K27me3), which is preferentially associated with promoters of inactive genes . SMYD3, a histone methyl transferase specific for H3K4, is over expressed in colorectal, hepatic and breast cancers, suggesting that H3K4 hypermethylation can occur at promoters of oncogenes . MLL, another H3K4 methyltransferase, is frequently translocated in various forms of leukaemia  and the polycomb protein EZH2, a H3K27 specific histone methyl-transferase, is altered in multiple types of cancer . Furthermore, recent reports suggest that H3K27me3 mediated silencing of tumor suppressor genes is a frequent event in prostate cancer [10, 11].
Human genome-wide binding profiles for H3K4me3, H3K27me3 and other histone modifications have been previously generated, but in most cases, cancer cell lines were used as biological material [10, 12]. Although this offers important insights into the functionality of the histone marks, it might not be optimal when investigating their importance in oncogenesis, as exemplified by the discrepancy between H3K27me3 profiles in normal tissue and cell lines of the same origin , or by increased DNA methylation in embryonic stem cells (ESCs) due to in vitro culture . Therefore, in order to establish how chromatin is altered during the oncogenic process, histone modification profiles from normal and tumor tissue samples are most desirable. However, and to the best of our knowledge, this has only been reported in liver  and pheochromocytoma .
Colon adenocarcinomas display numerous epigenetic alterations, including hypermethylation of tumor suppressor genes and loss of imprinting at IGF2/H19 . Interestingly, the major risk factor for CRC is age and epigenetic lesions accumulate with aging and contribute to cell transformation . Despite the evident role of histone modifications in colon cancer, histone modification profiles have not been analyzed in coupled normal and tumor colon samples. Here we present the first genome-wide maps of H3K4me3 and H3K27me3 in paired normal and tumor samples from CRC patients. We identify previously unknown changes in the methylation statuses between the normal and tumor samples and correlate these to CRC-related pathway and function. Taken together, these data could be an important resource in understanding the epigenetic alterations associated to CRC.
H3K4me3 and H3K27me3 binding profiles in paired normal and tumor samples from CRC patients
We generated enrichment profiles for H3K4me3 and H3K27me3 in tumor samples from two CRC patients, for which the profiles in normal mucosa have previously been generated  and for the colon adenocarcinoma cell line HT29. After determining the enriched regions (Figure 1A) for each modification in normal and tumor samples (Table 1 and Additional file 1), they were annotated to human genes using the UCSC knownGenes databases  (Table 1) (Methods).
Signal footprints were created around transcription start sites and both the H3K4me3 and H3K27me3 patterns were similar to those reported [20–22] in other cell types (Figure 1B). As a quality control, microarray expression data from public repositories for normal colon, CRC tumors and HT29 were used . Although the expression data comes from unrelated patient material, the objective was to study genome wide trends in transcription given the epigenetic status and not single gene effects where expression from the same patient material is essential. The expression of H3K4me3 enriched genes were higher than average, while H3K27me3 genes tended to be silent (Figure 2A). Bivalent and especially semi-bivalent (Methods) gene groups included fewer genes and were poorly expressed, although with higher variation in expression. The quality of our data is also supported by previous validation of H3K4me3 and H3K27me3 enriched regions by qPCR in normal colon  demonstrating that ChIP-chip analysis can be performed using human tissue as biological material.
For regions with associated transcripts, gene ontology (GO) annotations were obtained and an analysis was performed for over/under representation of terms compared to a whole genome background (Figure 2B, Methods). The H3K4me3 and H3K27me3 gene groups have opposite functions and bivalent and semi-bivalent gene groups lie closer in functional space to H3K27me3 than H3K4me3 genes, as also suggested by expression data. Therefore, bivalent and semi-bivalent genes were treated as part of H3K27me3 gene groups for most subsequent analysis. In general, H3K4me3 fractions displayed over-representation in metabolic processes while H3K27me3 fractions showed enrichment in developmental categories. The overrepresented categories in H3K4me3 and H3K27me3 genes were very similar in the HT29 cell line.
H3K4me3 and H3K27me3 profiles are similar between normal and tumor samples but distinct compared to CRC cell lines
Using all probes we calculated the Pearson's correlation between all combinations of histone modifications, patient material (normal/tumor) and cell line (Additional file 2, Table S1). In all cases the correlation between any patient materials was always strongest when that material was compared to another patient material. In some cases, e.g. patient 1 tumor H3K27me3/K4me3, the correlation was highest with patient 2 normal H3K27me3/K4me3 rather than the HT29 cancer cell line. This pattern is more evident when regarding enriched regions only and by counting overlapping enriched regions, a relatively high similarity between normal and tumor samples and between patients for both H3K4me3 (Figure 3) and H3K27me3 was noted (Table 1 and Additional file 2, Figure S6). Around 60% of the H3K4me3 and 35% of H3K27me3 regions in both patients are common between normal and tumor tissue. When annotated transcripts rather than individual regions were used these numbers rise to around 70% and 50% respectively. Considering H3K4me3 genes, HT29 targets showed large overlaps with both patients normal and tumor samples (Table 1-2). The overlaps for the adenocarcinoma cell line SW480  were quite poor (10-20%) both with tissue samples and HT29 (Table 2). We also noted that the number (around 1%, 80-431) of bivalent regions detected in promoters here is lower than the number detected in a previous study  on differentiated cells (CDC36+) where around 3% (693) of the genes were found to have bivalently marked promoters. Cui et al  performed their study using next generation sequencing and called more regions enriched throughout the genome suggesting that the main cause of difference lies in stringency of the downstream analysis.
Detection of genes and pathways with chromatin state differences between tumor and normal samples
To identify the genomic regions differing in their chromatin state between normal and tumor tissues (Methods), the results were split in fifteen so called transitions each representing a combined state of one or more marks in normal and tumor. Three major groups were created; i) gaining or loosing a single mark, ii) keeping a mark and iii) changes involving both marks (Figure 4A and Additional file 2, Figure S4A). The majority fall into category ii), that is no change, and category i) which is in agreement with the large overlaps reported in the previous section (Additional file 2,Table S2).
Category iii) implies chromatin transitions between normal and tumor that can result in a switch from transcriptionally active (e.g. H3K4me3) to inactive states (e.g. H3K27me3) or vice verse. Therefore, it is the most interesting category in terms of potentially functional consequences from an oncogenic point of view. However, only 4 genes were common to both patients in category iii) (Additional file 2, Table S2). KLF7 (N K27 to T K4), a transcription factor belonging to the KLF family that plays critical roles in differentiation, development, and maintenance of tissue homeostasis , but, to the best of our knowledge, without a specific role in colon cancer. The remaining three common genes (EBF3, DKFZp667I0324 and RBMS1) were all bivalent in normal tissue but lost all their histone methylation in tumor. Among these genes, 2 have previously been functionally characterized. EBF3 is a transcription factor reported to be a tumor suppressor gene  with observed silencing in CRC cell line HCT116, was bivalent in normal and lost both marks in tumor. RBMS1 is a transcription factor that has been suggested as an important factor in inducing apoptosis in mice , and has been shown to have predictive power on colorectal cancer recurrence in human patient cohorts . A number of transitions were specific to individual patients e.g. from category iii). Some genes had H3K4me3 in normal colon which changed to H3K27me3 in tumors, e.g. in patient 1 HOXB13 and KREMEN2, which are frequently down regulated in CRC [30, 31]. Among genes with a silent state (H3K27me3 mark) in normal colon that became active in CRC samples (gained H3K4me3) there were several well-characterized oncogenes [32–36], e.g. in patient 1 FLI1, WWTR1 & ZEB2 and in patient 2 TACSTD2 & TWIST1 (Figure 1A and Additional file 2, Figure S5 (patient 1)), the last two with proposed importance in CRC [32, 35].
Several genes in categories i) and ii) were found to have the same transitions in both patients. Both our patients have H3K4me3 on the PTGER2 gene in normal tissue and gained H3K27me3 in the tumor, which offers a mechanistic explanation to the previously reported expression pattern  of this gene in colorectal cancer. In addition, three genes common to both patients which carried bivalent marks in normal and lost H3K4me3 in tumor, LHX9, PKNOX2 and LBXCOR1 and one, PRDM8, lost the H3K27me3 in tumor, all of which are transcription factors with unknown function in colon.
The chromatin states of genes involved in biological pathways with relevance in CRC biology  was specifically investigated. The WNT-signalling pathway has been implicated in CRC [39, 40], with the APC gene frequently silenced or mutated. Patient 1 has tumor repressive transitions in the APC gene and tumor activating transitions in Axin, PP2A, β-catenin and in many downstream factors along the WNT-signalling pathway affecting transcription, cell cycle and CRC development, especially within the TCF/LEF family [41, 42] (Figure 4B). This suggests a repression on the β-catenin degradation complex and an up regulation of β-catenin dependent gene expression with downstream effects on the cell cycle . Recently, a dual role of APC in the WNT-signalling pathway has been suggested with both activating and repressing mechanisms  which could explain the patterns found in patient 2 where APC has gained activating markers (Additional file 2, Figure S4B) In addition, patient 1 shows an activation of several known repressors of the WNT-pathway. CtBP represses E-cadherin which in turn is a tumor repressor restricting tumor cell motility and invasion . The ubiquitin ligase Cul1 is part of a mechanism targeting proteins for degradation and thus an important factor in cell cycle control. Specifically, Cul1 acts as an oncogene that forms complexes that specifically target the cell cycle inhibitor p27/Kip1 for degradation .
Global expression changes of H3K4me3 and H3K27me3 targets in CRC
Although epigenetic transitions within category iii) are expected to have important functional consequences, few examples are detected in a given CRC patient. Therefore, these examples are less likely to have large-scale changes in gene expression. More global relationships between chromatin state and gene expression were analyzed using all epigenetic transitions involving more than 100 genes i.e. maintenance of H3K4me3 or H3K27me3, tumor loss of H3K4me3 or H3K27me3 as well as tumor gain of H3K4me3 or H3K27me3. Microarray expression data for paired normal and tumor samples from 24 CRC patients  were used, and the expression changes between normal and cancer samples for the genes included in each transition category (patient specific) was compared (Figure 5 and Additional file 2, Figures S1-S3).
Genes in transitions where H3K4me3 was maintained or specifically gained in tumors, showed significantly higher expression levels in cancer samples (Figure 5A and Additional file 2, Figure S2). Opposite expression shifts were observed for genes that kept, gained or lost H3K27me3 in the tumor samples (Figure 5B, Additional file 2, Figure S2). The H3K4me3 and H3K27me3 gene groups also differed in the baseline expression in normal samples as H3K4me3 genes were highly expressed, while most H3K27me3 genes showed low expression levels. Furthermore, genes in H3K4me3 transition groups are typically already activated in normal tissue and become "hyperactivated" in tumors, while genes in H3K27me3 transition groups are lowly expressed in normal colon and become "hypersilenced" in cancer samples. The generality of the observations is supported by (i) the similar expression shifts obtained when considering the epigenetic transition groups from either patient 1 and patient 2 (Additional file 2, Figure S2A); and (ii) from a second expression data set with paired normal and tumor samples from 9 CRC obtained in a different array platform  (Additional file 2, Figure S2 B-C).
The significantly up- or down- regulated genes in cancer samples in the two expression microarray data sets described above assuming a paired sample design was analyzed further. For both patients and microarray data sets (Figure 5C and Additional file 2, Figure S3) genes keeping or gaining H3K4me3 in tumors were significantly over-represented among up-regulated genes and significantly under-represented among down-regulated genes. Opposite patterns were observed for genes that kept, gained or lost H3K27me3 in tumor samples (Figure 5C and, Additional file 2, Figure S3). These results indicate that at least for a fraction of genes the reported expression changes may be of large magnitude and occur in most patients. We gathered microarray expression data for HT29  and considered the same epigenetic transitions as described above, and can report that the tendencies towards hyperactivation or hypersilencing of H3K4me3 and H3K27me3 gene groups became even more pronounced in HT29 cells (Figure 5D-E and Additional file 2, Figure S4).
Hypersilencing of some H3K27me3 genes may involve loss of H3K27me3 and DNA hypermethylation
The expression levels between normal colon, tumor samples and HT29 cell line for all H3K4me3 and H3K27me3 enriched genes were compared in the three samples types and for those genes that were shared or unique to any of them. The overall H3K4me3 or H3K27me3 gene groups showed progressive hyperactivation and hypersilencing in tumor and HT29 cells compared to normal colon, respectively (Figure 6A,B). This was also true for the genes with shared histone modification states, but less obvious in some of the unique groups. For the H3K27me3 unique groups, the hypersilencing pattern was still found for those genes carrying that histone mark only in normal colon samples (Figure 6B).
It has been recently reported  that genes enriched in H3K27me3 in normal prostate frequently become DNA hypermethylated and lose H3K27me3 in prostate cancer cell lines. Furthermore, DNA hypermethylation seems to be exacerbated by in vitro culture . In order to test whether genes with the H3K27me3 mark in normal colon overlapped more with genes hypermethylated in CRC than genes that carried the histone mark in tumor samples, the lists of H3K27me3 genes was compared with those previously reported as frequently hypermethylated in CRC . As shown in Figure 6C, both H3K27me3 genes in normal and tumor colon samples were significantly over-represented among DNA hypermethylated genes. This was not the case for the HT29 genes. Over-representation was only significant for normal colon, when considering H3K27me3 unique groups, and in fact even a trend towards under-representation was observed for H3K27me3 genes unique to HT29 (Figure 6C).
By comparing normal and tumor colon samples, it was observed that some regions, either shared or unique to one patient, showed transitions in the chromatin states that would suggest major functional consequences (i.e. from normal H3K4me3 to tumor H3K27me3 or vice verse). This resembles the results obtained after genetic mutation analysis of all coding sequences in CRC patients, where very few mutations are frequently shared between patients . These observations suggest that in a given CRC patient, there could be a few major alterations (genetic and/or epigenetic) that can be crucial for tumorogenesis e.g. key components that could disturb the WNT-pathway such as silencing of the APC gene and activation of the β-catenin gene or up-regulation of repressors of tumor suppressing genes found in patient 1.
These few major and likely causative alterations could, however, have major impacts on global gene expression programs and chromatin states. Our observations suggest that H3K4me3 genes, typically active in normal colon, become "hyperactivated" in CRC. Since many of these genes are involved in basic cellular processes, their increased expression could merely be the result of higher metabolic demands and proliferation rates in tumor cells . Mechanistically, a good candidate to mediate this "hyperactivation" response is SMYD3, a H3K4 methyltransferase frequently over-expressed in CRC and with oncogenic function . On the other hand, H3K27me3 genes are already silent or lowly expressed in normal colon mucosa, and seemed to become "hypersilenced" in CRC. Down-regulation of H3K27me3 genes could be explained by polycomb-premarking of cancer-specific DNA hypermethylation [13, 51] and subsequent loss of H3K27me3 and reduced epigenetic plasticity . Due to tumor heterogeneity between different CRC patients, it is not expected that all normal colon H3K27me3 genes will loose the mark and become hypermethylated in a given tumor, which could explain why even genes that kept H3K27me3 state in our two patients seem to be hypersilenced when considering larger patient cohorts. Finally, silencing of genes showing specific gain of H3K27me3 in tumors has been also reported in prostate cancer [10, 11]. Mechanistically, the silencing of all these H3K27me3 gene groups may involve EZH2, which is frequently over-expressed in CRC and other cancer types  and is able to interact with and recruit DNA methyltransferases .
Another important aspect of the work is the evaluation of how well cell lines may represent the chromatin states of normal or tumor tissues. According to our data, for H3K27me3, which is important in determining cellular identity , cell lines could give a misleading picture of chromatin states in tissues they are supposed to represent. Such discrepancy seemed to increase depending on the differentiation grade of the cell line, suggesting that low-grade cell lines could offer better cellular models than high grade poorly differentiated cell lines when studying chromatin states. Furthermore, our results and previous observations suggest that such alterations in H3K27me3 states can arise due to a combination of polycomb-premarking of DNA hypermethylation, subsequent lost of H3K27me3 and exacerbated DNA hypermethylation due to in vitro culture [12–14, 51].
One such example is the polycomb binding profiles that were previously generated for the CRC cell line SW480 and used as a model of differentiated cells. As compared to ESCs, SW480 polycomb (H3K27me3) targets were not enriched in developmental related processes . Here, and previously , we have shown that H3K27me3 marked genes, both in normal and tumor colon tissues, are largely enriched in developmental processes, and in fact, they significantly overlapped with ESC H3K27me3 targets. Therefore, we wanted to investigate how well CRC cell lines represent the chromatin states of colon tissue samples (Table 2). This might be expected since many H3K4me3 genes correspond to housekeeping genes typically active in any given cell type. On the other hand, when considering H3K27me3 gene groups, the overlaps between normal and tumor samples were quite high (40-60%), while somewhat lower for HT29 (20-30%), but as described above, H3K27me3 target genes in this cell line were still largely related to developmental processes. Moreover, the fact that the HT29 H3K27me3 profile, as compared to SW480, shows higher overlap with tissues and is enriched in developmental processes, can be explained by the fact that HT29 is a low grade CRC cell line (grade II), while SW480 is a high grade (grade III-IV), poorly differentiated and highly abnormal CRC cell line.
In this work we generated the first coupled normal-tumor histone modification profiles in CRC, one of the cancer types with highest incidence in developed countries, that should provide a valuable resource for future chromatin studies with CRC and/or gastrointestinal focus. By studying histone modification profiles we show that it is possible to get new insights into the epigenetic process and discover new biomarkers for cancer. Two such candidates discovered here are the PTGER2 and KLF7 genes. The former gained H3K37me3 in tumor in both our patients and the latter changed its modification status from H3K4me3 to H3K27me3 in both patients thus changing from activated to repressed state. PTGER2 has previously been implicated in CRC but further research will be needed to clarify the role of KLF7.
Normal colon mucosa  and tumor samples were obtained from two patients diagnosed with CRC at Uppsala University Hospital (Sweden). Immediately upon surgery, the removed material was kept on ice. Samples from the normal colon mucosa were taken at least 15 cm from the localized tumor, while the tumor sample was collected without containing parts of the normal mucosa. A gastrointestinal pathologist performed sample collection and tumor evaluation. The patients did not receive any anticancer treatment before surgery. All patients treated for CRC in Uppsala are registered in the Swedish Colorectal Cancer Registry. When admitted to the hospital they are asked to participate in that registry and also informed that biopsies will be saved in a biobank. After having read the information the informed consent is verbally obtained. Linked to this we have an ethical approval to perform studies like this, #2006/077.
ChIP and ChIP-chip
ChIP on patient samples and HT29 cells were performed as previously described [13, 54]. DNA amplification, fragmentation, labelling and hybridizations of ChIP and input DNAs were performed according to Affymetrix recommendations and basically as previously described , using Affymetrix GeneChip Human Promoter 1.0 arrays, which cover approximately from 7.5 kb upstream to 2.45 kb downstream of transcription start sites for over 25,500 human promoters. Raw array data has been deposited in ArrayExpress under accession number E-TABM-533
Data analysis pipeline
Affymetrix Tiling Analysis SDK revision 4 was used to quantile normalize each pair of replicate ChIP measurements together with the two input measurements and adjusted to have a median intensity of 200 [55, 56]. Replicates were combined using a sliding window of 301 bp (bandwidth = 150) assigning the median value of the replicates in the window to the center probe. From this point, all analysis was carried out using custom scripts in R . Log2-ratios of merged ChIP measurements over merged input measurements were used for further analysis. Enriched regions were defined using a Z-score based sliding window approach. For each probe on the array, a Z-score was calculated from the average of log2-ratios of probes within a window of 150 base pairs centred on the probe against the distribution of log2-ratios on the array.
Enriched regions were defined from enriched probes (Z-score >= 6) using the following criteria; (i) at least 2 enriched probes within the region, (ii) each enriched probe must be flanked on both sides by existing measurements within 110 base pairs (iii) the maximum gap between enriched probes within a region is 110 base pairs. Regions were prolonged with probes positioned within 75 base pairs of the probes defining the tails of the regions. Finally, regions were merged if the shortest distance was less than 1000 base pairs.
Bivalent regions were defined as regions sharing at least three probes between H3K4me3 and H3K27me3 regions out of which at least two were enriched. Regions classified as having changed their methylation status ('transitions') between normal and tumor where required to have at least three probes in common out of which two needed to be commonly labelled as enriched in the different states. In addition to this we required a Z-score of at least 3 between the normal and tumor in the non-changed methylation status ensuring that we are not selecting regions immediately on opposite sides of the chosen cut-off, i.e. clearly separating changes involving bivalent marks to single marks from instances where a single mark is changed to another single mark.
We associated enriched regions to UCSC knownGene transcripts if; (i) the regions were intragenic and positioned within 2.5 kb downstream of TSS or (ii) intergenic and positioned within 2.5 kb upstream of TSS but not associated to any other transcript through intragenic positioning. This window size (+/- 2.5 kb) was chosen based on the distribution of both H3K4me3 and H3K27me3 signals around transcription start sites. We defined transcripts to be marked by H3K4me3 if associated with at least one H3K4me3 region and no H3K27me3 region. The same was done for H3K27me3 associated transcripts. Transcripts associated with both H3K4me3 regions and H3K27me3 regions but not with bivalent regions were called semi-bivalent.
Gene Ontology (GO) analysis
GO terms were tested using a one-sided Fishers exact test. P-values were corrected for multiple hypotheses testing using Bonferroni correction. The hierarchical clustering was done using Euclidean distances and the Ward agglomeration method. All analysis was done using R 
Microarray Gene expression data sets
Microarray genes expression data sets were downloaded from Gene Expression Omnibus (GEO) public repositories. Expression data was obtained and analyzed using GEOquery and the 'limma' Bioconductor package in R . Expression data from paired tumor and normal colon mucosa samples from 24 and 9 different CRC patients corresponded to GSE10950 and GSE5364 datasets, respectively. HT29 expression data (GSM277543) was originally generated in parallel with the tissue samples from GSE10950 dataset.
histone H3 lysine 4 trimethylation
histone H3 lysine 27 trimethylation
chromatin immunoprecipitation and microarray hybridization.
Baylin SB, Ohm JE: Epigenetic gene silencing in cancer - a mechanism for early oncogenic pathway addiction?. Nature reviews. 2006, 6 (2): 107-116.
Esteller M: Epigenetics provides a new generation of oncogenes and tumour-suppressor genes. British journal of cancer. 2006, 94 (2): 179-183. 10.1038/sj.bjc.6602918.
Esteller M: Cancer epigenomics: DNA methylomes and histone-modification maps. Nat Rev Genet. 2007, 8 (4): 286-298. 10.1038/nrg2005.
Fraga MF, Ballestar E, Villar-Garea A, Boix-Chornet M, Espada J, Schotta G, Bonaldi T, Haydon C, Ropero S, Petrie K, Iyer NG, Perez-Rosado A, Calvo E, Lopez JA, Cano A, Calasanz MJ, Colomer D, Piris MA, Ahn N, Imhof A, Caldas C, Jenuwein T, Esteller M: Loss of acetylation at Lys16 and trimethylation at Lys20 of histone H4 is a common hallmark of human cancer. Nature genetics. 2005, 37 (4): 391-400. 10.1038/ng1531.
Seligson DB, Horvath S, Shi T, Yu H, Tze S, Grunstein M, Kurdistani SK: Global histone modification patterns predict risk of prostate cancer recurrence. Nature. 2005, 435 (7046): 1262-1266. 10.1038/nature03672.
Kouzarides T: Chromatin modifications and their function. Cell. 2007, 128 (4): 693-705. 10.1016/j.cell.2007.02.005.
Hamamoto R, Furukawa Y, Morita M, Iimura Y, Silva FP, Li M, Yagyu R, Nakamura Y: SMYD3 encodes a histone methyltransferase involved in the proliferation of cancer cells. Nature cell biology. 2004, 6 (8): 731-740. 10.1038/ncb1151.
Krivtsov AV, Armstrong SA: MLL translocations, histone modifications and leukaemia stem-cell development. Nature reviews. 2007, 7 (11): 823-833. 10.1038/nrc2253.
Simon JA, Lange CA: Roles of the EZH2 histone methyltransferase in cancer epigenetics. Mutat Res. 2008, 647 (1-2): 21-29. 10.1016/j.mrfmmm.2008.07.010.
Kondo Y, Shen L, Cheng AS, Ahmed S, Boumber Y, Charo C, Yamochi T, Urano T, Furukawa K, Kwabi-Addo B, Gold DL, Sekido Y, Huang TH, Issa JP: Gene silencing in cancer by histone H3 lysine 27 trimethylation independent of promoter DNA methylation. Nature genetics. 2008, 40 (6): 741-750. 10.1038/ng.159.
Yu J, Rhodes DR, Tomlins SA, Cao X, Chen G, Mehra R, Wang X, Ghosh D, Shah RB, Varambally S, Pienta KJ, Chinnaiyan AM: A polycomb repression signature in metastatic prostate cancer predicts cancer outcome. Cancer Res. 2007, 67 (22): 10657-10663. 10.1158/0008-5472.CAN-07-2498.
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.
Rada-Iglesias A, Enroth S, Andersson R, Wanders A, Pahlman L, Komorowski J, Wadelius C: Histone H3 lysine 27 trimethylation in adult differentiated colon associated to cancer DNA hypermethylation. Epigenetics. 2009, 4 (2): 107-113. 10.4161/epi.4.2.8038.
Meissner A, Mikkelsen TS, Gu H, Wernig M, Hanna J, Sivachenko A, Zhang X, Bernstein BE, Nusbaum C, Jaffe DB, Gnirke A, Jaenisch R, Lander ES: Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature. 2008, 454 (7205): 766-770.
Acevedo LG, Bieda M, Green R, Farnham PJ: Analysis of the mechanisms mediating tumor-specific changes in gene expression in human liver tumors. Cancer Res. 2008, 68 (8): 2641-2651. 10.1158/0008-5472.CAN-07-5590.
Sandgren J, Andersson R, Rada-Iglesias A, Enroth S, Akerstrom G, Dumanski JP, Komorowski J, Westin G, Wadelius C: Integrative epigenomic and genomic analysis of malignant pheochromocytoma. Exp Mol Med. 2010, 42 (7): 484-502. 10.3858/emm.2010.42.7.050.
Wong JJ, Hawkins NJ, Ward RL: Colorectal cancer: a model for epigenetic tumorigenesis. Gut. 2007, 56 (1): 140-148. 10.1136/gut.2005.088799.
Fraga MF, Agrelo R, Esteller M: Cross-talk between aging and cancer: the epigenetic language. Annals of the New York Academy of Sciences. 2007, 1100: 60-74. 10.1196/annals.1395.005.
Karolchik D, Kuhn RM, Baertsch R, Barber GP, Clawson H, Diekhans M, Giardine B, Harte RA, Hinrichs AS, Hsu F, Kober KM, Miller W, Pedersen JS, Pohl A, Raney BJ, Rhead B, Rosenbloom KR, Smith KE, Stanke M, Thakkapallayil A, Trumbower H, Wang T, Zweig AS, Haussler D, Kent WJ: The UCSC Genome Browser Database: 2008 update. Nucleic acids research. 2008, D773-779. 36 Database
Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K: High-resolution profiling of histone methylations in the human genome. Cell. 2007, 129 (4): 823-837. 10.1016/j.cell.2007.05.009.
Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, Margulies EH, Weng ZP, Snyder M, Dermitzakis ET, Stamatoyannopoulos JA, Thurman RE, Kuehn MS, Taylor CM, Neph S, Koch CM, Asthana S, Malhotra A, Adzhubei I, Greenbaum JA, Andrews RM, Flicek P, Boyle PJ, Cao H, Carter NP, Clelland GK, Davis S, Day N, Dhami P, Dillon SC, Dorschner MO, et al: Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007, 447 (7146): 799-816. 10.1038/nature05874.
Pan G, Tian S, Nie J, Yang C, Ruotti V, Wei H, Jonsdottir GA, Stewart R, Thomson JA: Whole-genome analysis of histone H3 lysine 4 and lysine 27 methylation in human embryonic stem cells. Cell stem cell. 2007, 1 (3): 299-312. 10.1016/j.stem.2007.08.003.
Jiang X, Tan J, Li J, Kivimae S, Yang X, Zhuang L, Lee PL, Chan MT, Stanton LW, Liu ET, Cheyette BN, Yu Q: DACT3 is an epigenetic regulator of Wnt/beta-catenin signaling in colorectal cancer and is a therapeutic target of histone modifications. Cancer cell. 2008, 13 (6): 529-541. 10.1016/j.ccr.2008.04.019.
Squazzo SL, O'Geen H, Komashko VM, Krig SR, Jin VX, Jang SW, Margueron R, Reinberg D, Green R, Farnham PJ: DACT3 is an epigenetic regulator of Wnt/beta-catenin signaling in colorectal cancer and is a therapeutic target of histone modifications. Genome Res. 2006, 16 (7): 890-900. 10.1101/gr.5306606.
Cui K, Zang C, Roh TY, Schones DE, Childs RW, Peng W, Zhao K: Chromatin signatures in multipotent human hematopoietic stem cells indicate the fate of bivalent genes during differentiation. Cell stem cell. 2009, 4 (1): 80-93. 10.1016/j.stem.2008.11.011.
McConnell BB, Yang VW: Mammalian Kruppel-Like Factors in Health and Diseases. Physiol Rev. 2010, 90 (4): 1337-1381. 10.1152/physrev.00058.2009.
Zhao LY, Niu Y, Santiago A, Liu J, Albert SH, Robertson KD, Liao D: An EBF3-mediated transcriptional program that induces cell cycle arrest and apoptosis. Cancer Res. 2006, 66 (19): 9445-9452. 10.1158/0008-5472.CAN-06-1713.
Nomura J, Matsumoto K, Iguchi-Ariga SMM, Ariga H: Positive regulation of Fas gene expression by MSSP and abrogation of Fas-mediated apoptosis induction in MSSP-deficient mice. Exp Cell Res. 2005, 305 (2): 324-332. 10.1016/j.yexcr.2005.01.006.
Reeve AE, Lin YH, Friederichs J, Black MA, Mages J, Rosenberg R, Guilford PJ, Phillips V, Thompson-Fawcett M, Kasabov N, Toro T, Merrie AE, van Rij A, Yoon HS, McCall JL, Siewert JR, Holzmann B: Multiple gene expression classifiers from different array platforms predict poor prognosis of colorectal cancer. Clinical Cancer Research. 2007, 13 (2): 498-507. 10.1158/1078-0432.CCR-05-2734.
Jung C, Kim RS, Zhang H, Lee SJ, Sheng H, Loehrer PJ, Gardner TA, Jeng MH, Kao C: HOXB13 is downregulated in colorectal cancer to confer TCF4-mediated transactivation. British journal of cancer. 2005, 92 (12): 2233-2239. 10.1038/sj.bjc.6602631.
Maehata T, Taniguchi H, Yamamoto H, Nosho K, Adachi Y, Miyamoto N, Miyamoto C, Akutsu N, Yamaoka S, Itoh F: Transcriptional silencing of Dickkopf gene family by CpG island hypermethylation in human gastrointestinal cancer. World J Gastroenterol. 2008, 14 (17): 2702-2714. 10.3748/wjg.14.2702.
Ohmachi T, Tanaka F, Mimori K, Inoue H, Yanaga K, Mori M: Clinical significance of TROP2 expression in colorectal cancer. Clin Cancer Res. 2006, 12 (10): 3057-3063. 10.1158/1078-0432.CCR-05-1961.
Chan SW, Lim CJ, Guo K, Ng CP, Lee I, Hunziker W, Zeng Q, Hong W: A role for TAZ in migration, invasion, and tumorigenesis of breast cancer cells. Cancer Res. 2008, 68 (8): 2592-2598. 10.1158/0008-5472.CAN-07-2696.
Truong AH, Ben-David Y: The role of Fli-1 in normal cell function and malignant transformation. Oncogene. 2000, 19 (55): 6482-6489. 10.1038/sj.onc.1204042.
Valdes-Mora F, Gomez Del Pulgar T, Bandres E, Cejas P, Ramirez de Molina A, Perez-Palacios R, Gallego-Ortega D, Garcia-Cabezas MA, Casado E, Larrauri J, Nistal M, Gonzalez-Baron M, Garcia-Foncillas J, Lacal JC: TWIST1 Overexpression is Associated with Nodal Invasion and Male Sex in Primary Colorectal Cancer. Ann Surg Oncol. 2008
Comijn J, Berx G, Vermassen P, Verschueren K, van Grunsven L, Bruyneel E, Mareel M, Huylebroeck D, van Roy F: The two-handed E box binding zinc finger protein SIP1 downregulates E-cadherin and induces invasion. Mol Cell. 2001, 7 (6): 1267-1278. 10.1016/S1097-2765(01)00260-X.
Gustafsson A, Hansson E, Kressner U, Nordgren S, Andersson M, Wang W, Lonnroth C, Lundholm K: EP1-4 subtype, COX and PPAR gamma receptor expression in colorectal cancer in prediction of disease-specific mortality. International journal of cancer. 2007, 121 (2): 232-240. 10.1002/ijc.22582.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic acids research. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.
Huang H, He X: Wnt/beta-catenin signaling: new (and old) players and new insights. Current opinion in cell biology. 2008, 20 (2): 119-125. 10.1016/j.ceb.2008.01.009.
Segditsas S, Tomlinson I: Colorectal cancer and genetic alterations in the Wnt pathway. Oncogene. 2006, 25 (57): 7531-7537. 10.1038/sj.onc.1210059.
Polakis P: Wnt signaling and cancer. Genes & development. 2000, 14 (15): 1837-1851.
Aoki M, Hecht A, Kruse U, Kemler R, Vogt PK: Nuclear endpoint of Wnt signaling: neoplastic transformation induced by transactivating lymphoid-enhancing factor 1. Proc Natl Acad Sci USA. 1999, 96 (1): 139-144. 10.1073/pnas.96.1.139.
Giles RH, van Es JH, Clevers H: Caught up in a Wnt storm: Wnt signaling in cancer. Biochimica et biophysica acta. 2003, 1653 (1): 1-24.
Takacs CM, Baird JR, Hughes EG, Kent SS, Benchabane H, Paik R, Ahmed Y: Dual Positive and Negative Regulation of Wingless Signaling by Adenomatous Polyposis Coli. Science. 2008, 319 (5861): 333-336. 10.1126/science.1151232.
Chinnadurai G: The Transcriptional Corepressor CtBP: A Foe of Multiple Tumor Suppressors. Cancer Research. 2009, 69 (3): 731-734. 10.1158/0008-5472.CAN-08-3349.
Lee J, Zhou P: Cullins and Cancer. Genes & Cancer. 2010, 1 (7): 690-699.
Yu K, Ganesan K, Tan LK, Laban M, Wu J, Zhao XD, Li H, Leung CH, Zhu Y, Wei CL, Hooi SC, Miller L, Tan P: A precisely regulated gene expression cassette potently modulates metastasis and survival in multiple solid cancers. PLoS genetics. 2008, 4 (7): e1000129-10.1371/journal.pgen.1000129.
Keshet I, Schlesinger Y, Farkash S, Rand E, Hecht M, Segal E, Pikarski E, Young RA, Niveleau A, Cedar H, Simon I: Evidence for an instructive mechanism of de novo methylation in cancer cells. Nature genetics. 2006, 38 (2): 149-153. 10.1038/ng1719.
Sjoblom T, Jones S, Wood LD, Parsons DW, Lin J, Barber TD, Mandelker D, Leary RJ, Ptak J, Silliman N, Szabo S, Buckhaults P, Farrell C, Meeh P, Markowitz SD, Willis J, Dawson D, Willson JK, Gazdar AF, Hartigan J, Wu L, Liu C, Parmigiani G, Park BH, Bachman KE, Papadopoulos N, Vogelstein B, Kinzler KW, Velculescu VE: The consensus coding sequences of human breast and colorectal cancers. Science. 2006, 314 (5797): 268-274. 10.1126/science.1133427.
Jin S, DiPaola RS, Mathew R, White E: Metabolic catastrophe as a means to cancer cell death. J Cell Sci. 2007, 120 (Pt 3): 379-383.
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. Nature genetics. 2007, 39 (2): 157-158. 10.1038/ng1941.
Vire E, Brenner C, Deplus R, Blanchon L, Fraga M, Didelot C, Morey L, Van Eynde A, Bernard D, Vanderwinden JM, Bollen M, Esteller M, Di Croce L, de Launoit Y, Fuks F: The Polycomb group protein EZH2 directly controls DNA methylation. Nature. 2006, 439 (7078): 871-874. 10.1038/nature04431.
Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, Alvarez P, Brockman W, Kim TK, Koche RP, Lee W, Mendenhall E, O'Donovan A, Presser A, Russ C, Xie X, Meissner A, Wernig M, Jaenisch R, Nusbaum C, Lander ES, Bernstein BE: Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 2007, 448 (7153): 553-560. 10.1038/nature06008.
Rada-Iglesias A, Wallerman O, Koch C, Ameur A, Enroth S, Clelland G, Wester K, Wilcox S, Dovey OM, Ellis PD, Wraight VL, James K, Andrews R, Langford C, Dhami P, Carter N, Vetrie D, Ponten F, Komorowski J, Dunham I, Wadelius C: Binding sites for metabolic disease related transcription factors inferred at base pair resolution by chromatin immunoprecipitation and genomic microarrays. Hum Mol Genet. 2005, 14 (22): 3435-3447. 10.1093/hmg/ddi378.
Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.
Tiling Array Software. [http://www.affymetrix.com/partners_programs/programs/developer/TilingArrayTools/index.affx]
R Development Core Team: R: A Language and Environment for Statistical Computing. 2011
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/11/450/prepub
AR-I was partially supported by Signe and Olof Wallenius Foundation and Lars Hiertas Foundation. JK was supported by the Polish Ministry of Science and Higher Education, grant number N301 239536 and Uppsala University. CW was supported by the Swedish Research Council for Science and Technology, the Cancer Foundation, the Markus Borgström Foundation and Uppsala University (ALF).
SE and ARI analysed the data, SE and RA designed and implemented the bioinformatics pipeline, ARI and OW prepared and performed the ChIP, AW and LP prepared the patient material, JK supervised the bioinformatics development. CW conceived and designed the study. All authors read and approved the final manuscript.
Stefan Enroth, Alvaro Rada-Iglesisas contributed equally to this work.