In order to verify the heterogeneity of BCSCs observed in mouse scRNA-seq data, we collected tumor samples from six breast cancer patients which were classified as luminal A, luminal B, TNBC and Her2+ subtypes and performed the scRNA-seq analysis (Table S7). The gene discussed is ERBB2; the disease is neoplasm.