- Open Access
HOXC8: a predictive glioma biomarker that induces epithelia-mesenchymal transition
Chinese Neurosurgical Journal volume 4, Article number: 24 (2018)
The transcription factor Homeobox C8 (HOXC8) is overexpressed and regulates many important genes involved in the proliferation and invasion of many malignant tumors. However, the function of HOXC8 in gliomas remains unclear.
Based on the Chinese Glioma Genome Atlas (CGGA) set, HOXC8 expression is negatively correlated with overall survival (OS). Small interfering RNA (si-HOXC8) was used to downregulate the mRNA and protein expression levels of HOXC8 to assess glioma cell proliferation, migration and invasion.
Patients with higher HOXC8 levels showed poorer prognosis. DAVID analysis results indicated that HOXC8 was related to cell cycle, cell adhesion and immune response. In U251 and LN229 glioma cells treated with small interfering RNA for HOXC8 (si-HOXC8) for gene knockdown, significantly lower cell capacity of growth, migration and invasion was observed. Moreover, HOXC8 knockdown could reduce the protein expression of classical epithelial mesenchymal transition (EMT) related markers.
HOXC8 may play an important role in glioma proliferation, migration and invasion. These findings indicated that HOXC8 may constitute a novel target for glioma treatment.
Gliomas account for most of malignant tumors of the central nervous system (CNS). Despite the recent progress in surgical resection, radiotherapy and chemotherapy, prognosis for glioma remains poor [1, 2]. Gliomas possessed strong proliferation, migration and invasion abilities, causing high mortality. Therefore, multiple studies have identified key genes that could help achieve a breakthrough in gliomas treatment [3, 4]. In this study, Homeobox C8 (HOXC8), a transcript factor, was selected for in-depth assessment.
HOXC8, located in a cluster on chromosome 12, belongs to the 39-member HOX family of proteins . Under normal circumstances, HOXC8 plays an indispensable role in embryonic morphogenesis and system morphogenesis [6, 7]. However, in various tumors, including ovarian cancer, hepatocellular carcinoma, and breast cancer, HOXC8 induces cell proliferation, migration and invasion, and is inversely correlated with overall survival (OS) through different mechanisms [8,9,10]. Although HOXC8 is important in tumorigenesis, its mechanism in glioma remains unclear.
Epithelial mesenchymal transition (EMT), a cellular alteration that confers a more invasive and drug-resistant phenotype, was initially observed in embryonic development [11, 12]. EMT was a complex process reflected by epithelial cells acquiring the mesenchymal phenotype and motility through a cascade of biological events . In gliomas, the mechanism by which HOXC8 affected EMT is unclear.
In summary, the aim of our study was to assess the oncogenic function of HOXC8. HOXC8 knockdown could inhibit glioma cell proliferation, migration and invasion. Therefore, HOXC8 could provide novel insights for glioma treatment.
Dataset and glioma samples
We obtained mRNA samples and follow-up information from the Chinese Glioma Genome Atlas (CGGA). OS was defined as the time from surgery to death. Glioma samples were collected from surgical resection, snap-frozen in liquid nitrogen and stored at − 80 °C until RNA extraction . This study was approved by the Ethics Committee of Beijing Tiantan Hospital, and written informed consent was obtained from all the patients.
Cell culture and siRNA transfection
Two glioma cell lines were cultured in Dulbecco’s Modified Eagle’s Medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and antibiotics (100 U/ml penicillin and 100 μg/mL streptomycin), in a humidified atmosphere with 5% CO2 at 37 °C . The siRNA segmentand a non-specific control siRNA sequence (NC) were purchased from Genepharma (Shanghai, China), and transfected into LN229 and U251 cells using Lipofectamine™ 2000 (Invitrogen, Shanghai, China), according to the manufacturer’s instructions. After 48 h of incubation, samples were collected for RT-PCR and Western blot assay .
Real-time polymerase chain reaction (RT-PCR) and western blot (WB)
RT-PCR and Western blot were performed as previously described [4, 17].The following primers were used for RT-PCR: HOXC8, Forward 5’-ACCGGCCTATTACGACTGC-3′ and Reverse 5′-TGCTGGTAGCCTGAGTTGGA-3′; GAPDH (Forward 5′- GGAGCGAGATCCCTCCAAAAT -3′ and Reverse 5′- GGCTGTTGTCATACTTCTCATGG-3′ was used as an internal control, and fold changes were calculated by the 2 − ΔΔCt method. For Western blotting, densities of specific protein bands were quantified after normalization to Tubulin or Actin level in the same sample.
Clone formation assay
A total of 1000 cells were seeded per well in a 6-well plate in DMEM and 10% FBS, and cultured for 2 weeks before Crystal violet stain. Then, clone numbers were recorded in siRNA and control groups, respectively.
Transwell migration assay
Transwell plates were used to study cell migration and invasion. A total of 1 × 105 LN229 or U251 cells were seeded in each well in DMEM only. Control medium was added to the lower chamber of the transwell plates. After 36 h (invasion assay) or 12 h (migration assay) of incubation, the cells in upper chambers were removed carefully before Crystal violet staining. Analysis was carried out by microscopy, counting cells in siRNA and control groups, respectively .
We used Kaplan-Meier analysis (log-rank test) to evaluate predictive value of HOXC8 expression for OS in different groups. Univariate and multivariate Cox regression analyses were used to determine that HOXC8 was an independent factor to predict OS. Pearson relation analysis was performed in CGGA array set using R language package. All values were considered statistically significant at p < 0.05.
HOXC8 was associated with overall survival in gliomas
We obtained the HOXC8 mRNA expression and prognosis information from the CGGA mRNA array and TCGA RNA-seq datasets, and assessed the association of HOXC8 expression with overall survival (OS). As shown in Fig. 1a and b, the level of HOXC8 mRNA expression was negatively related with OS in all grade gliomas and grade IV gliomas (p < 0.05). Moreover, we further validated the results in TCGA RNA-seq set (Fig. 1c and d, p < 0.05). Moreover, in Fig. 1e and f, in CGGA array set, progression-free survival time of was shorter in higher HOXC8 expression group (p < 0.05). To further validation our results, we performed a univariate Cox regression analysis, the results showed that HOXC8 was showed a risk factor for glioma patients (p < 0.0001, HR = 1.545, Table 1). What is more, age at diagnosis, IDH1 status, gliomas grade, TCGA subtype was significantly related with survival time of glioma patients (p < 0.0001, Table 1). Then, the results from multivariate Cox regression analysis including HOXC8 expression, age at diagnosis, IDH1 status, gliomas grade and TCGA subtype showed that HOXC8 was an independent prognostic factor for glioma patients (p < 0.0001, HR = 1.247, Table 1).
In conclusion, we hypothesized that HOXC8 could serve as a novel biomarker for predicting glioma OS.
Gene functional analysis
The function of HOXC8 in glioma remains unclear. We performed Pearson relation analysis in CGGA array set using R language package. To explore differential biological function of HOXC8, we selected top 1000 genes positively or negatively related with HOXC8 for DAVID analysis respectively. In Fig. 2a and b, Gene Ontology (GO) and KEGG pathways analysis showed that cell cycle, cell adhesion, immune response and inflammatory response were high expression in high-HOXC8 expression group. Meanwhile, taken together, gliomas in high-HOXC8 group had a greater capacity for proliferation, migration and invasion than these in low-HOXC8 group.
LN229 and U251 glioma cell lines were used to investigate biological processes
Then in Fig. 3a and b, seven (six glioma cell lines and one normal human astrocyte) frequently-studied cell lines (LN229, U87, U118, U251, SNB19, H4, and HA) were selected for HOXC8 mRNA and protein expression assessment. Interestingly, the six glioma cell lines showed higher HOXC8 expression levels compared with HA cells (human astrocytes). More importantly, LN229 and U251 cells had highest endogenous mRNA and protein levels of HOXC8. Based on the above results, LN229 and U251 were selected for subsequent experiments.
HOXC8 could significantly promote glioma proliferation
To assess the effect of HOXC8 on cell proliferation, HOXC8 was silenced in LN229 and U251 cells by siRNA; 48 h after transfection, siRNA successfully knocked down HOXC8 at the mRNA levels (Fig. 4a and b). Cell viability was measured by MTT and clonogenic assays, as shown in Fig. 5a-c, antagonism of HOXC8 expression suppressed the cell growth and the formation of cell clones (p < 0.05).
Downregulation of HOXC8 inhibits cell migration and invasion by blocking the EMT pathway
To assess whether the migration and invasion capacity of glioma cells could be inhibited by si-HOXC8, we performed migration and Matrigel invasion assays using LN229 and U251 cells. Crystal violet staining showed markedly less migration and invasion glioma cells after HOXC8 silencing compared with negative-control values (p < 0.01; Fig. 6a and b). In addition, EMT was a novel mechanism involved in metastasis [19, 20]. The protein expression levels of mesenchymal markers (N-Cadherin, Slug, Vimentin and MMP9) were reduced in si-HOXC8 LN229 and U251 cells compared with the control group, as assessed by Western blot (Fig. 6c and d). Taken together, si-HOXC8 suppressed cell migration and invasion by reversing EMT.
Massive evidence supports that HOXC8 plays an important role in malignant tumor formation and progression [21, 22]. As a transcription factor, HOXC8 is able to regulate and coordinate multiplevital genes (e.g. Mgl1, Embigin, Meis1, and Fyn) involved in cancer development and progression [8, 22, 23]. Previous findings demonstrated that HOXC8 depletion by small interfering RNA suppresses epithelial ovarian cancer proliferation and migration, and induces apoptosis by increasing ZAC1 expression . HOXC8, overexpressed in hepatocellular carcinoma (HCC) compared with adjacent non-tumor tissues, is associated with poor prognosis. Si-HOXC8 blocks G1-S phase transition, inhibits cell proliferation and renders cells more sensitive to oxaliplatin . Therefore, HOXC8 as a potential oncogenic driver, plays an important role in cancer cell proliferation, migration and invasion, by overregulating many oncogenes and downregulating tumor suppressor genes. In this study, we demonstrated that si-HOXC8 blunted glioma cell proliferation, migration and invasion by reversing EMT.
EMT is known to be implicated in cancer progression, metastasis and drug resistance . EMT reversibly enables polarized epithelial cells to lose their epithelial characteristics and to acquire mesenchymal properties . The availability of antibodies targeting epithelial (E-cadherin, β-catenin, and Claudin-1) and mesenchymal (N-cadherin, Slug, Vimentin and MMP9) markers makes it convenient to assess EMT . Multiple genes and pathways were reported to be involved in EMT [26,27,28]. However, the mechanisms by which HOXC8 affects EMT remains unclear. This study found that the HOXC8 was involved in EMT. Moreover, si-HOXC8 decreased the expression of classical mesenchymal markers (N-cadherin, Slug, Vimentin and MMP9). Hence, the present findings provide additional insights regarding the implication of HOXC8 and target genes in EMT as well as malignant tumor progression.
In summary, this study demonstrated that higher HOXC8 expression resulted in shorter OS. In agreement, HOXC8 knockdown inhibited glioma cell proliferation. Meanwhile, reduction of HOXC8 could also affected migration and invasion through blocking EMT pathway. The present findings suggested that HOXC8 should be considered a novel biomarker and target for glioma treatment.
Norden AD, Drappatz J, Wen PY. Antiangiogenic therapies for high-grade glioma. Nat Rev Neurol. 2009;5:610–20.
Zhou X, Ren Y, Zhang J, et al. HOTAIR is a therapeutic target in glioblastoma. Oncotarget. 2015;6:8353–65.
Chen PH, Chang CK, Shih CM, et al. The miR-204-3p-targeted IGFBP2 pathway is involved in xanthohumol-induced glioma cell apoptotic death. Neuropharmacology. 2016;110:362–75.
You G, Feng L, Yan W, et al. BCL2A1 is a potential biomarker for postoperative seizure control in patients with low-grade gliomas. CNS Neurosci Ther. 2013;19:882–8.
Gehring WJ, Hiromi Y. Homeotic genes and the homeobox. Annu Rev Genet. 1986;20:147–73.
Gao Y, Hyttel P, Hall VJ. Regulation of H3K27me3 and H3K4me3 during early porcine embryonic development. Mol Reprod Dev. 2010;77:540–9.
Kalyani R, Lee JY, Min H, Yoon H, Kim MH. Genes frequently Coexpressed with Hoxc8 provide insight into the discovery of target genes. Mol Cells. 2016;39:395–402.
Chao F, Zhang J, Zhang Y, et al. Embigin, regulated by HOXC8, plays a suppressive role in breast tumorigenesis. Oncotarget. 2015;6:23496–509.
Xu P, Zhang X, Ni W, et al. Upregulated HOXC8 expression is associated with poor prognosis and Oxaliplatin resistance in hepatocellular carcinoma. Dig Dis Sci. 2015;60:3351–63.
Zhou W, Ke SQ, Huang Z, et al. Periostin secreted by glioblastoma stem cells recruits M2 tumour-associated macrophages and promotes malignant growth. Nat Cell Biol. 2015;17:170–82.
Zeisberg M, Neilson EG. Biomarkers for epithelial-mesenchymal transitions. J Clin Invest. 2009;119:1429–37.
Zhang L, Zhang W, Li Y, et al. SHP-2-upregulated ZEB1 is important for PDGFRalpha-driven glioma epithelial-mesenchymal transition and invasion in mice and humans. Oncogene. 2016;35:5641–52.
Felipe Lima J, Nofech-Mozes S, Bayani J, Bartlett JM. EMT in breast carcinoma-a review. J Clin Med. 2016;5
Yang P, Liang T, Zhang C, et al. Clinicopathological factors predictive of postoperative seizures in patients with gliomas. Seizure. 2016;35:93–9.
Kaur H, Ali SZ, Huey L, et al. The transcriptional modulator HMGA2 promotes stemness and tumorigenicity in glioblastoma. Cancer Lett. 2016;377:55–64.
Wang X, Shi W, Shi H, et al. TRIM11 overexpression promotes proliferation, migration and invasion of lung cancer cells. J Exp Clin Cancer Res: CR. 2016;35:100.
Li MY, Yang P, Liu YW, et al. Low c-met expression levels are prognostic for and predict the benefits of temozolomide chemotherapy in malignant gliomas. Sci Rep. 2016;6:21141.
Zhao T, Yang H, Tian Y, et al. SOX7 is associated with the suppression of human glioma by HMG-box dependent regulation of Wnt/beta-catenin signaling. Cancer Lett. 2016;375:100–7.
Li J, Yuan J, Yuan X, et al. MicroRNA-200b inhibits the growth and metastasis of glioma cells via targeting ZEB2. Int J Oncol. 2016;48:541–50.
Thiyagarajan V, Tsai MJ, Weng CF. Antroquinonol targets FAK-signaling pathway suppressed cell migration, invasion, and tumor growth of C6 glioma. PLoS One. 2015;10:e0141285.
Du YB, Dong B, Shen LY, et al. The survival predictive significance of HOXC6 and HOXC8 in esophageal squamous cell carcinoma. J Surg Res. 2014;188:442–50.
Ruthala K, Gadi J, Lee JY, Yoon H, Chung HJ, Kim MH. Hoxc8 downregulates Mgl1 tumor suppressor gene expression and reduces its concomitant function on cell adhesion. Molecules Cells. 2011;32:273–9.
Lei H, Wang H, Juan AH, Ruddle FH. The identification of Hoxc8 target genes. Proc Natl Acad Sci U S A. 2005;102:2420–4.
Lu S, Liu R, Su M, et al. Overexpression of HOXC8 is associated with poor prognosis in epithelial ovarian Cancer. Reprod Sci. 2016;23:944–54.
Du B, Shim JS. Targeting epithelial-mesenchymal transition (EMT) to overcome drug resistance in Cancer. Molecules. 2016;21
Wang H, Tao T, Yan W, et al. Upregulation of miR-181s reverses mesenchymal transition by targeting KPNA4 in glioblastoma. Sci Rep. 2015;5:13072.
Ding Y, Li X, Hong D, Jiang L, He Y, Fang H. Silence of MACC1 decreases cell migration and invasion in human malignant melanoma through inhibiting the EMT. Biosci Trends. 2016;10:258–64.
Tang H, Fan D, Lei CT, et al. MAD2B promotes tubular epithelial-to-mesenchymal transition and renal tubulointerstitial fibrosis via Skp2. J Mol Med. 2016;94:1297–307.
Thanks for Chinese Glioma Genome Atlas dataset and The Cancer Genome Atlas dataset for data share.
This study was funded by National Natural Science Fund (Grant No. 81301112) and Yu Miao Plan of Beijing Ditan Hospital (Grant No. DTYM201603) and Beijing Medical Management Bureau Cultivation Plan (Grant No. PX2018079).
Availability of data and materials
Ethics approval and consent to participate
This study was approved by the Ethics Committee of Beijing Tiantan Hospital, and written informed consent was obtained from all the patients.
Consent for publication
The authors declare that they have no competing interest.