Signatures of mutational processes in human cancer

All cancers are caused by somatic mutations; however, understanding of the biological processes generating these mutations is limited. The catalogue of somatic mutations from a cancer genome bears the signatures of the mutational processes that have been operative. Here we analysed 4,938,362 mutations from 7,042 cancers and extracted more than 20 distinct mutational signatures. Some are present in many cancer types, notably a signature attributed to the APOBEC family of cytidine deaminases, whereas others are confined to a single cancer class. Certain signatures are associated with age of the patient at cancer diagnosis, known mutagenic exposures or defects in DNA maintenance, but many are of cryptic origin. In addition to these genome-wide mutational signatures, hypermutation localized to small genomic regions, 'kataegis', is found in many cancer types. The results reveal the diversity of mutational processes underlying the development of cancer, with potential implications for understanding of cancer aetiology, prevention and therapy.

Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis

Metastasis is a multistep process during which cancer cells disseminate from the site of primary tumors and establish secondary tumors in distant organs. In a search for key regulators of metastasis in a murine breast tumor model, we have found that the transcription factor Twist, a master regulator of embryonic morphogenesis, plays an essential role in metastasis. Suppression of Twist expression in highly metastatic mammary carcinoma cells specifically inhibits their ability to metastasize from the mammary gland to the lung. Ectopic expression of Twist results in loss of E-cadherin-mediated cell-cell adhesion, activation of mesenchymal markers, and induction of cell motility, suggesting that Twist contributes to metastasis by promoting an epithelial-mesenchymal transition (EMT). In human breast cancers, high level of Twist expression is correlated with invasive lobular carcinoma, a highly infiltrating tumor type associated with loss of E-cadherin expression. These results establish a mechanistic link between Twist, EMT, and tumor metastasis.

Stromal fibroblasts present in invasive human breast carcinomas promote tumor growth and angiogenesis through elevated SDF-1/CXCL12 secretion

Fibroblasts often constitute the majority of the stromal cells within a breast carcinoma, yet the functional contributions of these cells to tumorigenesis are poorly understood. Using a coimplantation tumor xenograft model, we demonstrate that carcinoma-associated fibroblasts (CAFs) extracted from human breast carcinomas promote the growth of admixed breast carcinoma cells significantly more than do normal mammary fibroblasts derived from the same patients. The CAFs, which exhibit the traits of myofibroblasts, play a central role in promoting the growth of tumor cells through their ability to secrete stromal cell-derived factor 1 (SDF-1); CAFs promote angiogenesis by recruiting endothelial progenitor cells (EPCs) into carcinomas, an effect mediated in part by SDF-1. CAF-secreted SDF-1 also stimulates tumor growth directly, acting through the cognate receptor, CXCR4, which is expressed by carcinoma cells. Our findings indicate that fibroblasts within invasive breast carcinomas contribute to tumor promotion in large part through the secretion of SDF-1.

Mesenchymal stem cells within tumour stroma promote breast cancer metastasis

Mesenchymal stem cells have been recently described to localize to breast carcinomas, where they integrate into the tumour-associated stroma. However, the involvement of mesenchymal stem cells (or their derivatives) in tumour pathophysiology has not been addressed. Here, we demonstrate that bone-marrow-derived human mesenchymal stem cells, when mixed with otherwise weakly metastatic human breast carcinoma cells, cause the cancer cells to increase their metastatic potency greatly when this cell mixture is introduced into a subcutaneous site and allowed to form a tumour xenograft. The breast cancer cells stimulate de novo secretion of the chemokine CCL5 (also called RANTES) from mesenchymal stem cells, which then acts in a paracrine fashion on the cancer cells to enhance their motility, invasion and metastasis. This enhanced metastatic ability is reversible and is dependent on CCL5 signalling through the chemokine receptor CCR5. Collectively, these data demonstrate that the tumour microenvironment facilitates metastatic spread by eliciting reversible changes in the phenotype of cancer cells.

The evaluation of tumor-infiltrating lymphocytes (TILs) in breast cancer: recommendations by an International TILs Working Group 2014

BACKGROUND

The morphological evaluation of tumor-infiltrating lymphocytes (TILs) in breast cancer (BC) is gaining momentum as evidence strengthens for the clinical relevance of this immunological biomarker. Accumulating evidence suggests that the extent of lymphocytic infiltration in tumor tissue can be assessed as a major parameter by evaluation of hematoxylin and eosin (H&E)-stained tumor sections. TILs have been shown to provide prognostic and potentially predictive value, particularly in triple-negative and human epidermal growth factor receptor 2-overexpressing BC.

DESIGN

A standardized methodology for evaluating TILs is now needed as a prerequisite for integrating this parameter in standard histopathological practice, in a research setting as well as in clinical trials. This article reviews current data on the clinical validity and utility of TILs in BC in an effort to foster better knowledge and insight in this rapidly evolving field, and to develop a standardized methodology for visual assessment on H&E sections, acknowledging the future potential of molecular/multiplexed approaches.

CONCLUSIONS

The methodology provided is sufficiently detailed to offer a uniformly applied, pragmatic starting point and improve consistency and reproducibility in the measurement of TILs for future studies.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Landscape of somatic mutations in 560 breast cancer whole-genome sequences

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

Mutational processes molding the genomes of 21 breast cancers

All cancers carry somatic mutations. The patterns of mutation in cancer genomes reflect the DNA damage and repair processes to which cancer cells and their precursors have been exposed. To explore these mechanisms further, we generated catalogs of somatic mutation from 21 breast cancers and applied mathematical methods to extract mutational signatures of the underlying processes. Multiple distinct single- and double-nucleotide substitution signatures were discernible. Cancers with BRCA1 or BRCA2 mutations exhibited a characteristic combination of substitution mutation signatures and a distinctive profile of deletions. Complex relationships between somatic mutation prevalence and transcription were detected. A remarkable phenomenon of localized hypermutation, termed “kataegis,” was observed. Regions of kataegis differed between cancers but usually colocalized with somatic rearrangements. Base substitutions in these regions were almost exclusively of cytosine at TpC dinucleotides. The mechanisms underlying most of these mutational signatures are unknown. However, a role for the APOBEC family of cytidine deaminases is proposed.

PaperClip

The landscape of cancer genes and mutational processes in breast cancer

All cancers carry somatic mutations in their genomes. A subset, known as driver mutations, confer clonal selective advantage on cancer cells and are causally implicated in oncogenesis, and the remainder are passenger mutations. The driver mutations and mutational processes operative in breast cancer have not yet been comprehensively explored. Here we examine the genomes of 100 tumours for somatic copy number changes and mutations in the coding exons of protein-coding genes. The number of somatic mutations varied markedly between individual tumours. We found strong correlations between mutation number, age at which cancer was diagnosed and cancer histological grade, and observed multiple mutational signatures, including one present in about ten per cent of tumours characterized by numerous mutations of cytosine at TpC dinucleotides. Driver mutations were identified in several new cancer genes including AKT2, ARID1B, CASP8, CDKN1B, MAP3K1, MAP3K13, NCOR1, SMARCD1 and TBX3. Among the 100 tumours, we found driver mutations in at least 40 cancer genes and 73 different combinations of mutated cancer genes. The results highlight the substantial genetic diversity underlying this common disease.

The life history of 21 breast cancers

Cancer evolves dynamically as clonal expansions supersede one another driven by shifting selective pressures, mutational processes, and disrupted cancer genes. These processes mark the genome, such that a cancer's life history is encrypted in the somatic mutations present. We developed algorithms to decipher this narrative and applied them to 21 breast cancers. Mutational processes evolve across a cancer's lifespan, with many emerging late but contributing extensive genetic variation. Subclonal diversification is prominent, and most mutations are found in just a fraction of tumor cells. Every tumor has a dominant subclonal lineage, representing more than 50% of tumor cells. Minimal expansion of these subclones occurs until many hundreds to thousands of mutations have accumulated, implying the existence of long-lived, quiescent cell lineages capable of substantial proliferation upon acquisition of enabling genomic changes. Expansion of the dominant subclone to an appreciable mass may therefore represent the final rate-limiting step in a breast cancer's development, triggering diagnosis.

PaperClip

Molecular definition of breast tumor heterogeneity

Cells with distinct phenotypes including stem-cell-like properties have been proposed to exist in normal human mammary epithelium and breast carcinomas, but their detailed molecular characteristics and clinical significance are unclear. We determined gene expression and genetic profiles of cells purified from cancerous and normal breast tissue using markers previously associated with stem-cell-like properties. CD24+ and CD44+ cells from individual tumors were clonally related but not always identical. CD44+ cell-specific genes included many known stem-cell markers and correlated with decreased patient survival. The TGF-beta pathway was specifically active in CD44+ cancer cells, where its inhibition induced a more epithelial phenotype. Our data suggest prognostic relevance of CD44+ cells and therapeutic targeting of distinct tumor cell populations.

Molecular characterization of the tumor microenvironment in breast cancer

Here we describe the comprehensive gene expression profiles of each cell type composing normal breast tissue and in situ and invasive breast carcinomas using serial analysis of gene expression. Based on these data, we determined that extensive gene expression changes occur in all cell types during cancer progression and that a significant fraction of altered genes encode secreted proteins and receptors. Despite the dramatic gene expression changes in all cell types, genetic alterations were detected only in cancer epithelial cells. The CXCL14 and CXCL12 chemokines overexpressed in tumor myoepithelial cells and myofibroblasts, respectively, bind to receptors on epithelial cells and enhance their proliferation, migration, and invasion. Thus, chemokines may play a role in breast tumorigenesis by acting as paracrine factors.

Sequence analysis of mutations and translocations across breast cancer subtypes

Breast carcinoma is the leading cause of cancer-related mortality in women worldwide, with an estimated 1.38 million new cases and 458,000 deaths in 2008 alone. This malignancy represents a heterogeneous group of tumours with characteristic molecular features, prognosis and responses to available therapy. Recurrent somatic alterations in breast cancer have been described, including mutations and copy number alterations, notably ERBB2 amplifications, the first successful therapy target defined by a genomic aberration. Previous DNA sequencing studies of breast cancer genomes have revealed additional candidate mutations and gene rearrangements. Here we report the whole-exome sequences of DNA from 103 human breast cancers of diverse subtypes from patients in Mexico and Vietnam compared to matched-normal DNA, together with whole-genome sequences of 22 breast cancer/normal pairs. Beyond confirming recurrent somatic mutations in PIK3CA, TP53, AKT1, GATA3 and MAP3K1, we discovered recurrent mutations in the CBFB transcription factor gene and deletions of its partner RUNX1. Furthermore, we have identified a recurrent MAGI3-AKT3 fusion enriched in triple-negative breast cancer lacking oestrogen and progesterone receptors and ERBB2 expression. The MAGI3-AKT3 fusion leads to constitutive activation of AKT kinase, which is abolished by treatment with an ATP-competitive AKT small-molecule inhibitor.

Phosphoglycerate dehydrogenase diverts glycolytic flux and contributes to oncogenesis

Most tumors exhibit increased glucose metabolism to lactate, however, the extent to which glucose-derived metabolic fluxes are used for alternative processes is poorly understood. Using a metabolomics approach with isotope labeling, we found that in some cancer cells a relatively large amount of glycolytic carbon is diverted into serine and glycine metabolism through phosphoglycerate dehydrogenase (PHGDH). An analysis of human cancers showed that PHGDH is recurrently amplified in a genomic region of focal copy number gain most commonly found in melanoma. Decreasing PHGDH expression impaired proliferation in amplified cell lines. Increased expression was also associated with breast cancer subtypes, and ectopic expression of PHGDH in mammary epithelial cells disrupted acinar morphogenesis and induced other phenotypic alterations that may predispose cells to transformation. Our findings show that the diversion of glycolytic flux into a specific alternate pathway can be selected during tumor development and may contribute to the pathogenesis of human cancer.

Homologous Recombination Deficiency (HRD) Score Predicts Response to Platinum-Containing Neoadjuvant Chemotherapy in Patients with Triple-Negative Breast Cancer

PURPOSE

BRCA1/2-mutated and some sporadic triple-negative breast cancers (TNBC) have DNA repair defects and are sensitive to DNA-damaging therapeutics. Recently, three independent DNA-based measures of genomic instability were developed on the basis of loss of heterozygosity (LOH), telomeric allelic imbalance (TAI), and large-scale state transitions (LST).

EXPERIMENTAL DESIGN

We assessed a combined homologous recombination deficiency (HRD) score, an unweighted sum of LOH, TAI, and LST scores, in three neoadjuvant TNBC trials of platinum-containing therapy. We then tested the association of HR deficiency, defined as HRD score ≥42 or BRCA1/2 mutation, with response to platinum-based therapy.

RESULTS

In a trial of neoadjuvant platinum, gemcitabine, and iniparib, HR deficiency predicted residual cancer burden score of 0 or I (RCB 0/I) and pathologic complete response (pCR; OR = 4.96, P = 0.0036; OR = 6.52, P = 0.0058). HR deficiency remained a significant predictor of RCB 0/I when adjusted for clinical variables (OR = 5.86, P = 0.012). In two other trials of neoadjuvant cisplatin therapy, HR deficiency predicted RCB 0/I and pCR (OR = 10.18, P = 0.0011; OR = 17.00, P = 0.0066). In a multivariable model of RCB 0/I, HR deficiency retained significance when clinical variables were included (OR = 12.08, P = 0.0017). When restricted to BRCA1/2 nonmutated tumors, response was higher in patients with high HRD scores: RCB 0/I P = 0.062, pCR P = 0.063 in the neoadjuvant platinum, gemcitabine, and iniparib trial; RCB 0/I P = 0.0039, pCR P = 0.018 in the neoadjuvant cisplatin trials.

CONCLUSIONS

HR deficiency identifies TNBC tumors, including BRCA1/2 nonmutated tumors more likely to respond to platinum-containing therapy. Clin Cancer Res; 22(15); 3764-73. ©2016 AACR.

Efficacy of neoadjuvant Cisplatin in triple-negative breast cancer

PURPOSE Cisplatin is a chemotherapeutic agent not used routinely for breast cancer treatment. As a DNA cross-linking agent, cisplatin may be effective treatment for hereditary BRCA1-mutated breast cancers. Because sporadic triple-negative breast cancer (TNBC) and BRCA1-associated breast cancer share features suggesting common pathogenesis, we conducted a neoadjuvant trial of cisplatin in TNBC and explored specific biomarkers to identify predictors of response. PATIENTS AND METHODS Twenty-eight women with stage II or III breast cancers lacking estrogen and progesterone receptors and HER2/Neu (TNBC) were enrolled and treated with four cycles of cisplatin at 75 mg/m(2) every 21 days. After definitive surgery, patients received standard adjuvant chemotherapy and radiation therapy per their treating physicians. Clinical and pathologic treatment response were assessed, and pretreatment tumor samples were evaluated for selected biomarkers. Results Six (22%) of 28 patients achieved pathologic complete responses, including both patients with BRCA1 germline mutations;18 (64%) patients had a clinical complete or partial response. Fourteen (50%) patients showed good pathologic responses (Miller-Payne score of 3, 4, or 5), 10 had minor responses (Miller-Payne score of 1 or 2), and four (14%) progressed. All TNBCs clustered with reference basal-like tumors by hierarchical clustering. Factors associated with good cisplatin response include young age (P = .001), low BRCA1 mRNA expression (P = .03), BRCA1 promoter methylation (P = .04), p53 nonsense or frameshift mutations (P = .01), and a gene expression signature of E2F3 activation (P = .03). CONCLUSION Single-agent cisplatin induced response in a subset of patients with TNBC. Decreased BRCA1 expression may identify subsets of TNBCs that are cisplatin sensitive. Other biomarkers show promise in predicting cisplatin response.

The JAK2/STAT3 signaling pathway is required for growth of CD44⁺CD24⁻ stem cell-like breast cancer cells in human tumors

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.

Paracrine and autocrine signals induce and maintain mesenchymal and stem cell states in the breast

The epithelial-mesenchymal transition (EMT) has been associated with the acquisition of motility, invasiveness, and self-renewal traits. During both normal development and tumor pathogenesis, this change in cell phenotype is induced by contextual signals that epithelial cells receive from their microenvironment. The signals that are responsible for inducing an EMT and maintaining the resulting cellular state have been unclear. We describe three signaling pathways, involving transforming growth factor (TGF)-β and canonical and noncanonical Wnt signaling, that collaborate to induce activation of the EMT program and thereafter function in an autocrine fashion to maintain the resulting mesenchymal state. Downregulation of endogenously synthesized inhibitors of autocrine signals in epithelial cells enables the induction of the EMT program. Conversely, disruption of autocrine signaling by added inhibitors of these pathways inhibits migration and self-renewal in primary mammary epithelial cells and reduces tumorigenicity and metastasis by their transformed derivatives.

HRDetect is a predictor of BRCA1 and BRCA2 deficiency based on mutational signatures

Approximately 1-5% of breast cancers are attributed to inherited mutations in BRCA1 or BRCA2 and are selectively sensitive to poly(ADP-ribose) polymerase (PARP) inhibitors. In other cancer types, germline and/or somatic mutations in BRCA1 and/or BRCA2 (BRCA1/BRCA2) also confer selective sensitivity to PARP inhibitors. Thus, assays to detect BRCA1/BRCA2-deficient tumors have been sought. Recently, somatic substitution, insertion/deletion and rearrangement patterns, or 'mutational signatures', were associated with BRCA1/BRCA2 dysfunction. Herein we used a lasso logistic regression model to identify six distinguishing mutational signatures predictive of BRCA1/BRCA2 deficiency. A weighted model called HRDetect was developed to accurately detect BRCA1/BRCA2-deficient samples. HRDetect identifies BRCA1/BRCA2-deficient tumors with 98.7% sensitivity (area under the curve (AUC) = 0.98). Application of this model in a cohort of 560 individuals with breast cancer, of whom 22 were known to carry a germline BRCA1 or BRCA2 mutation, allowed us to identify an additional 22 tumors with somatic loss of BRCA1 or BRCA2 and 47 tumors with functional BRCA1/BRCA2 deficiency where no mutation was detected. We validated HRDetect on independent cohorts of breast, ovarian and pancreatic cancers and demonstrated its efficacy in alternative sequencing strategies. Integrating all of the classes of mutational signatures thus reveals a larger proportion of individuals with breast cancer harboring BRCA1/BRCA2 deficiency (up to 22%) than hitherto appreciated (∼1-5%) who could have selective therapeutic sensitivity to PARP inhibition.

A pleiotropically acting microRNA, miR-31, inhibits breast cancer metastasis

MicroRNAs are well suited to regulate tumor metastasis because of their capacity to coordinately repress numerous target genes, thereby potentially enabling their intervention at multiple steps of the invasion-metastasis cascade. We identify a microRNA exemplifying these attributes, miR-31, whose expression correlates inversely with metastasis in human breast cancer patients. Overexpression of miR-31 in otherwise-aggressive breast tumor cells suppresses metastasis. We deploy a stable microRNA sponge strategy to inhibit miR-31 in vivo; this allows otherwise-nonaggressive breast cancer cells to metastasize. These phenotypes do not involve confounding influences on primary tumor development and are specifically attributable to miR-31-mediated inhibition of several steps of metastasis, including local invasion, extravasation or initial survival at a distant site, and metastatic colonization. Such pleiotropy is achieved via coordinate repression of a cohort of metastasis-promoting genes, including RhoA. Indeed, RhoA re-expression partially reverses miR-31-imposed metastasis suppression. These findings indicate that miR-31 uses multiple mechanisms to oppose metastasis.

Complex landscapes of somatic rearrangement in human breast cancer genomes

Multiple somatic rearrangements are often found in cancer genomes; however, the underlying processes of rearrangement and their contribution to cancer development are poorly characterized. Here we use a paired-end sequencing strategy to identify somatic rearrangements in breast cancer genomes. There are more rearrangements in some breast cancers than previously appreciated. Rearrangements are more frequent over gene footprints and most are intrachromosomal. Multiple rearrangement architectures are present, but tandem duplications are particularly common in some cancers, perhaps reflecting a specific defect in DNA maintenance. Short overlapping sequences at most rearrangement junctions indicate that these have been mediated by non-homologous end-joining DNA repair, although varying sequence patterns indicate that multiple processes of this type are operative. Several expressed in-frame fusion genes were identified but none was recurrent. The study provides a new perspective on cancer genomes, highlighting the diversity of somatic rearrangements and their potential contribution to cancer development.

X chromosomal abnormalities in basal-like human breast cancer

Sporadic basal-like cancers (BLC) are a distinct class of human breast cancers that are phenotypically similar to BRCA1-associated cancers. Like BRCA1-deficient tumors, most BLC lack markers of a normal inactive X chromosome (Xi). Duplication of the active X chromosome and loss of Xi characterized almost half of BLC cases tested. Others contained biparental but nonheterochromatinized X chromosomes or gains of X chromosomal DNA. These abnormalities did not lead to a global increase in X chromosome transcription but were associated with overexpression of a small subset of X chromosomal genes. Other, equally aneuploid, but non-BLC rarely displayed these X chromosome abnormalities. These results suggest that X chromosome abnormalities contribute to the pathogenesis of BLC, both inherited and sporadic.

Subclonal diversification of primary breast cancer revealed by multiregion sequencing

The sequencing of cancer genomes may enable tailoring of therapeutics to the underlying biological abnormalities driving a particular patient's tumor. However, sequencing-based strategies rely heavily on representative sampling of tumors. To understand the subclonal structure of primary breast cancer, we applied whole-genome and targeted sequencing to multiple samples from each of 50 patients' tumors (303 samples in total). The extent of subclonal diversification varied among cases and followed spatial patterns. No strict temporal order was evident, with point mutations and rearrangements affecting the most common breast cancer genes, including PIK3CA, TP53, PTEN, BRCA2 and MYC, occurring early in some tumors and late in others. In 13 out of 50 cancers, potentially targetable mutations were subclonal. Landmarks of disease progression, such as resistance to chemotherapy and the acquisition of invasive or metastatic potential, arose within detectable subclones of antecedent lesions. These findings highlight the importance of including analyses of subclonal structure and tumor evolution in clinical trials of primary breast cancer.

Reconstruction of functionally normal and malignant human breast tissues in mice

The study of normal breast epithelial morphogenesis and carcinogenesis in vivo has largely used rodent models. Efforts at studying mammary morphogenesis and cancer with xenotransplanted human epithelial cells have failed to recapitulate the full extent of development seen in the human breast. We have developed an orthotopic xenograft model in which both the stromal and epithelial components of the reconstructed mammary gland are of human origin. Genetic modification of human stromal cells before the implantation of ostensibly normal human mammary epithelial cells resulted in the outgrowth of benign and malignant lesions. This experimental model allows for studies of human epithelial morphogenesis and differentiation in vivo and underscores the critical role of heterotypic interactions in human breast development and carcinogenesis.