Identification of Immunological Characteristics and Immune Subtypes Based on Single-Sample Gene Set Enrichment Analysis Algorithm in Lower-Grade Glioma

Abstract

Few breakthroughs have been achieved in the treatment of lower-grade glioma (LGG) in recent decades. Apart from the conventional pathological and histological classifications, subtypes based on immunogenomics would provide reference for individualized treatment and prognosis prediction. Our study identified four immunotypes of lower-grade glioma (clusters A, B, C, and D) by bioinformatics methods in TCGA-LGG and two CGGA datasets. Cluster A was an “immune-cold” phenotype with the lowest immune infiltration and longest survival expectation, whereas cluster D was an “immune-rich” subtype with the highest immune infiltration and poor survival expectation. The expression of immune checkpoints increased along with immune infiltration degrees among the clusters. It was notable that immune clusters correlated with a variety of clinical and immunogenomic factors such as age, WHO grades, IDH1/2 mutation, PTEN, EGFR, ATRX, and TP53 status. In addition, LGGs in cluster D were sensitive to cisplatin, gemcitabine, and immune checkpoint PD-1 inhibitors. RTK-RAS and TP53 pathways were affected in cluster D. Functional pathways such as cytokine–cytokine receptor interaction, antigen processing and presentation, cell adhesion molecules (CAMs), and ECM–receptor interaction were also enriched in cluster D. Hub genes were selected by the Matthews correlation coefficient (MCC) algorithm in the blue module of a gene co-expression network. Our studies might provide an immunogenomics subtyping reference for immunotherapy in LGG.