- Open Access
SOX4 interacts with EZH2 and HDAC3 to suppress microRNA-31 in invasive esophageal cancer cells
Molecular Cancer volume 14, Article number: 24 (2015)
Tumor metastasis is responsible for 90% of cancer-related deaths. Recently, a strong link between microRNA dysregulation and human cancers has been established. However, the molecular mechanisms through which microRNAs regulate metastasis and cancer progression remain unclear.
We analyzed the reciprocal expression regulation of miR-31 and SOX4 in esophageal squamous and adenocarcinoma cell lines by qRT-PCR and Western blotting using overexpression and shRNA knock-down approaches. Furthermore, methylation studies were used to assess epigenetic regulation of expression. Functionally, we determined the cellular consequences using migration and invasion assays, as well as proliferation assays. Immunoprecipitation and ChIP were used to identify complex formation of SOX4 and co-repressor components.
Here, we report that SOX4 promotes esophageal tumor cell proliferation and invasion by silencing miR-31 via activation and stabilization of a co-repressor complex with EZH2 and HDAC3. We demonstrate that miR-31 is significantly decreased in invasive esophageal cancer cells, while upregulation of miR-31 inhibits growth, migration and invasion of esophageal adenocarcinoma (EAC) and squamous cell carcinoma (ESCC) cell lines. miR-31, in turn, targets SOX4 for degradation by directly binding to its 3′-UTR. Additionally, miR-31 regulates EZH2 and HDAC3 indirectly. SOX4, EZH2 and HDAC3 levels inversely correlate with miR-31 expression in ESCC cell lines. Ectopic expression of miR-31 in ESCC and EAC cell lines leads to down regulation of SOX4, EZH2 and HDAC3. Conversely, pharmacologic and genetic inhibition of SOX4 and EZH2 restore miR-31 expression. We show that SOX4, EZH2 and HDAC3 form a co-repressor complex that binds to the miR-31 promoter, repressing miR-31 through an epigenetic mark by H3K27me3 and by histone acetylation. Clinically, when compared to normal adjacent tissues, esophageal tumor samples show upregulation of SOX4, EZH2, and HDAC3, and EZH2 expression is significantly increased in metastatic ESCC tissues.
Thus, we identified a novel molecular mechanism by which the SOX4, EZH2 and miR-31 circuit promotes tumor progression and potential therapeutic targets for invasive esophageal carcinomas.
MicroRNAs (miRNAs) are a class of highly-conserved, noncoding 18-25- nucleotide RNAs that function as negative regulators of gene expression at the post-transcription level, binding to the 3′-untranslated region (3′-UTR) of mRNAs transcripts and targeting them for degradation . Though implicated in carcinogenesis, it is not clear how miRNAs promote tumorigenesis and metastasis or what networks regulate miRNAs expression. miRNA expression is commonly dysregulated in human cancers, [2,3] including esophageal cancers .
miR-31 expression is altered in multiple human cancers. Depending on the cellular context, miR-31 may be up- or downregulated, acting as an oncogene or tumor suppressor, respectively. Overexpression of miR-31 has been linked to disease progression in colorectal cancer , head-and-neck squamous cell carcinoma (HNSCC)  and lung cancer . miR-31 is downregulated in certain T-cell leukemias , breast cancer [9,10], melanoma , ovarian cancer  and prostate cancer . Downregulation and loss of miR-31 in esophageal adenocarcinoma (EAC) correlates with poor patient prognosis [14-16]. Additionally, miR-31 expression is reduced in EAC patients with poor histomorphologic response to neoadjuvant chemoradiation therapy . Conversely, miR-31 is upregulated in serum and tissue samples of esophageal squamous cell carcinoma (ESCC), with expression correlating to staging . Yet, in another ESCC cohort miR-31 expression was decreased, and low miR-31 expression correlated with poorly differentiated tumors and decreased survival . These reports emphasize the complexity of miR-31-associated phenotypes and the need to better define miR-31 targets, as well as pathways regulating miR-31 expression in different cancers.
SOX4 is a member of the highly conserved SoxC (SRY-related high-motility group box) transcription factors family, which contains two other members, SOX11 and SOX12 . SOX4 is a putative stem cell marker that plays a crucial role during cell fate determination [21,22]. SOX4-deficient mice suffer from multiple developmental defects, dying at embryonic day 14, secondary to ventricular outflow tract malformation . During embryogenesis, SoxC members are highly expressed, helping to maintain survival of pluripotent mesenchymal and neural progenitor cells . In adults, expression of SOX4 is restricted to certain cell types, including hematopoietic stem cells, mammary stem cells and hair follicle stem cells [25-27].
Meta-analysis has identified SOX4 as one of the 64 genes that constitute a general signature in all human cancers, and genome wide promoter analysis has shown that SOX4 regulates the transcription of genes involved in TGF-β, Wnt, Hedgehog, and Notch pathways and components of miRNA processing machinery such as Dicer, Argonaute 1 and RNA Helicase A [28,29]. SOX4 induces EMT and breast cancer progression by cooperating with oncogenic Ras. More recent work shows that SOX4 induces EMT via the polycomb epigenetic regulator EZH2 . miRNAs, such as miR-335, are known to target SOX4, suppressing metastasis and migration in breast cancer .
Polycomb group proteins have been linked to tumor progression in many cancers. The polycomb proteins can form at least two complexes: polycomb-repressive complexes 1 and 2 (PRC1 and PRC2). PRC2 contains three core proteins, EZH2, SUZ12, and EED. The histone methyltransferase EZH2 (enhancer of zeste homolog 2) epigenetically regulates genes involved in cell fate determination. Specifically, EZH2 trimethylates nucleosomal histone H3 at lysine 27 (H3K27me3). The H3K27me3 mark is associated with gene silencing and often found in the promoter of developmental genes [31,32]. However, it is unclear how EZH2 is recruited to the promoters it targets. Recent studies have shown that EZH2 interacts with various transcription factors such as androgen receptor (AR), GATA4, RORα and STAT3 and may directly activate or repress these genes independent of the H3K27me3 mark or chromatin modification [33-36]. With respect to miRNAs, prior work demonstrates that EZH2 interacts with AR to silence miR-31 in prostate cancers, and C-MYC recruits EZH2 to the miR-29 promoter in B-cell lymphomas .
EZH2 is upregulated in multiple cancers, promoting invasion and metastasis [38-40]. Genetic and epigenetic loss of miR-31 is associated with EZH2 overexpression in melanoma , suggesting that miR-31 directly or indirectly regulates EZH2 expression. Interestingly, studies show that polycomb complexes silence the CDKN2A and CDKN2B loci, which encode the tumor suppressors p14 (ARF), p15 (INK4B) and p16 (INK4A)  and contain the MIR31HG locus on chromosome 9 . In line with this observation, Yamagishi et al. reported that PRC2 binds the miR-31 coding region and directly represses transcription of miR-31 in adult T-cell leukemia . SOX4 positively regulates EZH2, indicating a potential functional link between miR-31, EZH2 and SOX4.
The roles of SOX4, HDAC3 and EZH2 in microRNA regulation are largely unknown and have been poorly defined so far. In this study, we explore the role of SOX4 and EZH2 in miR-31 repression and the contribution of miR-31 to survival, migration and invasion of aggressive esophageal cancers cells. We identify SOX4 as a direct target of miR-31. Expression of miR-31 inhibits SOX4, EZH2 and HDAC3 expression. We show that miR-31 is repressed in invasive esophageal cancers cell lines and that miR-31 levels inversely correlate with SOX4, EZH2 and HDAC3 expression. Co-immunoprecipitation demonstrates that SOX4 interacts with EZH2 and that HDAC3 may be important to bridge this interaction. We show that EZH2 and HDAC3 bind to the miR-31 promoter using chromatin immunoprecipitation. Altogether, our results identify a feed-forward loop that leads to the activation of SOX4, which in turn up-regulates and binds to EZH2, cooperating with HDAC3 to repress the miR-31 promoter and advance esophageal tumorigenesis.
miR-31 expression is downregulated in invasive esophageal cancer cells
To investigate the role of miR-31 in esophageal cancers, we examined the expression of miR-31 in ESCC, EAC and Barrett’s esophagus cell lines of differing invasive potential (Figure 1). Comparing esophageal squamous cell carcinoma cell lines, TE11 is less motile than TE8 and displays an epithelial phenotype (Figure 1A). The esophageal adenocarcinoma cell lines OE33 and FLO1 differ in that FLO1 is more mesenchymal and therefore more motile than the OE33 (Figure 1B). After miR-200a and 200b, which are known for their roles in EMT, miR-31 was the most downregulated miRNA in invasive FLO1 cells compared to their less invasive OE33 counterparts by qPCR screen (Figure 1C). Similarly, miR-31 downregulation was observed in TE8 ESCC cell lines compared to TE11 (Figure 1D). Likewise, miR-31 expression was higher in non-invasive cell lines such as the benign Barrett’s esophagus cell line CP-A compared to the metaplastic CP-B cell line (Figure 1E). Furthermore, we confirmed the elevated miR-31 expression in OE33 cells, which have an epithelial phenotype, compared to FLO1 cells (Figure 1F). Next, to focus on the biological significance and regulatory mechanisms of miR-31 expression in invasive adenocarcinoma and squamous cell carcinoma, we expressed miR-31 in invasive ESCC and EAC cell lines and analyzed the effects on cell migration and invasion.
miR-31 suppresses migration and invasion of aggressive ESCC and EAC cells
To examine the functional contribution of miR-31 in aggressive esophageal cancer, we transfected TE8 and FLO1 cells with vectors containing the precursor of miR-31 or an empty vector control. Ectopic expression of precursor and mature miR-31 in the respective cell lines was tested by quantitative RT-PCR (Figure 2A). In Boyden chamber migration and invasion assays, precursor miR-31 transfection significantly decreased cell migration and invasion in TE8 and FLO1 cells (Figure 2B, C, respectively). miR-31 expression had no significant effect on proliferation in TE8 cells and did not alter the number of colonies in colony formation assays (Figure 2D, E). In FLO1 cells, however, miR-31 suppressed proliferation and colony formation (Figure 2D, E, respectively), indicating miR-31 regulates esophageal carcinoma cell growth in some cell lines. These data suggest that miR-31 suppresses esophageal cancer cell motility and invasiveness, but cell growth depending on the cellular context.
miR-31 is epigenetically repressed in invasive esophageal cancer cells
Prior studies report that miR-31 expression is epigenetically silenced through promoter hypermethylation at CpG islands, as well as polycomb-mediated histone methylation [8,11]. We therefore speculated that loss of miR-31 in invasive esophageal cancer cells could be mediated, in part, by DNA and histone methylation. To determine the effect of PRC2 on miR-31 expression, we utilized the PRC2 inhibitor, 3-deazaneplanocin (DZNep). DZNep treatment of TE8 and FLO1 cells resulted in a decrease in EZH2 expression and caused a dose-dependent increase in miR-31 expression in both of our invasive cell lines (Figure 3A, B, respectively). DZNep treatment also led to a decrease in SOX4 protein level in FLO1 cells (Figure 3A). SOX4 was recently shown to upregulate EZH2 expression . Using Western blot and qRT-PCR, we evaluated the effect of a pan-HDAC inhibitor (SAHA) on miR-31 expression. We found that SAHA led to a decrease in HDAC3, EZH2 and EZH1 protein levels in TE8 and FLO1 cells (Figure 3C). Most importantly, miR-31 expression was significantly upregulated in both cell lines following SAHA treatment (Figure 3D). To test whether promoter methylation at CpG islands was involved in miR-31 silencing, we used the DNA methylation inhibitor 5′AZA-Deoxy-Cytidine (AZA). Treatment with AZA significantly increased the expression of miR-31 in TE8 cells and to a lesser extent in FLO1 cells (Figure 3E and F). These data suggest that PRC2, HDAC and DNA methylation are involved in miR-31 epigenetic silencing.
miR-31 directly targets SOX4 and indirectly targets EZH2 and HDAC3
A study by Asangani et al. recently showed that genetic and epigenetic loss of miR-31 leads to a feed forward upregulation of EZH2 . However, no mechanism was proposed. Previously, EZH2 was reported to interact with HDAC3 to repress miR-29 in lymphomas . More recent work shows that SOX4 binds to the EZH2 promoter, thereby upregulating EZH2 expression . We hypothesized that SOX4 initiates the feed forward activation of EZH2, which in turn represses miR-31. Analysis of the SOX4 3′-UTR using microrna.org (maintained at cBio, the Computational Biology Center at Memorial Sloan-Kettering Cancer Center) predicted a miR-31 binding site (Figure 4A). A sequence alignment search showed that the miR-31 target sequence in the SOX4 3′-UTR is conserved in humans and most great apes (Figure 4A). To test whether SOX4 is regulated by miR-31 through direct binding to its 3′UTR, we used psiCHECK2 SOX4 full length 3′-UTR plasmid (WT) , and constructed two derivatives, SOX4 WT 3′-UTR oligo plasmid (WT OLIGO) and SOX4 mutant 3′-UTR oligo plasmid (MUT OLIGO) (Figure 4A). The WT OLIGO plasmid contained a 71-nucleotide region including the miR-31 target sequence. In the SOX4 mutant 3′-UTR (MUT OLIGO), 4 nucleotides in the seed sequence were mutated . When co-transfected into HEK-293 cells, the luciferase reporter, SOX4 WT 3′-UTR and miR-31 plasmid showed reduced luciferase activity compared to co-transfection with miR-31 empty control vector (Figure 4B). This suppressive effect was reversed by the four-nucleotide substitution in the miR-31 binding sequence. Similarly, the suppressive effect of miR-31 on the SOX4 3′-UTR activity was observed in the esophageal tumor cell lines, TE8 and FLO1 (Figure 4B). In line with these results, overexpression of miR-31 in FLO1 cells suppressed the expression of SOX4 at both the protein (Figure 4C) and mRNA level (Figure 4D). As previously reported , our data confirm that miR-31 inhibits EZH2 expression (Figure 4C and D) whereas EZH1 expression was unchanged. Interestingly, miR-31 decreased HDAC3 on protein (Figure 4C) and mRNA levels (Figure 4D). However, target prediction algorithms do not detect any putative binding site for miR-31 in the 5′UTR, 3′UTR or coding sequence of HDAC3. Taken together, these results demonstrate that SOX4 is a direct target of miR-31, while EZH2 and HDAC3 are indirect targets.
To analyze the regulation of miR-31, EZH2, EZH1 and HDAC3 by SOX4, we chose TE11 cells, which express low SOX4 and EZH2 levels, had high expression of miR-31 as a model. Ectopic expression of SOX4 in TE11 cells (Figure 4E) decreased miR-31 expression, although this was not statistically significant (Figure 4F). Additionally, we observed that SOX4 induced EZH2 expression (Figure 4F) as previously reported , but not EZH1. HDAC3 expression increased as well. Consistent with this observation, overexpression of EZH2 in TE11 cells (Figure 4G) led to a significant decrease in miR-31 expression (Figure 4H). Interestingly, EZH2 induced SOX4 expression by more than a 2-fold (Figure 4H), suggesting that not only does SOX4 regulate EZH2, but EZH2 also regulates SOX4 thereby potentially repressing miR-31 and/or other miRNAs and transcription factors.
Depletion of SOX4 suppresses growth, migration and invasion of esophageal cancer cells
Based on the observation that miR-31 targets both SOX4 and EZH2 and given that SOX4 binds to the EZH2 promoter to activate its transcription, we examined whether SOX4 leads to feed forward activation of EZH2, and subsequent miR-31 silencing. TE8 and FLO1 cells, which express high levels of SOX4 and EZH2, showed the lowest expression level of miR-31. SOX4 knockdown (Figure 5A) using shRNA led to significant upregulation of miR-31 in TE8 and FLO1 cells, while miR-191 and miR-423-5p, used as controls, did not show any significant change (Figure 5B and C). We next investigated if loss of SOX4 functionally mimics overexpression of miR-31 in esophageal cancer cells. To assess the role of SOX4, we studied the migratory and invasive potential of TE8 and FLO1 cells following SOX4 knockdown with shRNA. Suppression of SOX4 inhibited transwell migration (Figure 5D and E) and invasion through Matrigel-coated Boyden chambers (Figure 5F and G) in TE8 and FLO1 cells. SOX4 knockdown cells also showed a significant reduction in proliferation (Figure 5H and I). Taken together, these data implicate SOX4 as a key mediator of the tumor-suppressive effects of miR-31 in this system.
SOX4, HDAC3 and EZH2 form a co-repressor complex to inhibit miR-31 expression
PRC2 is known to recruit DNA methyl-transferases (DNMTs), HDACs and other chromatin modifying enzymes to repress the transcription of developmental genes . Because SOX4 was recently shown to interact with other transcription factors, we tested whether SOX4 forms a co-repressor complex with EZH2 and HDAC3 to silence miR-31 expression. Using anti-SOX4 antibody for co-immunoprecipitation, cell lysates from TE8 (Figure 6A) and FLO1 cells (Figure 6B) showed interactions between SOX4 and EZH2, HDAC3 and H3K27me3 (Figure 6A, B). We next performed co-immunoprecipitation assays with anti-EZH2 antibody and showed that EZH2 equally interacts with SOX4 (Figure 6A, B). Immunoprecipitates obtained with EZH2 antibody detected also HDAC3, H3K27me3, and Suz12 (Figure 6A, B). To test whether HDAC3 interacts with SOX4 and EZH2, HDAC3-antibody was used for pull down. Lysates immunoprecipitated with HDAC3-specific antibody contained EZH2, SOX4, H3K27me3 in TE8 cells (Figure 6A), but only a faint band of EZH2 in FLO1 cells (Figure 6B). Suz12 was not detected after pull-down with HDAC3-specific antibody in either cell line.
Finally, we tested whether SOX4, EZH2 and HDAC3 tethered to the miR-31 promoter. Previous reports identify CpG islands and histone trimethylation mark sites upstream of the miR-31 promoter . Therefore, we tested whether SOX4, EZH2 and HDAC3 bind to those regions. ChIP analysis denoted H3K27me3 and HDAC3 enrichment in regions upstream of miR-31 (Figure 6C). However, we were unable to detect SOX4 at miR-31 promoter regions (data not shown). HDAC3 and H3K27me3 were equally enriched at the MYT1 promoter, which was used as a positive control as EZH2 has been shown to bind and regulate MYT1 promoter activity . We conclude that SOX4, HDAC3 and EZH2 function as a potential co-repressor complex to silence miR-31.
SOX4, EZH2 and HDAC3 inversely correlate with miR-31 expression in invasive esophageal cancer cells
To determine the expression of the complex-forming partners, SOX4, EZH2 and HDAC3, in esophageal cancer cell lines and tissues, we used qRT-PCR. We profiled SOX4, EZH2 and HDAC3 in ESCC (Figure 7A) and EAC (data not shown) cell lines and found a strong, inverse correlation between miR-31 and expression of SOX4, EZH2 and HDAC3 in invasive cancers of both histologies. Probing various publicly available datasets (GSE20347, GSE47404, GSE13937) to assess the expression of SOX4, EZH2 and HDAC3 in primary esophageal squamous and adenocarcinoma samples, we found SOX4 and EZH2 to be significantly upregulated in ESCC tumor samples (GSE20347), compared to normal tissues (Figure 7B). In another ESCC dataset (GSE47404), only EZH2 expression (P value = 0.008) was significantly upregulated in metastatic compared to non-metastatic patients. In EAC dataset (GSE13937), we found that SOX4 (P value = 0.003) and HDAC3 (P value = 0.042) were increased in tumor samples compared to adjacent normal tissues (Figure 7C). Taken together, these data identify a molecular circuit where SOX4, EZH2 and HDAC3 target miR-31 to promote esophageal malignancy. Inversely, in non-invasive tumor cells miR-31 targets SOX4, EZH2 and HDAC3 by direct and indirect means to inhibit tumor cell invasion (Figure 7D).
Differential microRNAs expression has been linked to tumor initiation and progression. Depending on cellular context, microRNA expression in tumors may be increased or decreased, with microRNAs behaving as tumor promoters or suppressors, respectively . Moreover, multiple microRNAs have been shown to promote or inhibit metastasis [3,45]. Because metastasis is responsible for more than 90% of cancers-related deaths, it is important to define molecular mechanisms by which microRNAs regulate metastasis and define new therapeutic targets. Studies on the role of miR-31 in esophageal tumors have produced conflicting results [16,18,19,46]. Hence, it is important to define not only the molecular pathways regulated by miR-31 but also the factors regulating miR-31 expression and functions in various tissues and cancers. The lack of consistency in the literature with respect to miR-31 expression in esophageal cancer may be due to platform choice or normalization methods. Variability may also result due to the limited number of patients in each study. Here, we show downregulation of miR-31 in invasive and aggressive esophageal cancer cells. miR-31 has a known role in prostate cancer and melanoma, suppressing key cell cycle regulators and pro-oncogenic genes such as CDK1, E2F2, EXO1, FOXM1, MCM2, Src or MET . We show that miR-31 significantly suppresses migration and invasion in vitro in both aggressive esophageal adenocarcinoma and squamous cell carcinoma. Ectopic expression of miR-31 did not significantly affect tumor cell proliferation, causing only a marginal decrease in colony formation in the adenocarcinoma cell line FLO1. Interestingly, Valastyan et al. also reported that miR-31 promoted metastasis but not cell proliferation in breast cancer . A number of studies argue that epithelial-mesenchymal transition (EMT) plays a crucial role in cancer metastasis and progression, and loss of miR-200 family members drives EMT in multiple cancers [47,48]. We examined expression of EMT markers in miR-31-overexpressing cell lines but found no significant alterations (data not shown). We conclude that miR-31 is not a strong inhibitor of EMT.
Several mechanisms could contribute to aberrant miR-31 expression in cancer. Based on data demonstrating epigenetic repression of miR-31 expression by DNA methylation and EZH2-mediated H3K27me3 epigenetic mark in melanoma, leukemia and prostate cancer [8,11], we examined the role of epigenetic regulation of miR-31 in invasive esophageal cancer. Our results indicate that DNA methylation and Polycomb-mediated histone methylation both contribute to miR-31 silencing, since treatment with the EZH2 inhibitor, DZNep, and the DNMT inhibitor, AZA, enhanced miR-31 expression. Additionally, we show for the first time that histone deacetylation contributes to miR-31 silencing as treatment with the pan-HDAC inhibitor SAHA restored miR-31 expression. This is in line with recent observations that EZH2 interacts with HDAC3 to downregulate miR-29  and that EZH2 can interact with DNMT to control DNA methylation . Therefore, we propose a key role for Polycomb/EZH2, HDAC and DNMT for survival and metastasis of esophageal cancers.
Recent work argues that SOX4 is a master regulator of EMT and cell invasion through its binding to the EZH2 promoter, inducing EZH2 transcription . Our results indicate that miR-31 down-regulates SOX4 by binding to its 3′-UTR in EAC and ESCC cells. Consistent with previous reports [8,11], we found that miR-31 negatively regulates EZH2 as well. Furthermore, we observed a significant decrease in HDAC3 expression with miR-31 overexpression. We, therefore, propose that one underlying mechanism by which miR-31 suppresses tumor cell invasion is by directly targeting SOX4 and indirectly targeting EZH2 and HDAC3. Notably, overexpression of SOX4 is observed in a variety of human cancers . Loss of SOX4 led to a significant increase in miR-31 expression and strongly inhibited tumor cells proliferation, migration and invasion. Additional studies will aim to identify a subset of patients with concomitant high SOX4, EZH2 and HDAC3 and low miR-31 to demonstrate the mechanistic and clinical correlation between these pathways.
Multiple reports show that the polycomb PRC2/EZH2 interacts with HDAC to repress the transcription of CDH1 . Moreover, MYC was recently shown to interact with EZH2 and HDAC3 to repress miR-29 in lymphomas . However, it is not clear how PRC2/EZH2 is recruited to the promoter of specific genes. Similarly, SOX4 has been shown to interact with multiple transcription factors. Despite the fact that our computational analysis did not detect any SOX4 regulatory elements upstream of the miR-31 promoter, co-IP assays show that SOX4 interacts with EZH2 and HDAC3. Based on the ChIP experiments, we propose that EZH2 and HDAC3 bind to a similar region on the miR-31 promoter, confirming that histone methylation and histone deacetylation contribute to miR-31 repression. It is also possible that other molecules such as PRC2-interacting YY1 or the SOX4-interracting protein GATA4 mediate SOX4 recruitment to the DNA. Future studies are needed to confirm such possibilities.
Taken together, our data identify a molecular circuit where SOX4, EZH2 and HDAC3 cooperate to repress miR-31. This may lead to tumor progression and the metastasis of EAC and ESCC; SOX4, EZH2, HDAC3 and miR-31 emerge as potential therapeutic targets.
Materials and methods
The squamous cell carcinoma cell lines TE8, TE1, TE11, TE7, TE12, TE2, TE3, TE5 and TE9 were developed and characterized by Kuroki et al. . OE33 and FLO1 were established by Dr. David Beer (University of Michigan). HEK293T, CP-A and CP-B were purchased from American Type Culture Collection (ATCC). KYSE140, KYSE180, KYSE150, KYSE520, KYSE70 were developed and characterized by Shimada et al. . All cells were cultured at 37°C in 5% CO2.
Chemicals and antibodies
Suberoyl anilide hydroxamic acid (SAHA, a pan-HDAC inhibitor) and 3-Deazaneplanocin A (DZNep, a polycomb EZH2 subunit inhibitor) were purchased from Millipore, Temecula, CA. The DNA methylation inhibitor 5-AZA-Cytidine (AZA) was from Sigma, St. Louis, MO. Cells were treated with indicated concentration and incubated for 72 hours before harvesting. Agarose A/G plus is from Santa Cruz. The following antibodies were used as primary antibodies: SOX4 (Santa Cruz; 1:1000); EZH2 (BD; 1:4000); HDAC3 (Millipore; 1:4000); total histone H3 (Abcam; 1:5000); Tri-methylated histone H3 H3K27me3 (Millipore; 1:3000); SUZ12 (Santa Cruz; 1:1,000); EZH1 (Santa Cruz; 1:1,000); β-tubulin (Sigma-Aldrich; 1:5,000).
Oligonucleotides and plasmids
The pBABE-miR-31 plasmid (Plasmid#26088), pWPXL-SOX4 (plasmid#36984), pCMVHA-hEZH2 (plasmid#24230), Psicheck2 SOX4 full-length 3′UTR (Plasmid#26989) were purchased from Addgene (Cambridge, MA).
A 71 bp WT fragment of the SOX4 3′UTR (SOX4 WT OLIGO) was created by overlapping extension PCR and cloned between the XhoI and NotI site of the psicheck2 plasmid. Similarly, the mutant construct of SOX4 3′UTR (SOX4 Mutant OLIGO) which carried a substitution of four nucleotides within the core seed sequence of miR-31, was carried out using overlapping extension PCR and cloned between the XhoI and NotI site of the psicheck2 plasmid.
The two set of plasmids containing shRNA specific to SOX4 and EZH2 were purchased from OriGene (Rockville, MD).
Primers used for SOX4, forward: 5′-AGCGACAAGATCCCTTTCATTC-3′, reverse: 5′-CGTTGCCGGACTTCACCTT-3′; EZH2, forward: 5′-GTACACGGGGATAGAGAATGTGG-3′, reverse: 5′GGTGGGCGGCTTTCTTTATCA-3′, for EZH1, forward: 5′-ATGCGACTTCGACAACTTAAACG-3′, reverse: 5′-GGCTTCATTGACTGAACAGGTT-3′, HDAC3, forward: 5′-CCTGGCATTGACCCATAGCC-3′, reverse: 5′-CTCTTGGTGAAGCCTTGCATA-3′, GAPDH, forward: 5′-GCGACACCCACTCCTCCAC-3′, reverse: 5′-TCCACCACCCTGTTGCTGTAG-3′.
3′UTR Luciferase Reporter assays
HEK293, TE8 and FLO1 cells were plated in triplicate in a 24 well plate. One day after plating, cells were transfected with the dual Renilla and Firefly luciferase reporter plasmid (psiCHECK-2) containing the full length 3′UTR of SOX4 (plasmid# 26989), the short WT oligo or mutant oligo along with a pBABE-miR-31 or pBABE Empty Vector expressing plasmid using FuGene HD (Promega). 48 hr post-transfection, cells were lysed using 1X passive lysis buffer and lysates were analyzed using the Dual-Glo Luciferase Reporter Assay System (Promega) on the Synergy4 multi-mode microplate reader (BioTeK).
Quantitative real time PCR
Total RNA was isolated from cells with Qiazol reagent (Qiagen) and reverse transcribed into cDNA using the miRNeasy mini kit and miScript miRNA Reverse Transcription kit (Qiagen). qRT-PCR was performed according to the manufacturer’s instructions (Applied Biosystems).
WST-1 reagent (Roche) was used according to manufacturers’ protocols to assess cell viability.
Cell migration and invasion assays
Migration and Invasion assays were performed as previously described . Migrating or invading cells were then photographed using the Zeiss Axioskop Plus or eVos microscope and quantified with ImageJ software.
Colony formation assay
For clonogenic assay, 500 transfected cells were seeded in six-well plates and maintained in complete medium for 2 weeks. Colonies were fixed with ice-cold methanol and stained with crystal violet. Colonies were photographed and counted using GelCount (Oxford Optronix).
Western blot analysis
Western blot was performed as previously described .
Cells were collected and lysed in IP lysis buffer (150 mM NaCL, 50 mM Tris pH8, 1% Triton X-100, 1% NP-40) supplemented with protease and phosphatase inhibitors, incubated on ice for 20 min, and cleared by centrifugation at 13,200 rpm at 4 C for 20 min. Total protein lysate (500 μg) was immunoprecipitated with the agarose-immobilized antibody (6 μg of anti-SOX4, EZH2, HDAC3 or isotype control antibodies) and incubated overnight at 4°C. Immune complexes were eluted from the agarose beads and analyzed by SDS-PAGE followed by immunoblot analysis. For co-IP in 293 T, cells were transfected with plasmids using FuGene HD (Promega). Cells were collected 48 hours post-transfection and analyzed as described above.
Chromatin Immunoprecipitation (ChIP)
The ChIP assays were adapted and performed according to previous publications . Briefly, cells were fixed using 1% formaldehyde for 15 minutes and quenched using 125 mM glycine for 5 minutes at room temperature. After centrifugation the cell pellet was re-suspended in the cell lysis buffer (150 mM NaCL, 50 mM Tris pH8, 1% Triton X-100, 1% NP-40, 0.01% SDS, 1.2 mM EDTA pH 8.0, 1 mM PMSF). Protein-bound chromatin was fragmented by sonication. Equal volumes of chromatin were immunoprecipitated with anti-HDAC3, anti-SOX4, anti-EZH2, anti-trimethyl-Histone H3 Lys27 or normal IgG as a negative control (Millipore). Following extensive washing the immunoprecipitated DNA was treated with RNase (Qiagen) for 30 min at 37°C and proteinase K (Roche) for an hour at 45°C. The DNA was eluted using 100 mM NaHCO3 and 1% SDS and the crosslinks were reversed using 300 mM NaCl at 65°C for 16 hours. Immunoprecipitated DNA and whole cell extract DNA were purified by Qiaquick PCR purification kit (Qiagen). The purified DNA was amplified by real-time quantitative PCR with Qiagen QuantiTech SYBR Green PCR master mix and analyzed for enrichment. Real-time qPCR amplification was performed with Applied Biosystems StepOnePlus real time PCR system.
Primers used as described in Lin et al. : -1,000bp forward: CCGATGACCTAGCCAGAAGT, reverse: CCCCACCCTTCAACTCGTAG; -500 bp, forward: TATCCTCAACCCTCCGTGTC, reverse: CATACACCTGAAGGGGCAGT; +500 bp, forward: CAATTTTGGCCCAGGAGATA, reverse: TTTCCGGGGACCTCTAGTTT; +42,500 bp, forward: TGGCCTATTTGCTGTTCTAATGAC, reverse: GCAAGCCAACCCCAACA; +45,000 bp, forward: AATGGGCCCTGCATTCTCT, reverse: AAAACCCACACCCTCACCAC; +47,500 bp, forward: CATCTTCAAAAGCGGACACTCT, reverse: ACAATACATAGCAGGACAGGAAG; MYT, forward: AGGCACCTTCTGTTGGCCGA, reverse: AGGCAGCTGCCTCCCGTACA; GAPDH, forward: CGGCTACTAGCGGTTTTACG, reverse: AAGAAGATGCGGCTGACTGT.
Datasets made publicly available from GEO Datasets (http://www.ncbi.nlm.nih.gov/gds/). The collected information from each dataset was analyzed and visualized in Prism version 6.00 for Mac (GraphPad software, La Jolla, California).
Each experiment was repeated at least three times. Numerical data are presented as mean ± standard deviation. The differences between two groups were analyzed using a Student’s t-test (two-tailed) or two-way ANOVA. Differences were considered statistically significant at p < 0.05. All statistical analysis was performed on GraphPad Prism 6.0c software (La Jolla, CA).
Esophageal squamous cell carcinoma
Head and neck squamous cell carcinoma
Suberoyl anilide hydroxamic acid
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, et al. MicroRNA expression profiles classify human cancers. Nature. 2005;435:834–8.
Zhang B, Pan X, Cobb GP, Anderson TA. microRNAs as oncogenes and tumor suppressors. Dev Biol. 2007;302:1–12.
Tavazoie SF, Alarcon C, Oskarsson T, Padua D, Wang Q, Bos PD, et al. Endogenous human microRNAs that suppress breast cancer metastasis. Nature. 2008;451:147–52.
Sakai NS, Samia-Aly E, Barbera M, Fitzgerald RC. A review of the current understanding and clinical utility of miRNAs in esophageal cancer. Semin Cancer Biol. 2013;23:512–21.
Bandres E, Cubedo E, Agirre X, Malumbres R, Zarate R, Ramirez N, et al. Identification by Real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoral tissues. Mol Cancer. 2006;5:29.
Liu CJ, Tsai MM, Hung PS, Kao SY, Liu TY, Wu KJ, et al. miR-31 ablates expression of the HIF regulatory factor FIH to activate the HIF pathway in head and neck carcinoma. Cancer Res. 2010;70:1635–44.
Liu X, Sempere LF, Ouyang H, Memoli VA, Andrew AS, Luo Y, et al. MicroRNA-31 functions as an oncogenic microRNA in mouse and human lung cancer cells by repressing specific tumor suppressors. J Clin Invest. 2010;120:1298–309.
Yamagishi M, Nakano K, Miyake A, Yamochi T, Kagami Y, Tsutsumi A, et al. Polycomb-mediated loss of miR-31 activates NIK-dependent NF-kappaB pathway in adult T cell leukemia and other cancers. Cancer Cell. 2012;21:121–35.
Valastyan S, Reinhardt F, Benaich N, Calogrias D, Szasz AM, Wang ZC, et al. A pleiotropically acting microRNA, miR-31, inhibits breast cancer metastasis. Cell. 2009;137:1032–46.
Augoff K, McCue B, Plow EF, Sossey-Alaoui K. miR-31 and its host gene lncRNA LOC554202 are regulated by promoter hypermethylation in triple-negative breast cancer. Mol Cancer. 2012;11:5.
Asangani IA, Harms PW, Dodson L, Pandhi M, Kunju LP, Maher CA, et al. Genetic and epigenetic loss of microRNA-31 leads to feed-forward expression of EZH2 in melanoma. Oncotarget. 2012;3:1011–25.
Creighton CJ, Fountain MD, Yu Z, Nagaraja AK, Zhu H, Khan M, et al. Molecular profiling uncovers a p53-associated role for microRNA-31 in inhibiting the proliferation of serous ovarian carcinomas and other cancers. Cancer Res. 2010;70:1906–15.
Fuse M, Kojima S, Enokida H, Chiyomaru T, Yoshino H, Nohata N, et al. Tumor suppressive microRNAs (miR-222 and miR-31) regulate molecular pathways based on microRNA expression signature in prostate cancer. J Hum Genet. 2012;57:691–9.
Chen Z, Saad R, Jia P, Peng D, Zhu S, Washington MK, et al. Gastric adenocarcinoma has a unique microRNA signature not present in esophageal adenocarcinoma. Cancer. 2013;119:1985–93.
Saad R, Chen Z, Zhu S, Jia P, Zhao Z, Washington MK, et al. Deciphering the unique microRNA signature in human esophageal adenocarcinoma. PLoS One. 2013;8:e64463.
Leidner RS, Ravi L, Leahy P, Chen Y, Bednarchik B, Streppel M, et al. The microRNAs, MiR-31 and MiR-375, as candidate markers in Barrett’s esophageal carcinogenesis. Genes Chromosomes Cancer. 2012;51:473–9.
Hotte GJ, Linam-Lennon N, Reynolds JV, Maher SG. Radiation sensitivity of esophageal adenocarcinoma: the contribution of the RNA-binding protein RNPC1 and p21-mediated cell cycle arrest to radioresistance. Radiation Res. 2012;177:272–9.
Zhang T, Wang Q, Zhao D, Cui Y, Cao B, Guo L, et al. The oncogenetic role of microRNA-31 as a potential biomarker in oesophageal squamous cell carcinoma. Clin Sci (Lond). 2011;121:437–47.
Lin RJ, Xiao DW, Liao LD, Chen T, Xie ZF, Huang WZ, et al. MiR-142-3p as a potential prognostic biomarker for esophageal squamous cell carcinoma. J Surg Oncol. 2012;105:175–82.
Bowles J, Schepers G, Koopman P. Phylogeny of the SOX family of developmental transcription factors based on sequence and structural indicators. Dev Biol. 2000;227:239–55.
Vervoort SJ, van Boxtel R, Coffer PJ. The role of SRY-related HMG box transcription factor 4 (SOX4) in tumorigenesis and metastasis: friend or foe? Oncogene. 2013;32:3397–409.
Dy P, Penzo-Mendez A, Wang H, Pedraza CE, Macklin WB, Lefebvre V. The three SoxC proteins–Sox4, Sox11 and Sox12–exhibit overlapping expression patterns and molecular properties. Nucleic Acids Res. 2008;36:3101–17.
Schilham MW, Oosterwegel MA, Moerer P, Ya J, de Boer PA, van de Wetering M, et al. Defects in cardiac outflow tract formation and pro-B-lymphocyte expansion in mice lacking Sox-4. Nature. 1996;380:711–4.
Bhattaram P, Penzo-Mendez A, Sock E, Colmenares C, Kaneko KJ, Vassilev A, et al. Organogenesis relies on SoxC transcription factors for the survival of neural and mesenchymal progenitors. Nat Commun. 2010;1:9.
Deneault E, Cellot S, Faubert A, Laverdure JP, Frechette M, Chagraoui J, et al. A functional screen to identify novel effectors of hematopoietic stem cell activity. Cell. 2009;137:369–79.
Pece S, Tosoni D, Confalonieri S, Mazzarol G, Vecchi M, Ronzoni S, et al. Biological and molecular heterogeneity of breast cancers correlates with their cancer stem cell content. Cell. 2010;140:62–73.
Kobielak K, Stokes N, de la Cruz J, Polak L, Fuchs E. Loss of a quiescent niche but not follicle stem cells in the absence of bone morphogenetic protein signaling. Proc Natl Acad Sci U S A. 2007;104:10063–8.
Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, et al. Large-scale meta-analysis of cancer microarray data identifies common transcriptional profiles of neoplastic transformation and progression. Proc Natl Acad Sci U S A. 2004;101:9309–14.
Scharer CD, McCabe CD, Ali-Seyed M, Berger MF, Bulyk ML, Moreno CS. Genome-wide promoter analysis of the SOX4 transcriptional network in prostate cancer cells. Cancer Res. 2009;69:709–17.
Tiwari N, Tiwari VK, Waldmeier L, Balwierz PJ, Arnold P, Pachkov M, et al. Sox4 is a master regulator of epithelial-mesenchymal transition by controlling Ezh2 expression and epigenetic reprogramming. Cancer Cell. 2013;23:768–83.
Sparmann A, van Lohuizen M. Polycomb silencers control cell fate, development and cancer. Nat Rev Cancer. 2006;6:846–56.
Bracken AP, Dietrich N, Pasini D, Hansen KH, Helin K. Genome-wide mapping of Polycomb target genes unravels their roles in cell fate transitions. Genes Dev. 2006;20:1123–36.
He A, Shen X, Ma Q, Cao J, von Gise A, Zhou P, et al. PRC2 directly methylates GATA4 and represses its transcriptional activity. Genes Dev. 2012;26:37–42.
Lee JM, Lee JS, Kim H, Kim K, Park H, Kim JY, et al. EZH2 generates a methyl degron that is recognized by the DCAF1/DDB1/CUL4 E3 ubiquitin ligase complex. Mol Cell. 2012;48:572–86.
Xu K, Wu ZJ, Groner AC, He HH, Cai C, Lis RT, et al. EZH2 oncogenic activity in castration-resistant prostate cancer cells is Polycomb-independent. Science. 2012;338:1465–9.
Kim E, Kim M, Woo DH, Shin Y, Shin J, Chang N, et al. Phosphorylation of EZH2 activates STAT3 signaling via STAT3 methylation and promotes tumorigenicity of glioblastoma stem-like cells. Cancer Cell. 2013;23:839–52.
Zhang X, Zhao X, Fiskus W, Lin J, Lwin T, Rao R, et al. Coordinated silencing of MYC-mediated miR-29 by HDAC3 and EZH2 as a therapeutic target of histone modification in aggressive B-Cell lymphomas. Cancer Cell. 2012;22:506–23.
Cao Q, Yu J, Dhanasekaran SM, Kim JH, Mani RS, Tomlins SA, et al. Repression of E-cadherin by the polycomb group protein EZH2 in cancer. Oncogene. 2008;27:7274–84.
Kleer CG, Cao Q, Varambally S, Shen R, Ota I, Tomlins SA, et al. EZH2 is a marker of aggressive breast cancer and promotes neoplastic transformation of breast epithelial cells. Proc Natl Acad Sci U S A. 2003;100:11606–11.
Varambally S, Cao Q, Mani RS, Shankar S, Wang X, Ateeq B, et al. Genomic loss of microRNA-101 leads to overexpression of histone methyltransferase EZH2 in cancer. Science. 2008;322:1695–9.
Bracken AP, Kleine-Kohlbrecher D, Dietrich N, Pasini D, Gargiulo G, Beekman C, et al. The Polycomb group proteins bind throughout the INK4A-ARF locus and are disassociated in senescent cells. Genes Dev. 2007;21:525–30.
Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42:D68–73.
Sharma S, Kelly TK, Jones PA. Epigenetics in cancer. Carcinogenesis. 2010;31:27–36.
Lin PC, Chui YL, Banerjee S, et al. Epigenetic Repression of miR-31 Disrupts Androgen Receptor Homeostasis and Contributes to Prostate Cancer Progression. Cancer Res. 2013;73:1232–44.
Taylor MA, Sossey-Alaoui K, Thompson CL, Danielpour D, Schiemann WP. TGF-beta upregulates miR-181a expression to promote breast cancer metastasis. J Clin Invest. 2013;123:150–63.
Lynam-Lennon N, Reynolds JV, Marignol L, Sheils OM, Pidgeon GP, Maher SG. MicroRNA-31 modulates tumour sensitivity to radiation in oesophageal adenocarcinoma. J Mol Med (Berl). 2012;90:1449–58.
Park SM, Gaur AB, Lengyel E, Peter ME. The miR-200 family determines the epithelial phenotype of cancer cells by targeting the E-cadherin repressors ZEB1 and ZEB2. Genes Dev. 2008;22:894–907.
Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, Farshid G, et al. The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nature Cell Biol. 2008;10:593–601.
Vire E, Brenner C, Deplus R, Blanchon L, Fraga M, Didelot C, et al. The Polycomb group protein EZH2 directly controls DNA methylation. Nature. 2006;439:871–4.
Nishihira T, Hashimoto Y, Katayama M, Mori S, Kuroki T. Molecular and cellular features of esophageal cancer cells. J Cancer Res Clin Oncol. 1993;119:441–9.
Shimada Y, Imamura M, Wagata T, Yamaguchi N, Tobe T. Characterization of 21 newly established esophageal cancer cell lines. Cancer. 1992;69:277–84.
Andl CD, McCowan KM, Allison GL, Rustgi AK. Cathepsin B is the driving force of esophageal cell invasion in a fibroblast-dependent manner. Neoplasia. 2010;12:485–98.
Tantin D, Voth WP, Shakya A. Efficient chromatin immunoprecipitation using limiting amounts of biomass. Journal of visualized experiments: JoVE. 2013:e50064.
This work was supported by the National Institute of Health (DK94900, DK091491). The use of Research Cores is supported by the Vanderbilt Ingram Cancer Center (P30 CA68485) and the Vanderbilt Digestive Disease Research Center (P30 DK058404). Further support is from the National Center for Research Resources, Grant UL1 RR024975-01, and is now at the National Center for Advancing Translational Sciences, Grant 2 UL1 TR000445-06. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH. We would like to thank Drs. Jeff Franklin and Christine Eischen for critical discussion of the manuscript.
The authors declare that they have no competing interests.
RBK carried out the experimental approach with assistance from KJT, STZ, and CJT. HAL performed the statistical analysis. RBK, TA and CA conceived of the study, and participated in its design and coordination and wrote the manuscript with editing assistance from KJT and CJT. All authors read and approved the final manuscript.