Abstract
Background: Glioblastoma (GBM) is the most common primary malignant intracranial tumor and closely related to metabolic alteration. However, few accepted prognostic models are currently available, especially models based on metabolic genes. Methods: The transcriptome data were obtained for all of the patients diagnosed with GBM from the Gene Expression Omnibus (GEO) (training cohort, n=369) and The Cancer Genome Atlas (TCGA) (validation cohort, n=152) with the following variables: age at diagnosis, sex, follow-up and overall survival (OS). Metabolic genes according to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were contracted, and a Lasso regression model was constructed. Survival was assessed by univariate or multivariate Cox proportional hazards regression and Kaplan-Meier analysis, and an independent external validation was also conducted to examine the model. Results: There were 341 metabolic genes showed significant differences between normal brain and GBM tissues in both the training and validation cohorts, among which 56 genes were dramatically correlated to the OS of patients. Lasso regression revealed that the metabolic prognostic model was composed of 18 genes, including COX10, COMT, and GPX2 with protective effects, as well as OCRL and RRM2 with unfavorable effects. Patients classified as high-risk by the risk score from this model had markedly shorter OS than low-risk patients (P<0.0001), and this significant result was also observed in independent external validation (P<0.001). Conclusions: The prognosis of GBM was dramatically related to metabolic pathways, and our metabolic prognostic model had high accuracy and application value in predicting the OS of GBM patients.