Abstract
Research indicated that Paclitaxel (PTX) can induce immunogenic cell death (ICD) through immunogenic modulation. However, the combination of PTX and ICD has not been extensively studied in breast cancer (BRCA). The TCGA-BRCA and GSE20685 datasets were enrolled in this study. Samples from the TCGA-BRCA dataset were consistently clustered based on selected immunogenic cell death-related genes (ICD-RGs). Next, candidate genes were obtained by overlapping differentially expressed genes (DEGs) between BRCA and normal groups, intersecting genes common to DEGs between cluster1 and cluster2 and hub module genes, and target genes of PTX from five databases. The univariate Cox algorithm and the least absolute shrinkage and selection operator (LASSO) were performed to obtain biomarkers and build a risk model. Following observing the immune microenvironment in differential risk subgroups, single-gene gene set enrichment analysis (GSEA) was carried out in all biomarkers. Finally, the expression of biomarkers was analyzed. Enrichment analysis showed that 626 intersecting genes were linked with inflammatory response. Further five biomarkers (CHI3L1, IL18, PAPLN, SH2D2A, and UBE2L6) were identified and a risk model was built. The model's performance was validated using GSE20685 dataset. Furthermore, the biomarkers were enriched with adaptive immune response. Lastly, the experimental results indicated that the alterations in IL18, SH2D2A, and CHI3L1 expression after treatment matched those in the public database. In this study, Five PTX-ICD-related biomarkers (CHI3L1, IL18, PAPLN, SH2D2A, and UBE2L6) were identified to aid in predicting BRCA treatment outcomes.