Populations of tumour cells display remarkable variability in almost every discernable phenotypic trait, including clinically important phenotypes such as ability to seed metastases and to survive therapy. This phenotypic diversity results from the integration of both genetic and non-genetic influences. Recent technological advances have improved the molecular understanding of cancers and the identification of targets for therapeutic interventions. However, it has become exceedingly apparent that the utility of profiles based on the analysis of tumours en masse is limited by intra-tumour genetic and epigenetic heterogeneity, as characteristics of the most abundant cell type might not necessarily predict the properties of mixed populations. In this Review, we discuss both genetic and non-genetic causes of phenotypic heterogeneity of tumour cells, with an emphasis on heritable phenotypes that serve as a substrate for clonal selection. We discuss the implications of intra-tumour heterogeneity in diagnostics and the development of therapeutic resistance.

Intratumor heterogeneity is a major clinical problem because tumor cell subtypes display variable sensitivity to therapeutics and may play different roles in progression. We previously characterized 2 cell populations in human breast tumors with distinct properties: CD44 + CD24 -cells that have stem cell-like characteristics, and CD44 -CD24 + cells that resemble more differentiated breast cancer cells. Here we identified 15 genes required for cell growth or proliferation in CD44 + CD24 -human breast cancer cells in a large-scale loss-of-function screen and found that inhibition of several of these (IL6, PTGIS, HAS1, CXCL3, and PFKFB3) reduced Stat3 activation. We found that the IL-6/JAK2/Stat3 pathway was preferentially active in CD44 + CD24 -breast cancer cells compared with other tumor cell types, and inhibition of JAK2 decreased their number and blocked growth of xenografts. Our results highlight the differences between distinct breast cancer cell types and identify targets such as JAK2 and Stat3 that may lead to more specific and effective breast cancer therapies.

SUMMARYCancers arise through a process of somatic evolution that can result in substantial sub-clonal heterogeneity within tumors. The mechanisms responsible for the coexistence of distinct sub-clones and the biological consequences of this coexistence remain poorly understood. Here we used a mouse xenograft model to investigate the impact of sub-clonal heterogeneity on tumor phenotypes and the competitive expansion of individual clones. We found that tumor growth can be driven by a minor cell subpopulation, which enhances the proliferation of all cells within a tumor by overcoming environmental constraints and yet can be outcompeted by faster proliferating competitors, resulting in tumor collapse. We then developed a mathematical modeling framework to identify the rules underlying the generation of intratumor clonal heterogeneity. We found that non-cell autonomous driving, together with clonal interference, stabilizes sub-clonal heterogeneity, thereby enabling inter-clonal interactions that can lead to new phenotypic traits.

Intratumor heterogeneity represents a major obstacle to effective cancer treatment and personalized medicine. However, investigators are now elucidating intratumor heterogeneity at the single-cell level due to improvements in technologies. Better understanding of the composition of tumors, and monitoring changes in cell populations during disease progression and treatment, will improve cancer diagnosis and therapeutic design. Measurements of intratumor heterogeneity may also be used as biomarkers to predict the risk of progression and therapeutic resistance. We summarize important considerations related to intratumor heterogeneity during tumor evolution. We also discuss experimental approaches that are commonly used to infer intratumor heterogeneity and describe how these methodologies can be translated into clinical practice.

The 8q24 gene desert contains risk loci for multiple epithelial cancers, including colon, breast, and prostate. Recent evidence suggests these risk loci contain enhancers. In this study, data are presented showing that each risk locus bears epigenetic marks consistent with enhancer elements and forms a long-range chromatin loop with the MYC proto-oncogene located several hundred kilobases telomeric and that these interactions are tissue-specific. We therefore propose that the 8q24 risk loci operate through a common mechanism-as tissue-specific enhancers of MYC.chromosome conformation capture | epigenetic | genetics

SUMMARY Cancer therapy exerts a strong selection pressure that shapes tumor evolution, yet our knowledge of how tumors change during treatment is limited. Here we report the analysis of cellular heterogeneity for genetic and phenotypic features and their spatial distribution in breast tumors pre- and post-neoadjuvant chemotherapy. We found that intratumor genetic diversity was tumor subtype-specific and it did not change during treatment in tumors with partial or no response. However, lower pre-treatment genetic diversity was significantly associated with complete pathologic response. In contrast, phenotypic diversity was different between pre- and post-treatment samples. We also observed significant changes in the spatial distribution of cells with distinct genetic and phenotypic features. We used these experimental data to develop a stochastic computational model to infer tumor growth patterns and evolutionary dynamics. Our results highlight the importance of integrated analysis of genotypes and phenotypes of single cells in intact tissues to predict tumor evolution.

SUMMARY Recurrent mutations in histone modifying enzymes imply key roles in tumorigenesis yet their functional relevance is largely unknown. Here we show that JARID1B, encoding a histone H3 lysine 4 (H3K4) demethylase, is frequently amplified and overexpressed in luminal breast tumors and a somatic mutation in a basal-like breast cancer results in the gain of unique chromatin binding and luminal expression and splicing patterns. Downregulation of JARID1B in luminal cells induces basal genes expression and growth arrest, which is rescued by TGFβ pathway inhibitors. Integrated JARID1B chromatin binding, H3K4 methylation, and expression profiles suggest a key function for JARID1B in luminal cell-specific expression programs. High luminal JARID1B activity is associated with poor outcome in patients with hormone receptor positive breast tumors.

BRCA1-associated breast tumors display loss of BRCA1 and frequent somatic mutations of PTEN and TP53. Here we describe the analysis of BRCA1, PTEN, and p53 at the single cell level in 55 BRCA1-associated breast tumors and computational methods to predict the relative temporal order of somatic events, on the basis of the frequency of cells with single or combined alterations. Although there is no obligatory order of events, we found that loss of PTEN is the most common first event and is associated with basal-like subtype, whereas in the majority of luminal tumors, mutation of TP53 occurs first and mutant PIK3CA is rarely detected. We also observed intratumor heterogeneity for the loss of wild-type BRCA1 and increased cell proliferation and centrosome amplification in the normal breast epithelium of BRCA1 mutation carriers. Our results have important implications for the design of chemopreventive and therapeutic interventions in this high-risk patient population. SIGNIFICANCE Defining the temporal order of tumor-driving somatic events is critical for early detection, risk stratification, and the design of chemopreventive therapies. Our combined experimental and computational approach reveal that the loss of wild-type BRCA1 may not be the first event in the majority of BRCA1-associated breast tumors and may not be present in all cancer cells within tumors.