Abstract
Glioblastoma (GBM) is the most common malignant tumor of the central nervous system. GBM causes poor clinical outcome and high mortality rate, mainly due to the lack of effective targeted therapy and prognostic biomarkers. Here, we developed a user-friendly Online Survival analysis web server for GlioBlastoMa, abbreviated OSgbm, to assess the prognostic value of candidate genes. Currently, OSgbm contains 684 samples with transcriptome profiles and clinical information from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) and Chinese Glioma Genome Atlas (CGGA). The survival analysis results can be graphically presented by Kaplan-Meier (KM) plot with Hazard ratio (HR) and log-rank p value. As demonstration, the prognostic value of 51 previously reported survival associated biomarkers, such as PROM1 (HR = 2.4120, p = 0.0071) and CXCR4 (HR = 1.5578, p < 0.001), were confirmed in OSgbm. In summary, OSgbm allows users to evaluate and develop prognostic biomarkers of GBM. The web server of OSgbm is available at http://bioinfo.henu.edu.cn/GBM/GBMList.jsp.
Introduction
Glioblastoma (GBM) is the most common malignant tumor of the central nervous system (CNS) and causes a high mortality rate (Nikiforova and Hamilton, 2011; Stoyanov et al., 2018). Although many new therapies have improved the clinical outcome and more clinical trials have demonstrated the high efficacy in treating GBM, the survival rate of GBM patients is still low. GBM is a complex disease to tackle with a median survival period of approximately 14 months, and a 5-year survival rate of 5% (Stupp et al., 2005; Johnson and O'Neill, 2012; Polivka et al., 2017). Prognostic biomarkers have been showing great roles in cancer patient management and may guide targeted therapies. Therefore, it is greatly needed to investigate prognostic biomarkers in GBM.
Previous studies have reported some prognostic biomarkers in GBM, such as gene mutation of gene IDH and PTEN, and expression variation of gene CD133 (Yang et al., 2016; Cai and Sughrue, 2017; Nguyen et al., 2018). However, these biomarkers have not been translated to clinical applications due to the lack of independent validation. In addition, due to the molecular heterogeneity among GBMs and limited patient samples (Nathanson et al., 2014; Aldape et al., 2015; Brown et al., 2017), the prognostic behavior of a certain biomarker may be inconsistent or even contradictory between different reports. In other words, cross population validation in a larger patient cohort is critical for evaluating the prognostic biomarker.
In current work, we collected the gene expression profiles and clinical information of 684 GBM patients from seven independent cohorts obtained from TCGA, GEO and CGGA. We developed a user-friendly web server, OSgbm, to analyze the prognostic value of genes of interests. With this web server, it would facilitate researchers and clinicians to screen, develop and validate new prognostic biomarkers in GBM.
Methods
Datasets Collection
GBM datasets are from three major data sources. First, level-3 gene expression profiling data (HiSeqV2) and clinical information of GBM samples were downloaded from TCGA on April 2018 (https://portal.gdc.cancer.gov/). Second, four cohorts (≥30 cases) with available gene expression profiles and clinical survival information were collected from GEO database (http://www.ncbi.nlm.nih.gov/geo/). Third, two GBM cohorts were gathered from CGGA (http://www.cgga.org.cn/). After an initial filtration and quality check (with available gene expression profiling data and clinical survival information), 153 samples from TCGA, 276 samples from GEO, and 255 samples from CGGA were included for the following database and web server construction. The histology of recurrent GBM (rGBM) were included in GSE7696 (10 samples), GSE42669 (11 samples), CGGAarray (9 samples) and CGGAseq (22 samples) datasets. Two CGGA datasets also included 20 samples of secondary GBM (sGBM).
System Implementation and Server Set-Up
OSgbm is a web-based tool which uses J2EE (Java 2 Platform Enterprise Edition) architecture as we previously described (Wang et al., 2019; Wang et al., 2019; Xie et al., 2019a; Zhang et al., 2019). The gene expression and clinical data were integrated in the background database, which was handled by a MySQL server. Dynamic web interfaces were written in HTML 5.0 and hosted by Tomcat on Windows Server. Using OSgbm requires a HTML 5.0-compliant browser with JavaScript enabled, but does not require any particular visual plug-in tool. Since the web server was designed for users with no specialized bioinformatics skills, we propose ‘out-of-the-box’ data. The input of OSgbm web server is official gene symbol. For the “Data Source: Combined” option, as all the datasets used in OSgbm already have been published, processed and normalized well, in order to avoid of the batch effect and platform biases among these datasets, we first stratify the patients into high- and low-expression group for the input gene in each dataset, and then merged relative patients from high- and low-expression group from each dataset into a combined high-expression group (Upper group in the Kaplan–Meier plot) and a combined low-expression group (Lower group in the Kaplan–Meier plot) for the analysis of Kaplan–Meier plot and log-rank test. The statistical analyses of input were performed with R package: KM curves with Hazard ratio (HR, 95% confidence interval) and log-rank p value were calculated by R package ‘survival’. OSgbm is available at http://bioinfo.henu.edu.cn/GBM/GBMList.jsp.
Validation of Previously Reported Prognostic Biomarkers
A PubMed search was performed to identify previously reported GBM prognostic biomarkers, using keywords ‘glioblastoma’, ‘survival’ and ‘biomarker’. Totally, 53 prognostic biomarkers were identified from 2013 publications. The flow chart of biomarker collection was showed in Figure S1. The prognostic values of these published biomarkers were analyzed in either a form of combined cohorts of all GBM patients or in a single cohort in our database.
Results
The Clinical Characteristics of GBM Datasets Used in OSgbm
In OSgbm, we included a total of 684 unique GBM samples from seven datasets, including one TCGA cohort, four GEO cohorts and two CGGA cohorts. The survival information includes overall survival (OS), disease specific survival (DSS), disease free interval (DFI) and progression free interval (PFI) (Liu et al., 2018). The confounding clinical factors, such as age, grade, gender, histology and treatment regimens were included as well. Clinical characteristics of these datasets in the OSgbm were presented in Table 1. All of the 684 patients have OS data, and the median OS time was 13.44 months, while 153 GBM patients from TCGA cohort have four above mentioned survival terms (OS, DSS, DFI and PFI). The median age of all the patients is 50 years. The death rate is 78.49%. A large proportion of the patients are in grade IV, especially in the two CGGA datasets (99.28% and 100%, respectively).
Table 1
| Data Source | Sample Size (n) | Median Age (years) | Death (%) | OS Median (years) | Gender (male, %) | Grade (I/II/III/IV, %) | Survival Terms |
|---|---|---|---|---|---|---|---|
| TCGA | 153 | 60 | 79.08 | 11.90 | 64.71 | – | OS, DSS, DFI, PFI |
| GSE7696 | 80 | 52 | 81.25 | 15.58 | 73.75 | – | OS |
| GSE4412 | 85 | 42 | 69.41 | 12.97 | 37.65 | 0/0/30.59/69.41 | OS |
| GSE42669 | 57 | 51 | 80.70 | 14.93 | 52.63 | – | OS |
| GSE30472 | 54 | – | 88.89 | 15.72 | – | 3.7/12.96/29.63/53.71 | OS |
| CGGAseq | 128 | 48 | 66.67 | 9.55 | 65.22 | 0/0/0.72/99.28 | OS |
| CGGAarray | 127 | 47 | 83.46 | 13.43 | 62.20 | 0/0/0/100 | OS |
| Total | 684 | 50 | 78.49 | 13.44 | 59.36 | – | – |
Clinical characteristics of each GBM dataset used in OSgbm.
Set-Up of OSgbm Web Server
The main function of OSgbm web server is to evaluate and determine the prognostic value of the quested genes. The users start by typing the gene symbol and choosing one dataset of interest or the combined dataset with pooling all the datasets together. To measure the association between a quested gene and survival, GBM samples are categorized according to the median (or other appropriate cutoff value, such as Trichotomy, Quartile) of the selected gene, and KM analysis is used to compare the outcomes between groups (Xie et al., 2019b). The user could limit the analysis in a subgroup of the patients by setting the age range, grade, gender and so on. Once the gene symbol is input and clinical characters are chosen, OS, DSS, DFI or PFI of each stratified group can be measured and analysis results will be available on the output web page. The prognostic value of each given gene is determined by HR (95% CI) and log-rank p value.
Validation of Previously Reported GBM Prognostic Biomarkers
To determine the performance of this online tool, 53 previously published GBM prognostic factors collected as the procedure shown in Figure S1 and then they were evaluated in OSgbm (Table 2, Figure 1) (Sano et al., 1999; Hung and Howng, 2003; Heimberger et al., 2005; Aaberg-Jessen et al., 2009; Shirai et al., 2009; Pu et al., 2011; Cui et al., 2013; De Tayrac et al., 2013; Lee et al., 2013; Rosati et al., 2013; Wang et al., 2013; Wu et al., 2013; Bai et al., 2014; Han et al., 2014; Meng et al., 2014; Olmez et al., 2014; Zupancic et al., 2014; Maris et al., 2015; Moutal et al., 2015; Yan et al., 2015; Zhao et al., 2015; Chaurasia et al., 2016; Nduom et al., 2016; Steponaitis et al., 2016; Steponaitis et al., 2016; Zhang et al., 2016a; Zhang et al., 2016b; Deng et al., 2017; Dong et al., 2017; Ge et al., 2017; Han et al., 2017; Haynes et al., 2017; Huang et al., 2017; Li et al., 2017; Luo and Zhuang, 2017; Ma et al., 2017; Mu et al., 2017; Roy et al., 2017; Wu et al., 2017; Zhang et al., 2017; Cai et al., 2018; Cetin et al., 2018; Cheng et al., 2018; Chou et al., 2018; Dahlrot et al., 2018; Gonçalves et al., 2018; Hayashi et al., 2018; Lv et al., 2018; Qian et al., 2018; Vasaikar et al., 2018; Hallal et al., 2019). OS was selected as the survival term. Among these prognostic genes, 51 of them showed significant prognostic ability in a large-scale combined cohort (33 genes) or in single cohort (18 genes), which were consistent with the prognostic value reported in the literature. The remaining two genes (IGF1R and PCBP2) display significant prognostic values in OSgbm, but is contradictory to what was reported in the literatures. Both of them were shown as favorable prognostic biomarkers in OSgbm but were reported to be unfavorable GBM prognostic biomarkers in previous reports (Table 2) (Maris et al., 2015; Luo and Zhuang, 2017).
Table 2
| Gene symbol | Validation results | Literature data | |||||||
|---|---|---|---|---|---|---|---|---|---|
| OS, HR (95% CI) | p Value | Cut Off | Osgbm | OS, HR (95% CI) | p Value | Sample (n) | Level | Reference | |
| PROM1 | 2.412 (1.040–4.174) | 0.007 | Upper 25% vs Lower 25% | GSE7679 | 2.39 (1.77–3.23) | <0.001 | 656 | mRNA | (Zhang et al., 2016) |
| SRGN | 2.371 (1.256–4.477) | 0.008 | Upper 25% vs Lower 25% | CGGAseq | – | 0.037 | 504 | mRNA | (Roy et al., 2017) |
| EDNRB | 2.272 (1.115–4.627) | 0.024 | Upper 25% vs Lower 75% | GSE30472 | 2.86 (1.12–7.34) | 0.031 | 25 | Protein | (Vasaikar et al., 2018) |
| PSMB4 | 2.074 (1.187–3.626) | 0.010 | Upper 25% vs Lower 25% | CGGAseq | – | <0.001 | 77 | Protein | (Cheng et al., 2018) |
| WNT6 | 2.035 (1.098–3.770) | 0.024 | Upper 25% vs Lower 25% | CGGAseq | – | 0.004 | 16 | Protein | (Gonçalves et al., 2018) |
| DPYSL5 | 2.023 (1.160–3.527) | 0.013 | Upper 25% vs Lower 25% | CGGAarray | – | 0.026 | 183 | Protein | (Moutal et al., 2015) |
| IL17A | 2.009 (1.107–3.646) | 0.022 | Upper 50% vs Lower 50% | GSE30472 | – | 0.007 | 41 | Protein | (Cui et al., 2013) |
| TLR9 | 1.976 (1.089–3.588) | 0.025 | Upper 25% vs Lower 25% | CGGAseq | – | 0.020 | 46 | Protein | (Mu et al., 2017) |
| ACKR3 | 1.974 (1.040–3.747) | 0.038 | Upper 30% vs Lower 30% | GSE7679 | 1.56 (1.04–2.51) | 0.03 | 146 | Protein | (Deng et al., 2017) |
| H19 | 1.864 (1.309–2.653) | <0.001 | Upper 25% vs Lower 25% | Combined | – | 0.034 | – | mRNA | (Wu et al., 2017) |
| EGFR | 1.845 (1.077–3.160) | 0.026 | Upper 25% vs Lower 75% | GSE7696 | – | <0.001 | 196 | Protein | (Heimberger et al., 2005) |
| NUSAP1 | 1.748 (1.006–3.040 | 0.048 | Upper 25% vs Lower 25% | CGGAarray | 0.65 (0.49–0.86)* | 0.003 | 518 | mRNA | (Qian et al., 2018) |
| CHAF1B | 1.707 (1.323–2.203) | <0.001 | Upper 30% vs Lower 30% | Combined | – | 0.004 | 96 | Protein | (De Tayrac et al., 2013) |
| TAGLN2 | 1.665 (1.282–2.161) | <0.001 | Upper 25% vs Lower 25% | Combined | – | <0.05 | 667 | mRNA | (Han et al., 2017) |
| BIRC1 | 1.658 (1.266–2.172) | <0.001 | Upper 25% vs Lower 25% | Combined | – | 0.0003 | 66 | Protein | (Shirai et al., 2009) |
| MGMT | 1.633 (1.260–2.115) | <0.001 | Upper 25% vs Lower 25% | Combined | 1.50 | 0.01 | 157 | Protein | (Dahlrot et al., 2018) |
| CD70 | 1.561 (1.180–2.065) | 0.002 | Upper 25% vs Lower 25% | Combined | 1.6 (0.98–2.51) | 0.046 | 107 | mRNA | (Ge et al., 2017) |
| CXCR4 | 1.558 (1.207–2.010) | <0.001 | Upper 25% vs Lower 25% | Combined | – | <0.05 | 156 | mRNA | (Ma et al., 2017) |
| CA9 | 1.556 (1.202–2.015) | <0.001 | Upper 25% vs Lower 75% | Combined | – | 0.004 | 66 | Protein | (Cetin et al., 2018) |
| PDCD1 | 1.508 (1.171–1.942) | 0.002 | Upper 30% vs Lower 30% | Combined | – | 0.028 | 149 | mRNA | (Nduom et al., 2016) |
| IDH1 | 1.490 (1.013–2.192) | 0.043 | Upper 50% vs Lower 50% | CGGAarray | – | 0.045 | 163 | Protein | (Chaurasia et al., 2016) |
| IGFBP2 | 1.467 (1.132–1.902) | 0.004 | Upper 25% vs Lower 25% | Combined | 1.04 (1.02–1.05) | 0.001 | 83 | Plasma | (Han et al., 2014) |
| PBK | 1.456 (1.131–1.875) | 0.004 | Upper 25% vs Lower 25% | Combined | – | 0.007 | 32 | Protein | (Hayashi et al., 2018) |
| EFEMP2 | 1.446 (1.117–1.871) | 0.005 | Upper 25% vs Lower 25% | Combined | – | <0.01 | 77 | mRNA | (Li et al., 2017) |
| MET | 1.434 (1.130–1.820) | 0.003 | Upper 30% vs Lower 30% | Combined | 1.7 (1.1–2.2) | <0.05 | 69 | Protein | (Olmez et al., 2014) |
| CHI3L1 | 1.438 (1.104–1.872) | 0.007 | Upper 25% vs Lower 25% | GSE30472 | – | <0.01 | 98 | mRNA | (Steponaitis et al., 2016) |
| TRAF2 | 1.443 (1.118–1.863) | 0.005 | Upper 25% vs Lower 25% | Combined | – | 0.03 | 105 | mRNA | (Zhang et al., 2017) |
| HMGB2 | 1.391 (1.099–1.759) | 0.006 | Upper 30% vs Lower 30% | Combined | 3.35 (1.25–9.02) | 0.017 | 51 | Protein | (Wu et al., 2013) |
| MCM6 | 1.387 (1.132–1.699) | 0.002 | Upper 25% vs Lower 75% | Combined | 1.19 | 0.006 | 325 | mRNA | (Cai et al., 2018) |
| CD44 | 1.386 (1.073–1.790) | 0.012 | Upper 25% vs Lower 25% | Combined | – | <0.001 | 28 | Protein | (Steponaitis et al., 2016) |
| TIMP1 | 1.342 (1.025–1.758) | 0.033 | Upper 25% vs Lower 25% | Combined | 3.2 (1.5–6.7) | 0.004 | 112 | Protein | (Aaberg-Jessen et al., 2009) |
| CD151 | 1.336 (1.023–1.746) | 0.034 | Upper 25% vs Lower 25% | Combined | 5.064 (1.427–17.969) | 0.012 | 211 | Protein | (Lee et al., 2013) |
| TWIST1 | 1.312 (1.013–1.699) | 0.039 | Upper 25% vs Lower 25% | Combined | 5.745 (1.331–1.89) | 0.017 | 86 | Protein | (Wang et al., 2013) |
| CCT6A | 1.316 (1.045–1.655) | 0.019 | Upper 30% vs Lower 30% | Combined | 3.21 (2.85–3.65) | 0.006 | 497 | Protein | (Hallal et al., 2019) |
| APC | 1.308 (1.093–1.566) | 0.004 | Upper 50% vs Lower 50% | Combined | – | <0.001 | 83 | Protein | (Rosati et al., 2013) |
| CD247 | 1.292 (1.022–1.633) | 0.032 | Upper 30% vs Lower 30% | Combined | 1.54 (1.05–2.28) | 0.023 | 149 | mRNA | (Nduom et al., 2016) |
| CXCR3 | 1.272 (1.027–1.575) | 0.028 | Upper 25% vs Lower 75% | Combined | 1.56 (1.04–2.51) | 0.03 | 146 | Protein | (Pu et al., 2011) |
| TCTN1 | 1.223 (1.011–1.493) | 0.039 | Upper 30% vs Lower 70% | Combined | 1.32 (1.08–1.61) | 0.006 | 518 | mRNA | (Meng et al., 2014) |
| BICD1 | 0.794 (0.644–0.978)# | 0.030 | Lower 25% vs Upper 75% | Combined | 1.577 (1.299–1.914) | <0.001 | 523 | mRNA | (Huang et al., 2017) |
| IFIT1 | 0.770 (0.609–0.973) | 0.029 | Upper 30% vs Lower 30% | Combined | 0.22 (0.10–0.52) | 0.001 | 70 | mRNA | (Zhang et al., 2016) |
| BRMS1L | 0.753 (0.587–0.966) | 0.026 | Upper 25% vs Lower 75% | Combined | – | <0.05 | 60 | mRNA | (Lv et al., 2018) |
| IGF1R | 0.745 (0.588–0.944) | 0.015 | Upper 30% vs Lower 30% | Combined | 1.65 (1.10–2.47) | 0.016 | 167 | Protein | (Maris et al., 2015) |
| GANO1 | 0.748 (0.585–0.957) | 0.021 | Upper 30% vs Lower 30% | Combined | – | 0.009 | 178 | Protein | (Zupancic et al., 2014) |
| PTEN | 0.729 (0.567–0.938) | 0.014 | Upper 25% vs Lower 25% | Combined | 3.3 (1.6–4.3)* | 0.0003 | 61 | mRNA | (Sano et al., 1999) |
| SEMA6A | 0.694 (0.556–0.867) | 0.001 | Upper 25% vs Lower 75% | Combined | 1.71 (1.01–2.65)* | 0.012 | 200 | Protein | (Zhao et al., 2015) |
| PHF3 | 0.683 (0.529–0.883) | 0.004 | Upper 25% vs Lower 25% | Combined | 0.44 (0.26–0.77) | 0.0031 | 35 | Protein | (Yan et al., 2015) |
| PPARα | 0.644 (0.503–0.825) | <0.001 | Upper 30% vs Lower 30% | Combined | 1.31 (1.05–1.63)* | 0.016 | 473 | mRNA | (Haynes et al., 2017) |
| PCBP2 | 0.632 (0.417–0.957) | 0.031 | Upper 25% vs Lower 75% | TCGA | – | <0.001 | 130 | mRNA | (Luo and Zhuang, 2017) |
| LAPTM4B | 0.626 (0.433–0.894) | 0.010 | Upper 50% vs Lower 50% | TCGA | – | <0.001 | 39 | Protein | (Dong et al., 2017) |
| ANXA7 | 0.619 (0.475–0.806) | <0.001 | Upper 25% vs Lower 25% | Combined | – | <0.001 | 99 | Protein | (Hung and Howng, 2003) |
| PHF20 | 0.557 (0.319–0.972) | 0.040 | Upper 50% vs Lower 50% | CGGAarray | 0.5 (0.29–0.86) | 0.012 | 62 | Protein | (Yan et al., 2015) |
| TES | 0.407 (0.173–0.958) | 0.040 | Upper 30% vs Lower 30% | GSE42669 | – | <0.05 | 37 | Protein | (Bai et al., 2014) |
| LGALS1 | 0.368 (0.157–0.863) | 0.022 | Upper 25% vs Lower 25% | GSE42669 | – | 0.009 | 45 | Protein | (Chou et al., 2018) |
Validation of previously reported prognostic biomarkers in OSgbm.
*: The lower gene expression compared with higher gene expression in the literature data.
#: The lower gene expression compared with higher gene expression in the OSgbm data.
Figure 1
Discussion
The development of prognostic biomarkers is important for guiding the treatments especially for therapy-resistant GBM patients. In our work, we developed a new web server, OSgbm, to help researchers to evaluate the prognostic value of a given gene for GBM patients. OSgbm is easy to use and requires no special skills (such as bioinformatics training). With filtering by one or several clinical confounding factors provided in OSgbm, users can also evaluate the prognostic value of their interested genes according to their special needs. The function and performance tests of OSgbm web server showed that 96% (51 out of 53) of previously reported prognostic biomarkers could be confirmed in OSgbm, which indicates that these biomarkers validated in independent cohorts have the potency of translating to clinical applications, and also indicates the well performance of OSgbm. Nevertheless, there are two genes including IGF1R and PCBP2 which showed different prognostic values to the literatures, the discrepancy of prognostic performance of IGF1R and PCBP2 between OSgbm and literatures may be caused by race, different cohort size, or analysis level and methods (mRNA vs. protein, gene microarray vs. immunohistochemistry) (Maris et al., 2015; Luo and Zhuang, 2017). For example, the race reported in literatures for PCBP2 is Asian, while that in validated cohort of OSgbm is mostly White. The mRNA level was analyzed in OSgbm for IGF1R, while IGF1R was determined by immunohistochemistry in literature. In addition, the race analyzed in OSgbm for IGF1R is Asian (Korea for GSE42669 and Chinese for CGGA), while the race reported in literature for IGF1R is European. As a result, it will be necessary to validate the prognostic performance of IGF1R and PCBP2 in a larger independent cohort of glioblastoma.
In conclusion, OSgbm is a user-friendly web server to help researchers and clinicians to identify suitable prognostic biomarkers in GBM. Furthermore, we will keep update the database of OSgbm to collect more and more GBM datasets when new GBM dataset is available, and will implement the multivariate cox proportional hazards model into OSgbm for the purpose of adjustment for the confounding clinical factors, and we also encourage users to contact us to upload their own data into OSgbm.
Statements
Data availability statement
All datasets for this study are included in the article/Supplementary Material.
Author contributions
XG conceived and directed the project. HD and QW collected data and developed the web server. HD, NL, JL, LG, MY, GZ, YA, FW, LX, and YL performed data analysis. WZ, HZ, and MZ contributed to data analysis and paper writing. XG and HD wrote the manuscript with the assistance and approval of all authors.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (No. 81602362), the program for Science and Technology Development in Henan Province (No. 162102310391), the supporting grants of Henan University (No. 2015YBZR048; No. B2015151), the program for Innovative Talents of Science and Technology in Henan Province (No. 18HASTIT048), and Yellow River Scholar Program (No. H2016012), Kaifeng Science and Technology Major Project (No. 18ZD008).
Conflict of interest
Author MZ is employed by company of Nanjing Jiliang Biotechnology Co., Ltd.
The remaining 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.01378/full#supplementary-material
References
1
Aaberg-JessenC.ChristensenK.OffenbergH.BartelsA.DreehsenT.HansenS.et al. (2009). Low expression of tissue inhibitor of metalloproteinases-1 (TIMP-1) in glioblastoma predicts longer patient survival. J. Neurooncol.95, 117–128. doi: 10.1007/s11060-009-9910-8
2
AldapeK.ZadehG.MansouriS.ReifenbergerG.von DeimlingA. (2015). Glioblastoma: pathology, molecular mechanisms and markers. Acta Neuropathol.129, 829–848. doi: 10.1007/s00401-015-1432-1
3
BaiY.ZhangQ.WangX. (2014). Downregulation of TES by hypermethylation in glioblastoma reduces cell apoptosis and predicts poor clinical outcome. Eur. J. Med. Res.19, 66. doi: 10.1186/s40001-014-0066-4
4
BrownD. V.FilizG.DanielP. M.HollandeF.DworkinS.AmiridisS.et al. (2017). Expression of CD133 and CD44 in glioblastoma stem cells correlates with cell proliferation, phenotype stability and intra-tumor heterogeneity. PloS One12, e0172791. doi: 10.1371/journal.pone.0172791
5
CaiX.SughrueM. E. (2017). Glioblastoma: new therapeutic strategies to address cellular and genomic complexity. Oncotarget20, 9540–9554. doi: 10.18632/oncotarget.23476
6
CaiH.ChengZ.ZhangH.WangP.ZhangY.HaoJ.et al. (2018). Overexpression of MCM6 predicts poor survival in patients with glioma. Hum. Pathol.78, 182–187. doi: 10.1016/j.humpath.2018.04.024
7
CetinB.GonulI. I.GumusayO.BilgetekinI.AlginE.OzetA.et al. (2018). Carbonic anhydrase IX is a prognostic biomarker in glioblastoma multiforme. Neuropathology38, 457–462. doi: 10.1111/neup.12485.
8
ChaurasiaA.ParkS. H.SeoJ. W.ParkC. K. (2016). Immunohistochemical analysis of ATRX, IDH1 and p53 in glioblastoma and their correlations with patient survival. J. Korean Med. Sci.31, 1208–1214. doi: 10.3346/jkms.2016.31.8.1208
9
ChengY.TsaiW.SungY.ChangH.ChenY. (2018). Interference with PSMB4 expression exerts an anti-tumor effect by decreasing the invasion and proliferation of human glioblastoma cells. Cell Physiol. Biochem.45, 819–831. doi: 10.1159/000487174
10
ChouS.YenS.HuangC.HuangE. (2018). Galectin-1 is a poor prognostic factor in patients with glioblastoma multiforme after radiotherapy. BMC Cancer18, 105. doi: 10.1186/s12885-018-4025-2
11
CuiX.XuZ.ZhaoZ.SuiD.RenX.HuangQ.et al. (2013). Analysis of CD137L and IL-17 expression in tumor tissue as prognostic indicators for gliblastoma. Int. J. Biol. Sci.9, 134–141. doi: 10.7150/ijbs.4891
12
DahlrotR. H.DowsettJ.FosmarkS.MalmströmA.HenrikssonR.BoldtH.et al. (2018). Prognostic value of O-6-methylguanine-DNA methyltransferase (MGMT) protein expression in glioblastoma excluding nontumour cells from the analysis. Neuropathol. Appl. Neurobiol.44, 172–184. doi: 10.1111/nan.12415
13
De TayracM.SaikaliS.AubryM.BellaudP.BonifaceR.QuillienV.et al. (2013). Prognostic significance of EDN/RB, HJURP, p60/CAF-1 and PDLI4, four new markers in high-grade gliomas. PloS One8, e73332. doi: 10.1371/journal.pone.0073332
14
DengL.ZhengW.DongX.LiuJ.ZhuC.LuD.et al. (2017). Chemokine receptor CXCR7 is an independent prognostic biomarker in glioblastoma. Cancer Biomark.20, 1–6. doi: 10.3233/CBM-151430
15
DongX.TamuraK.KobayashiD.AndoN.SumitaK.MaeharaT. (2017). LAPTM4B-35 is a novel prognostic factor for glioblastoma. J. Neurooncol.132, 295–303. doi: 10.1007/s11060-017-2369-0
16
GeH.MuL.JinL.YangC.ChangY.LongY.et al. (2017). Tumor associated CD70 expression is involved in promoting tumor migration and macrophage infiltration in GBM. Int. J. Cancer141, 1434–1444. doi: 10.1002/ijc.30830
17
GonçalvesC. S.Vieira de CastroJ.PojoM.MartinsE. P.QueirósS.ChautardE.et al. (2018). WNT6 is a novel oncogenic prognostic biomarker in human glioblastoma. Theranostics8, 4805–4823. doi: 10.7150/thno.25025
18
HallalS.RussellB. P.WeiH.LeeM. Y. T.ToonC. W.SyJ.et al. (2019). Extracellular Vesicles from neurosurgical aspirates identifies chaperonin containing TCP1 subunit 6A as a potential glioblastoma biomarker with prognostic significance. Proteomics19, e1800157. doi: 10.1002/pmic.201800157
19
HanS.MengL.HanS.WangY.WuA. (2014). Plasma IGFBP-2 levels after postoperative combined radiotherapy and chemotherapy predict prognosis in elderly glioblastoma patients. PloS One9, e93791. doi: 10.1371/journal.pone.0093791
20
HanM.XuR.XuY.ZhangX.NiS.HuangB.et al. (2017). TAGLN2 is a candidate prognostic biomarker promoting tumorigenesis in human gliomas. J. Exp. Clin. Cancer Res.36, 155. doi: 10.1186/s13046-017-0619-9
21
HayashiT.HayakawaY.KohM.TomitaT.NagaiS.KashiwazakiD.et al. (2018). Impact of a novel biomarker, T-LAK cell-originating protein kinase (TOPK) expression on outcome in malignant glioma. Neuropathology38, 144–153. doi: 10.1111/neup.12446
22
HaynesH. R.WhiteP.HaresK. M.RedondoJ.KempK. C.SingletonW. G. B.et al. (2017). The transcription factor PPARα is overexpressed and is associated with a favourable prognosis in IDH-wildtype primary glioblastoma. Histopathology70, 1030–1043. doi: 10.1111/his.13142
23
HeimbergerA. B.HlatkyR.SukiD.YangD.WeinbergJ.GilbertM.et al. (2005). Prognostic effect of epidermal growth factor receptor and EGFRvIII in glioblastoma multiforme patients. Clin. Cancer Res.11, 1462–1466. doi: 10.1158/1078-0432.CCR-04-1737
24
HuangS. P.ChangY. C.LowQ. H.WuA. T. H.ChenC. L.LinY. F.et al. (2017). BICD1 expression, as a potential biomarker for prognosis and predicting response to therapy in patients with glioblastomas. Oncotarget8, 113766–113791. doi: 10.18632/oncotarget.22667
25
HungK. S.HowngS. L. (2003). Prognostic significance of annexin VII expression in glioblastomas multiforme in humans. J. Neurosurg.99, 886–892. doi: 10.3171/jns.2003.99.5.0886
26
JohnsonD. R.O'NeillB. P. (2012). Glioblastoma survival in the United States before and during the temozolomide era. J. Neurooncol.107, 359–364. doi: 10.1002/cncr.26494
27
LeeD.SuhY. L.ParkT. I.DoI. G.SeolH. J.NamD. H.et al. (2013). Prognostic significance of tetraspanin CD151 in newly diagnosed glioblastomas. J. Surg. Oncol.107, 646–652. doi: 10.1002/jso.23249
28
LiF.LiY.ZhangK.LiY.HeP.LiuY.et al. (2017). FBLN4 as candidate gene associated with long-term and short-term survival with primary glioblastoma. Onco. Targets Ther.10, 387–395. doi: 10.2147/OTT.S117165
29
LiuJ.LichtenbergT.HoadleyK. A.PoissonL. M.LazarA. J.CherniackA. D.et al. (2018). An integrated tcga pan-cancer clinical data resource to drive high-quality survival outcome analytics. Cell173, 400–416. doi: 10.1016/j.cell.2018.02.052
30
LuoK.ZhuangK. (2017). High expression of PCBP2 is associated with progression and poor prognosis in patients with glioblastoma. BioMed. Pharmacother.94, 659–665. doi: 10.1016/j.biopha.2017.07.103
31
LvJ.YangH.WangX.HeR.DingL.SunX. (2018). Decreased BRMS1L expression is correlated with glioma grade and predicts poor survival in glioblastoma via an invasive phenotype. Cancer Biomark.22, 311–316. doi: 10.3233/CBM-171019
32
MaX.ShangF.ZhuW.LinQ. (2017). CXCR4 expression varies significantly among different subtypes of glioblastoma multiforme (GBM) and its low expression or hypermethylation might predict favorable overall survival. Expert Rev. Neurother7, 941–946. doi: 10.1080/14737175.2017.1351299
33
MarisC.D'HaeneN.TrépantA. L.Le MercierM.SauvageS.AllardJ. (2015). IGF-IR: a new prognostic biomarker for human glioblastoma. Br. J. Cancer113, 729–737. doi: 10.1038/bjc.2015.242
34
MengD.ChenY.ZhaoY.WangJ.YunD.YangS.et al. (2014). Expression and prognostic significance of TCTN1 in human glioblastoma. J. Transl. Med.12, 288. doi: 10.1186/s12967-014-0288-9
35
MoutalA.HonnoratJ.MassomaP.DésormeauxP.BertrandC.MallevalC.et al. (2015). CRMP5 controls glioblastoma cell proliferation and survival through notch-dependent signaling. Cancer Res.75, 3519–3528. doi: 10.1158/0008-5472.CAN-14-0631
36
MuL.WangY.WangY.ZhangH.ShangD.TanF.et al. (2017). Tumor location and survival outcomes in adult patients with supratentorial glioblastoma by levels of toll-like receptor 9 expression. World Neurosurg.97, 279–283. doi: 10.1016/j.wneu.2016.10.015
37
NathansonD. A.GiniB.MottahedehJ.VisnyeiK.KogaT.GomezG.et al. (2014). Targeted therapy resistance mediated by dynamic regulation of extrachromosomal mutant EGFR DNA. Science343, 72–76. doi: 10.1126/science.1241328
38
NduomE. K.WeiJ.YaghiN. K.HuangN.KongL. Y.GabrusiewiczK.et al. (2016). PD-L1 expression and prognostic impact in glioblastoma. Neuro Oncol.18, 195–205. doi: 10.1093/neuonc/nov172
39
NguyenH. S.ShabaniS.AwadA. J.KaushalM.DoanN. (2018). Molecular markers of therapy-resistant glioblastoma and potential strategy to combat resistance. Int. J. Mol. Sci.19. doi: 10.3390/ijms19061765
40
NikiforovaM. N.HamiltonR. L. (2011). Molecular diagnostics of gliomas. Arch. Pathol. Lab. Med.135, 558–568. doi: 10.1043/2010-0649-RAIR.1
41
OlmezO. F.CubukcuE.EvrenselT.KurtM.AvciN.TolunayS.et al. (2014). The immunohistochemical expression of c-Met is an independent predictor of survival in patients with glioblastoma multiforme. Clin. Transl. Oncol.16, 173–177. doi: 10.1007/s12094-013-1059-4
42
PolivkaJJrPolivkaJ.HolubecL.KubikovaT.PribanV.HesO.et al. (2017). Advances in experimental targeted therapy and immunotherapy for patients with glioblastoma multiforme. Anticancer Res.37, 21–33. doi: 10.21873/anticanres.11285
43
PuY.LiS.ZhangC.BaoZ.YangZ.SunL. (2011). High expression of CXCR3 is an independent prognostic factor in glioblastoma patients that promotes an invasive phenotype. J. Neurooncol.122, 43–51. doi: 10.1007/s11060-014-1692-y
44
QianZ.LiY.MaJ.XueY.XiY.et al. (2018). Prognostic value of NUSAP1 in progression and expansion of glioblastoma multiforme. J. Neurooncol.140, 199–208. doi: 10.1007/s11060-018-2942-1
45
RosatiA.PolianiP. L.TodeschiniA.CominelliM.MedicinaD.CenzatoM.et al. (2013). Glutamine synthetase expression as a valuable marker of epilepsy and longer survival in newly diagnosed glioblastoma multiforme. Neuro Oncol.15, 618–625. doi: 10.1093/neuonc/nos338.
46
RoyA.AttarhaS.WeishauptH.EdqvistP. H.SwartlingF. J.BergqvistM.et al. (2017). Serglycin as a potential biomarker for glioma: association of serglycin expression, extent of mast cell recruitment and glioblastoma progression. Oncotarget11, 24815–24827. doi: 10.18632/oncotarget.15820
47
SanoT.LinH.ChenX.LangfordL. A.KoulD.BondyM. L.et al. (1999). Differential expression of MMAC/PTEN in glioblastoma multiforme: relationship to localization and prognosis. Cancer Res.59, 1820–1824.
48
ShiraiK.SuzukiY.OkaK.NodaS. E.KatohH.SuzukiY.et al. (2009). Nuclear survivin expression predicts poorer prognosis in glioblastoma. J. Neurooncol.91, 353–358. doi: 10.1007/s11060-008-9720-4
49
SteponaitisG.SkiriutėD.KazlauskasA.GolubickaitėI.StakaitisR.TamašauskasA.et al. (2016). High CHI3L1 expression is associated with glioma patient survival. Diagn. Pathol.11, 42. doi: 10.1186/s13000-016-0492-4
50
StoyanovG. S.DzhenkovD.GhenevP.IlievB.EnchevY.TonchevA. B. (2018). Cell biology of glioblastoma multiforme: from basic science to diagnosis and treatment. Med. Oncol.35, 27. doi: 10.1007/s12032-018-1083-x
51
StuppR.MasonW. P.van den BentM. J.WellerM.FisherB.TaphoornM. J.et al. (2005). Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N Engl. J. Med.352, 987–996. doi: 10.1056/NEJMoa043330
52
VasaikarS.TsiprasG.LandázuriN.CostaH.WilhelmiV.SciclunaP.et al. (2018). Overexpression of endothelin B receptor in glioblastoma: a prognostic marker and therapeutic target?BMC Cancer18, 154. doi: 10.1074/jbc.M110.178012
53
WangS.KeY.LuG.SongZ.YuL.XiaoS.et al. (2013). Vasculogenic mimicry is a prognostic factor for postoperative survival in patients with glioblastoma. J. Neurooncol.112, 339–345. doi: 10.1007/s11060-013-1077-7
54
WangQ.XieL.DangY.SunX.XieT.GuoJ.et al. (2019). OSlms: A web server to evaluate the prognostic value of genes in Leiomyosarcoma. Front. Oncol.9, 190. doi: 10.3389/fonc.2019.00190
55
WangQ.WangF.LvJ.XinJ.XieL.ZhuW.et al. (2019). Interactive online consensus survival tool for esophageal squamous cell carcinoma prognosis analysis. Oncol. Lett.18, 1199–1206. doi: 10.3892/ol.2019.10440
56
WuZ.CaiL.LinS.XiongZ.LuJ.MaoY.et al. (2013). High-mobility group box 2 is associated with prognosis of glioblastoma by promoting cell viability, invasion, and chemotherapeutic resistance. Neuro Oncol.15, 1264–1275. doi: 10.1093/neuonc/not078
57
WuW.HuQ.NieE.YuT.WuY.ZhiT.et al. (2017). Hypoxia induces H19 expression through direct and indirect Hif-1α activity, promoting oncogenic effects in glioblastoma. Sci. Rep.7, 45029. doi: 10.1038/srep45029
58
XieL.WangQ.DangY.GeL.SunX.LiN.et al. (2019a). OSkirc: a web tool for identifying prognostic biomarkers in kidney renal clear cell carcinoma. Future Oncol15, 3103–3110. doi: 10.2217/fon-2019-0296
59
XieL.DangY.GuoJ.SunX.XieT.ZhangL.et al. (2019b). High KRT8 expression independently predicts poor prognosis for lung adenocarcinoma patients. Genes (Basel)10. doi: 10.3390/genes10010036
60
YanJ.KongL. Y.HuJ.GabrusiewiczK.DibraD.XiaX.et al. (2015). Autoantibodies against GLEA2 and PHF3 in glioblastoma: tumor-associated autoantibodies correlated with prolonged survival. J. Natl. Cancer Inst.107 (8). doi: 10.1002/ijc.20929
61
YangP.CaiJ.YanW.ZhangW.WangY.ChenB.et al. (2016). Classification based on mutations of TERT promoter and IDH characterizes subtypes in grade II/III gliomas. Neuro Oncol.18, 1099–1108. doi: 10.1093/neuonc/now021
62
ZhangW.ChenH.LvS.YangH. (2016a). High CD133 expression is associated with worse prognosis in patients with glioblastoma. Mol. Neurobiol.53, 2354–2360. doi: 10.1007/s12035-015-9187-1
63
ZhangJ.ChenY.LinG.ZhangJ.TangW.HuangJ.et al. (2016b). High IFIT1 expression predicts improved clinical outcome, and IFIT1 along with MGMT more accurately predicts prognosis in newly diagnosed glioblastoma. Hum. Pathol.52, 136–144. doi: 10.1016/j.humpath.2016.01.013
64
ZhangW.SunY.LiuL.LiZ. (2017). Prognostic significance of TNFR-associated factor 1 and 2 (TRAF1 and TRAF2) in glioblastoma23,4506–4512. Med. Sci. Monit.23, 4506–4512. doi: 10.12659/msm.903397
65
ZhangG.WangQ.YangM.YuanQ.DangY.SunX.et al. (2019). OSblca: A web server for investigating prognostic biomarkers of bladder cancer patients. Front. Oncol.9, 466. doi: 10.3389/fonc.2019.00466
66
ZhaoJ.TangH.ZhaoH.CheW.ZhangL.LiangP. (2015). SEMA6A is a prognostic biomarker in glioblastoma. Tumour Biol.36, 8333–8340. doi: 10.1007/s13277-015-3584-y
67
ZupancicK.BlejecA.HermanA.VeberM.VerbovsekU.KorsicM.et al. (2014). Identification of plasma biomarker candidates in glioblastoma using an antibody-array-based proteomic approach. Radiol. Oncol.48, 257–266. doi: 10.2478/raon-2014-0014
Summary
Keywords
glioblastoma, survival analysis, prognostic biomarker, OSgbm, transcriptome profiles, clinical information
Citation
Dong H, Wang Q, Li N, Lv J, Ge L, Yang M, Zhang G, An Y, Wang F, Xie L, Li Y, Zhu W, Zhang H, Zhang M and Guo X (2020) OSgbm: An Online Consensus Survival Analysis Web Server for Glioblastoma. Front. Genet. 10:1378. doi: 10.3389/fgene.2019.01378
Received
08 May 2019
Accepted
17 December 2019
Published
21 February 2020
Volume
10 - 2019
Edited by
Meng Zhou, Wenzhou Medical University, China
Reviewed by
Zhixiang Zuo, Sun Yat-sen University, China; Xuexin Yu, UT Southwestern Medical Center, United States
Updates
Copyright
© 2020 Dong, Wang, Li, Lv, Ge, Yang, Zhang, An, Wang, Xie, Li, Zhu, Zhang, Zhang and Guo.
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: Xiangqian Guo, xqguo@henu.edu.cn
†These authors have contributed equally to this work
This article was submitted to Cancer Genetics, a section of the journal Frontiers in Genetics
Disclaimer
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.