Higher cytolytic score correlates with an immunosuppressive tumor microenvironment and reduced survival in glioblastoma

Cytolytic score (CYT), calculated from mRNA expression levels of granzyme and perforin, positively correlates with CD8+ T cell infiltration/activity in a variety of cancers. Unlike other cancers, higher CYT has been associated with worse prognosis in glioblastoma (GBM). To address this discrepancy, we sought to investigate the relationship between CYT and immune checkpoint gene score (ICGscore), as well as their correlation with patient survival and tumor immune cell infiltration. Clinical and RNA-sequencing data for patients with newly diagnosed GBM were obtained from The Cancer Genome Atlas. Maximally-selected rank statistics was used to dichotomize subgroups. CIBERSORT was used to estimate abudence of immune cell-types. Spearman correlation was used to characterize the relationship between CYT and ICGscore. Kaplan–Meier curves were generated for survival analysis. Overall, 28/151 patients had high CYT. High CYT was associated with a mesenchymal subtype (p < 0.001) and worse survival (7.45 vs. 12.2 months, p < 0.001). There were no differences in patient demographics, IDH/MGMT mutation status, or treatment. On subgroup analysis, patients with high CYT/ICGscore had significantly increased CD8+ infiltration (p < 0.001), as expected, and worse survival (HR 0.445, p < 0.01). Furthermore, CYT strongly correlated with ICGscore (R(S) = 0.675, p < 0.001). The high CYT/ICGscore subgroup was associated with greater infiltration of M2 macrophages (p = 0.011) and neutrophils (p = 0.055). Our study highlights a multidimensional immunosuppressive GBM microenvironment in patients with higher CYT and potentially identifies patients with high CYT/ICGscore as a subgroup that may particularly benefit from multi-faceted immunotherapies, given their already elevated tumor CD8+ T cell levels.