FixNCut: single-cell genomics through reversible tissue fixation and dissociation

The use of single-cell technologies for clinical applications requires disconnecting sampling from downstream processing steps. Early sample preservation can further increase robustness and reproducibility by avoiding artifacts introduced during specimen handling. We present FixNCut, a methodology for the reversible fixation of tissue followed by dissociation that overcomes current limitations. We applied FixNCut to human and mouse tissues to demonstrate the preservation of RNA integrity, sequencing library complexity, and cellular composition, while diminishing stress-related artifacts. Besides single-cell RNA sequencing, FixNCut is compatible with multiple single-cell and spatial technologies, making it a versatile tool for robust and flexible study designs.

HAMSAB diet ameliorates dysfunctional signaling in pancreatic islets in autoimmune diabetes

An altered gut microbiota is associated with type 1 diabetes (T1D), affecting the production of short-chain fatty acids (SCFA) and glucose homeostasis. We previously demonstrated that enhancing serum acetate and butyrate using a dietary supplement (HAMSAB) improved glycemia in non-obese diabetic (NOD) mice and patients with established T1D. The effects of SCFA on immune-infiltrated islet cells remain to be clarified. Here, we performed single-cell RNA sequencing on islet cells from NOD mice fed an HAMSAB or control diet. HAMSAB induced a regulatory gene expression profile in pancreas-infiltrated immune cells. Moreover, HAMSAB maintained the expression of β-cell functional genes and decreased cellular stress. HAMSAB-fed mice showed preserved pancreatic endocrine cell identity, evaluated by decreased numbers of poly-hormonal cells. Finally, SCFA increased insulin levels in human β-like cells and improved transplantation outcome in NOD/SCID mice. Our findings support the use of metabolite-based diet as attractive approach to improve glucose control in T1D.

Ligand-dependent hedgehog signaling maintains an undifferentiated, malignant osteosarcoma phenotype

TP53 and RB1 loss-of-function mutations are common in osteosarcoma. During development, combined loss of TP53 and RB1 function leads to downregulation of autophagy and the aberrant formation of primary cilia, cellular organelles essential for the transmission of canonical Hedgehog (Hh) signaling. Excess cilia formation then leads to hypersensitivity to Hedgehog (Hh) ligand signaling. In mouse and human models, we now show that osteosarcomas with mutations in TP53 and RB1 exhibit enhanced ligand-dependent Hh pathway activation through Smoothened (SMO), a transmembrane signaling molecule required for activation of the canonical Hh pathway. This dependence is mediated by hypersensitivity to Hh ligand and is accompanied by impaired autophagy and increased primary cilia formation and expression of Hh ligand in vivo. Using a conditional genetic mouse model of Trp53 and Rb1 inactivation in osteoblast progenitors, we further show that deletion of Smo converts the highly malignant osteosarcoma phenotype to benign, well differentiated bone tumors. Conversely, conditional overexpression of SHH ligand, or a gain-of-function SMO mutant in committed osteoblast progenitors during development blocks terminal bone differentiation. Finally, we demonstrate that the SMO antagonist sonidegib (LDE225) induces growth arrest and terminal differentiation in vivo in osteosarcomas that express primary cilia and Hh ligand combined with mutations in TP53. These results provide a mechanistic framework for aberrant Hh signaling in osteosarcoma based on defining mutations in the tumor suppressor, TP53.

Fancm has dual roles in the limiting of meiotic crossovers and germ cell maintenance in mammals

Meiotic crossovers are required for accurate chromosome segregation and producing new allelic combinations. Meiotic crossover numbers are tightly regulated within a narrow range, despite an excess of initiating DNA double-strand breaks. Here, we reveal the tumor suppressor FANCM as a meiotic anti-crossover factor in mammals. We use unique large-scale crossover analyses with both single-gamete sequencing and pedigree-based bulk-sequencing datasets to identify a genome-wide increase in crossover frequencies in Fancm-deficient mice. Gametogenesis is heavily perturbed in Fancm loss-of-function mice, which is consistent with the reproductive defects reported in humans with biallelic FANCM mutations. A portion of the gametogenesis defects can be attributed to the cGAS-STING pathway after birth. Despite the gametogenesis phenotypes in Fancm mutants, both sexes are capable of producing offspring. We propose that the anti-crossover function and role in gametogenesis of Fancm are separable and will inform diagnostic pathways for human genomic instability disorders.

Systematic benchmarking of single-cell ATAC-sequencing protocols

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Inhibition of mutant IDH1 promotes cycling of acute myeloid leukemia stem cells

Approximately 20% of acute myeloid leukemia (AML) patients carry mutations in IDH1 or IDH2 that result in over-production of the oncometabolite D-2-hydroxyglutarate (2-HG). Small molecule inhibitors that block 2-HG synthesis can induce complete morphological remission; however, almost all patients eventually acquire drug resistance and relapse. Using a multi-allelic mouse model of IDH1-mutant AML, we demonstrate that the clinical IDH1 inhibitor AG-120 (ivosidenib) exerts cell-type-dependent effects on leukemic cells, promoting delayed disease regression. Although single-agent AG-120 treatment does not fully eradicate the disease, it increases cycling of rare leukemia stem cells and triggers transcriptional upregulation of the pyrimidine salvage pathway. Accordingly, AG-120 sensitizes IDH1-mutant AML to azacitidine, with the combination of AG-120 and azacitidine showing vastly improved efficacy in vivo. Our data highlight the impact of non-genetic heterogeneity on treatment response and provide a mechanistic rationale for the observed combinatorial effect of AG-120 and azacitidine in patients.

Epigenetic Activation of Plasmacytoid DCs Drives IFNAR-Dependent Therapeutic Differentiation of AML

Pharmacologic inhibition of epigenetic enzymes can have therapeutic benefit against hematologic malignancies. In addition to affecting tumor cell growth and proliferation, these epigenetic agents may induce antitumor immunity. Here, we discovered a novel immunoregulatory mechanism through inhibition of histone deacetylases (HDAC). In models of acute myeloid leukemia (AML), leukemia cell differentiation and therapeutic benefit mediated by the HDAC inhibitor (HDACi) panobinostat required activation of the type I interferon (IFN) pathway. Plasmacytoid dendritic cells (pDC) produced type I IFN after panobinostat treatment, through transcriptional activation of IFN genes concomitant with increased H3K27 acetylation at these loci. Depletion of pDCs abrogated panobinostat-mediated induction of type I IFN signaling in leukemia cells and impaired therapeutic efficacy, whereas combined treatment with panobinostat and IFNα improved outcomes in preclinical models. These discoveries offer a new therapeutic approach for AML and demonstrate that epigenetic rewiring of pDCs enhances antitumor immunity, opening the possibility of exploiting this approach for immunotherapies.

SIGNIFICANCE

We demonstrate that HDACis induce terminal differentiation of AML through epigenetic remodeling of pDCs, resulting in production of type I IFN that is important for the therapeutic effects of HDACis. The study demonstrates the important functional interplay between the immune system and leukemias in response to HDAC inhibition. This article is highlighted in the In This Issue feature, p. 1397.

Same-Cell Co-Occurrence of RAS Hotspot and BRAF V600E Mutations in Treatment-Naive Colorectal Cancer

PURPOSE

Mitogen-activated protein kinase pathway-activating mutations occur in the majority of colorectal cancer (CRC) cases and show mutual exclusivity. We identified 47 epidermal growth factor receptor/BRAF inhibitor-naive CRC patients with dual RAS hotspot/BRAF V600E mutations (CRC-DD) from a cohort of 4,561 CRC patients with clinical next-generation sequencing results. We aimed to define the molecular phenotypes of the CRC-DD and to test if the dual RAS hotspot/BRAF V600E mutations coexist within the same cell.

MATERIALS AND METHODS

We developed a single-cell genotyping method with a mutation detection rate of 96.3% and a genotype prediction accuracy of 92.1%. Mutations in the CRC-DD cohort were analyzed for clonality, allelic imbalance, copy number, and overall survival.

RESULTS

Application of single-cell genotyping to four CRC-DD revealed the co-occurrence of both mutations in the following percentages of cells per case: NRAS G13D/KRAS G12C, 95%; KRAS G12D/NRAS G12V, 48%; BRAF V600E/KRAS G12D, 44%; and KRAS G12D/NRAS G13V, 14%, respectively. Allelic imbalance favoring the oncogenic allele was less frequent in CRC-DD (24 of 76, 31.5%, somatic mutations) compared with a curated cohort of CRC with a single-driver mutation (CRC-SD; 119 of 232 mutations, 51.3%; P = .013). Microsatellite instability-high status was enriched in CRC-DD compared with CRC-SD (23% v 11.4%, P = .028). Of the seven CRC-DD cases with multiregional sequencing, five retained both driver mutations throughout all sequenced tumor sites. Both CRC-DD cases with discordant multiregional sequencing were microsatellite instability-high.

CONCLUSION

Our findings indicate that dual-driver mutations occur in a rare subset of CRC, often within the same tumor cells and across multiple tumor sites. Their presence and a lower rate of allelic imbalance may be related to dose-dependent signaling within the mitogen-activated protein kinase pathway.

Non-genetic determinants of malignant clonal fitness at single-cell resolution

All cancers emerge after a period of clonal selection and subsequent clonal expansion. Although the evolutionary principles imparted by genetic intratumour heterogeneity are becoming increasingly clear¹, little is known about the non-genetic mechanisms that contribute to intratumour heterogeneity and malignant clonal fitness². Here, using single-cell profiling and lineage tracing (SPLINTR)-an expressed barcoding strategy-we trace isogenic clones in three clinically relevant mouse models of acute myeloid leukaemia. We find that malignant clonal dominance is a cell-intrinsic and heritable property that is facilitated by the repression of antigen presentation and increased expression of the secretory leukocyte peptidase inhibitor gene (Slpi), which we genetically validate as a regulator of acute myeloid leukaemia. Increased transcriptional heterogeneity is a feature that enables clonal fitness in diverse tissues and immune microenvironments and in the context of clonal competition between genetically distinct clones. Similar to haematopoietic stem cells³, leukaemia stem cells (LSCs) display heritable clone-intrinsic properties of high, and low clonal output that contribute to the overall tumour mass. We demonstrate that LSC clonal output dictates sensitivity to chemotherapy and, although high- and low-output clones adapt differently to therapeutic pressure, they coordinately emerge from minimal residual disease with increased expression of the LSC program. Together, these data provide fundamental insights into the non-genetic transcriptional processes that underpin malignant clonal fitness and may inform future therapeutic strategies.

CDK4/6 inhibition promotes anti-tumor immunity through the induction of T cell memory

Pharmacological inhibitors of cyclin dependent kinases 4 and 6 (CDK4/6) are an approved treatment for hormone receptor-positive breast cancer and are currently under evaluation across hundreds of clinical trials for other cancer types. The clinical success of these inhibitors is largely attributed to well-defined tumor-intrinsic cytostatic mechanisms, while their emerging role as immunomodulatory agents is less understood. Using integrated epigenomic, transcriptomic and proteomic analyses, we demonstrated a novel action of CDK4/6 inhibitors in promoting the phenotypic and functional acquisition of immunological T cell memory. Short-term priming with a CDK4/6 inhibitor promoted long-term endogenous anti-tumor T cell immunity in mice, enhanced the persistence and therapeutic efficacy of chimeric antigen receptor (CAR)-T cells, and induced an RB-dependent T cell phenotype supportive of favorable responses to immune checkpoint blockade in melanoma patients. Together, these mechanistic insights significantly broaden the prospective utility of CDK4/6 inhibitors as clinical tools to boost anti-tumor T cell immunity.

γδ T Cells in Merkel Cell Carcinomas Have a Proinflammatory Profile Prognostic of Patient Survival

Merkel cell carcinomas (MCC) are immunogenic skin cancers associated with viral infection or UV mutagenesis. To study T-cell infiltrates in MCC, we analyzed 58 MCC lesions from 39 patients using multiplex-IHC/immunofluorescence (m-IHC/IF). CD4⁺ or CD8⁺ T cells comprised the majority of infiltrating T lymphocytes in most tumors. However, almost half of the tumors harbored prominent CD4/CD8 double-negative (DN) T-cell infiltrates (>20% DN T cells), and in 12% of cases, DN T cells represented the majority of T cells. Flow cytometric analysis of single-cell suspensions from fresh tumors identified DN T cells as predominantly Vδ2^- γδ T cells. In the context of γδ T-cell inflammation, these cells expressed PD-1 and LAG3, which is consistent with a suppressed or exhausted phenotype, and CD103, which indicates tissue residency. Furthermore, single-cell RNA sequencing (scRNA-seq) identified a transcriptional profile of γδ T cells suggestive of proinflammatory potential. T-cell receptor (TCR) analysis confirmed clonal expansion of Vδ1 and Vδ3 clonotypes, and functional studies using cloned γδ TCRs demonstrated restriction of these for CD1c and MR1 antigen-presenting molecules. On the basis of a 13-gene γδ T-cell signature derived from scRNA-seq analysis, gene-set enrichment on bulk RNA-seq data showed a positive correlation between enrichment scores and DN T-cell infiltrates. An improved disease-specific survival was evident for patients with high enrichment scores, and complete responses to anti-PD-1/PD-L1 treatment were observed in three of four cases with high enrichment scores. Thus, γδ T-cell infiltration may serve as a prognostic biomarker and should be explored for therapeutic interventions.See related Spotlight on p. 600.

Combined BRAF, MEK, and CDK4/6 Inhibition Depletes Intratumoral Immune-Potentiating Myeloid Populations in Melanoma

Combined inhibition of BRAF, MEK, and CDK4/6 is currently under evaluation in clinical trials for patients with melanoma harboring a BRAF^V600 mutation. While this triple therapy has potent tumor-intrinsic effects, the impact of this combination on antitumor immunity remains unexplored. Here, using a syngeneic Braf^V600ECdkn2a^-/-Pten^-/- melanoma model, we demonstrated that triple therapy promoted durable tumor control through tumor-intrinsic mechanisms and promoted immunogenic cell death and T-cell infiltration. Despite this, tumors treated with triple therapy were unresponsive to immune checkpoint blockade (ICB). Flow cytometric and single-cell RNA sequencing analyses of tumor-infiltrating immune populations revealed that triple therapy markedly depleted proinflammatory macrophages and cross-priming CD103⁺ dendritic cells, the absence of which correlated with poor overall survival and clinical responses to ICB in patients with melanoma. Indeed, immune populations isolated from tumors of mice treated with triple therapy failed to stimulate T-cell responses ex vivo While combined BRAF, MEK, and CDK4/6 inhibition demonstrates favorable tumor-intrinsic activity, these data suggest that collateral effects on tumor-infiltrating myeloid populations may impact antitumor immunity. These findings have important implications for the design of combination strategies and clinical trials that incorporate BRAF, MEK, and CDK4/6 inhibition with immunotherapy for the treatment of patients with melanoma.

Prevalence and potential biological role of TERT amplifications in ALK translocated adenocarcinoma of the lung

AIMS

The advent of specific ALK-targeting drugs has radically changed the outcome of patients with ALK translocated non-small-cell lung cancer (NSCLC). However, emerging resistance to treatment with ALK inhibitors in these patients remains a major concern. In previous studies, we analysed two ALK+ patient cohorts (TP53 wild-type/TP53 mutated) in terms of copy number alterations. All patients belonging to the TP53 wild-type group had mainly genetically stable genomes, with one exception showing chromosomal instability and amplifications of several gene loci, including TERT. Here, we aimed to determine the prevalence of TERT amplifications in these ALK+ lung cancer patients by analysing an independent cohort of 109 ALK translocated cases. We further analysed the copy numbers of numerous cancer-relevant genes and other genetic aberrations.

METHODS AND RESULTS

The prevalence of TERT amplifications was determined by means of FISH analyses. Copy numbers of 87 cancer-relevant genes were determined by NanoString nCounter^® technology, FoundationOne^® and lung-specific NGS panels in some of these TERT-amplified samples, and clinical data on patients with TERT-amplified tumours were collected. Our data revealed that five (4.6%) of all 109 analysed ALK+ patients harboured amplification of TERT and that these patients had genetically unstable genomes.

CONCLUSIONS

Our preliminary study shows that ALK+ adenocarcinomas should be evaluated in the context of their genomic background in order to more clearly understand and predict patients' individual course of disease.

From whole-mount to single-cell spatial assessment of gene expression in 3D

Unravelling spatio-temporal patterns of gene expression is crucial to understanding core biological principles from embryogenesis to disease. Here we review emerging technologies, providing automated, high-throughput, spatially resolved quantitative gene expression data. Novel techniques expand on current benchmark protocols, expediting their incorporation into ongoing research. These approaches digitally reconstruct patterns of embryonic expression in three dimensions, and have successfully identified novel domains of expression, cell types, and tissue features. Such technologies pave the way for unbiased and exhaustive recapitulation of gene expression levels in spatial and quantitative terms, promoting understanding of the molecular origin of developmental defects, and improving medical diagnostics.

Renal epithelial cells retain primary cilia during human acute renal allograft rejection injury

Objectives

Primary cilia are sensory organelles which co-ordinate several developmental/repair pathways including hedgehog signalling. Studies of human renal allografts suffering acute tubular necrosis have shown that length of primary cilia borne by epithelial cells doubles throughout the nephron and collecting duct, and then normalises as renal function returns. Conversely the loss of primary cilia has been reported in chronic allograft rejection and linked to defective hedgehog signalling. We investigated the fate of primary cilia in renal allografts suffering acute rejection.

Results

Here we observed that in renal allografts undergoing acute rejection, primary cilia were retained, with their length increasing 1 week after transplantation and remaining elevated. We used a mouse model of acute renal injury to demonstrate that elongated renal primary cilia in the injured renal tubule show evidence of smoothened accumulation, a biomarker for activation of hedgehog signalling. We conclude that primary cilium-mediated activation of hedgehog signalling is still possible during the acute phase of renal allograft rejection.

Inhibition of activin signaling in lung adenocarcinoma increases the therapeutic index of platinum chemotherapy

Resistance to platinum chemotherapy is a long-standing problem in the management of lung adenocarcinoma. Using a whole-genome synthetic lethal RNA interference screen, we identified activin signaling as a critical mediator of innate platinum resistance. The transforming growth factor-β (TGFβ) superfamily ligands activin A and growth differentiation factor 11 (GDF11) mediated resistance via their cognate receptors through TGFβ-activated kinase 1 (TAK1), rather than through the SMAD family of transcription factors. Inhibition of activin receptor signaling or blockade of activin A and GDF11 by the endogenous protein follistatin overcame this resistance. Consistent with the role of activin signaling in acute renal injury, both therapeutic interventions attenuated acute cisplatin-induced nephrotoxicity, its major dose-limiting side effect. This cancer-specific enhancement of platinum-induced cell death has the potential to dramatically improve the safety and efficacy of chemotherapy in lung cancer patients.

Targeting enhancer switching overcomes non-genetic drug resistance in acute myeloid leukaemia

Non-genetic drug resistance is increasingly recognised in various cancers. Molecular insights into this process are lacking and it is unknown whether stable non-genetic resistance can be overcome. Using single cell RNA-sequencing of paired drug naïve and resistant AML patient samples and cellular barcoding in a unique mouse model of non-genetic resistance, here we demonstrate that transcriptional plasticity drives stable epigenetic resistance. With a CRISPR-Cas9 screen we identify regulators of enhancer function as important modulators of the resistant cell state. We show that inhibition of Lsd1 (Kdm1a) is able to overcome stable epigenetic resistance by facilitating the binding of the pioneer factor, Pu.1 and cofactor, Irf8, to nucleate new enhancers that regulate the expression of key survival genes. This enhancer switching results in the re-distribution of transcriptional co-activators, including Brd4, and provides the opportunity to disable their activity and overcome epigenetic resistance. Together these findings highlight key principles to help counteract non-genetic drug resistance.

Genetic instability and recurrent MYC amplification in ALK-translocated NSCLC: a central role of TP53 mutations

The anaplastic lymphoma kinase (ALK) rearrangement defines a distinct molecular subtype of non-small cell lung cancer (NSCLC). Despite the excellent initial efficacy of ALK inhibitors in patients with ALK+ lung cancer, resistance occurs almost inevitably. To date, there is no reliable biomarker allowing the identification of patients at higher risk of relapse. Here, we analysed a subset of 53 ALK+ tumours with and without TP53 mutation and ALK+ NSCLC cell lines by NanoString nCounter technology. We found that the co-occurrence of early TP53 mutations in ALK+ NSCLC can lead to chromosomal instability: 24% of TP53-mutated patients showed amplifications of known cancer genes such as MYC (14%), CCND1 (10%), TERT (5%), BIRC2 (5%), ORAOV1 (5%), and YAP1 (5%). MYC-overexpressing ALK+ TP53-mutated cells had a proliferative advantage compared to wild-type cells. ChIP-Seq data revealed MYC-binding sites within the promoter region of EML4, and MYC overexpression in ALK+ TP53-mutated cells resulted in an upregulation of EML4-ALK, indicating a potential MYC-dependent resistance mechanism in patients with increased MYC copy number. Our study reveals that ALK+ NSCLC represents a more heterogeneous subgroup of tumours than initially thought, and that TP53 mutations in that particular cancer type define a subset of tumours that harbour chromosomal instability, leading to the co-occurrence of pathogenic aberrations. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Abstract P2-05-03: Novel driver genetic alterations in MYB-NFIB-negative breast adenoid cystic carcinomas

Introduction: Breast adenoid cystic carcinoma (AdCC) is a rare type of triple-negative breast cancer associated with an indolent clinical behavior. AdCCs provide a clear example of genotypic-phenotypic correlation with the majority harboring the MYB-NFIB fusion gene. In this study, we sought to identify alternative driver genetic alterations in breast AdCCs lacking the MYB-NFIB fusion gene.

Methods: Nucleic acids obtained from four breast AdCCs lacking the MYB-NFIB fusion gene as defined by reverse transcription (RT)-PCR and/or fluorescence in situ hybridization (FISH) were subjected to RNA-sequencing (n=3), whole-genome (n=2) and/or targeted (n=1) massively parallel sequencing. Sequencing data were analyzed using state-of-the-art bioinformatics algorithms, and potential alternative driver genetic alterations were validated using orthogonal sequencing and molecular pathology methods.

Results: RNA-sequencing revealed the presence of MYBL1-ACTN1 or MYBL1-NFIB fusion genes in two breast AdCCs, which were validated by whole-genome sequencing and/or MYBL1 FISH analysis. Both MYBL1 fusion gene-positive cases were found to overexpress MYBL1 as defined by quantitative RT-PCR analysis. In the third MYB-NFIB-negative breast AdCC studied, a high-level MYB gene amplification coupled with overexpression of MYB at the mRNA and protein levels was identified. In the fourth breast AdCC, which expressed high levels of MYB, whole-genome and RNA-sequencing revealed no definite alternative driver alteration, however, a MYBL2 intronic mutation was found in this case, which was associated with high levels of MYBL2 mRNA expression. In this case, single sample gene set enrichment analysis revealed activation of pathways similar to those activated in AdCCs harboring the MYB-NFIB or MYBL1 fusions genes.

Conclusion: We demonstrate that in breast AdCCs lacking the MYB-NFIB fusion gene MYBL1 rearrangements and MYB amplification are likely alternative driver genetic events. Given that activation of MYB/MYBL1 and their downstream targets can be driven by the MYB-NFIB fusion gene, MYBL1 rearrangements, MYB amplification or other yet to be validated mechanisms (e.g. MYBL2 non-coding mutations), our findings further suggest that breast AdCCs constitute a convergent phenotype.

Citation Format: Kim J, Geyer FC, Martelotto LG, Ng CKY, Lim RS, Selenica P, Li A, Pareja F, Fusco N, Edelweiss M, Mariani O, Badve S, Vincent-Salomon A, Norton L, Reis-Filho JS, Weigelt B. Novel driver genetic alterations in MYB-NFIB-negative breast adenoid cystic carcinomas [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P2-05-03.

Abstract PD4-13: Estrogen receptor-negative breast adenomyoepitheliomas are driven by co-occurring HRAS hotspot and PI3K pathway gene mutations: A genetic and functional analysis

Introduction:Adenomyoepithelioma (AME) of the breast is a rare biphasic tumor, characterized by epithelial and myoepithelial differentiation. Although AMEs have an indolent clinical course, a subset may progress to carcinoma and metastasize. We sought to define the mutational landscape of AMEs and investigate the functional impact of recurrent pathogenic mutations identified in these tumors.

Methods: Thirty-one AMEs were subjected to whole-exome sequencing (WES, n=8) or massively parallel sequencing targeting all coding regions of 410 key cancer genes and intronic and regulatory regions of selected genes (n=23). Somatic genetic alterations were defined using state-of-the-art bioinformatics algorithms. In an additional set of 12 AMEs, Sanger sequencing analysis of HRAS, PIK3CA and AKT1 was performed. Non-tumorigenic estrogen receptor (ER)-negative mammary epithelial cells (i.e. MCF10A, MCF10A with a PIK3CA H1047R mutation knock-in and MCF12A) were utilized for functional studies using both conventional monolayer and three-dimensional (3D) culture assays.

Results: 27 (63%) and 16 (37%) AMEs were ER-positive and ER-negative, respectively. ER-negativity was significantly associated with histologic features predictive of a more aggressive behavior, with a higher number of mutations and copy number alterations, and with a distinct mutational profile as compared to ER-positive AMEs. Of the 27 ER-positive AMEs, 12 cases (44%) harbored PIK3CA hotspot mutations, and 5 PIK3CA wild-type cases displayed E17K AKT1 hotspot mutations. By contrast, of the 16 ER-negative AMEs, 9 (56%), 9 (56%) and 3 (19%) harbored HRAS, PIK3CA (mostly E545K and H1047R hotspots) and PIK3R1 mutations, respectively. Strikingly, all HRAS mutations were restricted to ER-negative AMEs, affected the hotspot codon Q61 (Q61R/K), and all but one co-occurred with PIK3CA or PIK3R1 mutations. In addition, HRAS Q61 hotspot mutations were significantly associated with necrosis (p=0.01) and high mitotic rates (p=0.03). CDKN2A homozygous deletions were also detected only in ER-negative AMEs (19%) and found to be significantly associated with progression to carcinoma (p=0.001). Forced expression of HRAS Q61R in MCF10A and MCF12A cells resulted in i) increased proliferation and transformation, ii) an irregular growth pattern in 3D organotypic cell cultures, iii) partial loss of the epithelial phenotype, and iv) acquisition of myoepithelial differentiation, which was more overt in PIK3CA-mutant MCF10A cells. HRAS Q61Rinduced hyperactivation of the PI3K pathway, but both PI3K and MAPK pathways likely contributed to the RAS-mediated proliferation, which was completely arrested by combined AKT and MEK inhibition.

Conclusion: AMEs are phenotypically and genetically heterogeneous. Whilst pathogenic mutations in PI3K pathway-related genes occur across the spectrum of lesions, HRAS Q61 hotspot mutations are restricted to ER-negative AMEs. Our genomic and functional analyses indicate that HRAS Q61 mutations are driver events in the pathogenesis of ER-negative AMEs and, in conjunction with mutant PIK3CA, may lead to the acquisition of myoepithelial differentiation in breast epithelial cells.

Citation Format: Geyer FC, Li A, Papanastasiou AD, Smith A, Selenica P, Burke KA, Edelweiss M, Wen H-C, Piscuoglio S, Schultheis AM, Martelotto LG, Pareja F, Kumar R, Brandes A, Lozada J, Macedo GS, Muenst S, Terracciano LM, Jungbluth A, Foschini MP, Wen HY, Brogi E, Palazzo J, Rubin BP, Ng CKY, Norton L, Varga Z, Ellis IO, Rakha E, Chandarlapatty S, Weigelt B, Reis-Filho JS. Estrogen receptor-negative breast adenomyoepitheliomas are driven by co-occurring HRAS hotspot and PI3K pathway gene mutations: A genetic and functional analysis [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr PD4-13.

MYBL1 rearrangements and MYB amplification in breast adenoid cystic carcinomas lacking the MYB-NFIB fusion gene

Breast adenoid cystic carcinoma (AdCC), a rare type of triple-negative breast cancer, has been shown to be driven by MYB pathway activation, most often underpinned by the MYB-NFIB fusion gene. Alternative genetic mechanisms, such as MYBL1 rearrangements, have been reported in MYB-NFIB-negative salivary gland AdCCs. Here we report on the molecular characterization by massively parallel sequencing of four breast AdCCs lacking the MYB-NFIB fusion gene. In two cases, we identified MYBL1 rearrangements (MYBL1-ACTN1 and MYBL1-NFIB), which were associated with MYBL1 overexpression. A third AdCC harboured a high-level MYB amplification, which resulted in MYB overexpression at the mRNA and protein levels. RNA-sequencing and whole-genome sequencing revealed no definite alternative driver in the fourth AdCC studied, despite high levels of MYB expression and the activation of pathways similar to those activated in MYB-NFIB-positive AdCCs. In this case, a deletion encompassing the last intron and part of exon 15 of MYB, including the binding site of ERG-1, a transcription factor that may downregulate MYB, and the exon 15 splice site, was detected. In conclusion, we demonstrate that MYBL1 rearrangements and MYB amplification probably constitute alternative genetic drivers of breast AdCCs, functioning through MYBL1 or MYB overexpression. These observations emphasize that breast AdCCs probably constitute a convergent phenotype, whereby activation of MYB and MYBL1 and their downstream targets can be driven by the MYB-NFIB fusion gene, MYBL1 rearrangements, MYB amplification, or other yet to be identified mechanisms. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Diverse BRCA1 and BRCA2 Reversion Mutations in Circulating Cell-Free DNA of Therapy-Resistant Breast or Ovarian Cancer

Purpose: Resistance to platinum-based chemotherapy or PARP inhibition in germline BRCA1 or BRCA2 mutation carriers may occur through somatic reversion mutations or intragenic deletions that restore BRCA1 or BRCA2 function. We assessed whether BRCA1/2 reversion mutations could be identified in circulating cell-free DNA (cfDNA) of patients with ovarian or breast cancer previously treated with platinum and/or PARP inhibitors.Experimental Design: cfDNA from 24 prospectively accrued patients with germline BRCA1 or BRCA2 mutations, including 19 patients with platinum-resistant/refractory ovarian cancer and five patients with platinum and/or PARP inhibitor pretreated metastatic breast cancer, was subjected to massively parallel sequencing targeting all exons of 141 genes and all exons and introns of BRCA1 and BRCA2 Functional studies were performed to assess the impact of the putative BRCA1/2 reversion mutations on BRCA1/2 function.Results: Diverse and often polyclonal putative BRCA1 or BRCA2 reversion mutations were identified in cfDNA from four patients with ovarian cancer (21%) and from two patients with breast cancer (40%). BRCA2 reversion mutations were detected in cfDNA prior to PARP inhibitor treatment in a patient with breast cancer who did not respond to treatment and were enriched in plasma samples after PARP inhibitor therapy. Foci formation and immunoprecipitation assays suggest that a subset of the putative reversion mutations restored BRCA1/2 function.Conclusions: Putative BRCA1/2 reversion mutations can be detected by cfDNA sequencing analysis in patients with ovarian and breast cancer. Our findings warrant further investigation of cfDNA sequencing to identify putative BRCA1/2 reversion mutations and to aid the selection of patients for PARP inhibition therapy. Clin Cancer Res; 23(21); 6708-20. ©2017 AACR.

Diverse BRCA1 and BRCA2 Reversion Mutations in Circulating Cell-Free DNA of Therapy-Resistant Breast or Ovarian Cancer

Purpose: Resistance to platinum-based chemotherapy or PARP inhibition in germline BRCA1 or BRCA2 mutation carriers may occur through somatic reversion mutations or intragenic deletions that restore BRCA1 or BRCA2 function. We assessed whether BRCA1/2 reversion mutations could be identified in circulating cell-free DNA (cfDNA) of patients with ovarian or breast cancer previously treated with platinum and/or PARP inhibitors.Experimental Design: cfDNA from 24 prospectively accrued patients with germline BRCA1 or BRCA2 mutations, including 19 patients with platinum-resistant/refractory ovarian cancer and five patients with platinum and/or PARP inhibitor pretreated metastatic breast cancer, was subjected to massively parallel sequencing targeting all exons of 141 genes and all exons and introns of BRCA1 and BRCA2 Functional studies were performed to assess the impact of the putative BRCA1/2 reversion mutations on BRCA1/2 function.Results: Diverse and often polyclonal putative BRCA1 or BRCA2 reversion mutations were identified in cfDNA from four patients with ovarian cancer (21%) and from two patients with breast cancer (40%). BRCA2 reversion mutations were detected in cfDNA prior to PARP inhibitor treatment in a patient with breast cancer who did not respond to treatment and were enriched in plasma samples after PARP inhibitor therapy. Foci formation and immunoprecipitation assays suggest that a subset of the putative reversion mutations restored BRCA1/2 function.Conclusions: Putative BRCA1/2 reversion mutations can be detected by cfDNA sequencing analysis in patients with ovarian and breast cancer. Our findings warrant further investigation of cfDNA sequencing to identify putative BRCA1/2 reversion mutations and to aid the selection of patients for PARP inhibition therapy. Clin Cancer Res; 23(21); 6708-20. ©2017 AACR.

Abstract 3379: Massively parallel sequencing analysis of breast adenomyoepitheliomas reveals the heterogeneity of the disease and identifies a subset driven by HRAS hotspot mutations

Adenomyoepithelioma (AME) is a rare biphasic proliferative breast lesion, which may resemble salivary gland epithelial-myoepithelial carcinomas (EMCs). Most AMEs have an indolent clinical course, but malignant transformation and local and distant recurrences have been reported. We sought to define the mutational landscape of AMEs and investigate the functional impact of recurrent likely pathogenic mutations identified in these tumors. Nineteen AMEs were subjected to whole-exome massively parallel sequencing (MPS, n=7) or targeted capture MPS using MSK-IMPACT assay (n=12). Somatic genetic alterations and the cancer cell fraction of mutations were defined using state-of-the-art bioinformatics algorithms. Selected genes (i.e. HRAS and PIK3CA) were subjected to Sanger sequencing in a series of 17 additional AMEs (total n=36). Non-tumorigenic mammary epithelial cells (i.e. MCF10A, MCF10A with the PIK3CAH1047R mutation and MCF12A), which are estrogen receptor (ER)-negative, were utilized for 2D and 3D functional studies. Of 36 cases, 22 were ER-positive and 14 were ER-negative. MPS analysis revealed a low mutation burden and HRASQ61 and PIK3CA hotspot mutations in 6/19 (32%) and 11/19 (58%) AMEs, respectively. All HRASQ61 and all but one PIK3CA mutations were clonal. ER-positive and ER-negative AMEs were fundamentally histologically and genetically distinct. Whilst ER-positive AMEs displayed recurrent PIK3CA mutations (50%, 11/22) but lacked HRAS mutations, ER-negative AMEs displayed, in addition to PIK3CA mutations (57%, 8/14), recurrent HRASQ61 mutations (57%, 8/14). HRASQ61 mutations co-occurred with PIK3CA mutations (50%, 4/8), PIK3R1 deletions (12.5%, 1/8) and/or CDKN2A homozygous deletions (25%, 2/8). HRASQ61 mutations, but not PIK3CA mutations, were significantly associated with ER-negativity (100% vs 21%), concurrent carcinoma (50% vs 7%), axillary metastases (38% vs 0%), high proliferation (63% vs 4%), necrosis (63% vs 11%) and nuclear pleomorphism (75% vs 29%). In vitro forced HRASQ61R expression in MCF10A and MCF12A cells resulted in increased proliferation and transformation. In 3D organotypic cell cultures, forced HRASQ61R resulted in a highly disorganized growth pattern, a partial loss of epithelial phenotype and acquisition of aberrant myoepithelial differentiation, which was more overt in PIK3CA-mutant MCF10A cells. In conclusion, AMEs are phenotypically and genetically heterogeneous. Whilst PIK3CA hotspot mutations occur across the spectrum of lesions, HRASQ61 hotspot mutations are restricted to ER-negative AMEs, which should arguably be classified as breast EMCs. Our genomic and functional analyses are consistent with the notion that HRASQ61 mutations are driver events in the pathogenesis of ER-negative AMEs and may be sufficient for the acquisition of myoepithelial differentiation in breast cells.

Citation Format: Felipe C. Geyer, Kathleen A. Burke, Anqi Li, Anastasios D. Papanastatiou, Fresia Pareja, Anne S. Schulteis, Charlotte K. Ng, Salvatore Piscuoglio, Marcia Edelweiss, Luciano G. Martelotto, Pier Selenica, Maria R. Filippo, Gabriel S. Macedo, Achim Jungbluth, Hannah Y. Wen, Juan Palazzo, Zsuzsanna Varga, Emad Rakha, Ian O. Ellis, Brian Rubin, Britta Weigelt, Jorge S. Reis-Filho. Massively parallel sequencing analysis of breast adenomyoepitheliomas reveals the heterogeneity of the disease and identifies a subset driven by HRAS hotspot mutations [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 3379. doi:10.1158/1538-7445.AM2017-3379

Bi-allelic alterations in DNA repair genes underpin homologous recombination DNA repair defects in breast cancer

Homologous recombination (HR) DNA repair-deficient (HRD) breast cancers have been shown to be sensitive to DNA repair targeted therapies. Burgeoning evidence suggests that sporadic breast cancers, lacking germline BRCA1/BRCA2 mutations, may also be HRD. We developed a functional ex vivo RAD51-based test to identify HRD primary breast cancers. An integrated approach examining methylation, gene expression, and whole-exome sequencing was employed to ascertain the aetiology of HRD. Functional HRD breast cancers displayed genomic features of lack of competent HR, including large-scale state transitions and specific mutational signatures. Somatic and/or germline genetic alterations resulting in bi-allelic loss-of-function of HR genes underpinned functional HRD in 89% of cases, and were observed in only one of the 15 HR-proficient samples tested. These findings indicate the importance of a comprehensive genetic assessment of bi-allelic alterations in the HR pathway to deliver a precision medicine-based approach to select patients for therapies targeting tumour-specific DNA repair defects. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Generation of conditional oncogenic chromosomal translocations using CRISPR-Cas9 genomic editing and homology-directed repair

Chromosomal rearrangements encoding oncogenic fusion proteins are found in a wide variety of malignancies. The use of programmable nucleases to generate specific double-strand breaks in endogenous loci, followed by non-homologous end joining DNA repair, has allowed several of these translocations to be generated as constitutively expressed fusion genes within a cell population. Here, we describe a novel approach that combines CRISPR-Cas9 technology with homology-directed repair to engineer, capture, and modulate the expression of chromosomal translocation products in a human cell line. We have applied this approach to the genetic modelling of t(11;22)(q24;q12) and t(11;22)(p13;q12), translocation products of the EWSR1 gene and its 3' fusion partners FLI1 and WT1, present in Ewing's sarcoma and desmoplastic small round cell tumour, respectively. Our innovative approach allows for temporal control of the expression of engineered endogenous chromosomal rearrangements, and provides a means to generate models to study tumours driven by fusion genes. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Abstract S2-02: The landscape of somatic genetic alterations in BRCA1 and BRCA2 breast cancers

Introduction: Women carrying BRCA1 or BRCA2 germline mutations have a 45-80% lifetime risk of developing breast cancer (BC). BRCA1 and BRCA2 are perceived as bona fide tumor suppressor genes, whereby bi-allelic inactivation in tumor cells is required for tumorigenesis. Recent studies have indicated that loss of heterozygosity (LOH) of the wild-type allele of BRCA1 may be heterogeneous and constitute a late event. Therefore, additional somatic events prior to full BRCA1/2 inactivation may be required for tumorigenesis. Given that the somatic events that result in the development of BRCA1/2-BCs and their chronology are not understood, here we sought to define the genomic landscape of BRCA1/2-BCs and whether LOH of BRCA1/2 wild-type allele and/or mutations affecting additional tumor suppressor genes would be clonal or subclonal in these cancers.

Methods: We retrieved 29 BRCA1-BCs and 10 BRCA2-BCs from the Pathology Departments of the authors' institutions. DNA extracted from microdissected tumor and normal breast samples was subjected to targeted capture massively parallel sequencing using either the MSK-IMPACT assay or an assay targeting all exons of 254 genes recurrently mutated in BC or related to DNA repair. Somatic single nucleotide variants, small insertions and deletions and copy number alterations affecting genes present in both sequencing assays (111 genes) were defined using state-of-the-art bioinformatics algorithms. ABSOLUTE and FACETS were employed to define clonal (i.e. present virtually in 100% of the cancer cells of a given case) and subclonal mutations and the presence of LOH of the BRCA1 and BRCA2 wild-type alleles.

Results: Our analysis revealed bi-allelic inactivation of BRCA1 in 28 of 29 BRCA1-BCs (93% harbored LOH of the BRCA1 wild-type allele and 3% harbored a second somatic BRCA1 pathogenic mutation). The only BRCA1-BC lacking bi-allelic inactivation of BRCA1 was an estrogen receptor-positive lobular carcinoma, lacking genomic features consistent with homologous recombination DNA repair defects, diagnosed at 62 years of age. Bi-allelic inactivation of BRCA2 was found in all cases (100% of harbored LOH of the BRCA2 wild-type allele). A clonal somatic 'second hit' resulting in bi-allelic inactivation of BRCA1 or BRCA2 was detected in 76% and 100% of BRCA1-BCs and BRCA2-BCs, respectively. In BRCA1-BCs, TP53 mutations were detected in 76% of cases, and these mutations were found to be clonal in 58% of cases. The repertoire of somatic mutations affecting BRCA1-BCs included clonal somatic mutations or homozygous deletions of known tumor suppressor genes, such as PTEN, RB1, CDKN2A and NF1. In contrast, only 10% of the BRCA2-BCs harbored TP53 somatic mutations. Though clonal somatic mutations in several cancer genes were detected, 40% of BRCA2-BCs had no mutations affecting the cancer genes analyzed.

Conclusions: Bi-allelic inactivation of BRCA1 and BRCA2 are frequent events in BRCA1-BCs and BRCA2-BCs, respectively. In a subset of BRCA1-BCs, however, the second 'hit' appeared to be subclonal, whereas mutations affecting TP53 and other tumor suppressor genes were clonal, supporting the notion that at least in a subset of these tumors, loss of the wild-type allele of BRCA1 may be preceded by inactivation of another tumor suppressor gene.

Citation Format: Geyer FC, Burke KA, Macedo GS, Piscuoglio S, Ng CK, Martelotto LG, Papanastatiou AD, De Filippo MR, Schultheis AM, Brogi E, Robson M, Wen YH, Weigelt B, Schnitt SJ, Tung N, Reis-Filho JS. The landscape of somatic genetic alterations in BRCA1 and BRCA2 breast cancers [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr S2-02.

Whole-genome single-cell copy number profiling from formalin-fixed paraffin-embedded samples

A substantial proportion of tumors consist of genotypically distinct subpopulations of cancer cells. This intratumor genetic heterogeneity poses a substantial challenge for the implementation of precision medicine. Single-cell genomics constitutes a powerful approach to resolve complex mixtures of cancer cells by tracing cell lineages and discovering cryptic genetic variations that would otherwise be obscured in tumor bulk analyses. Because of the chemical alterations that result from formalin fixation, single-cell genomic approaches have largely remained limited to fresh or rapidly frozen specimens. Here we describe the development and validation of a robust and accurate methodology to perform whole-genome copy-number profiling of single nuclei obtained from formalin-fixed paraffin-embedded clinical tumor samples. We applied the single-cell sequencing approach described here to study the progression from in situ to invasive breast cancer, which revealed that ductal carcinomas in situ show intratumor genetic heterogeneity at diagnosis and that these lesions may progress to invasive breast cancer through a variety of evolutionary processes.

Genetic events in the progression of adenoid cystic carcinoma of the breast to high-grade triple-negative breast cancer

Adenoid cystic carcinoma of the breast is a rare histological type of triple-negative breast cancer with an indolent clinical behavior, often driven by the MYB-NFIB fusion gene. Here we sought to define the repertoire of somatic genetic alterations in two adenoid cystic carcinomas associated with high-grade triple-negative breast cancer. The different components of each case were subjected to copy number profiling and massively parallel sequencing targeting all exons and selected regulatory and intronic regions of 488 genes. Reverse transcription PCR and fluorescence in situ hybridization were employed to investigate the presence of the MYB-NFIB translocation. The MYB-NFIB fusion gene was detected in both adenoid cystic carcinomas and their associated high-grade triple-negative breast cancer components. Although the distinct components of both cases displayed similar patterns of gene copy number alterations, massively parallel sequencing analysis revealed intratumor genetic heterogeneity. In case 1, progression from the trabecular adenoid cystic carcinoma to the high-grade triple-negative breast cancer was found to involve clonal shifts with enrichment of mutations affecting EP300, NOTCH1, ERBB2 and FGFR1 in the high-grade triple-negative breast cancer. In case 2, a clonal KMT2C mutation was present in the cribriform adenoid cystic carcinoma, solid adenoid cystic carcinoma and high-grade triple-negative breast cancer components, whereas a mutation affecting MYB was present only in the solid and high-grade triple-negative breast cancer areas and additional three mutations targeting STAG2, KDM6A and CDK12 were restricted to the high-grade triple-negative breast cancer. In conclusion, adenoid cystic carcinomas of the breast with high-grade transformation are underpinned by the MYB-NFIB fusion gene and, akin to other forms of cancer, may be constituted by a mosaic of cancer cell clones at diagnosis. The progression from adenoid cystic carcinoma to high-grade triple-negative breast cancer of no special type may involve the selection of neoplastic clones and/or the acquisition of additional genetic alterations.

IDH2 Mutations Define a Unique Subtype of Breast Cancer with Altered Nuclear Polarity

Solid papillary carcinoma with reverse polarity (SPCRP) is a rare breast cancer subtype with an obscure etiology. In this study, we sought to describe its unique histopathologic features and to identify the genetic alterations that underpin SPCRP using massively parallel whole-exome and targeted sequencing. The morphologic and immunohistochemical features of SPCRP support the invasive nature of this subtype. Ten of 13 (77%) SPCRPs harbored hotspot mutations at R172 of the isocitrate dehydrogenase IDH2, of which 8 of 10 displayed concurrent pathogenic mutations affecting PIK3CA or PIK3R1 One of the IDH2 wild-type SPCRPs harbored a TET2 Q548* truncating mutation coupled with a PIK3CA H1047R hotspot mutation. Functional studies demonstrated that IDH2 and PIK3CA hotspot mutations are likely drivers of SPCRP, resulting in its reversed nuclear polarization phenotype. Our results offer a molecular definition of SPCRP as a distinct breast cancer subtype. Concurrent IDH2 and PIK3CA mutations may help diagnose SPCRP and possibly direct effective treatment. Cancer Res; 76(24); 7118-29. ©2016 AACR.

Abstract 100: The landscape of somatic genetic alterations in BRCA1 and BRCA2 breast cancers

Background: BRCA1 or BRCA2 germline mutations account for a substantial proportion of hereditary breast cancers. The majority of BRCA1-associated breast cancers are of triple-negative phenotype (estrogen receptor (ER)-, progesterone receptor (PR)- and HER2-negative) and harbor TP53 somatic mutations. BRCA2-associated cancers are less homogeneous and often ER-positive. Systematic analyses of the profiles of somatic genetic alterations in BRCA1- and BRCA2-breast cancers have yet to be reported. Here we sought to determine the repertoire of somatic mutations and copy number alterations (CNAs) in breast cancers occurring in patients with BRCA1 or BRCA2 germline mutations..

Methods: Eleven BRCA1- and five BRCA2-associated breast cancers were microdissected. DNA was extracted from microdissected frozen tumor-normal pairs and subjected to targeted capture massively parallel sequencing using the MSK-IMPACT platform, which includes all exons and selected introns and regulatory regions of 410 key cancer genes. Somatic single nucleotide variants, small insertions and deletions, CNAs and the cancer cell fraction (CCF) of each alteration were defined using state-of-the-art bioinformatics algorithms.

Results: 82% (9/11) and 18% (2/11) of BRCA1-breast cancers were triple-negative and ER-positive/HER2-negative, respectively. All BRCA2-cancers were ER-positive, of which one was HER2-positive. Sequencing analysis revealed a median of four (2-11) and two (0-6) non-synonymous somatic mutations in BRCA1- and BRCA2-breast cancers, respectively. Within BRCA1-breast cancers all but one (10/11, 91%) harbored TP53 clonal somatic mutations, and clonal somatic loss of the BRCA1 wild-type allele was found in 8 of these cases. The BRCA1-breast cancer lacking these alterations was ER-positive and the only case harboring a PIK3CA mutation. Additional clonal mutations identified in BRCA1-breast cancers included somatic mutations affecting EGFR and RB1. Subclonal mutations in known cancer genes (e.g. GATA3 and PTEN) were also identified, suggesting intra-tumor genetic heterogeneity. All BRCA2-breast cancers analyzed displayed clonal loss of the wild-type of BRCA2, but no gene was found to be recurrently mutated.

Conclusions: BRCA1- and BRCA2-breast cancers are both characterized by clonal inactivation of the BRCA1 and BRCA2 wild-type alleles, respectively, which likely constitute early genetic events. Within BRCA1-breast cancers, TP53 mutations are highly recurrent and clonal, and may precede somatic loss of the BRCA1 wild-type allele, as the latter was subclonal in one case harboring a clonal TP53 somatic mutation.

Citation Format: Gabriel S. Macedo, Kathleen A. Burke, Salvatore Piscuoglio, Charlotte K. Y. Ng, Felipe C. Geyer, Luciano G. Martelotto, Anastasios D. Papanastatiou, Maria R. De Filippo, Anne M. Schultheis, Edi Brogi, Mark E. Robson, Hannah Y. Wen, Britta Weigelt, Jorge S. Reis-Filho. The landscape of somatic genetic alterations in BRCA1 and BRCA2 breast cancers. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 100.

Abstract 138: Solid papillary carcinoma with reverse polarization are driven by IDH2 and PI3K pathway mutations

Background: Solid papillary carcinoma with reverse polarization (SPCRP) is a rare breast cancer subtype, whose most striking morphologic feature is the presence of back-to-back columnar epithelial cells with nuclei situated in an apical rather than in the normal basal location. We sought to characterize the morphologic and genetic landscape of SPCRP and determine whether it represents a distinct subtype of breast cancer underpinned by highly recurrent disease-specific genetic alterations.

Material and Methods: Archival sections of 13 SPCRPs were subjected to microdissection. DNA samples extracted from microdissected tumor and matched normal samples were subjected to whole exome sequencing (n = 2) on an Illumina HiSeq2000, and DNA samples from tumors were further subjected to targeted massively parallel sequencing (n = 1), SNaPshot profiling (n = 7) and/or Sanger sequencing. State-of-the-art bioinformatics methods were used to define somatic mutations. Non-malignant breast epithelial MCF10A cells with and without somatic PIK3CA H1047R knock-in were used to assess the functional impact of the mutations identified by sequencing of SPCRPs in two- and three-dimensional cell culture systems.

Results: Ten of 13 (77%) SPCRPs were found to harbor IDH2 R172 hotspot mutations, and eight of these cases had concurrent pathogenic mutations affecting PIK3CA or PIK3R1, in particular PIK3CA H1047R hotspot mutations. Functional monolayer and three-dimensional cell culture studies demonstrated that IDH2 R172 and PIK3CA mutations constitute likely drivers of this tumor and contribute to its unusual nuclear polarization phenotype.

Conclusions: Our results suggest that SPCRP is a distinct subtype of breast cancer underpinned by IDH2 hotspot mutations in conjunction with mutations affecting canonical genes of the PI3K pathway. Given that IDH2 hotspot mutations have not been described in breast cancer to date, these may be used as a diagnostic ancillary marker for SPCRPs. Furthermore, IDH2 mutations may serve as a potential therapeutic target in SPCRP patients with disseminated disease.

Citation Format: Sarah Chiang, Britta Weigelt, Huei-Chi Wen, Salvatore Piscuoglio, Luciano G. Martelotto, Charlotte K Y Ng, Jörg Balss, Gabrielle Baker, Kimberly S. Cole, Andreas von Deimling, Julie M. Batten, Jonathan D. Marotti, Hwei-Choo Soh, Benjamin L. McCalip, Raymond S. Lim, Kalliopi P. Siziopikou, Randi Burke, Song Lu, Xiaolong Liu, Tarek Hammour, Edi Brogi, A John Iafrate, Jorge S. Reis-Filho, Stuart J. Schnitt. Solid papillary carcinoma with reverse polarization are driven by IDH2 and PI3K pathway mutations. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 138.

The Genomic Landscape of Male Breast Cancers

PURPOSE

Male breast cancer is rare, and its genomic landscape has yet to be fully characterized. Lacking studies in men, treatment of males with breast cancer is extrapolated from results in females with breast cancer. We sought to define whether male breast cancers harbor somatic genetic alterations in genes frequently altered in female breast cancers.

EXPERIMENTAL DESIGN

All male breast cancers were estrogen receptor-positive, and all but two were HER2-negative. Fifty-nine male breast cancers were subtyped by immunohistochemistry, and tumor-normal pairs were microdissected and subjected to massively parallel sequencing targeting all exons of 241 genes frequently mutated in female breast cancers or DNA-repair related. The repertoires of somatic mutations and copy number alterations of male breast cancers were compared with that of subtype-matched female breast cancers.

RESULTS

Twenty-nine percent and 71% of male breast cancers were immunohistochemically classified as luminal A-like or luminal B-like, respectively. Male breast cancers displayed a heterogeneous repertoire of somatic genetic alterations that to some extent recapitulated that of estrogen receptor (ER)-positive/HER2-negative female breast cancers, including recurrent mutations affecting PIK3CA (20%) and GATA3 (15%). ER-positive/HER2-negative male breast cancers, however, less frequently harbored 16q losses, and PIK3CA and TP53 mutations than ER-positive/HER2-negative female breast cancers. In addition, male breast cancers were found to be significantly enriched for mutations affecting DNA repair-related genes.

CONCLUSIONS

Male breast cancers less frequently harbor somatic genetic alterations typical of ER-positive/HER2-negative female breast cancers, such as PIK3CA and TP53 mutations and losses of 16q, suggesting that at least a subset of male breast cancers are driven by a distinct repertoire of somatic changes. Given the genomic differences, caution may be needed in the application of biologic and therapeutic findings from studies of female breast cancers to male breast cancers. Clin Cancer Res; 22(16); 4045-56. ©2016 AACR.

Massively parallel sequencing of phyllodes tumours of the breast reveals actionable mutations, and TERT promoter hotspot mutations and TERT gene amplification as likely drivers of progression

Phyllodes tumours (PTs) are breast fibroepithelial lesions that are graded based on histological criteria as benign, borderline or malignant. PTs may recur locally. Borderline PTs and malignant PTs may metastasize to distant sites. Breast fibroepithelial lesions, including PTs and fibroadenomas, are characterized by recurrent MED12 exon 2 somatic mutations. We sought to define the repertoire of somatic genetic alterations in PTs and whether these may assist in the differential diagnosis of these lesions. We collected 100 fibroadenomas, 40 benign PTs, 14 borderline PTs and 22 malignant PTs; six, six and 13 benign, borderline and malignant PTs, respectively, and their matched normal tissue, were subjected to targeted massively parallel sequencing (MPS) using the MSK-IMPACT sequencing assay. Recurrent MED12 mutations were found in 56% of PTs; in addition, mutations affecting cancer genes (eg TP53, RB1, SETD2 and EGFR) were exclusively detected in borderline and malignant PTs. We found a novel recurrent clonal hotspot mutation in the TERT promoter (-124 C>T) in 52% and TERT gene amplification in 4% of PTs. Laser capture microdissection revealed that these mutations were restricted to the mesenchymal component of PTs. Sequencing analysis of the entire cohort revealed that the frequency of TERT alterations increased from benign (18%) to borderline (57%) and to malignant PTs (68%; p < 0.01), and TERT alterations were associated with increased levels of TERT mRNA (p < 0.001). No TERT alterations were observed in fibroadenomas. An analysis of TERT promoter sequencing and gene amplification distinguished PTs from fibroadenomas with a sensitivity and a positive predictive value of 100% (CI 95.38-100%) and 100% (CI 85.86-100%), respectively, and a sensitivity and a negative predictive value of 39% (CI 28.65-51.36%) and 68% (CI 60.21-75.78%), respectively. Our results suggest that TERT alterations may drive the progression of PTs, and may assist in the differential diagnosis between PTs and fibroadenomas. Copyright © 2015 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Lack of PRKD2 and PRKD3 kinase domain somatic mutations in PRKD1 wild-type classic polymorphous low-grade adenocarcinomas of the salivary gland

AIMS

Polymorphous low-grade adenocarcinoma (PLGA) is the second most common intra-oral salivary gland malignancy. The vast majority of PLGAs harbour a PRKD1 E710D hot-spot somatic mutation or somatic rearrangements of PRKD1, PRKD2 or PRKD3. Given the kinase domain homology among PRKD1, PRKD2 and PRKD3, we sought to define whether PLGAs lacking PRKD1 somatic mutations or PRKD gene family rearrangements would be driven by somatic mutations affecting the kinase domains of PRKD2 or PRKD3.

METHODS AND RESULTS

DNA was extracted from eight microdissected PLGAs lacking PRKD1 somatic mutations or PRKD gene family rearrangements. Samples were thoroughly centrally reviewed, microdissected and subjected to Sanger sequencing of the kinase domains of the PRKD2 and PRKD3 genes. None of the PLGAs lacking PRKD1 somatic mutations or PRKD gene family rearrangements harboured somatic mutations in the kinase domains of the PRKD2 or PRKD3 genes.

CONCLUSION

PLGAs lacking PRKD1 somatic mutations or PRKD gene family rearrangements are unlikely to harbour somatic mutations in the kinase domains of PRKD2 or PRKD3. Further studies are warranted to define the driver genetic events in this subgroup of PLGAs.

Uterine adenosarcomas are mesenchymal neoplasms

Uterine adenosarcomas (UAs) are biphasic lesions composed of a malignant mesenchymal (ie stromal) component and an epithelial component. UAs are generally low-grade and have a favourable prognosis, but may display sarcomatous overgrowth (SO), which is associated with a worse outcome. We hypothesized that, akin to breast fibroepithelial lesions, UAs are mesenchymal neoplasms in which clonal somatic genetic alterations are restricted to the mesenchymal component. To characterize the somatic genetic alterations in UAs and to test this hypothesis, we subjected 20 UAs to a combination of whole-exome (n = 6), targeted capture (n = 13) massively parallel sequencing (MPS) and/or RNA sequencing (n = 6). Only three genes, FGFR2, KMT2C and DICER1, were recurrently mutated, all in 2/19 cases; however, 26% (5/19) and 21% (4/19) of UAs harboured MDM2/CDK4/HMGA2 and TERT gene amplification, respectively, and two cases harboured fusion genes involving NCOA family members. Using a combination of laser-capture microdissection and in situ techniques, we demonstrated that the somatic genetic alterations detected by MPS were restricted to the mesenchymal component. Furthermore, mitochondrial DNA sequencing of microdissected samples revealed that epithelial and mesenchymal components of UAs were clonally unrelated. In conclusion, here we provide evidence that UAs are genetically heterogeneous lesions and mesenchymal neoplasms.

Abstract 2971: Whole exome sequencing reveals heterogeneity within lobular carcinoma in situ (LCIS) and clonal selection in the progression to malignant lesions

INTRODUCTION: Lobular carcinoma in situ (LCIS) is considered both a risk factor and non-obligate precursor of low-grade estrogen receptor-positive breast cancer. We sought to define the mutational repertoire, subclone complexity and heterogeneity of LCIS, and whether invasive lobular carcinomas (ILCs) would stem from specific subclones within a LCIS.

METHODS: Patients with a history of LCIS undergoing therapeutic or prophylactic mastectomy were prospectively enrolled in an IRB approved protocol. Frozen tissue blocks were collected, screened for lesions of interest (LCIS, ductal carcinoma in situ (DCIS), ILC, invasive ductal carcinoma (IDC)) and subject to microdissection and DNA extraction. Matched germline DNA was available for all cases. Whole exome sequencing was performed on a HiSeq2000 (Illumina) and data were aligned to the reference human genome hg19 and processed using GATK. SNVs were called using MuTect, and indels were called using a combination of Varscan and Strelka. Purity and ploidy estimates were calculated by ABSOLUTE. Clonal frequencies were estimated using Pyclone.

RESULTS: 30 LCIS, 10 ILCs, 6 IDCs and 7 DCIS from 15 patients qualified for data analysis, resulting in 18 LCIS-ILC pairs, 22 LCIS-LCIS pairs, 12 LCIS-DCIS pairs, and 14 LCIS-IDC pairs for comparison. 9/18 (50%) LCIS-ILC pairs and 8/22 (36%) LCIS-LCIS pairs were clonally related, supported by several shared mutations (median 18, range 7-81 for LCIS-ILC; median 14, range 5-22 for LCIS-LCIS). All related LCIS-ILC pairs and 6/8 related LCIS-LCIS pairs shared a pathogenic CDH1 mutation; 75% of related LCIS-ILC pairs also shared a PIK3CA hotspot mutation. 7/12 (58%) LCIS-DCIS pairs were found to be clonally related but the number of shared mutations was generally lower than that found in LCIS-ILC pairs (median 9, range 2-11). No evidence of a clonal relationship was found in any of the LCIS-IDC pairs tested. Clonal composition analysis revealed that samples of LCIS display intra-lesion genetic heterogeneity in the form of the presence of a minor clone in 70% of cases. In one case, the LCIS minor subclone (∼15%) constituted the major clone in the ILC and in another case the LCIS minor subclone constituted the major clone in the associated DCIS. The majority of the clonally related lesions were located in the same quadrant of the breast, however evidence of clonality was found in 5 LCIS-LCIS and 3 LCIS-DCIS pairs located in separate quadrants of the breast.

CONCLUSIONS: Intra-lesion genetic heterogeneity is a common phenomenon in LCIS. The dominant clone of a LCIS may not always be the clone directly involved in the progression to malignancy. The spatial relationships of clonally related lesions in this study suggest that anatomy does not always infer clonality, as lesions located in separate quadrants of the breast may be clonally related.

Citation Format: Michail Schizas, Rita A. Sakr, Britta Weigelt, Charlotte KY Ng, Jose Victor S. Carniello, Dilip Giri, Salvatore Piscuoglio, Luciano G. Martelotto, Russell Towers, Victor P. Andrade, Raymond Lim, David B. Solit, Jorge S. Reis-Filho, Tari A. King. Whole exome sequencing reveals heterogeneity within lobular carcinoma in situ (LCIS) and clonal selection in the progression to malignant lesions. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2971. doi:10.1158/1538-7445.AM2015-2971

Abstract 4735: Single-cell sequencing from formalin-fixed paraffin-embedded breast cancers: a powerful tool to address intratumor genetic heterogeneity

Background: At diagnosis, tumors often consist of multiple, genotypically distinct cell populations. Genome sequencing of single cells has opened new avenues of investigation, yet the application of this technology remains limited to fresh/frozen tissue samples. Here, we describe and validate a method to perform single-cell massively parallel sequencing using DNA extracted from individual nuclei from FFPE tumor samples.

Methods: Tissue sections are deparaffinized and the area of interest microdissected, DNA is reverse-crosslinked and the extracellular matrix is digested. Extracted single nuclei are then FACS-sorted, lysed and the whole genome repaired using a broad-spectrum DNA-repair enzyme cocktail. Repaired DNA is whole genome amplified (WGA) using GenomePlex WGA4 (Sigma-Aldrich) with modifications for heavily damaged genomic templates. Illumina sequencing libraries are generated using standard approaches followed by multiplex sequencing on an Illumina HiSeq2000. To obtain copy number (CN) profiles, single-cell sequencing data are mapped to the reference genome with PCR duplicates removed. Uniquely mapped reads are allocated and counted in genomic intervals of variable length (bins) with CN states roughly proportional to the number of allocated sequencing reads. Bin counts are then normalized on the basis of GC-content, segmented and transformed to CN values to identify long contiguous regions of equivalent CN.

Results: As a proof-of-principle we performed single-cell CN profiling on two aneuploid synchronous ductal carcinomas in situ (DCIS)/invasive ductal carcinomas (IDCs), where both FFPE and frozen material was available. Using our novel methodology, 24 single nuclei were multiplexed/lesion. An average of 3 million reads/cell (average coverage of 0.1X/cell) was obtained, providing sufficient data to infer CN profiles of single cells accurately. These data were successfully employed to identify non-neoplastic cells and distinct clonal lineages of neoplastic cells within each lesion. The data obtained from the analysis of FFPE samples were concordant with those obtained from the analysis of matched frozen samples.

Conclusions: We developed a robust procedure to perform single-cell massively parallel sequencing of individual nuclei isolated from FFPE samples, providing the opportunity to unlock pathology archives for studies aiming to catalog and dissect the biological and clinical relevance of intra-tumor genetic heterogeneity.

Citation Format: Luciano G. Martelotto, Rita A. Sakr, Timour Baslan, Linda Rodgers, Hilary Cox, Jude Kendall, Tari A. King, Britta Weigelt, James Hicks, Jorge S. Reis-Filho. Single-cell sequencing from formalin-fixed paraffin-embedded breast cancers: a powerful tool to address intratumor genetic heterogeneity. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4735. doi:10.1158/1538-7445.AM2015-4735

Abstract 2989: Intra-tumor heterogeneity and clonal changes in the progression of DCIS to invasiveness: Combined tumor bulk and single cell analysis

INTRODUCTION: Ductal carcinoma in situ (DCIS) is a clonal intraductal proliferation of epithelial cells which acts as a non-obligate precursor of invasive breast cancer (IBC), yet the genetic events leading to the acquisition of invasive behavior remain unclear. We hypothesize that DCIS is composed of mosaics of genetically diverse tumor cell clones, and that the process of invasion is an evolutionary bottleneck. To test this hypothesis we performed a detailed characterization of the repertoire of genetic alterations and intra-tumor genetic heterogeneity in synchronously diagnosed DCIS and IBC using NextGen sequencing of bulk tumor and single cells.

METHODS: DNA extracted from fresh frozen, microdissected DCIS, IBC and adjacent normal tissue was subjected to whole exome sequencing on an Illumina HiSeq2000. Reads were aligned to the reference human genome hg19. Single nucleotide variants (SNVs) were called by MuTect, and gene copy number alterations were determined using VarScan2. Single nuclei were isolated from 100μm serial sections microdissected to separate DCIS and IBC. Individual nuclei were FACS-sorted into a 96-well plate, lysed and whole genome amplified. Amplified DNA samples were barcoded, pooled and sequenced on a HiSeq2000. Single cell sequencing data were mapped to the reference genome and uniquely mapped reads were allocated into bins, normalized, segmented and CN values generated.

RESULTS: In 6 cases of synchronous DCIS and IBC, a median of 41 and 47 non-synonymous mutations were found in each component, respectively. The somatic mutations identified in both DCIS and adjacent IBC components affected known driver breast cancer genes, including AKT1, PIK3CA, GATA3, MAP2K4 and TP53. Interestingly, we also found mutations restricted to either DCIS or IBC: ATRX (IBC-3); ALK and PKD2 (IBC-5); ESR1 (DCIS-5). The gene copy number profiles of matched DCIS and IBC were similar in all 6 pairs, however we also identified gene copy number alterations restricted either to DCIS or IBC: 1q gain (IBC-5); 3p and 3q losses (DCIS-6); 12p homozygous deletion (DCIS-4). Single cell sequencing of two cases (3 and 4) revealed that the majority of cells from both DCIS and IBC were derived from a common precursor lineage with shared copy number losses and gains. Both sets of DCIS-IBC pairs in these cases displayed elements of a subclonal structure with dominant clones alongside genetically diverse derivatives as well as genetic heterogeneity reflected in variable copy number alterations and non-modal clones.

CONCLUSION: Synchronous DCIS and IBC share founder genetic events, but also harbor somatic genetic alterations restricted to either the DCIS or IBC components, demonstrating that although DCIS is a precursor of IBC, intra-tumor genetic heterogeneity is present at the DCIS stage. Changes in clonal composition likely take place in the progression from DCIS to IBC.

Citation Format: Rita A. Sakr, Luciano G. Martelotto, Timour Baslan, Charlotte KY Ng, Jude Kendall, Linda Rodgers, Hilary Cox, Mike Riggs, Sean D'Itali, Asya Stepansky, Narciso Olvera, Tari A. King, Britta Weigelt, Jorge S. Reis-Filho, James Hicks. Intra-tumor heterogeneity and clonal changes in the progression of DCIS to invasiveness: Combined tumor bulk and single cell analysis. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2989. doi:10.1158/1538-7445.AM2015-2989

The repertoire of somatic genetic alterations of acinic cell carcinomas of the breast: an exploratory, hypothesis-generating study

Acinic cell carcinoma (ACC) of the breast is a rare form of triple-negative (that is, oestrogen receptor-negative, progesterone receptor-negative, HER2-negative) salivary gland-type tumour displaying serous acinar differentiation. Despite its triple-negative phenotype, breast ACCs are reported to have an indolent clinical behaviour. Here, we sought to define whether ACCs have a mutational repertoire distinct from that of other triple-negative breast cancers (TNBCs). DNA was extracted from microdissected formalin-fixed, paraffin-embedded sections of tumour and normal tissue from two pure and six mixed breast ACCs. Each tumour component of the mixed cases was microdissected separately. Tumour and normal samples were subjected to targeted capture massively parallel sequencing targeting all exons of 254 genes, including genes most frequently mutated in breast cancer and related to DNA repair. Selected somatic mutations were validated by targeted amplicon resequencing and Sanger sequencing. Akin to other forms of TNBC, the most frequently mutated gene found in breast ACCs was TP53 (one pure and six mixed cases). Additional somatic mutations affecting breast cancer-related genes found in ACCs included PIK3CA, MTOR, CTNNB1, BRCA1, ERBB4, ERBB3, INPP4B, and FGFR2. Copy number alteration analysis revealed complex patterns of gains and losses similar to those of common forms of TNBCs. Of the mixed cases analysed, identical somatic mutations were found in the acinic and the high-grade non-acinic components in two out of four cases analysed, providing evidence of their clonal relatedness. In conclusion, breast ACCs display the hallmark somatic genetic alterations found in high-grade forms of TNBC, including complex patterns of gene copy number alterations and recurrent TP53 mutations. Furthermore, we provide circumstantial genetic evidence to suggest that ACCs may constitute the substrate for the development of more aggressive forms of triple-negative disease.

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.

Resolving quandaries: basaloid adenoid cystic carcinoma or breast cylindroma? The role of massively parallel sequencing

AIMS

The aims of this study were to perform a whole-exome sequencing analysis of a breast cylindroma and to investigate the role of molecular analyses in the differentiation between breast cylindroma, a benign tumour that displays MYB expression, and CYLD gene mutations, and its main differential diagnosis, the breast solid-basaloid adenoid cystic carcinoma, a malignant tumour that is characterized by the presence of the MYB-NFIB fusion gene and MYB overexpression.

METHODS AND RESULTS

A 66-year-old female underwent quadrantectomy after an irregular dense shadow was discovered in the right breast at the screening mammogram. Histologically, the tumour displayed features suggestive of a solid-basaloid variant of adenoid cystic carcinoma with a differential diagnosis of cylindroma. Fluorescence in situ hybridization, reverse transcription-polymerase chain reaction, immunohistochemistry and whole-exome sequencing revealed absence of the MYB-NFIB fusion gene, low levels of MYB protein expression and a clonal somatic CYLD splice site mutation associated with loss of heterozygosity of the wild-type allele.

CONCLUSIONS

The results of the histological, immunohistochemical and molecular analyses were consistent with a diagnosis of breast cylindroma, providing a proof-of-principle that the integration of histopathological and molecular approaches can help to differentiate between a low-malignant potential and a benign breast tumour of triple-negative phenotype.

MED12 somatic mutations in fibroadenomas and phyllodes tumours of the breast

AIMS

Somatic mutations in exon 2 of the mediator complex subunit 12 (MED12) gene have been identified in 60% of breast fibroadenomas (FAs). The aim of this study was to define whether phyllodes tumours (PTs) would harbour MED12 somatic mutations in a way akin to FAs.

METHODS AND RESULTS

A collection of 73 fibroepithelial tumours (including 26 FAs, 25 benign PTs, nine borderline PTs and 13 malignant PTs) from 64 patients was retrieved from the authors' institution. Sections from formalin-fixed paraffin-embedded (FFPE) blocks were microdissected to ensure an enrichment in neoplastic stromal elements of >70%. DNA samples extracted from tumour and matched normal tissues were subjected to Sanger sequencing of exon 2 of the MED12 gene. MED12 exon 2 somatic mutations, including 28 somatic single nucleotide variants and 19 insertions and deletions, were found in 65%, 88%, 78% and 8% of FAs, benign PTs, borderline PTs and malignant PTs, respectively. Malignant PTs harboured MED12 exon 2 somatic mutations significantly less frequently than FAs, benign and borderline PTs.

CONCLUSIONS

Although MED12 exon 2 somatic mutations probably constitute the driver genetic event of most FAs, benign and borderline PTs, our results suggest that the majority of malignant PTs may be driven by other genetic/epigenetic alterations.

Intra-tumor genetic heterogeneity and alternative driver genetic alterations in breast cancers with heterogeneous HER2 gene amplification

Background

HER2 is overexpressed and amplified in approximately 15% of invasive breast cancers, and is the molecular target and predictive marker of response to anti-HER2 agents. In a subset of these cases, heterogeneous distribution of HER2 gene amplification can be found, which creates clinically challenging scenarios. Currently, breast cancers with HER2 amplification/overexpression in just over 10% of cancer cells are considered HER2-positive for clinical purposes; however, it is unclear as to whether the HER2-negative components of such tumors would be driven by distinct genetic alterations. Here we sought to characterize the pathologic and genetic features of the HER2-positive and HER2-negative components of breast cancers with heterogeneous HER2 gene amplification and to define the repertoire of potential driver genetic alterations in the HER2-negative components of these cases.

Results

We separately analyzed the HER2-negative and HER2-positive components of 12 HER2 heterogeneous breast cancers using gene copy number profiling and massively parallel sequencing, and identified potential driver genetic alterations restricted to the HER2-negative cells in each case. In vitro experiments provided functional evidence to suggest that BRF2 and DSN1 overexpression/amplification, and the HER2 I767M mutation may be alterations that compensate for the lack of HER2 amplification in the HER2-negative components of HER2 heterogeneous breast cancers.

Conclusions

Our results indicate that even driver genetic alterations, such as HER2 gene amplification, can be heterogeneously distributed within a cancer, and that the HER2-negative components are likely driven by genetic alterations not present in the HER2-positive components, including BRF2 and DSN1 amplification and HER2 somatic mutations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0657-6) contains supplementary material, which is available to authorized users.

Abstract P2-03-08: Mutational landscape of metaplastic breast carcinomas

Introduction: Metaplastic breast carcinoma (MBC) is an aggressive histologic type of breast cancer, which preferentially displays a triple-negative phenotype. These tumors are characterized by the presence of malignant cells exhibiting differentiation towards squamous epithelium or mesenchymal elements, including spindle, chondroid, osseous and rhabdoid differentiation. Unlike other rare histologic types of breast cancer such as adenoid cystic and secretory carcinomas, which are underpinned by the MYB-NFIB and ETV6-NTRK3 fusion genes respectively, pathognomonic genetic alterations have not been identified in MBC. It has been suggested, however, that the frequency of PIK3CA somatic mutations would be significantly higher in MBCs than in other forms of triple-negative disease. Here we sought to characterize the mutational landscape of MBCs by means of high-depth whole exome sequencing analysis.

Material and Methods: Twenty-one triple-negative MBCs were retrieved from the authors’ institutions. Representative sections from frozen blocks were microdissected to ensure tumor cell content greater than 50%. DNA samples extracted from microdissected tumor and matched peripheral blood leukocytes were subjected to high-depth (250x) whole exome sequencing on an Illumina GAIIx or HiSeq2000. Somatic point mutations were called using MuTect and somatic insertions and deletions (indels) were called using Strelka, Varscan2 and Haplotype Caller. Potentially pathogenic mutations were predicted using computational algorithms including PolyPhen-2, Mutation Taster, Mutation Assessor, CHASM and FATHMM. Significantly mutated genes were identified using MutSigCV. Pathway and network enrichment analysis of mutations was performed with Ingenuity Pathway Analysis and HOTNET. The genomic landscape of MBCs was compared with that of triple-negative breast cancers (TNBCs) analyzed as part of The Cancer Genome Atlas project.

Results: A mean of 135 somatic non-synonymous point mutations and indels were identified per MBC. The most frequently mutated gene was TP53, found in 12/21 cases (57%), and the only significantly mutated gene as defined by MutSigCV (q&lt;0.01). The repertoire of somatic mutations found in MBCs was qualitatively similar to that of TNBCs of no special type, and recurrently mutated genes were altered at similar frequencies in MBCs and TNBCs of no special type. When somatic mutations were annotated in pathways and networks, MBCs were found to have potentially pathogenic mutations affecting genes directly related to the PI3K pathway, including pathogenic non-synonymous mutations affecting PIK3CA, PIK3R1, PIK3R2, PIK3C2B, PIK3C2G and PTEN, significantly more frequently than TNBCs of no special type (10 out of 21 MBCs vs. 11 out of 62 TNBCs; Fisher's exact test p-value=0.0099).

Conclusion: The majority (57%) of MBCs harbored non-synonymous mutations affecting TP53. While the frequencies of mutations affecting recurrently mutated genes in MBCs are similar to those found in other forms of TNBCs, MBCs significantly more frequently harbor mutations affecting PI3K pathway-related genes than TNBCs of no special type.

Citation Format: Charlotte KY Ng, Britta Weigelt, Salvatore Piscuoglio, Y Hannah Wen, Maria R De Filippo, Luciano G Martelotto, Rachael Natrajan, Raymond Lim, Edi Brogi, Larry Norton, Anne Vincent-Salomon, Jorge S Reis-Filho. Mutational landscape of metaplastic breast carcinomas [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr P2-03-08.

Benchmarking mutation effect prediction algorithms using functionally validated cancer-related missense mutations

BACKGROUND

Massively parallel sequencing studies have led to the identification of a large number of mutations present in a minority of cancers of a given site. Hence, methods to identify the likely pathogenic mutations that are worth exploring experimentally and clinically are required. We sought to compare the performance of 15 mutation effect prediction algorithms and their agreement. As a hypothesis-generating aim, we sought to define whether combinations of prediction algorithms would improve the functional effect predictions of specific mutations.

RESULTS

Literature and database mining of single nucleotide variants (SNVs) affecting 15 cancer genes was performed to identify mutations supported by functional evidence or hereditary disease association to be classified either as non-neutral (n = 849) or neutral (n = 140) with respect to their impact on protein function. These SNVs were employed to test the performance of 15 mutation effect prediction algorithms. The accuracy of the prediction algorithms varies considerably. Although all algorithms perform consistently well in terms of positive predictive value, their negative predictive value varies substantially. Cancer-specific mutation effect predictors display no-to-almost perfect agreement in their predictions of these SNVs, whereas the non-cancer-specific predictors showed no-to-moderate agreement. Combinations of predictors modestly improve accuracy and significantly improve negative predictive values.

CONCLUSIONS

The information provided by mutation effect predictors is not equivalent. No algorithm is able to predict sufficiently accurately SNVs that should be taken forward for experimental or clinical testing. Combining algorithms aggregates orthogonal information and may result in improvements in the negative predictive value of mutation effect predictions.

Hotspot activating PRKD1 somatic mutations in polymorphous low-grade adenocarcinomas of the salivary glands

Polymorphous low-grade adenocarcinoma (PLGA) is the second most frequent type of malignant tumor of the minor salivary glands. We identified PRKD1 hotspot mutations encoding p.Glu710Asp in 72.9% of PLGAs but not in other salivary gland tumors. Functional studies demonstrated that this kinase-activating alteration likely constitutes a driver of PLGA.

Integrative genomic and transcriptomic characterization of papillary carcinomas of the breast

Papillary carcinoma (PC) is a rare type of breast cancer, which comprises three histologic subtypes: encapsulated PC (EPC), solid PC (SPC) and invasive PC (IPC). Microarray‐based gene expression and Affymetrix SNP 6.0 gene copy number profiling, and RNA‐sequencing revealed that PCs are luminal breast cancers that display transcriptomic profiles distinct from those of grade‐ and estrogen receptor (ER)‐matched invasive ductal carcinomas of no special type (IDC‐NSTs), and that the papillary histologic pattern is unlikely to be underpinned by a highly recurrent expressed fusion gene or a highly recurrent expressed mutation. Despite displaying similar patterns of gene copy number alterations, significant differences in the transcriptomic profiles of EPCs, SPCs and IPCs were found, and may account for their different histologic features.

Papillary carcinomas of the breast display distinctive transcriptomic profiles.

Proliferation‐related genes are expressed at low levels in papillary carcinomas.

Papillary carcinomas are unlikely to be underpinned by a highly recurrent fusion gene.

Papillary carcinomas are unlikely to be underpinned by a highly recurrent expressed mutation.

Abstract P4-04-08: Genomic and transcriptomic characterization of papillary carcinomas of the breast

Introduction: Papillary carcinomas (PCs) are a rare (&lt;1%) special histological type of breast cancer that often affects postmenopausal women, and has an overall favorable outcome. Based on their histological characteristics, these tumors are classified into three different subtypes, namely encapsulated papillary carcinomas (EPC), solid papillary carcinomas (SPC) and invasive papillary carcinomas (IPC). In this study, we sought i) to investigate whether PCs constitute a molecular entity distinct from grade- and ER-matched invasive ductal carcinomas of no special type (IDC-NST) at the transcriptomic level, ii) to investigate whether EPC, SPC and IPC display distinct transcriptomic profiles, iii) to characterize the repertoire of copy number alterations in the different subtypes of PC, and iv) to identify recurrent fusion genes that may be potential drivers of this disease.

Material and methods: DNA and RNA were extracted from microdissected PCs (4 SPCs, 5 IPCs and 7 EPCs) and grade- and ER-matched IDC-NSTs (RNA only). 16 PCs and 16 grade- and ER-matched IDC-NSTs were subjected to gene expression profiling using the Illumina Human HT-12 v4 platform. Genes differentially expressed between the PC subtypes were identified using SAM, and functional annotation of these genes was performed using DAVID. Intrinsic molecular subtypes were determined using the PAM50 single sample predictor. Copy number profiling was performed using Affymetrix Human SNP 6.0 arrays with DNA extracted from 16 PCs. In addition, 8 PCs (3 IPCs, 3 EPCs, 2 SPCs) were subjected to paired-end massively parallel RNA sequencing (Illumina GAIIx). Putative expressed fusion transcripts were identified using validated algorithms (i.e. deFuse and Chimerascan), and confirmed by reverse transcription PCR.

Results: PCs were preferentially of histological grade I/II (82%) and ER-positive (100%). Unsupervised analysis revealed that PC subtypes show a high degree of similarity at the transcriptomic level, and form clusters distinct from grade- and ER-matched IDC-NSTs. Compared with IDC-NSTs, PCs displayed reduced expression of genes related to cell motility, adhesion and extracellular matrix. PAM50 subtyping classified 87.5%, 50% and 100% of EPCs, SPCs and IPCs as of luminal subtypes, respectively. 12.5% of EPCs were classified as of basal-like subtype, and 50% of SPCs as of HER2-enriched subtype. At the genomic level, PC subtypes displayed similar patterns of gene copy number aberrations. Five in-frame fusion genes, USF1-CCDC38, MDC1-SFTA2, DLD-LMBR1, PDCL-DENND1A and SUGT1-NOL6, were identified and validated in PCs, however none of these were recurrent in the cases included in this study.

Conclusion: Our results demonstrate that the majority of PCs are of luminal subtype, and support the contention that at the transcriptomic level, PCs are distinct from grade- and ER-matched IDC-NSTs. Our findings also demonstrate that unlike some other histological special types of breast cancer, PCs are not underpinned by a highly recurrent expressed fusion gene.

Citation Information: Cancer Res 2013;73(24 Suppl): Abstract nr P4-04-08.