| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Cell, Tumor, and Stem Cell Biology |
Divisions of 1 Cancer Biology and 2 Biometry and Clinical Trials, The Johns Hopkins Kimmel Cancer Center, Baltimore, Maryland; 3 Department of Molecular Biotechnology, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium; and 4 OncoMethylome Sciences, SA, Liege, Belgium
Requests for reprints: Stephen B. Baylin, The Sidney Kimmel Comprehensive Cancer Center at Johns Hopkins, Suite 541, 1650 Orleans Street, Baltimore, MD 21231. Phone: 410-955-8506; Fax: 410-614-9984; E-mail: sbaylin{at}jhmi.edu.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
We have previously shown in local promoter chromatin immunoprecipitation (ChIP) studies that several gene promoters retain many modifications that are typically associated with repression of transcription on DNA demethylation–induced gene expression in double-knockout cells (4). These genes, now actively transcribed, exist in what we have termed a semi-heterochromatic state (4). This state shares similarities with what others have described as bivalent promoter domains that contain both the active H3K4me2 mark and the repressive H3K27me3 mark (5). This chromatin status is thought to maintain a set of embryonic genes, which contain promoters with CpG islands, in a poised, low transcription state, which facilitates prevention of premature lineage commitment (5).
In the present study, we have sought to more globally establish whether what we have learned from the above chromatin analyses of subsets of colon cancer genes may define key features of the cancer epigenome. We have thus performed a promoter genome-wide localization, by ChIP-chip analyses, of the above key active (H3K4me2) and repressive (H3K27me3) marks in wild-type HCT116 colorectal cancer cells compared with double-knockout cells lacking DNA methylation. In this study, we are able to orient our global results around control genes from the previous local ChIP studies (4) mentioned above. Moreover, we have mapped our ChIP-chip results to previous microarray expression results used in an approach to randomly scan colon cancer cells for CpG island DNA-hypermethylated genes (6). This has allowed us to correlate the chromatin status of a large number of verified DNA-hypermethylated genes in HCT116 and double-knockout cells for which both reverse transcription-PCR (RT-PCR) and microarray expression analyses had previously been matched (6) to globally determine the chromatin signature of CpG island–containing genes before and after removal of DNA hypermethylation.
| Materials and Methods |
|---|
|
|
|---|
ChIP-chip. ChIP was combined with DNA microarray analysis results (6) on HCT116 colorectal cancer cells and double-knockout cells in duplicate using antibodies specific for H3K27me3 (7), H3K4me2 (Upstate), or an IgG control (Upstate). Using 5 x 107 cells per experiment, ChIP, DNA amplification, labeling, and array hybridization were done as previously described (8). Independent batches of HCT116 and double-knockout cells were used to perform independent ChIP experiments, and ChIPs for each antibody were hybridized to independent Agilent array sets. Whole genome expanded promoter 244K arrays from Agilent Technologies that span roughly 5 kb upstream to 2 kb downstream from gene transcription start sites were used.
Analyses. The data preprocessing for the tiling arrays was done according to the manufacturer's recommendations using the Agilent ChIP Analytics 1.3 software (i.e., a median blanks subtraction followed by an interarray median normalization and a dye-bias median normalization). The data were imported into R for downstream analysis. The expression arrays were preprocessed as previously described (6).
The expression data were analyzed as described in ref. 6, but some modifications were made to determine the basal expression of genes. Expression values were calculated using an Agilent whole human genome microarray with the double-knockout sample labeled with Cy5 and the original HCT116 sample labeled with Cy3. Expression status was defined as the natural logarithm of the normalized single channel signal (i.e., the Cy5 channel for double-knockout cells and the Cy3 channel for the HCT116 cells). The upper boundary of the silent expression range was chosen so that 90% (38-42) of previously verified DNA-hypermethylated genes (6) were contained within this zone. This cutoff was also chosen because it clearly separates the general expression trends of the verified genes from a set of false-positive genes (6) as shown in Supplementary Fig. S1.
Composite graphs were created to reflect the general chromatin trends of a particular group of genes. To retain maximal data quality, only genes that had a minimal amount of probes within the promoter-associated region around the transcription start site were allowed to contribute data points to the composite profile. A sliding window of 100 bp was normally used to compute the average H3K27me3 and H3K4me2 profiles for both HCT116 and double-knockout cells. An extended window size was used in case the number of input genes was lower than 50 to smooth the profiles and give a better view of the general trends. The profiles made with an extended window size were compared with the standard 100-bp window profiles to ensure that only the visible noise was reduced and nothing about the trends or gene characteristics was changed. Finally, the mean H3K27me3 enrichment in a broad region from –2,250 to +2,250 relative to the TSS was determined for all genes on the tiling array and was used to visualize the general H3K27me3 trends.
| Results and Discussion |
|---|
|
|
|---|
|
|
1.0 and ranging to 5.0) of H3K27me3 enrichment over the entire promoter region (Fig. 1D and E; Supplementary Figs. S3 and S4). These results match well with our previous studies linking, in local ChIP analyses (matches of such local ChIP data with patterns of ChIP-chip results for five individual genes shown in Supplementary Fig. S3), the association of H3K27me3 enrichment with CpG island DNA-hypermethylated genes in cancer (4). Interestingly, we found, however, that the H3K27me3 enrichment at DNA-hypermethylated genes is actually intermediate when compared with other PcG-marked genes. This is illustrated by a comparison of these DNA-hypermethylated genes to a group of 815 genes with very high H3K27me3 enrichment (Fig. 1F, for average patterns, and for individual genes see Supplementary Fig. S5—gene list will be included online). The difference in H3K27me3 enrichment at hypermethylated genes is highly statistically significant when compared with levels at either active or high PcG genes (P < 2e–16, Wilcoxon rank sum test). Interestingly, this high H3K27me3 enrichment selects for those genes lacking CpG islands (P < 2.2e–6, Fisher's exact test). Of the 815 genes selected for this highest H3K27me3 enrichment, 698 (86%) lack such islands in the promoter region. However, the pattern of broad distribution around the annotated start sites of the genes seems to be similar to that for the DNA-hypermethylated genes (Supplementary Fig. S5).
Having defined the above gene expression and chromatin parameters, we sought to characterize essentially all of the candidate promoter CpG island DNA-hypermethylated genes in HCT116 colon cancer cells. As shown schematically in Fig. 2A , our previous discovery approach using microarray expression analysis has discovered hypermethylated genes by defining a zone of increased expression following either genetic disruption of two DNA methyltransferases, DNMT1 and 3b (double-knockout cells), or treatment of the cells with the DNMT inhibitor 5-aza-2'-deoxycytidine (5-aza-dC; ref. 6). Further, we excluded those genes, which increased expression following treatment of the HCT116 cells with trichostatin A (TSA), an inhibitor of histone deacetylases (6). Such refractoriness to TSA is a well-defined feature of genes with densely DNA-methylated promoter CpG islands (18, 19). Finally, only genes with low basal expression on the expression microarray (previously using a cutoff that would include genes currently characterized in both the intermediate and very low expression ranges shown in Fig. 1A) were considered as possible candidates for silencing by DNA hypermethylation.
We now show, through our refined expression state classification in Fig. 1A, that many of the genes that ended up as being false positives for DNA hypermethylation in our previous study (6) actually exhibit chromatin characteristics of active genes. We examined a group of more than 1,792 genes in HCT116 cells (gene list will be included online) whose expression behavior on microarray analysis (increase with 5-aza-dC and no increase with TSA, regardless of basal expression) placed them in the top tier and next tier candidate zone (zones 1 and 2 in Fig. 2A). Surprisingly, the chromatin of 515 genes in this tier, which had intermediate expression on the microarray (Fig. 1A), have distinctly more H3K4me2 enrichment around the transcription start site when compared with a group of 610 genes within the very low expression region occupied by the verified DNA-hypermethylated genes (Fig. 2B). Thus, these genes with intermediate expression levels have an active chromatin pattern for this mark (Fig. 2B).
To place these findings into perspective for the purpose of enriching the efficiency of our expression array paradigm to randomly identify CpG island DNA-hypermethylated genes, we examined the H3K4me2 mark for 23 genes previously identified in our recent study (6) as being false positives for promoter CpG island DNA methylation. The leading characteristic of these false-positive genes in our previous study was that they had a basal level of expression by RT-PCR rather than the absent mRNA signal that characterized all of the genes verified to have dense promoter DNA hypermethylation (6). Indeed, we now find that these 23 genes have an intermediate level of expression on microarray analysis (Supplementary Fig. S1). Moreover, by ChIP-chip analyses, these false-positive genes have a H3K4me2 pattern similar to those for active genes (Fig. 2C).
Strikingly, having now eliminated the above 515 genes with intermediate basal expression and enriched H3K4me2 from the candidate list, the chromatin pattern for the remaining 610 genes within the silent zone (gene list in Supplementary Table S2) shows an identical chromatin pattern to that of the 42 verified DNA-hypermethylated genes. Individual patterns for selected genes from the 610 are shown in Supplementary Fig. S6. The pattern includes a broad distribution of H3K27me3 around the gene promoters and a low H3K4me2 peak positioned directly over the transcription start site without the dip for the nucleosome-free region (Fig. 2D). Thus, the level of H3K4me2, matched with a very stringent level for low basal transcription, can potentially eliminate false positives, eliminating some 28.75% of genes in the top tier (zone 1 in Fig. 2A) and 51.09% of the next tier (zone 2 in Fig. 2A) in the no TSA response zone that are not truly DNA methylated. We have previously identified that verification of genes as DNA methylated and silenced in the cell lines on which the microarray approach is used is
80% in the top tier and 50% in the next tier (6). Thus, knowledge of the chromatin status of the genes in cancer cells can markedly increase the efficiency of our expression array approach for identification of genes with promoter CpG island DNA methylation and complete transcriptional silencing.
Our previous chromatin studies at a limited number of DNA-hypermethylated genes suggested that even when these genes are induced to reexpress by DNA demethylation, they do not return to a highly active state but rather to one of low, poised transcription with exhibition of the bivalent chromatin pattern (4). In stem/progenitor cells, this pattern is thought to hold a key group of genes in a poised, low transcription state, necessary to maintain pluripotency (5). We were now in a position to test whether this is truly a universal property of CpG island DNA-hypermethylated genes in colon cancer cells. We examined this by comparing HCT116 cells to their isogenic derivative double-knockout cells, in which all tested DNA-hypermethylated genes are demethylated and reexpressed (3, 20).
We first find that the expression of the 42 verified DNA-hypermethylated genes is distinctly increased to an intermediate expression range in double-knockout cells (Fig. 1A), verifying our RT-PCR results and microarray analysis in our previous study (6). We then examined the chromatin profile of the genes in the double-knockout cells, shown for example for the individual SFRP1 gene promoter in Fig. 3A and B ; for other individual genes in Supplementary Figs. S3, S4, and S6; and for the composite profile of the validated hypermethylated genes in Fig. 3C (matching of previous local ChIP data to ChIP-chip data for five individual DNA-hypermethylated genes is shown in Supplementary Fig. S3). These promoters distinctly adopt a bivalent chromatin pattern on DNA demethylation, characterized by an increase in H3K27me3 enrichment near the gene start sites (P = 0.009873, Wilcoxon rank sum test; Fig. 1F), with a simultaneous increase in H3K4me2 levels (P = 5.865e–07, Wilcoxon rank sum test). We next examined the entire list of 610 best candidate DNA-hypermethylated genes from the TSA-negative zone (zones 1 and 2 in Fig. 2A), as defined by the stringent basal expression criteria above. These genes also assume an intermediate expression range in the double-knockout cells (Fig. 1A) and, again, adopt a bivalent pattern chromatin pattern (Fig. 3D). It is particularly interesting that the H3K4me2 mark not only increases but also redistributes the peak positions to either side of the direct transcription start site.
|
|
In summary, our global genome studies of two important histone modification marks in colon cancer cells have taught us much about a key component of abnormal epigenetic gene regulation in cancer cells, gene promoter CpG island DNA methylation, and associated tight transcriptional silencing. Our findings firmly suggest that DNA methylation is superimposed on a bivalent chromatin state, which has been best associated with the CpG islands of an important subset of genes in embryonic stem and committed progenitor cells. These global findings further link the origins of aberrantly DNA-hypermethylated genes in cancer cells to an underlying chromatin pattern similar to that of an important group of CpG island–containing, low-expression genes in stem/precursor cells. We can now provide additional support for these findings by the fact that 48.69% (P < 2.2e–16, Fisher's exact test) of the 610 best candidate DNA-hypermethylated genes, as we now define them in this study, are listed as PcG marked in published tiling studies of embryonic stem cells, embryonic fibroblasts, mouse embryonic fibroblasts, or neural progenitor cells (21–23). This accentuates the likelihood that this precursor cell state, which normally helps hold genes in a low-transcription poised state, may predispose such genes to DNA methylation and conversion of this state to very tight, heritable, transcriptional silencing in adult cancers. As we have discussed, this chromatin pattern may reflect the early cell compartments from which cancers arise (24). Loss of function of many of the genes involved may help cancers retain properties of these early cells at the expense of normal differentiation. Finally, our chromatin findings help explain how CpG island DNA methylation must be relieved to allow for reexpression via histone deacetylase inhibition. This finding may help refine the concepts underlying ongoing clinical trials for combined use of DNA-demethylating and histone deacetylase inhibitors in treating cancer.
| Disclosure of Potential Conflicts of Interest |
|---|
|
|
|---|
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| Footnotes |
|---|
K.M. McGarvey and L. Van Neste contributed equally to this work.
Received 2/26/08. Revised 4/10/08. Accepted 5/16/08.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. Zeng and J. Li Maximization of negative correlations in time-course gene expression data for enhancing understanding of molecular pathways Nucleic Acids Res., October 23, 2009; (2009) gkp822v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Jin, B. Yao, J.-L. Li, C. R. Fields, A. L. Delmas, C. Liu, and K. D. Robertson DNMT1 and DNMT3B Modulate Distinct Polycomb-Mediated Histone Modifications in Colon Cancer Cancer Res., September 15, 2009; 69(18): 7412 - 7421. [Abstract] [Full Text] [PDF] |
||||
![]() |
P.-Y. Hsu, D. E. Deatherage, B. A.T. Rodriguez, S. Liyanarachchi, Y.-I Weng, T. Zuo, J. Liu, A. S.L. Cheng, and T. H-M. Huang Xenoestrogen-Induced Epigenetic Repression of microRNA-9-3 in Breast Epithelial Cells Cancer Res., July 15, 2009; 69(14): 5936 - 5945. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. C. Glockner, M. Dhir, J. M. Yi, K. E. McGarvey, L. Van Neste, J. Louwagie, T. A. Chan, W. Kleeberger, A. P. de Bruine, K. M. Smits, et al. Methylation of TFPI2 in Stool DNA: A Potential Novel Biomarker for the Detection of Colorectal Cancer Cancer Res., June 1, 2009; 69(11): 4691 - 4699. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Hussain, M. Rao, A. E. Humphries, J. A. Hong, F. Liu, M. Yang, D. Caragacianu, and D. S. Schrump Tobacco Smoke Induces Polycomb-Mediated Repression of Dickkopf-1 in Lung Cancer Cells Cancer Res., April 15, 2009; 69(8): 3570 - 3578. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Rodriguez, M. Munoz, L. Vives, C. G. Frangou, M. Groudine, and M. A. Peinado Bivalent domains enforce transcriptional memory of DNA methylated genes in cancer cells PNAS, December 16, 2008; 105(50): 19809 - 19814. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |