Abstract
Breast cancer is now the leading cause of cancer morbidity and mortality among women worldwide. Paclitaxel and anthracycline-based neoadjuvant chemotherapy is widely used for the treatment of breast cancer, but its sensitivity remains difficult to predict for clinical use. In our study, a LASSO logistic regression method was applied to develop a genomic classifier for predicting pathologic complete response (pCR) to neoadjuvant chemotherapy in breast cancer. The predictive accuracy of the signature classifier was further evaluated using four other independent test sets. Also, functional enrichment analysis of genes in the signature was performed, and the correlations between the prediction score of the signature classifier and immune characteristics were explored. We found a 25-gene signature classifier through the modeling, which showed a strong ability to predict pCR to neoadjuvant chemotherapy in breast cancer. For T/FAC-based training and test sets, and a T/AC-based test set, the AUC of the signature classifier is 1.0, 0.9071, 0.9683, 0.9151, and 0.7350, respectively, indicating that it has good predictive ability for both T/FAC and T/AC schemes. The multivariate model showed that 25-gene signature was far superior to other clinical parameters as independent predictor. Functional enrichment analysis indicated that genes in the signature are mainly enriched in immune-related biological processes. The prediction score of the classifier was significantly positively correlated with the immune score. There were also significant differences in immune cell types between pCR and residual disease (RD) samples. Conclusively, we developed a 25-gene signature classifier that can effectively predict pCR to paclitaxel and anthracycline-based neoadjuvant chemotherapy in breast cancer. Our study also suggests that the immune ecosystem is actively involved in modulating clinical response to neoadjuvant chemotherapy and is beneficial to patient outcomes.