Abstract Background: Neoadjuvant pertuzumab and trastuzumab in combination with chemotherapy can achieve pathological complete response (pCR, defined as ypT0/is ypN0) in approximately 40–60% of patients with human epidermal growth factor receptor 2-positive (HER2+) breast cancer (BC). Patients who achieve pCR have significantly improved survival. Therefore, the prediction of pCR is crucial for optimizing neoadjuvant therapy. Methods: To identify a gene expression signature that could predict pCR in patients with HER2+ BC, we used the nCounter Breast Cancer 360TM V2 panel to quantify the expression of 758 genes in 104 HER2+ BC samples from patients who received neoadjuvant docetaxel, carboplatin, trastuzumab, and pertuzumab (TCHP) chemotherapy and subsequent curative breast surgery. Data normalization, differentially expressed gene (DEG) analysis, and gene set enrichment analysis (GSEA) were performed using the NanoTube Bioconductor package. Using the elastic net logistic regression model with optimal hyperparameters determined by 5-fold cross-validation (5-CV), we developed a chemosensitivity score based on the expression of significant DEGs (genes with absolute log2-fold-change ≥0.5 and q-value < 0.1). We tested its predictive value for pCR using multivariable logistic regression analysis. Results: The median age at diagnosis was 49.5 years. The clinical stage was II in 53 (51%) patients and III in 51 (49%) patients at diagnosis. Ninety (86.5%) patients had a HER2 immunohistochemistry (IHC) score of 3+, while the remaining 14 (13.5%) patients had HER2 IHC 2+ but HER2 amplification detected by in situ hybridization. Pathological evaluation of surgical specimens revealed that 55 (52.9%) patients achieved pCR. Menopausal status (p = 0.01), PAM50 intrinsic subtype (p = 0.002), hormone receptor (HR) status (p < 0.001), and HER2 IHC score (p = 0.019) were significantly correlated with the pCR rate. DEG analysis using the PAM50 intrinsic subtype, HR status, and HER2 IHC score as covariates revealed 6 downregulated DEGs and 33 upregulated DEGs in patients who achieved pCR compared to patients who did not (Table 1). The 39 DEGs were significantly enriched in nine preselected MSigDB hallmark gene sets representing immunological processes (p < 0.001). The optimal elastic net logistic regression model containing 39 DEGs achieved the mean 5-CV area under the receiver operating characteristic curve of 0.81. The chemosensitivity score calculated using the optimal model parameters was significantly correlated with pCR after adjustment for the PAM50 intrinsic subtype, HR status, and HER2 IHC score (p < 0.001). GSEA revealed eight hallmark gene sets that were significantly up- or downregulated in samples with pCR, including epithelial-mesenchymal transition (normalized enrichment score [NES] = -2.16; q < 0.001), late estrogen response (NES = -2.01; q = 0.004), early estrogen response (NES = -1.972; q = 0.005), and inflammatory response (NES = 1.53; q = 0.06). Conclusion: Our findings demonstrate that gene expression patterns can predict the response to neoadjuvant TCHP chemotherapy in HER2+ BC. The functional characterization of 39 DEGs suggests that the tumor microenvironment, including adjacent immune cells, plays a significant role in modulating the response to neoadjuvant chemotherapy in HER2+ BC. Table 1. Significant DEGs Citation Format: Junghoon Shin, Ji-Yeon Kim, Eun Yoon Cho, Seok Won Kim, Jeong Eon Lee, Hyunwoo Lee, Yeon Hee Park, Jin Seok Ahn, Young-Hyuck Im. Gene expression signature as a predictor of pathological complete response to neoadjuvant docetaxel, carboplatin, trastuzumab, and pertuzumab treatment in locally advanced HER2-positive breast cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO1-15-05.
Read full abstract