- Open Access
TET1 regulates hypoxia-induced epithelial-mesenchymal transition by acting as a co-activator
Genome Biology volume 15, Article number: 513 (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
TET enzymes have been shown to convert 5mC to 5hmC to regulate gene expression -. In spite of the various epigenetic mechanisms demonstrated to regulate hypoxia-responsive gene expression ,, the role of TET enzymes and DNA demethylation in regulating hypoxia-responsive genes remain largely unknown. We tested whether TET1 could be regulated by hypoxia to mediate hypoxia-regulated gene expression. Exposure of various cell lines to hypoxia showed the activation of TET1 mRNA expression (Additional file 1a). Western blot analysis confirmed the upregulation of TET1 protein levels by hypoxia (Figure 1a). Overexpression of HIF-2α, but not HIF-1α, activated the expression of TET1 (Figure 1b and data not shown). Knockdown of HIF-2α abolished the activation of TET1 under hypoxia in two different cell lines (Figure 1c and Additional file 1b and c), indicating that HIF-2α was the major regulator of TET1 expression under hypoxia. We further identified the region in the proximal promoter of TET1 gene that responded to hypoxia and HIF-2α. Reporter gene assay showed that the promoter region (-381 to +17 bp upstream of the transcription start site, TSS) of the TET1 gene was activated by hypoxia/HIF-2α (Figure 1d and Additional file 2). The hypoxia response region was further narrowed down (-158 to +17 bp upstream of the TSS) (Figure 1d). The construct containing the promoter region of -91 to +17 bp upstream of the TSS did not respond to hypoxia/HIF-2α stimulation (data not shown), further narrowing down the hypoxia/HIF-2α responsive region to -158 to -91 bp upstream of the TSS of TET1 gene. Chromatin immunoprecipitation experiments showed that HIF-2α directly bound to the proximal promoter region (Additional file 1d). Positive control also showed that HIF-2α bound to the promoter region of WDR5 gene (Additional file 1d) . All these results indicate the direct regulation of TET1 gene expression by HIF-2α.
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
Since TET1 demethylates 5mC to regulate gene expression ,, we wanted to determine the genes whose expression was regulated by TET1 under hypoxia. RNA sequencing experiments using FADU control vs. FADU-TET1 knockdown cells under normoxia or hypoxia showed that a total of 1,044 genes were regulated by TET1 under hypoxia since their expression decreased under TET1 knockdown (Figure 2a and b). Transcriptome profiles were shown in Figure 2a. Gene Ontology (GO) analysis of these genes showed that they belonged to various biological processes (Additional file 5a). To further narrow down the genes that were regulated by TET1 and also showed changes in the levels of 5hmC in their promoter regions, 5hmC sequencing was performed on these FADU clones under normoxia or hypoxia (control vs. TET1 knockdown) . Analysis of 5hmC sequencing results showed the distribution of 5hmC peaks (Additional file 5b). There were no genome-wide changes of 5hmC labeling with regards to their genomic locations under normoxia or hypoxia (Additional file 6a), indicating that 5hmC regulated by TET1 only influenced the expression of certain specific genes. To search for genes that were regulated by TET1 from RNA expression analysis and also had increased levels of 5hmC in their promoters, Venn diagram showed that 98 genes had such patterns (Figure 2b). GO analysis of the 98 genes categorized them into different groups (Figure 2c). It is interesting that the genes belonging to cholesterol metabolic process represented the most statistically significant group (Figure 2c), which was also one of the prominently represented groups from RNA-Seq analysis (Additional file 5a). INSIG1 (insulin induced gene 1), a major regulator of cholesterol biosynthesis ,, was inside this group (Figure 2a and c). INSIG1 sequesters SREBP in ER to inhibit cholesterol synthesis and promotes the degradation of HMG CoA reductase ,. We further focused on the analysis of INSIG1 gene since HIF-1α inhibits sterol regulatory element binding protein-1c (SREBP-1c) and suppresses excessive lipid accumulation ,. Analysis of INSIG1 RNA levels confirmed the activation of INSIG1 expression by hypoxia and knockdown of TET1 decreased the expression of INSIG1 under hypoxia in two different cell lines using quantitative real-time PCR analysis (Figure 2d and e and Additional file 6b and c).
Demethylation of the INSIG1promoter regions by TET1
To correlate the RNA expression and 5hmC levels in the promoter region regulated by TET1, the group of cholesterol metabolic process genes including INSIG1, FDPS (farnesyl diphosphate synthase), APOC1 (apolipoprotein C1), and SQLE (squalene epoxidase) were analyzed. There were increased 5hmC peaks in their promoters (Figure 3a and Additional file 7). Since the expression of INSIG1 was confirmed (Figure 2d and e and Additional file 6b and c), we further analyzed the 5hmC peaks in the promoter of INSIG1. The 5hmC peaks in two different promoter regions of INSIG1 gene were confirmed to be increased by hypoxia and decreased under TET1 knockdown in two different cell lines using hMeDIP assays (Figure 3b and Additional file 8). In addition, the levels of 5mC also decreased in the promoter region of INSIG1 under hypoxia in two different cell lines using MeDIP assays (Figure 3c and Additional file 9). Finally, to test the role of INSIG1 in hypoxia-induced EMT, knockdown of INSIG1 was performed. The result showed that knockdown of INSIG1 abolished hypoxia-induced EMT in FADU and H1299 cells, mainly abolishing mesenchymal gene activation (Figure 3d and Additional file 10a). Increased in vitro migration and invasion activity in H1299 cells induced by hypoxia was abolished under INSIG1 knockdown (Additional file 10b). Further analysis showed that the abolishment of hypoxia-induced EMT by knockdown of INSIG1 was not mediated through repression of EMT regulators including Twist1 and Snail (Additional file 10c). All the above results indicate that INSIG1 plays a significant role in hypoxia-induced EMT, which is regulated by TET1.
Regulation of glucose metabolism gene expression by TET1 or INSIG1
We tested whether TET1 or INSIG1 regulates glucose metabolism pathway to possibly contribute to hypoxia-induced Warburg effect ,. The genes in the glucose metabolism pathway were tested for their expression under hypoxia and knockdown of TET1 or INSIG1. The results showed that the genes including glucose transporter 3 (GLUT3), hexokinase 1 (HK1), phosphoglycerate kinase 2 (PGK2), pyruvate kinase M (PKM), and lactate dehydrogenase A (LDHA) were induced by hypoxia and knockdown of TET1 or INSIG1 abolished or decreased their expression levels in two different cell lines (Figure 4a,b, and Additional file 11a,b). In contrast, phosphoglycerate kinase 1 (PGK1) was induced by hypoxia but did not respond to TET1 or INSIG1 knockdown (Figure 4c,d and Additional file 11c,d). These results indicate that TET1 and INSIG1 regulate the expression of hypoxia-induced glycolytic genes and may contribute to Warburg effect induced by hypoxia.
TET1 acts as a transcription co-activator to regulate hypoxia-induced gene expression
Since there were TET1-regulated genes that did not show difference in the 5hmC levels in their promoters (Figure 2b), we speculated that TET1 may serve as a co-regulator together with HIFs to regulate gene expression. To test this hypothesis, reporter gene assays using HIF-1α, TET1, CBP (CREB binding protein), and Twist1 promoter-driven reporter construct were performed in 293 T cells. Our previous results indicated that wild-type HIF-1α could be stably expressed in 293 T cells due to the overexpression of HIF-1α that is beyond the HIF-1α degradation activity under normoxia ,. Indeed TET1 further enhanced the reporter gene activity similar to CBP (Figure 5a). In addition, the TET1 enzymatically inactive mutant (TET1-CDmt) also had similar activity to increase the ability of HIF-1α to promote Twist1 promoter activity (Figure 5a) . Similar results were shown using the WDR5 promoter-driven reporter construct (Additional file 12a). In addition, using the INSIG1 promoter-driven reporter gene, the ability of HIF-2α to promote its activity was further enhanced by either TET1 wild-type or TET1-CDmt protein (Additional file 12b and c). These results indicate that TET1 may serve as a co-activator to enhance the transactivation activity of either HIF-1α or HIF-2α.
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
Since TET1 knockdown mitigated hypoxia-induced EMT and TET1 can act as a transcriptional co-activator without the enzymatic activity, we tested whether a TET1 catalytically inactive mutant could rescue hypoxia-induced EMT in TET1 knockdown cells. The TET1-CDmt expression vector was transfected into two different TET1 knockdown cell lines to test its ability to rescue hypoxia-induced EMT. The results showed that TET1-CDmt was capable of restoring hypoxia-induced EMT in two different TET1 knockdown cell lines (Figure 6a,b), indicating that the catalytic activity of TET1 is dispensable for hypoxia-induced EMT.
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).
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.
Semenza GL: Hypoxia-inducible factors in physiology and medicine. Cell. 2012, 148: 399-408. 10.1016/j.cell.2012.01.021.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Shaw RJ: Glucose metabolism and cancer. Curr Opin Cell Biol. 2006, 18: 598-608. 10.1016/j.ceb.2006.10.005.
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.
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.
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.
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.
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.
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.
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.
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.
Rose MD, Winston F, Hieter P: Methods in Yeast Genetics. 1990, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY
Gietz RD, Schiestl RH: Transforming yeast with DNA. Methods Mol Cell Biol. 1995, 1995: 255-269.
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.
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.
The authors declare that they have no competing interests.
KJW conceived and initiated the study. YPT, HFC, and SYC performed the experiments. CWC and HWW analyzed the bioinformatics data. ZJS performed the yeast experiments. CS performed the 5hmC pull down assays. SCT, CH, and KJW analyzed the data. SCT and KJW wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 4: Figure S3.: Knockdown of TET1 mitigated hypoxia-induced epithelial-mesenchymal transition in FADU cells, immunofluorescence staining of E-cadherin and vimentin, and in vitro migration and invasion activity of FADU cell lines with scrambled or TET1 knockdown under normoxia or hypoxia. (DOC 202 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.