Bioinformatic and Statistical Analysis

The human glioblastoma analysis was first performed on the TCGA dataset, available as ‘Tumour Glioblastoma - TCGA - MAS 5.0 - u133a’ (Network, C. G. A. R 2008) in the R project 3.5.0, and source data are available at https://www.cbioportal.org. The TCGA database used comprises information about 540 patients, including 85 samples sub-classified in classical (n = 17), mesenchymal (n = 27), neural (n = 17) and proneural (n = 24) GB. While,for CGGA datasets, data were downloaded from http://www.cgga.org.cn. The mRNA expression and suvival analysis was performed only on primary or “de novo” GB cohort (n=220). The survival statistical analysis was performed using R2 webtool (36) and Survminer package (37). The Survminer package provides Kaplan-Meyer plots based on a Log Rank Scale p-value for the comparison of the subgroups. The optimal cut off was calculated by maximally selected rank statistics using Maxstat R-package (38).