Abstract
Background: Immunogenic cell death (ICD) remodels the tumor immune microenvironment, plays an inherent role in tumor cell apoptosis, and promotes durable protective antitumor immunity. Currently, appropriate biomarker-based ICD immunotherapy for breast cancer (BC) is under active exploration. Methods: To determine the potential link between ICD genes and the clinical risk of BC, TCGA-BC was used as the training set and GSE58812 was used as the validation set. Gene expression, consistent clustering, enrichment analysis, and mutation omics analyses were performed to analyze the potential biological pathways of ICD genes involved in BC. Furthermore, a risk and prognosis model of ICD was constructed to evaluate the correlation between risk grade and immune infiltration, clinical stage, and survival prognosis. Results: We identified two ICD-related subtypes by consistent clustering and found that the C2 subtype was associated with good survival prognosis, abundant immune cell infiltration, and high activity of immune biological processes. Based on this, we constructed and validated an ICD risk and prognosis model of BC, including ATG5, HSP90AA1, PIK3CA, EIF2AK3, MYD88, IL1R1, and CD8A. This model can effectively predict the survival rate of patients with BC and is negatively correlated with the immune microenvironment and clinical stage. Conclusion: This study provides new insights into the role of ICD in BC. The novel classification risk model based on ICD in BC established in this study can aid in estimating the potential prognosis of patients with BC and the clinical outcomes of immunotherapy and postulates targets that are more useful in comprehensive treatment strategies.