Molecular subtyping of acute myeloid leukemia through ferroptosis signatures predicts prognosis and deciphers the immune microenvironment

Abstract

Acute myeloid leukemia (AML) is one of the most aggressive hematological malignancies with a low 5-year survival rate and high rate of relapse. Developing more efficient therapies is an urgent need for AML treatment. Accumulating evidence showed that ferroptosis, an iron-dependent form of programmed cell death, is closely correlated with cancer initiation and clinical outcome through reshaping the tumor microenvironment. However, understanding of AML heterogeneity based on extensive profiling of ferroptosis signatures remains to be investigated yet. Herein, five independent AML transcriptomic datasets (TCGA-AML, GSE37642, GSE12417, GSE10358, and GSE106291) were obtained from the GEO and TCGA databases. Then, we identified two ferroptosis-related molecular subtypes (C1 and C2) with distinct prognosis and tumor immune microenvironment (TIME) by consensus clustering. Patients in the C1 subtype were associated with favorable clinical outcomes and increased cytotoxic immune cell infiltration, including CD8+/central memory T cells, natural killer (NK) cells, and non-regulatory CD4+ T cells while showing decreased suppressive immune subsets such as M2 macrophages, neutrophils, and monocytes. Functional enrichment analysis of differentially expressed genes (DEGs) implied that cell activation involved in immune response, leukocyte cell–cell adhesion and migration, and cytokine production were the main biological processes. Phagosome, antigen processing and presentation, cytokine–cytokine receptor interaction, B-cell receptor, and chemokine were identified as the major pathways. To seize the distinct landscape in C1 vs. C2 subtypes, a 5-gene prognostic signature (LSP1, IL1R2, MPO, CRIP1, and SLC24A3) was developed using LASSO Cox stepwise regression analysis and further validated in independent AML cohorts. Patients were divided into high- and low-risk groups, and decreased survival rates were observed in high- vs. low-risk groups. The TIME between high- and low-risk groups has a similar scenery in C1 vs. C2 subtypes. Single-cell-level analysis verified that LSP1 and CRIP1 were upregulated in AML and exhausted CD8+ T cells. Dual targeting of these two markers might present a promising immunotherapeutic for AML. In addition, potential effective chemical drugs for AML were predicted. Thus, we concluded that molecular subtyping using ferroptosis signatures could characterize the TIME and provide implications for monitoring clinical outcomes and predicting novel therapies.