Abstract Molecular mechanisms underlying breast cancer lymph node metastasis remain unclear. Using single-cell sequencing, we investigated the transcriptome profile of 96,796 single cells from 15 paired samples of primary tumors and axillary lymph nodes. We identified nine cancer cell subclusters including CD44 + / ALDH2 + /ALDH6A1 + breast cancer stem cells (BCSCs), which had a copy-number variants profile similar to that of normal breast tissue. Importantly, BCSCs existed only in primary tumors and evolved into metastatic clusters infiltrating into lymph nodes. Furthermore, transcriptome data suggested that NECTIN2-TIGIT-mediated interactions between metastatic breast cancer cells and tumor microenvironment (TME) cells, which promoted immune escape and lymph node metastasis. This study is the first to delineate the transcriptome profile of breast cancer lymph node metastasis using single-cell RNA sequencing. Our findings offer novel insights into the mechanisms underlying breast cancer metastasis and have implications in developing novel therapies to inhibit the initiation of breast cancer metastasis.
