Search

Cancer immunotherapies produce remarkable results in B cell malignancies; however, optimal cell surface targets for many solid cancers remain elusive. Here, we present an integrative proteomic, transcriptomic, and epigenomic analysis of tumor and normal tissues to identify biologically relevant cell surface immunotherapeutic targets for neuroblastoma, an often-fatal childhood cancer. Proteogenomic analyses reveal sixty high-confidence candidate immunotherapeutic targets, and we prioritize delta-like canonical notch ligand 1 (DLK1) for further study. High expression of DLK1 directly correlates with a super-enhancer. Immunofluorescence, flow cytometry, and immunohistochemistry show robust cell surface expression of DLK1. Short hairpin RNA mediated silencing of DLK1 in neuroblastoma cells results in increased cellular differentiation. ADCT-701, a DLK1-targeting antibody-drug conjugate (ADC), shows potent and specific cytotoxicity in DLK1-expressing neuroblastoma xenograft models. Since high DLK1 expression is found in several adult and pediatric cancers, our study demonstrates the utility of a proteogenomic approach and credentials DLK1 as an immunotherapeutic target., Competing Interests: Declaration of interests F. Zammarchi, K. Havenith, and P.H.v.B. are or were employed by ADC Therapeutics at the time the work was conducted and hold or previously held shares/stocks in ADC Therapeutics. The following patent is held: WO2018146199A1., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Pediatric brain cancer is the leading cause of disease-related mortality in children, and many aggressive tumors still lack effective treatment strategies. We characterized aberrant alternative splicing across pediatric brain tumors, identifying pediatric high-grade gliomas (HGGs) among the most heterogeneous. Annotating these events with UniProt, we identified 11,940 splice events in 5,368 genes leading to potential protein function changes. We discovered CDC-like kinase 1 ( CLK1 ) is aberrantly spliced to include exon 4, resulting in a gain of two phosphorylation sites and subsequent activation. Inhibition of CLK1 with Cirtuvivint significantly decreased both cell viability and proliferation in the pediatric HGG KNS-42 cell line. Morpholino-mediated depletion of CLK1 exon 4 splicing reduced RNA expression, protein abundance, and cell viability with concurrent differential expression of 78 cancer genes and differential splicing at functional sites in 193 cancer genes. Our findings highlight a dependency of pediatric HGGs on CLK1 and represent a promising therapeutic strategy., Competing Interests: Conflicts of Interest The authors declare no conflicts of interest.

Background: In 2019, the Open Pediatric Brain Tumor Atlas (OpenPBTA) was created as a global, collaborative open-science initiative to genomically characterize 1,074 pediatric brain tumors and 22 patient-derived cell lines. Here, we extend the OpenPBTA to create the Open Pediatric Cancer (OpenPedCan) Project, a harmonized open-source multi-omic dataset from 6,112 pediatric cancer patients with 7,096 tumor events across more than 100 histologies. Combined with RNA-Seq from the Genotype-Tissue Expression (GTEx) and The Cancer Genome Atlas (TCGA), OpenPedCan contains nearly 48,000 total biospecimens (24,002 tumor and 23,893 normal specimens)., Findings: We utilized Gabriella Miller Kids First (GMKF) workflows to harmonize WGS, WXS, RNA-seq, and Targeted Sequencing datasets to include somatic SNVs, InDels, CNVs, SVs, RNA expression, fusions, and splice variants. We integrated summarized CPTAC whole cell proteomics and phospho-proteomics data, miRNA-Seq data, and have developed a methylation array harmonization workflow to include m-values, beta-vales, and copy number calls. OpenPedCan contains reproducible, dockerized workflows in GitHub, CAVATICA, and Amazon Web Services (AWS) to deliver harmonized and processed data from over 60 scalable modules which can be leveraged both locally and on AWS. The processed data are released in a versioned manner and accessible through CAVATICA or AWS S3 download (from GitHub), and queryable through PedcBioPortal and the NCI's pediatric Molecular Targets Platform. Notably, we have expanded PBTA molecular subtyping to include methylation information to align with the WHO 2021 Central Nervous System Tumor classifications, allowing us to create research- grade integrated diagnoses for these tumors., Conclusions: OpenPedCan data and its reproducible analysis module framework are openly available and can be utilized and/or adapted by researchers to accelerate discovery, validation, and clinical translation., Competing Interests: Declarations of Interest The authors declare no conflicts.

Summary: With the increasing rates of exome and whole genome sequencing, the ability to classify large sets of germline sequencing variants using up-to-date American College of Medical Genetics-Association for Molecular Pathology (ACMG-AMP) criteria is crucial. Here, we present Automated Germline Variant Pathogenicity (AutoGVP), a tool that integrates germline variant pathogenicity annotations from ClinVar and sequence variant classifications from a modified version of InterVar (PVS1 strength adjustments, removal of PP5/BP6). This tool facilitates large-scale, clinically focused classification of germline sequence variants in a research setting., Availability and Implementation: AutoGVP is an open source dockerized workflow implemented in R and freely available on GitHub at https://github.com/diskin-lab-chop/AutoGVP., (© The Author(s) 2024. Published by Oxford University Press.)

Purpose: The PNOC001 phase II single-arm trial sought to estimate progression-free survival (PFS) associated with everolimus therapy for progressive/recurrent pediatric low-grade glioma (pLGG) on the basis of phosphatidylinositol 3-kinase (PI3K)/AKT/mammalian target of rapamycin (mTOR) pathway activation as measured by phosphorylated-ribosomal protein S6 and to identify prognostic and predictive biomarkers., Patients and Methods: Patients, age 3-21 years, with progressive/recurrent pLGG received everolimus orally, 5 mg/m 2 once daily. Frequency of driver gene alterations was compared among independent pLGG cohorts of newly diagnosed and progressive/recurrent patients. PFS at 6 months (primary end point) and median PFS (secondary end point) were estimated for association with everolimus therapy., Results: Between 2012 and 2019, 65 subjects with progressive/recurrent pLGG (median age, 9.6 years; range, 3.0-19.9; 46% female) were enrolled, with a median follow-up of 57.5 months. The 6-month PFS was 67.4% (95% CI, 60.0 to 80.0) and median PFS was 11.1 months (95% CI, 7.6 to 19.8). Hypertriglyceridemia was the most common grade ≥3 adverse event. PI3K/AKT/mTOR pathway activation did not correlate with clinical outcomes (6-month PFS, active 68.4% v nonactive 63.3%; median PFS, active 11.2 months v nonactive 11.1 months; P = .80). Rare/novel KIAA1549::BRAF fusion breakpoints were most frequent in supratentorial midline pilocytic astrocytomas, in patients with progressive/recurrent disease, and correlated with poor clinical outcomes (median PFS, rare/novel KIAA1549::BRAF fusion breakpoints 6.1 months v common KIAA1549::BRAF fusion breakpoints 16.7 months; P < .05). Multivariate analysis confirmed their independent risk factor status for disease progression in PNOC001 and other, independent cohorts. Additionally, rare pathogenic germline variants in homologous recombination genes were identified in 6.8% of PNOC001 patients., Conclusion: Everolimus is a well-tolerated therapy for progressive/recurrent pLGGs. Rare/novel KIAA1549::BRAF fusion breakpoints may define biomarkers for progressive disease and should be assessed in future clinical trials.

Background: High-risk neuroblastoma is a complex genetic disease that is lethal in more than 50% of patients despite intense multimodal therapy. Through genome-wide association studies (GWAS) and next-generation sequencing, we have identified common single nucleotide polymorphisms and rare, pathogenic or likely pathogenic germline loss-of-function variants in BARD1 enriched in neuroblastoma patients. The functional implications of these findings remain poorly understood., Methods: We correlated BARD1 genotype with expression in normal tissues and neuroblastomas, along with the burden of DNA damage in tumors. To validate the functional consequences of germline pathogenic or likely pathogenic BARD1 variants, we used CRISPR-Cas9 to generate isogenic neuroblastoma (IMR-5) and control (RPE1) cellular models harboring heterozygous BARD1 loss-of-function variants (R112*, R150*, E287fs, and Q564*) and quantified genomic instability in these cells via next-generation sequencing and with functional assays measuring the efficiency of DNA repair., Results: Both common and rare neuroblastoma-associated BARD1 germline variants were associated with lower levels of BARD1 mRNA and an increased burden of DNA damage. Using isogenic heterozygous BARD1 loss-of-function variant cellular models, we functionally validated this association with inefficient DNA repair. BARD1 loss-of-function variant isogenic cells exhibited reduced efficiency in repairing Cas9-induced DNA damage, ineffective RAD51 focus formation at DNA double-strand break sites, and enhanced sensitivity to cisplatin and poly (ADP-ribose) polymerase (PARP) inhibition both in vitro and in vivo., Conclusions: Taken together, we demonstrate that germline BARD1 variants disrupt DNA repair fidelity. This is a fundamental molecular mechanism contributing to neuroblastoma initiation that may have important therapeutic implications., (© The Author(s) 2023. Published by Oxford University Press. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Background: To achieve replicative immortality, most cancers develop a telomere maintenance mechanism, such as reactivation of telomerase or alternative lengthening of telomeres (ALT). There are limited data on the prevalence and clinical significance of ALT in pediatric brain tumors, and ALT-directed therapy is not available., Methods: We performed C-circle analysis (CCA) on 579 pediatric brain tumors that had corresponding tumor/normal whole genome sequencing through the Open Pediatric Brain Tumor Atlas (OpenPBTA). We detected ALT in 6.9% (n = 40/579) of these tumors and completed additional validation by ultrabright telomeric foci in situ on a subset of these tumors. We used CCA to validate TelomereHunter for computational prediction of ALT status and focus subsequent analyses on pediatric high-grade gliomas (pHGGs) Finally, we examined whether ALT is associated with recurrent somatic or germline alterations., Results: ALT is common in pHGGs (n = 24/63, 38.1%), but occurs infrequently in other pediatric brain tumors (<3%). Somatic ATRX mutations occur in 50% of ALT+ pHGGs and in 30% of ALT- pHGGs. Rare pathogenic germline variants in mismatch repair (MMR) genes are significantly associated with an increased occurrence of ALT., Conclusions: We demonstrate that ATRX is mutated in only a subset of ALT+ pHGGs, suggesting other mechanisms of ATRX loss of function or alterations in other genes may be associated with the development of ALT in these patients. We show that germline variants in MMR are associated with the development of ALT in patients with pHGG., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Society for Neuro-Oncology.)

Pediatric brain and spinal cancers are collectively the leading disease-related cause of death in children; thus, we urgently need curative therapeutic strategies for these tumors. To accelerate such discoveries, the Children's Brain Tumor Network (CBTN) and Pacific Pediatric Neuro-Oncology Consortium (PNOC) created a systematic process for tumor biobanking, model generation, and sequencing with immediate access to harmonized data. We leverage these data to establish OpenPBTA, an open collaborative project with over 40 scalable analysis modules that genomically characterize 1,074 pediatric brain tumors. Transcriptomic classification reveals universal TP53 dysregulation in mismatch repair-deficient hypermutant high-grade gliomas and TP53 loss as a significant marker for poor overall survival in ependymomas and H3 K28-mutant diffuse midline gliomas. Already being actively applied to other pediatric cancers and PNOC molecular tumor board decision-making, OpenPBTA is an invaluable resource to the pediatric oncology community., Competing Interests: C.S.G.’s spouse was an employee of Alex’s Lemonade Stand Foundation, which was a sponsor of this research. J.A.S., C.L.S., C.J.B., S.J.S., and J.N.T. are or were employees of Alex’s Lemonade Stand Foundation, a sponsor of this research. A.J.W. is a member of the Scientific Advisory boards for Alexion and DayOne Biopharmaceuticals., (© 2023 The Author(s).)

Nicotine is the reinforcing ingredient in tobacco. Following chronic exposure, sudden cessation of nicotine use produces negative symptoms of withdrawal that contribute to dependence. The molecular mechanisms underlying nicotine withdrawal behaviors, however, are poorly understood. Using recombinant inbred mice, chronic nicotine was delivered by minipump and withdrawal induced using mecamylamine. Somatic signs of withdrawal, and anxiety-like behavior using elevated plus maze, were then assessed. Interval mapping was used to identify associations between genetic variation and withdrawal behaviors, and with basal gene expression. Differential gene expression following nicotine exposure and withdrawal was also assessed in progenitor mice using microarrays. Quantitative trait loci mapping identified chromosome intervals with significant genetic associations to somatic signs of withdrawal or withdrawal-induced anxiety-like behavior. Using bioinformatics, and association with basal gene expression in nucleus accumbens, we implicated Rb1, Bnip3l, Pnma2, Itm2b, and Kif13b as candidate genes for somatic signs of withdrawal, and Galr1, which showed trans-regulation from a region of chromosome 14 that was associated with somatic signs of withdrawal. Candidate genes within the chromosome 9 region associated with anxiety-like withdrawal behavior included Dixdc1, Ncam1, and Sorl1. Bioinformatics identified six genes that were also significantly associated with nicotine or alcohol traits in recent human genome-wide association studies. Withdrawal-associated somatic signs and anxiety-like behavior had strong non-overlapping genetic associations, respectively, with regions of chromosome 14 and chromosome 9. Genetic, behavioral and gene expression correlations, and bioinformatics analysis identified several candidate genes that may represent novel molecular targets for modulating nicotine withdrawal symptoms., (© 2023 The Authors. Genes, Brain and Behavior published by International Behavioural and Neural Genetics Society and John Wiley & Sons Ltd.)

Pediatric brain tumors are the leading cause of cancer-related death in children in the United States and contribute a disproportionate number of potential years of life lost compared to adult cancers. Moreover, survivors frequently suffer long-term side effects, including secondary cancers. The Children's Brain Tumor Network (CBTN) is a multi-institutional international clinical research consortium created to advance therapeutic development through the collection and rapid distribution of biospecimens and data via open-science research platforms for real-time access and use by the global research community. The CBTN's 32 member institutions utilize a shared regulatory governance architecture at the Children's Hospital of Philadelphia to accelerate and maximize the use of biospecimens and data. As of August 2022, CBTN has enrolled over 4700 subjects, over 1500 parents, and collected over 65,000 biospecimen aliquots for research. Additionally, over 80 preclinical models have been developed from collected tumors. Multi-omic data for over 1000 tumors and germline material are currently available with data generation for > 5000 samples underway. To our knowledge, CBTN provides the largest open-access pediatric brain tumor multi-omic dataset annotated with longitudinal clinical and outcome data, imaging, associated biospecimens, child-parent genomic pedigrees, and in vivo and in vitro preclinical models. Empowered by NIH-supported platforms such as the Kids First Data Resource and the Childhood Cancer Data Initiative, the CBTN continues to expand the resources needed for scientists to accelerate translational impact for improved outcomes and quality of life for children with brain and spinal cord tumors., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. David S. Ziegler is a consultant, or on the advisory board, of Bayer, AstraZeneca, Accendatech, Novartis, Day One, FivePhusion, Amgen, Alexion, and Norgine. Angela J. Waanders is on the advisory board of Alexion and Day One., (Copyright © 2022. Published by Elsevier Inc.)

Purpose: PNOC003 is a multicenter precision medicine trial for children and young adults with newly diagnosed diffuse intrinsic pontine glioma (DIPG)., Patients and Methods: Patients (3-25 years) were enrolled on the basis of imaging consistent with DIPG. Biopsy tissue was collected for whole-exome and mRNA sequencing. After radiotherapy (RT), patients were assigned up to four FDA-approved drugs based on molecular tumor board recommendations. H3K27M-mutant circulating tumor DNA (ctDNA) was longitudinally measured. Tumor tissue and matched primary cell lines were characterized using whole-genome sequencing and DNA methylation profiling. When applicable, results were verified in an independent cohort from the Children's Brain Tumor Network (CBTN)., Results: Of 38 patients enrolled, 28 patients (median 6 years, 10 females) were reviewed by the molecular tumor board. Of those, 19 followed treatment recommendations. Median overall survival (OS) was 13.1 months [95% confidence interval (CI), 11.2-18.4] with no difference between patients who followed recommendations and those who did not. H3K27M-mutant ctDNA was detected at baseline in 60% of cases tested and associated with response to RT and survival. Eleven cell lines were established, showing 100% fidelity of key somatic driver gene alterations in the primary tumor. In H3K27-altered DIPGs, TP53 mutations were associated with worse OS (TP53mut 11.1 mo; 95% CI, 8.7-14; TP53wt 13.3 mo; 95% CI, 11.8-NA; P = 3.4e-2), genome instability (P = 3.1e-3), and RT resistance (P = 6.4e-4). The CBTN cohort confirmed an association between TP53 mutation status, genome instability, and clinical outcome., Conclusions: Upfront treatment-naïve biopsy provides insight into clinically relevant molecular alterations and prognostic biomarkers for H3K27-altered DIPGs., (©2022 The Authors; Published by the American Association for Cancer Research.)

Background: Pediatric brain tumors are the leading cause of cancer death in children with an urgent need for innovative therapies. Glypican 2 (GPC2) is a cell surface oncoprotein expressed in neuroblastoma for which targeted immunotherapies have been developed. This work aimed to characterize GPC2 expression in pediatric brain tumors and develop an mRNA CAR T cell approach against this target., Methods: We investigated GPC2 expression across a cohort of primary pediatric brain tumor samples and cell lines using RNA sequencing, immunohistochemistry, and flow cytometry. To target GPC2 in the brain with adoptive cellular therapies and mitigate potential inflammatory neurotoxicity, we used optimized mRNA to create transient chimeric antigen receptor (CAR) T cells. We developed four mRNA CAR T cell constructs using the highly GPC2-specific fully human D3 single chain variable fragment for preclinical testing., Results: We identified high GPC2 expression across multiple pediatric brain tumor types including medulloblastomas, embryonal tumors with multilayered rosettes, other central nervous system embryonal tumors, as well as definable subsets of highly malignant gliomas. We next validated and prioritized CAR configurations using in vitro cytotoxicity assays with GPC2-expressing neuroblastoma cells, where the light-to-heavy single chain variable fragment configurations proved to be superior. We expanded the testing of the two most potent GPC2-directed CAR constructs to GPC2-expressing medulloblastoma and high-grade glioma cell lines, showing significant GPC2-specific cell death in multiple models. Finally, biweekly locoregional delivery of 2-4 million GPC2-directed mRNA CAR T cells induced significant tumor regression in an orthotopic medulloblastoma model and significantly prolonged survival in an aggressive orthotopic thalamic diffuse midline glioma xenograft model. No GPC2-directed CAR T cell related neurologic or systemic toxicity was observed., Conclusion: Taken together, these data show that GPC2 is a highly differentially expressed cell surface protein on multiple malignant pediatric brain tumors that can be targeted safely with local delivery of mRNA CAR T cells, laying the framework for the clinical translation of GPC2-directed immunotherapies for pediatric brain tumors., Competing Interests: Competing interests: TS is currently employed by Spark Therapeutics. KK is currently employed by BioNTech and is an inventor on a patent related to use of nucleoside-modified mRNA. DMB is currently employed by Tmunity Therapeutics. JBF, DMB, JMM, and KRB hold patents for the discovery and development of immunotherapies for cancer, including patents related to glypican 2 (GPC2)-directed immunotherapies. KRB and JMM receive research funding from Tmunity for research on GPC2-directed immunotherapies and JBF, DMB, JMM, and KRB receive royalties from Tmunity for licensing of GPC2-related intellectual property. JMM is a founder of both Tantigen Bio and Hula Therapeutics, focused on cellular therapies for childhood cancers, but neither are working on GPC2-directed therapeutics. All other authors have nothing to disclose., (© Author(s) (or their employer(s)) 2022. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

One of the greatest barriers to curative treatment of neuroblastoma is its frequent metastatic outgrowth prior to diagnosis, especially in cases driven by amplification of the MYCN oncogene. However, only a limited number of regulatory proteins that contribute to this complex MYCN -mediated process have been elucidated. Here we show that the growth arrest-specific 7 ( GAS7 ) gene, located at chromosome band 17p13.1, is preferentially deleted in high-risk MYCN-driven neuroblastoma. GAS7 expression was also suppressed in MYCN -amplified neuroblastoma lacking 17p deletion. GAS7 deficiency led to accelerated metastasis in both zebrafish and mammalian models of neuroblastoma with overexpression or amplification of MYCN . Analysis of expression profiles and the ultrastructure of zebrafish neuroblastoma tumors with MYCN overexpression identified that GAS7 deficiency led to (i) downregulation of genes involved in cell-cell interaction, (ii) loss of contact among tumor cells as critical determinants of accelerated metastasis, and (iii) increased levels of MYCN protein. These results provide the first genetic evidence that GAS7 depletion is a critical early step in the cascade of events culminating in neuroblastoma metastasis in the context of MYCN overexpression. SIGNIFICANCE: Heterozygous deletion or MYCN -mediated repression of GAS7 in neuroblastoma releases an important brake on tumor cell dispersion and migration to distant sites, providing a novel mechanism underlying tumor metastasis in MYCN-driven neuroblastoma. See related commentary by Menard, p. 2815 ., (©2021 American Association for Cancer Research.)

Tumor-associated macrophages (TAMs) play an important role in tumor immunity and comprise of subsets that have distinct phenotype, function, and ontology. Transcriptomic analyses of human medulloblastoma, the most common malignant pediatric brain cancer, showed that medulloblastomas (MBs) with activated sonic hedgehog signaling (SHH-MB) have significantly more TAMs than other MB subtypes. Therefore, we examined MB-associated TAMs by single-cell RNA sequencing of autochthonous murine SHH-MB at steady state and under two distinct treatment modalities: molecular-targeted inhibitor and radiation. Our analyses reveal significant TAM heterogeneity, identify markers of ontologically distinct TAM subsets, and show the impact of brain microenvironment on the differentiation of tumor-infiltrating monocytes. TAM composition undergoes dramatic changes with treatment and differs significantly between molecular-targeted and radiation therapy. We identify an immunosuppressive monocyte-derived TAM subset that emerges with radiation therapy and demonstrate its role in regulating T cell and neutrophil infiltration in MB., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

We report a comprehensive proteogenomics analysis, including whole-genome sequencing, RNA sequencing, and proteomics and phosphoproteomics profiling, of 218 tumors across 7 histological types of childhood brain cancer: low-grade glioma (n = 93), ependymoma (32), high-grade glioma (25), medulloblastoma (22), ganglioglioma (18), craniopharyngioma (16), and atypical teratoid rhabdoid tumor (12). Proteomics data identify common biological themes that span histological boundaries, suggesting that treatments used for one histological type may be applied effectively to other tumors sharing similar proteomics features. Immune landscape characterization reveals diverse tumor microenvironments across and within diagnoses. Proteomics data further reveal functional effects of somatic mutations and copy number variations (CNVs) not evident in transcriptomics data. Kinase-substrate association and co-expression network analysis identify important biological mechanisms of tumorigenesis. This is the first large-scale proteogenomics analysis across traditional histological boundaries to uncover foundational pediatric brain tumor biology and inform rational treatment selection., Competing Interests: Declaration of Interests E.E.S. serves as chief executive officer for Sema4 and has an equity interest in this company., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Background: Gene fusion events are significant sources of somatic variation across adult and pediatric cancers and are some of the most clinically-effective therapeutic targets, yet low consensus of RNA-Seq fusion prediction algorithms makes therapeutic prioritization difficult. In addition, events such as polymerase read-throughs, mis-mapping due to gene homology, and fusions occurring in healthy normal tissue require informed filtering, making it difficult for researchers and clinicians to rapidly discern gene fusions that might be true underlying oncogenic drivers of a tumor and in some cases, appropriate targets for therapy., Results: We developed annoFuse, an R package, and shinyFuse, a companion web application, to annotate, prioritize, and explore biologically-relevant expressed gene fusions, downstream of fusion calling. We validated annoFuse using a random cohort of TCGA RNA-Seq samples (N = 160) and achieved a 96% sensitivity for retention of high-confidence fusions (N = 603). annoFuse uses FusionAnnotator annotations to filter non-oncogenic and/or artifactual fusions. Then, fusions are prioritized if previously reported in TCGA and/or fusions containing gene partners that are known oncogenes, tumor suppressor genes, COSMIC genes, and/or transcription factors. We applied annoFuse to fusion calls from pediatric brain tumor RNA-Seq samples (N = 1028) provided as part of the Open Pediatric Brain Tumor Atlas (OpenPBTA) Project to determine recurrent fusions and recurrently-fused genes within different brain tumor histologies. annoFuse annotates protein domains using the PFAM database, assesses reciprocality, and annotates gene partners for kinase domain retention. As a standard function, reportFuse enables generation of a reproducible R Markdown report to summarize filtered fusions, visualize breakpoints and protein domains by transcript, and plot recurrent fusions within cohorts. Finally, we created shinyFuse for algorithm-agnostic interactive exploration and plotting of gene fusions., Conclusions: annoFuse provides standardized filtering and annotation for gene fusion calls from STAR-Fusion and Arriba by merging, filtering, and prioritizing putative oncogenic fusions across large cancer datasets, as demonstrated here with data from the OpenPBTA project. We are expanding the package to be widely-applicable to other fusion algorithms and expect annoFuse to provide researchers a method for rapidly evaluating, prioritizing, and translating fusion findings in patient tumors.

Medulloblastoma is a highly heterogeneous pediatric brain tumor with five molecular subtypes, Sonic Hedgehog TP53-mutant, Sonic Hedgehog TP53-wildtype, WNT, Group 3, and Group 4, defined by the World Health Organization. The current mechanism for classification into these molecular subtypes is through the use of immunostaining, methylation, and/or genetics. We surveyed the literature and identified a number of RNA-Seq and microarray datasets in order to develop, train, test, and validate a robust classifier to identify medulloblastoma molecular subtypes through the use of transcriptomic profiling data. We have developed a GPL-3 licensed R package and a Shiny Application to enable users to quickly and robustly classify medulloblastoma samples using transcriptomic data. The classifier utilizes a large composite microarray dataset (15 individual datasets), an individual microarray study, and an RNA-Seq dataset, using gene ratios instead of gene expression measures as features for the model. Discriminating features were identified using the limma R package and samples were classified using an unweighted mean of normalized scores. We utilized two training datasets and applied the classifier in 15 separate datasets. We observed a minimum accuracy of 85.71% in the smallest dataset and a maximum of 100% accuracy in four datasets with an overall median accuracy of 97.8% across the 15 datasets, with the majority of misclassification occurring between the heterogeneous Group 3 and Group 4 subtypes. We anticipate this medulloblastoma transcriptomic subtype classifier will be broadly applicable to the cancer research and clinical communities., Competing Interests: The authors have declared that no competing interests exist.

Neuroblastoma is a childhood cancer with heterogeneous clinical outcomes. To comprehensively assess the impact of telomere maintenance mechanism (TMM) on clinical outcomes in high-risk neuroblastoma, we integrated the C-circle assay [a marker for alternative lengthening of telomeres (ALT)], TERT mRNA expression by RNA-sequencing, whole-genome/exome sequencing, and clinical covariates in 134 neuroblastoma patient samples at diagnosis. In addition, we assessed TMM in neuroblastoma cell lines ( n = 104) and patient-derived xenografts ( n = 28). ALT was identified in 23.4% of high-risk neuroblastoma tumors and genomic alterations in ATRX were detected in 60% of ALT tumors; 40% of ALT tumors lacked genomic alterations in known ALT-associated genes. Patients with high-risk neuroblastoma were classified into three subgroups (TERT-high, ALT + , and TERT-low/non-ALT) based on presence of C-circles and TERT mRNA expression (above or below median TERT expression). Event-free survival was similar among TERT-high, ALT + , or TERT-low/non-ALT patients. However, overall survival (OS) for TERT-low/non-ALT patients was significantly higher relative to TERT-high or ALT patients (log-rank test; P < 0.01) independent of current clinical and molecular prognostic markers. Consistent with the observed higher OS in patients with TERT-low/non-ALT tumors, continuous shortening of telomeres and decreasing viability occurred in low TERT-expressing, non-ALT patient-derived high-risk neuroblastoma cell lines. These findings demonstrate that assaying TMM with TERT mRNA expression and C-circles provides precise stratification of high-risk neuroblastoma into three subgroups with substantially different OS: a previously undescribed TERT-low/non-ALT cohort with superior OS (even after relapse) and two cohorts of patients with poor survival that have distinct molecular therapeutic targets. SIGNIFICANCE: These findings assess telomere maintenance mechanisms with TERT mRNA and the ALT DNA biomarker C-circles to stratify neuroblastoma into three groups, with distinct overall survival independent of currently used clinical risk classifiers., (©2020 American Association for Cancer Research.)

Background: The treatment of high-risk neuroblastoma continues to present a formidable challenge to pediatric oncology. Previous studies have shown that Bromodomain and extraterminal (BET) inhibitors can inhibit MYCN expression and suppress MYCN-amplified neuroblastoma in vivo. Furthermore, alterations within RAS-MAPK (mitogen-activated protein kinase) signaling play significant roles in neuroblastoma initiation, maintenance, and relapse, and mitogen-activated extracellular signal-regulated kinase (MEK) inhibitors demonstrate efficacy in subsets of neuroblastoma preclinical models. Finally, hyperactivation of RAS-MAPK signaling has been shown to promote resistance to BET inhibitors. Therefore, we examined the antitumor efficacy of combined BET/MEK inhibition utilizing I-BET726 or I-BET762 and trametinib in high-risk neuroblastoma., Procedure: Utilizing a panel of genomically annotated neuroblastoma cell line models, we investigated the in vitro effects of combined BET/MEK inhibition on cell proliferation and apoptosis. Furthermore, we evaluated the effects of combined inhibition in neuroblastoma xenograft models., Results: Combined BET and MEK inhibition demonstrated synergistic effects on the growth and survival of a large panel of neuroblastoma cell lines through augmentation of apoptosis. A combination therapy slowed tumor growth in a non-MYCN-amplified, NRAS-mutated neuroblastoma xenograft model, but had no efficacy in an MYCN-amplified model harboring a loss-of-function mutation in NF1., Conclusions: Combinatorial BET and MEK inhibition was synergistic in the vast majority of neuroblastoma cell lines in the in vitro setting but showed limited antitumor activity in vivo. Collectively, these data do not support clinical development of this combination in high-risk neuroblastoma., (© 2020 Wiley Periodicals, Inc.)

Neuroblastoma is an aggressive pediatric malignancy of the neural crest with suboptimal cure rates and a striking predilection for widespread metastases, underscoring the need to identify novel therapeutic vulnerabilities. We recently identified the RNA binding protein LIN28B as a driver in high-risk neuroblastoma and demonstrated it promotes oncogenic cell proliferation by coordinating a RAN-Aurora kinase A network. Here, we demonstrate that LIN28B influences another key hallmark of cancer, metastatic dissemination. Using a murine xenograft model of neuroblastoma dissemination, we show that LIN28B promotes metastasis. We demonstrate that this is in part due to the effects of LIN28B on self-renewal and migration, providing an understanding of how LIN28B shapes the metastatic phenotype. Our studies reveal that the let-7 family, which LIN28B inhibits, decreases self-renewal and migration. Next, we identify PDZ Binding Kinase (PBK) as a novel LIN28B target. PBK is a serine/threonine kinase that promotes the proliferation and self-renewal of neural stem cells and serves as an oncogenic driver in multiple aggressive malignancies. We demonstrate that PBK is both a novel direct target of let-7i and that MYCN regulates PBK expression, thus elucidating two oncogenic drivers that converge on PBK. Functionally, PBK promotes self-renewal and migration, phenocopying LIN28B. Taken together, our findings define a role for LIN28B in neuroblastoma metastasis and define the targetable kinase PBK as a potential novel vulnerability in metastatic neuroblastoma., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Understanding the aberrant transcriptional landscape of neuroblastoma is necessary to provide insight to the underlying influences of the initiation, progression and persistence of this developmental cancer. Here, we present chromatin immunoprecipitation sequencing (ChIP-Seq) data for the oncogenic transcription factors, MYCN and MYC, as well as regulatory histone marks H3K4me1, H3K4me3, H3K27Ac, and H3K27me3 in ten commonly used human neuroblastoma-derived cell line models. In addition, for all of the profiled cell lines we provide ATAC-Seq as a measure of open chromatin. We validate specificity of global MYCN occupancy in MYCN amplified cell lines and functional redundancy of MYC occupancy in MYCN non-amplified cell lines. Finally, we show with H3K27Ac ChIP-Seq that these cell lines retain expression of key neuroblastoma super-enhancers (SE). We anticipate this dataset, coupled with available transcriptomic profiling on the same cell lines, will enable the discovery of novel gene regulatory mechanisms in neuroblastoma.

We developed a computational pipeline designed to use RNA sequencing ( n = 136) and gene expression profiling ( n = 250) data from neuroblastoma tumors to identify cell surface proteins predicted to be highly expressed in MYCN amplified neuroblastomas and with little or no expression in normal human tissues. We then performed ChIP-seq in the MYCN amplified cell lines KELLY, NB-1643, and NGP to identify gene promoters that are occupied by MYCN protein to define the intersection with the differentially-expressed gene list. We initially identified 116 putative immunotherapy targets with predicted transmembrane domains, with the most significant differentially-expressed of these being the calmodulin kinase-like vesicle-associated gene (CAMKV, p = 2 × 10 -6 ). CAMKV encodes a protein that binds calmodulin in the presence of calcium, but lacks the kinase activity of other calmodulin kinase family members. We confirmed that CAMKV is selectively expressed in 7/7 MYCN amplified neuroblastoma cell lines and showed that the transcription of CAMKV is directly controlled by MYCN. From membrane fractionation and immunohistochemistry, we verified that CAMKV is membranous in MYCN amplified neuroblastoma cell lines and patient-derived xenografts. Finally, immunohistochemistry showed that CAMKV is not expressed on normal tissues outside of the central nervous system. Together, these data demonstrate that CAMKV is a differentially-expressed cell surface protein that is transcriptionally regulated by MYCN, making it a candidate for targeting with antibodies or antibody-drug conjugates that do not cross the blood brain barrier., (Copyright © 2020 Sussman, Rokita, Huang, Raman, Rathi, Martinez, Bosse, Lane, Hart, Bhatti, Pawel and Maris.)

Relapsed neuroblastomas are enriched with activating mutations of the RAS-MAPK signaling pathway. The MEK1/2 inhibitor trametinib delays tumor growth but does not sustain regression in neuroblastoma preclinical models. Recent studies have implicated the Hippo pathway transcriptional coactivator protein YAP1 as an additional driver of relapsed neuroblastomas, as well as a mediator of trametinib resistance in other cancers. Here, we used a highly annotated set of high-risk neuroblastoma cellular models to modulate YAP1 expression and RAS pathway activation to test whether increased YAP1 transcriptional activity is a mechanism of MEK1/2 inhibition resistance in RAS-driven neuroblastomas. In NLF (biallelic NF1 inactivation) and SK-N-AS (NRAS Q61K) cell lines, trametinib caused a near-complete translocation of YAP1 protein into the nucleus. YAP1 depletion sensitized neuroblastoma cells to trametinib, while overexpression of constitutively active YAP1 protein induced trametinib resistance. Mechanistically, significant enhancement of G 1 -S cell-cycle arrest, mediated by depletion of MYC/MYCN and E2F transcriptional output, sensitized RAS-driven neuroblastomas to trametinib following YAP1 deletion. These findings underscore the importance of YAP activity in response to trametinib in RAS-driven neuroblastomas, as well as the potential for targeting YAP in a trametinib combination. SIGNIFICANCE: High-risk neuroblastomas with hyperactivated RAS signaling escape the selective pressure of MEK inhibition via YAP1-mediated transcriptional reprogramming and may be sensitive to combination therapies targeting both YAP1 and MEK., (©2019 American Association for Cancer Research.)

A heritable polymorphism within regulatory sequences of the LMO1 gene is associated with its elevated expression and increased susceptibility to develop neuroblastoma, but the oncogenic pathways downstream of the LMO1 transcriptional co-regulatory protein are unknown. Our ChIP-seq and RNA-seq analyses reveal that a key gene directly regulated by LMO1 and MYCN is ASCL1, which encodes a basic helix-loop-helix transcription factor. Regulatory elements controlling ASCL1 expression are bound by LMO1, MYCN and the transcription factors GATA3, HAND2, PHOX2B, TBX2 and ISL1-all members of the adrenergic (ADRN) neuroblastoma core regulatory circuitry (CRC). ASCL1 is required for neuroblastoma cell growth and arrest of differentiation. ASCL1 and LMO1 directly regulate the expression of CRC genes, indicating that ASCL1 is a member and LMO1 is a coregulator of the ADRN neuroblastoma CRC.

Background: Unlike some adult cancers, most pediatric cancers are considered immunologically cold and generally less responsive to immunotherapy. While immunotherapy has already been incorporated into standard of care treatment for pediatric patients with high-risk neuroblastoma, overall survival remains poor. In a mouse melanoma model, we found that radiation and tumor-specific immunocytokine generate an in situ vaccination response in syngeneic mice bearing large tumors. Here, we tested whether a novel immunotherapeutic approach utilizing radiation and immunocytokine together with innate immune stimulation could generate a potent antitumor response with immunologic memory against syngeneic murine neuroblastoma., Methods: Mice bearing disialoganglioside (GD2)-expressing neuroblastoma tumors (either NXS2 or 9464D-GD2) were treated with radiation and immunotherapy (including anti-GD2 immunocytokine with or without anti-CTLA-4, CpG and anti-CD40 monoclonal antibody). Tumor growth, animal survival and immune cell infiltrate were analyzed in the tumor microenvironment in response to various treatment regimens., Results: NXS2 had a moderate tumor mutation burden (TMB) while N-MYC driven 9464D-GD2 had a low TMB, therefore the latter served as a better model for high-risk neuroblastoma (an immunologically cold tumor). Radiation and immunocytokine induced a potent in situ vaccination response against NXS2 tumors, but not in the 9464D-GD2 tumor model. Addition of checkpoint blockade with anti-CTLA-4 was not effective alone against 9464D-GD2 tumors; inclusion of CpG and anti-CD40 achieved a potent antitumor response with decreased T regulatory cells within the tumors and induction of immunologic memory., Conclusions: These data suggest that a combined innate and adaptive immunotherapeutic approach can be effective against immunologically cold syngeneic murine neuroblastoma. Further testing is needed to determine how these concepts might translate into development of more effective immunotherapeutic approaches for the treatment of clinically high-risk neuroblastoma.

Accelerating cures for children with cancer remains an immediate challenge as a result of extensive oncogenic heterogeneity between and within histologies, distinct molecular mechanisms evolving between diagnosis and relapsed disease, and limited therapeutic options. To systematically prioritize and rationally test novel agents in preclinical murine models, researchers within the Pediatric Preclinical Testing Consortium are continuously developing patient-derived xenografts (PDXs)-many of which are refractory to current standard-of-care treatments-from high-risk childhood cancers. Here, we genomically characterize 261 PDX models from 37 unique pediatric cancers; demonstrate faithful recapitulation of histologies and subtypes; and refine our understanding of relapsed disease. In addition, we use expression signatures to classify tumors for TP53 and NF1 pathway inactivation. We anticipate that these data will serve as a resource for pediatric oncology drug development and will guide rational clinical trial design for children with cancer., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Neuroblastoma is a cancer of the developing sympathetic nervous system. It is diagnosed in 600-700 children per year in the United States and accounts for 12% of pediatric cancer deaths. Despite recent advances in our understanding of this malignancy's complex genetic architecture, the contribution of rare germline variants remains undefined. Here, we conducted a genome-wide analysis of large (>500 kb), rare (<1%) germline copy number variants (CNVs) in two independent, multi-ethnic cohorts totaling 5,585 children with neuroblastoma and 23,505 cancer-free control children. We identified a 550-kb deletion on chromosome 16p11.2 significantly enriched in neuroblastoma cases (0.39% of cases and 0.03% of controls; p = 3.34 × 10 -9 ). Notably, this CNV corresponds to a known microdeletion syndrome that affects approximately one in 3,000 children and confers risk for diverse developmental phenotypes including autism spectrum disorder and other neurodevelopmental disorders. The CNV had a substantial impact on neuroblastoma risk, with an odds ratio of 13.9 (95% confidence interval = 5.8-33.4). The association remained significant when we restricted our analysis to individuals of European ancestry in order to mitigate potential confounding by population stratification (0.42% of cases and 0.03% of controls; p = 4.10 × 10 -8 ). We used whole-genome sequencing (WGS) to validate the deletion in paired germline and tumor DNA from 12 cases. Finally, WGS of four parent-child trios revealed that the deletion primarily arose de novo without maternal or paternal bias. This finding expands the clinical phenotypes associated with 16p11.2 microdeletion syndrome to include cancer, and it suggests that disruption of the 16p11.2 region may dysregulate neurodevelopmental pathways that influence both neurological phenotypes and neuroblastoma., (Copyright © 2019 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Single-cell gene expression analyses of mammalian tissues have uncovered profound stage-specific molecular regulatory phenomena that have changed the understanding of unique cell types and signaling pathways critical for lineage determination, morphogenesis, and growth. We discuss here the case for a Pediatric Cell Atlas as part of the Human Cell Atlas consortium to provide single-cell profiles and spatial characterization of gene expression across human tissues and organs. Such data will complement adult and developmentally focused HCA projects to provide a rich cytogenomic framework for understanding not only pediatric health and disease but also environmental and genetic impacts across the human lifespan., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)