Focusing on scRNA-seq-Derived T Cell-Associated Genes to Identify Prognostic Signature and Immune Microenvironment...

The clinical outcomes of low-grade glioma (LGG) are associated with T cell infiltration, but the specific contribution of heterogeneous T cell types remains unclear.

To study the different functions of T cells in LGG, we mapped the single-cell RNA sequencing results of 10 LGG samples to obtain T cell marker genes. In addition, bulk RNA data of 975 LGG samples were collected for model construction. Algorithms such as TIMER, CIBERSORT, QUANTISEQ, MCPCOUTER, XCELL, and EPIC were used to depict the tumor microenvironment landscape. Subsequently, three immunotherapy cohorts, PRJEB23709, GSE78820, and IMvigor210, were used to explore the efficacy of immunotherapy.

The Human Primary Cell Atlas was used as a reference dataset to identify each cell cluster; a total of 15 cell clusters were defined and cells in cluster 12 were defined as T cells. According to the distribution of T cell subsets (CD4+ T cell, CD8+ T cell, Naïve T cell, and Treg cell), we selected the differentially expressed genes. Among the CD4+ T cell subsets, we screened 3 T cell-related genes, and the rest were 28, 4, and 13, respectively. Subsequently, according to the T cell marker genes, we screened six genes for constructing the model, namely, RTN1, HERPUD1, MX1, SEC61G, HOPX, and CHI3L1. The ROC curve showed that the predictive ability of the prognostic model for 1, 3, and 5 years was 0.881, 0.817, and 0.749 in the TCGA cohort, respectively. In addition, we found that risk scores were positively correlated with immune infiltration and immune checkpoints. To this end, we obtained three immunotherapy cohorts to verify their predictive ability of immunotherapy effects and found that high-risk patients had better clinical effects of immunotherapy.

This single-cell RNA sequencing combined with bulk RNA sequencing may elucidate the composition of the tumor microenvironment and pave the way for the treatment of low-grade gliomas.