Phyllodes tumors (PTs) and fibroadenomas (FAs) are fibroepithelial lesions (FELs) of the breast. Although mutations affecting exon 2 of MED12 are highly recurrent in FAs and PTs, TERT promoter hotspot mutations are frequently found in PTs but are vanishingly rare in FAs. Malignant transformation of benign PTs is well-documented, but the progression from FA to PT remains a matter of contention. Here we report on the somatic genetic alterations in multiple ipsilateral synchronous FELs (three FAs, one benign PT, and one malignant PT) occurring in the same patient. DNA samples extracted from each tumor and matched normal tissue were subjected to targeted massively parallel sequencing using the Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets (MSK-IMPACT) assay. This analysis revealed MED12 mutations in all lesions. One FA and the benign PT harbored a MED12Gly44Val mutation, whereas another FA and the malignant PT displayed a MED12Gly44Asp mutation. The remaining FA had an independent distinct MED12Gly44Cys mutation. A formal clonality analysis suggested a clonal relationship between the FELs with identical MED12 mutations (P<0.05). A clonal TERT promoter hotspot mutation was identified exclusively in the malignant PT. The identification of distinct MED12 mutations in multifocal ipsilateral and synchronous FELs supports the notion that co-existing mammary fibroepithelial tumors can arise independently. Conversely, the co-existence of identical MED12 mutations indicates clonal relatedness among FAs and PTs, corroborating the hypothesis that FAs may constitute the substrate from which PTs develop. Our findings also support the notion that acquisition of TERT promoter mutations may drive the progression of FELs.

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

Metaplastic breast carcinoma (MBC) and uterine carcinosarcoma (UCS) are rare aggressive cancers, characterized by an admixture of adenocarcinoma and areas displaying mesenchymal/sarcomatoid differentiation. We sought to define whether MBCs and UCSs harbor similar patterns of genetic alterations, and whether the different histologic components of MBCs and UCSs are clonally related. Whole-exome sequencing (WES) data from MBCs (n = 35) and UCSs (n = 57, The Cancer Genome Atlas) were reanalyzed to define somatic genetic alterations, altered signaling pathways, mutational signatures, and genomic features of homologous recombination DNA repair deficiency (HRD). In addition, the carcinomatous and sarcomatous components of an additional cohort of MBCs (n = 11) and UCSs (n = 6) were microdissected separately and subjected to WES, and their clonal relatedness was assessed. MBCs and UCSs harbored recurrent genetic alterations affecting TP53, PIK3CA, and PTEN, similar patterns of gene copy number alterations, and an enrichment in alterations affecting the epithelial-to-mesenchymal transition (EMT)-related Wnt and Notch signaling pathways. Differences were observed, however, including a significantly higher prevalence of FAT3 and FAT1 somatic mutations in MBCs compared to UCSs, and conversely, UCSs significantly more frequently harbored somatic mutations affecting FBXW7 and PPP2R1A as well as HER2 amplification than MBCs. Genomic features of HRD and biallelic alterations affecting bona fide HRD-related genes were found to be more prevalent in MBCs than in UCSs. The distinct histologic components of MBCs and UCSs were clonally related in all cases, with the sarcoma component likely stemming from a minor subclone of the carcinoma component in the samples with interpretable chronology of clonal evolution. Despite the similar histologic features and pathways affected by genetic alterations, UCSs differ from MBCs on the basis of FBXW7 and PPP2R1A mutations, HER2 amplification, and lack of HRD, supporting the notion that these entities are more than mere phenocopies of the same tumor type in different anatomical sites.

Mono-allelic germline pathogenic variants in the Partner And Localizer of BRCA2 (PALB2) gene predispose to a high-risk of breast cancer development, consistent with the role of PALB2 in homologous recombination (HR) DNA repair. Here, we sought to define the repertoire of somatic genetic alterations in PALB2-associated breast cancers (BCs), and whether PALB2-associated BCs display bi-allelic inactivation of PALB2 and/or genomic features of HR-deficiency (HRD). Twenty-four breast cancer patients with pathogenic PALB2 germline mutations were analyzed by whole-exome sequencing (WES, n = 16) or targeted capture massively parallel sequencing (410 cancer genes, n = 8). Somatic genetic alterations, loss of heterozygosity (LOH) of the PALB2 wild-type allele, large-scale state transitions (LSTs) and mutational signatures were defined. PALB2-associated BCs were found to be heterogeneous at the genetic level, with PIK3CA (29%), PALB2 (21%), TP53 (21%), and NOTCH3 (17%) being the genes most frequently affected by somatic mutations. Bi-allelic PALB2 inactivation was found in 16 of the 24 cases (67%), either through LOH (n = 11) or second somatic mutations (n = 5) of the wild-type allele. High LST scores were found in all 12 PALB2-associated BCs with bi-allelic PALB2 inactivation sequenced by WES, of which eight displayed the HRD-related mutational signature 3. In addition, bi-allelic inactivation of PALB2 was significantly associated with high LST scores. Our findings suggest that the identification of bi-allelic PALB2 inactivation in PALB2-associated BCs is required for the personalization of HR-directed therapies, such as platinum salts and/or PARP inhibitors, as the vast majority of PALB2-associated BCs without PALB2 bi-allelic inactivation lack genomic features of HRD.

Tall cell carcinoma with reverse polarity is a rare subtype of breast carcinoma with solid and papillary growth and nuclear features reminiscent of those of the tall cell variant of papillary thyroid carcinoma. These tumors harbor recurrent IDH2 R172 hotspot mutations or TET2 mutations, co-occurring with mutations in PI3K pathway genes. Diagnosis of tall cell carcinomas with reverse polarity is challenging in view of their rarity and the range of differential diagnosis. We sought to determine the sensitivity and specificity of IDH2 R172 immunohistochemistry for the detection of IDH2 R172 hotspot mutations in this entity. We evaluated 14 tall cell carcinomas with reverse polarity (ten excision and five core needle biopsy specimens), 13 intraductal papillomas, 16 solid papillary carcinomas, and 5 encapsulated papillary carcinomas by Sanger sequencing of the IDH2 R172 hotspot locus and of exons 9 and 20 of PIK3CA, and by immunohistochemistry using monoclonal antibodies (11C8B1) to the IDH2 R172S mutation. The 14 tall cell carcinomas with reverse polarity studied harbored IDH2 R172 hotspot mutations, which co-occurred with PIK3CA hotspot mutations in 50% of cases. None of the other papillary neoplasms analyzed displayed IDH2 R172 mutations, however PIK3CA hotspot mutations were detected in 54% of intraductal papillomas, 6% of solid papillary carcinomas, and 20% of encapsulated papillary carcinomas tested. Immunohistochemical analysis with anti-IDH2 R172S antibodies (11C8B1) detected IDH2 R172 mutated protein in 93% (14/15) of tall cell carcinomas with reverse polarity samples including excision (n = 9/10) and core needle biopsy specimens (n = 5), whereas the remaining papillary neoplasms (n = 34) were negative. Our findings demonstrate that immunohistochemical analysis of IDH2 R172 is highly sensitive and specific for the detection of IDH2 R172 hotspot mutations, and likely suitable as a diagnostic tool in the evaluation of excision and core needle biopsy material of tall cell carcinomas with reverse polarity.

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.

Breast fibroepithelial lesions (fibroadenomas and phyllodes tumors) are underpinned by recurrent MED12 exon 2 mutations, which are more common in fibroadenomas and benign phyllodes tumors. TERT promoter hotspot mutations have been documented in phyllodes tumors, and found to be more frequent in borderline and malignant lesions. Several lines of evidence suggest that a subset of phyllodes tumors might arise from fibroadenomas. Here we sought to investigate the genetic differences between phyllodes tumors with fibroadenoma-like areas vs. those without. We retrieved data for 16 borderline/ malignant phyllodes tumors, including seven phyllodes tumors with fibroadenoma-like areas and nine phyllodes tumors without fibroadenoma-like areas, which had been previously subjected to targeted capture massively parallel sequencing. Whilst MED12 exon 2 mutations were significantly more frequent in tumors with fibroadenoma-like areas (71 vs. 11%), an enrichment in genetic alterations targeting bona fide cancer genes was found in those without fibroadenoma-like areas, in particular in EGFR mutations and amplifications (78 vs. 14%). No significant difference in the frequency of TERT genetic alterations was observed (71% in cases with fibroadenoma-like areas vs 56% in those without fibroadenoma-like areas). Our data suggest that the development of phyllodes tumors might follow two different evolutionary pathways: a MED12-mutant pathway that involves the progression from a fibroadenoma to a malignant phyllodes tumor; and a MED12-wild-type pathway, where malignant phyllodes tumors arise de novo through the acquisition of genetic alterations targeting cancer genes. Additional studies are warranted to confirm our observations and define whether the outcome differs between both pathways.

The American Society of Clinical Oncology/College of American Pathologists (ASCO/CAP) 2013 guidelines for HER2 assessment have increased the number of HER2 equivocal breast carcinomas following in situ hybridization reflex testing, that is, HER2 "double equivocal" (equivocal protein expression and equivocal gene copy number). Forty-five double-equivocal carcinomas were subjected to Prosigna analysis. Twenty-seven cases were investigated for the expression of genes found to be differentially expressed between estrogen receptor (ER)-positive/HER2-positive (N=22) and ER-positive/HER2-negative (N=22) control cases. Twenty-nine of the 45 cases were also analyzed by targeted sequencing using a panel of 14 genes. We then explored the pathologic complete response rates in an independent series of double-equivocal carcinoma patients treated with trastuzumab-containing chemotherapy. All cases were ER-positive, with a mean Ki67 of 28%. Double-equivocal carcinomas were predominantly luminal B (76%); 9 cases (20%) were luminal A, and 2 cases (4%) HER2-enriched. The majority (73%) showed a high risk of recurrence by Prosigna, even when the carcinomas were small (<2 cm), node-negative/micrometastatic, and/or grade 2. Double-equivocal carcinomas showed TP53 (6/29, 20%), PIK3CA (3/29, 10%), HER2 (1/29, 3%), and MAP2K4 (1/29, 3%) mutations. Compared with grade-matched ER-positive/HER2-negative breast carcinomas from METABRIC, double-equivocal carcinomas harbored more frequently TP53 mutations and less frequently PIK3CA mutations (P<0.05). No significant differences were observed with grade-matched ER-positive/HER2-positive carcinomas. Lower pathologic complete response rates were observed in double-equivocal compared with HER2-positive patients (10% vs. 60%, P=0.009). Double-equivocal carcinomas are preferentially luminal B and show a high risk of recurrence. A subset of these tumors can be labeled as HER2-enriched by transcriptomic analysis. HER2 mutations can be identified in HER2 double-equivocal cases.

Granular cell tumors (GCTs) are rare tumors that can arise in multiple anatomical locations, and are characterized by abundant intracytoplasmic granules. The genetic drivers of GCTs are currently unknown. Here, we apply whole-exome sequencing and targeted sequencing analysis to reveal mutually exclusive, clonal, inactivating somatic mutations in the endosomal pH regulators ATP6AP1 or ATP6AP2 in 72% of GCTs. Silencing of these genes in vitro results in impaired vesicle acidification, redistribution of endosomal compartments, and accumulation of intracytoplasmic granules, recapitulating the cardinal phenotypic characteristics of GCTs and providing a novel genotypic-phenotypic correlation. In addition, depletion of ATP6AP1 or ATP6AP2 results in the acquisition of oncogenic properties. Our results demonstrate that inactivating mutations of ATP6AP1 and ATP6AP2 are likely oncogenic drivers of GCTs and underpin the genesis of the intracytoplasmic granules that characterize them, providing a genetic link between endosomal pH regulation and tumorigenesis.

Uterine adenosarcomas and breast phyllodes tumors (PTs) are morphologically similar, being composed of stromal projections in a leaf-like fashion lined by epithelial cells. Here, we investigated whether their histologic similarities would be mirrored at the genetic level. The previously reported repertoires of somatic genetic alterations found in 19 adenosarcomas and 22 PTs (six benign, six borderline, and 10 malignant) were compared. PTs significantly more frequently displayed mutations affecting MED12, the TERT gene promoter and bona fide cancer genes, whereas adenosarcomas harbored a higher rate of MDM2/CDK4 and TERT gene amplifications. Pathway analyses based on the genes affected by somatic genetic alterations in these tumors indicated that Wnt signaling likely plays a role in the biology of adenosarcomas and benign/borderline PTs. In conclusion, despite the differences at the gene level, PTs and adenosarcomas share remarkable morphologic similarities and enrichment for somatic genetic alterations affecting Wnt pathway-related genes.

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

Adenomyoepithelioma of the breast is a rare tumor characterized by epithelial-myoepithelial differentiation, whose genetic underpinning is largely unknown. Here we show through whole-exome and targeted massively parallel sequencing analysis that whilst estrogen receptor (ER)-positive adenomyoepitheliomas display PIK3CA or AKT1 activating mutations, ER-negative adenomyoepitheliomas harbor highly recurrent codon Q61 HRAS hotspot mutations, which co-occur with PIK3CA or PIK3R1 mutations. In two- and three-dimensional cell culture models, forced expression of HRASQ61R in non-malignant ER-negative breast epithelial cells with or without a PIK3CAH1047R somatic knock-in results in transformation and the acquisition of the cardinal features of adenomyoepitheliomas, including the expression of myoepithelial markers, a reduction in E-cadherin expression, and an increase in AKT signaling. Our results demonstrate that adenomyoepitheliomas are genetically heterogeneous, and qualify mutations in HRAS, a gene whose mutations are vanishingly rare in common-type breast cancers, as likely drivers of ER-negative adenomyoepitheliomas.

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.

Acinic cell carcinoma is an indolent form of invasive breast cancer, whereas microglandular adenosis has been shown to be a neoplastic proliferation. Both entities display a triple-negative phenotype, and may give rise to and display somatic genomic alterations typical of high-grade triple-negative breast cancers. Here we report on a comparison of previously published data on eight carcinoma-associated microglandular adenosis and eight acinic cell carcinomas subjected to targeted massively parallel sequencing targeting all exons of 236 genes recurrently mutated in breast cancer and/or DNA repair-related. Somatic mutations, insertions/ deletions, and copy number alterations were detected using state-of-the-art bioinformatic algorithms. All cases were of triple-negative phenotype. A median of 4.5 (1-13) and 4.0 (1-7) non-synonymous somatic mutations per carcinoma-associated microglandular adenosis and acinic cell carcinoma were identified, respectively. TP53 was the sole highly recurrently mutated gene (75% in microglandular adenosis versus 88% in acinic cell carcinomas), and TP53 mutations were consistently coupled with loss of heterozygosity of the wild-type allele. Additional somatic mutations shared by both groups included those in BRCA1, PIK3CA, and INPP4B. Recurrent (n=2) somatic mutations restricted to microglandular adenosis or acinic cell carcinomas included those affecting PTEN and MED12 or ERBB4, respectively. No significant differences in the repertoire of somatic mutations were detected between microglandular adenosis and acinic cell carcinomas, and between this group of lesions and 77 triple-negative carcinomas from The Cancer Genome Atlas. Microglandular adenosis and acinic cell carcinomas, however, were genetically distinct from estrogen receptor-positive and/or HER2-positive breast cancers from The Cancer Genome Atlas. Our findings support the contention that microglandular adenosis and acinic cell carcinoma are part of the same spectrum of lesions harboring frequent TP53 somatic mutations, and likely represent low-grade forms of triple-negative disease with no/minimal metastatic potential, of which a subset has the potential to progress to high-grade triple-negative breast cancer.

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.

Breast cancers in BRCA1 mutation carriers frequently have a distinctive basal-like phenotype. It has been suggested that this results from an origin in basal breast epithelial stem cells. Here, we demonstrate that deleting Brca1 in mouse mammary epithelial luminal progenitors produces tumors that phenocopy human BRCA1 breast cancers. They also resemble the majority of sporadic basal-like breast tumors. However, directing Brca1 deficiency to basal cells generates tumors that express molecular markers of basal breast cancers but do not histologically resemble either human BRCA1 or the majority of sporadic basal-like breast tumors. These findings support a derivation of the majority of human BRCA1-associated and sporadic basal-like tumors from luminal progenitors rather than from basal stem cells. They also demonstrate that when target cells for transformation have the potential for phenotypic plasticity, tumor phenotypes may not directly reflect histogenesis. This has important implications for cancer prevention strategies.

Primary pleomorphic adenomas (PAs) and mucoepidermoid carcinomas (MECs) of the breast are vanishingly rare. Here we sought to determine whether breast PAs and MECs would be underpinned by the fusion genes reported to occur in their salivary gland counterparts. Our study included three breast PAs and one breast MEC, which were subjected to RNA sequencing (PAs, n = 2; MEC, n = 1) or to Archer FusionPlex sequencing (PA, n = 1). Our analyses revealed the presence of the HMGA2-WIF1 fusion gene in breast PA3, the CTNNB1-PLAG1 fusion gene in breast PA2, and the CRTC1-MAML2 fusion gene in the breast MEC analyzed (1/1). No oncogenic fusion genes were detected in breast PA1, and no additional oncogenic fusion genes were detected in the cases studied. The presence of the fusion genes identified was validated by fluorescence in situ hybridization (n = 1), reverse transcription-PCR (n = 1), or by both methods (n = 1). Taken together, our findings indicate that PAs and MECs arising in the breast resemble their salivary gland counterparts not only phenotypically but also at the genetic level. Furthermore, our data suggest that the molecular analysis of breast PAs and MECs might constitute a useful tool to aid in their differential diagnosis.

Polymorphous adenocarcinoma (PAC) and cribriform adenocarcinoma of (minor) salivary gland (CASG) are salivary gland tumors with overlapping spectrum of morphology. Whether these represent distinct entities or a histologic spectrum of the same tumor remains contentious. PACs harbor recurrent PRKD1 E710D hotspot mutations in >70% of cases, whereas 80% of CASGs display rearrangements involving PRKD1, PRKD2, or PRKD3 (PRKD1/2/3). We studied the molecular and morphologic features of 37 PACs/CASGs, seeking to identify the associations among genotype, histologic phenotype, and classification. DNA was subjected to Sanger sequencing analysis of the PRKD1 hotspot locus. Fluorescence in situ hybridization (FISH) analysis for PRKD1/2/3 was performed using dual-color break-apart probes. Tumors were classified into four categories as described previously: PAC, CASG, tumor with indeterminate features (TIF), and tumor with a predominant papillary pattern (TPPP). PRKD1 E710D hotspot mutations were identified in 56%, 20%, 43% and 0% of PACs, CASGs, TIFs, and TPPPs, respectively. FISH demonstrated PRKD1/2/3 rearrangements in 13%, 78%, 36%, and 75% of PACs, CASGs, TIFs, and TPPPs, respectively. Histologically, fusion-positive tumors were associated with a high percentage of papillary growth, low percentage of single filing arrangement, a propensity of base of tongue location, and frequent (50%) lymph node metastasis, compared with the mutation-related tumors which had negligible nodal metastasis risk. Our results demonstrated that (1) PACs/CASGs are underpinned by genetic alterations affecting PRKD genes; (2) despite the associations between PAC and PRKD1 hotspot mutations and CASG and PRKD1/2/3 fusion, such distinction is not absolute; and (3) there is of a novel genotypic-phenotypic association whereby fusion-positive tumors are usually located in the base of the tongue, show papillary architecture and have a high risk of nodal metastasis. Genetic analysis of PRKD genes appears to be useful characterizing this spectrum of tumors, not only histologically but also clinically identifying those tumors with high risk of nodal metastasis.