New Brain Tumor Entities Emerge from Molecular Classification of CNS-PNETs.

Primitive neuroectodermal tumors of the central nervous system (CNS-PNETs) are highly aggressive, poorly differentiated embryonal tumors occurring predominantly in young children but also affecting adolescents and adults. Herein, we demonstrate that a significant proportion of institutionally diagnosed CNS-PNETs display molecular profiles indistinguishable from those of various other well-defined CNS tumor entities, facilitating diagnosis and appropriate therapy for patients with these tumors. From the remaining fraction of CNS-PNETs, we identify four new CNS tumor entities, each associated with a recurrent genetic alteration and distinct histopathological and clinical features. These new molecular entities, designated "CNS neuroblastoma with FOXR2 activation (CNS NB-FOXR2)," "CNS Ewing sarcoma family tumor with CIC alteration (CNS EFT-CIC)," "CNS high-grade neuroepithelial tumor with MN1 alteration (CNS HGNET-MN1)," and "CNS high-grade neuroepithelial tumor with BCOR alteration (CNS HGNET-BCOR)," will enable meaningful clinical trials and the development of therapeutic strategies for patients affected by poorly differentiated CNS tumors.

Accurate and efficient detection of gene fusions from RNA sequencing data.

The identification of gene fusions from RNA sequencing data is a routine task in cancer research and precision oncology. However, despite the availability of many computational tools, fusion detection remains challenging. Existing methods suffer from poor prediction accuracy and are computationally demanding. We developed Arriba, a novel fusion detection algorithm with high sensitivity and short runtime. When applied to a large collection of published pancreatic cancer samples (n = 803), Arriba identified a variety of driver fusions, many of which affected druggable proteins, including ALK, BRAF, FGFR2, NRG1, NTRK1, NTRK3, RET, and ROS1. The fusions were significantly associated with KRAS wild-type tumors and involved proteins stimulating the MAPK signaling pathway, suggesting that they substitute for activating mutations in KRAS In addition, we confirmed the transforming potential of two novel fusions, RRBP1-RAF1 and RASGRP1-ATP1A1, in cellular assays. These results show Arriba's utility in both basic cancer research and clinical translation.

Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas.

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

T-cell prolymphocytic leukemia is associated with deregulation of oncogenic microRNAs on transcriptional and epigenetic level.

Deregulation of micro(mi)-RNAs is a common mechanism in tumorigenesis. We investigated the expression of 2083 miRNAs in T-cell prolymphocytic leukemia (T-PLL). Compared to physiologic CD4+ and CD8+ T-cell subsets, 111 miRNAs were differentially expressed in T-PLL. Of these, 33 belonged to miRNA gene clusters linked to cancer. Genomic variants affecting miRNAs were infrequent with the notable exception of copy number aberrations. Remarkably, we found strong upregulation of the miR-200c/-141 cluster in T-PLL to be associated with DNA hypomethylation and active promoter marks. Our findings suggest that copy number aberrations and epigenetic changes could contribute to miRNA deregulation in T-PLL.

Comprehensive genomic analysis of refractory multiple myeloma reveals a complex mutational landscape associated with drug resistance and novel therapeutic vulnerabilities.

The outcomes of patients with multiple myeloma (MM) refractory to immunomodulatory agents (IMiDs) and proteasome inhibitors (PIs) remain poor. In this study, we performed whole genome and transcriptome sequencing of 39 heavily pretreated relapsed/refractory MM (RRMM) patients to identify mechanisms of resistance and potential therapeutic targets. We observed a high mutational load and indications of increased genomic instability. Recurrently mutated genes in RRMM, which had not been previously reported or only observed at a lower frequency in newly diagnosed MM, included NRAS, BRAF, TP53, SLC4A7, MLLT4, EWSR1, HCFC2, and COPS3. We found multiple genomic regions with bi-allelic events affecting tumor suppressor genes and demonstrated a significant adverse impact of bi-allelic TP53 alterations on survival. With regard to potentially resistance conferring mutations, recurrently mutated gene networks included genes with relevance for PI and IMiD activity; the latter particularly affecting members of the Cereblon and the COP9 signalosome complex. We observed a major impact of signatures associated with exposure to melphalan or impaired DNA double-strand break homologous recombination repair in RRMM. The latter coincided with mutations in genes associated with PARP inhibitor sensitivity in 49% of RRMM patients; a finding with potential therapeutic implications. In conclusion, this comprehensive genomic characterization revealed a complex mutational and structural landscape in RRMM and highlights potential implications for therapeutic strategies.

Reconstruction of rearranged T-cell receptor loci by whole genome and transcriptome sequencing gives insights into the initial steps of T-cell prolymphocytic leukemia.

T-cell prolymphocytic leukemia (T-PLL) is an aggressive tumor with leukemic presentation of mature T-lymphocytes. Here, we aimed at characterizing the initial events in the molecular pathogenesis of T-PLL and particularly, at determining the point in T-cell differentiation when the hallmark oncogenic events, that is, inv(14)(q11q32)/t(14;14)(q11;q32) and t(X;14)(q28;q11) occur. To this end, we mined whole genome and transcriptome sequencing data of 17 and 11 T-PLL cases, respectively. Mapping of the 14q32.1 locus breakpoints identified only TCL1A, which was moreover significantly overexpressed in T-PLL as compared to benign CD4+ and CD8+ T-cells, as the only common oncogenic target of aberrations. In cases with t(14;14), the breakpoints mapped telomeric and in cases with inv(14) centromeric or in the 3'-untranslated region of TCL1A. Regarding the T-cell receptor alpha (TRA) locus-TCL1A breakpoint junctions, all 17 breakpoints involved recombination signal sequences and 15 junctions contained nontemplated (N-) nucleotides. All T-PLL cases studied carried in-frame TRA rearrangements on the intact allele, which skewed significantly toward usage of distal/central TRAV/TRAJ gene segments as compared to the illegitimate TRA rearrangements. Our findings suggest that the oncogenic TRA-TCL1A/MTCP1 rearrangements in T-PLL occur during opening of the TRA locus, that is, during the progression from CD4+ immature single positive to early double positive thymocyte stage, just before physiologic TCL1A expression is silenced. The cell carrying such an oncogenic event continues maturation and rearranges the second TRA allele to achieve a functional T-cell receptor. Thereafter, it switches off RAG and DNTT expression in line with the mature T-cell phenotype at presentation of T-PLL.

IG-MYC + neoplasms with precursor B-cell phenotype are molecularly distinct from Burkitt lymphomas.

The WHO Classification of Tumours of Haematopoietic and Lymphoid Tissue notes instances of Burkitt lymphoma/leukemia (BL) with IG-MYC rearrangement displaying a B-cell precursor immunophenotype (termed herein "preBLL"). To characterize the molecular pathogenesis of preBLL, we investigated 13 preBLL cases (including 1 cell line), of which 12 were analyzable using genome, exome, and targeted sequencing, imbalance mapping, and DNA methylation profiling. In 5 patients with reads across the IG-MYC breakpoint junctions, we found evidence that the translocation derived from an aberrant VDJ recombination, as is typical for IG translocations arising in B-cell precursors. Genomic changes like biallelic IGH translocations or VDJ rearrangements combined with translocation into the VDJ region on the second allele, potentially preventing expression of a productive immunoglobulin, were detected in 6 of 13 cases. We did not detect mutations in genes frequently altered in BL, but instead found activating NRAS and/or KRAS mutations in 7 of 12 preBLLs. Gains on 1q, recurrent in BL and preB lymphoblastic leukemia/lymphoma (pB-ALL/LBL), were detected in 7 of 12 preBLLs. DNA methylation profiling showed preBLL to cluster with precursor B cells and pB-ALL/LBL, but apart from BL. We conclude that preBLL genetically and epigenetically resembles pB-ALL/LBL rather than BL. Therefore, we propose that preBLL be considered as a pB-ALL/LBL with recurrent genetic abnormalities.

Mutational patterns and regulatory networks in epigenetic subgroups of meningioma.

DNA methylation patterns delineate clinically relevant subgroups of meningioma. We previously established the six meningioma methylation classes (MC) benign 1-3, intermediate A and B, and malignant. Here, we set out to identify subgroup-specific mutational patterns and gene regulation. Whole genome sequencing was performed on 62 samples across all MCs and WHO grades from 62 patients with matched blood control, including 40 sporadic meningiomas and 22 meningiomas arising after radiation (Mrad). RNA sequencing was added for 18 of these cases and chromatin-immunoprecipitation for histone H3 lysine 27 acetylation (H3K27ac) followed by sequencing (ChIP-seq) for 16 samples. Besides the known mutations in meningioma, structural variants were found as the mechanism of NF2 inactivation in a small subset (5%) of sporadic meningiomas, similar to previous reports for Mrad. Aberrations of DMD were found to be enriched in MCs with NF2 mutations, and DMD was among the most differentially upregulated genes in NF2 mutant compared to NF2 wild-type cases. The mutational signature AC3, which has been associated with defects in homologous recombination repair (HRR), was detected in both sporadic meningioma and Mrad, but widely distributed across the genome in sporadic cases and enriched near genomic breakpoints in Mrad. Compared to the other MCs, the number of single nucleotide variants matching the AC3 pattern was significantly higher in the malignant MC, which also exhibited higher genomic instability, determined by the numbers of both large segments affected by copy number alterations and breakpoints between large segments. ChIP-seq analysis for H3K27ac revealed a specific activation of genes regulated by the transcription factor FOXM1 in the malignant MC. This analysis also revealed a super enhancer near the HOXD gene cluster in this MC, which, together with general upregulation of HOX genes in the malignant MC, indicates a role of HOX genes in meningioma aggressiveness. This data elucidates the biological mechanisms rendering different epigenetic subgroups of meningiomas, and suggests leveraging HRR as a novel therapeutic target.

Screening drug effects in patient-derived cancer cells links organoid responses to genome alterations.

Cancer drug screening in patient-derived cells holds great promise for personalized oncology and drug discovery but lacks standardization. Whether cells are cultured as conventional monolayer or advanced, matrix-dependent organoid cultures influences drug effects and thereby drug selection and clinical success. To precisely compare drug profiles in differently cultured primary cells, we developed DeathPro, an automated microscopy-based assay to resolve drug-induced cell death and proliferation inhibition. Using DeathPro, we screened cells from ovarian cancer patients in monolayer or organoid culture with clinically relevant drugs. Drug-induced growth arrest and efficacy of cytostatic drugs differed between the two culture systems. Interestingly, drug effects in organoids were more diverse and had lower therapeutic potential. Genomic analysis revealed novel links between drug sensitivity and DNA repair deficiency in organoids that were undetectable in monolayers. Thus, our results highlight the dependency of cytostatic drugs and pharmacogenomic associations on culture systems, and guide culture selection for drug tests.

Conserved peptide fragmentation as a benchmarking tool for mass spectrometers and a discriminating feature for targeted proteomics.

Quantifying the similarity of spectra is an important task in various areas of spectroscopy, for example, to identify a compound by comparing sample spectra to those of reference standards. In mass spectrometry based discovery proteomics, spectral comparisons are used to infer the amino acid sequence of peptides. In targeted proteomics by selected reaction monitoring (SRM) or SWATH MS, predetermined sets of fragment ion signals integrated over chromatographic time are used to identify target peptides in complex samples. In both cases, confidence in peptide identification is directly related to the quality of spectral matches. In this study, we used sets of simulated spectra of well-controlled dissimilarity to benchmark different spectral comparison measures and to develop a robust scoring scheme that quantifies the similarity of fragment ion spectra. We applied the normalized spectral contrast angle score to quantify the similarity of spectra to objectively assess fragment ion variability of tandem mass spectrometric datasets, to evaluate portability of peptide fragment ion spectra for targeted mass spectrometry across different types of mass spectrometers and to discriminate target assays from decoys in targeted proteomics. Altogether, this study validates the use of the normalized spectral contrast angle as a sensitive spectral similarity measure for targeted proteomics, and more generally provides a methodology to assess the performance of spectral comparisons and to support the rational selection of the most appropriate similarity measure. The algorithms used in this study are made publicly available as an open source toolset with a graphical user interface.

Modulation of the chromatin phosphoproteome by the Haspin protein kinase.

Recent discoveries have highlighted the importance of Haspin kinase activity for the correct positioning of the kinase Aurora B at the centromere. Haspin phosphorylates Thr(3) of the histone H3 (H3), which provides a signal for Aurora B to localize to the centromere of mitotic chromosomes. To date, histone H3 is the only confirmed Haspin substrate. We used a combination of biochemical, pharmacological, and mass spectrometric approaches to study the consequences of Haspin inhibition in mitotic cells. We quantified 3964 phosphorylation sites on chromatin-associated proteins and identified a Haspin protein-protein interaction network. We determined the Haspin consensus motif and the co-crystal structure of the kinase with the histone H3 tail. The structure revealed a unique bent substrate binding mode positioning the histone H3 residues Arg(2) and Lys(4) adjacent to the Haspin phosphorylated threonine into acidic binding pockets. This unique conformation of the kinase-substrate complex explains the reported modulation of Haspin activity by methylation of Lys(4) of the histone H3. In addition, the identification of the structural basis of substrate recognition and the amino acid sequence preferences of Haspin aided the identification of novel candidate Haspin substrates. In particular, we validated the phosphorylation of Ser(137) of the histone variant macroH2A as a target of Haspin kinase activity. MacroH2A Ser(137) resides in a basic stretch of about 40 amino acids that is required to stabilize extranucleosomal DNA, suggesting that phosphorylation of Ser(137) might regulate the interactions of macroH2A and DNA. Overall, our data suggest that Haspin activity affects the phosphorylation state of proteins involved in gene expression regulation and splicing.

The genomic and transcriptional landscape of primary central nervous system lymphoma.

Primary lymphomas of the central nervous system (PCNSL) are mainly diffuse large B-cell lymphomas (DLBCLs) confined to the central nervous system (CNS). Molecular drivers of PCNSL have not been fully elucidated. Here, we profile and compare the whole-genome and transcriptome landscape of 51 CNS lymphomas (CNSL) to 39 follicular lymphoma and 36 DLBCL cases outside the CNS. We find recurrent mutations in JAK-STAT, NFkB, and B-cell receptor signaling pathways, including hallmark mutations in MYD88 L265P (67%) and CD79B (63%), and CDKN2A deletions (83%). PCNSLs exhibit significantly more focal deletions of HLA-D (6p21) locus as a potential mechanism of immune evasion. Mutational signatures correlating with DNA replication and mitosis are significantly enriched in PCNSL. TERT gene expression is significantly higher in PCNSL compared to activated B-cell (ABC)-DLBCL. Transcriptome analysis clearly distinguishes PCNSL and systemic DLBCL into distinct molecular subtypes. Epstein-Barr virus (EBV)+ CNSL cases lack recurrent mutational hotspots apart from IG and HLA-DRB loci. We show that PCNSL can be clearly distinguished from DLBCL, having distinct expression profiles, IG expression and translocation patterns, as well as specific combinations of genetic alterations.

FOXR2 Stabilizes MYCN Protein and Identifies Non-MYCN-Amplified Neuroblastoma Patients With Unfavorable Outcome.

PURPOSE: Clinical outcomes of patients with neuroblastoma range from spontaneous tumor regression to fatality. Hence, understanding the mechanisms that cause tumor progression is crucial for the treatment of patients. In this study, we show that FOXR2 activation identifies a subset of neuroblastoma tumors with unfavorable outcome and we investigate the mechanism how FOXR2 relates to poor outcome in patients. MATERIALS AND METHODS: We analyzed three independent transcriptional data sets of in total 1030 primary neuroblastomas with full clinical annotation. We performed immunoprecipitation for FOXR2 and MYCN and silenced FOXR2 expression in two neuroblastoma cell lines to examine the effect on cellular processes, transcriptome, and MYCN protein levels. Tumor samples were analyzed for protein levels of FOXR2 and MYCN. RESULTS: In three combined neuroblastoma data sets, 9% of tumors show expression of FOXR2 but have low levels of MYCN mRNA. FOXR2 expression identifies a group of patients with unfavorable outcome, showing 10-year overall survival rates of 53%-59%, and proves to be an independent prognostic factor compared with established risk factors. Transcriptionally, FOXR2-expressing tumors are very similar to MYCN-amplified tumors, suggesting that they might share a common mechanism of tumor initiation. FOXR2 knockdown in FOXR2-expressing neuroblastoma cell lines resulted in cell cycle arrest, reduced cell growth, cell death, and reduced MYCN protein levels, all indicating that FOXR2 is essential for these tumors. Finally, we show that FOXR2 binds and stabilizes MYCN protein and MYCN protein levels are highly increased in FOXR2-expressing tumors, in several cases comparable with MYCN-amplified samples. CONCLUSION: The stabilization of MYCN by FOXR2 represents an alternative mechanism to MYCN amplification to increase MYCN protein levels. As such, FOXR2 expression identifies another subset of neuroblastoma patients with unfavorable clinical outcome.

Kalirin-RAC controls nucleokinetic migration in ADRN-type neuroblastoma.

The migrational propensity of neuroblastoma is affected by cell identity, but the mechanisms behind the divergence remain unknown. Using RNAi and time-lapse imaging, we show that ADRN-type NB cells exhibit RAC1- and kalirin-dependent nucleokinetic (NUC) migration that relies on several integral components of neuronal migration. Inhibition of NUC migration by RAC1 and kalirin-GEF1 inhibitors occurs without hampering cell proliferation and ADRN identity. Using three clinically relevant expression dichotomies, we reveal that most of up-regulated mRNAs in RAC1- and kalirin-GEF1-suppressed ADRN-type NB cells are associated with low-risk characteristics. The computational analysis shows that, in a context of overall gene set poverty, the upregulomes in RAC1- and kalirin-GEF1-suppressed ADRN-type cells are a batch of AU-rich element-containing mRNAs, which suggests a link between NUC migration and mRNA stability. Gene set enrichment analysis-based search for vulnerabilities reveals prospective weak points in RAC1- and kalirin-GEF1-suppressed ADRN-type NB cells, including activities of H3K27- and DNA methyltransferases. Altogether, these data support the introduction of NUC inhibitors into cancer treatment research.

Single-cell transcriptomic analyses provide insights into the developmental origins of neuroblastoma.

Neuroblastoma is a pediatric tumor of the developing sympathetic nervous system. However, the cellular origin of neuroblastoma has yet to be defined. Here we studied the single-cell transcriptomes of neuroblastomas and normal human developing adrenal glands at various stages of embryonic and fetal development. We defined normal differentiation trajectories from Schwann cell precursors over intermediate states to neuroblasts or chromaffin cells and showed that neuroblastomas transcriptionally resemble normal fetal adrenal neuroblasts. Importantly, neuroblastomas with varying clinical phenotypes matched different temporal states along normal neuroblast differentiation trajectories, with the degree of differentiation corresponding to clinical prognosis. Our work highlights the roles of oncogenic MYCN and loss of TFAP2B in blocking differentiation and may provide the basis for designing therapeutic interventions to overcome differentiation blocks.

Super enhancers define regulatory subtypes and cell identity in neuroblastoma.

Half of the children diagnosed with neuroblastoma (NB) have high-risk disease, disproportionately contributing to overall childhood cancer-related deaths. In addition to recurrent gene mutations, there is increasing evidence supporting the role of epigenetic deregulation in disease pathogenesis. Yet, comprehensive cis-regulatory network descriptions from NB are lacking. Here, using genome-wide H3K27ac profiles across 60 NBs, covering the different clinical and molecular subtypes, we identified four major super-enhancer-driven epigenetic subtypes and their underlying master regulatory networks. Three of these subtypes recapitulated known clinical groups; namely, MYCN-amplified, MYCN non-amplified high-risk and MYCN non-amplified low-risk NBs. The fourth subtype, exhibiting mesenchymal characteristics, shared cellular identity with multipotent Schwann cell precursors, was induced by RAS activation and was enriched in relapsed disease. Notably, CCND1, an essential gene in NB, was regulated by both mesenchymal and adrenergic regulatory networks converging on distinct super-enhancer modules. Overall, this study reveals subtype-specific super-enhancer regulation in NBs.

Alternative lengthening of telomeres in childhood neuroblastoma from genome to proteome.

Telomere maintenance by telomerase activation or alternative lengthening of telomeres (ALT) is a major determinant of poor outcome in neuroblastoma. Here, we screen for ALT in primary and relapsed neuroblastomas (n = 760) and characterize its features using multi-omics profiling. ALT-positive tumors are molecularly distinct from other neuroblastoma subtypes and enriched in a population-based clinical sequencing study cohort for relapsed cases. They display reduced ATRX/DAXX complex abundance, due to either ATRX mutations (55%) or low protein expression. The heterochromatic histone mark H3K9me3 recognized by ATRX is enriched at the telomeres of ALT-positive tumors. Notably, we find a high frequency of telomeric repeat loci with a neuroblastoma ALT-specific hotspot on chr1q42.2 and loss of the adjacent chromosomal segment forming a neo-telomere. ALT-positive neuroblastomas proliferate slowly, which is reflected by a protracted clinical course of disease. Nevertheless, children with an ALT-positive neuroblastoma have dismal outcome.

Mutational mechanisms shaping the coding and noncoding genome of germinal center derived B-cell lymphomas.

B cells have the unique property to somatically alter their immunoglobulin (IG) genes by V(D)J recombination, somatic hypermutation (SHM) and class-switch recombination (CSR). Aberrant targeting of these mechanisms is implicated in lymphomagenesis, but the mutational processes are poorly understood. By performing whole genome and transcriptome sequencing of 181 germinal center derived B-cell lymphomas (gcBCL) we identified distinct mutational signatures linked to SHM and CSR. We show that not only SHM, but presumably also CSR causes off-target mutations in non-IG genes. Kataegis clusters with high mutational density mainly affected early replicating regions and were enriched for SHM- and CSR-mediated off-target mutations. Moreover, they often co-occurred in loci physically interacting in the nucleus, suggesting that mutation hotspots promote increased mutation targeting of spatially co-localized loci (termed hypermutation by proxy). Only around 1% of somatic small variants were in protein coding sequences, but in about half of the driver genes, a contribution of B-cell specific mutational processes to their mutations was found. The B-cell-specific mutational processes contribute to both lymphoma initiation and intratumoral heterogeneity. Overall, we demonstrate that mutational processes involved in the development of gcBCL are more complex than previously appreciated, and that B cell-specific mutational processes contribute via diverse mechanisms to lymphomagenesis.