Identification and Validation of an 11-Ferroptosis Related Gene Signature and Its Correlation With Immune Checkpoint Molecules in Glioma

Front Cell Dev Biol. 2021 Jun 23:9:652599. doi: 10.3389/fcell.2021.652599. eCollection 2021.

Abstract

Background: Glioma is the most common primary malignant brain tumor with significant mortality and morbidity. Ferroptosis, a novel form of programmed cell death (PCD), is critically involved in tumorigenesis, progression and metastatic processes.

Methods: We revealed the relationship between ferroptosis-related genes and glioma by analyzing the mRNA expression profiles from The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), GSE16011, and the Repository of Molecular Brain Neoplasia Data (REMBRANDT) datasets. The least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed to construct a ferroptosis-associated gene signature in the TCGA cohort. Glioma patients from the CGGA, GSE16011, and REMBRANDT cohorts were used to validate the efficacy of the signature. Receiver operating characteristic (ROC) curve analysis was applied to measure the predictive performance of the risk score for overall survival (OS). Univariate and multivariate Cox regression analyses of the 11-gene signature were performed to determine whether the ability of the prognostic signature in predicting OS was independent. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to identify the potential biological functions and pathways of the signature. Subsequently, we performed single sample gene set enrichment analysis (ssGSEA) to explore the correlation between risk scores and immune status. Finally, seven putative small molecule drugs were predicted by Connectivity Map.

Results: The 11-gene signature was identified to divide patients into two risk groups. ROC curve analysis indicated the 11-gene signature as a potential diagnostic factor in glioma patients. Multivariate Cox regression analyses showed that the risk score was an independent predictive factor for overall survival. Functional analysis revealed that genes were enriched in iron-related molecular functions and immune-related biological processes. The results of ssGSEA indicated that the 11-gene signature was correlated with the initiation and progression of glioma. The small molecule drugs we selected showed significant potential to be used as putative drugs.

Conclusion: we identified a novel ferroptosis-related gene signature for prognostic prediction in glioma patients and revealed the relationship between ferroptosis-related genes and immune checkpoint molecules.

Keywords: ferroptosis; gene signature; glioma; prognosis; tumor immunity.