Extra-cranial malignant rhabdoid tumors (MRTs) and cranial atypical teratoid RTs (ATRTs) are heterogeneous pediatric cancers driven primarily by SMARCB1 loss. To understand the genome-wide molecular relationships between MRTs and ATRTs, we analyze multi-omics data from 140 MRTs and 161 ATRTs. We detect similarities between the MYC subgroup of ATRTs (ATRT-MYC) and extra-cranial MRTs, including global DNA hypomethylation and overexpression of HOX genes and genes involved in mesenchymal development, distinguishing them from other ATRT subgroups that express neural-like features. We identify five DNA methylation subgroups associated with anatomical sites and SMARCB1 mutation patterns. Groups 1, 3, and 4 exhibit cytotoxic T cell infiltration and expression of immune checkpoint regulators, consistent with a potential role for immunotherapy in rhabdoid tumor patients.

Esthesioneuroblastoma (ENB) is a rare cancer of the olfactory mucosa, with no established molecular stratification to date. We report similarities of ENB with tumors arising in the neural crest and perform integrative analysis of these tumors. We propose a molecular-based subtype classification of ENB as basal or neural, both of which have distinct pathological, transcriptomic, proteomic, and immune features. Among the basal subtype, we uncovered an IDH2 R172 mutant-enriched subgroup (∼35%) harboring a CpG island methylator phenotype reminiscent of IDH2 mutant gliomas. Compared with the basal ENB methylome, the neural ENB methylome shows genome-wide reprogramming with loss of DNA methylation at the enhancers of axonal guidance genes. Our study reveals insights into the molecular pathogenesis of ENB and provides classification information of potential therapeutic relevance.

This integrated, multiplatform PanCancer Atlas study co-mapped and identified distinguishing molecular features of squamous cell carcinomas (SCCs) from five sites associated with smoking and/or human papillomavirus (HPV). SCCs harbor 3q, 5p, and other recurrent chromosomal copy-number alterations (CNAs), DNA mutations, and/or aberrant methylation of genes and microRNAs, which are correlated with the expression of multi-gene programs linked to squamous cell stemness, epithelial-to-mesenchymal differentiation, growth, genomic integrity, oxidative damage, death, and inflammation. Low-CNA SCCs tended to be HPV(+) and display hypermethylation with repression of TET1 demethylase and FANCF, previously linked to predisposition to SCC, or harbor mutations affecting CASP8, RAS-MAPK pathways, chromatin modifiers, and immunoregulatory molecules. We uncovered hypomethylation of the alternative promoter that drives expression of the ΔNp63 oncogene and embedded miR944. Co-expression of immune checkpoint, T-regulatory, and Myeloid suppressor cells signatures may explain reduced efficacy of immune therapy. These findings support possibilities for molecular classification and therapeutic approaches.

PBRM1 is the second most commonly mutated gene after VHL in clear cell renal cell carcinoma (ccRCC). However, the biological consequences of PBRM1 mutations for kidney tumorigenesis are unknown. Here, we find that kidney-specific deletion of Vhl and Pbrm1, but not either gene alone, results in bilateral, multifocal, transplantable clear cell kidney cancers. PBRM1 loss amplified the transcriptional outputs of HIF1 and STAT3 incurred by Vhl deficiency. Analysis of mouse and human ccRCC revealed convergence on mTOR activation, representing the third driver event after genetic inactivation of VHL and PBRM1. Our study reports a physiological preclinical ccRCC mouse model that recapitulates somatic mutations in human ccRCC and provides mechanistic and therapeutic insights into PBRM1 mutated subtypes of human ccRCC.

RNA sequencing (RNA-seq) detects estrogen receptor alpha gene (ESR1) fusion transcripts in estrogen receptor-positive (ER+) breast cancer, but their role in disease pathogenesis remains unclear. We examined multiple ESR1 fusions and found that two, both identified in advanced endocrine treatment-resistant disease, encoded stable and functional fusion proteins. In both examples, ESR1-e6>YAP1 and ESR1-e6>PCDH11X, ESR1 exons 1-6 were fused in frame to C-terminal sequences from the partner gene. Functional properties include estrogen-independent growth, constitutive expression of ER target genes, and anti-estrogen resistance. Both fusions activate a metastasis-associated transcriptional program, induce cellular motility, and promote the development of lung metastasis. ESR1-e6>YAP1- and ESR1-e6>PCDH11X-induced growth remained sensitive to a CDK4/6 inhibitor, and a patient-derived xenograft (PDX) naturally expressing the ESR1-e6>YAP1 fusion was also responsive. Transcriptionally active ESR1 fusions therefore trigger both endocrine therapy resistance and metastatic progression, explaining the association with fatal disease progression, although CDK4/6 inhibitor treatment is predicted to be effective.

Malignant rhabdoid tumors (MRTs) represent one of the most aggressive childhood malignancies. No effective treatment options are available, and prognosis is, therefore, dismal. Previous studies have demonstrated that tumor organoids capture the heterogeneity of patient tumors and can be used to predict patient response to therapy. Here, we perform drug screening on patient-derived normal and tumor organoids to identify MRT-specific therapeutic vulnerabilities. We identify neddylation inhibitor MLN4924 as a potential therapeutic agent. Mechanistically, we find increased neddylation in MRT organoids and tissues and show that MLN4924 induces a cytotoxic response via upregulation of the unfolded protein response. Lastly, we demonstrate in vivo efficacy in an MRT PDX mouse model, in which single-agent MLN4924 treatment significantly extends survival. Our study demonstrates that organoids can be used to find drugs selectively targeting tumor cells while leaving healthy cells unharmed and proposes neddylation inhibition as a therapeutic strategy in MRT.

Lung cancer remains the leading cause of cancer death. Genome sequencing of lung tumors from patients with squamous cell carcinoma has identified SMAD4 to be frequently mutated. Here, we use a mouse model to determine the molecular mechanisms by which Smad4 loss leads to lung cancer progression. Mice with ablation of Pten and Smad4 in airway epithelium develop metastatic adenosquamous tumors. Comparative transcriptomic and in vivo cistromic analyses determine that loss of PTEN and SMAD4 results in ELF3 and ErbB2 pathway activation due to decreased expression of ERRFI1, a negative regulator of ERBB2 in mouse and human cells. The combinatorial inhibition of ErbB2 and Akt signaling attenuate tumor progression and cell invasion, respectively. Expression profile analysis of human lung tumors substantiated the importance of the ErbB2/Akt/ELF3 signaling pathway as both a prognostic biomarker and a therapeutic drug target for treating lung cancer.

The ampulla of Vater is a complex cellular environment from which adenocarcinomas arise to form a group of histopathologically heterogenous tumors. To evaluate the molecular features of these tumors, 98 ampullary adenocarcinomas were evaluated and compared to 44 distal bile duct and 18 duodenal adenocarcinomas. Genomic analyses revealed mutations in the WNT signaling pathway among half of the patients and in all three adenocarcinomas irrespective of their origin and histological morphology. These tumors were characterized by a high frequency of inactivating mutations of ELF3, a high rate of microsatellite instability, and common focal deletions and amplifications, suggesting common attributes in the molecular pathogenesis are at play in these tumors. The high frequency of WNT pathway activating mutation, coupled with small-molecule inhibitors of β-catenin in clinical trials, suggests future treatment decisions for these patients may be guided by genomic analysis.

Cancer cells function as primary architects of the tumor microenvironment. However, the molecular features of cancer cells that govern stromal cell phenotypes remain unclear. Here, we show that cancer-associated fibroblast (CAF) heterogeneity is driven by lung adenocarcinoma (LUAD) cells at either end of the epithelial-to-mesenchymal transition (EMT) spectrum. LUAD cells that have high expression of the EMT-activating transcription factor ZEB1 reprogram CAFs through a ZEB1-dependent secretory program and direct CAFs to the tips of invasive projections through a ZEB1-driven CAF repulsion process. The EMT, in turn, sensitizes LUAD cells to pro-metastatic signals from CAFs. Thus, CAFs respond to contextual cues from LUAD cells to promote metastasis.

Most human cancers arise from stem and progenitor cells by the sequential accumulation of genetic and epigenetic alterations, while cancer modeling typically requires simultaneous multiple oncogenic events. Here, we show that a single p53 mutation, despite causing no defect in the mouse brain, promoted neural stem and progenitor cells to spontaneously accumulate oncogenic alterations, including loss of multiple chromosomal (chr) regions syntenic to human chr10 containing Pten, forming malignant gliomas with PI3K/Akt activation. Rictor/mTORC2 loss inhibited Akt signaling, greatly delaying and reducing glioma formation by suppressing glioma precursors within the subventricular zone stem cell niche. Rictor/mTORC2 loss delayed timely differentiation of granule cell precursors (GCPs) during cerebellar development, promoting sustained GCP proliferation and medulloblastoma formation, which recapitulated critical features of TP53 mutant sonic hedgehog (SHH) medulloblastomas with GLI2 and/or N-MYC amplification. Our study demonstrates that Rictor/mTORC2 has opposing functions in neural stem cells and GCPs in the adult and the developing brain, promoting malignant gliomas and suppressing SHH-medulloblastoma formation, respectively.

On the basis of multidimensional and comprehensive molecular characterization (including DNA methalylation and copy number, RNA, and protein expression), we classified 894 renal cell carcinomas (RCCs) of various histologic types into nine major genomic subtypes. Site of origin within the nephron was one major determinant in the classification, reflecting differences among clear cell, chromophobe, and papillary RCC. Widespread molecular changes associated with TFE3 gene fusion or chromatin modifier genes were present within a specific subtype and spanned multiple subtypes. Differences in patient survival and in alteration of specific pathways (including hypoxia, metabolism, MAP kinase, NRF2-ARE, Hippo, immune checkpoint, and PI3K/AKT/mTOR) could further distinguish the subtypes. Immune checkpoint markers and molecular signatures of T cell infiltrates were both highest in the subtype associated with aggressive clear cell RCC. Differences between the genomic subtypes suggest that therapeutic strategies could be tailored to each RCC disease subset.

Hippo signaling has been recognized as a key tumor suppressor pathway. Here, we perform a comprehensive molecular characterization of 19 Hippo core genes in 9,125 tumor samples across 33 cancer types using multidimensional "omic" data from The Cancer Genome Atlas. We identify somatic drivers among Hippo genes and the related microRNA (miRNA) regulators, and using functional genomic approaches, we experimentally characterize YAP and TAZ mutation effects and miR-590 and miR-200a regulation for TAZ. Hippo pathway activity is best characterized by a YAP/TAZ transcriptional target signature of 22 genes, which shows robust prognostic power across cancer types. Our elastic-net integrated modeling further reveals cancer-type-specific pathway regulators and associated cancer drivers. Our results highlight the importance of Hippo signaling in squamous cell cancers, characterized by frequent amplification of YAP/TAZ, high expression heterogeneity, and significant prognostic patterns. This study represents a systems-biology approach to characterizing key cancer signaling pathways in the post-genomic era.

The global impact of somatic structural variants (SVs) on gene regulation in advanced tumors with complex treatment histories has been mostly uncharacterized. Here, using whole-genome and RNA sequencing from 570 recurrent or metastatic tumors, we report the altered expression of hundreds of genes in association with nearby SV breakpoints, including oncogenes and G-protein-coupled receptor-related genes such as PLEKHG2. A significant fraction of genes with SV-expression associations correlate with worse patient survival in primary and advanced cancers, including SRD5A1. In many instances, SV-expression associations involve retrotransposons being translocated near genes. High overall SV burden is associated with treatment with DNA alkylating agents or taxanes and altered expression of metabolism-associated genes. SV-expression associations within tumors from topoisomerase I inhibitor-treated patients include chromatin-related genes. Within anthracycline-treated tumors, SV breakpoints near chromosome 1p genes include PDE4B. Patient treatment and history can help understand the widespread SV-mediated cis-regulatory alterations found in cancer.

Functional implication of stromal heterogeneity in the prostate remains incompletely understood. Using lineage tracing and light-sheet imaging, we show that some fibroblast cells at the mouse proximal prostatic ducts and prostatic urethra highly express Lgr5. Genetic ablation of these anatomically restricted stromal cells, but not nonselective ablation of prostatic stromal cells, rapidly induces prostate epithelial turnover and dedifferentiation that are reversed following spontaneous restoration of the Lgr5+ stromal cells. RNA sequencing (RNA-seq) analysis indicates that ablating the Lgr5+ stromal cells activates a mechanosensory response. Ablating the Lgr5+ stromal cells impairs the control of prostatic ductal outlet, increases prostate tissue stiffness, and activates the mitogen-activated protein kinase (MAPK). Suppressing MAPK overrides the elevated epithelial proliferation. In summary, the Lgr5+ stromal cells regulate prostate tissue homeostasis and maintain its functional integrity in a long-distance manner. Our study implies that the cells at organ junctions most likely control organ homeostasis by sustaining a balanced mechanoforce.