- Department of Neurosurgery, Ningbo City First Hospital, Ningbo, China
Glioblastoma (GBM), an aggressive primary tumor, is common in humans, accounting for 12–15% of all intracranial tumors, and has median survival of fewer than 15 months. Since a growing body of evidence suggests that conventional drugs are ineffective against GBM, our goal is to find emerging therapies that play a role in its treatment. This research constructs a risk model to predict the prognosis of GBM patients. A set of genes associated with GBM was taken from a GBM gene data bank, and clinical information on patients with GBM was retrieved from the Cancer Genome Atlas (TCGA) data bank. One-way Cox and Kaplan–Meier analyses were performed to identify genes in relation to prognosis. Groups were classified into high and low expression level of PTEN expression. Prognosis-related genes were further identified, and multi-factor Cox regression analysis was used to build risk score equations for the prognostic model to construct a survival prognostic model. The area under the ROC curve suggested that the pattern had high accuracy. When combined with nomogram analysis, GJB2 was considered an independent predictor of GBM prognosis. This study provides a potential prognostic predictive biological marker for GBM patients and confirms that GJB2 is a key gene for GBM progression.
Introduction
Glioblastoma (GBM) is an aggressive primary tumor that is the most aggressive and common type of brain tumor in humans (1), with a median survival of fewer than 15 months (2). It consists of highly malignant cells with metastatic and angiogenic properties that lead to resistance to agents such as temozolomide, although none of these agents significantly improves overall survival (3). The conventional treatment is a combination of chemotherapy and surgical resection followed by radiotherapy and adjuvant chemotherapy. This regimen has been effective in improving overall survival but largely fails to prevent recurrence because surgical treatment fails to completely eradicate GBM cells, which remain surrounded by scattered GBM-infiltrating cells (4). Therefore, the identification of reliable prognostic markers has become crucial for GBM treatment.
Phosphatase and tensin homolog (PTEN), a tumor suppressor gene, is closely involved in cell translation, proliferation, and tumorigenesis. PTEN gene mutations are frequently found in the genetic landscape of high-grade gliomas, and is hallmarks of glioma malignancy,they influence cell proliferation, proangiogenetic pathways, and antitumoral immune response (5).This study evaluated the prognosis and expression of GBM patients by analyzing the mutation spectrum of PTEN-related microenvironment.
We collected a set of genes related to GBM. Candidate genes were obtained by differentially expressed genes. One-way Cox and Kaplan–Meier analyses were performed to identify the genes. Combined with LASSO regression analysis to establish prognostic features and further identify prognosis-related genes, a survival prognostic model was constructed, and the area under the ROC curve suggested that the model had high accuracy. Combined with nomogram analysis, the prognostic significance of genetic features in GBM was assessed.
Materials and method
Data acquisition
Expression matrices of all tumors, including the clinically relevant pathological features of tumor tissue samples, were obtained from the Cancer Genome Atlas (TCGA) database and evaluated based on the data completeness of clinical samples and degree of matching with sequenced samples. Duplicate and censored samples and cases without clinical findings were excluded.
Variance analysis
The Limma package of R software can screen for differentially expressed genes between GBM and para cancer with |logFC| = 0.3785 (adj. P-value <0.05) and plotted volcanoes.To further confirm the potential functions of potential targets,the data were analyzed by functional enrichment. ClusterProfiler program package in R software was used to analyze the GO Function of potential mRNA and KEGG pathway gene ontology (GO) containing Molecular functions (MF), Biological process (BP) and Cellular Component (CC).
Kaplan–Meier survival analysis
The effect of genetic characteristics on prognosis was verified by One-way Cox. Kaplan–Meier survival curves were plotted to compare the survival of patients in the high and low expression level groups.
Establishment and analysis of risk prediction models
LASSO regression analysis was applied to narrow the range of prognosis-related genes and ensure the stability of the results. The median risk score of each sample was calculated using the risk score formula as the threshold, and patients were classified into high expression level and low expression level groups according to the expression level of PTEN. Survival curves were plotted. In addition, the working curves (ROC) of the subjects were plotted, and the area under the curve (AUC) was calculated to assess the predictive validity of the model.
Nomogram analysis
Independent prognostic factors for GBM were determined by univariate and multifactorial analyses. Column plots and calibration curve plots predicted the predictive power of survival at 1, 2, and 3 years.
Pan-cancer analysis
Data from normal and TCGA tumor tissues in the GTEx database were combined to analyze the differences in their gene expression. The amount of tumor mutations in each tumor sample was counted separately, and the relationship between gene expression and tumor mutational burden was analyzed.
Results
Mutation landscape
PTEN gene mutation data, clinical data, and transcriptome data from the TCGA data bank and somatic mutations in GBM patients were downloaded. Mutation analysis revealed that PTEN had a high mutation rate of 30.03% in GBM, which was the highest among TCGA tumors (Figure 1A). PTEN mutation frequency ranked first (34.1%) among all mutated genes in GBM (Figure 1B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/952521/fonc-12-952521-HTML/image_m/fonc-12-952521-g001.jpg)
Figure 1 Somatic mutation rate and mutation landscape map of PTEN in GBM. (A) Lollipop plot of PTEN mutation distribution in GBM; (B) oncoplot showing the somatic landscape of the GBM cohort.
Differentially expressed gene screening
In | logFC | = 0.58 under the condition of the filter, identified the genes from the TCGA - GBM data too little number, change | logFC | = 0.3785, identified 45 differentially expressed genes (adj. P-value <0.05). There were 5 down-regulated genes and 40 up-regulated genes (Figures 2A, B). The results of differentially up-regulated genes and down-regulated genes, KEGG pathway enrichment, and GO term enrichment are shown in Figure 2C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/952521/fonc-12-952521-HTML/image_m/fonc-12-952521-g002.jpg)
Figure 2 Screening for differentially expressed genes. (A, B) Volcano plot showing differentially expressed genes in GBM; (C) KEGG pathway enrichment results and GO term enrichment results.
Prognostic analysis
To understand the effect of each gene on prognosis, one-way Cox regression was used to validate the patients into high and low expression level groups by scoring each gene’s median expression level as the median value. It was further observed that only KIRREL2, TENT5B, DIRAS3, SDC1, GJB2, DDIT4L, HOXA10, and H2AW were prognostically significant in GBM (Figure 3A). The relationship between the prognosis of GBM and the expression levels of the above genes was further analyzed by plotting Kaplan–Meier survival curves. It was found that the higher the expression of DDIT4L, GJB2, KIRREL2, DIRAS3, HOXA10, SDC1, and TENT5B, the worse the prognosis, and the lower the expression of H2AW, the worse the prognosis (Figure 3B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/952521/fonc-12-952521-HTML/image_m/fonc-12-952521-g003.jpg)
Figure 3 Prognostic analysis, (A) Validation of gene effects on prognosis using one-way Cox regression, presented using forest plots; (B) Kaplan–Meier survival analysis of the relationship between high and low expression levels of DDIT4L, DIRAS3, GJB2, H2AW, HOXA10, KIRREL2, SDC1, and TENT5B genes and GBM prognosis.
Expression analysis
The expression of DDIT4L, DIRAS3, GJB2, H2AW, HOXA10, KIRREL2, SDC1, and TENT5B in the PTEN mutant group, PTEN wild group, and normal group were observed, and the box plot identified that DDIT4L, DIRAS3, GJB2, HOXA10, SDC1, and TENT5B were up-regulated in PTEN mutant group and PTEN wild group. H2AW and KIRREL2 were down-regulated in PTEN mutant group and PTEN wild group, compared with the normal group (Figure 4).
Construction of the LASSO prognostic model
LASSO regression analysis was performed to further narrow the range of prognosis-related genes and ensure the stability of the results (Figures 5A, B). The risk score was calculated for each sample using the following risk score formula:
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/952521/fonc-12-952521-HTML/image_m/fonc-12-952521-g005.jpg)
Figure 5 LASSO regression analysis. (A) Coefficients of selected features expressed as λ parameters; (B) partial likelihood deviation plotted with log(l) using LASSO regression analysis model; (C) Kaplan–Meier survival curves for patients in the high and low expression level groups; (D) time-dependent ROC curve plot.
The median expression level score was used as a threshold for dividing the patients into high and low expression level groups and plotting the survival curves, which showed that the survival prognosis of the high expression level group was significantly worse than that of the low expression level group (Figure 5C). The ROC curves of 1, 3, and 5-year survival of GBM patients predicted by this risk score model were plotted (Figure 5D), and their AUC areas were 0.694,0.651, and 0.763, respectively.
Nomogram analysis
The GJB2 gene was shown to be a prognostic factor for GBM by univariate and multifactorial analyses (Figures 6A, B). Column and calibration curve plots indicated the predictive power of survival at 1, 2, and 3 years. The C-index was 0.583 and the P-value was 0.02.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/952521/fonc-12-952521-HTML/image_m/fonc-12-952521-g006.jpg)
Figure 6 Nomogram analysis. (A) Univariate and (B) multifactorial Cox analyses showing that the GJB2 gene is an independent prognostic factor for GBM; (C) columnar plots predicting 1, 2, and 3-year overall survival in GBM patients; (D) calibration curves for the overall survival columnar plot model.
Expression of GJB2 in GBM
The expression of GJB2 in tumors was analyzed by integrating the data from TCGA, and the results showed that GJB2 was highly expressed in BLCA, BRCA, CESC, COAD, DLBC, ESCA, GBM, HNSC, KIRC, KIRP, LGG, LIHC, LUAD, LUSC, OV, PAAD, PRAD, STAD, TGCT, THCA, UCEC, and UCS, and there was low expression in CHOL, KICH, and READ (Figure 7A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/952521/fonc-12-952521-HTML/image_m/fonc-12-952521-g007.jpg)
Figure 7 Pan-cancer analysis of GJB2 in multiple tumors. (A) Expression of GJB2 in multiple tumors; (B) prognostic analysis of GJB2 in multiple tumors; (C) correlation of GJB2 with immune checkpoints in multiple tumors; (D) correlation of GJB2 with immune mutation load in multiple tumors; (E) correlation of GJB2 with microsatellites in multiple tumors. *P<0.05, **P<0.01,***P<0.001.
One-way Cox regression showed prognostic significance of GJB2 in KICH and PCPG (Figure 7B).
The relationship between GJB2 gene expression and immune checkpoint gene expression was analyzed, and CD276 was highly correlated with multiple cancers (Figure 7C).
Immune mutation load correlation indicated that GJB2 was significantly associated with immune mutation load in CESC, COAD, KIRC, KIRP, LAML, LUAD, PAAD, SKCM, STAD, and UCEC (Figure 7D).
GJB2 was correlated with microsatellite instability in COAD, DLBC, ESCA, HNSC, KIRC, LIHC, and PRAD (Figure 7E).
Discussion
As contemporary living standards and the material conditions of people continue to improve, health issues are receiving more and more attention. However, the incidence of cancer is increasing year by year, threatening human health. Among the cancers, GBM, which is classified as a grade IV diffuse glioma, is the primary tumor in adults and has a terrible prognosis. It has a high rate of recurrence and can spread rapidly to other parts of the brain (6), causing thousands of deaths worldwide each year. GBM consists of highly malignant cells that have metastatic and angiogenic properties that lead to resistance to agents such as temozolomide (3). The average survival of GBM patients remains below 20 months using current therapies (7), and conventional systemic chemotherapeutic agents used to treat GBM are ineffective. Several studies have found that the extremely high recurrence rate of GBM is associated with the expression of strongly proliferative genes in cells (8). Since this process usually involves multiple genes (9), we aimed to explore the prognostic genes that may play a role in the treatment of GBM.
PTEN, one of the most frequently mutated genes in human cancers (10), is a tumor suppressor with growth and survival regulatory functions (11). It controls many processes, including survival, proliferation, differentiation, energy metabolism, and deregulation of cell structure and mobility (12, 13). Deletion or mutation of PTEN leads to increased cell proliferation and decreased cell death and tumor development (14). Somatic loss-of-function mutations in PTEN can cause or drive the malignant phenotype of various human cancers (13), which is consistent with our findings. GJB2 is considered an oncogene and is related to tumor growth, EMT, and lymph node metastasis in a variety of cancers (15–18). Mutations in the GJB2 gene are a major cause of autosomal recessive hereditary nonsyndromic hearing loss (ARNSHL) in many populations (19). GJB2 expression is elevated in many tumor cell lines, tumor tissues (20), and breast cancer (21). It has been shown that GJB2 is an independent prognostic biomarker for LUAD, and patients with GJB2 overexpression have shorter overall survival. GJB2 may also be a potential prognostic factor for KIRC, as shown by pan-cancer analysis (22).
We identified genes related to the prognosis of GBM, namely GJB2. We first screened differentially expressed genes by differential analysis and performed one-way Cox regression and Kaplan–Meier analyses to identify genes related to prognosis. The groups were classified into high and low expression level, based on PTEN expression. Combined with LASSO regression analysis to establish prognostic characteristics and further identify prognosis-related genes, multi-factor Cox regression analysis was used to establish risk score equations for the prognostic model to construct the survival prognostic model, and the area under the ROC curve proved that the model had high accuracy. Combined with nomogram analysis, GJB2 was considered an independent predictor of GBM prognosis.
Many studies on the genetic correlation of GBM prognosis have been conducted. For example, Yanxin Li et al. found that the prognosis of GBM patients was poorer when the HOXD10 gene was highly expressed, and HOXD10 may play different roles at different stages of GBM development (23). To the best of our knowledge, the GJB2 gene that we screened is a new GBM biomarker, and no previous reports of this gene associated with the development and progression of GBM have appeared.
In summary, this study provides a predictive biological marker for GBM patients and confirms GJB2 as a key gene for GBM progression. These conclusions may provide a direction for prognosis prediction and treatment of GBM patients.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.
Funding
This study was supported by the grants from the Ningbo medical and health brand discipline (PPXK2018-04) and Ningbo Science and Technology Innovation 2025 Major Project (2020Z094, 2022Z125), Key Laboratory of Precision Medicine for Atherosclerotic Diseases of Zhejiang Province (2022E10026).
Author contributions
We contributed equally for this work.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Yang Y, Tu Y, Lu J, Chen Q, Zhu Z, Peng W, et al. PT109, a novel multi-kinase inhibitor, suppresses glioblastoma multiforme through cell reprogramming: Involvement of PTBP1/PKM1/2 pathway. Eur J Pharmacol (2022) 23:174837. doi: 10.1016/j.ejphar.2022.174837
2. Shergalis A, Bankhead A 3rd, Luesakul U, Muangsin N, Neamati N. Current challenges and opportunities in treating glioblastoma. Pharmacol Rev (2018) 703:412–45. doi: 10.1124/pr.117.014944
3. Carlsson SK, Brothers SP, Wahlestedt C. Emerging cellular therapies for glioblastoma multiforme. Cureus (2018) 10(3):e2305. doi: 10.7759/cureus.2305
4. Vasilev A, Sofi R, Rahman R, Smith SJ, Teschemacher AG, Kasparov S, et al. Using light for therapy of glioblastoma multiforme (GBM). Brain Sci (2020) 10(3):75. doi: 10.3390/brainsci10020075
5. Giotta Lucifero A, Luzzi S. Immune landscape in PTEN-related glioma microenvironment: A bioinformatic analysis. Brain Sci (2022) 12(4):501. doi: 10.3390/brainsci12040501
6. Batash R, Asna N, Schaffer P, Francis N, Schaffer M, et al. Glioblastoma multiforme, diagnosis and treatment; recent literature review. Curr Medicin Chem (2017) 24(27):3002–9. doi: 10.2174/0929867324666170516123206
7. Wang S, Yao F, Lu X, Li Q, Su Z, Lee JH, et al. Temozolomide promotes immune escape of GBM cells via upregulating PD-L1. Am J Cancer Res (2019) 9(6):1161–71.
8. Lara-Velazquez M, Zarco N, Carrano A, Phillipps J, Norton ES, Schiapparelli P, et al. Alpha 1-antichymotrypsin contributes to stem cell characteristics and enhances tumorigenicity of glioblastoma. Neuro-oncology (2021) 23(4):599–610. doi: 10.1093/neuonc/noaa264
9. Malik V, Garg S, Afzal S, Dhanjal JK, Yun CO, Kaul SC, et al. Bioinformatics and molecular insights to anti-metastasis activity of triethylene glycol derivatives. Int J Mol Sci (2020) 21(15):5463. doi: 10.3390/ijms21155463
10. Zheng L, Liang H, Zhang Q, Shen Z, Sun Y, Zhao X, et al. circPTEN1, a circular RNA generated from PTEN, suppresses cancer progression through inhibition of TGF-β/Smad signaling. Mol Cancer (2022) 21(1):41. doi: 10.1186/s12943-022-01495-y
11. Chen C-Y, Chen J, He L, Stiles BL. PTEN: Tumor suppressor and metabolic regulator. Front Endocrinol (2018) 9:338. doi: 10.3389/fendo.2018.00338
12. Chow JT, Salmena L. Recent advances in PTEN signalling axes in cancer. Facult Rev (2020) 9:31. doi: 10.12703/r/9-31
13. Chen L, Deyin G. The functions of tumor suppressor PTEN in innate and adaptive immunity. Cell Mol Immunol (2017) 14(7):581–9. doi: 10.1038/cmi.2017.30
14. Song MS, Salmena L, Pandolfi PP. The functions and regulation of the PTEN tumour suppressor. Nat Rev Mol Cell Biol (2012) 13(5):283–96. doi: 10.1038/nrm3330
15. Ezumi K, Yamamoto H, Murata K, Higashiyama M, Damdinsuren B, Nakamura Y, et al. Aberrant expression of connexin 26 is associated with lung metastasis of colorectal cancer. Clin Cancer Res (2008) 14(3):677–84. doi: 10.1158/1078-0432.CCR-07-1184
16. Inose T, Kato H, Kimura H, Faried A, Tanaka N, Sakai M, et al. Correlation between connexin 26 expression and poor prognosis of esophageal squamous cell carcinoma. Ann Surg Oncol (2009) 16(6):1704–10. doi: 10.1245/s10434-009-0443-3
17. Naoi Y, Miyoshi Y, Taguchi T, Kim SJ, Arai T, Maruyama N, et al. Connexin26 expression is associated with aggressive phenotype in human papillary and follicular thyroid cancers. Cancer Lett (2008) 262(2):248–56. doi: 10.1016/j.canlet.2007.12.008
18. Wu J-I, Wang LH. Emerging roles of gap junction proteins connexins in cancer metastasis, chemoresistance and clinical application. J Biomed Sci (2019) 26(1):8. doi: 10.1186/s12929-019-0497-x
19. Koohiyan M, Koohian F, Azadegan-Dehkordi F. GJB2-related hearing loss in central Iran: Review of the spectrum and frequency of gene mutations. Ann Hum Genet (2020) 84(2):107–13. doi: 10.1111/ahg.12354
20. Zhu T, Gao YF, Chen YX, Wang ZB, Yin JY, Mao XY, et al. Genome-scale analysis identifies GJB2 and ERO1LB as prognosis markers in patients with pancreatic cancer. Oncotarget (2017) 8(13):21281–9. doi: 10.18632/oncotarget.15068
21. Shettar A, Damineni S, Mukherjee G, Kondaiah P. Gap junction β- 2 expression is negatively associated with the estrogen receptor status in breast cancer tissues and is a regulator of breast tumorigenesis. Oncol Rep (2018) 40(6):3645–53. doi: 10.3892/or.2018.6764
22. Lu Ao, Shi Y, Liu Y, Lin J, Zhang H, Guo Y, et al. Integrative analyses identified ion channel genes GJB2 and SCNN1B as prognostic biomarkers and therapeutic targets for lung adenocarcinoma. Lung Cancer (Amsterd Netherl) (2021) 158:29–39. doi: 10.1016/j.lungcan.2021.06.001
Keywords: prognosis, LASSO, pan-cancer analysis, GBM, glioblastoma
Citation: Zhou S, Wang H, Huang Y, Wu Y and Lin Z (2022) The global change of gene expression pattern caused by PTEN mutation affects the prognosis of glioblastoma. Front. Oncol. 12:952521. doi: 10.3389/fonc.2022.952521
Received: 25 May 2022; Accepted: 30 June 2022;
Published: 09 August 2022.
Edited by:
Ye Wang, The Second Affiliated Hospital of Medical College of Qingdao University, ChinaReviewed by:
Limin Xu, Huzhou Central Hospital, ChinaShengbao Wang, Lanzhou University Second Hospital, China
Copyright © 2022 Zhou, Wang, Huang, Wu and Lin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhiqing Lin, zhiqinglinly@outlook.com
†These authors share first authorship