In ER+/HER2- breast cancer, multiple measures of intra-tumor heterogeneity are associated with worse response to endocrine therapy. To investigate heterogeneity in response to treatment, we developed an operating room-to-laboratory pipeline for the collection of live human tumors and normal breast specimens immediately after surgical resection for processing into single-cell workflows for experimentation and genomic analyses. We demonstrate differences in tamoxifen response by cell type and identify distinctly responsive and resistant subpopulations within the malignant cell compartment of human tumors. Tamoxifen resistance signatures from 3 distinct resistant subpopulations are prognostic in large cohorts of ER+ breast cancer patients and enriched in endocrine therapy resistant tumors. This novel ex vivo model system now provides a foundation to define responsive and resistant sub-populations within heterogeneous tumors, to develop precise single cell-based predictors of response to therapy, and to identify genes and pathways driving resistance to therapy.

The ErbB lineage of oncogenic receptor tyrosine kinases is frequently overexpressed in head and neck squamous cell carcinomas. A common co-regulon triggered by the ErbB proteins; involving shared signaling circuitries; may harbor co-druggable targets or response biomarkers for potential future multimodal precision therapy in ErbB-driven head and neck squamous cell carcinoma. We here present a cohort-based; genome-wide analysis of 488 head and neck squamous cell carcinomas curated as part of The Cancer Genome Atlas Project to characterize genes that are significantly positively co-regulated with the four ErbB proteins and those that are shared among all ErbBs denoting a common ErbB co-regulon. Significant positive gene correlations involved hundreds of genes that were co-expressed with the four ErbB family members (q < 0.05). A common; overlapping co-regulon consisted of a core set of 268 genes that were uniformly co-regulated with all four ErbB genes and highly enriched for functions in chromatin organization and histone modifications. This high-priority set of genes contained ten putative antineoplastic drug-gene interactions. The nature and directionality of these ten drug-gene associations was an inhibiting interaction for seven (PIK3CB; PIK3C2B; HDAC4; FRK; PRKCE; EPHA4; and DYRK1A) of them in which the drug decreases the biological activity or expression of the gene target. For three (CHD4; ARID1A; and PBRM1) of the associations; the directionality of the interaction was such that the gene predicted sensitivit y to the drug suggesting utility as potential response biomarkers. Drug-gene interactions that predicted the gene product to be reduced by the drug included a variety of potential targeted molecular agent classes. This unbiased genome-wide analysis identified a target-rich environment for multimodal therapeutic approaches in tumors that are putatively ErbB-driven. The results of this study require preclinical validation before ultimately devising lines of combinatorial treatment strategies for ErbB-dependent head and neck squamous cell carcinomas that incorporate these findings.