Penile squamous cell carcinoma is a rare disease, in which somatic genetic aberrations have yet to be characterized. We hypothesized that gene copy aberrations might correlate with human papillomavirus status and clinico-pathological features. We sought to determine the spectrum of gene copy number aberrations in a large series of PSCCs and to define their correlations with human papillomavirus, histopathological subtype, and tumor grade, stage and lymph node status. Seventy formalin-fixed, paraffin embedded penile squamous cell carcinomas were centrally reviewed by expert uropathologists. DNA was extracted from micro-dissected samples, subjected to PCR-based human papillomavirus assessment and genotyping (INNO-LiPA human papillomavirus Genotyping Extra Assay) and microarray-based comparative genomic hybridization using a 32K Bacterial Artificial Chromosome array platform. Sixty-four samples yielded interpretable results. Recurrent gains were observed in chromosomes 1p13.3-q44 (88%), 3p12.3-q29 (86%), 5p15.33-p11 (67%) and 8p12-q24.3 (84%). Amplifications of 5p15.33-p11 and 11p14.1-p12 were found in seven (11%) and four (6%) cases, respectively. Losses were observed in chromosomes 2q33-q37.3 (86%), 3p26.3-q11.1 (83%) and 11q12.2-q25 (81%). Although many losses and gains were similar throughout the cohort, there were small significant differences observed at specific loci, between human papillomavirus positive and negative tumors, between tumor types, and tumor grade and nodal status. These results demonstrate that despite the diversity of genetic aberrations in penile squamous cell carcinomas, there are significant correlations between the clinico-pathological data and the genetic changes that may play a role in disease natural history and progression and highlight potential driver genes, which may feature in molecular pathways for existing therapeutic agents.

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

Triple-negative breast cancers (BC) represent a heterogeneous subtype of BCs, generally associated with an aggressive clinical course and where targeted therapies are currently limited. Target validation studies for all BC subtypes have largely employed established BC cell lines, which have proven to be effective tools for drug discovery.

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.

Oesophageal cancer is the seventh most common cause of cancer-related death worldwide. Disease relapse is frequent and treatment options are limited.

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.

Mutations in the pioneer transcription factor FOXA1 are a hallmark of estrogen receptor-positive (ER+) breast cancers. Examining FOXA1 in ∼5,000 breast cancer patients identifies several hotspot mutations in the Wing2 region and a breast cancer-specific mutation SY242CS, located in the third β strand. Using a clinico-genomically curated cohort, together with breast cancer models, we find that FOXA1 mutations associate with a lower response to aromatase inhibitors. Mechanistically, Wing2 mutations display increased chromatin binding at ER loci upon estrogen stimulation, and an enhanced ER-mediated transcription without changes in chromatin accessibility. In contrast, SY242CS shows neomorphic properties that include the ability to open distinct chromatin regions and activate an alternative cistrome and transcriptome. Structural modeling predicts that SY242CS confers a conformational change that mediates stable binding to a non-canonical DNA motif. Taken together, our results provide insights into how FOXA1 mutations perturb its function to dictate cancer progression and therapeutic response.

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.

Traditional genetic testing for patients with gastrointestinal stromal tumors (GISTs) focus on those with syndromic features. To assess whether expanded genetic testing of GIST patients could identify hereditary cancer predisposition, we analyzed matched tumor-germline sequencing results from 103 patients with GISTs over a 6-year period. Germline pathogenic/likely pathogenic (P/LP) variants in GIST-associated genes (SDHA, SDHB, SDHC, NF1, KIT) were identified in 69% of patients with KIT/PDGFRA-wildtype GISTs, 63% of whom did not have any personal or family history of syndromic features. To evaluate the frequency of somatic versus germline variants identified in tumor-only sequencing of GISTs, we analyzed 499 de-identified tumor-normal pairs. P/LP variants in certain genes (e.g., BRCA1/2, SDHB) identified in tumor-only sequencing of GISTs were almost exclusively germline in origin. Our results provide guidance for genetic testing of GIST patients and indicate that germline testing should be offered to all patients with KIT/PDGFRA-wildtype GISTs regardless of their history of syndromic features.

Small cell carcinoma (SCC) of the uterine cervix is a rare and aggressive form of neuroendocrine carcinoma, which resembles small cell lung cancer (SCLC) in its histology and poor survival rate. Here, we sought to define the genetic underpinning of SCCs of the uterine cervix and compare their mutational profiles with those of human papillomavirus (HPV)-positive head and neck squamous cell carcinomas, HPV-positive cervical carcinomas, and SCLCs using publicly available data. Using a combination of whole-exome and targeted massively parallel sequencing, we found that the nine uterine cervix SCCs, which were HPV18-positive (n = 8) or HPV16-positive (n = 1), harbored a low mutation burden, few copy number alterations, and other than TP53 in two cases no recurrently mutated genes. The majority of mutations were likely passenger missense mutations, and only few affected previously described cancer-related genes. Using RNA-sequencing, we identified putative viral integration sites on 18q12.3 and on 8p22 in two SCCs of the uterine cervix. The overall nonsilent mutation rate of uterine cervix SCCs was significantly lower than that of SCLCs, HPV-driven cervical adeno- and squamous cell carcinomas, or HPV-positive head and neck squamous cell carcinomas. Unlike SCLCs, which are reported to harbor almost universal TP53 and RB1 mutations and a dominant tobacco smoke-related signature 4, uterine cervix SCCs rarely harbored mutations affecting these genes (2/9, 22% TP53; 0% RB1) and displayed a dominant aging (67%) or APOBEC mutational signature (17%), akin to HPV-driven cancers, including cervical adeno- and squamous cell carcinomas and head and neck squamous cell carcinomas. Taken together, in contrast to SCLCs, which are characterized by highly recurrent TP53 and RB1 alterations, uterine cervix SCCs were positive for HPV leading to inactivation of the suppressors p53 and RB, suggesting that these SCCs are convergent phenotypes.

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.

Sclerosing stromal tumor (SST) of the ovary is a rare type of sex cord-stromal tumor (SCST), whose genetic underpinning is currently unknown. Here, using whole-exome, targeted capture and RNA-sequencing, we report recurrent FHL2-GLI2 fusion genes in 65% (17/26) of SSTs and other GLI2 rearrangements in additional 15% (4/26) SSTs, none of which are detected in other types of SCSTs (n = 48) or common cancer types (n = 9,950). The FHL2-GLI2 fusions result in transcriptomic activation of the Sonic Hedgehog (SHH) pathway in SSTs. Expression of the FHL2-GLI2 fusion in vitro leads to the acquisition of phenotypic characteristics of SSTs, increased proliferation, migration and colony formation, and SHH pathway activation. Targeted inhibition of the SHH pathway results in reversal of these oncogenic properties, indicating its role in the pathogenesis of SSTs. Our results demonstrate that the FHL2-GLI2 fusion is likely the oncogenic driver of SSTs, defining a genotypic-phenotypic correlation in ovarian neoplasms.

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.

Late-stage relapse (LSR) in patients with breast cancer (BC) occurs more than five years and up to 10 years after initial treatment and has less than 30% 5-year relative survival rate. Long non-coding RNAs (lncRNAs) play important roles in BC yet have not been studied in LSR BC. Here, we identify 1127 lncRNAs differentially expressed in LSR BC via transcriptome sequencing and analysis of 72 early-stage and 24 LSR BC patient tumors. Decreasing expression of the most up-regulated lncRNA, LINC00355, in BC and MCF7 long-term estrogen deprived cell lines decreases cellular invasion and proliferation. Subsequent mechanistic studies show that LINC00355 binds to MENIN and changes occupancy at the CDKN1B promoter to decrease p27Kip. In summary, this is a key study discovering lncRNAs in LSR BC and LINC00355 association with epigenetic regulation and proliferation in BC.

Cyclin-dependent kinases 4 and 6 (CDK4/6) represent a major therapeutic vulnerability for breast cancer. The kinases are clinically targeted via ATP competitive inhibitors (CDK4/6i); however, drug resistance commonly emerges over time. To understand CDK4/6i resistance, we surveyed over 1,300 breast cancers and identified several genetic alterations (e.g., FAT1, PTEN, or ARID1A loss) converging on upregulation of CDK6. Mechanistically, we demonstrate CDK6 causes resistance by inducing and binding CDK inhibitor INK4 proteins (e.g., p18INK4C). In vitro binding and kinase assays together with physical modeling reveal that the p18INK4C-cyclin D-CDK6 complex occludes CDK4/6i binding while only weakly suppressing ATP binding. Suppression of INK4 expression or its binding to CDK6 restores CDK4/6i sensitivity. To overcome this constraint, we developed bifunctional degraders conjugating palbociclib with E3 ligands. Two resulting lead compounds potently degraded CDK4/6, leading to substantial antitumor effects in vivo, demonstrating the promising therapeutic potential for retargeting CDK4/6 despite CDK4/6i resistance. SIGNIFICANCE: CDK4/6 kinase activation represents a common mechanism by which oncogenic signaling induces proliferation and is potentially targetable by ATP competitive inhibitors. We identify a CDK6-INK4 complex that is resilient to current-generation inhibitors and develop a new strategy for more effective inhibition of CDK4/6 kinases.This article is highlighted in the In This Issue feature, p. 275.

Mucinous carcinomas can arise in any organ with epithelial cells that produce mucus. While mucinous tumors from different organs are histologically similar, it remains to be elucidated whether they share molecular alterations. Here we analyzed a total of 902 patients across six cancer types by comparing mucinous and non-mucinous samples, integrating text mining of pathology reports, gene expression, methylation, mutational and copy-number profiling. We found that, in addition to genes involved in mucin processing and secretion, MUC2 up-regulation is a multi-cancer biomarker of mucinous histology and is regulated by DNA methylation in colorectal, breast and stomach cancer. The majority of carcinomas with mucinous differentiation had fewer DNA copy-number alterations than non-mucinous tumors. The tumor mutational burden was lower in breast and lung with mucinous differentiation compared to their non-mucinous counterparts. We found several differences in the frequency of oncogenic gene and pathway alterations between mucinous and non-mucinous carcinomas, including a lower frequency of p53 pathway alterations in colorectal and lung cancer, and a lower frequency of PI-3-Kinase/Akt pathway alterations in breast and stomach cancer with mucinous differentiation. This study shows that carcinomas with mucinous differentiation originating from different organs share transcriptomic and genomic similarities. These results might pave the way for a more biologically relevant taxonomy for these rare cancers.

Physiologic activation of estrogen receptor α (ERα) is mediated by estradiol (E2) binding in the ligand-binding pocket of the receptor, repositioning helix 12 (H12) to facilitate binding of coactivator proteins in the unoccupied coactivator binding groove. In breast cancer, activation of ERα is often observed through point mutations that lead to the same H12 repositioning in the absence of E2. Through expanded genetic sequencing of breast cancer patients, we identified a collection of mutations located far from H12 but nonetheless capable of promoting E2-independent transcription and breast cancer cell growth. Using machine learning and computational structure analyses, this set of mutants was inferred to act distinctly from the H12-repositioning mutants and instead was associated with conformational changes across the ERα dimer interface. Through both in vitro and in-cell assays of full-length ERα protein and isolated ligand-binding domain, we found that these mutants promoted ERα dimerization, stability, and nuclear localization. Point mutations that selectively disrupted dimerization abrogated E2-independent transcriptional activity of these dimer-promoting mutants. The results reveal a distinct mechanism for activation of ERα function through enforced receptor dimerization and suggest dimer disruption as a potential therapeutic strategy to treat ER-dependent cancers.

Cyclin dependent kinase 4/6 (CDK4/6) inhibitors (CDK4/6i) are effective in breast cancer; however, drug resistance is frequently encountered and poorly understood. We conducted a genomic analysis of 348 estrogen receptor-positive (ER+) breast cancers treated with CDK4/6i and identified loss-of-function mutations affecting FAT1 and RB1 linked to drug resistance. FAT1 loss led to marked elevations in CDK6, the suppression of which restored sensitivity to CDK4/6i. The induction of CDK6 was mediated by the Hippo pathway with accumulation of YAP and TAZ transcription factors on the CDK6 promoter. Genomic alterations in other Hippo pathway components were also found to promote CDK4/6i resistance. These findings uncover a tumor suppressor function of Hippo signaling in ER+ breast cancer and establish FAT1 loss as a mechanism of resistance to CDK4/6i.

Aromatase inhibitors (AIs) are a vital component of estrogen receptor positive (ER+) breast cancer treatment. De novo and acquired resistance, however, is common. The aims of this study were to relate patterns of copy number aberrations to molecular and proliferative response to AIs, to study differences in the patterns of copy number aberrations between breast cancer samples pre- and post-AI neoadjuvant therapy, and to identify putative biomarkers for resistance to neoadjuvant AI therapy using an integrative analysis approach.

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.