CEACAM6 as a Machine Learning Derived Immune Biomarker for Predicting Neoadjuvant Chemotherapy Response in HR+/HER2-Breast Cancer

Dalang Fang^1*Jie Lin¹Wang Jin¹Qingxiao Nong¹Shouwen Tao¹Bimin Lu¹Yanrong Yu¹Hao Peng¹Yingying Tian¹Qunying Su¹Yanfei Ma¹Yuanlu Huang²

• ¹Affiliated Hospital of Youjiang Medical University for Nationalities, Baise, China
• ²Baise People's Hospital, Baise, China

Background: Hormone receptor-positive/human epidermal growth factor receptor 2negative (HR+/HER2-) breast cancer is the most common subtype, characterized by heterogeneous responses to neoadjuvant chemotherapy (NAC) and a low pathological complete response (pCR) rate. Existing biomarkers have limited predictive accuracy, hindering personalized treatment. This study aimed to identify predictive biomarkers for NAC response and explore their therapeutic potential in HR+/HER2-breast cancer. Methods: We integrated 497 HR+/HER2-samples from TCGA and 956 from nine GEO datasets (training set: n=708; test set: n=248). Differentially expressed genes (DEGs) between tumors and normal tissues (TCGA) and between pCR and residual disease (RD) groups (GEO) were identified. Overlapping DEGs were further screened using LASSO, random forest, and SVM-RFE algorithms. Predictive models were constructed with 10 machine learning algorithms and interpreted using SHAP. Gene set enrichment analysis (GSEA), CIBERSORT-based immune infiltration, and drug sensitivity prediction using oncoPredict and GDSC2 were performed. Immunohistochemistry (IHC) was conducted on paired pre/post-NAC samples (n=9). Clinical correlation was analyzed in a retrospective cohort of 106 HR+/HER2-NAC patients. Results: Thirty-eight overlapping DEGs were identified, and four key genes (CEACAM6, MELK, RARRES1, BIRC5) were selected. NeuralNet showed the best model performance (AUC=0.816). CEACAM6 was the most important SHAP feature; its high expression predicted RD and was associated with poor survival (p=0.014). GSEA revealed CEACAM6-high tumors were enriched in drug resistance pathways (e.g., oxidative phosphorylation), while low expression correlated with immune activation. Immune analysis showed pCR tumors had more effector cells (Tfh, γδ T cells, M1 macrophages), whereas RD tumors were enriched in Tregs and resting mast cells. CEACAM6 positively correlated with Tregs and naïve CD4+ T cells, and negatively with CD8+ T cells and M1 macrophages. CEACAM6-high tumors had higher IC50s for six NAC-related drugs. IHC confirmed persistent CEACAM6 expression in RD tumors post-NAC. Clinically, pCR patients had higher lymphocyte counts and more frequent N2-N3 nodal status. Conclusion: CEACAM6 is a promising predictive biomarker in HR+/HER2-breast cancer, associated with chemoresistance and immune suppression. Machine learning models incorporating immune and pathway features may enhance personalized NAC strategies.