Abstract
The tumor micro-environment (TME) plays an important role in the biology of cancer, including gliomas. Single cell studies have highlighted the role of specific TME components in gliomas, and the methods to deconvolve bulk profiling data may serve to complement these studies on clinically annotated tumors. In this study, we estimated cell type proportions in 3 large glioma datasets (TCGA, CGGA-325, CGGA-693) using CIBERSORTx. Using a signature matrix comprising 22 immune cell types, we identified IDH mutation status-specific immune cell distributions and found that the proportions of 10 cell types were significantly different between IDHmut and IDHwt tumors across the 3 datasets. Looking further within IDHmut tumors, we found that monocytes were enriched in 1p/19q non-co-deleted tumors across the 3 glioma datasets, consistent with prior single cell studies. We then examined estimated gene expression among immune cell types relative to IDH mutation status and found clear separation of gene expression in 15 of 22 cell types in all 3 datasets. When we applied these 22 gene expression signatures in each tumor sample onto cluster-of-cluster analyses to identify tumor groups with distinct immune signature patterns, we found that samples were distributed largely according to the IDH status in all 3 datasets, confirming that immune cell expression is distinct based on IDH status. Among IDH-specific groups, cluster-of-cluster analyses showed that immune cell-based cluster groups had distinct survival outcomes, and that IDHwt samples were distributed significantly based on tumor grades as well as based on EGFR overexpression. Among IDHmut tumors, the distributions of tumor grade and 1p/19q co-deletion status were significantly different in the immune-based clusters in 2 of the 3 datasets examined. Overall, these results highlight the biological and clinical significance of the immune cell environment in gliomas, including distinctions based on IDH mutation status as well as prognosis within IDH-specific groups.