# Derivation of a novel multi‐gene prognostic model based on regulated cell death pathways in acute myeloid leukemia: A comprehensive bioinformatic analysis integrating gene expression, mutation profiling, and immune infiltration

**Authors:** Ali Ahmadi, Amir Abas Navidinia, Davood Bashash, Behzad Poopak, Shadi Esmaeili

PMC · DOI: 10.1371/journal.pone.0328412 · PLOS One · 2025-08-01

## TL;DR

A new four-gene model based on cell death pathways helps predict survival in acute myeloid leukemia patients and reveals genetic and immune patterns linked to prognosis.

## Contribution

A novel multi-gene prognostic model based on regulated cell death pathways in AML, validated across cohorts and linked to genomic and immune features.

## Key findings

- A four-gene model (ARHGEF35, GSN, ELANE, AKT3) accurately predicts AML patient survival with high AUC values.
- High-risk patients show enrichment of DNMT3A and RUNX1 mutations and increased M2 macrophage infiltration.
- Low-risk patients are associated with KIT mutations and reduced resting mast cells.

## Abstract

Acute myeloid leukemia (AML) is a highly aggressive hematologic malignancy with dismal survival outcomes, where dysregulation of regulated cell death (RCD) pathways plays a pivotal role in leukemogenesis and therapeutic resistance.

Differential expression analyses were performed comparing AML samples with healthy bone marrow. Diagnostic differentially expressed genes (DEGs) were then intersected with curated gene sets representing apoptosis, pyroptosis, autophagy, necroptosis, and ferroptosis to derive an RCD-based gene signature. Prognostic markers were identified by univariate Cox regression, and these markers were refined using LASSO regression to construct a multi-gene prognostic model that generated an individual risk score (RS) for each patient. The performance of the model was validated internally through Kaplan–Meier survival analyses and receiver operating characteristic (ROC) curves for 1-, 3-, and 5-year survival, and externally confirmed in an independent TARGET-AML cohort. In addition, mutation analysis was conducted using the maftools package, and immune infiltration profiling was performed with CIBERSORT and xCell to characterize the molecular landscape of the risk groups.

Our integrative approach yielded a four-gene prognostic model incorporating ARHGEF35, GSN, ELANE, and AKT3. High RS was strongly associated with adverse overall survival, with Kaplan–Meier analyses showing p-value < 0.0001 in the training cohort and p-value = 0.0026 in the testing cohort. The model demonstrated robust predictive accuracy with AUC values of 82%, 87%, and 91% for 1-, 3-, and 5-year survival in the training set, and 65%, 81%, and 94% in the testing set. Mutation analysis revealed that DNMT3A and RUNX1 mutations were significantly enriched in high-RS patients (p-value = 0.0015 and p-value = 0.0086, respectively), whereas KIT mutations were more prevalent in low-RS patients (p-value = 0.0058). Immune profiling indicated that high-RS patients had increased M2 macrophage infiltration (p-value = 0.0027) and reduced resting mast cells (p-value = 0.0033).

These findings establish that an RCD-based multi-gene risk model can robustly stratify AML patients by prognosis and illuminate underlying genomic and immunologic mechanisms, thereby offering promising avenues for personalized therapeutic strategies.

## Linked entities

- **Genes:** ARHGEF35 (Rho guanine nucleotide exchange factor 35) [NCBI Gene 445328], GSN (gelsolin) [NCBI Gene 2934], ELANE (elastase, neutrophil expressed) [NCBI Gene 1991], AKT3 (AKT serine/threonine kinase 3) [NCBI Gene 10000], DNMT3A (DNA methyltransferase 3 alpha) [NCBI Gene 1788], RUNX1 (RUNX family transcription factor 1) [NCBI Gene 861], KIT (KIT proto-oncogene, receptor tyrosine kinase) [NCBI Gene 3815]
- **Diseases:** acute myeloid leukemia (MONDO:0015667)

## Full-text entities

- **Genes:** EGFR (epidermal growth factor receptor) [NCBI Gene 1956] {aka ERBB, ERBB1, ERRP, HER1, NISBD2, NNCIS}, DNMT3A (DNA methyltransferase 3 alpha) [NCBI Gene 1788] {aka DNMT3A2, HESJAS, M.HsaIIIA, TBRS}, PTEN (phosphatase and tensin homolog) [NCBI Gene 5728] {aka 10q23del, BZS, CWS1, DEC, GLM2, MHAM}, RUNX1 (RUNX family transcription factor 1) [NCBI Gene 861] {aka AML1, AML1-EVI-1, AMLCR1, CBF2alpha, CBFA2, EVI-1}, IDH2 (isocitrate dehydrogenase (NADP(+)) 2) [NCBI Gene 3418] {aka D2HGA2, ICD-M, IDH, IDH-2, IDHM, IDP}, GSN (gelsolin) [NCBI Gene 2934] {aka ADF, AGEL, AMYLD4}, WT1 (WT1 transcription factor) [NCBI Gene 7490] {aka AWT1, GUD, NPHS4, WAGR, WIT-2, WT-1}, NPM1 (nucleophosmin 1) [NCBI Gene 4869] {aka B23, NPM}, H1-0 (H1.0 linker histone) [NCBI Gene 3005] {aka H1.0, H10, H1F0, H1FV}, ARHGEF35 (Rho guanine nucleotide exchange factor 35) [NCBI Gene 445328] {aka ARHGEF5L}, CEBPA (CCAAT enhancer binding protein alpha) [NCBI Gene 1050] {aka C/EBP-alpha, CEBP}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, BAX (BCL2 associated X, apoptosis regulator) [NCBI Gene 581] {aka BCL2L4}, AKT1 (AKT serine/threonine kinase 1) [NCBI Gene 207] {aka AKT, PKB, PKB-ALPHA, PRKBA, RAC, RAC-ALPHA}, FLT3 (fms related receptor tyrosine kinase 3) [NCBI Gene 2322] {aka CD135, FLK-2, FLK2, STK1}, TET2 (tet methylcytosine dioxygenase 2) [NCBI Gene 54790] {aka IMD75, KIAA1546, MDS}, BCL2 (BCL2 apoptosis regulator) [NCBI Gene 596] {aka Bcl-2, PPP1R50}, ELANE (elastase, neutrophil expressed) [NCBI Gene 1991] {aka ELA2, GE, HLE, HNE, NE, PMN-E}, TP53 (tumor protein p53) [NCBI Gene 7157] {aka BCC7, BMFS5, LFS1, P53, TRP53}, PIK3CB (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit beta) [NCBI Gene 5291] {aka P110BETA, PI3K, PI3KBETA, PIK3C1}, CD4 (CD4 molecule) [NCBI Gene 920] {aka CD4mut, IMD79, Leu-3, OKT4D, T4}, KIT (KIT proto-oncogene, receptor tyrosine kinase) [NCBI Gene 3815] {aka C-Kit, CD117, MASTC, PBT, SCFR}, BAK1 (BCL2 antagonist/killer 1) [NCBI Gene 578] {aka BAK, BAK-LIKE, BCL2L7, CDN1}, CEBPZ (CCAAT enhancer binding protein zeta) [NCBI Gene 10153] {aka CBF, CBF2, HSP-CBF, NOC1}, IDH1 (isocitrate dehydrogenase (NADP(+)) 1) [NCBI Gene 3417] {aka HEL-216, HEL-S-26, IDCD, IDH, IDP, IDPC}, AKT3 (AKT serine/threonine kinase 3) [NCBI Gene 10000] {aka MPPH, MPPH2, PKB-GAMMA, PKBG, PRKBG, RAC-PK-gamma}, MTOR (mechanistic target of rapamycin kinase) [NCBI Gene 2475] {aka FRAP, FRAP1, FRAP2, RAFT1, RAPT1, SKS}, FAS (Fas cell surface death receptor) [NCBI Gene 355] {aka ALPS1A, APO-1, APT1, CD95, FAS1, FASTM}, MIR497 (microRNA 497) [NCBI Gene 574456] {aka MIRN497, hsa-mir-497, mir-497}, MCL1 (MCL1 apoptosis regulator, BCL2 family member) [NCBI Gene 4170] {aka BCL2L3, EAT, MCL1-ES, MCL1L, MCL1S, Mcl-1}
- **Diseases:** RCD (MESH:D003643), sepsis (MESH:D018805), metastasis (MESH:D009362), hematologic malignancies (MESH:D019337), chronic lymphocytic leukemia (MESH:D015451), myelodysplastic syndromes (MESH:D009190), cytopenias (MESH:D006402), bladder cancer (MESH:D001749), Tumor (MESH:D009369), inflammatory (MESH:D007249), FAB (OMIM:176500), Leukemia (MESH:D007938), lung cancer (MESH:D008175), PDA (MESH:D004374), AML (MESH:D015470), gastric carcinoma (MESH:D013274), mitochondrial dysfunction (MESH:D028361)
- **Chemicals:** lipid (MESH:D008055), GDP (MESH:D006153), sphingolipid (MESH:D013107), ceramide (MESH:D002518), GTP (MESH:D006160), S1P (MESH:C060506), iron (MESH:D007501), calcium (MESH:D002118), PIP2 (MESH:D019269)
- **Species:** Homo sapiens (human, species) [taxon 9606]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12316299/full.md

## Figures

7 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12316299/full.md

## References

51 references — full list in the complete paper: https://tomesphere.com/paper/PMC12316299/full.md

---
Source: https://tomesphere.com/paper/PMC12316299