Research Paper Volume 15, Issue 12 pp 5412—5425
Tumor-suppressive function and mechanism of miR-873-5p in glioblastoma: evidence based on bioinformatics analysis and experimental validation
- 1 Department of Neurosurgery, General Hospital of the Northern Theater Command of Chinese People’s Liberation Army, Shenyang 110000, China
- 2 Department of Neurosurgery, Jinqiu Hospital of Liaoning Province, Shenyang 110000, China
Received: February 27, 2023 Accepted: May 16, 2023 Published: June 28, 2023
https://doi.org/10.18632/aging.204800How to Cite
Copyright: © 2023 Zhang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Abstract
This study aims to clarify the mechanistic actions of microRNA-873-5p (miR-873-5p) on glioblastoma (GBM) progression. The most differentially expressed miRNAs were retrieved from the GEO database. It was established that miR-873-5p was downregulated in GBM tissues and cells. Based on in silico prediction and experimental data, HMOX1 was demonstrated to be a target gene of miR-873-5p. Further, miR-873-5p was then ectopically expressed in GBM cells to examine its effect on the malignant behaviors of GBM cells. Overexpression of miR-873-5p inhibited GBM cell proliferation and invasion by targeting HMOX1. HMOX1 promoted SPOP expression by increasing HIF1α expression, thus stimulating GBM cell malignant phenotypes. miR-873-5p suppressed the malignant phenotypes of GBM cells and tumorigenesis in vitro and in vivo by inhibiting the HMOX1/HIF1α/SPOP signaling axis. This study uncovers a novel miR-873-5p/HMOX1/HIF1α/SPOP axis in GBM, providing new insights into GBM progression and therapeutic targets for GBM treatment.
Introduction
Glioblastoma (GBM) represents the most aggressive form of brain tumor [1]. Partly owing to the plastic and heterogenous nature, GBM continues to exhibit very high mortality and morbidity rates [2, 3]. Therefore, gaining a comprehensive understanding of the molecular mechanisms of GBM pathogenesis is urgently needed to develop novel therapeutic strategies for GBM treatment.
Of note, microRNA (miR)-873-5p is a main member of the miR-873 family [4]. A recent study has proved the inhibiting effect of miR-873-5p in glioma cell proliferation and invasion and the promoting effect on cell apoptosis [5]. Meanwhile, miR-873-5p can repress GBM cell growth and migration, but inhibition of miR-873-5p shows opposite effects [6]. miR-873-5p has been demonstrated to target heme oxygenase 1 (HMOX1) and regulate its expression in PC12 cells [7]. HMOX1 is an inducible enzyme capable of catalyzing the degradation of heme [8]. HMOX1 has demonstrated a correlation with the stemness and metastasis of GBM, substantiating the histopathological hallmark of GBM [9]. In human glioma tumors, HMOX1 is robustly induced compared with normal brain tissues and relates to the molecular mechanism driving GBM cell proliferation survival and metastasis [10].
Furthermore, hypoxia-induced factor 1α (HIF1α) transcription is lowered in bone marrow-derived macrophages with HMOX1 knockdown [11]. HIF1α is a subunit of HIF1 activated by hypoxia and a key driver of tumor progression in GBM patients [12, 13]. Elevation of HIF1α promotes the expression of angiogenic factors and elicits vascular abnormality and GBM progression [14]. SPOP encodes an E3 ubiquitin ligase component, and in clear cell renal cell carcinoma, it is a transcriptional target of HIFs [15]. SPOP exerts essential roles in physiological and pathological processes [16]. Notably, a previous study has suggested SPOP as a tumor suppressor in glioma to inhibit cell viability, migration, and invasion in vitro [17]. Therefore, it is reasonable to hypothesize that miR-873-5p may serve as a tumor suppressor in GBM with the involvement of HMOX1, HIF1α and SPOP. Consequently, this study employed bioinformatics analysis combined with experimental validation with the aim to validate the role of miR-873-5p in GBM progression and the potential mechanism.
Results
miR-873-5p is downregulated in GBM tissues and cell lines, which targets HMOX1
First, we retrieved a GBM-related miRNA expression dataset, GSE103229, from the GEO database and conducted the differential analysis with |logFC| > 1 and p < 0.05 as the threshold. miR-873-5p was found to be the most significantly differentially expressed miRNA, with a significant down-regulation in brain tissue samples from GBM patients compared to normal brain tissue samples (Figure 1A).
Figure 1. miR-873-5p, downregulated in GBM tissues or cells, targets HMOX1 and reduces its expression. (A) A volcano plot of the differentially expressed genes in GBM and normal samples from the GSE103229 dataset. Red dots indicate significantly up-regulated genes, and blue dots indicate significantly downregulated genes. The top hsa-miR-873-5p with the most significant differential expression (the lowest p value) was labeled. (B) Venn diagram of the downstream genes of miR-873-5p predicted using the miRDB, TargetScan, and mirDIP databases with the significantly up-regulated genes in GBM samples from the GSE104267 dataset. (C) HMOX1 expression in normal brain tissue samples and GBM samples analyzed by the GEPIA database. The red box represents the GBM tissue samples (n = 163), and the gray box represents the normal brain tissue samples (n = 207). (D) Binding sites of miR-873-5p in HMOX1 mRNA predicted by the miRDB database. (E) Expression of miR-873-5p and HMOX1 in tumor tissues of GBM patients (n = 20) and normal brain tissues (n = 20) measured by RT-qPCR. (F) Expression of miR-873-5p and HMOX1 in human normal brain glial cell line HEB and GBM cell lines T98G, U251, LN-229, and A172 determined by RT-qPCR. (G) Binding of miR-873-5p to HMOX1 validated by dual-luciferase reporter assay. (H) Expression of miR-873-5p in T98G cells transfected with miR-873-5p mimic or miR-873-5p inhibitor determined by RT-qPCR. (I) Western blot of HMOX1 protein in T98G cells transfected with miR-873-5p mimic or miR-873-5p inhibitor. *p < 0.05, compared to normal brain tissues, HEB cells or NC mimic-transfected T98G cells. #p < 0.05, compared to NC inhibitor-transfected T98G cells.
miRDB, TargetScan, and mirDIP databases predicted 123, 72 and 196 downstream target genes of miR-873-5p, respectively. Meanwhile, 816 significantly up-regulated genes were obtained following the differential analysis of the GBM-related gene expression dataset GSE104267, and then intersected with the predicted downstream genes, with one candidate gene determined, namely HMOX1 (Figure 1B). According to the analysis results of the GEPIA database, HMOX1 expression was amplified in GBM samples relative to normal brain tissue samples (Figure 1C). In addition, miR-873-5p binding sites in HMOX1 mRNA were found following analysis using the miRDB database (Figure 1D).
Based on RT-qPCR results, miR-873-5p expression was lowered while HMOX1 expression was elevated in tumor tissues of GBM patients compared with normal tissues (Figure 1E). Furthermore, similar results were obtained in the T98G, U251, LN-229, and A172 cells compared to HEB cells. Of note, the T98G cell line showed the lowest miR-873-5p expression and the highest HMOX1 expression (Figure 1F), which was thus selected for subsequent experiments.
Dual-luciferase reporter assay results demonstrated a decline of HMOX1-Wt luciferase activity upon miR-873-5p mimic transfection, while no changes were noted in HMOX1-Mut luciferase activity (Figure 1G), indicating that miR-873-5p could bind to HMOX1.
To test miR-873-5p regulation on HMOX1 expression, we transfected T98G cells with miR-873-5p mimic and miR-873-5p inhibitor. RT-qPCR determined the transfection efficiency (Figure 1H). In addition, western blot results showed reduced HMOX1 expression in T98G cells with miR-873-5p mimic, while opposite results were found upon miR-873-5p inhibitor (Figure 1I).
miR-873-5p inhibits GBM cell malignant characteristics via targeting HMOX1
Then, we shifted to testing the effect of miR-873-5p and HMOX1 on GBM cell functions. HMOX1 protein expression was decreased in T98G cells transfected with miR-873-5p mimic + oe-NC while it was increased in response to oe-HMOX1 alone or combined with miR-873-5p mimic (Figure 2A).
Figure 2. miR-873-5p impairs GBM cell proliferation and invasion by targeting HMOX1. T98G cells were transfected with miR-873-5p mimic, oe-HMOX1 or in combination. (A) Western blot of HMOX1 protein in T98G cells. (B) Proliferation of T98G cells measured by CCK-8 assay. (C) Invasion of T98G cells measured by Transwell assay. (D) Western blot of the CCND1, E-cadherin, MYC, MMP-7, and N-cadherin proteins in T98G cells. *p < 0.05, compared to NC-mimic + oe-NC-transfected T98G cells. #p < 0.05, compared to miR-873-5p mimic + oe-NC-transfected T98G cells. The cell experiment was run in triplicate independently.
Moreover, a reduction of cell proliferation and invasion was noted upon miR-873-5p overexpression, the effect of which was undermined by HMOX1 overexpression (Figure 2B, 2C). Meanwhile, the protein expression of CCND1, MYC, MMP-7, and N-cadherin was reduced while that of E-cadherin was increased in T98G cells with miR-873-5p mimic + oe-NC; conversely, oe-HMOX1 alone or combined with miR-873-5p mimic caused opposite results (Figure 2D).
HMOX1 up-regulates SPOP expression by augmenting HIF1α expression
Analysis of the GEPIA database indicated that HIF1α expression was augmented in the GBM samples versus the normal brain tissue samples (Figure 3A). In addition, Western blot results showed reduced expression of HMOX1 and HIF1α in T98G cells treated with sh-HMOX1 (Figure 3B).
Figure 3. HMOX1 up-regulates SPOP expression by promoting HIF1α expression. (A) Expression of HIF1α in normal brain tissue samples and GBM tissue samples analyzed by the GEPIA database. The red box represents the GBM tissue samples (n = 163) and the gray box represents the normal brain tissue samples (n = 207). (B) Transfection efficiency of sh-HMOX1 in T98G cells determined by Western blot. (C) Expression of SPOP in normal brain tissue samples and GBM tissue samples analyzed by the GEPIA database. The red box represents the GBM tissue samples (n = 163), and the gray box represents the normal brain tissue samples (n = 207). (D) Western blot of SPOP protein in tumor tissues of GBM patients (n = 20) and normal brain tissues (n = 20). (E) Western blot of SPOP protein in HEB and T98G cells. (F) Enrichment of SPOP promoter using anti-HIF1α antibody detected by ChIP assay. (G) Transfection efficiency of sh-HIF1α in T98G cells determined by Western blot. (H) Western blot of HMOX1, HIF1α and SPOP proteins in T98G cells transfected with oe-HMOX1 alone or combined with sh-HIF1α. *p < 0.05, compared to normal brain tissues, HEB cells or T98G cells transfected with sh-NC or oe-NC + sh-NC. #p < 0.05, compared to oe-HMOX1 + sh-NC-transfected T98G cells. The cell experiment was run in triplicate independently.
Previous literature has reported that HIF1α can transcriptionally activate the expression of SPOP [15], which has tumor-promoting effects in several types of cancer [15, 18, 19]. Here, we found that SPOP was significantly up-regulated in GBM samples by GEPIA analysis (Figure 3C). Meanwhile, SPOP protein expression was amplified in tumor tissues of GBM patients (Figure 3D), and in T98G cells (Figure 3E) compared with normal tissues and HEB cells, respectively.
We observed decreased enrichment of SPOP promoter in ChIP assays using specific antibody for HIF1α (Figure 3F). Furthermore, Western blot verified the transfection efficiency of sh-HIF1α in T98G cells, as shown by reduced HIF1α protein expression in T98G cells transfected with sh-HIF1α-1 or sh-HIF1α-2. The sh-HIF1α-1 sequence showed better efficiency (Figure 3G) and was used for the subsequent experiments.
In response to oe-HMOX1 + sh-NC transfection, HMOX1, HIF1α, and SPOP protein expression was increased, while further sh-HIF1α transfection resulted in decreased protein expression of HIF1α and SPOP without affecting that of HMOX1 (Figure 3H).
HMOX1 promotes GBM cell malignant characteristics through activation of the HIF1α/SPOP signaling axis
Next, this study focused on the effect of HMOX1 on GBM cell function by regulating the HIF1α/SPOP signaling axis. Western blot results presented increased expression of HMOX1, HIF1α, and SPOP in the T98G cells transfected with oe-HMOX1 + sh-NC while further transfection with sh-HIF1α reduced HIF1α and SPOP expression. Upon transfection with oe-HMOX1 + sh-SPOP, HMOX1 and HIF1α expression was unchanged in addition to diminished SPOP expression (Figure 4A).
Figure 4. HMOX1 promotes GBM cell proliferation and invasion through activation of the HIF1α/SPOP signaling axis. T98G cells were transfected with oe-HMOX1, oe-HMOX1 + sh-HIF1α or oe-HMOX1 + sh-SPOP. (A) Western blot of HMOX1, HIF1α and SPOP proteins in T98G cells. (B) Proliferation of T98G cells measured by CCK-8 assay. (C) Invasion of T98G cells measured by Transwell assay. (D) Western blot of the CCND1, E-cadherin, MYC, MMP-7, and N-cadherin proteins in T98G cells. *p < 0.05, compared to oe-NC + sh-NC-transfected T98G cells. #p < 0.05, compared to oe-HMOX1 + sh-NC-transfected T98G cells. The cell experiment was run in triplicate independently.
As shown in Figure 4B, 4C, cell proliferation and invasion were augmented upon HMOX1 overexpression, whereas contrary results were found following HIF1α silencing or SPOP silencing. Meanwhile, HMOX1 overexpression up-regulated CCND1, MYC, MMP-7 and N-cadherin protein expression while downregulating E-cadherin protein expression; in contrast, HIF1α silencing or SPOP silencing abolished these effects of HMOX1 overexpression (Figure 4D).
miR-873-5p inhibits GBM cell malignant characteristics by disrupting the HMOX1/HIF1α/SPOP signaling axis
The aforementioned results allowed us to determine whether miR-873-5p affected GBM cell function through the HMOX1/HIF1α/SPOP signaling. Transfection with miR-873-5p mimic + oe-NC in T98G cells amplified miR-873-5p expression and reduced the expression of HIF1α and SPOP. Further transfection with oe-SPOP failed to affect miR-873-5p and HIF1α expression but elevated SPOP expression (Figure 5A, 5B).
Figure 5. miR-873-5p inhibits GBM cell proliferation and invasion by inhibiting the HMOX1/HIF1α/SPOP signaling axis. T98G cells were transfected with miR-873-5p mimic alone or combined with oe-SPOP. (A) Expression of miR-873-5p, HIF1α and SPOP in T98G cells determined by RT-qPCR. (B) Western blot of HIF1α and SPOP proteins in T98G cells. (C) Western blot of the CCND1, E-cadherin, MYC, MMP-7, and N-cadherin proteins in T98G cells. (D) Proliferation of T98G cells measured by CCK-8 assay. (E) Invasion of T98G cells measured by Transwell assay. (F) Western blot of HMOX1, HIF1α, and SPOP proteins in T98G cells transfected with miR-873-5p mimic alone or combined with oe-HMOX1 or oe-HIF1α. *p < 0.05, compared to NC-mimic + oe-NC-transfected T98G cells. #p < 0.05, compared to miR-873-5p mimic + oe-NC-transfected T98G cells. The cell experiment was run in triplicate independently.
Furthermore, Western blot data exhibited a decline of CCND1, MYC, MMP-7, and N-cadherin expression, yet an increase of E-cadherin expression upon abundantly expressed miR-873-5p. However, oe-SPOP alone or combined with miR-873-5p mimic led to contrasting results (Figure 5C). In addition, decreased cell proliferation and invasion were noted upon miR-873-5p overexpression, the effect of which was rescued by oe-SPOP alone or combined with miR-873-5p mimic (Figure 5D, 5E). To further confirm that miR-873-5p regulates SPOP expression via HMOX1/HIF1α, we transfected T98G cells with miR-873-5p mimic + oe-NC, miR-873-5p mimic + oe-HMOX1 and miR-873-5p mimic +oe-HIF1α. The results of Western blot suggested higher expression of HMOX1, HIF1α and SPOP in T98G cells transfected with miR-873-5p mimic + oe-HMOX1 than those with miR-873-5p mimic + oe-NC. In addition, transfection with miR-873-5p mimic + oe-HIF1α did not change the expression of HMOX1 while augmenting that of HIF1α and SPOP versus transfection with miR-873-5p mimic + oe-NC (Figure 5F).
Overall, miR-873-5p can repress GBM cell malignant characteristics by inhibiting the HMOX1/HIF1α/SPOP signaling axis.
miR-873-5p curtails the tumorigenesis of GBM cells by inhibiting the HMOX1/HIF1α/SPOP signaling axis in vivo
Finally, the effect of miR-873-5p on the tumorigenesis of GBM cells in vivo was tested in the in vivo animal experiments. Injection of miR-873-5p agomir + oe-NC-treated T98G cells into the nude mice decreased tumor volume and weight, while miR-873-5p agomir + oe-SPOP led to opposite results (Figure 6A, 6B).
Figure 6. miR-873-5p inhibits the tumorigenesis of GBM cells by downregulating the HMOX1/HIF1α/SPOP signaling axis in vivo. (A) Tumor volume of nude mice. (B) Tumor weight of nude mice. (C) Expression of miR-873-5p, HMOX1, HIF1α, and SPOP in the tumor tissues of mice detected by RT-qPCR. (D) Western blot of the HMOX1, HIF1α, and SPOP proteins in the tumor tissues of mice. (E) Western blot of the CCND1, E-cadherin, MYC, MMP-7, and N-cadherin proteins in the tumor tissues of mice. n = 10. *p < 0.05, compared to NC agomir-treated mice. #p < 0.05, compared to miR-873-5p agomir + oe-NC-treated mice.
Furthermore, expression of miR-873-5p was increased, but that of HMOX1, HIF1α, and SPOP was reduced in the tumor tissue of mice treated with miR-873-5p agomir + oe-NC. However, further treatment with oe-SPOP caused no changes in miR-873-5p, HMOX1, and HIF1α expression and enhanced SPOP expression (Figure 6C, 6D).
A decline of CCND1, MYC, MMP-7, and N-cadherin expression along with an increase of E-cadherin expression was evident following treatment with miR-873-5p agomir + oe-NC; miR-873-5p agomir + oe-SPOP contributed to contrasting results (Figure 6E). Collectively, miR-873-5p can target HMOX1 and consequently downregulate the HIF1α/SPOP signaling axis, delaying the tumorigenesis of GBM cells in vivo.
Discussion
This study revealed that miR-873-5p can target HMOX1 and reduce its expression, and leads to consequent downregulation of the HIF1α/SPOP signaling axis, ultimately suppressing GBM cell growth and tumorigenesis, based on bioinformatics analysis together with experimental validation.
The results of this study showed that miR-873-5p was poorly expressed in GBM tissues and cell lines; in addition, the abundant expression of miR-873-5p could repress the malignant functions of GBM cells in vitro along with cell tumorigenesis in vivo. Indeed, miR-873-5p has been found to be down-regulated in GBM tissues; and overexpression of miR-873-5p exerts tumor-inhibiting roles to arrest the growth and migration of GBM cells [20]. Additionally, miR-873 has been observed to suppress the expression of IGF2BP1 and inhibit GBM cell tumorigenesis and metastasis [21].
Subsequent results showed that the inhibiting effect of miR-873-5p on GBM cell malignant behaviors was associated with the reduction of the miR-873-5p target gene HMOX1. In line with this finding, a previous study has demonstrated that miR-873-5p can target HMOX1 and negatively regulate its expression in PC12 cells [7]. HMOX1 is obviously up-regulated in GBM cells and tissues and can serve as a potential biopsy marker for GBM [22, 23]. In addition, evidence has shown that the knockdown of HMOX1 contributes to the impairment of GBM invasion [24]. Thus, miR-873-5p-mediated inhibition of HMOX1 may represent a useful target for inhibiting GBM tumor progression.
Notably, the present study revealed the possible downstream mechanism of HMOX1 promoting GBM cell proliferation and invasion, which was due to the activation of the HIF1α/SPOP signaling. In specific terms, HMOX1 promoted SPOP expression by increasing HIF1α expression and thus stimulated the malignant phenotypes of GBM cells. HMOX1 acts as upstream of HIF1α and can promote the stabilization of HIF1α [8]. Previous data indicate that SPOP protein can be regulated by HIF1α at the post-transcription level and accumulation of HIF1α augments the expression of SPOP in trophoblasts [25]. Meanwhile, amplified expression of HIF1α was noted in GBM brain tissues and cell lines; diminishing the expression of HIF1α leads to a decrease in U87MG cell viability [26]. In addition, inhibition of HIF1α can retard tumor development in a murine U251-xenograft model of GBM [27]. Another study has highlighted that the knockout of HIF1α in GBM cells eradicates cell invasion in vitro and inhibits tumor growth in vivo [28]. SPOP can suppress tumorigenesis in various human malignancies, including lung, prostate, colon, gastric, and liver cancers; however, there is also evidence revealing that SPOP exhibits oncogenic functions in kidney cancer [18]. The possible function of SPOP in GBM, however, has not yet been reported. Further investigation is still required to validate the promoting property of SPOP on GBM cell malignant phenotype along with tumorigenesis in vivo.
In addition to the above results, the current research also found that ectopically expressed miR-873-5p or HIF1α silencing could reduce the expression of N-cadherin while increasing that of E-cadherin in T98G cells. Conversely, overexpression of HMOX1 or SPOP led to opposite results. E-cadherin and N-cadherin are proteins related to EMT, which is critical in tumorigenesis due to its enhancement in chemoresistance, metastasis, and tumor stemness. Moreover, the up-regulation of N-cadherin and the down-regulation of E-cadherin are the hallmarks of EMT [29, 30]. Furthermore, existing evidence has shown that miR-873-5p overexpression can inhibit colorectal cancer cell migration, invasion, and EMT formation [31]. Another study also showed that miR-873-5p suppresses EMT and colorectal cancer cell proliferation via directly targeting JMJD8 through the NF-κB pathway [32]. In addition, HMOX1 expression has been found to significantly enhance EMT and promote lung adenocarcinoma metastasis [33]. In clear cell renal cell carcinoma, SPOP drives EMT and promotes cell invasion via activation of β-catenin/TCF4 complex [34]. Mounting studies have demonstrated that the EMT of cancer cells can be enhanced by overexpression of HIF1α [35, 36]. Therefore, miR-873-5p may contribute to the suppression of GBM EMT and invasiveness by downregulating the HMOX1/HIF1α/SPOP signaling axis.
In conclusion, the present study suggested that miR-873-5p may act as a tumor suppressor in GBM progression by downregulating the HMOX1/HIF1α/SPOP signaling axis (Figure 7). These findings provide a mechanistic understanding for developing alternative strategies for GBM treatment. Nonetheless, HIF1α is directly related to the hypoxia environment and signaling pathways [37, 38], and SPOP is a commonly studied E3 ubiquitin ligase [18, 39]. Meanwhile, another study has revealed that hypoxia-induced GLT8D1 can impede CD133 degradation via N-linked glycosylation. Direct blockade of the GLT8D1/CD133 complex formation through CD133N1-108, or inhibiting GLT8D1 expression through lercanidipine can suppress Wnt/β-catenin-dependent tumorigenesis in vitro and in a xenograft mouse model of GBM patient-derived tumors [40]. In addition, ivacaftor has been identified as an effective inhibitor of GBM stem cell maintenance, which can promote cell apoptosis and delay the growth of patient-derived xenografts from GBM. Besides, ivacaftor reduces the expression of stem cell marker genes, including CD133, CD44, and Sox2, thereby inhibiting the progression of GBM [41]. Therefore, in the context of GBM, the involvement of other hypoxia-related pathways cannot be excluded in the pronounced protection of miR-873-5p against GBM.
Figure 7. Schematic representation of the mechanism of miR-873-5p in GBM. miR-873-5p targets HMOX1 and inhibits its expression, thus downregulating the HIF1α/SPOP signaling axis, reducing the proliferation and invasion as well as the tumorigenesis of GBM cells.
Materials and Methods
In silico prediction
GBM-related miRNA expression dataset GSE103229 and mRNA expression dataset GSE104267 were retrieved from the GEO database. The GSE103229 dataset contains 5 normal brain tissue samples and 5 brain tissue samples from GBM patients. The GSE104267 dataset contains 3 control brain tissue samples and 9 brain tissue samples from GBM patients. Differentially expressed miRNAs and genes were screened using R “limma” package, and |logFC| > 1 and p value < 0.05 served as the threshold.
Downstream target genes of miR-873-5p were predicted using miRDB (Target Score > 80), TargetScan (Total context++ score < −0.20), and mirDIP databases (Score Class = Very High). The predicted genes underwent intersection analysis with the DEGs from the GSE104267 dataset. The GEPIA database was used to analyze HMOX1, HIF1α, and SPOP expression in normal brain tissue samples (n = 207) and GBM tissue samples (n = 163).
Tissue sample collection
Tumor tissue samples were collected from 20 patients with GBM receiving tumor resection between August 2020 and August 2021. All patients (aged 20-68 years with a mean age of 46 years; 12 males and 8 females) had no anti-tumor treatment before surgery. The normal brain tissues of 20 cases of accidental injuries requiring brain surgery were used as the control. No patients underwent chemotherapy or radiotherapy before surgery. The included patients were confirmed as GBM as per the 2016 WHO Classification of Tumors of the Central Nervous System and pathological examination report. Tissue specimens were stored in liquid nitrogen.
Cell culture and treatment
Human normal brain glial cell line HEB (27142244, Shanghai Aolu Biotechnology Co. Ltd., Shanghai, China) was cultured with DMEM (11965092, Gibco; consisted of 100 U/mL penicillin, 100 μg/mL streptomycin, 10% FBS). In addition, GBM cell lines (Procell Life Science and Technology Co. Ltd., Wuhan, China) T98G (CL-0583), U251 (CL-0237), LN-229 (CL-0578) and A172 (CL-0012) were cultured with EMEM (30-2003, ATCC, Manassas, VA, USA; encompassing 100 U/mL penicillin, 100 μg/mL streptomycin, 10% FBS). Cells were cultured at 37°C, 5% CO2, and 95% saturated humidity [42].
Using Lipofectamine 2000 reagent (11668019, Invitrogen), T98G cells were transfected with plasmids (GenePharma, Shanghai, China) of miR-873-5p inhibitor, miR-873-5p mimic, NC mimic, NC inhibitor, NC-mimic + oe-NC, miR-873-5p mimic + oe-NC, NC-mimic + oe-HMOX1, miR-873-5p mimic + oe-HMOX1, miR-873-5p mimic + oe-HIF1α, sh-NC, sh-HMOX1, oe-HMOX1 + sh-NC, oe-HMOX1 + sh-HIF1α and oe-HMOX1 + sh-SPOP. The transfection sequences are described in Supplementary Table 1.
RT-qPCR and Western blot
Total RNA was extracted by TRIzol (15596026, Invitrogen, USA) and reverse transcribed into cDNA by Tailing method using EZ-press miRNA Reverse Transcription Kit (EZB-miRT2-S, for miRNA detection; EZBioscience, USA) or PrimeScript RT reagent Kit (Takara, Japan; RR047A, for mRNA detection). ABI 7500 quantitative PCR instrument (Applied Biosystems, Waltham, MA, USA) was applied for RT-qPCR. 2−ΔΔCT was employed to calculate levels of target genes, normalized to U6 or GAPDH. Primer sequences are shown in Supplementary Table 2.
Total protein was extracted, separated by SDS-PAGE, and transferred onto PVDF membranes. The membrane was blocked with 5% skim milk powder and probed overnight at 4°C with primary rabbit antibodies (Supplementary Table 3) and then with secondary antibody (ab97051, 1: 2000, Abcam, HRP-labeled goat anti-rabbit IgG H&L). ECL reagent (Amersham, Little Chalfont, UK, BB-3501) was employed to visualize protein bands. Quantity One v4.6.2 software was applied for band intensity quantification, normalized to GAPDH [43].
Dual-luciferase reporter assay
The cDNA fragment of the HMOX1 3′-UTR containing miR-873-5p binding sites was cloned into the pmirGLO vector. Through the site-directed mutagenesis, the HMOX1 3′-UTR fragment (HMOX1-Mut) was cloned into the pmirGLO vector. The two plasmids were transfected with miR-873-5p mimic or NC mimic into HEK293T cells (CL-0005, Procell; cultured with DMEM containing 10% FBS at 37°C and 5% CO2) for 48 h. Next, luciferase activities were detected by a multi-mode microplate reader (SpectraMax M5, Molecular Devices, Sunnyvale, CA, USA).
CCK-8 assay
CCK-8 kit (CK04, Dojindo Laboratories, Kumamoto, Japan) was used for this experiment. GBM cells were seeded into 96-well plates (1 × 104 cells/well) and cultured for 24 h. Then cells were incubated at 37°C for 3 h with CCK-8 reagent (10 μL). OD values at 450 nm were determined.
Transwell assay
The 50 μL of Matrigel (YB356234, Yubo Biotechnology Co. Ltd., Shanghai, China) was spread in each Transwell chamber. The cells were digested, counted, and resuspended. The upper chamber was added with 200 μL cell suspension, and 800 μL conditioned medium with 20% FBS was supplemented to the lower chamber, followed by incubation in a 37°C incubator for 20–24 h. Afterward, cells were stained with 0.1% crystal violet. The invaded cells were counted under an inverted optical microscope.
ChIP assay
GBM cells were fixed to produce DNA-protein crosslinks. Then, cells were sonicated to produce chromatin fragments of an appropriate size. Lysates were then reacted with magnetic protein A beads with antibodies. Anti-HIF1α (ab2185, 1: 100, Abcam) and the NC IgG (rabbit, ab109489, 1: 100, Abcam) antibodies were used, followed by incubation with Protein Agarose/Sepharose. Finally, the precipitated SPOP (Forward: 5′-AATGAGCCTGTTGCTCCTGT-3′; Reverse: 5′-GGGCAACACCAGATGAGAAT-3′) was analyzed by RT-qPCR [44].
Construction of a mouse model of GBM
Thirty SPF BALB/c nude mice (4–6 weeks; 16–20 g) were housed in an SPF laboratory (60–65% humidity, 22–25°C, a 12-h light/dark cycle) individually for one week, with free access to water and food. The mice were acclimatized for one week before experiment. The health of the mice was observed before the experiment.
The T98G cell suspension (1 × 107 cells/mL) harboring NC-agomir + oe-NC, miR-873-5p agomir + oe-NC or miR-873-5p agomir + oe-SPOP (10 nmol/mouse, 100 μL, once 5 days, for 5 weeks) was subcutaneously inoculated into the left armpit of mice by 1 mL syringe. All overexpression vectors were constructed in the EcoRI site of the pLV vector. All agomir and overexpression plasmids were purchased from GenePharma. After 5 weeks, mice were euthanized, with the tumor removed and weighed.
Statistical analysis
Measurement data, presented as mean ± standard deviation, are analyzed using SPSS 19.0 software (IBM Corp. Armonk, NY, USA). Data between two groups were compared by unpaired t-test. One-way ANOVA with Tukey’s test was applied for multi-group data comparison. Bonferroni-corrected two-way ANOVA or repeated measures ANOVA was employed for multi-group data comparison at varied time points. p < 0.05 was considered significant.
Data availability statement
The data underlying this article will be shared on reasonable request to the corresponding author.
Supplementary Materials
Supplementary Tables
Author Contributions
X.B.Z. and F.K.J. designed the study. X.B.Z., F.K.J., C.G. and X.N.L. collated the data, carried out data analyses and produced the initial draft of the manuscript. J.N.L. and G.B.L. contributed to drafting the manuscript. All authors have read and approved the final submitted manuscript.
Conflicts of Interest
The authors declare no conflicts of interest related to this study.
Ethical Statement and Consent
The current study was approved by the Animal Ethics Committee of General Hospital of the Northern Theater Command of Chinese People’s Liberation Army and performed according to the Guide for the Care and Use of Laboratory Animals published by the US National Institutes of Health. Tissue sample collection was approved by the Clinical Ethics Committee of General Hospital of the Northern Theater Command of Chinese People’s Liberation Army and conducted in accordance with the Declaration of Helsinki. All participants signed informed consents before sample collection.
Funding
This work was funded by General Program of Liaoning Natural Science Foundation (2020-MS-056; 2022-MS-074).
References
- 1. Venkataramani V, Yang Y, Schubert MC, Reyhan E, Tetzlaff SK, Wißmann N, Botz M, Soyka SJ, Beretta CA, Pramatarov RL, Fankhauser L, Garofano L, Freudenberg A, et al. Glioblastoma hijacks neuronal mechanisms for brain invasion. Cell. 2022; 185:2899–917.e31. https://doi.org/10.1016/j.cell.2022.06.054 [PubMed]
- 2. Preddy I, Nandoliya K, Miska J, Ahmed AU. Checkpoint: Inspecting the barriers in glioblastoma immunotherapies. Semin Cancer Biol. 2022; 86:473–81. https://doi.org/10.1016/j.semcancer.2022.02.012 [PubMed]
- 3. Khan I, Baig MH, Mahfooz S, Imran MA, Khan MI, Dong JJ, Cho JY, Hatiboglu MA. Nanomedicine for glioblastoma: Progress and future prospects. Semin Cancer Biol. 2022; 86:172–86. https://doi.org/10.1016/j.semcancer.2022.06.007 [PubMed]
- 4. Zou Y, Zhong C, Hu Z, Duan S. MiR-873-5p: A Potential Molecular Marker for Cancer Diagnosis and Prognosis. Front Oncol. 2021; 11:743701. https://doi.org/10.3389/fonc.2021.743701 [PubMed]
- 5. Li K, Zhang Q, Niu D, Xing H. Mining miRNAs' Expressions in Glioma Based on GEO Database and Their Effects on Biological Functions. Biomed Res Int. 2020; 2020:5637864. https://doi.org/10.1155/2020/5637864 [PubMed]
- 6. Wang X, Li L, Zhao K, Lin Q, Li H, Xue X, Ge W, He H, Liu D, Xie H, Wu Q, Hu Y. A novel LncRNA HITT forms a regulatory loop with HIF-1α to modulate angiogenesis and tumor growth. Cell Death Differ. 2020; 27:1431–46. https://doi.org/10.1038/s41418-019-0449-8 [PubMed]
- 7. Shi R, Zhang S, Cheng G, Yang X, Zhao N, Chen C. Ginsenoside Rg1 and Acori Graminei Rhizoma Attenuates Neuron Cell Apoptosis by Promoting the Expression of miR-873-5p in Alzheimer's Disease. Neurochem Res. 2018; 43:1529–38. https://doi.org/10.1007/s11064-018-2567-y [PubMed]
- 8. Dunn LL, Kong SMY, Tumanov S, Chen W, Cantley J, Ayer A, Maghzal GJ, Midwinter RG, Chan KH, Ng MKC, Stocker R. Hmox1 (Heme Oxygenase-1) Protects Against Ischemia-Mediated Injury via Stabilization of HIF-1α (Hypoxia-Inducible Factor-1α). Arterioscler Thromb Vasc Biol. 2021; 41:317–30. https://doi.org/10.1161/ATVBAHA.120.315393 [PubMed]
- 9. Ghosh D, Ulasov IV, Chen L, Harkins LE, Wallenborg K, Hothi P, Rostad S, Hood L, Cobbs CS. TGFβ-Responsive HMOX1 Expression Is Associated with Stemness and Invasion in Glioblastoma Multiforme. Stem Cells. 2016; 34:2276–89. https://doi.org/10.1002/stem.2411 [PubMed]
- 10. Castruccio Castracani C, Longhitano L, Distefano A, Di Rosa M, Pittalà V, Lupo G, Caruso M, Corona D, Tibullo D, Li Volti G. Heme Oxygenase-1 and Carbon Monoxide Regulate Growth and Progression in Glioblastoma Cells. Mol Neurobiol. 2020; 57:2436–46. https://doi.org/10.1007/s12035-020-01869-7 [PubMed]
- 11. Zhao H, Narasimhan P, Kalish F, Wong RJ, Stevenson DK. Dysregulation of hypoxia-inducible factor-1α (Hif1α) expression in the Hmox1-deficient placenta. Placenta. 2020; 99:108–16. https://doi.org/10.1016/j.placenta.2020.07.015 [PubMed]
- 12. Domènech M, Hernández A, Plaja A, Martínez-Balibrea E, Balañà C. Hypoxia: The Cornerstone of Glioblastoma. Int J Mol Sci. 2021; 22:12608. https://doi.org/10.3390/ijms222212608 [PubMed]
- 13. Lee JW, Bae SH, Jeong JW, Kim SH, Kim KW. Hypoxia-inducible factor (HIF-1)alpha: its protein stability and biological functions. Exp Mol Med. 2004; 36:1–12. https://doi.org/10.1038/emm.2004.1 [PubMed]
- 14. Fan Y, Potdar AA, Gong Y, Eswarappa SM, Donnola S, Lathia JD, Hambardzumyan D, Rich JN, Fox PL. Profilin-1 phosphorylation directs angiocrine expression and glioblastoma progression through HIF-1α accumulation. Nat Cell Biol. 2014; 16:445–56. https://doi.org/10.1038/ncb2954 [PubMed]
- 15. Li G, Ci W, Karmakar S, Chen K, Dhar R, Fan Z, Guo Z, Zhang J, Ke Y, Wang L, Zhuang M, Hu S, Li X, et al. SPOP promotes tumorigenesis by acting as a key regulatory hub in kidney cancer. Cancer Cell. 2014; 25:455–68. https://doi.org/10.1016/j.ccr.2014.02.007 [PubMed]
- 16. Yuan D, Yang Z, Chen Y, Li S, Tan B, Yu Q. Hypoxia-induced SPOP attenuates the mobility of trophoblast cells through inhibition of the PI3K/AKT/GSK3β pathway. Cell Biol Int. 2021; 45:599–611. https://doi.org/10.1002/cbin.11501 [PubMed]
- 17. Ding D, Song T, Jun W, Tan Z, Fang J. Decreased expression of the SPOP gene is associated with poor prognosis in glioma. Int J Oncol. 2015; 46:333–41. https://doi.org/10.3892/ijo.2014.2729 [PubMed]
- 18. Song Y, Xu Y, Pan C, Yan L, Wang ZW, Zhu X. The emerging role of SPOP protein in tumorigenesis and cancer therapy. Mol Cancer. 2020; 19:2. https://doi.org/10.1186/s12943-019-1124-x [PubMed]
- 19. Clark A, Burleson M. SPOP and cancer: a systematic review. Am J Cancer Res. 2020; 10:704–26. [PubMed]
- 20. Lin YH, Guo L, Yan F, Dou ZQ, Yu Q, Chen G. Long non-coding RNA HOTAIRM1 promotes proliferation and inhibits apoptosis of glioma cells by regulating the miR-873-5p/ZEB2 axis. Chin Med J (Engl). 2020; 133:174–82. https://doi.org/10.1097/CM9.0000000000000615 [PubMed]
- 21. Wang RJ, Li JW, Bao BH, Wu HC, Du ZH, Su JL, Zhang MH, Liang HQ. MicroRNA-873 (miRNA-873) inhibits glioblastoma tumorigenesis and metastasis by suppressing the expression of IGF2BP1. J Biol Chem. 2015; 290:8938–48. https://doi.org/10.1074/jbc.M114.624700 [PubMed]
- 22. Gao Y, Zhang E, Liu B, Zhou K, He S, Feng L, Wu G, Cao M, Wu H, Cui Y, Zhang X, Liu X, Wang Y, et al. Integrated analysis identified core signal pathways and hypoxic characteristics of human glioblastoma. J Cell Mol Med. 2019; 23:6228–37. https://doi.org/10.1111/jcmm.14507 [PubMed]
- 23. Ozyerli-Goknar E, Sur-Erdem I, Seker F, Cingöz A, Kayabolen A, Kahya-Yesil Z, Uyulur F, Gezen M, Tolay N, Erman B, Gönen M, Dunford J, Oppermann U, Bagci-Onder T. The fungal metabolite chaetocin is a sensitizer for pro-apoptotic therapies in glioblastoma. Cell Death Dis. 2019; 10:894. https://doi.org/10.1038/s41419-019-2107-y [PubMed]
- 24. Ghosh D, Funk CC, Caballero J, Shah N, Rouleau K, Earls JC, Soroceanu L, Foltz G, Cobbs CS, Price ND, Hood L. A Cell-Surface Membrane Protein Signature for Glioblastoma. Cell Syst. 2017; 4:516–29.e7. https://doi.org/10.1016/j.cels.2017.03.004 [PubMed]
- 25. Duan A, Shen K, Li B, Li C, Zhou H, Kong R, Shao Y, Qin J, Yuan T, Ji J, Guo W, Wang X, Xue T, et al. Extracellular vesicles derived from LPS-preconditioned human synovial mesenchymal stem cells inhibit extracellular matrix degradation and prevent osteoarthritis of the knee in a mouse model. Stem Cell Res Ther. 2021; 12:427. https://doi.org/10.1186/s13287-021-02507-2 [PubMed]
- 26. Nie W, Luo X, Lu D, Yuan P, Liu B, Xu H, Ye M. Casein kinase 1α 1 is involved in the progression of glioblastoma through HIF-1α-mediated autophagy. J Neurophysiol. 2022; 128:910–8. https://doi.org/10.1152/jn.00316.2022 [PubMed]
- 27. Lee J, Kim E, Chong K, Ryu SW, Kim C, Choi K, Kim JH, Choi C. Atypical induction of HIF-1α expression by pericellular Notch1 signaling suffices for the malignancy of glioblastoma multiforme cells. Cell Mol Life Sci. 2022; 79:537. https://doi.org/10.1007/s00018-022-04529-2 [PubMed]
- 28. Peng G, Wang Y, Ge P, Bailey C, Zhang P, Zhang D, Meng Z, Qi C, Chen Q, Chen J, Niu J, Zheng P, Liu Y, Liu Y. The HIF1α-PDGFD-PDGFRα axis controls glioblastoma growth at normoxia/mild-hypoxia and confers sensitivity to targeted therapy by echinomycin. J Exp Clin Cancer Res. 2021; 40:278. https://doi.org/10.1186/s13046-021-02082-7 [PubMed]
- 29. Loh CY, Chai JY, Tang TF, Wong WF, Sethi G, Shanmugam MK, Chong PP, Looi CY. The E-Cadherin and N-Cadherin Switch in Epithelial-to-Mesenchymal Transition: Signaling, Therapeutic Implications, and Challenges. Cells. 2019; 8:1118. https://doi.org/10.3390/cells8101118 [PubMed]
- 30. Canesin G, Cuevas EP, Santos V, López-Menéndez C, Moreno-Bueno G, Huang Y, Csiszar K, Portillo F, Peinado H, Lyden D, Cano A. Lysyl oxidase-like 2 (LOXL2) and E47 EMT factor: novel partners in E-cadherin repression and early metastasis colonization. Oncogene. 2015; 34:951–64. https://doi.org/10.1038/onc.2014.23 [PubMed]
- 31. Li G, Xu Y, Wang S, Yan W, Zhao Q, Guo J. MiR-873-5p inhibits cell migration, invasion and epithelial-mesenchymal transition in colorectal cancer via targeting ZEB1. Pathol Res Pract. 2019; 215:34–9. https://doi.org/10.1016/j.prp.2018.10.008 [PubMed]
- 32. Wang L, Jiang F, Ma F, Zhang B. MiR-873-5p suppresses cell proliferation and epithelial-mesenchymal transition via directly targeting Jumonji domain-containing protein 8 through the NF-κB pathway in colorectal cancer. J Cell Commun Signal. 2019; 13:549–60. https://doi.org/10.1007/s12079-019-00522-w [PubMed]
- 33. Chen B, Zhang L, Zhou H, Ye W, Luo C, Yang L, Fang N, Tang A. HMOX1 promotes lung adenocarcinoma metastasis by affecting macrophages and mitochondrion complexes. Front Oncol. 2022; 12:978006. https://doi.org/10.3389/fonc.2022.978006 [PubMed]
- 34. Zhao W, Zhou J, Deng Z, Gao Y, Cheng Y. SPOP promotes tumor progression via activation of β-catenin/TCF4 complex in clear cell renal cell carcinoma. Int J Oncol. 2016; 49:1001–8. https://doi.org/10.3892/ijo.2016.3609 [PubMed]
- 35. Kang MK, Kim SI, Oh SY, Na W, Kang YH. Tangeretin Ameliorates Glucose-Induced Podocyte Injury through Blocking Epithelial to Mesenchymal Transition Caused by Oxidative Stress and Hypoxia. Int J Mol Sci. 2020; 21:8577. https://doi.org/10.3390/ijms21228577 [PubMed]
- 36. Zhang J, Zhang Q, Lou Y, Fu Q, Chen Q, Wei T, Yang J, Tang J, Wang J, Chen Y, Zhang X, Zhang J, Bai X, Liang T. Hypoxia-inducible factor-1α/interleukin-1β signaling enhances hepatoma epithelial-mesenchymal transition through macrophages in a hypoxic-inflammatory microenvironment. Hepatology. 2018; 67:1872–89. https://doi.org/10.1002/hep.29681 [PubMed]
- 37. Rezaeian AH, Li CF, Wu CY, Zhang X, Delacerda J, You MJ, Han F, Cai Z, Jeong YS, Jin G, Phan L, Chou PC, Lee MH, et al. A hypoxia-responsive TRAF6-ATM-H2AX signalling axis promotes HIF1α activation, tumorigenesis and metastasis. Nat Cell Biol. 2017; 19:38–51. https://doi.org/10.1038/ncb3445 [PubMed]
-
38.
Reglero C, Lafarga V, Rivas V, Albitre Á, Ramos P, Berciano SR, Tapia O, Martínez-Chantar ML, Mayor F
Jr , Penela P. GRK2-Dependent HuR Phosphorylation Regulates HIF1α Activation under Hypoxia or Adrenergic Stress. Cancers (Basel). 2020; 12:1216. https://doi.org/10.3390/cancers12051216 [PubMed] - 39. Zhang H, Jin X, Huang H. Deregulation of SPOP in Cancer. Cancer Res. 2023; 83:489–99. https://doi.org/10.1158/0008-5472.CAN-22-2801 [PubMed]
- 40. Liu K, Jiang L, Shi Y, Liu B, He Y, Shen Q, Jiang X, Nie Z, Pu J, Yang C, Chen Y. Hypoxia-induced GLT8D1 promotes glioma stem cell maintenance by inhibiting CD133 degradation through N-linked glycosylation. Cell Death Differ. 2022; 29:1834–49. https://doi.org/10.1038/s41418-022-00969-2 [PubMed]
- 41. Liu K, Pu J, Nie Z, Shi Y, Jiang L, Wu Q, Chen Y, Yang C. Ivacaftor Inhibits Glioblastoma Stem Cell Maintenance and Tumor Progression. Front Cell Dev Biol. 2021; 9:678209. https://doi.org/10.3389/fcell.2021.678209 [PubMed]
- 42. Lu J, Zhang M, Yang X, Cui T, Dai J. MicroRNA-548c-3p inhibits T98G glioma cell proliferation and migration by downregulating c-Myb. Oncol Lett. 2017; 13:3866–72. https://doi.org/10.3892/ol.2017.5870 [PubMed]
-
43.
Huan J, Hornick NI, Shurtleff MJ, Skinner AM, Goloviznina NA, Roberts CT
Jr , Kurre P. RNA trafficking by acute myelogenous leukemia exosomes. Cancer Res. 2013; 73:918–29. https://doi.org/10.1158/0008-5472.CAN-12-2184 [PubMed] - 44. Zhu P, Wang Y, Wu J, Huang G, Liu B, Ye B, Du Y, Gao G, Tian Y, He L, Fan Z. LncBRM initiates YAP1 signalling activation to drive self-renewal of liver cancer stem cells. Nat Commun. 2016; 7:13608. https://doi.org/10.1038/ncomms13608 [PubMed]