Construction and validation of a TAMRGs prognostic signature for gliomas by integrated analysis of scRNA and bulk RNA sequencing data

Background

This study aimed to construct and validate a prognostic model based on tumor associated macrophage-related genes (TAMRGs) by integrating single-cell RNA sequencing (scRNA-seq) and bulk RNA sequencing (bulk RNA-seq) data.

Methods

The scRNA-seq data of three inhouse glioma tissues were used to identify the tumor-associated macrophages (TAMs) marker genes, the DEGs from the The Cancer Genome Atlas (TCGA) − Genotype-Tissue Expression (GTEx) dataset were used to further select TAMs marker genes. Subsequently, a TAMRG-score was constructed by Least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox regression analysis in the TCGA dataset and validated in the Chinese Glioma Genome Atlas (CGGA) dataset.

Results

We identified 186 TAMs marker genes, and a total of 6 optimal prognostic genes including CKS2, LITAF, CTSB, TWISTNB, PPIF and G0S2 were selected to construct a TAMRG-score. The high TAMRG-score was significantly associated with worse prognosis (log-rank test, P<0.001). Moreover, the TAMRG-score outperformed the other three models with AUC of 0.808. Immune cell infiltration, TME scores, immune checkpoints, TMB and drug susceptibility were significantly different between TAMRG-score groups. In addition, a nomogram were constructed by combing the TAMRG-score and clinical information (Age, Grade, IDH mutation and 1p19q codeletion) to predict the survival of glioma patients with AUC of 0.909 for 1-year survival.

Conclusion

The high TAMRG-score group was associated with a poor prognosis. A nomogram by incorporating TMARG-score could precisely predict glioma survival, and provide evidence for personalized treatment of glioma.