Diffuse midline gliomas and posterior fossa type A ependymomas contain the recurrent histone H3 lysine 27 (H3 K27M) mutation and express the H3 K27M-mimic EZHIP (CXorf67), respectively. H3 K27M and EZHIP are competitive inhibitors of Polycomb Repressive Complex 2 (PRC2) lysine methyltransferase activity. In vivo, these proteins reduce overall H3 lysine 27 trimethylation (H3K27me3) levels; however, residual peaks of H3K27me3 remain at CpG islands (CGIs) through an unknown mechanism. Here, we report that EZHIP and H3 K27M preferentially interact with PRC2 that is allosterically activated by H3K27me3 at CGIs and impede its spreading. Moreover, H3 K27M oncohistones reduce H3K27me3 in trans, independent of their incorporation into the chromatin. Although EZHIP is not found outside placental mammals, expression of human EZHIP reduces H3K27me3 in Drosophila melanogaster through a conserved mechanism. Our results provide mechanistic insights for the retention of residual H3K27me3 in tumors driven by H3 K27M and EZHIP.

Myxopapillary ependymoma (MPE) is a heterogeneous disease regarding histopathology and outcome. The underlying molecular biology is poorly understood, and markers that reliably predict the patients' clinical course are unknown.

Ependymomas encompass a heterogeneous group of central nervous system (CNS) neoplasms that occur along the entire neuroaxis. In recent years, extensive (epi-)genomic profiling efforts have identified several molecular groups of ependymoma that are characterized by distinct molecular alterations and/or patterns. Based on unsupervised visualization of a large cohort of genome-wide DNA methylation data, we identified a highly distinct group of pediatric-type tumors (n = 40) forming a cluster separate from all established CNS tumor types, of which a high proportion were histopathologically diagnosed as ependymoma. RNA sequencing revealed recurrent fusions involving the pleomorphic adenoma gene-like 1 (PLAGL1) gene in 19 of 20 of the samples analyzed, with the most common fusion being EWSR1:PLAGL1 (n = 13). Five tumors showed a PLAGL1:FOXO1 fusion and one a PLAGL1:EP300 fusion. High transcript levels of PLAGL1 were noted in these tumors, with concurrent overexpression of the imprinted genes H19 and IGF2, which are regulated by PLAGL1. Histopathological review of cases with sufficient material (n = 16) demonstrated a broad morphological spectrum of tumors with predominant ependymoma-like features. Immunohistochemically, tumors were GFAP positive and OLIG2- and SOX10 negative. In 3/16 of the cases, a dot-like positivity for EMA was detected. All tumors in our series were located in the supratentorial compartment. Median age of the patients at the time of diagnosis was 6.2 years. Median progression-free survival was 35 months (for 11 patients with data available). In summary, our findings suggest the existence of a novel group of supratentorial neuroepithelial tumors that are characterized by recurrent PLAGL1 fusions and enriched for pediatric patients.

Large-scale molecular profiling studies in recent years have shown that central nervous system (CNS) tumors display a much greater heterogeneity in terms of molecularly distinct entities, cellular origins and genetic drivers than anticipated from histological assessment. DNA methylation profiling has emerged as a useful tool for robust tumor classification, providing new insights into these heterogeneous molecular classes. This is particularly true for rare CNS tumors with a broad morphological spectrum, which are not possible to assign as separate entities based on histological similarity alone. Here, we describe a molecularly distinct subset of predominantly pediatric CNS neoplasms (n = 60) that harbor PATZ1 fusions. The original histological diagnoses of these tumors covered a wide spectrum of tumor types and malignancy grades. While the single most common diagnosis was glioblastoma (GBM), clinical data of the PATZ1-fused tumors showed a better prognosis than typical GBM, despite frequent relapses. RNA sequencing revealed recurrent MN1:PATZ1 or EWSR1:PATZ1 fusions related to (often extensive) copy number variations on chromosome 22, where PATZ1 and the two fusion partners are located. These fusions have individually been reported in a number of glial/glioneuronal tumors, as well as extracranial sarcomas. We show here that they are more common than previously acknowledged, and together define a biologically distinct CNS tumor type with high expression of neural development markers such as PAX2, GATA2 and IGF2. Drug screening performed on the MN1:PATZ1 fusion-bearing KS-1 brain tumor cell line revealed preliminary candidates for further study. In summary, PATZ1 fusions define a molecular class of histologically polyphenotypic neuroepithelial tumors, which show an intermediate prognosis under current treatment regimens.

Methylation of histone 3 lysine 36 (H3K36me) has emerged as an essential epigenetic component for the faithful regulation of gene expression. Despite its importance in development, disease, and cancer, how the molecular agents collectively shape the H3K36me landscape is unclear.

As the progression of low-grade diffuse astrocytomas into grade 4 tumors significantly impacts patient prognosis, a better understanding of this process is of paramount importance for improved patient care. In this project, we analyzed matched IDH-mutant astrocytomas before and after progression to grade 4 from six patients (discovery cohort) with genome-wide sequencing, 21 additional patients with targeted sequencing, and 33 patients from Glioma Longitudinal AnalySiS cohort for validation. The Cancer Genome Atlas data from 595 diffuse gliomas provided supportive information. All patients in our discovery cohort received radiation, all but one underwent chemotherapy, and no patient received temozolomide (TMZ) before progression to grade 4 disease. One case in the discovery cohort exhibited a hypermutation signature associated with the inactivation of the MSH2 and DNMT3A genes. In other patients, the number of chromosomal rearrangements and deletions increased in grade 4 tumors. The cell cycle checkpoint gene CDKN2A, or less frequently RB1, was most commonly inactivated after receiving both chemo- and radiotherapy when compared to other treatment groups. Concomitant activating PDGFRA/MET alterations were detected in tumors that acquired a homozygous CDKN2A deletion. NRG3 gene was significantly downregulated and recurrently altered in progressed tumors. Its decreased expression was associated with poorer overall survival in both univariate and multivariate analysis. We also detected progression-related alterations in RAD51B and other DNA repair pathway genes associated with the promotion of error-prone DNA repair, potentially facilitating tumor progression. In our retrospective analysis of patient treatment and survival timelines (n = 75), the combination of postoperative radiation and chemotherapy (mainly TMZ) outperformed radiation, especially in the grade 3 tumor cohort, in which it was typically given after primary surgery. Our results provide further insight into the contribution of treatment and genetic alterations in cell cycle, growth factor signaling, and DNA repair-related genes to tumor evolution and progression.

The large diversity of central nervous system (CNS) tumor types in children and adolescents results in disparate patient outcomes and renders accurate diagnosis challenging. In this study, we prospectively integrated DNA methylation profiling and targeted gene panel sequencing with blinded neuropathological reference diagnostics for a population-based cohort of more than 1,200 newly diagnosed pediatric patients with CNS tumors, to assess their utility in routine neuropathology. We show that the multi-omic integration increased diagnostic accuracy in a substantial proportion of patients through annotation to a refining DNA methylation class (50%), detection of diagnostic or therapeutically relevant genetic alterations (47%) or identification of cancer predisposition syndromes (10%). Discrepant results by neuropathological WHO-based and DNA methylation-based classification (30%) were enriched in histological high-grade gliomas, implicating relevance for current clinical patient management in 5% of all patients. Follow-up (median 2.5 years) suggests improved survival for patients with histological high-grade gliomas displaying lower-grade molecular profiles. These results provide preliminary evidence of the utility of integrating multi-omics in neuropathology for pediatric neuro-oncology.

BRAF genomic alterations are the most common oncogenic drivers in pediatric low-grade glioma (pLGG). Arm 1 (n = 77) of the ongoing phase 2 FIREFLY-1 (PNOC026) trial investigated the efficacy of the oral, selective, central nervous system-penetrant, type II RAF inhibitor tovorafenib (420 mg m-2 once weekly; 600 mg maximum) in patients with BRAF-altered, relapsed/refractory pLGG. Arm 2 (n = 60) is an extension cohort, which provided treatment access for patients with RAF-altered pLGG after arm 1 closure. Based on independent review, according to Response Assessment in Neuro-Oncology High-Grade Glioma (RANO-HGG) criteria, the overall response rate (ORR) of 67% met the arm 1 prespecified primary endpoint; median duration of response (DOR) was 16.6 months; and median time to response (TTR) was 3.0 months (secondary endpoints). Other select arm 1 secondary endpoints included ORR, DOR and TTR as assessed by Response Assessment in Pediatric Neuro-Oncology Low-Grade Glioma (RAPNO) criteria and safety (assessed in all treated patients and the primary endpoint for arm 2, n = 137). The ORR according to RAPNO criteria (including minor responses) was 51%; median DOR was 13.8 months; and median TTR was 5.3 months. The most common treatment-related adverse events (TRAEs) were hair color changes (76%), elevated creatine phosphokinase (56%) and anemia (49%). Grade ≥3 TRAEs occurred in 42% of patients. Nine (7%) patients had TRAEs leading to discontinuation of tovorafenib. These data indicate that tovorafenib could be an effective therapy for BRAF-altered, relapsed/refractory pLGG. ClinicalTrials.gov registration: NCT04775485 .

The response of triple-negative breast cancer (TNBC) patients to pre-operative (neoadjuvant chemotherapy) is a critical factor of their outcome. To determine the effects of chemotherapy on the tumor genome and to identify mutations associated with chemoresistance and sensitivity, we performed whole exome sequencing on pre/post-chemotherapy tumors and matched lymphocytes from 26 patients. We observed great inter-tumoral heterogeneity with no gene mutated recurrently in more than four tumors besides TP53. Although the degree of response to chemotherapy in residual tumors was associated with more subclonal changes during chemotherapy, there was minimal evolution between pre/post-tumors. Indeed, gene sets enriched for mutations in pre- and post-chemotherapy tumors were very similar and reflected genes involved in the biological process of neurogenesis. Somatically mutated genes present in chemosensitive tumors included COL1A2, PRMD15, APOBEC3B, PALB2 and histone protein encoding genes, while BRCA1, ATR, ARID1A, XRCC3 and genes encoding for tubulin-associated proteins were present in the chemoresistant tumors. We also found that the mutational spectrum of post-chemotherapy tumors was more reflective of matching metastatic tumor biopsies than pre-chemotherapy samples. These findings support a portrait of modest ongoing genomic instability with respect to single-nucleotide variants induced by or selected for by chemotherapy in TNBCs.

Posterior fossa type A (PF-EPN-A, PFA) ependymoma are aggressive tumors that mainly affect children and have a poor prognosis. Histopathology shows significant intratumoral heterogeneity, ranging from loose tissue to often sharply demarcated, extremely cell-dense tumor areas. To determine molecular differences in morphologically different areas and to understand their clinical significance, we analyzed 113 PF-EPN-A samples, including 40 corresponding relapse samples. Cell-dense areas ranged from 0 to 100% of the tumor area and displayed a higher proportion of proliferating tumor cells (p < 0.01). Clinically, cell density was associated with poor progression-free and overall survival (pPFS = 0.0026, pOS < 0.01). Molecularly, tumor areas with low and high cell density showed diverging DNA methylation profiles regarding their similarity to distinct previously discovered PF-EPN-A subtypes in 9/21 cases. Prognostically relevant chromosomal changes at 1q and 6q showed spatial heterogeneity within single tumors and were significantly enriched in cell-dense tumor areas as shown by single-cell RNA (scRNA)-sequencing as well as copy number profiling and fluorescence in situ hybridization (FISH) analyses of different tumor areas. Finally, spatial transcriptomics revealed cell-dense areas of different tumors to be more similar than various different areas of the same tumor. High-density areas distinctly overexpressed genes encoding histone proteins, WNT5A, TGFB1, or IGF2. Relapsing tumors displayed a higher proportion of cell-dense areas (p = 0.036), a change in PF-EPN-A methylation subtypes (13/32 patients), and novel chromosome 1q gains and 6q losses (12/32 cases) compared to corresponding primary tumors. Our data suggest that PF-EPN-A ependymomas habor a previously unrecognized intratumoral heterogeneity with clinical implications, which has to be accounted for when selecting diagnostic material, inter alia, by histological evaluation of the proportion of cell-dense areas.

Diffuse intrinsic pontine gliomas (DIPGs) are deadly pediatric brain tumors, non-resectable due to brainstem localization and diffusive growth. Over 80% of DIPGs harbor a mutation in histone 3 (H3.3 or H3.1) resulting in a lysine-to-methionine substitution (H3K27M). Patients with DIPG have a dismal prognosis with no effective therapy. We show that histone deacetylase (HDAC) inhibitors lead to a significant reduction in the H3.3K27M protein (up to 80%) in multiple glioma cell lines. We discover that the SB939-mediated H3.3K27M loss is partially blocked by a lysosomal inhibitor, chloroquine. The H3.3K27M loss is facilitated by co-occurrence of H2A.Z, as evidenced by the knockdown of H2A.Z isoforms. Chromatin immunoprecipitation sequencing (ChIP-seq) analysis confirms the occupancy of H3.3K27M and H2A.Z at the same SB939-inducible genes. We discover a mechanism showing that HDAC inhibition in DIPG leads to pharmacological modulation of the oncogenic H3.3K27M protein levels. These findings show the possibility of directly targeting the H3.3K27M oncohistone.

The KCC2 cotransporter establishes the low neuronal Cl(-) levels required for GABAA and glycine (Gly) receptor-mediated inhibition, and KCC2 deficiency in model organisms results in network hyperexcitability. However, no mutations in KCC2 have been documented in human disease. Here, we report two non-synonymous functional variants in human KCC2, R952H and R1049C, exhibiting clear statistical association with idiopathic generalized epilepsy (IGE). These variants reside in conserved residues in the KCC2 cytoplasmic C-terminus, exhibit significantly impaired Cl(-)-extrusion capacities resulting in less hyperpolarized Gly equilibrium potentials (EG ly), and impair KCC2 stimulatory phosphorylation at serine 940, a key regulatory site. These data describe a novel KCC2 variant significantly associated with a human disease and suggest genetically encoded impairment of KCC2 functional regulation may be a risk factor for the development of human IGE.

The cure rate for childhood intracranial ependymoma is approximately 70% in the setting of a gross total resection followed by radiation, but management remains challenging in patients with residual disease. Therefore, robust biomarkers are needed to guide the development of new targeted therapy. The authors evaluated the expression of several biomarkers in pediatric intracranial ependymoma and observed that the expression of enhancer of zeste homolog 2 (EZH2), a polycomb complex protein involved in epigenetic regulation of gene expression, was independently associated with poor survival.

Diffuse Intrinsic Pontine Gliomas (DIPGs) are deadly paediatric brain tumours where needle biopsies help guide diagnosis and targeted therapies. To address spatial heterogeneity, here we analyse 134 specimens from various neuroanatomical structures of whole autopsy brains from nine DIPG patients. Evolutionary reconstruction indicates histone 3 (H3) K27M--including H3.2K27M--mutations potentially arise first and are invariably associated with specific, high-fidelity obligate partners throughout the tumour and its spread, from diagnosis to end-stage disease, suggesting mutual need for tumorigenesis. These H3K27M ubiquitously-associated mutations involve alterations in TP53 cell-cycle (TP53/PPM1D) or specific growth factor pathways (ACVR1/PIK3R1). Later oncogenic alterations arise in sub-clones and often affect the PI3K pathway. Our findings are consistent with early tumour spread outside the brainstem including the cerebrum. The spatial and temporal homogeneity of main driver mutations in DIPG implies they will be captured by limited biopsies and emphasizes the need to develop therapies specifically targeting obligate oncohistone partnerships.

Encephalocraniocutaneous lipomatosis (ECCL) is a sporadic condition characterized by ocular, cutaneous, and central nervous system anomalies. Key clinical features include a well-demarcated hairless fatty nevus on the scalp, benign ocular tumors, and central nervous system lipomas. Seizures, spasticity, and intellectual disability can be present, although affected individuals without seizures and with normal intellect have also been reported. Given the patchy and asymmetric nature of the malformations, ECCL has been hypothesized to be due to a post-zygotic, mosaic mutation. Despite phenotypic overlap with several other disorders associated with mutations in the RAS-MAPK and PI3K-AKT pathways, the molecular etiology of ECCL remains unknown. Using exome sequencing of DNA from multiple affected tissues from five unrelated individuals with ECCL, we identified two mosaic mutations, c.1638C>A (p.Asn546Lys) and c.1966A>G (p.Lys656Glu) within the tyrosine kinase domain of FGFR1, in two affected individuals each. These two residues are the most commonly mutated residues in FGFR1 in human cancers and are associated primarily with CNS tumors. Targeted resequencing of FGFR1 in multiple tissues from an independent cohort of individuals with ECCL identified one additional individual with a c.1638C>A (p.Asn546Lys) mutation in FGFR1. Functional studies of ECCL fibroblast cell lines show increased levels of phosphorylated FGFRs and phosphorylated FRS2, a direct substrate of FGFR1, as well as constitutive activation of RAS-MAPK signaling. In addition to identifying the molecular etiology of ECCL, our results support the emerging overlap between mosaic developmental disorders and tumorigenesis.

The emergence of isoform-sensitive microarrays has helped fuel in-depth studies of the human transcriptome. The Affymetrix GeneChip Human Exon 1.0 ST Array (Exon Array) has been previously shown to be effective in profiling gene expression at the isoform level. More recently, the Affymetrix GeneChip Human Gene 1.0 ST Array (Gene Array) has been released for measuring gene expression and interestingly contains a large subset of probes from the Exon Array. Here, we explore the potential of using Gene Array probes to assess expression variation at the sub-transcript level. Utilizing datasets of the high quality Microarray Quality Control (MAQC) RNA samples previously assayed on the Exon Array and Gene Array, we compare the expression measurements of the two platforms to determine the performance of the Gene Array in detecting isoform variations.

We recently described a novel autosomal recessive neurodevelopmental disorder with intellectual disability in four patients from two related Hutterite families. Identity-by-descent mapping localized the gene to a 5.1 Mb region at chromosome 16p13.3 containing more than 170 known or predicted genes. The objective of this study was to identify the causative gene for this rare disorder.

Herpes simplex encephalitis (HSE) is a lethal neurological disease resulting from infection with Herpes Simplex Virus 1 (HSV-1). Loss-of-function mutations in the UNC93B1, TLR3, TRIF, TRAF3, and TBK1 genes have been associated with a human genetic predisposition to HSE, demonstrating the UNC93B-TLR3-type I IFN pathway as critical in protective immunity to HSV-1. However, the TLR3, UNC93B1, and TRIF mutations exhibit incomplete penetrance and represent only a minority of HSE cases, perhaps reflecting the effects of additional host genetic factors. In order to identify new host genes, proteins and signaling pathways involved in HSV-1 and HSE susceptibility, we have implemented the first genome-wide mutagenesis screen in an in vivo HSV-1 infectious model. One pedigree (named P43) segregated a susceptible trait with a fully penetrant phenotype. Genetic mapping and whole exome sequencing led to the identification of the causative nonsense mutation L3X in the Receptor-type tyrosine-protein phosphatase C gene (Ptprc(L3X)), which encodes for the tyrosine phosphatase CD45. Expression of MCP1, IL-6, MMP3, MMP8, and the ICP4 viral gene were significantly increased in the brain stems of infected Ptprc(L3X) mice accounting for hyper-inflammation and pathological damages caused by viral replication. Ptprc(L3X) mutation drastically affects the early stages of thymocytes development but also the final stage of B cell maturation. Transfer of total splenocytes from heterozygous littermates into Ptprc(L3X) mice resulted in a complete HSV-1 protective effect. Furthermore, T cells were the only cell population to fully restore resistance to HSV-1 in the mutants, an effect that required both the CD4⁺ and CD8⁺ T cells and could be attributed to function of CD4⁺ T helper 1 (Th1) cells in CD8⁺ T cell recruitment to the site of infection. Altogether, these results revealed the CD45-mediated T cell function as potentially critical for infection and viral spread to the brain, and also for subsequent HSE development.

Recently, thanks to the increasing throughput of new technologies, we have begun to explore the full extent of alternative pre-mRNA splicing (AS) in the human transcriptome. This is unveiling a vast layer of complexity in isoform-level expression differences between individuals. We used previously published splicing sensitive microarray data from lymphoblastoid cell lines to conduct an in-depth analysis on splicing efficiency of known and predicted exons. By combining publicly available AS annotation with a novel algorithm designed to search for AS, we show that many real AS events can be detected within the usually unexploited, speculative majority of the array and at significance levels much below standard multiple-testing thresholds, demonstrating that the extent of cis-regulated differential splicing between individuals is potentially far greater than previously reported. Specifically, many genes show subtle but significant genetically controlled differences in splice-site usage. PCR validation shows that 42 out of 58 (72%) candidate gene regions undergo detectable AS, amounting to the largest scale validation of isoform eQTLs to date. Targeted sequencing revealed a likely causative SNP in most validated cases. In all 17 incidences where a SNP affected a splice-site region, in silico splice-site strength modeling correctly predicted the direction of the micro-array and PCR results. In 13 other cases, we identified likely causative SNPs disrupting predicted splicing enhancers. Using Fst and REHH analysis, we uncovered significant evidence that 2 putative causative SNPs have undergone recent positive selection. We verified the effect of five SNPs using in vivo minigene assays. This study shows that splicing differences between individuals, including quantitative differences in isoform ratios, are frequent in human populations and that causative SNPs can be identified using in silico predictions. Several cases affected disease-relevant genes and it is likely some of these differences are involved in phenotypic diversity and susceptibility to complex diseases.

There is tremendous potential for genome sequencing to improve clinical diagnosis and care once it becomes routinely accessible, but this will require formalizing research methods into clinical best practices in the areas of sequence data generation, analysis, interpretation and reporting. The CLARITY Challenge was designed to spur convergence in methods for diagnosing genetic disease starting from clinical case history and genome sequencing data. DNA samples were obtained from three families with heritable genetic disorders and genomic sequence data were donated by sequencing platform vendors. The challenge was to analyze and interpret these data with the goals of identifying disease-causing variants and reporting the findings in a clinically useful format. Participating contestant groups were solicited broadly, and an independent panel of judges evaluated their performance.