# Multi-cohort and single-cell profiling of aging genes reveals prognostic and therapeutic targets in breast cancer

**Authors:** Li Huang, Lei Zhang, Xiaoyu Shi, Chun Wang, Xin Chen, Miao Li, Ni Ni, Ge Gao, Tao Wang, Xiaonan Zhang

PMC · DOI: 10.1016/j.isci.2026.114847 · 2026-01-29

## TL;DR

A machine learning model based on aging-related genes helps predict breast cancer outcomes and treatment responses.

## Contribution

A novel aging gene signature (MLAG) is developed that outperforms existing models in predicting breast cancer prognosis.

## Key findings

- The MLAG score consistently stratifies breast cancer patients by survival risk across multiple datasets.
- High MLAG tumors exhibit genomic instability and immunosuppressive features, while low MLAG tumors show enhanced immune infiltration.
- MLAG predicts immunotherapy response and identifies potential therapeutic strategies like panobinostat.

## Abstract

Aging-related transcriptional programs shape breast cancer progression, immune regulation, and therapeutic response. We integrated curated aging-associated genes with 108 survival modeling strategies to derive a machine learning-based aging gene signature. Trained in the TCGA-BRCA cohort and evaluated across twelve independent datasets, the eight-gene MLAG score consistently stratified patients by survival risk across clinical contexts. High MLAG exhibited increased genomic instability, including elevated tumor mutational burden and copy number alterations, together with distinct regulatory and intercellular communication patterns. In contrast, low MLAG showed enhanced immune infiltration, coordinated microenvironment signaling, and higher immune checkpoint expression. Integration of bulk and single-cell transcriptomic analyses localized MLAG-associated aging programs primarily to malignant epithelial cells and linked them to aneuploidy and stress-response pathways. The MLAG score associated with immune checkpoint blockade response and supported identification of candidate therapeutic strategies. These findings highlight the relevance of aging-associated transcriptional states in breast cancer prognosis and treatment stratification.

•A machine learning aging gene signature robustly stratifies breast cancer prognosis•The MLAG model outperforms 100 published signatures across multiple cohorts•High MLAG tumors show genomic instability and immunosuppressive microenvironments•MLAG predicts immunotherapy response and identifies panobinostat as a candidate therapy

A machine learning aging gene signature robustly stratifies breast cancer prognosis

The MLAG model outperforms 100 published signatures across multiple cohorts

High MLAG tumors show genomic instability and immunosuppressive microenvironments

MLAG predicts immunotherapy response and identifies panobinostat as a candidate therapy

Computational bioinformatics; Cancer; Omics

## Linked entities

- **Chemicals:** panobinostat (PubChem CID 6918837)
- **Diseases:** breast cancer (MONDO:0004989)

## Full-text entities

- **Genes:** GADD45A (growth arrest and DNA damage inducible alpha) [NCBI Gene 1647] {aka DDIT1, GADD45}, TTN (titin) [NCBI Gene 7273] {aka CMD1G, CMH9, CMPD4, CMYO5, CMYP5, EOMFC}, TWIST1 (twist family bHLH transcription factor 1) [NCBI Gene 7291] {aka ACS3, BPES2, BPES3, CRS, CRS1, CSO}, COL5A3 (collagen type V alpha 3 chain) [NCBI Gene 50509], ATM (ATM serine/threonine kinase) [NCBI Gene 472] {aka AT1, ATA, ATC, ATD, ATDC, ATE}, SPI1 (Spi-1 proto-oncogene) [NCBI Gene 6688] {aka AGM10, OF, PU.1, SFPI1, SPI-1, SPI-A}, ETV1 (ETS variant transcription factor 1) [NCBI Gene 2115] {aka ER81}, TLE2 (TLE family member 2, transcriptional corepressor) [NCBI Gene 7089] {aka ESG, ESG2, GRG2}, MYC (MYC proto-oncogene, bHLH transcription factor) [NCBI Gene 4609] {aka MRTL, MYCC, bHLHe39, c-Myc}, EREG (epiregulin) [NCBI Gene 2069] {aka EPR, ER, Ep}, PVT1 (Pvt1 oncogene) [NCBI Gene 5820] {aka LINC00079, MIR1204HG, NCRNA00079, TP53LC09, onco-lncRNA-100}, EPRS1 (glutamyl-prolyl-tRNA synthetase 1) [NCBI Gene 2058] {aka EARS, EPRS, GLUPRORS, HLD15, PARS, PIG32}, GSDMC (gasdermin C) [NCBI Gene 56169] {aka MLZE}, MMP14 (matrix metallopeptidase 14) [NCBI Gene 4323] {aka MMP-14, MMP-X1, MT-MMP, MT-MMP 1, MT1-MMP, MT1MMP}, CD99 (CD99 molecule (Xg blood group)) [NCBI Gene 4267] {aka HBA71, MIC2, MIC2X, MIC2Y, MSK5X}, NF1 (neurofibromin 1) [NCBI Gene 4763] {aka NFNS, VRNF, WSS}, F3 (coagulation factor III, tissue factor) [NCBI Gene 2152] {aka CD142, TF, TFA}, ZNF580 (zinc finger protein 580) [NCBI Gene 51157], CREBBP (CREB binding lysine acetyltransferase) [NCBI Gene 1387] {aka CBP, KAT3A, MKHK1, RSTS, RSTS1}, ABL2 (ABL proto-oncogene 2, non-receptor tyrosine kinase) [NCBI Gene 27] {aka ABLL, ARG}, HOXC10 (homeobox C10) [NCBI Gene 3226] {aka HOX3I}, MTOR (mechanistic target of rapamycin kinase) [NCBI Gene 2475] {aka FRAP, FRAP1, FRAP2, RAFT1, RAPT1, SKS}, MUC16 (mucin 16, cell surface associated) [NCBI Gene 94025] {aka CA125}, MITF (melanocyte inducing transcription factor) [NCBI Gene 4286] {aka CMM8, COMMAD, MI, MITF-A, WS2, WS2A}, CTLA4 (cytotoxic T-lymphocyte associated protein 4) [NCBI Gene 1493] {aka ALPS5, CD, CD152, CELIAC3, CTLA-4, GRD4}, INPP4B (inositol polyphosphate-4-phosphatase type II B) [NCBI Gene 8821], MDK (midkine) [NCBI Gene 4192] {aka ARAP, MK, NEGF2}, HOXA3 (homeobox A3) [NCBI Gene 3200] {aka HOX1, HOX1E}, LAG3 (lymphocyte activating 3) [NCBI Gene 3902] {aka CD223}, PDCD1 (programmed cell death 1) [NCBI Gene 5133] {aka ADMIO4, AIMTBS, CD279, PD-1, PD1, SLEB2}, HDAC9 (histone deacetylase 9) [NCBI Gene 9734] {aka HD7, HD7b, HD9, HDAC, HDAC7B, HDAC9B}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, CCND2 (cyclin D2) [NCBI Gene 894] {aka KIAK0002, MPPH3}, ZFR2 (zinc finger RNA binding protein 2) [NCBI Gene 23217] {aka KIAA1086}, FOXM1 (forkhead box M1) [NCBI Gene 2305] {aka FKHL16, FOXM1A, FOXM1B, FOXM1C, HFH-11, HFH11}, NR4A1 (nuclear receptor subfamily 4 group A member 1) [NCBI Gene 3164] {aka GFRP1, HMR, N10, NAK-1, NGFIB, NP10}, SDC4 (syndecan 4) [NCBI Gene 6385] {aka SYND4}, RUNX3 (RUNX family transcription factor 3) [NCBI Gene 864] {aka AML2, CBFA3, PEBP2aC}, PGR (progesterone receptor) [NCBI Gene 5241] {aka NR3C3, PR}, BCR (BCR activator of RhoGEF and GTPase) [NCBI Gene 613] {aka ALL, BCR1, CML, D22S11, D22S662, PHL}, SELE (selectin E) [NCBI Gene 6401] {aka CD62E, ELAM, ELAM1, ESEL, LECAM2, selectin-e}, ERBB2 (erb-b2 receptor tyrosine kinase 2) [NCBI Gene 2064] {aka CD340, HER-2, HER-2/neu, HER2, MLN 19, MLN-19}, HMOX1 (heme oxygenase 1) [NCBI Gene 3162] {aka HMOX1D, HO-1, HSP32, bK286B10}, PIK3CB (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit beta) [NCBI Gene 5291] {aka P110BETA, PI3K, PI3KBETA, PIK3C1}, ENO1 (enolase 1) [NCBI Gene 2023] {aka ENO1-IT1, ENO1L1, HEL-S-17, MPB1, NNE, PPH}, ATCAY (ATCAY kinesin light chain interacting caytaxin) [NCBI Gene 85300] {aka BNIP-H, CLAC}, MAP2K6 (mitogen-activated protein kinase kinase 6) [NCBI Gene 5608] {aka CRCMSL, MAPKK6, MEK6, MKK6, PRKMK6, SAPKK-3}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, CD4 (CD4 molecule) [NCBI Gene 920] {aka CD4mut, IMD79, Leu-3, OKT4D, T4}, CD274 (CD274 molecule) [NCBI Gene 29126] {aka ADMIO5, B7-H, B7H1, PD-L1, PDCD1L1, PDCD1LG1}, MLXIP (MLX interacting protein) [NCBI Gene 22877] {aka MIR, MONDOA, bHLHe36}, CCNA2 (cyclin A2) [NCBI Gene 890] {aka CCN1, CCNA}, CCDC26 (CCDC26 long non-coding RNA) [NCBI Gene 137196] {aka GLM7, RAM}, HMGB1 (high mobility group box 1) [NCBI Gene 3146] {aka HMG-1, HMG1, HMG3, SBP-1}, TBX2 (T-box transcription factor 2) [NCBI Gene 6909] {aka VETD}, DDR1 (discoidin domain receptor tyrosine kinase 1) [NCBI Gene 780] {aka CAK, CD167, DDR, EDDR1, HGK2, MCK10}, MAZ (MYC associated zinc finger protein) [NCBI Gene 4150] {aka PUR1, Pur-1, SAF-1, SAF-2, SAF-3, ZF87}, TP53 (tumor protein p53) [NCBI Gene 7157] {aka BCC7, BMFS5, LFS1, P53, TRP53}, RELB (RELB proto-oncogene, NF-kB subunit) [NCBI Gene 5971] {aka I-REL, IMD53, IREL, REL-B}, TJP3 (tight junction protein 3) [NCBI Gene 27134] {aka ZO-3, ZO3}
- **Diseases:** fibrosis (MESH:D005355), chronic inflammation (MESH:D007249), tumorigenicity (MESH:D002471), death (MESH:D003643), Tumor Immune Dysfunction and (MESH:D007154), Tumor (MESH:D009369), epithelial tumor (MESH:D002277), nodal (MESH:D013611), TNBC (MESH:D064726), Breast cancer (MESH:D001943), oncogenesis (MESH:D063646), Dysfunction (MESH:D006331), MLAG (MESH:D007859), hematologic malignancies (MESH:D019337), chronic (MESH:D002908), aneuploidy (MESH:D000782)
- **Chemicals:** H&amp;E (MESH:D006371), paraffin (MESH:D010232), hematoxylin (MESH:D006416), WZ-4002 (MESH:C571455), sapitinib (MESH:C548875), Panobinostat (MESH:D000077767), eosin (MESH:D004801), navitoclax (MESH:C528561), formalin (MESH:D005557), BI-2536 (MESH:C518477)
- **Species:** Homo sapiens (human, species) [taxon 9606]
- **Cell lines:** MCF7 — Homo sapiens (Human), Invasive breast carcinoma of no special type, Cancer cell line (CVCL_0031), HeLa — Homo sapiens (Human), Human papillomavirus-related endocervical adenocarcinoma, Cancer cell line (CVCL_0030)

## Figures

10 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12924738/full.md

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