GBMdeconvoluteR accurately infers proportions of neoplastic and immune cell populations from bulk glioblastoma transcriptomics data

AbstractThe biological and clinical impact of neoplastic and immune cell type ratios in the glioblastoma (GBM) tumour microenvironment is being realised. Characterising and quantifying cell types within GBMs at scale will facilitate a better understanding of the association between the cellular landscape and tumour phenotypes or clinical correlates. This study aimed to develop a tool that can deconvolute immune and neoplastic cells within the GBM tumour microenvironment from bulk RNA sequencing data. We developed an IDH wild-type (IDHwt) GBM specific single immune cell reference dataset, from four independent studies, consisting of B cells, T cells, NK cells, microglia, tumour associated macrophages, monocytes, mast and DC cells. We used this alongside an existing neoplastic single cell-type dataset consisting of astrocyte-like, oligodendrocyte- and neuronal-progenitor like and mesenchymal GBM cancer cells to create both marker and gene signature matrix-based deconvolution tools. We then applied single-cell resolution imaging mass cytometry (IMC) to ten IDHwt GBM samples, five paired primary and recurrent tumours, in parallel with these tools to determine which performed best. Marker based gene expression deconvolution using GBM tissue specific markers, which we have packaged as GBMdeconvoluteR, gave the most accurate results. The correlation between immune cell quantification by IMC and by GBMdeconvoluteR for primary IDHwt GBM samples was 0.52 (Pearson’s P=7.8×10−3) and between neoplastic cell quantification by IMC and by GBMdeconvoluteR was 0.75 (Pearson’s P=1.2×10−3). We applied GBMdeconvoluteR to bulk GBM RNAseq data from The Cancer Genome Atlas (TCGA) and were able to recapitulate recent findings from multi-omics single cell studies with regards associations between mesenchymal GBM cancer cells and both lymphoid and myeloid cells. Furthermore, we were able to expand upon this to show that these associations are stronger in patients with worse prognosis. GBMdeconvoluteR is accessible online athttps://gbmdeconvoluter.leeds.ac.uk.Key pointsGBMdeconvoluteR is a glioblastoma-specific cellular deconvolution tool. When applied to bulk GBM RNAseq data, it accurately quantifies the neoplastic and immune cells in that tumour. It is available online athttps://gbmdeconvoluter.leeds.ac.uk