Accurate pathological diagnosis is crucial for optimal management of patients with cancer. For the approximately 100 known tumour types of the central nervous system, standardization of the diagnostic process has been shown to be particularly challenging-with substantial inter-observer variability in the histopathological diagnosis of many tumour types. Here we present a comprehensive approach for the DNA methylation-based classification of central nervous system tumours across all entities and age groups, and demonstrate its application in a routine diagnostic setting. We show that the availability of this method may have a substantial impact on diagnostic precision compared to standard methods, resulting in a change of diagnosis in up to 12% of prospective cases. For broader accessibility, we have designed a free online classifier tool, the use of which does not require any additional onsite data processing. Our results provide a blueprint for the generation of machine-learning-based tumour classifiers across other cancer entities, with the potential to fundamentally transform tumour pathology.

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

Recent exon-sequencing studies of human tumours have revealed that subunits of BAF (mammalian SWI/SNF) complexes are mutated in more than 20% of all human malignancies, but the mechanisms involved in tumour suppression are unclear. BAF chromatin-remodelling complexes are polymorphic assemblies that use energy provided by ATP hydrolysis to regulate transcription through the control of chromatin structure and the placement of Polycomb repressive complex 2 (PRC2) across the genome. Several proteins dedicated to this multisubunit complex, including BRG1 (also known as SMARCA4) and BAF250a (also known as ARID1A), are mutated at frequencies similar to those of recognized tumour suppressors. In particular, the core ATPase BRG1 is mutated in 5-10% of childhood medulloblastomas and more than 15% of Burkitt's lymphomas. Here we show a previously unknown function of BAF complexes in decatenating newly replicated sister chromatids, a requirement for proper chromosome segregation during mitosis. We find that deletion of Brg1 in mouse cells, as well as the expression of BRG1 point mutants identified in human tumours, leads to anaphase bridge formation (in which sister chromatids are linked by catenated strands of DNA) and a G2/M-phase block characteristic of the decatenation checkpoint. Endogenous BAF complexes interact directly with endogenous topoisomerase IIα (TOP2A) through BAF250a and are required for the binding of TOP2A to approximately 12,000 sites across the genome. Our results demonstrate that TOP2A chromatin binding is dependent on the ATPase activity of BRG1, which is compromised in oncogenic BRG1 mutants. These studies indicate that the ability of TOP2A to prevent DNA entanglement at mitosis requires BAF complexes and suggest that this activity contributes to the role of BAF subunits as tumour suppressors.

Genomic structural variants (SVs) are abundant in humans, differing from other forms of variation in extent, origin and functional impact. Despite progress in SV characterization, the nucleotide resolution architecture of most SVs remains unknown. We constructed a map of unbalanced SVs (that is, copy number variants) based on whole genome DNA sequencing data from 185 human genomes, integrating evidence from complementary SV discovery approaches with extensive experimental validations. Our map encompassed 22,025 deletions and 6,000 additional SVs, including insertions and tandem duplications. Most SVs (53%) were mapped to nucleotide resolution, which facilitated analysing their origin and functional impact. We examined numerous whole and partial gene deletions with a genotyping approach and observed a depletion of gene disruptions amongst high frequency deletions. Furthermore, we observed differences in the size spectra of SVs originating from distinct formation mechanisms, and constructed a map of SV hotspots formed by common mechanisms. Our analytical framework and SV map serves as a resource for sequencing-based association studies.

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and 'enhancer hijacking' events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Genomic sequencing has driven precision-based oncology therapy; however, the genetic drivers of many malignancies remain unknown or non-targetable, so alternative approaches to the identification of therapeutic leads are necessary. Ependymomas are chemotherapy-resistant brain tumours, which, despite genomic sequencing, lack effective molecular targets. Intracranial ependymomas are segregated on the basis of anatomical location (supratentorial region or posterior fossa) and further divided into distinct molecular subgroups that reflect differences in the age of onset, gender predominance and response to therapy. The most common and aggressive subgroup, posterior fossa ependymoma group A (PF-EPN-A), occurs in young children and appears to lack recurrent somatic mutations. Conversely, posterior fossa ependymoma group B (PF-EPN-B) tumours display frequent large-scale copy number gains and losses but have favourable clinical outcomes. More than 70% of supratentorial ependymomas are defined by highly recurrent gene fusions in the NF-κB subunit gene RELA (ST-EPN-RELA), and a smaller number involve fusion of the gene encoding the transcriptional activator YAP1 (ST-EPN-YAP1). Subependymomas, a distinct histologic variant, can also be found within the supratetorial and posterior fossa compartments, and account for the majority of tumours in the molecular subgroups ST-EPN-SE and PF-EPN-SE. Here we describe mapping of active chromatin landscapes in 42 primary ependymomas in two non-overlapping primary ependymoma cohorts, with the goal of identifying essential super-enhancer-associated genes on which tumour cells depend. Enhancer regions revealed putative oncogenes, molecular targets and pathways; inhibition of these targets with small molecule inhibitors or short hairpin RNA diminished the proliferation of patient-derived neurospheres and increased survival in mouse models of ependymomas. Through profiling of transcriptional enhancers, our study provides a framework for target and drug discovery in other cancers that lack known genetic drivers and are therefore difficult to treat.

Embryonal tumours with multilayered rosettes (ETMRs) are aggressive paediatric embryonal brain tumours with a universally poor prognosis1. Here we collected 193 primary ETMRs and 23 matched relapse samples to investigate the genomic landscape of this distinct tumour type. We found that patients with tumours in which the proposed driver C19MC2-4 was not amplified frequently had germline mutations in DICER1 or other microRNA-related aberrations such as somatic amplification of miR-17-92 (also known as MIR17HG). Whole-genome sequencing revealed that tumours had an overall low recurrence of single-nucleotide variants (SNVs), but showed prevalent genomic instability caused by widespread occurrence of R-loop structures. We show that R-loop-associated chromosomal instability can be induced by the loss of DICER1 function. Comparison of primary tumours and matched relapse samples showed a strong conservation of structural variants, but low conservation of SNVs. Moreover, many newly acquired SNVs are associated with a mutational signature related to cisplatin treatment. Finally, we show that targeting R-loops with topoisomerase and PARP inhibitors might be an effective treatment strategy for this deadly disease.

Mutated isocitrate dehydrogenase 1 (IDH1) defines a molecularly distinct subtype of diffuse glioma1-3. The most common IDH1 mutation in gliomas affects codon 132 and encodes IDH1(R132H), which harbours a shared clonal neoepitope that is presented on major histocompatibility complex (MHC) class II4,5. An IDH1(R132H)-specific peptide vaccine (IDH1-vac) induces specific therapeutic T helper cell responses that are effective against IDH1(R132H)+ tumours in syngeneic MHC-humanized mice4,6-8. Here we describe a multicentre, single-arm, open-label, first-in-humans phase I trial that we carried out in 33 patients with newly diagnosed World Health Organization grade 3 and 4 IDH1(R132H)+ astrocytomas (Neurooncology Working Group of the German Cancer Society trial 16 (NOA16), ClinicalTrials.gov identifier NCT02454634). The trial met its primary safety endpoint, with vaccine-related adverse events restricted to grade 1. Vaccine-induced immune responses were observed in 93.3% of patients across multiple MHC alleles. Three-year progression-free and death-free rates were 0.63 and 0.84, respectively. Patients with immune responses showed a two-year progression-free rate of 0.82. Two patients without an immune response showed tumour progression within two years of first diagnosis. A mutation-specificity score that incorporates the duration and level of vaccine-induced IDH1(R132H)-specific T cell responses was associated with intratumoral presentation of the IDH1(R132H) neoantigen in pre-treatment tumour tissue. There was a high frequency of pseudoprogression, which indicates intratumoral inflammatory reactions. Pseudoprogression was associated with increased vaccine-induced peripheral T cell responses. Combined single-cell RNA and T cell receptor sequencing showed that tumour-infiltrating CD40LG+ and CXCL13+ T helper cell clusters in a patient with pseudoprogression were dominated by a single IDH1(R132H)-reactive T cell receptor.

The International Cancer Genome Consortium (ICGC) was launched to coordinate large-scale cancer genome studies in tumours from 50 different cancer types and/or subtypes that are of clinical and societal importance across the globe. Systematic studies of more than 25,000 cancer genomes at the genomic, epigenomic and transcriptomic levels will reveal the repertoire of oncogenic mutations, uncover traces of the mutagenic influences, define clinically relevant subtypes for prognosis and therapeutic management, and enable the development of new cancer therapies.

Medulloblastoma, the most common malignant paediatric brain tumour, arises in the cerebellum and disseminates through the cerebrospinal fluid in the leptomeningeal space to coat the brain and spinal cord. Dissemination, a marker of poor prognosis, is found in up to 40% of children at diagnosis and in most children at the time of recurrence. Affected children therefore are treated with radiation to the entire developing brain and spinal cord, followed by high-dose chemotherapy, with the ensuing deleterious effects on the developing nervous system. The mechanisms of dissemination through the cerebrospinal fluid are poorly studied, and medulloblastoma metastases have been assumed to be biologically similar to the primary tumour. Here we show that in both mouse and human medulloblastoma, the metastases from an individual are extremely similar to each other but are divergent from the matched primary tumour. Clonal genetic events in the metastases can be demonstrated in a restricted subclone of the primary tumour, suggesting that only rare cells within the primary tumour have the ability to metastasize. Failure to account for the bicompartmental nature of metastatic medulloblastoma could be a major barrier to the development of effective targeted therapies.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals1. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample.