TET1 regulates hypoxia-induced epithelial-mesenchymal transition by acting as a co-activator
© Tsai et al.; licensee BioMed Central Ltd. 2014
Received: 18 April 2014
Accepted: 28 October 2014
Published: 3 December 2014
Hypoxia induces the epithelial-mesenchymal transition, EMT, to promote cancer metastasis. In addition to transcriptional regulation mediated by hypoxia-inducible factors, HIFs, other epigenetic mechanisms of gene regulation, such as histone modifications and DNA methylation, are utilized under hypoxia. However, whether DNA demethylation mediated by TET1, a DNA dioxygenase converting 5-methylcytosine, 5mC, into 5-hydroxymethylcytosine, 5hmC, plays a role in hypoxia-induced EMT is largely unknown.
We show that TET1 regulates hypoxia-responsive gene expression. Hypoxia/HIF-2α regulates the expression of TET1. Knockdown of TET1 mitigates hypoxia-induced EMT. RNA sequencing and 5hmC sequencing identified the set of TET1-regulated genes. Cholesterol metabolic process genes are among the genes that showed high prevalence and statistical significance. We characterize one of the genes, INSIG1 (insulin induced gene 1), to confirm its expression and the 5hmC levels in its promoter. Knockdown of INSIG1 also mitigates hypoxia-induced EMT. Finally, TET1 is shown to be a transcriptional co-activator that interacts with HIF-1α and HIF-2α to enhance their transactivation activity independent of its enzymatic activity. TET1 acts as a co-activator to further enhance the expression of INSIG1 together with HIF-2α. We define the domain in HIF-1α that interacts with TET1 and map the domain in TET1 that confers transactivation to a 200 amino acid region that contains a CXXC domain. The TET1 catalytically inactive mutant is capable of rescuing hypoxia-induced EMT in TET1 knockdown cells.
These findings demonstrate that TET1 serves as a transcription co-activator to regulate hypoxia-responsive gene expression and EMT, in addition to its role in demethylating 5mC.
Cells develop various mechanisms to cope with hypoxia and survive ,. Hypoxia induces the epithelial-mesenchymal transition (EMT) to promote cancer metastasis -. In addition to transcriptional regulation mediated by hypoxia-inducible factors (HIFs), other epigenetic mechanisms of gene regulation, such as histone modifications and DNA methylation, are utilized under hypoxia -. Certain chromatin changes have been observed during EMT. For example in Snail-induced EMT, loss of H3K4Me3, H3K4Ac, and H3K27Ac, and gain of H3K27Me3 were observed for genes repressed, whereas gain of H3K4Me3, H3K4Me1, and loss of H3K27Me3 were observed for genes activated . Other specific chromatin changes have also been observed during hypoxia or TGF-β-induced EMT ,.
DNA demethylation is recently shown to be an important epigenetic mechanism that regulates gene expression due to the discovery of TET (ten-eleven translocation) enzymes that demethylate DNA ,. TET enzymes have been shown to convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) to regulate gene expression -. Tet proteins can also convert 5-methylcytosine to 5-formylcytosine and 5-carboxycytosine ,. Tet1 and Tet2 are regulated by Oct4 during somatic cell reprogramming into induced pluripotent stem cells . However, other mechanisms regulating expression of Tet1 genes remain largely unknown.
In this report, we explored the mechanism of TET1 regulation by hypoxia. The role TET1 in regulating the process of EMT induced by hypoxia was investigated. The set of genes that was regulated by TET1 under hypoxia and their role in hypoxia-induced EMT were delineated. Finally, we showed the extra role of TET1 in serving as a transcription co-activator. These results provide a fresh insight into regulation of hypoxia-responsive gene expression by TET1 and further expand the role of TET1.
Results and discussion
Hypoxia/HIF-2α activates TET1 expression and knockdown of TET1 mitigates hypoxia-induced epithelial-mesenchymal transition
Since EMT is induced by hypoxia to promote tumor metastasis -, we tested whether TET1 plays a significant role in hypoxia-induced EMT. TET1 knockdown was performed in a H1299 lung cancer cell line. Exposure of H1299 cells to hypoxia induced EMT as characterized by repression of E-cadherin/plakoglobin and upregulation of vimentin/N-cadherin expression , (Figure 1e). Knockdown of TET1 abolished the upregulation of mesenchymal gene expression including vimentin and N-cadherin (Figure 1e). Immunofluorescence assays showed the no change of E-cadherin expression and abolishment of vimentin upregulation in H1299 cells under TET1 knockdown (Additional file 3a). TET1 knockdown also significantly decreased the in vitro migration and invasion activity induced by hypoxia in H1299 cells (Additional file 3b). The same results were obtained using the head and neck cancer FADU cell line (Additional file 4). These results indicate that TET1 plays a crucial role in hypoxia-induced EMT and induction of in vitro migration/invasion activity.
Cholesterol metabolic process is regulated by TET1 and knockdown of INSIG1 mitigates hypoxia-induced epithelial-mesenchymal transition
Demethylation of the INSIG1promoter regions by TET1
Regulation of glucose metabolism gene expression by TET1 or INSIG1
TET1 acts as a transcription co-activator to regulate hypoxia-induced gene expression
For TET1 to be a co-activator of HIFs, there should be interaction between TET1 and HIFs. Co-immunoprecipitation assays were performed by co-transfecting HIF-1α and TET1 into 293 T cells. The results showed that the anti-HIF-1α antibody pulled down both HIF-1α and TET1 (Figure 5b). Reverse assays also showed that the anti-TET1 antibody pulled down TET1 and HIF-1α (Additional file 13a). HIF-1α also interacted with the enzymatically inactive TET1 (TET1-CDmt) by co-immunoprecipitation assays (Figure 5b and Additional file 13a). The interaction also occurred between HIF-2α and TET1 as well as between HIF-2α and TET1-CDmt (Additional file 13b and c). To further map the domain in HIF-1α that interacts with TET1, co-immunoprecipitation assays using various truncation mutants of HIF-1α and TET1 showed that the amino terminal domain (a.a. 1-400) of HIF-1α interacted with TET1 (Additional file 14). Further mapping experiments showed that the a.a. 175-305 region of HIF-1α interacted with TET1 (Figure 5c and Additional file 15a and b). To map the domain in TET1 that has transactivation activity, a yeast one hybrid assay was performed. The a.a. 451-674 region in TET1 displayed transactivation activity in contrast to the other regions that did not provide transactivation activity (Figure 5d). Further dissection of the transactivation domain showed that dividing this domain into two regions disrupted the transactivation activity (Additional file 15c). All the above results indicate that TET1 serves as a co-activator to interact with HIFs and activate the expression of hypoxia downstream targets.
Since O-linked N-glucosamine (O-GlcNAc) transferase (OGT) interacts with TET1 to modulate TET1 levels and regulate CpG island methylation ,, we tested whether OGT could also interact with HIF-1α and TET1 to regulate HIF-1α target gene expression. Co-immunoprecipitation experiments using an anti-TET1 antibody to pull down extracts from H1299 cells showed that both HIF-1α and OGT can be pulled down by an anti-TET1 antibody (Figure 5e), indicating that endogenous TET1 interacts with HIF-1α and OGT. However, only TET1 and CBP can be pulled down by an anti-HIF-1α antibody, but not OGT (Additional file 16a). This result suggested that there may not be direct interaction of HIF-1α with OGT. Other co-immunoprecipitation experiments showed that only TET1 can be pulled down by an anti-OGT antibody, but not HIF-1α or CBP (Additional file 16b). In addition, only HIF-1α can be pulled down by an anti-CBP antibody, but not OGT or TET1 (Additional file 16c). Co-expression of OGT with HIF-1α and CBP did not further increase the HIF-1α target gene expression using transient transfection assays (data not shown). These results indicate that TET1 may form different complexes with HIF-1α/CBP or OGT to regulate different target gene expression.
Rescue of hypoxia-induced EMT in TET1 knockdown cells by a TET1 catalytically inactive mutant
This report demonstrated that hypoxia-induced TET1 uses another important mechanism by serving as a transcription co-activator to regulate hypoxia-responsive gene expression and EMT. It is interesting that INSIG1, the master regulator of cholesterol biosynthesis ,, is involved in the regulation of hypoxia-induced EMT, highlighting the novelty and importance of lipid metabolism in regulating EMT and metastasis. Due to the ability of INSIG1 to inhibit cholesterol biosynthesis , that may shun down lipid synthesis to favor glucose utilization similar to the role of AMPK that also inhibits cholesterol synthesis , the activation of INSIG1 expression through hypoxia-induced TET1 may also contributes to Warburg effect observed under hypoxia as supported by our results (Figure 4) ,. It is interesting that INSIG1 did not regulate hypoxia-induced EMT through Twist1 or Snail, indicating that other EMT regulators or mechanisms (for example, chromatin modifiers) may be involved. Since we readily observed the regulation of various glucose metabolic genes by INSIG1, the cross talk between metabolism and other mechanisms may explain the role of INISG1 in hypoxia-induced EMT.
Our results show that TET1 serves as a 5mC-specific dioxygenase to mediate DNA demethylation of promoters ,, while it also simultaneously interacts with transcription factors as a co-activator to promote gene transcription. This mechanism is a unique model of gene transcription and a model is presented (Figure 6c). TET1 therefore represents a new class of co-activator compared to the recent demonstration of PKM2 (pyruvate kinase M2) as a HIF-1α co-activator . The domain in TET1 that confers transactivation activity is mapped (a.a. 451-674) by co-immunoprecipitation and yeast one-hybrid assays. It is interesting that this region contains the CXXC domain that is located in the TET enzymes . In contrast to the recruitment of TET1 by NANOG that requires the catalytic activity of TET1 to establish pluripotency , synergistic activation by TET1 and HIFs is independent of its enzymatic activity and supports that TET1 is a bona fide transcription co-activator. Our results demonstrate that TET1 plays a crucial role in the regulation of hypoxia-responsive gene expression by acting as a transcriptional co-activator since the catalytically inactive TET1 could rescue hypoxia-induced EMT in TET1 knockdown cells. It is possible that TET1 also serves as a co-repressor since Tet1 has dual functions on gene transcription and also interacts with SIN3A ,. Although TET1 also interacts with OGT to regulate gene expression ,, our results showed that TET1/OGT may form a different complex from the HIF-1α/CBP/TET1 complex. It will be important to search for other transcription factors that can interact with TET1 to regulate gene expression. Our results further expand the role of TET1 beyond epigenetic regulation of gene expression.
Materials and methods
The cell lines used were described , including the human hypopharyngeal squamous carcinoma cell line FADU, human tongue squamous carcinoma cell line SAS, and lung cancer cell lines H1299 and A549 (obtained from ATCC). Human embryonic kidney 293 T cell line was used in transient transfection experiments. SNU398 is a human hepatoma cell line obtained from Dr. YS Jou (Academia Sinica, Taiwan).
Protein extraction, western blot analysis, RNA extraction, quantitative real-time PCR, in vitromigration/invasion assay
For extraction of proteins from cell lines, cell lysis buffer (50 mM Tris, pH 7.5, 150 mM NaCl, 0.5% sodium deoxycholate, 1% NP-40, 0.1% SDS) containing protease inhibitors was used. Cell lysates were clarified by centrifugation at 13,000 rpm, 4°C for 10 min. The protein content was determined by Bradford method (Bio-Rad Laboratories, Hercules, CA, USA). For western blot analysis, 50 μg to 100 μg protein extracts from each clone were loaded to 10% SDS-PAGE gels and transferred to nitrocellulose filters. The filters were probed with different antibodies (Additional file 17), and an anti-ß-actin antibody was selected as a loading control. Signals were developed using an ECL chemiluminescence kit (Millipore Corporation, Billerica, MA, USA). Total RNA was isolated using Trizol reagent (Life Technologies Corporation, Carlsbad, CA, USA) according to manufacturer’s recommendations. Single-stranded cDNA was synthesized by the RevertAid™ First Strand cDNA Synthesis Kit (Fermentas International, Inc., Burlington, Canada). Real-time PCR was performed on a StepOnePlus™ Real-Time PCR System (ABI: Applied Biosystems, Foster City, CA USA) according to the manufacturer’s instructions. The 2-ΔΔCt method of relative quantification was used to estimate the copy number of gene expression, and 18S was selected as an internal control. The oligonucleotides used for real-time PCR were shown in Additional file 18. For in vitro migration/invasion assay, eight-μm pore size Boyden chamber was used. Cells (4 × 104) in 0.5% serum-containing RPMI were plated in the upper chamber and 15% fetal bovine serum was added to RPMI 1640 in the lower chamber as a chemo-attractant. For invasion assay, the upper side of the filter was covered with Matrigel (BD Biosciences, San Jose, CA, USA) (1:2 dilution with RPMI). After 12 h for migration assay or 24 h for invasion assay, cells on the upper side of the filter were removed, and cells that remained adherent to the underside of membrane were fixed in 4% formaldehyde and stained with Hoechst 33342 dye. The number of migrated cells was counted using a fluorescence microscope. Ten contiguous fields of each sample were examined using a 40× objective to obtain a representative number of cells which migrated/invaded across the membrane.
Plasmids, transfection, and luciferase assays
The plasmids used and transfection methods were described using calcium phosphate or lipofection transfection method ,. The expression vectors of HIF-1α, HIF-2α, CBP, Twist1 promoter-driven luciferase construct, and WDR5 promoter-driven luciferase construct were described ,. The TET1 and INSIG1 promoter regions were cloned and the reporter constructs were shown in Additional file 2. The oligonucleotides used for plasmis constructions were shown in Additional file 19. The reporter constructs were co-transfected into 293 T cells with different expression vectors and an internal control plasmid. Luciferase assays were performed using the same amount of cell extracts and corrected for transfection efficiency using an internal control (GFP) ,.
Quantitative chromatin immunoprecipitation (Chip)
ChIP assay was performed as described . Briefly, cells were cross-linked with 1% formaldehyde for 10 min and stopped by adding glycine to a final concentration of 0.125 M. Fixed cells were washed twice with TBS (20 mM Tris, pH 7.5, 150 mM NaCl) and harvested in 5 mL of SDS buffer (50 mM Tris, pH 8.0, 0.5% SDS, 100 mM NaCL, 5 mM EDTA, and protease inhibitors). Cells were pelleted by centrifugation and suspended in 2 mL of IP buffer (100 mM Tris, pH 8.6, 0.3% SDS, 1.7% Triton X-100, 5 mM EDTA). Cells were sonicated with a 0.25-inch diameter probe for 15 s twice using an MSE-soniprep 1500 sonicator (setting 18). For each immunoprecipitation, 1 mL of lysate was precleared by adding 50 μL of blocked protein A beads (50% protein A-Sepharose, Amersham Biosciences; 0.5 mg/mL bovine serum albumin, 0.2 mg/mL salmon sperm DNA) at 4°C for 1 h. Samples were spun, and the supernatants were incubated at 4°C for overnight with no antibody, IgG, anti-HIF-2α antibody. Immune complexes were recovered by adding 50 μL of blocked protein A beads and incubated overnight at 4°C. Beads were successively washed with: (1) mixed micelle buffer (20 mM Tris, pH 8.1, 150 mM NaCl, 5 mM EDTA, 5% w/v sucrose, 1% Triton X-100, 0.2% SDS); (2) buffer 500 (50 mM Hepes, pH 7.5, 0.1% w/v deoxycholic acid, 1% Triton X-100, 500 mM NaCl, 1 mM EDTA); (3) LiCl detergent wash buffer (10 mM Tris, pH 8,0.5% deoxycholic acid, 0.5% Nonidet P-40, 250 mM LiCl, 1 mM EDTA); and (4) TE buffer (10 mM Tris, 1 mM EDTA) and then eluted with 1% SDS and 0.1 M NaHCO3. Twenty milliliters of 5 M NaCl was added to the elutes, and the mixture was incubated at 65°C for 5 h to reverse the cross-linking. After digestion with proteinase K, the solution was phenol/chloroform-extracted and ethanol-precipitated. For qChIP analysis , DNA fragments were resuspended in 400 μL of water and 5 μL was used by real-time PCR. Each sample was calculated as the percentage of input sample. The sequences of PCR regions and primers used in ChIP assay are listed (Additional file 20). The antibodies used are listed (Additional file 17).
For immunofluorescence staining ,, cells on glass coverslips or chamber slides were fixed with 4% paraformaldehyde and permeabilized with 0.5% Triton X-100. After washing three times for 10 min with PBS, fixed cells or slides were blocked with blocking buffer (PBS, 0.1% Tween-20, with 3% goat serum) for 1 h and incubated with primary antibody diluted in blocking buffer overnight at 4°C. After washing three times with PBS for 10 min, the fixed cells were treated with the appropriate secondary antibody (fluorescein isothiocyanate (FITC)-conjugated anti-mouse IgG or rhodamine-conjugated anti-rabbit IgG) (Sigma) that was diluted in blocking buffer for 1 h at room temperature. Finally, the fixed cells were washed three times for 10 min with TBS, and their nuclei were counterstained, mounted, and observed by using fluorescence microscope or confocal microscope.
Lentivirus containing short hairpin RNAs (shRNAs) expressed in a lentiviral vector (pLKO.1-puro) were generated in 293 T cells as previously described . Packaging plasmid pCMVΔR8.91 was obtained from SC Teng (National Taiwan University, Taiwan). Various pLKO plasmids to knockdown HIF-1α, HIF-2α, TET1, INSIG1, and scrambled control were provided by National RNAi Core Facility of Academia Sinica, Taipei, Taiwan. For lentivirus production, 293 T cells were transfected with 15 μg pLKO.1-puro lentiviral vectors expressing different shRNAs along with 1.5 μg of envelope plasmid pMD.G and 15 μg of packaging plasmid pCMVΔR8.91. Virus was collected 48 h after transfection. To prepare various knockdown cells, FADU or H1299 cells were infected with lentivirus for 24 h, and used for various assays after 4 to 5 days. The sequence of the lentiviral siRNA vectors were shown in Additional file 21.
RNA-Seq data analysis
The RNA-Seq libraries were prepared according to the standard Illumina protocol with the mRNAseq Illumina TruSeq and were sequenced using Illumina Hiseq2000 to obtain 100 bp paired-end reads. The reads were aligned to the hg19 reference assembly with the TopHat/Cufflinks  alignment package using Ensembl annotations. Transcript abundance was measured in fragments per kb of exon per million fragments mapped (FPKM).
5hmC chemical labeling
5hmC labeling reactions were performed according to the previous protocol . Briefly, sonicated genomic DNA (average 400 bp, 500 ng/μL) was incubated with 50 mM HEPES buffer (pH 7.9), 25 mM MgCl2, 100 mM UDP-6-N3-Glc, and 2 mM βGT for 1 h at 37°C. The labeled DNA was purified by the QIAquick Nucleotide Removal kit (QIAGEN) and eluted in H2O. The click chemistry was performed with the addition of 150 mM of disulfide-biotin, and the mixture was incubated for 2 h at 37°C. The labeled DNA fragments were then purified by the QIAquick Nucleotide Removal kit (QIAGEN) and enriched by Dynabeads Streptavidin C1 (Invitrogen), and subsequently released by DTT treatment. The enriched DNA fragments were first purified by Micro Bio-Spin 6 spin columns (Bio-Rad) followed by MinElute PCR Purification Kit (QIAGEN).
5hmC-Seq data analysis
All the 5hmc-enriched DNAs were sequenced using Illumina Hiseq2000 containing more than 2.5 × 107, 50 bp single-end reads per sample. Raw sequence reads were mapped onto the reference human genome (NCBI Build UCSC hg19) using the Bowtie v0.12.7 algorithm (-m 1 -v 3 --best --strata) . Unique and monoclonal reads were used for further analysis. The distribution of 5hmC reads at promoters or in gene body regions were analyzed by PAVIS  (Additional file 6a). Further analysis was done using Bioconductor using packages ChIPpeakAnno, biomaRt, MEDIPS and in-house scripts. The genes with 5hmC peaks at promoter region (approximately -3000 to 0 bp upstream of the transcription start site) were chosen for further analysis.
Selection of overlapping groups of genes from RNA-Seq and 5hmC-Seq data followed by Gene Ontology (GO) analysis
The genes that have differential 5hmC-enriched promoter regions from 5hmC-Seq analysis and are differentially regulated from RNA-Seq analysis are selected for functional enrichment analysis. For 5hmC-Seq data, the criteria for genes with 5hmC-enriched promoter regions are: the peak levels of Hypoxia (scrambled control) >4X Hypoxia (TET1-si), Hypoxia (scrambled control) >2X Normoxia (scrambled control), and Normoxia (scrambled control) >2X Normoxia (TET1-si). The differentially expressed genes from RNA-Seq are identified by the criteria: Hypoxia (scrambled control) >2X Hypoxia (TET1-si), Hypoxia (scrambled control) > Normoxia (scrambled control), and Normoxia (scrambled control) > or = Normoxia (TET1-si). The 98 genes that were chosen through the overlapping criteria by 5hmC-Seq and RNA-Seq analyses were further analyzed by GO-based functional enrichment analysis. GO-based functional enrichment analysis (Fisher’s Exact test) was used to measure the gene-enrichment in annotation terms for selected genes. The genes in the GO terms that passed the criteria of P value <0.01 and at least two genes in each GO term were considered for further analysis.
Detection of 5hmC levels within specific gene locus (hMeDIP)
Genomic DNA was prepared using a genomic DNA extraction kit (QIAGEN, Germany). Detection of 5-hydroxymethylcytosine (5hmC) within specific gene locus was performed by using EpiJET 5-hmC Analysis Kit (Thermo Scientific, USA). The DNA was glucosylated and digested with the Glc-5-hmC sensitive restriction endonuclease Epi MspI. Glucosylation/digestion levels were analyzed by qPCR. The sequences of oligonucleotides used for hMeDIP were shown in Additional file 22.
Methylated DNA immunoprecipitation (MeDIP)
Genomic DNA was prepared using a genomic DNA extraction kit (QIAGEN, Germany) and sonicated with Bioruptor (Diagenode Inc., USA) to produce random fragments ranging in mean size from 300 to 1,000 bp. Two micrograms of fragmented DNA was denatured for 10 min at 95°C and immunoprecipitated for 16 h at 4°C with 5 μL of 5-methylcytidine antibody (Eurogentec, Germany) in a final volume of 250 μL IP buffer (20 mM Tri-HCl pH 7.5, 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 2.5 mM sodium pyrophosphate, 1 mM β-glycerophosphate, 0.05% Triton X-100). The mixture was incubated with 50 μL protein A beads for another 4 h at 4°C and washed three times with 1 mL of IP buffer. Beads were resuspended with 500 μL digestion buffer (50 mM Tris (pH 8.0), 10 mM EDTA, 1% SDS) containing 10 μL proteinase K (20 mg/mL stock) overnight at 65°C. DNA was phenol/chloroform-extracted, ethanol-precipitated, and resuspended in 200 μL of water for real-time PCR analysis. The sequences of oligonucleotides used for MeDIP were shown in Additional file 22.
Yeast strains, transformations, and media
PJ69-4A (MATa trp1-901 leu2-3,112 ura3-52 his3-200 gal4Δ gal80Δ LYS2∷GAL1-HIS3 GAL2-ADE2 met2∷GAL7-lacZ) was used as the yeast strain. Media was prepared as described . All yeast transformations were conducted using the high efficiency method described .
Co-immunoprecipitation was performed by incubating antibody with 500 μL of whole cell extracts in IP buffer (20 mM Tri-HCl pH 7.5, 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 1% Triton X-100, 2.5 mM sodium pyrophosphate, 1 mM β-glycerophosphate, and protease inhibitor) with 50 μL of protein A/G mix magnetic beads (Millipore) for 4 h at 4°C. The immunoprecipitates were washed several times with IP buffer to removing any non-bound proteins. Finally, components of the bound immune complex are eluted, and then analyzed by western blot detection to verify the identity of proteins.
Measurement of transcriptional activation strength
The activation strength was measured by two different assays which measure the expression of the reporter genes: HIS3, which encodes imidazoleglycerol-phosphate dehydratase and catalyzes the sixth step in yeast histidine biosynthesis. The pGBDU constructs containing HIF-1α, TET1, or TET1 truncation mutants were generated according to the length and restriction sites mentioned in Additional file 3. Plasmids were transformed into the yeast strain PJ69-4A . Transformed clones were selected on SC-Ura plates and confirmed for HIS3 expression on SC-His plates. Five to 10 freshly transformed colonies were mixed and spotted in five-fold serial dilutions onto plates. Plates were kept at 30°C until colonies formed.
The independent Student’s t-test was used to compare the continuous variables between two groups, and the χ 2 test was applied for comparison of dichotomous variables. The control groups of all the statistical analyses were specified in the figure legends. The level of statistical significance was set at 0.05 for all tests.
Sequencing data have been deposited to the Gene Expression Omnibus (GEO) under accession number GSE59990.
Correspondence and requests for materials should be addressed to KJW (firstname.lastname@example.org).
We are grateful to Dr. LJ Juan at Academia Sinica for the gift of TET1 wild-type and mutant cDNA expression vectors. We thank Dr. Ting-Yu Chang, WM Kung, and YH Lee for technical assistance. We declare no competing financial interests. This work was supported in part to KJW by Ministry of Science and Technology Summit grant (MOST 103-2745-B-039-001-ASP), National Science Council Frontier grant (NSC102-2321-B-010-001), National Research Program for Biopharmaceuticals (NSC102-2325-B-010-004), a grant from Ministry of Education, Aim for the Top University Plan (103 AC-T301), center of excellence for cancer research at Taipei Veterans General Hospital (MOHW103-TD-B-111-02), Taichung Veterans General Hospital (TCVGH-YM1010301), Kee-Lung Chang-Gung Memorial Hospital (CMRPG2D0031), and National Health Research Institutes (NHRI-EX102-10230SI, NHRI-EX103-10230SI). HC is an investigator of the Howard Hughes Medical Institute.
- Majmundar AJ, Wong WJ, Simon MC: Hypoxia-inducible factors and the response to hypoxic stress. Mol Cell. 2010, 40: 294-309. 10.1016/j.molcel.2010.09.022.PubMedPubMed CentralView ArticleGoogle Scholar
- Semenza GL: Hypoxia-inducible factors in physiology and medicine. Cell. 2012, 148: 399-408. 10.1016/j.cell.2012.01.021.PubMedPubMed CentralView ArticleGoogle Scholar
- Semenza GL: Hypoxia-inducible factors: mediators of cancer progression and targets for cancer therapy. Trends Pharmacol Sci. 2012, 33: 207-214. 10.1016/j.tips.2012.01.005.PubMedPubMed CentralView ArticleGoogle Scholar
- Lu X, Kang Y: Hypoxia and hypoxia-inducible factors: master regulators of metastasis. Clin Cancer Res. 2011, 16: 5928-5935. 10.1158/1078-0432.CCR-10-1360.View ArticleGoogle Scholar
- Thiery JP, Acloque H, Huang RY, Nieto MA: Epithelial-mesenchymal transitions in development and disease. Cell. 2009, 139: 871-890. 10.1016/j.cell.2009.11.007.PubMedView ArticleGoogle Scholar
- Yang MH, Wu MZ, Chiou SH, Chen PM, Chang SY, Liu CJ, Teng SC, Wu KJ: Direct regulation of TWIST by HIF-1α promotes metastasis. Nat Cell Biol. 2008, 10: 295-305. 10.1038/ncb1691.PubMedView ArticleGoogle Scholar
- Johnson AB, Denko N, Barton MC: Hypoxia induces a novel signature of chromatin modifications and global repression of transcription. Mutat Res. 2008, 640: 174-179. 10.1016/j.mrfmmm.2008.01.001.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsai YP, Wu KJ: Epigenetic regulation of hypoxia-responsive gene expression: focusing on chromatin and DNA modifications. Int J Cancer. 2014, 134: 249-256. 10.1002/ijc.28190.PubMedView ArticleGoogle Scholar
- Wu CY, Tsai YP, Wu MZ, Teng SC, Wu KJ: Epigenetic reprogramming and post-transcriptional regulation during the epithelial-mesenchymal transition. Trends Genet. 2012, 28: 454-463. 10.1016/j.tig.2012.05.005.PubMedView ArticleGoogle Scholar
- Javaid S, Zhang J, Anderssen E, Black JC, Wittner BS, Tajima K, Ting DT, Smolen GA, Zubrowski M, Desai R, Maheswaran S, Ramaswamy S, Whetstine JR, Haber DA: Dynamic chromatin modification sustains epithelial-mesenchymal transition following inducible expression of Snail-1. Cell Rep. 2013, 5: 1679-1689. 10.1016/j.celrep.2013.11.034.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu MZ, Tsai YP, Yang MH, Huang CH, Chang SY, Chang CC, Teng SC, Wu KJ: Interplay between HDAC3 and WDR5 is essential for hypoxia-induced epithelial-mesenchymal transition. Mol Cell. 2011, 43: 811-822. 10.1016/j.molcel.2011.07.012.PubMedView ArticleGoogle Scholar
- McDonald OG, Wu H, Timp W, Doi A, Feinberg AP: Genome-scale epigenetic reporgramming during epithelial-mesenchymal transition. Nat Struc Mol Biol. 2011, 18: 867-874. 10.1038/nsmb.2084.View ArticleGoogle Scholar
- Tahiliani M, Koh KP, Shen Y, Pastor WA, Bandukwala H, Brudno Y, Agarwal S, Iyer LM, Liu DR, Aravind L, Rao A: Conversion of 5-methylcytosine to 5-hdroxymethylcytosine in mammalian DNA by MLL partner TET1. Science. 2009, 324: 930-935. 10.1126/science.1170116.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu H, Zhang Y: Mechanisms and functions of Tet protein-mediated 5-methycytosine oxidation. Genes Dev. 2011, 25: 2436-2452. 10.1101/gad.179184.111.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu H, D’Alessio AC, Ito S, Xia K, Wang Z, Cui K, Zhao K, Sun YE, Zhang Y: Dual functions of Tet1 in transcriptional regulation in mouse embryonic stem cells. Nature. 2011, 473: 389-393. 10.1038/nature09934.PubMedPubMed CentralView ArticleGoogle Scholar
- Ito S, Shen L, Dai Q, Wu SC, Collins LB, Swenberg JA, He C, Zhang Y: Tet proteins can convert 5-methylcytosine to 5-formylcytosine and 5-carboxycytosine. Science. 2011, 333: 1300-1303. 10.1126/science.1210597.PubMedPubMed CentralView ArticleGoogle Scholar
- He YF, Li BZ, Li Z, Liu P, Wang Y, Tang Q, Ding J, Jia Y, Chen Z, Li L, Sun Y, Li X, Dai Q, Song CX, Zhang K, He C, Xu GL: Tet-mediated formation of 5-carboxycytosine and its excision by TDG in mammalian DNA. Science. 2011, 333: 1303-1307. 10.1126/science.1210944.PubMedPubMed CentralView ArticleGoogle Scholar
- Koh KP, Yabuuchi A, Rao S, Huang Y, Cunniff K, Nardone J, Laiho A, Tahiliani M, Sommer CA, Mostoslavsky G, Lahesmaa R, Orkin SH, Rodig SJ, Daley GQ, Rao A: Tet1 and Tet1 regulates 5-hydroxymethylcytosine production and cell lineage specification in mouse embryonic cells. Cell Stem Cell. 2011, 8: 200-213. 10.1016/j.stem.2011.01.008.PubMedPubMed CentralView ArticleGoogle Scholar
- Song CX, Szulwach KE, Fu Y, Dai Q, Yi C, Li X, Li Y, Chen CH, Zhang W, Jian X, Wang J, Zhang L, Looney TJ, Zhang B, Godley LA, Hicks LM, Lahn BT, Jin P, He C: Selective chemical labeling reveals the genome-wide distribution of 5-hydroxymethylcytosine. Nat Biotechnol. 2011, 29: 68-72. 10.1038/nbt.1732.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang T, Espenshade PJ, Wright ME, Yabe D, Gong Y, Aebersold R, Goldstein JL, Brown MS: Crucial step in cholesterol homeostasis: sterols promotes binding of SCAP to INSIG-1, a membrane protein that facilitates retention of SREBPs in ER. Cell. 2002, 110: 489-500. 10.1016/S0092-8674(02)00872-3.PubMedView ArticleGoogle Scholar
- Sever N, Yang T, Brown MS, Goldstein JL, DeBose-Boyd RA: Accelerated degradation of HMG CoA reductase mediated by binding of insig-1 to its sterol-sensing domain. Mol Cell. 2003, 11: 25-33. 10.1016/S1097-2765(02)00822-5.PubMedView ArticleGoogle Scholar
- Choi SM, Cho HJ, Cho H, Kim KH, Kim JB, Park H: Stra13/DEC1 and DEC2 inhibit sterol regulatory element binding protein-1c in a hypoxia-inducible factor-dependent mechanism. Nucl Acids Res. 2008, 36: 6372-6385. 10.1093/nar/gkn620.PubMedPubMed CentralView ArticleGoogle Scholar
- Nishiyama Y, Goda N, Kanai M, Niwa D, Osanai K, Yamamoto Y, Senoo-Matsuda N, Johnson RS, Miura S, Kabe Y: Suematsu M.HIF-1α induction suppresses excessive lipid accumulation in alcoholic fatty liver in mice. J Hepatol. 2012, 56: 441-447. 10.1016/j.jhep.2011.07.024.PubMedView ArticleGoogle Scholar
- Semenza GL, Roth PH, Fang HM, Wang GL: Transcriptional regulation of genes encoding glycolytic enzymes by hypoxia-inducible factor 1. J Biol Chem. 1994, 269: 23757-23763.PubMedGoogle Scholar
- Meijer TW, Kaanders JH, Span PN, Bussink J: Targeting hypoxia, HIF-1, and tumor glucose metabolism to improve radiation efficacy. Clin Cancer Res. 2012, 18: 1327-1338. 10.1158/1078-0432.CCR-12-0858.View ArticleGoogle Scholar
- Hsu CH, Peng KL, Kang ML, Chen YR, Yang YC, Tsai CH, Chu CS, Jeng YM, Chen YT, Lin FM, Huang HD, Lu YY, Teng YC, Lin ST, Lin RK, Tang FM, Lee SB, Hsu HM, Yu JC, Hsiao PW, Juan LJ: TET1 suppresses cancer invasion by activating the tissue inhibitors of metalloproteinases. Cell Rep. 2012, 2: 568-579. 10.1016/j.celrep.2012.08.030.PubMedView ArticleGoogle Scholar
- Vella P, Scelfo A, Jammula S, Chiacchiera F, Williams K, Cuomo A, Roberto A, Christensen J, Bonaldi T, Helin K, Pasini D: Tet proteins connec the O-linked N-acetylglucosamine transferase Ogt to chromatin in embryonic stem cells. Mol Cell. 2013, 49: 645-656. 10.1016/j.molcel.2012.12.019.PubMedView ArticleGoogle Scholar
- Shi FT, Kim H, Lu W, He Q, Liu D, Goodell MA, Wan M, Songyang Z: Ten-eleven translocation 1 (Tet1) is regulated by O-linked N-acetylglucosamine transferase (Ogt) for target gene repression in mouse embryonic stem cells. J Biol Chem. 2013, 288: 20776-20784. 10.1074/jbc.M113.460386.PubMedPubMed CentralView ArticleGoogle Scholar
- Shaw RJ: Glucose metabolism and cancer. Curr Opin Cell Biol. 2006, 18: 598-608. 10.1016/j.ceb.2006.10.005.PubMedView ArticleGoogle Scholar
- Luo W, Hu H, Chang R, Zhong J, Knabel M, O’Meally R, Cole RN, Pandey A, Semenza GL: Pyruvate kinase M2 is a PHD3-stimulated coactivator for hypoxia-inducible factor 1. Cell. 2011, 145: 732-744. 10.1016/j.cell.2011.03.054.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu Y, Xu C, Kato A, Tempel W, Abreu JG, Bian C, Hu Y, Hu D, Zhao B, Cerovina T, Diao J, Wu F, He HH, Cui Q, Clark E, Ma C, Barbara A, Veenstra GJ, Xu G, Kaiser UB, Liu XA, Sugrue SP, He X, Min J, Kato Y, Shi YG: Tet3 CXXC domain and dioxygenase activity cooperatively regulate key genes for Xenopus eye and neural development. Cell. 2012, 151: 1200-1213. 10.1016/j.cell.2012.11.014.PubMedPubMed CentralView ArticleGoogle Scholar
- Costa Y, Ding J, Theunissen TW, Faiola F, Hore TA, Shliaha PV, Fidalgo M, Saunders A, Lawrence M, Dietmann S, Das S, Levasseur DN, Li Z, Xu M, Reik W, Silva JC, Wang J: NANOG-dependent function of TET1 and TET2 in establishment of pluripotency. Nature. 2013, 495: 370-374. 10.1038/nature11925.PubMedPubMed CentralView ArticleGoogle Scholar
- Williams K, Christensen J, Pedersen MT, Johansen JV, Cloos PA, Rappsilber J, Helin K: TET1 and hydroxymethylcytosine in transcription and DNA methylation fidelity. Nature. 2011, 473: 343-348. 10.1038/nature10066.PubMedPubMed CentralView ArticleGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L: Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012, 7: 562-578. 10.1038/nprot.2012.016.PubMedPubMed CentralView ArticleGoogle Scholar
- Song C, Clark TA, Lu XY, Kislyuk A, Dai Q, Turner SW, He C: Sensitive and specific single molecule sequencing of 5-hydroxymethylcytosine. Nat Methods. 2011, 9: 75-77. 10.1038/nmeth.1779.PubMedPubMed CentralView ArticleGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang W, Loganantharaj R, Schroeder B, Fargo D, Li L: PAVIS: a tool for Peak Annotation and Visualization. Bioinformatics. 2013, 29: 3097-3309. 10.1093/bioinformatics/btt520.PubMedPubMed CentralView ArticleGoogle Scholar
- Rose MD, Winston F, Hieter P: Methods in Yeast Genetics. 1990, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NYGoogle Scholar
- Gietz RD, Schiestl RH: Transforming yeast with DNA. Methods Mol Cell Biol. 1995, 1995: 255-269.Google Scholar
- James P, Halladay J, Craig EA: Genomic libraries and a host strain designed for highly efficient two-hybrid selection in yeast. Genetics. 1996, 144: 1425-1436.PubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.