Primary triple negative breast cancers (TNBC) represent approximately 16% of all breast cancers1 and are a tumour type defined by exclusion, for which comprehensive landscapes of somatic mutation have not been determined. Here we show in 104 early TNBC cases, that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some exhibiting only a handful of somatic aberrations in a few pathways, whereas others contain hundreds of somatic events and multiple pathways implicated. Integration with matched whole transcriptome sequence data revealed that only ~36% of mutations are expressed. By examining single nucleotide variant (SNV) allelic abundance derived from deep re-sequencing (median >20,000 fold) measurements in 2414 somatic mutations, we determine for the first time in an epithelial tumour, the relative abundance of clonal genotypes among cases in the population. We show that TNBC vary widely and continuously in their clonal frequencies at the time of diagnosis, with basal subtype TNBC2,3 exhibiting more variation than non-basal TNBC. Although p53 and PIK3CA/PTEN somatic mutations appear clonally dominant compared with other pathways, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal and cell shape/motility proteins occurred at lower clonal frequencies, suggesting they occurred later during tumour progression. Taken together our results show that future attempts to dissect the biology and therapeutic responses of TNBC will require the determination of individual tumour clonal genotypes.

Recent advances in next generation sequencing have made it possible to precisely characterize all somatic coding mutations that occur during the development and progression of individual cancers. Here we used these approaches to sequence the genomes (>43-fold coverage) and transcriptomes of an oestrogen-receptor-alpha-positive metastatic lobular breast cancer at depth. We found 32 somatic non-synonymous coding mutations present in the metastasis, and measured the frequency of these somatic mutations in DNA from the primary tumour of the same patient, which arose 9 years earlier. Five of the 32 mutations (in ABCB11, HAUS3, SLC24A4, SNX4 and PALB2) were prevalent in the DNA of the primary tumour removed at diagnosis 9 years earlier, six (in KIF1C, USP28, MYH8, MORC1, KIAA1468 and RNASEH2A) were present at lower frequencies (1-13%), 19 were not detected in the primary tumour, and two were undetermined. The combined analysis of genome and transcriptome data revealed two new RNA-editing events that recode the amino acid sequence of SRP9 and COG3. Taken together, our data show that single nucleotide mutational heterogeneity can be a property of low or intermediate grade primary breast cancers and that significant evolution can occur with disease progression.

The telomeric amplicon at 8p12 is common in oestrogen receptor-positive (ER+) breast cancers. Array-CGH and expression analyses of 1172 primary breast tumours revealed that ZNF703 was the single gene within the minimal amplicon and was amplified predominantly in the Luminal B subtype. Amplification was shown to correlate with increased gene and protein expression and was associated with a distinct expression signature and poor clinical outcome. ZNF703 transformed NIH 3T3 fibroblasts, behaving as a classical oncogene, and regulated proliferation in human luminal breast cancer cell lines and immortalized human mammary epithelial cells. Manipulation of ZNF703 expression in the luminal MCF7 cell line modified the effects of TGFβ on proliferation. Overexpression of ZNF703 in normal human breast epithelial cells enhanced the frequency of in vitro colony-forming cells from luminal progenitors. Taken together, these data strongly point to ZNF703 as a novel oncogene in Luminal B breast cancer.

No abstract