ARTICLE ABSTRACTThere is a need to individualize assays for tumor molecular phenotyping, given variations in the differentiation status of tumor and normal tissues in different patients. To address this, we performed single-cell genomics of breast tumors and adjacent normal cells propagated for a short duration under growth conditions that enable epithelial reprogramming. Cells analyzed were either unselected for a specific subpopulation or phenotypically defined as undifferentiated and highly clonogenic ALDH+/CD49f+/EpCAM+ luminal progenitors, which express both basal cell and luminal cell–enriched genes. We analyzed 420 tumor cells and 284 adjacent normal cells for expression of 93 genes that included a PAM50-intrinsic subtype classifier and stemness-related genes. ALDH+/CD49f+/EpCAM+ tumor and normal cells clustered differently compared with unselected tumor and normal cells. PAM50 gene-set analyses of ALDH+/CD49f+/EpCAM+ populations efficiently identified major and minor clones of tumor cells, with the major clone resembling clinical parameters of the tumor. Similarly, a stemness-associated gene set identified clones with divergent stemness pathway activation within the same tumor. This refined expression profiling technique distinguished genes truly deregulated in cancer from genes that identify cellular precursors of tumors. Collectively, the assays presented here enable more precise identification of cancer-deregulated genes, allow for early identification of therapeutically targetable tumor cell subpopulations, and ultimately provide a refinement of precision therapeutics for cancer treatment. Cancer Res; 77(10); 2759–69. ©2017 AACR.