The initiation and progression of breast cancer from the transformation of the normal epithelium to ductal carcinoma in situ (DCIS) and invasive disease is a complex process involving the acquisition of genetic alterations and changes in gene expression, alongside microenvironmental and recognized histological alterations. Here, we sought to comprehensively characterise the genomic and transcriptomic features of the MCF10 isogenic model of breast cancer progression, and to functionally validate potential driver alterations in three-dimensional (3D) spheroids that may provide insights into breast cancer progression, and identify targetable alterations in conditions more similar to those encountered in vivo. We performed whole genome, exome and RNA sequencing of the MCF10 progression series to catalogue the copy number and mutational and transcriptomic landscapes associated with progression. We identified a number of predicted driver mutations (including PIK3CA and TP53) that were acquired during transformation of non-malignant MCF10A cells to their malignant counterparts that are also present in analysed primary breast cancers from The Cancer Genome Atlas (TCGA). Acquisition of genomic alterations identified MYC amplification and previously undescribed RAB3GAP1-HRAS and UBA2-PDCD2L expressed in-frame fusion genes in malignant cells. Comparison of pathway aberrations associated with progression showed that, when cells are grown as 3D spheroids, they show perturbations of cancer-relevant pathways. Functional interrogation of the dependency on predicted driver events identified alterations in HRAS, PIK3CA and TP53 that selectively decreased cell growth and were associated with progression from preinvasive to invasive disease only when cells were grown as spheroids. Our results have identified changes in the genomic repertoire in cell lines representative of the stages of breast cancer progression, and demonstrate that genetic dependencies can be uncovered when cells are grown in conditions more like those in vivo. The MCF10 progression series therefore represents a good model with which to dissect potential biomarkers and to evaluate therapeutic targets involved in the progression of breast cancer. © 2016 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

The intra-tumor diversity of cancer cells is under intense investigation; however, little is known about the heterogeneity of the tumor microenvironment that is key to cancer progression and evolution. We aimed to assess the degree of microenvironmental heterogeneity in breast cancer and correlate this with genomic and clinical parameters.

Metaplastic breast cancer (MBC) is a rare special histologic type of triple-negative breast cancer, characterized by the presence of neoplastic cells showing differentiation towards squamous epithelium and/or mesenchymal elements. Here we sought to define whether histologically distinct subgroups of MBCs would be underpinned by distinct genomic and/or transcriptomic alterations. Microarray-based copy number profiling identified limited but significant differences between the distinct MBC subtypes studied here, despite the limited sample size (n = 17). In particular, we found that, compared to MBCs with chondroid or squamous cell metaplasia, MBCs with spindle cell differentiation less frequently harbored gain of 7q11.22-23 encompassing CLDN3 and CLDN4, consistent with their lower expression of claudins and their association with the claudin-low molecular classification. Microarray-based and RNA-sequencing-based gene expression profiling revealed that MBCs with spindle cell differentiation differ from MBCs with chondroid or squamous cell metaplasia on the expression of epithelial-to-mesenchymal transition-related genes, including down-regulation of CDH1 and EPCAM. In addition, RNA-sequencing revealed that the histologic patterns observed in MBCs are unlikely to be underpinned by a highly recurrent expressed fusion gene or a pathognomonic expressed mutation in cancer genes. Loss of PTEN expression or mutations affecting PIK3CA or TSC2 observed in 8/17 MBCs support the contention that PI3K pathway activation plays a role in the development of MBCs. Our data demonstrate that despite harboring largely similar patterns of gene copy number alterations, MBCs with spindle cell, chondroid and squamous differentiation are distinct at the transcriptomic level but are unlikely to be defined by specific pathognomonic genetic alterations.

Triple-negative breast cancers (TNBC) are resistant to standard-of-care chemotherapy and lack known targetable driver gene alterations. Identification of novel drivers could aid the discovery of new treatment strategies for this hard-to-treat patient population, yet studies using high-throughput and accurate models to define the functions of driver genes in TNBC to date have been limited. Here, we employed unbiased functional genomics screening of the 200 most frequently mutated genes in breast cancer, using spheroid cultures to model in vivo-like conditions, and identified the histone acetyltransferase CREBBP as a novel tumor suppressor in TNBC. CREBBP protein expression in patient tumor samples was absent in 8% of TNBCs and at a high frequency in other tumors, including squamous lung cancer, where CREBBP-inactivating mutations are common. In TNBC, CREBBP alterations were associated with higher genomic heterogeneity and poorer patient survival and resulted in upregulation and dependency on a FOXM1 proliferative program. Targeting FOXM1-driven proliferation indirectly with clinical CDK4/6 inhibitors (CDK4/6i) selectively impaired growth in spheroids, cell line xenografts, and patient-derived models from multiple tumor types with CREBBP mutations or loss of protein expression. In conclusion, we have identified CREBBP as a novel driver in aggressive TNBC and identified an associated genetic vulnerability in tumor cells with alterations in CREBBP and provide a preclinical rationale for assessing CREBBP alterations as a biomarker of CDK4/6i response in a new patient population. SIGNIFICANCE: This study demonstrates that CREBBP genomic alterations drive aggressive TNBC, lung cancer, and lymphomas and may be selectively treated with clinical CDK4/6 inhibitors.

Breast cancers can be classified by hierarchical clustering using an "intrinsic" gene list into one of at least five molecular subtypes: basal-like, HER2, luminal A, luminal B, and normal breast-like. Five different intrinsic gene lists composed of varying numbers of genes have been used for molecular subtype identification and classification of breast cancers. The aim of this study was to determine the objectivity and interobserver reproducibility of the assignment of molecular subtype classes by hierarchical cluster analysis.

Genome-wide association studies have identified more than 70 common variants that are associated with breast cancer risk. Most of these variants map to non-protein-coding regions and several map to gene deserts, regions of several hundred kilobases lacking protein-coding genes. We hypothesized that gene deserts harbor long-range regulatory elements that can physically interact with target genes to influence their expression. To test this, we developed Capture Hi-C (CHi-C), which, by incorporating a sequence capture step into a Hi-C protocol, allows high-resolution analysis of targeted regions of the genome. We used CHi-C to investigate long-range interactions at three breast cancer gene deserts mapping to 2q35, 8q24.21, and 9q31.2. We identified interaction peaks between putative regulatory elements ("bait fragments") within the captured regions and "targets" that included both protein-coding genes and long noncoding (lnc) RNAs over distances of 6.6 kb to 2.6 Mb. Target protein-coding genes were IGFBP5, KLF4, NSMCE2, and MYC; and target lncRNAs included DIRC3, PVT1, and CCDC26. For one gene desert, we were able to define two SNPs (rs12613955 and rs4442975) that were highly correlated with the published risk variant and that mapped within the bait end of an interaction peak. In vivo ChIP-qPCR data show that one of these, rs4442975, affects the binding of FOXA1 and implicate this SNP as a putative functional variant.

The ability to utilize formalin-fixed, paraffin-embedded (FFPE) archival specimens reliably for high-resolution molecular genetic analysis would be of immense practical application in the study of human disease. We have evaluated the ability of the GenomePlex whole genome amplification (WGA) kit to amplify frozen and FFPE tissue for use in array CGH (aCGH). GenomePlex gave highly representative data compared with unamplified controls both from frozen material (Pearson's R(2) = 0.898) and from FFPE (R(2) = 0.883). Artifactual amplification observed using DOP-PCR at chromosomes 1p, 3, 13q, and 16p was not seen with GenomePlex. Highly reproducible aCGH profiles were obtained using as little as 5 ng starting material from FFPE (R(2) = 0.918). This WGA method should readily lend itself to the determination of DNA copy number alterations from small fresh-frozen and FFPE clinical tumor specimens, although some care must be taken to optimize the DNA extraction procedure.

Understanding how breast cancer (BC) grows in axillary lymph nodes (ALNs), and refining how therapies might halt that process, is clinically important. However, modelling the complex ALN microenvironment is difficult, and no human models exist at present. We harvested ALNs from ten BC patients, and perfused them at 37 °C ex vivo for up to 24 h. Controlled autologous testing showed that ALNs remain viable after 24 h of ex vivo perfusion: haematoxylin and eosin-stained histological appearance and proliferation (by Ki67 immunohistochemistry) did not change significantly over time for any perfused ALN compared with a control from time-point zero. Furthermore, targeted gene expression analysis (NanoString PanCancer IO360 panel) showed that only 21/750 genes were differentially expressed between control and perfused ALNs (|log2 FC| > 1 and q < 0.1): none were involved in apoptosis and metabolism, but rather all 21 genes were involved in immune function and angiogenesis. During perfusion, tissue acid-base balance remained stable. Interestingly, the flow rate increased (p < 0.001) in cancer-replaced (i.e. metastasis occupied more than 90% of the surface area on multiple levels) compared to cancer-free nodes (i.e. nodes with no metastasis on multiple sections). CXCL11 transcripts were significantly more abundant in cancer-replaced nodes, while CXCL12 transcripts were significantly more abundant in cancer-free nodes. These cytokines were also detected in the circulating perfusate. Monoclonal antibodies (nivolumab and trastuzumab) were administered into a further three ALNs to confirm perfusion efficacy. These drugs saturated the nodes; nivolumab even induced cancer cell death. Normothermic ALN perfusion is not only feasible but sustains the tumour microenvironment ex vivo for scientific investigation. This model could facilitate the identification of actionable immuno-oncology targets. © 2019 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Using genome-wide clustered regularly interspaced short palindromic repeats (CRISPR) screens to understand endocrine drug resistance, we discovered ARID1A and other SWI/SNF complex components as the factors most critically required for response to two classes of estrogen receptor-alpha (ER) antagonists. In this context, SWI/SNF-specific gene deletion resulted in drug resistance. Unexpectedly, ARID1A was also the top candidate in regard to response to the bromodomain and extraterminal domain inhibitor JQ1, but in the opposite direction, with loss of ARID1A sensitizing breast cancer cells to bromodomain and extraterminal domain inhibition. We show that ARID1A is a repressor that binds chromatin at ER cis-regulatory elements. However, ARID1A elicits repressive activity in an enhancer-specific, but forkhead box A1-dependent and active, ER-independent manner. Deletion of ARID1A resulted in loss of histone deacetylase 1 binding, increased histone 4 lysine acetylation and subsequent BRD4-driven transcription and growth. ARID1A mutations are more frequent in treatment-resistant disease, and our findings provide mechanistic insight into this process while revealing rational treatment strategies for these patients.

Ovarian clear cell carcinoma (OCCC) is a rare subtype of epithelial ovarian cancer characterised by a high frequency of loss-of-function ARID1A mutations and a poor response to chemotherapy. Despite their generally low mutational burden, an intratumoural T cell response has been reported in a subset of OCCC, with ARID1A purported to be a biomarker for the response to the immune checkpoint blockade independent of micro-satellite instability (MSI). However, assessment of the different immune cell types and spatial distribution specifically within OCCC patients has not been described to date. Here, we characterised the immune landscape of OCCC by profiling a cohort of 33 microsatellite stable OCCCs at the genomic, gene expression and histological level using targeted sequencing, gene expression profiling using the NanoString targeted immune panel, and multiplex immunofluorescence to assess the spatial distribution and abundance of immune cell populations at the protein level. Analysis of these tumours and subsequent independent validation identified an immune-related gene expression signature associated with risk of recurrence of OCCC. Whilst histological quantification of tumour-infiltrating lymphocytes (TIL, Salgado scoring) showed no association with the risk of recurrence or ARID1A mutational status, the characterisation of TILs via multiplexed immunofluorescence identified spatial differences in immunosuppressive cell populations in OCCC. Tumour-associated macrophages (TAM) and regulatory T cells were excluded from the vicinity of tumour cells in low-risk patients, suggesting that high-risk patients have a more immunosuppressive microenvironment. We also found that TAMs and cytotoxic T cells were also excluded from the vicinity of tumour cells in ARID1A-mutated OCCCs compared to ARID1A wild-type tumours, suggesting that the exclusion of these immune effectors could determine the host response of ARID1A-mutant OCCCs to therapy. Overall, our study has provided new insights into the immune landscape and prognostic associations in OCCC and suggest that tailored immunotherapeutic approaches may be warranted for different subgroups of OCCC patients.

In breast cancer (BC), detecting low volumes of axillary lymph node (ALN) metastasis pre-operatively is difficult and novel biomarkers are needed. We recently showed that patient-derived ALNs can be sustained ex-vivo using normothermic perfusion. We now compare reactive (tumour-free; n = 5) and macrometastatic (containing tumour deposits >2 mm; n = 4) ALNs by combining whole section multiplex immunofluorescence with TMT-labelled LC-MS/MS of the circulating perfusate. Macrometastases contained significantly fewer B cells and T cells (CD4+/CD8+/regulatory) than reactive nodes (p = 0.02). Similarly, pathway analysis of the perfusate proteome (119/1453 proteins significantly differentially expressed) showed that immune function was diminished in macrometastases in favour of 'extracellular matrix degradation'; only 'neutrophil degranulation' was preserved. Qualitative comparison of the perfusate proteome to that of node-positive pancreatic and prostatic adenocarcinoma also highlighted 'neutrophil degranulation' as a contributing factor to nodal metastasis. Thus, metastasis-induced changes in the REPLICANT perfusate proteome are detectable, and could facilitate biomarker discovery.

Breast cancers with BRCA1 germline mutation have a characteristic DNA copy number (CN) pattern. We developed a test that assigns CN profiles to be 'BRCA1-like' or 'non-BRCA1-like', which refers to resembling a BRCA1-mutated tumor or resembling a tumor without a BRCA1 mutation, respectively. Approximately one third of the BRCA1-like breast cancers have a BRCA1 mutation, one third has hypermethylation of the BRCA1 promoter and one third has an unknown reason for being BRCA1-like. This classification is indicative of patients' response to high dose alkylating and platinum containing chemotherapy regimens, which targets the inability of BRCA1 deficient cells to repair DNA double strand breaks. We investigated whether this classification can be reliably obtained with next generation sequencing and copy number platforms other than the bacterial artificial chromosome (BAC) array Comparative Genomic Hybridization (aCGH) on which it was originally developed. We investigated samples from 230 breast cancer patients for which a CN profile had been generated on two to five platforms, comprising low coverage CN sequencing, CN extraction from targeted sequencing panels (CopywriteR), Affymetrix SNP6.0, 135K/720K oligonucleotide aCGH, Affymetrix Oncoscan FFPE (MIP) technology, 3K BAC and 32K BAC aCGH. Pairwise comparison of genomic position-mapped profiles from the original aCGH platform and other platforms revealed concordance. For most cases, biological differences between samples exceeded the differences between platforms within one sample. We observed the same classification across different platforms in over 80% of the patients and kappa values of at least 0.36. Differential classification could be attributed to CN profiles that were not strongly associated to one class. In conclusion, we have shown that the genomic regions that define our BRCA1-like classifier are robustly measured by different CN profiling technologies, providing the possibility to retro- and prospectively investigate BRCA1-like classification across a wide range of CN platforms.

Poly(ADP-ribose) polymerase (PARP) is an enzyme involved in DNA repair. PARP inhibitors can act as chemosensitizers, or operate on the principle of synthetic lethality when used as single agent. Clinical trials have shown drugs in this class to be promising for BRCA mutation carriers. We postulated that inability to demonstrate response in non-BRCA carriers in which BRCA is inactivated by other mechanisms or with deficiency in homologous recombination for DNA repair is due to lack of molecular markers that define a responding subpopulation. We identified candidate markers for this purpose for olaparib (AstraZeneca) by measuring inhibitory effects of nine concentrations of olaparib in 22 breast cancer cell lines and identifying features in transcriptional and genome copy number profiles that were significantly correlated with response. We emphasized in this discovery process genes involved in DNA repair. We found that the cell lines that were sensitive to olaparib had a significant lower copy number of BRCA1 compared to the resistant cell lines (p value 0.012). In addition, we discovered seven genes from DNA repair pathways whose transcriptional levels were associated with response. These included five genes (BRCA1, MRE11A, NBS1, TDG, and XPA) whose transcript levels were associated with resistance and two genes (CHEK2 and MK2) whose transcript levels were associated with sensitivity. We developed an algorithm to predict response using the seven-gene transcription levels and applied it to 1,846 invasive breast cancer samples from 8 U133A/plus 2 (Affymetrix) data sets and found that 8-21 % of patients would be predicted to be responsive to olaparib. A similar response frequency was predicted in 536 samples analyzed on an Agilent platform. Importantly, tumors predicted to respond were enriched in basal subtype tumors. Our studies support clinical evaluation of the utility of our seven-gene signature as a predictor of response to olaparib.

The majority of borderline ovarian tumours (BOTs) behave in a benign fashion, but some may show aggressive behavior. The reason behind this has not been elucidated. The epidermal growth factor receptor (EGFR) is known to contribute to cell survival signals as well as metastatic potential of some tumours. EGFR expression and gene status have not been thoroughly investigated in BOTs as it has in ovarian carcinomas. In this study we explore protein expression as well as gene mutations and amplifications of EGFR in BOTs in comparison to a subset of other epithelial ovarian tumours. We studied 85 tumours, including 61 BOTs, 10 low grade serous carcinomas (LGSCs), 9 high grade serous carcinomas (HGSCs) and 5 benign epithelial tumours. EGFR protein expression was studied using immunohistochemistry. Mutations were investigated by Sanger sequencing exons 18-21 of the tyrosine kinase domain of EGFR. Cases with comparatively higher protein expression were examined for gene amplification by chromogenic in situ hybridization. We also studied the tumours for KRAS and BRAF mutations. Immunohistochemistry results revealed both cytoplasmic and nuclear EGFR expression with variable degrees between tumours. The level of nuclear localization was relatively higher in BOTs and LGSCs as compared to HGSCs or benign tumours. The degree of nuclear expression of BOTs showed no significant difference from that in LGSCs (mean ranks 36.48, 33.05, respectively, p=0.625), but was significantly higher than in HGSCs (mean ranks: 38.88, 12.61 respectively, p< 0.001) and benign tumours (mean ranks: 35.18, 13.00 respectively, p= 0.010). Cytoplasmic expression level was higher in LGSCs. No EGFR gene mutations or amplification were identified, yet different polymorphisms were detected. Five different types of point mutations in the KRAS gene and the V600E BRAF mutation were detected exclusively in BOTs and LGSCs. Our study reports for the first time nuclear localization of EGFR in BOTs. The nuclear localization similarities between BOTs and LGSCs and not HGSCs support the hypothesis suggesting evolution of LGSCs from BOTs. We also confirm that EGFR mutations and amplifications are not molecular events in the pathogenesis of BOTs.

The 19q12 locus is amplified in a subgroup of oestrogen receptor (ER)-negative grade III breast cancers. This amplicon comprises nine genes, including cyclin E1 (CCNE1), which has been proposed as its 'driver'. The aim of this study was to identify the genes within the 19q12 amplicon whose expression is required for the survival of cancer cells harbouring their amplification.

ARID1A (AT-rich interactive domain containing protein 1A) loss-of-function mutations have been reported in gynecological cancers, including rarer subtypes such as clear cell carcinoma. Preclinical studies indicate that ARID1A mutant cancers display sensitivity to ATR inhibition while tumors without ARID1A mutations may be sensitive to Ataxia telangiectasia and Rad3 related (ATR) inhibitors in combination with poly-ADP ribose polymerase (PARP) inhibitors.

Cancer cells are known to display varying degrees of metastatic propensity, but the molecular basis underlying such heterogeneity remains unclear. Our aims in this study were to (i) elucidate prognostic subtypes in primary tumors based on an epithelial-to-mesenchymal-to-amoeboid transition (EMAT) continuum that captures the heterogeneity of metastatic propensity and (ii) to more comprehensively define biologically informed subtypes predictive of breast cancer metastasis and survival in lymph node-negative (LNN) patients.

One approach to identifying cancer-specific vulnerabilities and therapeutic targets is to profile genetic dependencies in cancer cell lines. Here, we describe data from a series of siRNA screens that identify the kinase genetic dependencies in 117 cancer cell lines from ten cancer types. By integrating the siRNA screen data with molecular profiling data, including exome sequencing data, we show how vulnerabilities/genetic dependencies that are associated with mutations in specific cancer driver genes can be identified. By integrating additional data sets into this analysis, including protein-protein interaction data, we also demonstrate that the genetic dependencies associated with many cancer driver genes form dense connections on functional interaction networks. We demonstrate the utility of this resource by using it to predict the drug sensitivity of genetically or histologically defined subsets of tumor cell lines, including an increased sensitivity of osteosarcoma cell lines to FGFR inhibitors and SMAD4 mutant tumor cells to mitotic inhibitors.

RAD21 is a component of the cohesin complex, which is essential for chromosome segregation and error-free DNA repair. We assessed its prognostic and predictive power in a cohort of in situ and invasive breast cancers, and its effect on chemosensitivity in vitro.

Gene fusions arising from chromosomal translocations have been implicated in cancer. However, the role of gene fusions in BRCA1-related breast cancers is not well understood. Mutations in BRCA1 are associated with an increased risk for breast cancer (up to 80% lifetime risk) and ovarian cancer (up to 50%). We sought to identify putative gene fusions in the transcriptomes of these cancers using high-throughput RNA sequencing (RNA-Seq).