A Multielement Prognostic Nomogram Based on a Peripheral Blood Test, Conventional MRI and Clinical Factors for Glioblastoma

Abstract

Background:

Glioblastoma (GBM) is one of the most malignant types of tumors in the central nervous system, and the 5-year survival remains low. Several studies have shown that preoperative peripheral blood tests and preoperative conventional Magnetic Resonance Imaging (MRI) examinations affect the prognosis of GBM patients. Therefore, it is necessary to construct a risk score based on a preoperative peripheral blood test and conventional MRI and develop a multielement prognostic nomogram for GBM.

Methods:

This study retrospectively analyzed 131 GBM patients. Determination of the association between peripheral blood test variables and conventional MRI variables and prognosis was performed by univariate Cox regression. The nomogram model, which was internally validated using a cohort of 56 GBM patients, was constructed by multivariate Cox regression. RNA sequencing data from Gene Expression Omnibus (GEO) and Chinese Glioma Genome Atlas (CGGA datasets were used to determine peripheral blood test-related genes based on GBM prognosis.

Results:

The constructed risk score included the neutrophil/lymphocyte ratio (NLR), lymphocyte/monocyte ratio (LMR), albumin/fibrinogen (AFR), platelet/lymphocyte ratio (PLR), and center point–to-ventricle distance (CPVD). A final nomogram was developed using factors associated with prognosis, including age, sex, the extent of tumor resection, IDH mutation status, radiotherapy status, chemotherapy status, and risk. The Area Under Curve (AUC) values of the receiver operating characteristic curve (ROC) curve were 0.876 (12-month ROC), 0.834 (24-month ROC) and 0.803 (36-month ROC) in the training set and 0.906 (12-month ROC), 0.800 (18-month ROC) and 0.776 (24-month ROC) in the validation set. In addition, vascular endothelial growth factor A (VEGFA) was closely associated with NLR and LMR and identified as the most central negative gene related to the immune microenvironment and influencing immune activities.

Conclusion:

The risk score was established as an independent predictor of GBM prognosis, and the nomogram model exhibit appropriate predictive power. In addition, VEGFA is the key peripheral blood test-related gene that is significantly associated with poor prognosis.

Background

Glioblastomas (GBMs) are prevalent malignant tumors of the brain and the central nervous system (CNS), causing most deaths among primary brain tumor patients (1). Despite advances in cancer diagnosis and treatment, no significant improvements in patient outcomes have been reported. GBM patients have poor clinical outcomes with 5-year relative survival rates and median overall survival times of 5.6% and 12–15 months, respectively (2, 3). Therefore, it is necessary to determine appropriate prognostic factors for suitable risk classification. Studies have found several molecular markers, including the methylation state of the O6-methylguanine-DNA methyltransferase (MGMT) gene promoter and isocitrate dehydrogenase enzyme 1/2 (IDH1/2) mutation. However, these markers are not preoperatively available. Therefore, there is a need to develop a reliable and simple preoperative scoring system for the prediction of GBM prognosis (4–9).

Peripheral blood tests can significantly impact tumor growth (10). Peripheral blood tests, such as NLR, PLR, LMR, red cell volume distribution width (RDW), albumin/globulin ratio (AGR), albumin/gamma-glutamyl transferase ratio (AGgR) and fibrinogen, can refine patient stratification to therapy and predict survival outcomes in gliomas (11–14). Moreover, conventional MRI offers diagnostic value for genotype classification (15). We aimed to build a risk score based on peripheral blood tests and conventional MRI for GBM prognosis. A nomogram was also developed to predict individualized survival outcomes for GBM patients using this scoring system, age, sex, IDH mutation, and other important prognostic factors. RNA sequencing data from the Chinese Glioma Genome Atlas (CGGA) were used to determine key peripheral blood test-related genes based on GBM prognosis.

Materials and Methods

Study Population

This study retrospectively reviewed 233 patients pathologically diagnosed with WHO grade IV GBM at the First Affiliated Hospital of Wenzhou Medical University from March 2011 to March 2020. The inclusion criteria for this study were as follows: (1) Age ≥ 18 years. (2) Patients with first-time preoperative peripheral blood cell data for GBM samples. (3) Patients without acute conditions, such as microbial infections, or those patients who underwent drug treatments that could have affected the immune system. (4) Patients with complete clinical data as well as follow-up data. A total of 187 patients met our inclusion criteria and were included in our study. A total of 131 patients treated from March 2011 to September 2018 were included in the training set, whereas 56 patients identified from October 2018 to March 2020 were included in the validation set. This study followed the Declaration of Helsinki and was permitted by the review board of the First Affiliated Hospital of WMU. Patients were required to sign informed consent for the use of their clinical data in future studies.

RNA expression profiles of immune cell subsets of humans were retrieved from Gene GEO (https://www.ncbi.nlm.nih.gov/geo) datasets (GSE28491). Molecular pathology and RNA sequencing data as well as clinical information for 693 glioma patients were retrieved from the CGGA database (http://www.cgga.org.cn), and we also retrieved RNA sequencing data of GBM from The Cancer Genome Atlas (TCGA; http://tcga-data.nci.nih.gov/tcga/), GSE16011 (16). RNA sequencing data of normal tissues were acquired from the Genotype-Tissue Expression Project (GTEx) (http://www.gtexportal.org) and log2(X+1) transformed for subsequent analyses. Single-cell sequencing (ScRNA-seq) data of key genes were obtained from Cancer Single-cell State Atlas (CancerSEA) (http://biocc.hrbmu.edu.cn/CancerSEA/), Patel AP. Science. 2014 (Brain) (EXP0058) (GSE57872) (17).

Data Collection

The patients' clinical parameters, including sex, tumor location, age, tumor size, resection extent, postoperative adjuvant chemoradiotherapy, preoperative fibrinogen, preoperative neutrophil and lymphocyte counts, preoperative albumin and globulin levels, and preoperative gamma-glutamyl transferase, were collected from the electronic medical records system.

Two blinded neuroradiologists with ≥7 years of MR as well as CT neuroimaging experience independently reviewed MRI scans. The assessed imaging features were (1) ring-shaped gadolinium enhancement (RE) (18) (presence or absence), (2) peritumoural oedema (presence or absence), (3) involvement of the subventricular zone (SVZ) ion on T2WI (presence or absence), (4) involvement of the insula on T2WI (presence or absence), (5) center point–to-ventricle distance (CPVD), and (6) tumor volume (D1^*D2^*D3/2) (19– 21).

Statistical Analysis

A log-rank test was used to establish the optimal threshold values for NLR, LMR, PLR, CPVD, and AFR. Continuous variables are expressed as numbers and means ± SD and were compared with unpaired Student's t-tests. Categorical data were assessed using the chi-square test. Correlations were calculated using Spearman's rank correlation. Kaplan–Meier survival curves were used for survival analysis (22–24). Cox proportional hazards regression models were used for univariate and multivariate analyses of clinical variables for the determination of independent prognostic factors. Finally, differentially expressed genes were examined using the Wilcoxon test. Univariate Cox regression, LASSO regression and multivariate Cox regression were performed to select the key genes. p ≤ 0.05 denoted significance, and all p values were two-sided. SPSS 22.0 was used for all statistical analyses.

Establishment and Validation of the Nomogram

The risk score system was calculated as follows: GBM patients with high NLR (>4.16), high PLR (>193.75), low LMR (≤2.52), low AFR (≤9.43), and short CPVD (≤32.56) had a score of 5 (five abnormalities). A risk score of 0 was defined as low risk; risk scores 1 and 2 were defined as moderate risk; risk scores 3 and 4 were defined as high risk; and a risk score of 5 was defined as ultra-high risk. With the exception of the number of chemotherapies in the training set, the four risk groups did not differ in participant characteristics at baseline in either the training set or validation set. Then, independent prognostic factors, such as age, sex, extent of resection (EOR), IDH mutation, postoperative adjuvant chemoradiotherapy and risk score, were screened by multivariate Cox regression analysis and used to construct a nomogram to predict the 12-, 24-, and 36-month overall survival (OS) outcomes of GBM. A “nomogramEx” package was used to obtain each patient's point in the nomogram. Receiver operating characteristic (ROC) curves were drawn, and the area under the curve (AUC) was determined to test the predictive significance of the nomogram. Moreover, the nomogram calibration curve was plotted to approximate the predictive ability of the nomogram. The models were trained and validated using the training and validation set data, respectively.

Exploration of Peripheral Blood Test-Related Genes

The downloaded RNA expression profiles of human immune cell subsets (GSE28491, 5 neutrophil samples, 5 monocyte samples and 20 lymphocyte samples) were used to explore the peripheral blood test-related genes that also influence the immune infiltrates of GBM. The edgeR package in R was used for the analysis of differential expression between neutrophils and lymphocytes as well as monocytes and lymphocytes. Genes with p < 0.05 were determined to be differentially expressed genes. Signatures related to NLR and LMR were denoted as significantly differentially expressed genes between neutrophils and lymphocytes and between monocytes and lymphocytes. The differentially expressed genes were also screened between CGGA and GTEx, representing normal and GBM tissues. Finally, various immune-associated genes were obtained from ImmPort (www.immport.org/shared/genelists). Thirty-eight genes were expressed in the final intersection and were used for the subsequent analysis. Genes that were markedly associated with prognosis were assessed by univariate Cox regression analysis at a cut-off p value < 0.05. The training set was further subjected to LASSO regression analysis and multivariate Cox regression analysis to identify the best prognostic genes via the glmnet package in R (22). Based on GBM prognosis, the key genes (SLC11A1 and VEGF) were selected by multivariable Cox regression. Therefore, VEGFA was determined to be a key prognostic-associated immune-related gene (IRG) by combining univariate Cox regression and PPI network analyses with clinical factors. Finally, RNA sequencing data from CGGA, TCGA and GSE16011 were used to calculate the correlation between NLR, LMR and expression of VEGFA. This was achieved with the CIBERSORT algorithm based on the deconvolution, using the ‘CIBERSORT' R package. The thresholds for inclusion was P ≤ 0.05 (25).

Multivariate Analysis of VEGFA

Correlations between VEGFA expression and clinical factors, including sample type, gender, age, IDH mutation status,1p19q codeletion status and MGMT promoter methylation status were evaluated. GO and KEGG analyses of VEGFA were performed via gene set enrichment analysis (GSEA) (26). Single-sample gene set enrichment analyses (ssGSEA) were performed using “gsva” in R. Infiltration scores for 16 immune cells as well as the activities of 13 immune-associated pathways were determined for each sample in the 693 CGGA dataset. Subsequently, the correlation among the expression level of the key genes, immune-related pathways, immune cell infiltration and special immune checkpoints in GBM were assessed. Finally, the relationship between VEGFA and gene functional states were analyzed by CancerSEA (17).

Results

Patient Characteristics

Patient demographics for the imputed, training set (N = 131), and validation set (N = 56) are shown in Table 1. Patients in the training set were significantly younger than those in the validation set (p = 0.016). Most patients in the training set harbored IDH mutations, whereas most patients in the validation set had gross total resection. Compared to the training set, a higher number of patients in the validation set received current radiation and TMZ. Finally, compared to the validation set, more patients died in the training set (P < 0.001).

Calculation of the Scoring System

The point with the minimum P value of the log-rank test was considered to be the optimal cut-off value. The cut-off values were 4.16, 2.52, 193.75, 9.43, 1.61 and 2.61 for NLR, LMR, PLR, AFR, AGR and AGgR, respectively. The AGR and AGgR values were removed given the lack of significant results in the KM survival curve. In addition, the differences in overall survival caused by the presence or absence of ring-shaped gadolinium enhancement (RE), peritumoural and subventricular zone (SVZ) involvement on T2WI and insular involvement on T2WI on MRI were not statistically significant. The cut-off values were 32.56 mm and 48.44 cm³ for CPVD and the tumor volume, respectively. However, no statistically significant results for volume were obtained. Therefore, the scoring system contained the following predictor variables: NLR, PLR, LMR, AFR, and CPVD. The risk score was defined as follows: patients with high NLR (>4.16), high PLR (>193.75), low LMR (≤2.52), low AFR (≤9.43), and short CPVD (≤32.56) received a score of 5 (five abnormalities). Each abnormality was assigned one point. If all 5 parameters did not meet the standards, the patient was given a score of 0 (no abnormality). Approximately 18.3, 29.0, 28.2, 13.7, 8.4, and 2.3% of the patients had risk scores of 0, 1, 2, 3, 4, and 5, respectively. The median overall patient survival was 23.6 months, 15.4 months, 12.6 months, 6.1 months, 6.6 months and 6.0 months for risk scores of 0, 1, 2, 3, 4 and 5, respectively (P < 0.001) (Table 2).

Univariate and Multivariate Analyses for OS

These analyses were conducted using clinical data of the training set. Univariate analysis revealed that age (hazard ratio (HR) = 1.034; 95% CI = 1.020–1.049; P < 0.001), extent of resection (EOR) (HR = 0.374; 95% CI = 1.060–2.704; P = 0.027), IDH mutation (HR = 2.848; 95% CI = 1.438–5.643; P = 0.003), radiotherapy (HR = 3.139; 95% CI = 2.138–4.609; P < 0.001), chemotherapy (HR = 4.531; 95% CI = 2.892–7.100; P < 0.001) and risk score (P < 0.001) were markedly correlated with OS. Furthermore, multivariate analysis showed that age (HR = 1.030; 95% CI = 1.013–1.047; P < 0.001), EOR (HR = 0.543; 95% CI = 0.321–0.919; P = 0.023), IDH mutation (HR = 2.275; 95% CI = 1.090–4.749; P = 0.029), radiotherapy (HR = 2.072; 95% CI = 1.269–3.385; P = 0.014), chemotherapy (HR = 2.044; 95% CI = 1.142–3.660; P = 0.006) and risk score (P = 0.002) were independent predictive factors for OS (Table 3).

Nomogram and Validation

Patients with risk scores of 1 and 2 had slightly different OS. Additionally, patients with risk scores 3 and 4 had slightly different OS. A risk score 0 was defined as low risk; risk scores of 1 and 2 were defined as moderate risk; risk scores of 3 and 4 were defined as high risk; and a risk score of 5 was defined as ultrahigh risk. The baseline characteristics of the patients in the four groups are shown in Supplementary Material. Although sex was not significantly related to survival, this variable was retained in the multivariable models due to its clinical importance (Figure 1A). A nomogram was built using the training set to estimate the 12-, 24-, and 36-month survival probabilities (Figure 1B). This model showed good predictive ability with a global p value (log–rank) of 7.8321e−16, AIC of 872.71, and concordance index of 0.77. Calibration curves for the training set were established to predict 12-, 24-, and 36-month survival. Similarly, calibration curves for the validation set were established to predict 12-, 18- and 24-month survival since no patient in the validation set had an overall survival of greater than 36 months. The gray and red lines indicate the ideal survival rates and the observed survival rates, respectively. All the calibration curves were closely aligned with the 45-degree line, demonstrating good calibration. The scoring standard of the nomogram was obtained, and each patient's points were calculated. ROC curves were drawn, and the corresponding AUCs were calculated to compare the precision at 12 months (0.876), 24 months (0.834) and 36 months (0.803) in the training set. Moreover, the 12-, 18- and 24-month AUCs (0.906, 0.800, and 0.776, respectively) were calculated in the validation set. These findings validated the nomogram's good predictive ability (Figures 1C,D).

Figure 1

Peripheral Blood Test-Related Genes

NLR and LMR were associated with GBM patient prognosis. Nevertheless, the clinical significance of neutrophil, monocyte, and lymphocyte signatures from RNA sequencing data in GBM was not comprehensively evaluated. The final intersection included 38 genes (28 upregulated and 10 downregulated) (Figures 2A,B). Through univariate Cox regression analyses, six genes (VEGFA, SLC11A1, TNFRSF12A, PLAU, PTX3, and PLAUR) were significantly correlated with the OS of GBM patients in CGGA (Figure 2C). LASSO (Figure 2D) and multivariate Cox regression analyses further verified that two genes (SLC11A1 and VEGFA) represented the optimal combination for evaluating GBM patient prognosis (Figure 2E). SLC11A1 and VEGFA were more highly expressed in GBM samples compared with normal samples. The independent factors, including age, sex, radiotherapy status, 1p19q codeletion status, chemotherapy status, IDH mutation status, and MGMT promoter methylation status, were combined with the two genes. Only VEGFA (p = 0.022) was established to be an independent prognostic factor for GBM (Figure 2F). Patients with high VEGFA (HR = 1.103; CI = 1.014–1.199) expression had a significantly poorer prognosis. Cytoscape software was used to build a PPI network with 38 nodes and 94 edges based on the STRING database to identify the probable interaction network among the 38 genes (Figure 2G). The bar plot indicates the top 10 genes based on nodal numbers. These data reconfirmed that VEGFA was the core gene associated with high NLR and MLR (Figure 2H). This study also assessed the correlation between VEGFA and clinical factors. The results showed that VEGFA expression levels were markedly elevated in the IDH wild type compared with the IDH mutant type (Figure 3). Finally, the NLR and LMR of samples in CGGA, TCGA and GSE16011 were estimated by using CIBERSORT. The results showed NLR was statistically positively associated with expression of VEGFA in three datasets. However, the association between LMR and VEGFA were not statistically significant (Supplementary Figure 1).

Figure 2

Figure 3

Multivariate Analysis of VEGFA

GSEA was applied to identify enriched features as well as functional differences between the low and high VEGFA expression groups using the top five entries of GO and KEGG terms. The high VEGFA expression group was found to be enriched in “ACTIN_CYTOSKELETON, ADAPTIVE_ IMMUNE_ RESPONSE, AMEBOIDAL_ TYPE_ CELL_ MIGRATION, APICAL_ PART_OF_ CELL and CANONICAL_ WNT_ SIGNALING_ PATHWAY” in GO and “CYTOKINE_ RECEPTOR_ INTERACTION, FOCAL_ ADHESION, NEUROACTIVE_ LIGAND_ RECEPTOR_ INTERACTION, PATHWAYS_ IN_CANCER and REGULATION_OF_ACTIN_ CYTOSKELETON” in KEGG. The enrichment scores of different immune cells and immune-associated functions as well as pathways were further evaluated based on the ssGSEA algorithm to assess the correlations between VEGFA expression and immune infiltration. Seven immune cells (Tregs, Th2 cells, Th1 cells, T helper cells, pDCs, iDCs, and macrophages) were markedly correlated with VEGFA expression. Moreover, nine immune-related functions (type II IFN response, parainflammation, CCR, MHC class I, inflammation promotion, checkpoint, T cell costimulation, APC costimulation, and cytolytic activity) were significantly correlated with VEGFA expression. In addition, six immune checkpoints (PD-1, PD-L1, TIM-3, B7-H3, CD40 and CD28) were closely associated with VEGFA, suggesting a potential association between VEGFA expression and immune checkpoint inhibitor treatment. ssGSEA confirmed that VEGFA levels might affect the immune status in the tumor microenvironment. Finally, single-cell RNA sequencing results showed that hypoxia increased VEGFA expression. The expression of VEGFA was significantly positively associated with hypoxia (R = 0.36) and angiogenesis (R = 0.27) (Figure 4).

Figure 4

Discussion

Despite considerable advances in brain cancer therapy in the past decade, the prognosis of GBM is extremely poor (2). Comprehensive preoperative assessment tools are needed to classify GBM risk, thereby appropriately guiding further individualized treatment. Herein, peripheral blood test indicators and imaging indicators of 187 GBM patients were used to define a scoring system based on NLR, PLR, LMR, AFR, and CPVD. An individual survival nomogram for GBM patients was also developed and validated. The CPH survival model included basic patient information (age and sex), preoperative scoring system, extent of tumor resection, IDH status, and subsequent treatment after surgery. The AUC values of the ROC curve were 0.876 (12-month ROC), 0.834 (24-month ROC) and 0.803 (36-month ROC) in the training set and 0.906 (12-month ROC), 0.800 (18-month ROC) and 0.776 (24-month ROC) in the validation set. Compared to other models, our prediction models provide better prediction results and enhance the evaluation of the entire course of the disease. However, our risk score combined with multielement data included not only peripheral blood tests but also imaging factors and was finally combined with clinical risk factors to construct a more comprehensive and more individualized prognostic model in GBM. Many studies have constructed a series of risk score systems to evaluate the prognosis of gliomas. For example, Wu et al. constructed a risk score based on preoperative fibrinogen (F), the neutrophil to lymphocyte ratio (NLR), and albumin to globulin ratio (AGR) and found that a high F-NLR-AGR score was an independent predictor of poor prognosis; however, this study did not incorporate imaging analysis (27). Similarly, Zhang constructed a radiomics nomogram from 4,000 radiomics features, demonstrating a good predictive ability in both the training and validation sets without incorporating the peripheral blood analysis results (28). Finally, Wang et al. constructed a model based on radiomic signatures and blood urea nitrogen (BUN). However, our study used a larger sample size, and our model showed better predictive ability (29).

Recent studies have reported that GBM tumor cells may increase the neutrophil count both in the peripheral blood and around the tumor by secreting chemotactic factors, including G-CSF, VEGF, IL-1β, and IL-6, and neutrophils may promote tumor progression by secreting important cytokines, such as VEGF, IL-6, IL-8, elastases, and matrix metalloproteinases (30, 31). In addition, an increasing neutrophil count may regulate lymphocyte function by releasing reactive oxygen species and arginase, thereby inhibiting lymphocyte survival and normal cytotoxic function (32). It has also been reported that elevated platelet levels may promote tumor growth, angiogenesis, and dissemination by secreting crucial factors, such as VEGF (33). In addition, platelets release increased levels of soluble CD40, which is a known inhibitor of regulatory T cell recruitment that may promote the immunosuppressive microenvironment and eventually form an environment conducive to tumor growth (34). Abundant macrophage infiltration is a common feature of GBMs, but these tumor-associated macrophages (TAMs) lack apparent phagocytic activity in GBMs. Recent studies have shown that macrophages can be recruited and induced to become the M2 type by a wide variety of factors secreted by glioma cells, including IL-10, IL-4, IL-6, macrophage-colony stimulating factor (M-CSF), TGF-β, and prostaglandin E2 (PGE2), and promote cancer cell proliferation in GBMs (35, 36). In addition, M2-like macrophages also significantly promote angiogenesis and induce immunosuppression (37). Therefore, we found that neutrophil platelet count and macrophage count may be negatively correlated with lymphocyte count, and high NLR, PLR and MLR may promote angiogenesis in GBM. In our study, patients with a higher NLR (>4.16), higher PLR (>193.75) and lower LMR (≤2.52) had a poorer prognosis, and the outcomes were consistent with past studies.

Regarding the AFR, mounting evidence indicates that fibrinogen is an important regulator of tumor progression and systemic inflammatory responses in several malignant cancers and that hyperfibrinogenaemia is associated with a high invasiveness of GBM (38). Fibrinogen promotes tumor angiogenesis and increases the adhesive, migratory, and invasive abilities of tumor cells (39–41). Moreover, the physical barrier formed by platelet-fibrin deposition surrounding tumor cells can prevent NK cell destruction (42). Malignant tumor cells can synthesize fibrinogen, thus promoting tumor cell growth and angiogenesis through interactions with vascular endothelial growth factor and fibroblast growth factor-2 (43, 44). As a malnutrition and inflammation marker, a reduced albumin level is a risk factor for malignant tumors. Therefore, a lower AFR may imply inadequate antitumour immunity, malnutrition, and cancer-associated inflammation, which are unfavorable for cancer prognosis.

Matarredona et al. (20) indicated that the tumor location and the subventricular zone (SVZ) are related to a patient's prognosis. The SVZ of the adult human brain has several neural stem cells (NSCs), which can undergo multilineage differentiation. In the SVZ, NSCs are potential cells of origin containing driver mutations of human GBM (45). This study demonstrated that GBM patients with elevated CPVD have shorter OS than those with low CPVD. Thus, CPVD may be a marker for tumor invasion in GBM.

VEGFA is essential for physiological and pathological angiogenesis, and bevacizumab, a molecular-targeted drug, binds and neutralizes human VEGFA to suppress the VEGF signaling pathway. Given its antiangiogenic effects, bevacizumab has become the standard GBM treatment and positively affects the quality of life and survival of recurrent GBM patients. However, bevacizumab increases PFS but has no effect on OS. N. García-Romero et al. also reported that tumor growth of greater than 20% in GBM patients is independent of the VEGFA pathway and thus does not benefit from antiangiogenic therapy (46–48). Herein, differentially expressed genes between neutrophils and lymphocytes as well as monocytes and lymphocytes were screened by differential gene expression analysis. VEGFA was not only significantly related to survival but also ranked first in the protein–protein interaction (PPI) network map. Therefore, we confirmed that VEGFA expression was closely related to our risk score level, which was manifested through NLR, LMR and PLR. Based on the current achievements, we propose to perform prospective clinical trials in which we obtain VEGFA expression levels using ddPCR to detect and quantify ctDNA levels in GBM patient blood and cerebrospinal fluid to confirm the close correlation between VEGFA expression and peripheral blood tests. We also attempted to explore the prognostic predictive ability of the ctDNA level of VEGFA and whether ctDNA can be a potential molecular indicator of sensitivity or resistance to bevacizumab.

Conclusions

In summary, the nomogram is a novel tool for predicting the prognosis of GBM patients, thereby informing individualized treatment. In addition, VEGFA may influence the immune microenvironment of GBM, leading to poor prognosis in GBM patients.

Funding

This research was supported by grants from followings: Wenzhou Science and Technology Project under Grants, Y20190144, Key Research Project of Traditional Chinese Medicine of Zhejiang Province of China, 2019ZZ015, and Medical Health Science and Technology Research Project of Zhejiang Province of China, 2018KY515.

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.

  OstromQTCioffiGGittlemanHPatilNWaiteKKruchkoCet al. Statistical report: primary brain and other central nervous system tumors diagnosed in the United States in 2012-2016. Neuro Oncol. (2019) 21:v1–v100. 10.1093/neuonc/noz150

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2.

  ZinnPOColenRRKasperEMBurkhardtJK. Extent of resection and radiotherapy in GBM: A 1973 to 2007 surveillance, epidemiology and end results analysis of 21,783 patients. Int J Oncol. (2013) 42:929–34. 10.3892/ijo.2013.1770

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3.

  DarefskyASKingJTDubrowR. Adult glioblastoma multiforme survival in the temozolomide era: a population-based analysis of Surveillance, Epidemiology, and End Results registries. Cancer. (2012) 118:2163–72. 10.1002/cncr.26494

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4.

  ThakkarJPDolecekTAHorbinskiCOstromQTLightnerDDBarnholtz-SloanJSet al. Epidemiologic and molecular prognostic review of glioblastoma. Cancer Epidemiol Biomarkers Prev. (2014) 23:1985–96. 10.1158/1055-9965.EPI-14-0275

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5.

  GaistDHallasJFriisSHansenSSørensenHT. Statin use and survival following glioblastoma multiforme. Cancer Epidemiol. (2014) 38:722–7. 10.1016/j.canep.2014.09.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6.

  OstromQTGittlemanHTruittGBosciaAKruchkoCBarnholtz-SloanJSet al. Statistical report: primary brain and other central nervous system tumors diagnosed in the United States in 2011-2015. Neuro Oncol. (2018) 20:iv1–iv86. 10.1093/neuonc/noy131

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7.

  LouisDNPerryAReifenbergerGvon DeimlingAFigarella-BrangerDCaveneeWKet al. The 2016 world health organization classification of tumors of the central nervous system: a summary. Acta Neuropathol. (2016) 131:803–20. 10.1007/s00401-016-1545-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8.

  StuppRBradaMvan den BentMJTonnJCPentheroudakisG. High-grade glioma: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol. (2014) 3:iii93–101. 10.1093/annonc/mdu050

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9.

  CeccarelliMBarthelFPMaltaTMSabedotTSSalamaSRMurrayBAet al. Molecular profiling reveals biologically discrete subsets and pathways of progression in diffuse glioma. Cell. (2016) 164:550–63. 10.1016/j.cell.2015.12.028

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10.

  MarshallHTDjamgozMBA. Immuno-oncology: emerging targets and combination therapies. Front Oncol. (2018) 8:315. 10.3389/fonc.2018.00315

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11.

  KayhanAKorkmazTSBaranOKemerdereRYeniSNTanriverdiT. Preoperative systemic inflammatory markers in different brain pathologies: an analysis of 140 patients. Turk Neurosurg. (2019) 29:799–803. 10.5137/1019-5149.JTN.24244-18.2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12.

  AuezovaRRyskeldievNDoskaliyevAKuanyshevYZhetpisbaevBAldiyarovaNet al. Association of preoperative levels of selected blood inflammatory markers with prognosis in gliomas. Onco Targets Ther. (2016) 9:6111–7. 10.2147/OTT.S113606

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13.

  XuWWangDZhengXOuQHuangL. Sex-dependent association of preoperative hematologic markers with glioma grade and progression. J Neurooncol. (2018) 137:279–87. 10.1007/s11060-017-2714-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14.

  PerisanidisCPsyrriACohenEEEngelmannJHeinzeGPerisanidisBet al. Prognostic role of pretreatment plasma fibrinogen in patients with solid tumors: A systematic review and meta-analysis. Cancer Treat Rev. (2015) 41:960–70. 10.1016/j.ctrv.2015.10.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15.

  PatelSHPoissonLMBratDJZhouYCooperLSnuderlMet al. T2-FLAIR mismatch, an imaging biomarker for IDH and 1p/19q status in lower-grade gliomas: a TCGA/TCIA project. Clin Cancer Res. (2017) 23:6078–85. 10.1158/1078-0432.CCR-17-0560

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16.

  VerhaakRGHoadleyKAPurdomEWangVQiYWilkersonMDet al. Cancer Genome Atlas Research Network. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell. (2010) 17:98–110. 10.1016/j.ccr.2009.12.020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17.

  YuanHYanMZhangGLiuWDengCLiaoGXuLLuoTYanHLongZShiAZhaoTXiaoYLiX. CancerSEA: a cancer single-cell state atlas. Nucleic Acids Res. (2019) 47:D900–D908. 10.1093/nar/gky939

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18.

  PieralliniABonaminiMBozzaoAPantanoPStefanoDDFeroneEet al. Supratentorial diffuse astrocytic tumours: proposal of an MRI classification. Eur Radiol. (1997) 7:395–9. 10.1007/s003300050173

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19.

  ZhaoKSunGWangQXueZLiuGXiaYet al. The diagnostic value of conventional MRI and CT features in the identification of the IDH1-Mutant and 1p/19q co-deletion in WHO grade II gliomas. Acad Radiol. (2021) 28:e189–98. 10.1016/j.acra.2020.03.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20.

  MatarredonaERPastorAM. Neural stem cells of the subventricular zone as the origin of human glioblastoma stem cells. Therapeutic implications. Front Oncol. (2019) 9:779. 10.3389/fonc.2019.00779

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21.

  PalludJBlonskiMMandonnetEAudureauEFontaineDSanaiNet al. Velocity of tumor spontaneous expansion predicts long-term outcomes for diffuse low-grade gliomas. Neuro Oncol. (2013) 15:595–606. 10.1093/neuonc/nos331

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22.

  FriedmanJHastieTTibshiraniR. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. (2010) 33:1–22. 10.18637/jss.v033.i01

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23.

  IshwaranHKogalurUB. Fast unified random forests for survival, regression, and classification (RF-SRC), R package version 2.9.0. (2019).

  • Google Scholar

• 24.

  ThulPJÅkessonLWikingMMahdessianDGeladakiAAit BlalHet al. A subcellular map of the human proteome. Science. (2017) 356:eaal3321. 10.1126/science.aal3321

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25.

  NewmanAMLiuCLGreenMRGentlesAJFengWXuYet al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. 10.1038/nmeth.3337

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26.

  RooneyMSShuklaSAWuCJGetzGHacohenN. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. (2015) 160:48–61. 10.1016/j.cell.2014.12.033

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27.

  WuYSongZSunKRongSGaoPWangFet al. novel scoring system based on peripheral blood test in predicting grade and prognosis of patients with glioma. Onco Targets Ther. (2019) 12:11413–23. 10.2147/OTT.S236598

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28.

  ZhangXLuHTianQFengNYinLXuXet al. radiomics nomogram based on multiparametric MRI might stratify glioblastoma patients according to survival. Eur Radiol. (2019) 29:5528–38. 10.1007/s00330-019-06069-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29.

  WangJYiXFuYPangPDengHTangHet al. Preoperative Magnetic Resonance Imaging Radiomics for Predicting Early Recurrence of Glioblastoma. Front Oncol. (2021) 11:769188. 10.3389/fonc.2021.769188

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30.

  ScapiniPNesiLMoriniMTanghettiEBelleriMNoonanDet al. Generation of biologically active angiostatin kringle 1-3 by activated human neutrophils. J Immunol. (2002) 168:5798–804. 10.4049/jimmunol.168.11.5798

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31.

  LechnerMGLiebertzDJEpsteinAL. Characterization of cytokine-induced myeloid-derived suppressor cells from normal human peripheral blood mononuclear cells. J Immunol. (2010) 185:2273–84. 10.4049/jimmunol.1000901

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32.

  GregoryADHoughtonAM. Tumor-associated neutrophils: new targets for cancer therapy. Cancer Res. (2011) 71:2411–6. 10.1158/0008-5472.CAN-10-2583

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33.

  MenterDGTuckerSCKopetzSSoodAKCrissmanJDHonnKV. Platelets and cancer: a casual or causal relationship: revisited. Cancer Metastasis Rev. (2014) 33:231–69. 10.1007/s10555-014-9498-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34.

  McCourtMWangJHSookhaiSRedmondHP. Proinflammatory mediators stimulate neutrophil-directed angiogenesis. Arch Surg. (1999) 134:1325–31. 10.1001/archsurg.134.12.1325

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35.

  WachJApallasSSchneiderMWellerJSchussPVatterHet al. Mean platelet volume/platelet count ratio and risk of progression in glioblastoma. Front Oncol. (2021) 11:695316. 10.3389/fonc.2021.695316

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36.

  CuiXMoralesRTQianWWangHGagnerJPDolgalevIet al. Hacking macrophage-associated immunosuppression for regulating glioblastoma angiogenesis. Biomaterials. (2018) 161:164–78. 10.1016/j.biomaterials.2018.01.053

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37.

  ZhouWKeSQHuangZFlavahanWFangXPaulJet al. Periostin secreted by glioblastoma stem cells recruits M2 tumour-associated macrophages and promotes malignant growth. Nat Cell Biol. (2015) 17:170–82. 10.1038/ncb3090

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38.

  PalumboJSTalmageKEMassariJVLa JeunesseCMFlickMJKombrinckKWet al. Platelets and fibrin(ogen) increase metastatic potential by impeding natural killer cell-mediated elimination of tumor cells. Blood. (2005) 105:178–85. 10.1182/blood-2004-06-2272

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39.

  Simpson-HaidarisPJRybarczykB. Tumors and fibrinogen. The role of fibrinogen as an extracellular matrix protein. Ann N Y Acad Sci. (2001) 936:406–25. 10.1111/j.1749-6632.2001.tb03525.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40.

  LiuXLiuZLinEChenYSunXZhouZet al. cumulative score based on preoperative fibrinogen and the neutrophil-lymphocyte ratio to predict outcomes in resectable gastric cancer. Cancer Manag Res. (2018) 10:3007–14. 10.2147/CMAR.S174656

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41.

  NieswandtBHafnerMEchtenacherBMännelDN. Lysis of tumor cells by natural killer cells in mice is impeded by platelets. Cancer Res. (1999) 59:1295–300.

  • Pubmed Abstract
  • Google Scholar

• 42.

  SahniAFrancisCW. Vascular endothelial growth factor binds to fibrinogen and fibrin and stimulates endothelial cell proliferation. Blood. (2000) 96:3772–8. 10.1182/blood.V96.12.3772.h8003772_3772_3778

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43.

  SahniASimpson-HaidarisPJSahniSKVadayGGFrancisCW. Fibrinogen synthesized by cancer cells augments the proliferative effect of fibroblast growth factor-2 (FGF-2). J Thromb Haemost. (2008) 6:176–83. 10.1111/j.1538-7836.2007.02808.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44.

  ChiangGCPisapiaDJLiechtyBMaggeRRamakrishnaRKniselyJet al. The prognostic value of MRI subventricular zone involvement and tumor genetics in lower grade gliomas. J Neuroimaging. (2020) 30:901–9. 10.1111/jon.12763

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45.

  FengHLiuKShenXLiangJWangCQiuWet al. Targeting tumor cell-derived CCL2 as a strategy to overcome Bevacizumab resistance in ETV5+ colorectal cancer. Cell Death Dis. (2020) 11:916. 10.1038/s41419-020-03208-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46.

  SimonTPiniotiSSchellenbergerPRajeeveVWendlerFCutillasPRet al. Shedding of bevacizumab in tumour cells-derived extracellular vesicles as a new therapeutic escape mechanism in glioblastoma. Mol Cancer. (2018) 17:132. 10.1186/s12943-018-0878-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47.

  GilbertMRDignamJJArmstrongTSWefelJSBlumenthalDTVogelbaumMAet al. A randomized trial of bevacizumab for newly diagnosed glioblastoma. N Engl J Med. (2014) 370:699–708. 10.1056/NEJMoa1308573

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48.

  García-RomeroNPalacín-AlianaIMadurgaRCarrión-NavarroJEsteban-RubioSJiménezBet al. Bevacizumab dose adjustment to improve clinical outcomes of glioblastoma. BMC Med. (2020) 18:142. 10.1186/s12916-020-01610-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar