Pathway-based signatures predict patient outcome, chemotherapy benefit and synthetic lethal dependencies in invasive lobular breast cancer.

BACKGROUND: Invasive Lobular Carcinoma (ILC) is a morphologically distinct breast cancer subtype that represents up to 15% of all breast cancers. Compared to Invasive Breast Carcinoma of No Special Type (IBC-NST), ILCs exhibit poorer long-term outcome and a unique pattern of metastasis. Despite these differences, the systematic discovery of robust prognostic biomarkers and therapeutically actionable molecular pathways in ILC remains limited. METHODS: Pathway-centric multivariable models using statistical machine learning were developed and tested in seven retrospective clinico-genomic cohorts (n = 996). Further external validation was performed using a new RNA-Seq clinical cohort of aggressive ILCs (n = 48). RESULTS AND CONCLUSIONS: mRNA dysregulation scores of 25 pathways were strongly prognostic in ILC (FDR-adjusted P < 0.05). Of these, three pathways including Cell-cell communication, Innate immune system and Smooth muscle contraction were also independent predictors of chemotherapy response. To aggregate these findings, a multivariable machine learning predictor called PSILC was developed and successfully validated for predicting overall and metastasis-free survival in ILC. Integration of PSILC with CRISPR-Cas9 screening data from breast cancer cell lines revealed 16 candidate therapeutic targets that were synthetic lethal with high-risk ILCs. This study provides interpretable prognostic and predictive biomarkers of ILC which could serve as the starting points for targeted drug discovery for this disease.

Implications of tumour heterogeneity on cancer evolution and therapy resistance: lessons from breast cancer.

Tumour heterogeneity is pervasive amongst many cancers and leads to disease progression, and therapy resistance. In this review, using breast cancer as an exemplar, we focus on the recent advances in understanding the interplay between tumour cells and their microenvironment using single cell sequencing and digital spatial profiling technologies. Further, we discuss the utility of lineage tracing methodologies in pre-clinical models of breast cancer, and how these are being used to unravel new therapeutic vulnerabilities and reveal biomarkers of breast cancer progression. © 2023 The Pathological Society of Great Britain and Ireland.

Definitive study shows no association between ARID1A mutation status and clinical outcome in endometriosis-related ovarian cancers‡.

The ARID1A tumour suppressor protein is a component of the SWI/SNF chromatin remodelling complex, which is mutated in approximately 20% of all human cancers. ARID1A mutational status is considered to hold prognostic significance in a range of solid malignancies, yet in endometriosis-related ovarian carcinomas there has been a lack of clarity of its prognostic role. Moreover, the relationship between ARID1A status and immune infiltrate is also poorly understood. In a recent issue of The Journal of Pathology, a large comprehensive study by Heinze, Nazeran et al addressed these areas by reviewing 1,623 endometriosis-associated ovarian carcinomas and correlating ARID1A status using standardised immunohistochemistry to infer mutation status, with comprehensive clinicopathological features, mismatch repair status and CD8+ tumour infiltrating lymphocytes. The study definitively showed that ARID1A status does not provide any independent prognostic value in endometriosis-associated ovarian carcinomas. ARID1A loss was, however, shown to be associated with mismatch repair deficiency and increased CD8+ tumour infiltrating lymphocytes in endometrioid ovarian carcinoma, which may be relevant for future studies. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Translational genomics of ovarian clear cell carcinoma.

Ovarian clear cell carcinomas (OCCC) are rare aggressive, chemo-resistant tumours comprising approximately 13% of all epithelial ovarian cancers, which have distinct clinical and molecular features, when compared to other gynaecological malignancies. At present, there are no specific licensed targeted therapies for OCCC, although a number of candidate targets have been identified. This review focuses on recent knowledge underpinning our understanding of the pathogenesis of OCCC including direct and synthetic-lethal therapeutic strategies in particular focussing on ARID1A deficiency. We also discuss current targeted clinical trials and immunotherapeutic approaches.

Real-time ex vivo perfusion of human lymph nodes invaded by cancer (REPLICANT): a feasibility study.

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.

ATARI trial: ATR inhibitor in combination with olaparib in gynecological cancers with ARID1A loss or no loss (ENGOT/GYN1/NCRI).

BACKGROUND: 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. PRIMARY OBJECTIVE: To determine whether the ATR inhibitor, ceralasertib, has clinical activity as a single agent and in combination with the PARP inhibitor, olaparib, in patients with ARID1A 'loss' and 'no loss' clear cell carcinomas and other relapsed gynecological cancers. STUDY HYPOTHESIS: ARID1A deficient clear cell carcinoma of the ovary or endometrium is sensitive to ATR inhibition, while the combination of ATR and PARP inhibition has activity in other gynecological tumors, irrespective of ARID1A status. TRIAL DESIGN: ATARI (ENGOT/GYN1/NCRI) is a multicenter, international, proof-of-concept, phase II, parallel cohort trial assessing ceralasertib activity as a single agent and in combination with olaparib in ARID1A stratified gynecological cancers. Patients with relapsed ovarian/endometrial clear cell carcinoma with ARID1A loss will receive ceralasertib monotherapy (cohort 1A). Relapsed ovarian/endometrial clear cell carcinoma patients with no ARID1A loss (cohort 2) or patients with other histological subtypes (endometrioid, carcinosarcoma, cervical) (cohort 3) will receive combination therapy (olaparib/ceralasertib). Treatment will continue until disease progression. MAJOR INCLUSION/EXCLUSION CRITERIA: Patients with histologically confirmed recurrent clear cell (ovarian, endometrial, or endometriosis related), endometrioid (ovarian, endometrial, or endometriosis related), cervical (adenocarcinomas or squamous), or carcinosarcomas (ovarian or endometrial) are eligible. Patients progressing after ≥1 prior platinum with evidence of measurable (RECIST v1.1) radiological disease progression since last systemic anticancer therapy and prior to trial entry are eligible. Previous ATR or PARP inhibitor treatment is not permissible. PRIMARY ENDPOINT: Best overall objective response rate (RECIST v1.1). SAMPLE SIZE: A minimum of 40 and a maximum of 116. ESTIMATED DATES FOR COMPLETING ACCRUAL AND PRESENTING RESULTS: Accrual is anticipated to be complete by the second quarter of 2022, with reporting of results by the fourth quarter of 2022. Overall accrual targets and reporting timelines are dependent on individual cohort progression to stage 2. TRIAL REGISTRATION NUMBER: NCT0405269.

Superior breast cancer metastasis risk stratification using an epithelial-mesenchymal-amoeboid transition gene signature.

BACKGROUND: 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. METHODS: We constructed a novel metastasis biology-based gene signature (EMAT) derived exclusively from cancer cells induced to undergo either epithelial-to-mesenchymal transition (EMT) or mesenchymal-to-amoeboid transition (MAT) to gauge their metastatic potential. Genome-wide gene expression data obtained from 913 primary tumors of lymph node-negative breast cancer (LNNBC) patients were analyzed. EMAT gene signature-based prognostic stratification of patients was performed to identify biologically relevant subtypes associated with distinct metastatic propensity. RESULTS: Delineated EMAT subtypes display a biologic range from less stem-like to more stem-like cell states and from less invasive to more invasive modes of cancer progression. Consideration of EMAT subtypes in combination with standard clinical parameters significantly improved survival prediction. EMAT subtypes outperformed prognosis accuracy of receptor or PAM50-based BC intrinsic subtypes even after adjusting for treatment variables in 3 independent, LNNBC cohorts including a treatment-naïve patient cohort. CONCLUSIONS: EMAT classification is a biologically informed method that provides prognostic information beyond that which can be provided by traditional cancer staging or PAM50 molecular subtype status and may improve metastasis risk assessment in early stage, LNNBC patients, who may otherwise be perceived to be at low metastasis risk.

Unbiased analysis of potential targets of breast cancer susceptibility loci by Capture Hi-C.

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.

Three-dimensional modelling identifies novel genetic dependencies associated with breast cancer progression in the isogenic MCF10 model.

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.

Comparative proteomic assessment of matrisome enrichment methodologies.

The matrisome is a complex and heterogeneous collection of extracellular matrix (ECM) and ECM-associated proteins that play important roles in tissue development and homeostasis. While several strategies for matrisome enrichment have been developed, it is currently unknown how the performance of these different methodologies compares in the proteomic identification of matrisome components across multiple tissue types. In the present study, we perform a comparative proteomic assessment of two widely used decellularisation protocols and two extraction methods to characterise the matrisome in four murine organs (heart, mammary gland, lung and liver). We undertook a systematic evaluation of the performance of the individual methods on protein yield, matrisome enrichment capability and the ability to isolate core matrisome and matrisome-associated components. Our data find that sodium dodecyl sulphate (SDS) decellularisation leads to the highest matrisome enrichment efficiency, while the extraction protocol that comprises chemical and trypsin digestion of the ECM fraction consistently identifies the highest number of matrisomal proteins across all types of tissue examined. Matrisome enrichment had a clear benefit over non-enriched tissue for the comprehensive identification of matrisomal components in murine liver and heart. Strikingly, we find that all four matrisome enrichment methods led to significant losses in the soluble matrisome-associated proteins across all organs. Our findings highlight the multiple factors (including tissue type, matrisome class of interest and desired enrichment purity) that influence the choice of enrichment methodology, and we anticipate that these data will serve as a useful guide for the design of future proteomic studies of the matrisome.

Microenvironmental Heterogeneity Parallels Breast Cancer Progression: A Histology-Genomic Integration Analysis.

BACKGROUND: 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. METHODS AND FINDINGS: We developed a quantitative measure of microenvironmental heterogeneity along three spatial dimensions (3-D) in solid tumors, termed the tumor ecosystem diversity index (EDI), using fully automated histology image analysis coupled with statistical measures commonly used in ecology. This measure was compared with disease-specific survival, key mutations, genome-wide copy number, and expression profiling data in a retrospective study of 510 breast cancer patients as a test set and 516 breast cancer patients as an independent validation set. In high-grade (grade 3) breast cancers, we uncovered a striking link between high microenvironmental heterogeneity measured by EDI and a poor prognosis that cannot be explained by tumor size, genomics, or any other data types. However, this association was not observed in low-grade (grade 1 and 2) breast cancers. The prognostic value of EDI was superior to known prognostic factors and was enhanced with the addition of TP53 mutation status (multivariate analysis test set, p = 9 × 10-4, hazard ratio = 1.47, 95% CI 1.17-1.84; validation set, p = 0.0011, hazard ratio = 1.78, 95% CI 1.26-2.52). Integration with genome-wide profiling data identified losses of specific genes on 4p14 and 5q13 that were enriched in grade 3 tumors with high microenvironmental diversity that also substratified patients into poor prognostic groups. Limitations of this study include the number of cell types included in the model, that EDI has prognostic value only in grade 3 tumors, and that our spatial heterogeneity measure was dependent on spatial scale and tumor size. CONCLUSIONS: To our knowledge, this is the first study to couple unbiased measures of microenvironmental heterogeneity with genomic alterations to predict breast cancer clinical outcome. We propose a clinically relevant role of microenvironmental heterogeneity for advanced breast tumors, and highlight that ecological statistics can be translated into medical advances for identifying a new type of biomarker and, furthermore, for understanding the synergistic interplay of microenvironmental heterogeneity with genomic alterations in cancer cells.

APRIN is a cell cycle specific BRCA2-interacting protein required for genome integrity and a predictor of outcome after chemotherapy in breast cancer.

Mutations in BRCA2 confer an increased risk of cancer development, at least in part because the BRCA2 protein is required for the maintenance of genomic integrity. Here, we use proteomic profiling to identify APRIN (PDS5B), a cohesion-associated protein, as a BRCA2-associated protein. After exposure of cells to hydroxyurea or aphidicolin, APRIN and other cohesin components associate with BRCA2 in early S-phase. We demonstrate that APRIN expression is required for the normal response to DNA-damaging agents, the nuclear localisation of RAD51 and BRCA2 and efficient homologous recombination. The clinical significance of these findings is indicated by the observation that the BRCA2/APRIN interaction is compromised by BRCA2 missense variants of previously unknown significance and that APRIN expression levels are associated with histological grade in breast cancer and the outcome of breast cancer patients treated with DNA-damaging chemotherapy.

Breast cancer heterogeneity: parallel evolution or conscious uncoupling?

Breast cancer is known to display considerable inter- and intra-tumour genetic heterogeneity. It is now widely accepted that no two breast cancers harbour the same complement of genomic alterations, and that both primary and metastatic breast cancers are composed of multiple genetically diverse subclones that evolve under different selective pressures. Recent work published in the Journal of Pathology by Desmedt and colleagues questions the evolutionary dynamics of multi-focal breast cancer with similar pathological features by studying the mutational repertoire of different lesions. Whilst the majority of the lesions showed some common driver alterations, one-third lacked any common mutations, suggesting very early clonal divergence. These and other recent observations underscore the need for a fundamental understanding of the rules governing breast cancer evolution, and highlight the need for in-depth assessment of driver alterations for appropriate patient management and selective treatment.

SF3B1 mutations constitute a novel therapeutic target in breast cancer.

Mutations in genes encoding proteins involved in RNA splicing have been found to occur at relatively high frequencies in several tumour types including myelodysplastic syndromes, chronic lymphocytic leukaemia, uveal melanoma, and pancreatic cancer, and at lower frequencies in breast cancer. To investigate whether dysfunction in RNA splicing is implicated in the pathogenesis of breast cancer, we performed a re-analysis of published exome and whole genome sequencing data. This analysis revealed that mutations in spliceosomal component genes occurred in 5.6% of unselected breast cancers, including hotspot mutations in the SF3B1 gene, which were found in 1.8% of unselected breast cancers. SF3B1 mutations were significantly associated with ER-positive disease, AKT1 mutations, and distinct copy number alterations. Additional profiling of hotspot mutations in a panel of special histological subtypes of breast cancer showed that 16% and 6% of papillary and mucinous carcinomas of the breast harboured the SF3B1 K700E mutation. RNA sequencing identified differentially spliced events expressed in tumours with SF3B1 mutations including the protein coding genes TMEM14C, RPL31, DYNL11, UQCC, and ABCC5, and the long non-coding RNA CRNDE. Moreover, SF3B1 mutant cell lines were found to be sensitive to the SF3b complex inhibitor spliceostatin A and treatment resulted in perturbation of the splicing signature. Albeit rare, SF3B1 mutations result in alternative splicing events, and may constitute drivers and a novel therapeutic target in a subset of breast cancers.

Genomic landscape of adenoid cystic carcinoma of the breast.

Adenoid cystic carcinoma (AdCC) is a rare type of triple-negative breast cancer (TNBC) characterized by the presence of the MYB-NFIB fusion gene. The molecular underpinning of breast AdCCs other than the MYB-NFIB fusion gene remains largely unexplored. Here we sought to define the repertoire of somatic genetic alterations of breast AdCCs. We performed whole-exome sequencing, followed by orthogonal validation, of 12 breast AdCCs to determine the landscape of somatic mutations and gene copy number alterations. Fluorescence in situ hybridization and reverse-transcription PCR were used to define the presence of MYB gene rearrangements and MYB-NFIB chimeric transcripts. Unlike common forms of TNBC, we found that AdCCs have a low mutation rate (0.27 non-silent mutations/Mb), lack mutations in TP53 and PIK3CA and display a heterogeneous constellation of known cancer genes affected by somatic mutations, including MYB, BRAF, FBXW7, SMARCA5, SF3B1 and FGFR2. MYB and TLN2 were affected by somatic mutations in two cases each. Akin to salivary gland AdCCs, breast AdCCs were found to harbour mutations targeting chromatin remodelling, cell adhesion, RNA biology, ubiquitination and canonical signalling pathway genes. We observed that, although breast AdCCs had rather simple genomes, they likely display intra-tumour genetic heterogeneity at diagnosis. Taken together, these findings demonstrate that the mutational burden and mutational repertoire of breast AdCCs are more similar to those of salivary gland AdCCs than to those of other types of TNBCs, emphasizing the importance of histological subtyping of TNBCs. Furthermore, our data provide direct evidence that AdCCs harbour a distinctive mutational landscape and genomic structure, irrespective of the disease site of origin.

An ecological measure of immune-cancer colocalization as a prognostic factor for breast cancer.

INTRODUCTION: Abundance of immune cells has been shown to have prognostic and predictive significance in many tumor types. Beyond abundance, the spatial organization of immune cells in relation to cancer cells may also have significant functional and clinical implications. However there is a lack of systematic methods to quantify spatial associations between immune and cancer cells. METHODS: We applied ecological measures of species interactions to digital pathology images for investigating the spatial associations of immune and cancer cells in breast cancer. We used the Morisita-Horn similarity index, an ecological measure of community structure and predator-prey interactions, to quantify the extent to which cancer cells and immune cells colocalize in whole-tumor histology sections. We related this index to disease-specific survival of 486 women with breast cancer and validated our findings in a set of 516 patients from different hospitals. RESULTS: Colocalization of immune cells with cancer cells was significantly associated with a disease-specific survival benefit for all breast cancers combined. In HER2-positive subtypes, the prognostic value of immune-cancer cell colocalization was highly significant and exceeded those of known clinical variables. Furthermore, colocalization was a significant predictive factor for long-term outcome following chemotherapy and radiotherapy in HER2 and Luminal A subtypes, independent of and stronger than all known clinical variables. CONCLUSIONS: Our study demonstrates how ecological methods applied to the tumor microenvironment using routine histology can provide reproducible, quantitative biomarkers for identifying high-risk breast cancer patients. We found that the clinical value of immune-cancer interaction patterns is highly subtype-specific but substantial and independent to known clinicopathologic variables that mostly focused on cancer itself. Our approach can be developed into computer-assisted prediction based on histology samples that are already routinely collected.

Metaplastic breast carcinomas display genomic and transcriptomic heterogeneity [corrected]. .

Metaplastic breast carcinoma is a rare and aggressive histologic type of breast cancer, preferentially displaying a triple-negative phenotype. We sought to define the transcriptomic heterogeneity of metaplastic breast cancers on the basis of current gene expression microarray-based classifiers, and to determine whether these tumors display gene copy number profiles consistent with those of BRCA1-associated breast cancers. Twenty-eight consecutive triple-negative metaplastic breast carcinomas were reviewed, and the metaplastic component present in each frozen specimen was defined (ie, spindle cell, squamous, chondroid metaplasia). RNA and DNA extracted from frozen sections with tumor cell content >60% were subjected to gene expression (Illumina HumanHT-12 v4) and copy number profiling (Affymetrix SNP 6.0), respectively. Using the best practice PAM50/claudin-low microarray-based classifier, all metaplastic breast carcinomas with spindle cell metaplasia were of claudin-low subtype, whereas those with squamous or chondroid metaplasia were preferentially of basal-like subtype. Triple-negative breast cancer subtyping using a dedicated website (http://cbc.mc.vanderbilt.edu/tnbc/) revealed that all metaplastic breast carcinomas with chondroid metaplasia were of mesenchymal-like subtype, spindle cell carcinomas preferentially of unstable or mesenchymal stem-like subtype, and those with squamous metaplasia were of multiple subtypes. None of the cases was classified as immunomodulatory or luminal androgen receptor subtype. Integrative clustering, combining gene expression and gene copy number data, revealed that metaplastic breast carcinomas with spindle cell and chondroid metaplasia were preferentially classified as of integrative clusters 4 and 9, respectively, whereas those with squamous metaplasia were classified into six different clusters. Eight of the 26 metaplastic breast cancers subjected to SNP6 analysis were classified as BRCA1-like. The diversity of histologic features of metaplastic breast carcinomas is reflected at the transcriptomic level, and an association between molecular subtypes and histology was observed. BRCA1-like genomic profiles were found only in a subset (31%) of metaplastic breast cancers, and were not associated with a specific molecular or histologic subtype.

Characterization of the genomic features and expressed fusion genes in micropapillary carcinomas of the breast.

Micropapillary carcinoma (MPC) is a rare histological special type of breast cancer, characterized by an aggressive clinical behaviour and a pattern of copy number aberrations (CNAs) distinct from that of grade- and oestrogen receptor (ER)-matched invasive carcinomas of no special type (IC-NSTs). The aims of this study were to determine whether MPCs are underpinned by a recurrent fusion gene(s) or mutations in 273 genes recurrently mutated in breast cancer. Sixteen MPCs were subjected to microarray-based comparative genomic hybridization (aCGH) analysis and Sequenom OncoCarta mutation analysis. Eight and five MPCs were subjected to targeted capture and RNA sequencing, respectively. aCGH analysis confirmed our previous observations about the repertoire of CNAs of MPCs. Sequencing analysis revealed a spectrum of mutations similar to those of luminal B IC-NSTs, and recurrent mutations affecting mitogen-activated protein kinase family genes and NBPF10. RNA-sequencing analysis identified 17 high-confidence fusion genes, eight of which were validated and two of which were in-frame. No recurrent fusions were identified in an independent series of MPCs and IC-NSTs. Forced expression of in-frame fusion genes (SLC2A1-FAF1 and BCAS4-AURKA) resulted in increased viability of breast cancer cells. In addition, genomic disruption of CDK12 caused by out-of-frame rearrangements was found in one MPC and in 13% of HER2-positive breast cancers, identified through a re-analysis of publicly available massively parallel sequencing data. In vitro analyses revealed that CDK12 gene disruption results in sensitivity to PARP inhibition, and forced expression of wild-type CDK12 in a CDK12-null cell line model resulted in relative resistance to PARP inhibition. Our findings demonstrate that MPCs are neither defined by highly recurrent mutations in the 273 genes tested, nor underpinned by a recurrent fusion gene. Although seemingly private genetic events, some of the fusion transcripts found in MPCs may play a role in maintenance of a malignant phenotype and potentially offer therapeutic opportunities.

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.

Resistance to therapy caused by intragenic deletion in BRCA2.

Cells with loss of BRCA2 function are defective in homologous recombination (HR) and are highly sensitive to inhibitors of poly(ADP-ribose) polymerase (PARP), which provides the basis for a new therapeutic approach. Here we show that resistance to PARP inhibition can be acquired by deletion of a mutation in BRCA2. We derived PARP-inhibitor-resistant (PIR) clones from the human CAPAN1 pancreatic cancer cell line, which carries the protein-truncating c.6174delT frameshift mutation. PIR clones could form DNA-damage-induced RAD51 nuclear foci and were able to limit genotoxin-induced genomic instability, both hallmarks of a competent HR pathway. New BRCA2 isoforms were expressed in the resistant lines as a result of intragenic deletion of the c.6174delT mutation and restoration of the open reading frame (ORF). Reconstitution of BRCA2-deficient cells with these revertant BRCA2 alleles rescued PARP inhibitor sensitivity and HR deficiency. Most of the deletions in BRCA2 were associated with small tracts of homology, and possibly arose from error-prone repair caused by BRCA2 deficiency. Similar ORF-restoring mutations were present in carboplatin-resistant ovarian tumours from c.6174delT mutation carriers. These observations have implications for understanding drug resistance in BRCA mutation carriers as well as in defining functionally important domains within BRCA2.