Postoperative progression and chemotherapy resistance is the major cause of treatment failure in patients with triple-negative breast cancer (TNBC). Currently, there is a lack of an ideal predictive model for the progression and drug sensitivity of postoperative TNBC patients. Diverse programmed cell death (PCD) patterns play an important role in tumor progression, which has the potential to be a prognostic and drug sensitivity indicator for TNBC after surgery. Twelve PCD patterns (apoptosis, necroptosis, pyroptosis, ferroptosis, cuproptosis, entotic cell death, netotic cell death, parthanatos, lysosome-dependent cell death, autophagy-dependent cell death, alkaliptosis, and oxeiptosis) were analyzed for model construction. Bulk transcriptome, single-cell transcriptome, genomics, and clinical information were collected from TCGA-BRCA, METABRIC, GSE58812, GSE21653, GSE176078, GSE75688, and KM-plotter cohorts to validate the model. The machine learning algorithm established a cell death index (CDI) with a 12-gene signature. Validated in five independent datasets, TNBC patients with high CDI had a worse prognosis after surgery. Two molecular subtypes of TNBC with distinct vital biological processes were identified by an unsupervised clustering model. A nomogram with high predictive performance was constructed by incorporating CDI with clinical features. Furthermore, CDI was associated with immune checkpoint genes and key tumor microenvironment components by integrated analysis of bulk and single-cell transcriptome. TNBC patients with high CDI are resistant to standard adjuvant chemotherapy regimens (docetaxel, oxaliplatin, etc.); however, they might be sensitive to palbociclib (an FDA-approved drug for luminal breast cancer). Generally, we established a novel CDI model by comprehensively analyzing diverse cell death patterns, which can accurately predict clinical prognosis and drug sensitivity of TNBC after surgery. A user-friendly website was created to facilitate the application of this prediction model (https://tnbc.shinyapps.io/CDI_Model/).

Distant metastasis remains the major cause of morbidity for breast cancer. Individuals with liver or brain metastasis have an extremely poor prognosis and low response rates to anti‐PD‐1/L1 immune checkpoint therapy compared to those with metastasis at other sites. Therefore, it is urgent to investigate the underlying mechanism of anti‐PD‐1/L1 resistance and develop more effective immunotherapy strategies for these patients. Using single‐cell RNA sequencing, a high‐resolution map of the entire tumor ecosystem based on 44 473 cells from breast cancer liver and brain metastases is depicted. Identified by canonical markers and confirmed by multiplex immunofluorescent staining, the metastatic ecosystem features remarkable reprogramming of immunosuppressive cells such as FOXP3+ regulatory T cells, LAMP3+ tolerogenic dendritic cells, CCL18+ M2‐like macrophages, RGS5+ cancer‐associated fibroblasts, and LGALS1+ microglial cells. In addition, PD‐1 and PD‐L1/2 are barely expressed in CD8+ T cells and cancer/immune/stromal cells, respectively. Interactions of the immune checkpoint molecules LAG3‐LGALS3 and TIGIT‐NECTIN2 between CD8+ T cells and cancer/immune/stromal cells are found to play dominant roles in the immune escape. In summary, this study dissects the intratumoral heterogeneity and immunosuppressive microenvironment in liver and brain metastases of breast cancer for the first time, providing insights into the most appropriate immunotherapy strategies for these patients.

Intrinsic and acquired anti-HER2 resistance remains a major hurdle for treating HER2-positive breast cancer. Using genome-wide CRISPR/Cas9 screening in vitro and in vivo, we identify FGFR4 as an essential gene following anti-HER2 treatment. FGFR4 inhibition enhances susceptibility to anti-HER2 therapy in resistant breast cancer. Mechanistically, m6A-hypomethylation regulated FGFR4 phosphorylates GSK-3β and activates β-catenin/TCF4 signaling to drive anti-HER2 resistance. Notably, suppression of FGFR4 dramatically diminishes glutathione synthesis and Fe²⁺ efflux efficiency via the β-catenin/TCF4-SLC7A11/FPN1 axis, resulting in excessive ROS production and labile iron pool accumulation. Ferroptosis, a unique iron-dependent form of oxidative cell death, is triggered after FGFR4 inhibition. Experiments involving patient-derived xenografts and organoids reveals a synergistic effect of anti-FGFR4 with anti-HER2 therapy in breast cancer with either intrinsic or acquired resistance. Together, these results pinpoint a mechanism of anti-HER2 resistance and provide a strategy for overcoming resistance via FGFR4 inhibition in recalcitrant HER2-positive breast cancer.

This study aims to explore the function and mechanism of G Protein-coupled receptor class C group 5 member A (GPRC5A) in docetaxel-resistance and liver metastasis of breast cancer. Single-cell RNA transcriptomic analysis and bioinformatic analysis are used to screen relevant genes in breast cancer metastatic hepatic specimens. MeRIP, dual-luciferase analysis and bioinformation were used to detect m6A modulation. Mass spectrometry (MS), co-inmunoprecipitation (co-IP) and immunofluorescence colocalization were executed to explore the mechanism of GPRC5A in breast cancer cells. GPRC5A was upregulated in triple-negative breast cancer (TNBC) and was associated with a poor prognosis. In vitro and in vivo experiments demonstrated that knockdown of GPRC5A alleviated metastasis and resistance to docetaxel in TNBC. Overexpression of GPRC5A had the opposite effects. The m6A methylation of GPRC5A mRNA was modulated by METTL3 and YTHDF1, which facilitates its translation. GPRC5A inhibited the ubiquitination-dependent degradation of LAMTOR1, resulting in the recruitment of mTORC1 to lysosomes and activating the mTORC1/p70s6k signaling pathway. METTL3/YTHDF1 axis up-regulates GPRC5A expression by m6A methylation. GPRC5A activates mTORC1/p70s6k signaling pathway by recruiting mTORC1 to lysosomes, consequently promotes docetaxel-resistance and liver metastasis.

The heterogeneity of macrophages influences the response to immune checkpoint inhibitor (ICI) therapy. However, few studies explore the impact of APOE ⁺ macrophages on ICI therapy using single‐cell RNA sequencing (scRNA‐seq) and machine learning methods. The scRNA‐seq and bulk RNA‐seq data are Integrated to construct an M.Sig model for predicting ICI response based on the distinct molecular signatures of macrophage and machine learning algorithms. Comprehensive single‐cell analysis as well as in vivo and in vitro experiments are applied to explore the potential mechanisms of the APOE ⁺ macrophage in affecting ICI response. The M.Sig model shows clear advantages in predicting the efficacy and prognosis of ICI therapy in pan‐cancer patients. The proportion of APOE ⁺ macrophages is higher in ICI non‐responders of triple‐negative breast cancer compared with responders, and the interaction and longer distance between APOE ⁺ macrophages and CD8 ⁺ exhausted T (Tex) cells affecting ICI response is confirmed by multiplex immunohistochemistry. In a 4T1 tumor‐bearing mice model, the APOE inhibitor combined with ICI treatment shows the best efficacy. The M.Sig model using real‐world immunotherapy data accurately predicts the ICI response of pan‐cancer, which may be associated with the interaction between APOE ⁺ macrophages and CD8 ⁺ Tex cells.

Disruption of disulfide homeostasis during biological processes can have fatal consequences. Excess disulfides induce cell death in a novel manner, termed as “disulfidptosis.” However, the specific mechanism of disulfidptosis has not yet been elucidated. To determine the cancer types sensitive to disulfidptosis and outline the corresponding treatment strategies, we firstly investigated the crucial functions of disulfidptosis regulators pan‐cancer at multi‐omics levels. We found that different tumor types expressed dysregulated levels of disulfidptosis regulators, most of which had an impact on tumor prognosis. Moreover, we calculated the disulfidptosis activity score in tumors and validated it using multiple independent datasets. Additionally, we found that disulfidptosis activity was correlated with classic biological processes and pathways in various cancers. Disulfidptosis activity was also associated with tumor immune characteristics and could predict immunotherapy outcomes. Notably, the disulfidptosis regulator, glycogen synthase 1 ( GYS1 ), was identified as a promising target for triple‐negative breast cancer and validated via in vitro and in vivo experiments. In conclusion, our study elucidated the complex molecular phenotypes and clinicopathological correlations of disulfidptosis regulators in tumors, laying a solid foundation for the development of disulfidptosis‐targeting strategies for cancer treatment.