Despite advances in molecular characterization of glioblastoma multiforme (GBM), only a handful of predictive biomarkers exist with limited clinical relevance. We aimed to identify differentially expressed genes in tumor samples collected at surgery associated with response to subsequent treatment, including temozolomide (TMZ) and nitrosoureas.

Gene expression was collected from multiple independent datasets. Patients were categorized as responders/nonresponders based on their survival status at 16 months post-surgery. For each gene, the expression was compared between responders and nonresponders with a Mann-Whitney U test and receiver operating characteristic. The package „roc” was used to calculate the area under the curve (AUC).

The integrated database comprises 454 GBM patients from three independent datasets and 10,103 genes. The highest proportion of responders (68%) were among patients treated with TMZ combined with nitrosoureas, where FCGR2B upregulation provided the strongest predictive value (AUC=0.72, p < 0.001). Elevated expression of CSTA and MRPS17 was associated with a lack of response to multiple treatment strategies. DLL3 upregulation was present in subsequent responders to any treatment combination containing TMZ. Three genes (PLSCR1, MX1, and MDM2) upregulated both in the younger cohort and in patients expressing low MGMT delineate a subset of patients with worse prognosis within a population generally associated with a favorable outcome.

The identified transcriptomic changes provide biomarkers of responsiveness, offer avenues for preclinical studies, and may enhance future GBM patient stratifications. The described methodology provides a reliable pipeline for the initial testing of potential biomarker candidates for future validation studies.