Sarcomas are malignant soft tissue and bone tumours affecting adults, adolescents and children. They represent a morphologically heterogeneous class of tumours and some entities lack defining histopathological features. Therefore, the diagnosis of sarcomas is burdened with a high inter-observer variability and misclassification rate. Here, we demonstrate classification of soft tissue and bone tumours using a machine learning classifier algorithm based on array-generated DNA methylation data. This sarcoma classifier is trained using a dataset of 1077 methylation profiles from comprehensively pre-characterized cases comprising 62 tumour methylation classes constituting a broad range of soft tissue and bone sarcoma subtypes across the entire age spectrum. The performance is validated in a cohort of 428 sarcomatous tumours, of which 322 cases were classified by the sarcoma classifier. Our results demonstrate the potential of the DNA methylation-based sarcoma classification for research and future diagnostic applications.

Medulloblastoma is the most common malignant brain tumor of childhood. Group 3 medulloblastoma, the most aggressive molecular subtype, frequently disseminates through the leptomeningeal cerebral spinal fluid (CSF) spaces in the brain and spinal cord. The mechanism of dissemination through the CSF remains poorly understood, and the molecular pathways involved in medulloblastoma metastasis and self-renewal are largely unknown. Here we show that NOTCH1 signaling pathway regulates both the initiation of metastasis and the self-renewal of medulloblastoma. We identify a mechanism in which NOTCH1 activates BMI1 through the activation of TWIST1. NOTCH1 expression and activity are directly related to medulloblastoma metastasis and decreased survival rate of tumor-bearing mice. Finally, medulloblastoma-bearing mice intrathecally treated with anti-NRR1, a NOTCH1 blocking antibody, present lower frequency of spinal metastasis and higher survival rate. These findings identify NOTCH1 as a pivotal driver of Group 3 medulloblastoma metastasis and self-renewal, supporting the development of therapies targeting this pathway.

How the genetic landscape governs a tumor's response to immunotherapy remains poorly understood. To assess the immune-modulatory capabilities of 573 genes associated with altered cytotoxicity in human cancers, here we perform CRISPR/Cas9 screens directly in mouse lung cancer models. We recover the known immune evasion factors Stat1 and Serpinb9 and identify the cancer testis antigen Adam2 as an immune modulator, whose expression is induced by KrasG12D and further elevated by immunotherapy. Using loss- and gain-of-function experiments, we show that ADAM2 functions as an oncogene by restraining interferon and TNF cytokine signaling causing reduced presentation of tumor-associated antigens. ADAM2 also restricts expression of the immune checkpoint inhibitors PDL1, LAG3, TIGIT and TIM3 in the tumor microenvironment, which might explain why ex vivo expanded and adoptively transferred cytotoxic T-cells show enhanced cytotoxic efficacy in ADAM2 overexpressing tumors. Together, direct in vivo CRISPR/Cas9 screens can uncover genetic alterations that control responses to immunotherapies.

We present a novel hierarchical Bayes statistical model, xseq, to systematically quantify the impact of somatic mutations on expression profiles. We establish the theoretical framework and robust inference characteristics of the method using computational benchmarking. We then use xseq to analyse thousands of tumour data sets available through The Cancer Genome Atlas, to systematically quantify somatic mutations impacting expression profiles. We identify 30 novel cis-effect tumour suppressor gene candidates, enriched in loss-of-function mutations and biallelic inactivation. Analysis of trans-effects of mutations and copy number alterations with xseq identifies mutations in 150 genes impacting expression networks, with 89 novel predictions. We reveal two important novel characteristics of mutation impact on expression: (1) patients harbouring known driver mutations exhibit different downstream gene expression consequences; (2) expression patterns for some mutations are stable across tumour types. These results have critical implications for identification and interpretation of mutations with consequent impact on transcription in cancer.

Multi-omics datasets represent distinct aspects of the central dogma of molecular biology. Such high-dimensional molecular profiles pose challenges to data interpretation and hypothesis generation. ActivePathways is an integrative method that discovers significantly enriched pathways across multiple datasets using statistical data fusion, rationalizes contributing evidence and highlights associated genes. As part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancers across 38 tumor types, we integrated genes with coding and non-coding mutations and revealed frequently mutated pathways and additional cancer genes with infrequent mutations. We also analyzed prognostic molecular pathways by integrating genomic and transcriptomic features of 1780 breast cancers and highlighted associations with immune response and anti-apoptotic signaling. Integration of ChIP-seq and RNA-seq data for master regulators of the Hippo pathway across normal human tissues identified processes of tissue regeneration and stem cell regulation. ActivePathways is a versatile method that improves systems-level understanding of cellular organization in health and disease through integration of multiple molecular datasets and pathway annotations.

The most common events in breast cancer (BC) involve chromosome arm losses and gains. Here we describe identification of 1089 gene-centric common insertion sites (gCIS) from transposon-based screens in 8 mouse models of BC. Some gCIS are driver-specific, others driver non-specific, and still others associated with tumor histology. Processes affected by driver-specific and histology-specific mutations include well-known cancer pathways. Driver non-specific gCIS target the Mediator complex, Ca++ signaling, Cyclin D turnover, RNA-metabolism among other processes. Most gCIS show single allele disruption and many map to genomic regions showing high-frequency hemizygous loss in human BC. Two gCIS, Nf1 and Trps1, show synthetic haploinsufficient tumor suppressor activity. Many gCIS act on the same pathway responsible for tumor initiation, thereby selecting and sculpting just enough and just right signaling. These data highlight ~1000 genes with predicted conditional haploinsufficient tumor suppressor function and the potential to promote chromosome arm loss in BC.

Sonic hedgehog medulloblastoma encompasses a clinically and molecularly diverse group of cancers of the developing central nervous system. Here, we use unbiased sequencing of the transcriptome across a large cohort of 250 tumors to reveal differences among molecular subtypes of the disease, and demonstrate the previously unappreciated importance of non-coding RNA transcripts. We identify alterations within the cAMP dependent pathway (GNAS, PRKAR1A) which converge on GLI2 activity and show that 18% of tumors have a genetic event that directly targets the abundance and/or stability of MYCN. Furthermore, we discover an extensive network of fusions in focally amplified regions encompassing GLI2, and several loss-of-function fusions in tumor suppressor genes PTCH1, SUFU and NCOR1. Molecular convergence on a subset of genes by nucleotide variants, copy number aberrations, and gene fusions highlight the key roles of specific pathways in the pathogenesis of Sonic hedgehog medulloblastoma and open up opportunities for therapeutic intervention.

The catalog of cancer driver mutations in protein-coding genes has greatly expanded in the past decade. However, non-coding cancer driver mutations are less well-characterized and only a handful of recurrent non-coding mutations, most notably TERT promoter mutations, have been reported. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancer across 38 tumor types, we perform multi-faceted pathway and network analyses of non-coding mutations across 2583 whole cancer genomes from 27 tumor types compiled by the ICGC/TCGA PCAWG project that was motivated by the success of pathway and network analyses in prioritizing rare mutations in protein-coding genes. While few non-coding genomic elements are recurrently mutated in this cohort, we identify 93 genes harboring non-coding mutations that cluster into several modules of interacting proteins. Among these are promoter mutations associated with reduced mRNA expression in TP53, TLE4, and TCF4. We find that biological processes had variable proportions of coding and non-coding mutations, with chromatin remodeling and proliferation pathways altered primarily by coding mutations, while developmental pathways, including Wnt and Notch, altered by both coding and non-coding mutations. RNA splicing is primarily altered by non-coding mutations in this cohort, and samples containing non-coding mutations in well-known RNA splicing factors exhibit similar gene expression signatures as samples with coding mutations in these genes. These analyses contribute a new repertoire of possible cancer genes and mechanisms that are altered by non-coding mutations and offer insights into additional cancer vulnerabilities that can be investigated for potential therapeutic treatments.

Pediatric medulloblastoma (MB) is the most common solid malignant brain neoplasm, with Group 3 (G3) MB representing the most aggressive subgroup. MYC amplification is an independent poor prognostic factor in G3 MB, however, therapeutic targeting of the MYC pathway remains limited and alternative therapies for G3 MB are urgently needed. Here we show that the RNA-binding protein, Musashi-1 (MSI1) is an essential mediator of G3 MB in both MYC-overexpressing mouse models and patient-derived xenografts. MSI1 inhibition abrogates tumor initiation and significantly prolongs survival in both models. We identify binding targets of MSI1 in normal neural and G3 MB stem cells and then cross referenced these data with unbiased large-scale screens at the transcriptomic, translatomic and proteomic levels to systematically dissect its functional role. Comparative integrative multi-omic analyses of these large datasets reveal cancer-selective MSI1-bound targets sharing multiple MYC associated pathways, providing a valuable resource for context-specific therapeutic targeting of G3 MB.

Pediatric high-grade gliomas of the subclass MYCN (HGG-MYCN) are highly aggressive tumors frequently carrying MYCN amplifications, TP53 mutations, or both alterations. Due to their rarity, such tumors have only recently been identified as a distinct entity, and biological as well as clinical characteristics have not been addressed specifically. To gain insights into tumorigenesis and molecular profiles of these tumors, and to ultimately suggest alternative treatment options, we generated a genetically engineered mouse model by breeding hGFAP-cre::Trp53Fl/Fl::lsl-MYCN mice. All mice developed aggressive forebrain tumors early in their lifetime that mimic human HGG-MYCN regarding histology, DNA methylation, and gene expression. Single-cell RNA sequencing revealed a high intratumoral heterogeneity with neuronal and oligodendroglial lineage signatures. High-throughput drug screening using both mouse and human tumor cells finally indicated high efficacy of Doxorubicin, Irinotecan, and Etoposide as possible therapy options that children with HGG-MYCN might benefit from.