BackgroundNatural killer (NK) cells play crucial roles in mediating anti-cancer activity in breast cancer (BRCA). However, the potential of NK cell-related molecules in predicting BRCA outcomes and guiding personalized therapy remains largely unexplored. This study focused on developing a prognostic and therapeutic prediction model for BRCA by incorporating NK cell-related genes.MethodsThe data analyzed primarily originated from the TCGA and GEO databases. The prognostic role of NK cells was evaluated, and marker genes of NK cells were identified via single-cell analysis. Module genes closely associated with immunotherapy resistance were identified by bulk transcriptome-based weighted correlation network analysis (WGCNA). Following taking intersection and LASSO regression, NK-related genes (NKRGs) relevant to BRCA prognosis were screened, and the NK-related prognostic signature was subsequently constructed. Analyses were further expanded to clinicopathological relevance, GSEA, tumor microenvironment (TME) analysis, immune function, immunotherapy responsiveness, and chemotherapeutics. Key NKRGs were screened by machine learning and validated by spatial transcriptomics (ST) and immunohistochemistry (IHC).ResultsTumor-infiltrating NK cells are a favorable prognostic factor in BRCA. By combining scRNA-seq and bulk transcriptomic analyses, we identified 7 NK-related prognostic NKRGs (CCL5, EFHD2, KLRB1, C1S, SOCS3, IRF1, and CCND2) and developed an NK-related risk scoring (NKRS) system. The prognostic reliability of NKRS was verified through survival and clinical relevance analyses across multiple cohorts. NKRS also demonstrated robust predictive power in various aspects, including TME landscape, immune functions, immunotherapy responses, and chemotherapeutic sensitivity. Additionally, KLRB1 and CCND2 emerged as key prognostic NKRGs identified through machine learning and external validation, with their expression correlation with NK cells confirmed in BRCA specimens by ST and IHC.ConclusionsWe developed a novel NK-related gene signature that has proven valuable for evaluating prognosis and treatment response in BRCA, expecting to advance precision medicine of BRCA.