Cancer metastasis requires that primary tumour cells evolve the capacity to intravasate into the lymphatic system or vasculature, and extravasate into and colonize secondary sites1. Others have demonstrated that individual cells within complex populations show heterogeneity in their capacity to form secondary lesions2–5. Here we develop a polyclonal mouse model of breast tumour heterogeneity, and show that distinct clones within a mixed population display specialization, for example, dominating the primary tumour, contributing to metastatic populations, or showing tropism for entering the lymphatic or vasculature systems. We correlate these stable properties to distinct gene expression profiles. Those clones that efficiently enter the vasculature express two secreted proteins, Serpine2 and Slpi, which were necessary and sufficient to program these cells for vascular mimicry. Our data indicate that these proteins not only drive the formation of extra-vascular networks but also ensure their perfusion by acting as anticoagulants. We propose that vascular mimicry drives the ability of some breast tumour cells to contribute to distant metastases while simultaneously satisfying a critical need of the primary tumour to be fed by the vasculature. Enforced expression of SERPINE2 and SLPI in human breast cancer cell lines also programmed them for vascular mimicry, and SERPINE2 and SLPI were overexpressed preferentially in human patients that had lung-metastatic relapse. Thus, these two secreted proteins, and the phenotype they promote, may be broadly relevant as drivers of metastatic progression in human cancer.

Using a functional model of breast cancer heterogeneity, we previously showed that clonal sub-populations proficient at generating circulating tumour cells were not all equally capable of forming metastases at secondary sites. A combination of differential expression and focused in vitro and in vivo RNA interference screens revealed candidate drivers of metastasis that discriminated metastatic clones. Among these, asparagine synthetase expression in a patient's primary tumour was most strongly correlated with later metastatic relapse. Here we show that asparagine bioavailability strongly influences metastatic potential. Limiting asparagine by knockdown of asparagine synthetase, treatment with l-asparaginase, or dietary asparagine restriction reduces metastasis without affecting growth of the primary tumour, whereas increased dietary asparagine or enforced asparagine synthetase expression promotes metastatic progression. Altering asparagine availability in vitro strongly influences invasive potential, which is correlated with an effect on proteins that promote the epithelial-to-mesenchymal transition. This provides at least one potential mechanism for how the bioavailability of a single amino acid could regulate metastatic progression.

Inappropriate activation of developmental pathways is a well-recognized tumor-promoting mechanism. Here we show that overexpression of the homeoprotein Six1, normally a developmentally restricted transcriptional regulator, increases TGF-β signaling in human breast cancer cells and induces an epithelial-mesenchymal transition (EMT) that is in part dependent on its ability to increase TGF-β signaling. TGF-β signaling and EMT have been implicated in metastatic dissemination of carcinoma. Accordingly, we used spontaneous and experimental metastasis mouse models to demonstrate that Six1 overexpression promotes breast cancer metastasis. In addition, we show that, like its induction of EMT, Six1-induced experimental metastasis is dependent on its ability to activate TGF-β signaling. Importantly, in human breast cancers Six1 correlated with nuclear Smad3 and thus increased TGF-β signaling. Further, breast cancer patients whose tumors overexpressed Six1 had a shortened time to relapse and metastasis and an overall decrease in survival. Finally, we show that the effects of Six1 on tumor progression likely extend beyond breast cancer, since its overexpression correlated with adverse outcomes in numerous other cancers including brain, cervical, prostate, colon, kidney, and liver. Our findings indicate that Six1, acting through TGF-β signaling and EMT, is a powerful and global promoter of cancer metastasis.

The ability to predict metastatic potential could be of great clinical importance, however, it is uncertain if predicting metastasis to specific vital organs is feasible. As a first step in evaluating metastatic predictions, we analyzed multiple primary tumors and metastasis pairs and determined that >90% of 298 gene expression signatures were found to be similarly expressed between matched pairs of tumors and metastases; therefore, primary tumors may be a good predictor of metastatic propensity. Next, using a dataset of >1,000 human breast tumor gene expression microarrays we determined that HER2-enriched subtype tumors aggressively spread to the liver, while basal-like and claudin-low subtypes colonize the brain and lung. Correspondingly, brain and lung metastasis signatures, along with embryonic stem cell, tumor initiating cell, and hypoxia signatures, were also strongly expressed in the basal-like and claudin-low tumors. Interestingly, low “Differentiation Scores,” or high expression of the aforementioned signatures, further predicted for brain and lung metastases. In total, these data identify that depending upon the organ of relapse, a combination of gene expression signatures most accurately predicts metastatic behavior.Electronic supplementary materialThe online version of this article (doi:10.1007/s10549-011-1619-7) contains supplementary material, which is available to authorized users.

BackgroundHuman breast cancer is a heterogeneous disease consisting of multiple molecular subtypes. Genetically engineered mouse models are a useful resource for studying mammary cancers in vivo under genetically controlled and immune competent conditions. Identifying murine models with conserved human tumor features will facilitate etiology determinations, highlight the effects of mutations on pathway activation, and should improve preclinical drug testing.ResultsTranscriptomic profiles of 27 murine models of mammary carcinoma and normal mammary tissue were determined using gene expression microarrays. Hierarchical clustering analysis identified 17 distinct murine subtypes. Cross-species analyses using three independent human breast cancer datasets identified eight murine classes that resemble specific human breast cancer subtypes. Multiple models were associated with human basal-like tumors including TgC3(1)-Tag, TgWAP-Myc and Trp53-/-. Interestingly, the TgWAPCre-Etv6 model mimicked the HER2-enriched subtype, a group of human tumors without a murine counterpart in previous comparative studies. Gene signature analysis identified hundreds of commonly expressed pathway signatures between linked mouse and human subtypes, highlighting potentially common genetic drivers of tumorigenesis.ConclusionsThis study of murine models of breast carcinoma encompasses the largest comprehensive genomic dataset to date to identify human-to-mouse disease subtype counterparts. Our approach illustrates the value of comparisons between species to identify murine models that faithfully mimic the human condition and indicates that multiple genetically engineered mouse models are needed to represent the diversity of human breast cancers. The reported trans-species associations should guide model selection during preclinical study design to ensure appropriate representatives of human disease subtypes are used.

Accurate classification is essential for understanding the pathophysiology of a disease and can inform therapeutic choices. For hematopoietic malignancies, a classification scheme based on the phenotypic similarity between tumor cells and normal cells has been successfully used to define tumor subtypes; however, use of normal cell types as a reference by which to classify solid tumors has not been widely emulated, in part due to more limited understanding of epithelial cell differentiation compared with hematopoiesis. To provide a better definition of the subtypes of epithelial cells comprising the breast epithelium, we performed a systematic analysis of a large set of breast epithelial markers in more than 15,000 normal breast cells, which identified 11 differentiation states for normal luminal cells. We then applied information from this analysis to classify human breast tumors based on normal cell types into 4 major subtypes, HR0-HR3, which were differentiated by vitamin D, androgen, and estrogen hormone receptor (HR) expression. Examination of 3,157 human breast tumors revealed that these HR subtypes were distinct from the current classification scheme, which is based on estrogen receptor, progesterone receptor, and human epidermal growth factor receptor 2. Patient outcomes were best when tumors expressed all 3 hormone receptors (subtype HR3) and worst when they expressed none of the receptors (subtype HR0). Together, these data provide an ontological classification scheme associated with patient survival differences and provides actionable insights for treating breast tumors.

Five molecular subtypes (luminal A, luminal B, HER2-enriched, basal-like, and claudin-low) with clinical implications exist in breast cancer. Here, we evaluated the molecular and phenotypic relationships of (1) a large in vitro panel of human breast cancer cell lines (BCCLs), human mammary fibroblasts (HMFs), and human mammary epithelial cells (HMECs); (2) in vivo breast tumors; (3) normal breast cell subpopulations; (4) human embryonic stem cells (hESCs); and (5) bone marrow-derived mesenchymal stem cells (hMSC). First, by integrating genomic data of 337 breast tumor samples with 93 cell lines we were able to identify all the intrinsic tumor subtypes in the cell lines, except for luminal A. Secondly, we observed that the cell lines recapitulate the differentiation hierarchy detected in the normal mammary gland, with claudin-low BCCLs and HMFs cells showing a stromal phenotype, HMECs showing a mammary stem cell/bipotent progenitor phenotype, basal-like cells showing a luminal progenitor phenotype, and luminal B cell lines showing a mature luminal phenotype. Thirdly, we identified basal-like and highly migratory claudin-low subpopulations of cells within a subset of triple-negative BCCLs (SUM149PT, HCC1143, and HCC38). Interestingly, both subpopulations within SUM149PT were enriched for tumor-initiating cells, but the basal-like subpopulation grew tumors faster than the claudin-low subpopulation. Finally, claudin-low BCCLs resembled the phenotype of hMSCs, whereas hESCs cells showed an epithelial phenotype without basal or luminal differentiation. The results presented here help to improve our understanding of the wide range of breast cancer cell line models through the appropriate pairing of cell lines with relevant in vivo tumor and normal cell counterparts.Electronic supplementary materialThe online version of this article (doi:10.1007/s10549-013-2743-3) contains supplementary material, which is available to authorized users.

Currently available human tumour cell line panels consist of a small number of lines in each lineage that generally fail to retain the phenotype of the original patient tumour. Here we develop a cell culture medium that enables us to routinely establish cell lines from diverse subtypes of human ovarian cancers with >95% efficiency. Importantly, the 25 new ovarian tumour cell lines described here retain the genomic landscape, histopathology and molecular features of the original tumours. Furthermore, the molecular profile and drug response of these cell lines correlate with distinct groups of primary tumours with different outcomes. Thus, tumour cell lines derived using this methodology represent a significantly improved platform to study human tumour pathophysiology and response to therapy.