IL4I1 Is a Metabolic Immune Checkpoint that Activates the AHR and Promotes Tumor Progression.

Aryl hydrocarbon receptor (AHR) activation by tryptophan (Trp) catabolites enhances tumor malignancy and suppresses anti-tumor immunity. The context specificity of AHR target genes has so far impeded systematic investigation of AHR activity and its upstream enzymes across human cancers. A pan-tissue AHR signature, derived by natural language processing, revealed that across 32 tumor entities, interleukin-4-induced-1 (IL4I1) associates more frequently with AHR activity than IDO1 or TDO2, hitherto recognized as the main Trp-catabolic enzymes. IL4I1 activates the AHR through the generation of indole metabolites and kynurenic acid. It associates with reduced survival in glioma patients, promotes cancer cell motility, and suppresses adaptive immunity, thereby enhancing the progression of chronic lymphocytic leukemia (CLL) in mice. Immune checkpoint blockade (ICB) induces IDO1 and IL4I1. As IDO1 inhibitors do not block IL4I1, IL4I1 may explain the failure of clinical studies combining ICB with IDO1 inhibition. Taken together, IL4I1 blockade opens new avenues for cancer therapy.

Interleukin-10 receptor signaling promotes the maintenance of a PD-1int TCF-1+ CD8+ T cell population that sustains anti-tumor immunity.

T cell exhaustion limits anti-tumor immunity and responses to immunotherapy. Here, we explored the microenvironmental signals regulating T cell exhaustion using a model of chronic lymphocytic leukemia (CLL). Single-cell analyses identified a subset of PD-1hi, functionally impaired CD8+ T cells that accumulated in secondary lymphoid organs during disease progression and a functionally competent PD-1int subset. Frequencies of PD-1int TCF-1+ CD8+ T cells decreased upon Il10rb or Stat3 deletion, leading to accumulation of PD-1hi cells and accelerated tumor progression. Mechanistically, inhibition of IL-10R signaling altered chromatin accessibility and disrupted cooperativity between the transcription factors NFAT and AP-1, promoting a distinct NFAT-associated program. Low IL10 expression or loss of IL-10R-STAT3 signaling correlated with increased frequencies of exhausted CD8+ T cells and poor survival in CLL and in breast cancer patients. Thus, balance between PD-1hi, exhausted CD8+ T cells and functional PD-1int TCF-1+ CD8+ T cells is regulated by cell-intrinsic IL-10R signaling, with implications for immunotherapy.

Genome sequencing of pediatric medulloblastoma links catastrophic DNA rearrangements with TP53 mutations.

Genomic rearrangements are thought to occur progressively during tumor development. Recent findings, however, suggest an alternative mechanism, involving massive chromosome rearrangements in a one-step catastrophic event termed chromothripsis. We report the whole-genome sequencing-based analysis of a Sonic-Hedgehog medulloblastoma (SHH-MB) brain tumor from a patient with a germline TP53 mutation (Li-Fraumeni syndrome), uncovering massive, complex chromosome rearrangements. Integrating TP53 status with microarray and deep sequencing-based DNA rearrangement data in additional patients reveals a striking association between TP53 mutation and chromothripsis in SHH-MBs. Analysis of additional tumor entities substantiates a link between TP53 mutation and chromothripsis, and indicates a context-specific role for p53 in catastrophic DNA rearrangements. Among these, we observed a strong association between somatic TP53 mutations and chromothripsis in acute myeloid leukemia. These findings connect p53 status and chromothripsis in specific tumor types, providing a genetic basis for understanding particularly aggressive subtypes of cancer.

T-bet suppresses proliferation of malignant B cells in chronic lymphocytic leukemia.

The T-box transcription factor T-bet is known as a master regulator of T-cell response but its role in malignant B cells is not sufficiently explored. Here, we conducted single-cell resolved multi-omics analyses of malignant B cells from patients with chronic lymphocytic leukemia (CLL) and studied a CLL mouse model with genetic knockout of TBX21. We found that T-bet acts as a tumor suppressor in malignant B cells by decreasing their proliferation rate. NF-κB activity induced by inflammatory signals provided by the microenvironment, triggered T-bet expression which impacted on promoter proximal and distal chromatin co-accessibility and controlled a specific gene signature by mainly suppressing transcription. Gene set enrichment analysis identified a positive regulation of interferon signaling, and a negative control of proliferation by T-bet. In line, we showed that T-bet represses cell cycling and is associated with longer overall survival of CLL patients. Our study uncovers a novel tumor suppressive role of T-bet in malignant B cells via its regulation of inflammatory processes and cell cycling which has implications for stratification and therapy of CLL patients. Linking T-bet activity to inflammation explains the good prognostic role of genetic alterations in inflammatory signaling pathways in CLL.

Conumee 2.0: enhanced copy-number variation analysis from DNA methylation arrays for humans and mice.

MOTIVATION: Copy-number variations (CNVs) are common genetic alterations in cancer and their detection may impact tumor classification and therapeutic decisions. However, detection of clinically relevant large and focal CNVs remains challenging when sample material or resources are limited. This has motivated us to create a software tool to infer CNVs from DNA methylation arrays which are often generated as part of clinical routines and in research settings. RESULTS: We present our R package, conumee 2.0, that combines tangent normalization, an adjustable genomic binning heuristic, and weighted circular binary segmentation to utilize DNA methylation arrays for CNV analysis and mitigate technical biases and batch effects. Segmentation results were validated in a lung squamous cell carcinoma dataset from TCGA (n = 367 samples) by comparison to segmentations derived from genotyping arrays (Pearson's correlation coefficient of 0.91). We further introduce a segmented block bootstrapping approach to detect focal alternations that achieved 60.9% sensitivity and 98.6% specificity for deletions affecting CDKN2A/B (60.0% and 96.9% for RB1, respectively) in a low-grade glioma cohort from TCGA (n = 239 samples). Finally, our tool provides functionality to detect and summarize CNVs across large sample cohorts. AVAILABILITY AND IMPLEMENTATION: Conumee 2.0 is available under open-source license at: https://github.com/hovestadtlab/conumee2.

Therapeutic targeting of ependymoma as informed by oncogenic enhancer profiling.

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.

scSNPdemux: a sensitive demultiplexing pipeline using single nucleotide polymorphisms for improved pooled single-cell RNA sequencing analysis.

BACKGROUND: Here we present scSNPdemux, a sample demultiplexing pipeline for single-cell RNA sequencing data using natural genetic variations in humans. The pipeline requires alignment files from Cell Ranger (10× Genomics), a population SNP database and genotyped single nucleotide polymorphisms (SNPs) per sample. The tool works on sparse genotyping data in VCF format for sample identification. RESULTS: The pipeline was tested on both single-cell and single-nuclei based RNA sequencing datasets and showed superior demultiplexing performance over the lipid-based CellPlex and Multi-seq sample multiplexing technique which incurs additional single cell library preparation steps. Specifically, our pipeline demonstrated superior sensitivity and specificity in cell-identity assignment over CellPlex, especially on immune cell types with low RNA content. CONCLUSIONS: We designed a streamlined pipeline for single-cell sample demultiplexing, aiming to overcome common problems in multiplexing samples using single cell libraries which might affect data quality and can be costly.

The whole-genome landscape of medulloblastoma subtypes.

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.

A synergistic interaction between HDAC- and PARP inhibitors in childhood tumors with chromothripsis.

Chromothripsis is a form of genomic instability characterized by the occurrence of tens to hundreds of clustered DNA double-strand breaks in a one-off catastrophic event. Rearrangements associated with chromothripsis are detectable in numerous tumor entities and linked with poor prognosis in some of these, such as Sonic Hedgehog medulloblastoma, neuroblastoma and osteosarcoma. Hence, there is a need for therapeutic strategies eliminating tumor cells with chromothripsis. Defects in DNA double-strand break repair, and in particular homologous recombination repair, have been linked with chromothripsis. Targeting DNA repair deficiencies by synthetic lethality approaches, we performed a synergy screen using drug libraries (n = 375 compounds, 15 models) combined with either a PARP inhibitor or cisplatin. This revealed a synergistic interaction between the HDAC inhibitor romidepsin and PARP inhibition. Functional assays, transcriptome analyses and in vivo validation in patient-derived xenograft mouse models confirmed the efficacy of the combinatorial treatment.

Clonal evolution in chronic lymphocytic leukemia is scant in relapsed but accelerated in refractory cases after chemo(immune) therapy.

Clonal evolution is involved in the progression of chronic lymphocytic leukemia (CLL). In order to link evolutionary patterns to different disease courses, we performed a long-term longitudinal mutation profiling study of CLL patients. Tracking somatic mutations and their changes in allele frequency over time and assessing the underlying cancer cell fraction revealed highly distinct evolutionary patterns. Surprisingly, in long-term stable disease and in relapse after long-lasting clinical response to treatment, clonal shifts are minor. In contrast, in refractory disease major clonal shifts occur although there is little impact on leukemia cell counts. As this striking pattern in refractory cases is not linked to a strong contribution of known CLL driver genes, the evolution is mostly driven by treatment-induced selection of sub-clones, underlining the need for novel, non-genotoxic treatment regimens.

Active medulloblastoma enhancers reveal subgroup-specific cellular origins.

Medulloblastoma is a highly malignant paediatric brain tumour, often inflicting devastating consequences on the developing child. Genomic studies have revealed four distinct molecular subgroups with divergent biology and clinical behaviour. An understanding of the regulatory circuitry governing the transcriptional landscapes of medulloblastoma subgroups, and how this relates to their respective developmental origins, is lacking. Here, using H3K27ac and BRD4 chromatin immunoprecipitation followed by sequencing (ChIP-seq) coupled with tissue-matched DNA methylation and transcriptome data, we describe the active cis-regulatory landscape across 28 primary medulloblastoma specimens. Analysis of differentially regulated enhancers and super-enhancers reinforced inter-subgroup heterogeneity and revealed novel, clinically relevant insights into medulloblastoma biology. Computational reconstruction of core regulatory circuitry identified a master set of transcription factors, validated by ChIP-seq, that is responsible for subgroup divergence, and implicates candidate cells of origin for Group 4. Our integrated analysis of enhancer elements in a large series of primary tumour samples reveals insights into cis-regulatory architecture, unrecognized dependencies, and cellular origins.

A versatile system to introduce clusters of genomic double-strand breaks in large cell populations.

In vitro assays for clustered DNA lesions will facilitate the analysis of the mechanisms underlying complex genome rearrangements such as chromothripsis, including the recruitment of repair factors to sites of DNA double-strand breaks (DSBs). We present a novel method generating localized DNA DSBs using UV irradiation with photomasks. The size of the damage foci and the spacing between lesions are fully adjustable, making the assay suitable for different cell types and targeted areas. We validated this setup with genomically stable epithelial cells, normal fibroblasts, pluripotent stem cells, and patient-derived primary cultures. Our method does not require a specialized device such as a laser, making it accessible to a broad range of users. Sensitization by 5-bromo-2-deoxyuridine incorporation is not required, which enables analyzing the DNA damage response in post-mitotic cells. Irradiated cells can be cultivated further, followed by time-lapse imaging or used for downstream biochemical analyses, thanks to the high throughput of the system. Importantly, we showed genome rearrangements in the irradiated cells, providing a proof of principle for the induction of structural variants by localized DNA lesions.

Unraveling most abundant mutational signatures in head and neck cancer.

Genomic alterations are a driving force in the multistep process of head and neck cancer (HNC) and result from the interaction of exogenous environmental exposures and endogenous cellular processes. Each of these processes leaves a characteristic pattern of mutations on the tumor genome providing the unique opportunity to decipher specific signatures of mutational processes operative during HNC pathogenesis and to address their prognostic value. Computational analysis of whole exome sequencing data of the HIPO-HNC (Heidelberg Center for Personalized Oncology-head and neck cancer) (n = 83) and TCGA-HNSC (The Cancer Genome Atlas-Head and Neck Squamous Cell Carcinoma) (n = 506) cohorts revealed five common mutational signatures (Catalogue of Somatic Mutations in Cancer [COSMIC] Signatures 1, 2, 3, 13 and 16) and demonstrated their significant association with etiological risk factors (tobacco, alcohol and HPV16). Unsupervised hierarchical clustering identified four clusters (A, B, C1 and C2) of which Subcluster C2 was enriched for cases with a higher frequency of signature 16 mutations. Tumors of Subcluster C2 had significantly lower p16INK4A expression accompanied by homozygous CDKN2A deletion in almost one half of cases. Survival analysis revealed an unfavorable prognosis for patients with tumors characterized by a higher mutation burden attributed to signature 16 as well as cases in Subcluster C2. Finally, a LASSO-Cox regression model was applied to prioritize clinically relevant signatures and to establish a prognostic risk score for head and neck squamous cell carcinoma patients. In conclusion, our study provides a proof of concept that computational analysis of somatic mutational signatures is not only a powerful tool to decipher environmental and intrinsic processes in the pathogenesis of HNC, but could also pave the way to establish reliable prognostic patterns.

Linking aberrant chromatin features in chronic lymphocytic leukemia to transcription factor networks.

In chronic lymphocytic leukemia (CLL), a diverse set of genetic mutations is embedded in a deregulated epigenetic landscape that drives cancerogenesis. To elucidate the role of aberrant chromatin features, we mapped DNA methylation, seven histone modifications, nucleosome positions, chromatin accessibility, binding of EBF1 and CTCF, as well as the transcriptome of B cells from CLL patients and healthy donors. A globally increased histone deacetylase activity was detected and half of the genome comprised transcriptionally downregulated partially DNA methylated domains demarcated by CTCF CLL samples displayed a H3K4me3 redistribution and nucleosome gain at promoters as well as changes of enhancer activity and enhancer linkage to target genes. A DNA binding motif analysis identified transcription factors that gained or lost binding in CLL at sites with aberrant chromatin features. These findings were integrated into a gene regulatory enhancer containing network enriched for B-cell receptor signaling pathway components. Our study predicts novel molecular links to targets of CLL therapies and provides a valuable resource for further studies on the epigenetic contribution to the disease.

Molecular subgrouping of primary pineal parenchymal tumors reveals distinct subtypes correlated with clinical parameters and genetic alterations.

Tumors of the pineal region comprise several different entities with distinct clinical and histopathological features. Whereas some entities predominantly affect adults, pineoblastoma (PB) constitutes a highly aggressive malignancy of childhood with a poor outcome. PBs mainly arise sporadically, but may also occur in the context of cancer predisposition syndromes including DICER1 and RB1 germline mutation. With this study, we investigate clinico-pathological subgroups of pineal tumors and further characterize their biological features. We performed genome-wide DNA methylation analysis in 195 tumors of the pineal region and 20 normal pineal gland controls. Copy-number profiles were obtained from DNA methylation data; gene panel sequencing was added for 93 tumors and analysis was further complemented by miRNA sequencing for 22 tumor samples. Unsupervised clustering based on DNA methylation profiling separated known subgroups, like pineocytoma, pineal parenchymal tumor of intermediate differentiation, papillary tumor of the pineal region and PB, and further distinct subtypes within these groups, including three subtypes within the core PB subgroup. The novel molecular subgroup Pin-RB includes cases of trilateral retinoblastoma as well as sporadic pineal tumors with RB1 alterations, and displays similarities with retinoblastoma. Distinct clinical associations discriminate the second novel molecular subgroup PB-MYC from other PB cases. Alterations within the miRNA processing pathway (affecting DROSHA, DGCR8 or DICER1) are found in about two thirds of cases in the three core PB subtypes. Methylation profiling revealed biologically distinct groups of pineal tumors with specific clinical and molecular features. Our findings provide a foundation for further clinical as well as molecular and functional characterization of PB and other pineal tumors, including the role of miRNA processing defects in oncogenesis.

Enhancer hijacking activates GFI1 family oncogenes in medulloblastoma.

Medulloblastoma is a highly malignant paediatric brain tumour currently treated with a combination of surgery, radiation and chemotherapy, posing a considerable burden of toxicity to the developing child. Genomics has illuminated the extensive intertumoral heterogeneity of medulloblastoma, identifying four distinct molecular subgroups. Group 3 and group 4 subgroup medulloblastomas account for most paediatric cases; yet, oncogenic drivers for these subtypes remain largely unidentified. Here we describe a series of prevalent, highly disparate genomic structural variants, restricted to groups 3 and 4, resulting in specific and mutually exclusive activation of the growth factor independent 1 family proto-oncogenes, GFI1 and GFI1B. Somatic structural variants juxtapose GFI1 or GFI1B coding sequences proximal to active enhancer elements, including super-enhancers, instigating oncogenic activity. Our results, supported by evidence from mouse models, identify GFI1 and GFI1B as prominent medulloblastoma oncogenes and implicate 'enhancer hijacking' as an efficient mechanism driving oncogene activation in a childhood cancer.

Decoding the regulatory landscape of medulloblastoma using DNA methylation sequencing.

Epigenetic alterations, that is, disruption of DNA methylation and chromatin architecture, are now acknowledged as a universal feature of tumorigenesis. Medulloblastoma, a clinically challenging, malignant childhood brain tumour, is no exception. Despite much progress from recent genomics studies, with recurrent changes identified in each of the four distinct tumour subgroups (WNT-pathway-activated, SHH-pathway-activated, and the less-well-characterized Group 3 and Group 4), many cases still lack an obvious genetic driver. Here we present whole-genome bisulphite-sequencing data from thirty-four human and five murine tumours plus eight human and three murine normal controls, augmented with matched whole-genome, RNA and chromatin immunoprecipitation sequencing data. This comprehensive data set allowed us to decipher several features underlying the interplay between the genome, epigenome and transcriptome, and its effects on medulloblastoma pathophysiology. Most notable were highly prevalent regions of hypomethylation correlating with increased gene expression, extending tens of kilobases downstream of transcription start sites. Focal regions of low methylation linked to transcription-factor-binding sites shed light on differential transcriptional networks between subgroups, whereas increased methylation due to re-normalization of repressed chromatin in DNA methylation valleys was positively correlated with gene expression. Large, partially methylated domains affecting up to one-third of the genome showed increased mutation rates and gene silencing in a subgroup-specific fashion. Epigenetic alterations also affected novel medulloblastoma candidate genes (for example, LIN28B), resulting in alternative promoter usage and/or differential messenger RNA/microRNA expression. Analysis of mouse medulloblastoma and precursor-cell methylation demonstrated a somatic origin for many alterations. Our data provide insights into the epigenetic regulation of transcription and genome organization in medulloblastoma pathogenesis, which are probably also of importance in a wider developmental and disease context.

Analysis of Epstein-Barr Virus Genomes and Expression Profiles in Gastric Adenocarcinoma.

Epstein-Barr virus (EBV) is a causative agent of a variety of lymphomas, nasopharyngeal carcinoma (NPC), and ∼9% of gastric carcinomas (GCs). An important question is whether particular EBV variants are more oncogenic than others, but conclusions are currently hampered by the lack of sequenced EBV genomes. Here, we contribute to this question by mining whole-genome sequences of 201 GCs to identify 13 EBV-positive GCs and by assembling 13 new EBV genome sequences, almost doubling the number of available GC-derived EBV genome sequences and providing the first non-Asian EBV genome sequences from GC. Whole-genome sequence comparisons of all EBV isolates sequenced to date (85 from tumors and 57 from healthy individuals) showed that most GC and NPC EBV isolates were closely related although American Caucasian GC samples were more distant, suggesting a geographical component. However, EBV GC isolates were found to contain some consistent changes in protein sequences regardless of geographical origin. In addition, transcriptome data available for eight of the EBV-positive GCs were analyzed to determine which EBV genes are expressed in GC. In addition to the expected latency proteins (EBNA1, LMP1, and LMP2A), specific subsets of lytic genes were consistently expressed that did not reflect a typical lytic or abortive lytic infection, suggesting a novel mechanism of EBV gene regulation in the context of GC. These results are consistent with a model in which a combination of specific latent and lytic EBV proteins promotes tumorigenesis.IMPORTANCE Epstein-Barr virus (EBV) is a widespread virus that causes cancer, including gastric carcinoma (GC), in a small subset of individuals. An important question is whether particular EBV variants are more cancer associated than others, but more EBV sequences are required to address this question. Here, we have generated 13 new EBV genome sequences from GC, almost doubling the number of EBV sequences from GC isolates and providing the first EBV sequences from non-Asian GC. We further identify sequence changes in some EBV proteins common to GC isolates. In addition, gene expression analysis of eight of the EBV-positive GCs showed consistent expression of both the expected latency proteins and a subset of lytic proteins that was not consistent with typical lytic or abortive lytic expression. These results suggest that novel mechanisms activate expression of some EBV lytic proteins and that their expression may contribute to oncogenesis.

Driver mutations in histone H3.3 and chromatin remodelling genes in paediatric glioblastoma.

Glioblastoma multiforme (GBM) is a lethal brain tumour in adults and children. However, DNA copy number and gene expression signatures indicate differences between adult and paediatric cases. To explore the genetic events underlying this distinction, we sequenced the exomes of 48 paediatric GBM samples. Somatic mutations in the H3.3-ATRX-DAXX chromatin remodelling pathway were identified in 44% of tumours (21/48). Recurrent mutations in H3F3A, which encodes the replication-independent histone 3 variant H3.3, were observed in 31% of tumours, and led to amino acid substitutions at two critical positions within the histone tail (K27M, G34R/G34V) involved in key regulatory post-translational modifications. Mutations in ATRX (α-thalassaemia/mental retardation syndrome X-linked) and DAXX (death-domain associated protein), encoding two subunits of a chromatin remodelling complex required for H3.3 incorporation at pericentric heterochromatin and telomeres, were identified in 31% of samples overall, and in 100% of tumours harbouring a G34R or G34V H3.3 mutation. Somatic TP53 mutations were identified in 54% of all cases, and in 86% of samples with H3F3A and/or ATRX mutations. Screening of a large cohort of gliomas of various grades and histologies (n = 784) showed H3F3A mutations to be specific to GBM and highly prevalent in children and young adults. Furthermore, the presence of H3F3A/ATRX-DAXX/TP53 mutations was strongly associated with alternative lengthening of telomeres and specific gene expression profiles. This is, to our knowledge, the first report to highlight recurrent mutations in a regulatory histone in humans, and our data suggest that defects of the chromatin architecture underlie paediatric and young adult GBM pathogenesis.