Prior studies have highlighted that novel programmed cell death (PCD) modalities, including ferroptosis, pyroptosis, and necroptosis, are correlated with tumor progression and antitumor immunity. Nonetheless, comprehensive analysis of tumor microenvironment (TME) profiles mediated by the crosstalk of distinct PCD forms has not been conducted in breast cancer (BC).
Here, we curated 34 identified PCD-associated genes (PCDAGs) and applied the consensus clustering algorithm to establish PCD-mediated tumor patterns in BC. Subsequently, based on prognostic differentially expressed genes extracted from distinct PCD-mediated patterns, we applied the LASSO algorithm to construct CD_Score. Furthermore, the correlation analysis between CD_Score and TME features, molecular subtypes, clinicopathological characteristics, drug response, and immunotherapeutic efficacy was performed.
Three distinct PCD-clusters were determined among 2,038 BC samples, which did not only display different clinical outcomes but highly correlated to the established immunological tumor phenotypes: “desert,” “excluded,” and “inflamed” immune profiles. Based on the CD_Score derived from the PCD-related gene signature, BC patients could be stratified into CD_Score-low and -high group, of which the former displayed satisfactory survival outcome and enhanced immune infiltration. Further exploration identified that the CD_Score-high group significantly correlated with elevated neoantigen load and higher mutation frequency in SMGs (e.g., TP53 and MAP3K1) and reduced expression of immune checkpoint proteins.
This research is the first to emphasize the close relationship between distinct cell death modalities and the diversity and complexity of immune infiltration in TME. We established the CD_Score, which could help enhance our cognition of TME features and facilitate the clinical application of immunotherapy.