Abstract
The formation of differentiated cell types from pluripotent progenitors involves epigenetic regulation of gene expression. DNA hydroxymethylation results from the enzymatic oxidation of 5-methylcytosine (5-mC) to 5-hydroxymethylcytosine (5-hmC) by the ten-eleven translocation (TET) 5-mC dioxygenase enzymes. Previous work has mapped changes in 5-mC during differentiation of intestinal stem cells. However, whether or not 5-hmC regulates colonocyte differentiation is unknown. Here we show that 5-hmC regulates gene expression during colonocyte differentiation and controls gene expression in human colon cancers. Genome-wide profiling of 5-hmC during in vitro colonic differentiation demonstrated that 5-hmC is gained at highly expressed and induced genes and is associated with intestinal transcription factor binding sites, including those for HNF4A and CDX2. TET1 induction occurred during differentiation, and TET1 knockdown altered gene expression and inhibited barrier formation of colonocytes. We find that the 5-hmC distribution in primary human colonocytes parallels the distribution found in differentiated cells in vitro, and that gene-specific 5-hmC changes in human colon cancers are directly correlated with changes in gene expression. Our results support a model in which 5-hmC regulates differentiation of adult human intestine and 5-hmC alterations contribute to the disrupted gene expression in colon cancer.
Differentiation of the intestinal epithelium along the spatially distinct crypt-luminal axis requires gene expression changes mediated in part by epigenetic pathways1. Modified cytosine bases include 5-methylcytosine (5-mC), 5-hydroxymethylcytosine (5-hmC), 5-formylcytosine (5-fC), and 5-carboxylcytosine (5-caC) and are one mechanism of epigenetic regulation of gene expression2,3,4,5. Changes to 5-mC have previously been implicated in regulation of intestinal differentiation. Specifically, differentially methylated regions have been identified between the stem and differentiated cells of the mouse small intestine. Co-localization of these regions with lineage-specific transcription factor binding sites and enhancer regions suggests that these 5-mC changes play an important role in intestinal differentiation6,7. Similar observations in hematopoietic and skin cell differentiation suggest that this is a common feature of lineage commitment8,9.
5-mC is converted to 5-hmC by the ten-eleven-translocation 5-methylcytosine dioxygenase (TET) enzymes2,3. Although accumulation of 5-mC at gene promoters and enhancers has been associated with gene repression10, 5-hmC has distinct biological functions. 5-hmC has been associated with active genes, has been proposed to play a role in numerous demethylation pathways, and may serve as a bona fide epigenetic mark with a distinct set of binding proteins11. Recently, mapping of 5-hmC in differentiation of erythroid cells and studies of Tet2 knockout mice have demonstrated that TET activity is critical for normal hematopoietic differentiation12,13. However, the importance of TET activity in the differentiation of other cell types remains unclear.
To determine if 5-hmC has a functional role in regulating colonocyte differentiation, we mapped 5-hmC changes during cell-cell adhesion-initiated differentiation of T84 colon adenocarcinoma cells since a similar system had previously been used to map chromatin regulatory regions of small intestinal differentiation14. When seeded at low density, these contact-naïve cells proliferate to form a confluent monolayer consisting of polarized cells with high transepithelial electrical resistance and morphological, structural, functional, and transcriptional features of colonocytes in vivo15,16. We extend our studies to human samples by comparing our in vitro 5-hmC maps to the 5-hmC profile of primary human colonocytes. Finally, since developmental pathways are frequently dysregulated in cancer17, we define regions losing and gaining 5-hmC in human colon cancers, and correlate these alterations with changes in gene expression.
Results
5-hmC is increased during T84 cell differentiation and is associated with epithelial pathways and transcription factor binding sites
T84 cells were seeded at low density and transepithelial electrical resistance was used to monitor T84 cell monolayer formation as cells differentiated (Fig. S1a). Total 5-hmC levels increased in differentiated cells (day 15) relative to proliferating cells (day 0) by dot blot assay (Fig. 1a). We used the hMe-Seal method to isolate and sequence 5-hmC-enriched DNA from cells at days 0, 4, 12, and 15 to determine how 5-hmC distribution changed during differentiation (Fig. S1b–e)18. Consistent with our dot blot results, we found that 5-hmC covered an increasing amount of the genome and that hMe-Seal peaks became more intense as differentiation progressed (Fig. 1b and Fig. S1f–j). We found no enrichment of 5-hmC at various genomic elements (CpG islands, CpG shores, promoters, 5′ UTRs, exons, introns, 3′ UTRs, and intergenic regions) at day 0, but significant enrichment at CpG shores and promoter regions by day 4. By day 12 and day 15, a strong 5-hmC signature was observed, with significant enrichment for 5-hmC over CpG islands, CpG shores, promoters, and gene bodies (Fig. 1c). We observed a relative preference for 5-hmC at CpG shores relative to CpG islands (Fig. S1k). To visualize 5-hmC changes over genes, we plotted the 5-hmC profile of an average gene at each time point. This demonstrated that 5-hmC was gained over promoters and gene bodies (Fig. 1d). KEGG pathway analysis demonstrated that 5-hmC peaks were enriched at genes involved in epithelial barrier function, including focal adhesion, adherens junctions, regulation of actin cytoskeleton, and endocytosis (Fig. 1e).
5-hmC frequently colocalizes with transcription factor binding sites12,19. Therefore, we analyzed genomic sequences covered by 5-hmC and found that they were predicted to bind the HNF4A, RXRA, and CDX2 transcription factors, which are known to regulate intestinal development (Table S1)20. We validated this result against ENCODE HNF4A ChIP-seq data acquired from HepG2 cells (Fig. S1l)21.
Previous work mapped early and late binding sites of HNF4A, CDX2, and GATA6 as well as the active enhancer mark H3K4me2 during differentiation of the Caco2 colon cancer cell line14. We calculated enrichment of 5-hmC at these regions and found that 5-hmC becomes especially enriched at the late binding sites of HNF4A and CDX2. We also observed enrichment of 5-hmC at enhancer regions (Figs 1f,g and S1m). GATA6 binding sites, served as a negative control and showed only weak overlap with 5-hmC. Furthermore, we examined HNF4A binding sites by Tet-assisted bisulfite sequencing (TAB-seq)22, which allows for quantification of cytosine, 5-mC, and 5-hmC at single base resolution. We performed TAB-seq at the HNF4A binding sites of RXRA, GNA12, VAV2 and at auto-regulatory sites of HNF4A (Figs 1h,i and S1n–t). In addition to gain of 5-hmC, we observed demethylation at the binding sites of VAV2 and GNA12. HNF4A binding sites have been associated with regions losing methylation in intestinal differentiation7. Since 5-hmC is implicated in demethylation pathways11, our data indicate that TET activity at these positions is a mechanism for this observation.
5-hmC is associated with highly expressed and induced genes
The increase of 5-hmC in differentiation, association of 5-hmC with important epithelial pathways, and colocalization of 5-hmC with lineage-specific intestinal transcription factor binding sites implicated 5-hmC gains as a regulator of colonocyte differentiation. To identify the relationship between 5-hmC and gene expression, we used RNA-seq to quantify transcript levels at each time point of differentiation. KEGG pathway analysis showed that genes associated with MAPK signaling were induced, and numerous metabolic and disease-associated pathways were repressed (Fig. S2a). We separated genes by expression level at day 15 (Fig. 2a), and found a positive correlation between 5-hmC and expression (Fig. 2b). To test the relationship between gene induction over the time course and 5-hmC levels, we identified genes with changed expression and plotted 5-hmC levels. We found that up-regulated genes had higher 5-hmC levels than genes that were down-regulated or had unchanged expression over the time course (Fig. 2c). To confirm these correlations, we performed the inverse analyses and examined expression as a function of 5-hmC. Highly hydroxymethylated genes were more highly expressed (Figs 2d and S2b), and genes with the highest levels of 5-hmC were more likely to be induced (Fig. 2C), and had a greater median fold change than genes with lower 5-hmC levels (Fig. S2c). Examples of genes that gain both 5-hmC and expression include MUC20, SLC11A2, SLC26A2, SLC6A8, FXYD1, and NDRG1 (Figs 2f,g and S2d–g). Association of 5-hmC with highly expressed and induced genes suggested that 5-hmC had an important role in regulation of gene expression during differentiation.
TET1 regulates gene expression in differentiating colonocytes
To test for functional significance of TET activity during differentiation, we performed TET1 knockdown experiments since TET1 was induced during differentiation (Fig. 3a, S3a,b) and since TET1 has been implicated as a tumor suppressor and regulator of the WNT pathway in carcinomas, including colon cancer23. Differential expression analysis of cells expressing a control versus TET1 targeting hairpin at day 15 (Fig. 3b) revealed that differentially regulated genes were enriched for genes coding for proteins targeted to the cell membrane and extracellular space (Fig. 3c). To identify potential direct TET1 targets, we identified differentially expressed genes that were hydroxymethylated, including: SLC26A3, SLC20A1, SLC26A2, AXIN2, ST3GAL4, RARG, and PHGR1 (Figs 3d–f, S3c,d). Given the exquisite sensitivity of SLC26A3 to TET1 knockdown, we integrated the 5-hmC profile of this gene with ENCODE data from HepG2 cells to infer the epigenetic context of 5-hmC at this locus. This analysis revealed that 5-hmC is gained at a region marked by H3K4me1, H3K4me2, and H3K27ac and is defined as an enhancer region by the ChromHMM model of chromatin segmentation (Fig. 3g)24. To determine if the gene expression changes observed on TET1 knockdown had a functional consequence on colonocyte function, transepithelial electrical resistance (TER) measurements were performed on cells depleted of TET1. We found that TET1 knock down attenuated TER at late time points (Fig. 3h), demonstrating that gene expression changes introduced on TET1 knockdown compromised barrier formation in vitro.
The 5-hmC distribution of primary human colonocytes is similar to that of differentiated T84 cells and is altered in human cancers in association with gene expression changes
To compare the 5-hmC profile observed in vitro with the distribution found in vivo, we performed hMe-Seal on isolated colonocyte fractions from two patients. We found that more regions of the genome were covered by 5-hmC peaks in primary human colonocytes than in differentiated T84 cells, consistent with a previous report that 5-hmC levels are lower in cell lines25. Despite the larger amount of 5-hmC found in vivo, sixty-five percent of regions covered by 5-hmC in vitro were also hydroxymethylated in vivo (Fig. 4a). Moreover, nearly all genes with a 5-hmC peak in vitro also had a 5-hmC peak in vivo (Fig. 4b). This analysis demonstrated that genes regulated by 5-hmC in vitro are also regulated by 5-hmC in vivo. Like our in vitro data, patient samples had high 5-hmC levels at HNF4A binding sites and at genes in similar KEGG pathways (Fig. S4a–e). These results show that while 5-hmC may regulate more genes in vivo the biological function is likely similar.
Altered regulation of developmental pathways such as the WNT, Hedgehog, Notch, and BMP pathways is recognized as a hallmark of cancer17. Since our results delineate a role for 5-hmC in regulating normal differentiation, we tested whether or not 5-hmC distribution is altered in primary human colon cancers by performing cell-type specific hMe-Seal on isolated colonocytes from matched colon cancer and adjacent normal samples to identify genes that lose or gain 5-hmC in cancer (Fig. S4f). Using The Cancer Genome Atlas (TCGA) dataset, we found that 163 of 222 genes losing 5-hmC in cancer also lost expression, and that 96 of 175 genes gaining 5-hmC in cancer gained expression (Fig. S4g). Examples of genes with hydroxymethylation and gene expression changes include: CA2, FMN2, PDCD4, PKIB, SLC26A2, BMP7, NKD2, TESC, and TGFBI (Figs 4c–h and S4h–l). Gene expression changes for these genes were further validated using qPCR (Fig. S4m). We further quantified changes in 5-hmC, 5-mC, and cytosine in additional independent patient samples at PKIB and FMN2 using TAB-seq. Consistent loss of 5-hmC at these loci demonstrates that this is a common feature of colon cancers (Fig. 4i,j).
Given that we observed gene-specific changes in 5-hmC, we performed pathway analysis to determine if certain KEGG pathways were enriched for gaining or losing 5-hmC in cancer. Surprisingly, we could not identify pathways that were significantly enriched for genes gaining or losing 5-hmC in the colon cancer samples relative to adjacent normals. Since we had identified the focal adhesion, adherens junctions, regulation of actin cytoskeleton, and endocytosis pathways as being enriched for 5-hmC in differentiated T84 cells, we determined if these pathways had alterations to 5-hmC levels in cancer. Plotting 5-hmC levels over genes in these pathways revealed that there were no consistent changes of 5-hmC levels over these pathways (Fig. S5), demonstrating that while gene-specific changes may contribute to altered expression in cancer, consistent pathway-wide changes in 5-hmC could not be identified.
Discussion
Our results identify a previously unappreciated role for 5-hmC in regulating epithelial differentiation. We observe accumulation of 5-hmC at intestine specific genes and transcription factor binding sites. Integration of hMe-Seal and RNA-seq data demonstrates that 5-hmC is particularly enriched at highly expressed and/or induced genes, and TET1 knockdown shows functional significance of TET activity in colonocyte differentiation. Mapping of 5-hmC in primary colonocytes and cancers demonstrates that perturbations to TET activity are a mechanism of altered expression in cancer.
Previous studies have profiled cytosine modifications during differentiation of small intestinal epithelium using bisulfite-based techniques unable to distinguish 5-mC and 5-hmC6,7,26,27. These experiments revealed relatively few changes between the stem and differentiated populations. In contrast, our studies reveal that gains in 5-hmC during differentiation are relatively common. This suggests that transitions between 5-mC and 5-hmC, which would not be detected by bisulfite-based techniques, may be more frequent than transitions between modified (5-mC and 5-hmC) and unmodified cytosines.
Co-localization of 5-hmC with lineage-specific transcription factor binding sites suggests either that multiple transcription factors direct TET enzymes to target loci or that TET activity is required for subsequent transcription factor binding. Previous work has demonstrated that CDX2 functions as a pioneer transcription factor triggering chromatin remodeling28. Although we find greater enrichment of 5-hmC at HNF4A binding sites, these transcription factors frequently regulate the same regions14. Future work should explore the molecular interactions between these transcription factors and the TET enzymes, and determine how TET activity at these regions is coordinated with other aspects of chromatin remodeling. Our data at SLC26A3 and regions marked by H3K4me2 suggest that TET activity is likely integrated with other chromatin modifying enzymes at enhancer regions. Future work should interrogate the role of each TET enzyme in regulating distinct chromosomal regions and gene sets during epithelial differentiation.
Previous reports have identified 5-hmC loss as a hallmark of cancer, including colon cancer23,29,30,31,32,33,34. However, many of these studies measured only global 5-hmC content without genome-wide mapping. Our studies reveal both locus-specific losses and gains of 5-hmC in cancer. These 5-hmC changes occur at genes coding for proteins involved in regulating intestinal development and cancer cell proliferation. For example, we demonstrate significant loss of 5-hmC at FMN2, a gene whose protein product promotes cell cycle arrest by inhibiting the degradation of p2135. 5-hmC gains occur at genes including NKD2, a gene coding for a regulator of Wnt and TGFa/EGFR signaling36,37,38, and BMP7, a member of the TGF-beta superfamily39. Gain of 5-hmC also occurs at TGFBI, a gene whose over-expression in colon cancer has been considered as a prognostic biomarker and target for therapy since TGFBI promotes metastasis40. Similarly, TESC codes for a member of the calcineurin homologous protein family that has been reported to activate NF-κB and possibly Akt signaling pathways to control the survival and proliferation of colorectal cancer cells41.
Correlation of 5-hmC alterations with tumor-associated changes in gene expression suggests that 5-hmC changes can contribute to both gene activation and repression in cancer. These epigenetic lesions are of particular interest since they are potentially reversible with epigenetic targeting agents, and future work should determine whether or not any of these locus-specific 5-hmC changes predict response to these drugs. Finally, by defining regions with altered 5-hmC in colon cancer, our work provides a starting point for identifying key genomic regions that could be epigenetically altered in other epithelial malignancies and advances our understanding of the molecular mechanisms underlying tumorigenesis.
Methods
T84 in vitro differentiation cell culture
T84 human colon carcinoma epithelial cells were grown in biologic triplicate on tissue culture-treated 6-well Transwell filter support units (0.4-μm pore size). The culture medium employed was a 1:1 mixture of Dulbecco’s modified Eagle’s medium (DMEM) and Ham’s F-12 nutrient mixture supplemented with heat-inactivated fetal bovine serum (FBS) and 100 U/mL penicillin/0.1 mg/mL streptomycin. The cells were grown in a humidified 5% CO2 atmosphere at 37°C. Culture media was replaced approximately every 48 hours. The monolayers were grown until achieving differentiation with transepithelial resistance (TER) measurements of >2,000 Ω ⋅ cm2.
Transepithelial electrical resistance (TER) measurement
TER was measured with a MilliCell-ERS volt-ohmmeter (Millipore) and chopstick current electrodes approximately every 48 hours prior to culture medium exchange. All TER values were normalized by subtracting the baseline resistance measurement across the Transwell filter support unit in the absence of cells.
Quantification of 5-hmC by dot blot
Dot blot assays were performed as previously described19. Briefly, genomic DNA was sonicated, 5-hmC was conjugated to UDP-6-N3-Glu18, and biotinylated by addition of 4.5 nm Click-IT Biotin DIBO Alkyne (Life Technologies). The biotin labeled DNA was purified and spotted on HyBond N + membrane (GE Healthcare Life Sciences). Membranes were probed with Avidin-HRP and dots were visualized with Western Lightning ECL reagent (PerkinElmer).
5-hydroxymethylcytosine selective chemical labeling (hMe-Seal)
hMe-Seal followed previously published protocols18,19. Briefly, 20 μg genomic DNA was sonicated to 250 base pairs (Covaris), conjugated to UDP-6-N3-Glu with β-GT, and purified with Micro Bio-Spin Columns (Bio-Rad). 5-hmC was subsequently biotinylated by addition of DMCO-S-S-PEG3-Biotin Conjugate (Click Chemistry Tools). Biotinylated DNA was isolated using streptavidin-coated magnetic beads (Invitrogen). After washing beads, labeled DNA was released with 50 mM DTT and purified sequentially with Micro Bio-Spin Columns (Bio-Rad) and MinElute reaction cleanup kits (Qiagen). The concentration of the isolated DNA was measured by Qubit (Qiagen).
hMe-Seal read mapping and peak calling
Reads were mapped to the hg19 genome using the BWA- backtrack algorithm. Peak calling was performed using MACS 1.4.
Enrichment calculations for 5-hmC at genomic elements
Enrichment for 5-hmC at genomic regions used previously described permutation analyses19.
Genomic coordinates of 5′ UTRs, 3′ UTRs, exons, introns, and CpG islands were obtained from the UCSC Table Browser42. Promoter regions were defined by downloading TSS sites from the UCSC Table Browser and using the bedtools suite to generate a bed file of coordinates 2 kb downstream to 500 bp upstream of TSSs. Genomic coordinates of CpG shores were deduced from the locations of CpG islands by using the bedtools suite to define regions 2 kb adjacent to each CpG island.
To calculate enrichment of 5-hmC over each annotation, each peak identified by MACS 1.4 was annotated with the number of tags mapping to that region (also outputted by MACS 1.4). The number of tags mapping to each peak was divided by the length of that peak to determine tags/bp over that peak. The peak bed file was then intersected against genomic annotations. For each intersection of a peak and a genomic annotation, the length of the intersection in bases was multiplied with the number of tags per base for that peak island in order to obtain the number of sequencing tags at each intersection:
(Length of the intersection between peak and genomic annotation in bps) × (sequencing tags per base at peak island) = sequencing tags over intersection
Summing these values for all intersections genome-wide gives the absolute number of sequencing tags over each annotation being analyzed. This value served as an “observed” value for enrichment calculations.
To calculate enrichment, the observed value was divided by an expected value. The expected value was determined by scrambling 5-hmC peaks across the genome using the bedtools shuffle command and repeating the calculation described above. Shuffling was performed 1,000 times for each enrichment analysis, and the average value obtained from these random permutations was used as the expected value in enrichment calculations.
For all enrichment calculations, each biological replicate was treated independently, and results averaged to give enrichment values at each time point and over each annotation.
KEGG pathway analysis
A list of genes with 5-hmC peaks within 2 kb of their gene body was generated by first downloading the positions of all RefSeq genes from the UCSC table browser. The bedtools slop command was used to lengthen the positions of these genes by 2 kb. This file was then intersected with the output of MACS 1.4 peaks at day 15. The gene list outputted by this command was used in KEGG pathway analysis performed on the DAVID Bioinformatics Resources 6.7 website (https://david.ncifcrf.gov/).
Genomic positions of GATA6, HNF4A, CDX2, and H3K4me2 in proliferating and differentiated cells
Positions of these transcription factor binding sites and histone modification have previously been determined in Caco2 cells14. The positions of these features were acquired from GEO (GSE23436). Peaks were converted from hg18 to hg19 coordinates using the UCSC liftover tool.
Quantification of 5-hmC over genes
To quantify 5-hmC over gene bodies, a gtf file was created with positions of all RefSeq genes. For genes with a MACS 1.4 called peak over their gene body, Ht-seq count was used to sum hMe-Seal reads over these genes. To control for gene length, the number of 5-hmC sequencing tags mapping to each gene was then divided by the length of the gene to give 5-hmC density over each gene. Plotting the distribution of gene density yielded the histogram in Fig. S2b. The histogram demonstrated a normal distribution with a right tail. Therefore, genes were binned by 0-25, 25-75, 75-98, and 98–100 percentile.
Plotting hMe-Seal data around genomic elements
To visualize 5-hMe-Seal sequencing data around a given genomic feature (e.g., gene bodies, CpG islands, and transcription factor binding sites), the ngs-plot software was used43. Positions of CpG islands were obtained from the UCSC table browser. Positions of transcription factor binding sites were acquired from previously published datasets as described above. Sam files aligned by bwa were used as input.
Tet-assisted bisulfite sequencing (TAB-seq)
TAB-Seq was performed using the Wisegene 5-hmC TAB-Seq kit according to the manufacturer’s protocol. Bisulfite treatment was performed using the MethylCode Bisulfite Conversion Kit (Life Technologies, MECOV-50). PCR reactions used ZymoTaq (Zymo Research). PCR products were cloned using the NEB PCR cloning kit (NEB E1202S).
RNA-seq read mapping and expression analysis
RNA-seq analysis followed a published protocol44. Reads were mapped to the hg19 genome using Bowtie2 and TopHat2. Gene expression quantification was performed using Cufflinks to calculate fragments per kilobase of transcript per million mapped reads (FPKM), and Cuffdiff to perform differential expression analysis. All gene analyses used a GTF file obtained from Ensembl.
Generation of T84 cells stably expressing shRNAs
Production of lentiviral particles and subsequent infection of target cells followed the pLKO.1 protocol found at http://www.addgene.org/tools/protocols/plko/. Infection was performed with a final polybrene concentration of 8 μg/mL. Cells expressing the desired shRNAs were selected using puromycin (8 μg/mL, Santa Cruz). pLKO.1 plasmids coding for TET1 targeting hairpins were previously described45.
Analysis of the cancer genome atlas (TCGA) data
Level 3 TCGA data were downloaded from the TCGA data matrix. Gene level RSEM data, calculated by the TCGA, were used to measure expression. Heat maps were generated by the UCSC cancer browser46. Wilcoxon-Rank sum tests were used to test for differential expression on TCGA data. Test statistics were Bonferonni corrected.
Preparation of isolated colonocytes from surgically resected mucosa
We adapted a method to isolate a nearly pure population of colonocytes47. Colon crypts from surgically resected mucosa stripped from muscle of four patients with moderately differentiated sporadic colon cancer and adjacent normal colon were isolated by EDTA/EGTA at 4 °C. Crypts were released by gentle shaking and separated from stromal cells by centrifugation at 40 g for 2 min. Isolated crypts stained for positive (CK20) and negative (vimentin) markers.
Additional Information
Accession Numbers: Next generation sequencing datasets are available at http://www.ncbi.nlm.nih.gov/geo: GSE55391.
How to cite this article: Chapman, C. G. et al. TET-catalyzed 5-hydroxymethylcytosine regulates gene expression in differentiating colonocytes and colon cancer. Sci. Rep. 5, 17568; doi: 10.1038/srep17568 (2015).
Supplementary Material
Acknowledgments
We would like to thank Julie-Aurore Losman (Dana Farber Cancer Institute) for plasmids coding for control and TET1-targeting shRNAs. This project was supported by the Gastrointestinal Research Foundation Associates’ Board and the National Center for Advancing Translational Sciences of the National Institutes of Health through Grant Number UL1 TR000430.
The authors declare no competing financial interests.
Author Contributions C.G.C. and C.J.M. performed the experiments, analyzed the data, and wrote the manuscript. F.W., A.C. and J.M. performed experiments. K.M., F.B. and M.B.B. provided colonocyte fractions from primary patient samples. J.H.K. and L.A.G. conceived the study and edited the manuscript.
04/28/2016
A correction has been published and is appended to both the HTML and PDF versions of this paper. The error has been fixed in the paper.
References
- Hammoud S. S., Cairns B. R. & Jones D. A. Epigenetic regulation of colon cancer and intestinal stem cells. Current opinion in cell biology 10.1016/j.ceb.2013.01.007 (2013). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Tahiliani M. et al. Conversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1. Science 324, 930–935 (2009). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kriaucionis S. & Heintz N. The nuclear DNA base 5-hydroxymethylcytosine is present in Purkinje neurons and the brain. Science 324, 929–930 (2009). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Ito S. et al. Tet Proteins Can Convert 5-Methylcytosine to 5-Formylcytosine and 5-Carboxylcytosine. Science 333, 1300–1303 (2011). [DOI] [PMC free article] [PubMed] [Google Scholar]
- He Y. F. et al. Tet-Mediated Formation of 5-Carboxylcytosine and Its Excision by TDG in Mammalian DNA. Science 333, 1303–1307 (2011). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kaaij L. T. J. et al. DNA methylation dynamics during intestinal stem cell differentiation reveals enhancers driving gene expression in the villus. Genome Biol 14, R50 (2013). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Sheaffer K. L. et al. DNA methylation is required for the control of stem cell differentiation in the small intestine. Genes Dev. 28, 652–664 (2014). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Bock C. et al. DNA methylation dynamics during in vivo differentiation of blood and skin stem cells. Mol Cell 47, 633–647 (2012). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Shearstone J. R. et al. Global DNA demethylation during mouse erythropoiesis in vivo. Science 334, 799–802 (2011). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Moen E. L. et al. New themes in the biological functions of 5-methylcytosine and 5-hydroxymethylcytosine. Immunological reviews 263, 36–49 (2015). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Mariani C. J., Madzo J., Moen E. L., Yesilkanal A. & Godley L. A. Alterations of 5-hydroxymethylcytosine in human cancers. Cancers (Basel) 5, 786–814 (2013). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Madzo J. et al. Hydroxymethylation at Gene Regulatory Regions Directs Stem/Early Progenitor Cell Commitment during Erythropoiesis. Cell Rep 6, 231–244 (2014). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Madzo J., Vasanthakumar A. & Godley L. A. Perturbations of 5-Hydroxymethylcytosine Patterning in Hematologic Malignancies. Seminars in hematology 50, 61–69 (2013). [DOI] [PubMed] [Google Scholar]
- Verzi M. P. et al. Differentiation-Specific Histone ModificationsReveal Dynamic Chromatin Interactions and Partners for the Intestinal Transcription Factor CDX2. Dev Cell 19, 713–726 (2010). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Christensen J. et al. Defining new criteria for selection of cell-based intestinal models using publicly available databases. BMC Genomics 13, 274 (2012). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Madara J. L., Stafford J., Dharmsathaphorn K. & Carlson S. Structural analysis of a human intestinal epithelial cell line. Gastroenterology 92, 1133–1145 (1987). [DOI] [PubMed] [Google Scholar]
- Bertrand F. E., Angus C. W., Partis W. J. & Sigounas G. Developmental pathways in colon cancer: crosstalk between WNT, BMP, Hedgehog and Notch. Cell Cycle 11, 4344–4351 (2012). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Song C.-X. et al. Selective chemical labeling reveals the genome-wide distribution of 5-hydroxymethylcytosine. Nat Biotechnol 29, 68–72 (2010). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Mariani C. J. et al. TET1-Mediated Hydroxymethylation Facilitates Hypoxic Gene Induction in Neuroblastoma. Cell Rep 7, 1343–1352 (2014). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Heinz S. et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell 38, 576–589 (2010). [DOI] [PMC free article] [PubMed] [Google Scholar]
- ENCODE Project Consortium et al. An integrated encyclopedia of DNA elements in the human genome. Nature 489, 57–74 (2012). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Yu M. et al. Base-resolution analysis of 5-hydroxymethylcytosine in the mammalian genome. Cell 149, 1368–1380 (2012). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Neri F. et al. TET1 is a tumour suppressor that inhibits colon cancer growth by derepressing inhibitors of the WNT pathway. Oncogene 1–9 10.1038/onc.2014.356 (2014). [DOI] [PubMed] [Google Scholar]
- Ernst J. & Kellis M. ChromHMM: automating chromatin-state discovery and characterization. Nat Meth 9, 215–216 (2012). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Nestor C. E. et al. Rapid reprogramming of epigenetic and transcriptional profiles in mammalian culture systems. Genome Biol 16, 11 (2015). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Jin S.-G., Kadam S. & Pfeifer G. P. Examination of the specificity of DNA methylation profiling techniques towards 5-methylcytosine and 5-hydroxymethylcytosine. Nucleic Acids Res 38, e125 (2010). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Huang Y. et al. The Behaviour of 5-Hydroxymethylcytosine in Bisulfite Sequencing. PLoS ONE 5, e8888 (2010). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Verzi M. P., Shin H., San Roman A. K., Liu X. S. & Shivdasani R. A. Intestinal Master Transcription Factor CDX2 Controls Chromatin Access for Partner Transcription Factor Binding. Mol. Cell. Biol. 33, 281–292 (2012). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Haffner M. C. et al. Global 5-hydroxymethylcytosine content is significantly reduced in tissue stem/progenitor cell compartments and in human cancers. Oncotarget 2, 627–637 (2011). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Jin S. G. et al. 5-Hydroxymethylcytosine Is Strongly Depleted in Human Cancers but Its Levels Do Not Correlate with IDH1 Mutations. Cancer Res 71, 7360–7365 (2011). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Yang H. et al. Tumor development is associated with decrease of TET gene expression and 5-methylcytosine hydroxylation. Oncogene 32, 663–669 (2013). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Chen M. L. et al. Quantification of 5-Methylcytosine and 5-Hydroxymethylcytosine in Genomic DNA from Hepatocellular Carcinoma Tissues by Capillary Hydrophilic-Interaction Liquid Chromatography/Quadrupole TOF Mass Spectrometry. Clinical Chemistry 59, 824–832 (2013). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kudo Y. et al. Loss of 5-hydroxymethylcytosine is accompanied with malignant cellular transformation. Cancer Science 103, 670–676 (2012). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Lian C. G. et al. Loss of 5-Hydroxymethylcytosine Is an Epigenetic Hallmark of Melanoma. Cell 150, 1135–1146 (2012). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Yamada K., Ono M., Perkins N. D., Rocha S. & Lamond A. I. Identification and Functional Characterization of FMN2, a Regulatorof the Cyclin-Dependent Kinase Inhibitor p21. Mol Cell 49, 922–933 (2013). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Zeng W. et al. naked cuticle encodes an inducible antagonist of Wnt signalling. Nature 403, 789–795 (2000). [DOI] [PubMed] [Google Scholar]
- Hu T. & Li C. Convergence between Wnt-β-catenin and EGFR signaling in cancer. Mol Cancer 9, 236 (2010). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Li C. et al. Myristoylated Naked2 escorts transforming growth factor alpha to the basolateral plasma membrane of polarized epithelial cells. Proc. Natl. Acad. Sci. USA 101, 5571–5576 (2004). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Motoyama K. et al. Clinical Significance of BMP7 in Human Colorectal Cancer. Ann Surg Oncol 15, 1530–1537 (2008). [DOI] [PubMed] [Google Scholar]
- Ma C. et al. Extracellular matrix protein ig-h3/TGFBI promotes metastasis of colon cancer by enhancing cell extravasation. Genes Dev. 22, 308–321 (2008). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kang Y. H. et al. The EF-hand calcium-binding protein tescalcin is a potential oncotarget in colorectal cancer. Oncotarget 5, 2149–2160 (2014). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Karolchik D. et al. The UCSC Table Browser data retrieval tool. Nucleic Acids Res 32, D493–6 (2004). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Shen L., Shao N., Liu X. & Nestler E. ngs.plot: Quick mining and visualization of next-generation sequencing data by integrating genomic databases. BMC Genomics 15, 1–14 (2014). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Roberts A. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc 7, 562–578 (2012). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Losman J.-A. et al. (R)-2-Hydroxyglutarate Is Sufficient to Promote Leukemogenesis and Its Effects Are Reversible. Science 10.1126/science.1231677 (2013). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Zhu J. et al. The UCSC Cancer Genomics Browser. Nat Meth 6, 239–240 (2009). [DOI] [PMC free article] [PubMed] [Google Scholar]
- Seidelin J. B., Horn T. & Nielsen O. H. Simple and efficient method for isolation and cultivation of endoscopically obtained human colonocytes. Am. J. Physiol. Gastrointest. Liver Physiol. 285, G1122–G1128 (2003). [DOI] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.