# Exploring the relationship between metabolism and immune microenvironment in breast cancer bone metastasis based on metabolic pathways

**Authors:** Changyou Yang, Zhaofeng Li, Guojian Li, Houlin Mi, Yi Qin

PMC · DOI: 10.1371/journal.pone.0341270 · 2026-01-29

## TL;DR

This study explores how metabolism and the immune environment interact in breast cancer bone metastasis, identifying key pathways and genes that affect patient outcomes.

## Contribution

The study introduces novel metabolic pathway subtypes and risk models that predict prognosis in breast cancer bone metastasis.

## Key findings

- Three metabolic pathways (Amino Acid, Cofactor/Vitamin, and Secondary Metabolite) are significantly linked to patient prognosis.
- Two subtypes (C1 and C2) and three gene clusters (MRG1–3) show distinct prognostic outcomes and immune infiltration differences.
- A metabolism-based risk model accurately predicts prognosis, and scRNA-seq reveals ALDH1A1 and macrophages as key players in BCBM.

## Abstract

Bone metastasis is a significant contributor to mortality in patients with advanced breast cancer. Its progression is deeply intertwined with tumor metabolic reprogramming and the remodeling of the immune microenvironment. However, the dynamic interplay between metabolic pathways and immune regulation remains incompletely elucidated.

In this study, leveraging RNA-seq data and clinical information from breast cancer bone metastasis (BCBM) patients sourced from the GEO database, integrated bioinformatic analyses were employed to determine the activity of metabolic pathways significantly associated with prognosis. Metabolism-related genes were identified, and different metabolism-related gene clusters (MRGs) were subsequently identified by unsupervised clustering. Furthermore, a risk model was constructed based on hub prognostic genes, and differences in immune cell infiltration and drug sensitivity were compared between different subgroups. Finally, through single-cell RNA sequencing (scRNA-seq) analysis, we elucidated cellular heterogeneity and cell-cell communication within the tumor microenvironment (TME).

This study identified three metabolic pathways (Amino Acid, Cofactor/Vitamin, and Secondary Metabolite) significantly associated with patient prognosis. Two metabolic pathway-related subtypes (C1 and C2) were defined, which exhibited differing prognostic outcomes. Concurrently, MRG1–3 were also identified, and there were significant differences in prognosis and immune infiltration levels between the three clusters, with MRG2 having a significantly better prognosis than MRG1 and MRG3. In addition, metabolism-related risk models based on risk scores were developed. The risk model had strong prognostic predictive ability. Subsequently, scRNA-seq analysis revealed that ALDH1A1 and macrophages may play a key role in BCBM.

This study reveals the prognostic metabolic pathways and important prognostic target genes in BCBM from the perspective of metabolism-immunity interaction. MRGs can well distinguish the prognosis of different patients, and metabolism-related risk modeling can be used as a good prognostic predictor, which provides valuable insights into the “metabolic-immune” perspective of treatment.

## Linked entities

- **Genes:** ALDH1A1 (aldehyde dehydrogenase 1 family member A1) [NCBI Gene 216]
- **Diseases:** breast cancer (MONDO:0004989)

## Full-text entities

- **Genes:** IL2RB (interleukin 2 receptor subunit beta) [NCBI Gene 3560] {aka CD122, IL15RB, IMD63, P70-75}, MEIS3 (Meis homeobox 3) [NCBI Gene 56917] {aka MRG2}, TLR4 (toll like receptor 4) [NCBI Gene 7099] {aka ARMD10, CD284, TLR-4, TOLL}, CD163 (CD163 molecule) [NCBI Gene 9332] {aka M130, MM130, SCARI1}, IL7 (interleukin 7) [NCBI Gene 3574] {aka IL-7, IMD130}, MIF (macrophage migration inhibitory factor) [NCBI Gene 4282] {aka GIF, GLIF, MMIF}, CTLA4 (cytotoxic T-lymphocyte associated protein 4) [NCBI Gene 1493] {aka ALPS5, CD, CD152, CELIAC3, CTLA-4, GRD4}, TDG (thymine DNA glycosylase) [NCBI Gene 6996] {aka hTDG}, BCL2 (BCL2 apoptosis regulator) [NCBI Gene 596] {aka Bcl-2, PPP1R50}, GRM1 (glutamate metabotropic receptor 1) [NCBI Gene 2911] {aka GPRC1A, MGLU1, MGLUR1, PPP1R85, SCA44, SCAR13}, PGR (progesterone receptor) [NCBI Gene 5241] {aka NR3C3, PR}, PTHLH (parathyroid hormone like hormone) [NCBI Gene 5744] {aka BDE2, HHM, PLP, PTHR, PTHRP}, TLR2 (toll like receptor 2) [NCBI Gene 7097] {aka CD282, TIL4}, MAP3K7 (mitogen-activated protein kinase kinase kinase 7) [NCBI Gene 6885] {aka CSCF, FMD2, MEKK7, TAK1, TGF1a}, RBM39 (RNA binding motif protein 39) [NCBI Gene 9584] {aka CAPER, CAPERalpha, FSAP59, HCC1, RNPC2}, TGFB1 (transforming growth factor beta 1) [NCBI Gene 7040] {aka CAEND1, CED, DPD1, IBDIMDE, LAP, TGF-beta1}, CDK9 (cyclin dependent kinase 9) [NCBI Gene 1025] {aka C-2k, CDC2L4, CTK1, PITALRE, TAK}, IL10 (interleukin 10) [NCBI Gene 3586] {aka CSIF, GVHDS, IL-10, IL10A, TGIF}, ASNS (asparagine synthetase (glutamine-hydrolyzing)) [NCBI Gene 440] {aka ASNSD, TS11}, CSF3 (colony stimulating factor 3) [NCBI Gene 1440] {aka C17orf33, CSF3OS, GCSF}, MAS1L (MAS1 proto-oncogene like, G protein-coupled receptor) [NCBI Gene 116511] {aka MAS-L, MRG, dJ994E9.2}, CSF2 (colony stimulating factor 2) [NCBI Gene 1437] {aka CSF, GMCSF}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, MRC1 (mannose receptor C-type 1) [NCBI Gene 4360] {aka CD206, CLEC13D, CLEC13DL, MMR, MRC1L1, bA541I19.1}, CCR7 (C-C motif chemokine receptor 7) [NCBI Gene 1236] {aka BLR2, CC-CKR-7, CCR-7, CD197, CDw197, CMKBR7}, MSR1 (macrophage scavenger receptor 1) [NCBI Gene 4481] {aka CD204, SCARA1, SR-A, SR-AI, SR-AII, SR-AIII}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, TLE3 (TLE family member 3, transcriptional corepressor) [NCBI Gene 7090] {aka ESG, ESG3, GRG3, HsT18976}, CD44 (CD44 molecule (IN blood group)) [NCBI Gene 960] {aka CDW44, CSPG8, ECM-III, ECMR-III, H-CAM, HCELL}, CD80 (CD80 molecule) [NCBI Gene 941] {aka B7, B7-1, B7.1, BB1, CD28LG, CD28LG1}, CXCR4 (C-X-C motif chemokine receptor 4) [NCBI Gene 7852] {aka CD184, D2S201E, FB22, HM89, HSY3RR, LCR1}, CITED4 (Cbp/p300 interacting transactivator with ED-rich tail 4) [NCBI Gene 163732], ITGAM (integrin subunit alpha M) [NCBI Gene 3684] {aka CD11B, CR3A, HNA-4, MAC-1, MAC1A, MO1A}, MORF4L1P1 (mortality factor 4 like 1 pseudogene 1) [NCBI Gene 326591] {aka MORF4LP1, MRG1}, CD4 (CD4 molecule) [NCBI Gene 920] {aka CD4mut, IMD79, Leu-3, OKT4D, T4}, TNFSF11 (TNF superfamily member 11) [NCBI Gene 8600] {aka CD254, ODF, OPGL, OPTB2, RANKL, TNLG6B}, CD86 (CD86 molecule) [NCBI Gene 942] {aka B7-2, B7.2, B70, BU63, CD28LG2, CD86 v6}, PHGDH (phosphoglycerate dehydrogenase) [NCBI Gene 26227] {aka 3-PGDH, 3PGDH, HEL-S-113, NLS, NLS1, PDG}, CD74 (CD74 molecule) [NCBI Gene 972] {aka CLIP, DHLAG, HLADG, II, Ia-GAMMA, p33}, ADAM12 (ADAM metallopeptidase domain 12) [NCBI Gene 8038] {aka ADAM12-OT1, CAR10, MCMP, MCMPMltna, MLTN, MLTNA}, ALDH1A1 (aldehyde dehydrogenase 1 family member A1) [NCBI Gene 216] {aka ALDC, ALDH-E1, ALDH1, ALDH11, HEL-9, HEL-S-53e}, CXCL8 (C-X-C motif chemokine ligand 8) [NCBI Gene 3576] {aka GCP-1, GCP1, IL8, LECT, LUCT, LYNAP}, TNF (tumor necrosis factor) [NCBI Gene 7124] {aka DIF, IMD127, TNF-alpha, TNFA, TNFSF2, TNLG1F}
- **Diseases:** osteosarcoma (MESH:D012516), bone metastatic cancers (MESH:D001859), hypoxic (MESH:D002534), hypercalcemia (MESH:D006934), lung cancer (MESH:D008175), Cancer (MESH:D009369), liver cancer (MESH:D006528), osteolytic (MESH:D030981), MRGs (MESH:D003027), pathological fractures (MESH:D005598), MDSCs (OMIM:601308), Bone metastasis (MESH:D009362), BCBM (MESH:D001943), -negative (MESH:D064726)
- **Chemicals:** Sabutoclax (MESH:C550162), Amino Acid (MESH:D000596), glucose (MESH:D005947), NO (MESH:D009614), MG-132 (MESH:C072553), polyamines (MESH:D011073), glutamate (MESH:D018698), lactate (MESH:D019344), CDK.9.1 (-), glutamine (MESH:D005973), CTX (MESH:C552428), acid (MESH:D000143), L-serine (MESH:D012694), Staurosporine (MESH:D019311), NAD(P)+ (MESH:D009249)
- **Species:** Homo sapiens (human, species) [taxon 9606]

## Figures

9 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12854414/full.md

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