AutoGVP: a dockerized workflow integrating ClinVar and InterVar germline sequence variant classification

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.

Everolimus for Children With Recurrent or Progressive Low-Grade Glioma: Results From the Phase II PNOC001 Trial

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/m2 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.

BARD1 germline variants induce haploinsufficiency and DNA repair defects in neuroblastoma

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.

A proteogenomic surfaceome study identifies DLK1 as an immunotherapeutic target in neuroblastoma

Cancer immunotherapies have produced 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 specimens along with normal tissues to identify biologically relevant cell surface proteins that can serve as immunotherapeutic targets for neuroblastoma, an often-fatal childhood cancer of the developing nervous system. We apply this approach to human-derived cell lines (N=9) and cell/patient-derived xenograft (N=12) models of neuroblastoma. Plasma membrane-enriched mass spectrometry identified 1,461 cell surface proteins in cell lines and 1,401 in xenograft models, respectively. Additional proteogenomic analyses revealed 60 high-confidence candidate immunotherapeutic targets and we prioritized Delta-like canonical notch ligand 1 (DLK1) for further study. High expression of DLK1 directly correlated with the presence of a super-enhancer spanning the DLK1 locus. Robust cell surface expression of DLK1 was validated by immunofluorescence, flow cytometry, and immunohistochemistry. Short hairpin RNA mediated silencing of DLK1 in neuroblastoma cells resulted in increased cellular differentiation. ADCT-701, a DLK1-targeting antibody-drug conjugate (ADC), showed potent and specific cytotoxicity in DLK1-expressing neuroblastoma xenograft models. Moreover, DLK1 is highly expressed in several adult cancer types, including adrenocortical carcinoma (ACC), pheochromocytoma/paraganglioma (PCPG), hepatoblastoma, and small cell lung cancer (SCLC), suggesting potential clinical benefit beyond neuroblastoma. Taken together, our study demonstrates the utility of comprehensive cancer surfaceome characterization and credentials DLK1 as an immunotherapeutic target.

Highlights

Plasma membrane enriched proteomics defines surfaceome of neuroblastomaMulti-omic data integration prioritizes DLK1 as a candidate immunotherapeutic target in neuroblastoma and other cancersDLK1 expression is driven by a super-enhancer DLK1 silencing in neuroblastoma cells results in cellular differentiation ADCT-701, a DLK1-targeting antibody-drug conjugate, shows potent and specific cytotoxicity in DLK1-expressing neuroblastoma preclinical models.

AutoGVP: a dockerized workflow integrating ClinVar and InterVar germline sequence variant classification

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.

OpenPBTA: The Open Pediatric Brain Tumor Atlas

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.

Alternative lengthening of telomeres (ALT) in pediatric high-grade gliomas can occur without ATRX mutation and is enriched in patients with pathogenic germline mismatch repair (MMR) variants

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.

Abstract 3566: Expansion of the Pediatric Brain Tumor Atlas: Children's Brain Tumor Network, Kids First Data Resource and Childhood Cancer Data Initiative Open Science effort

Pediatric central nervous system (CNS) cancers are the leading disease-related cause of death in children and there is urgent need for curative therapeutic strategies for these tumors. To address the urgency, Children’s Brain Tumor Network (CBTN) has advanced an open science model to accelerate the research discovery for pediatric brain tumors. In first phase of Open Pediatric Brain Tumor Atlas (OpenPBTA) effort CBTN together with Pacific Pediatric Neuro-Oncology Consortium (PNOC) with support of Gabriella Miller Kids First Data Resource Center (KFDRC) created and comprehensively characterized over 1000 clinically annotated pediatric brain tumors.

In the second phase of the OpenPBTA effort, through resource awards and collaboration across KFDRC, the NCI Childhood Cancer Data Initiative (CCDI), NCI’s Clinical Proteomic Tumor Analysis Consortium (CPTAC), NCI Center for Cancer Research and additional partnered institutions and foundations, CBTN has expanded OpenPBTA to support high throughput molecular characterization for an additional 1900 pediatric brain tumor patients and their families. This includes the processing and characterization of over 8000 specimens across &gt;50 brain tumor diagnoses. The cohort expansion builds on &gt;1000 previously characterized samples with a portfolio of multimodal data including whole genome sequencing, RNA sequencing, miRNA sequencing, methylation sequencing, proteomics, lipidomics and/or metabolomics. Molecular data is linked to patient longitudinal clinical data, imaging data (MRIs and radiology reports), histology slide images, and pathology reports. To inform novel discovery and clinical implementation of genomic approaches for diagnostic/therapeutic purposes, the data is deposited the cloud-based research environment of the NCI’s CCDI and the KFDRC to provide near real-time integration, dissemination, processing, and sharing of associated petabyte-scale harmonized data. The approach leverages the DRC platform’s cloud-based computational environment in CAVATICA. Processed annotations are facilitated via CAVATICA-enabled shareable pipelines and can be explored through PedcBioPortal, a data visualization/analysis application further integrating additional public and deposited datasets. This expansion phase of OpenPBTA is released with no embargo period and provides one of the largest deeply characterized cohorts of pediatric brain tumor samples and associated clinical data for &gt;3000 pediatric brain tumor patients. CBTN’s open-science, rapid-release model aims to advance novel biomarkers and therapeutic exploratory research, supporting new clinical trial development and accelerated discovery on behalf of changing the outcome for kids with brain tumors.

Citation Format: Mateusz P. Koptyra, Komal Rahti, Yuankun Zhu, Bailey Farrow, Daniel Miller, Adam Kraya, Yiran Guo, Peter Madsen, Nicholas Van Kuren, Xiaoyan Huang, Miguel A. Brown, Jennifer L. Mason, Meen Chul Kim, Allison P. Heath, Brian M. Ennis, Bo Zhang, Jena V. Lilly, Jo Lynne Rokita, Christopher Friedman, Ximena P. Cuellar, Catherine A. Sullivan, Noel Coleman, Trang Duros, Thinh Q. Nguyen, Emmett C. Drake, Zeinab Helili, Beth A. Frenkel, Gerri R. Trooskin, Ariana Familiar, Karthik Viswanathan, Christopher M. Beck, Madison L. Hollawell, Valerie P. Baubet, Cassie Kline, Mariarita Santi, Tatiana S. Patton, Stephanie Stefankiewicz, Arya Kamnaa, Ryan A. Velasco, Dani Cardona, Phillip J. Storm, Adam C. Resnick, o/b/o Children's Brain Tumor Network. Expansion of the Pediatric Brain Tumor Atlas: Children's Brain Tumor Network, Kids First Data Resource and Childhood Cancer Data Initiative Open Science effort. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3566.

BARD1 germline variants induce haploinsufficiency and DNA repair defects in neuroblastoma

Importance

High-risk neuroblastoma is a complex genetic disease that is lethal in 50% of patients despite intense multimodal therapy. Our genome-wide association study (GWAS) identified single-nucleotide polymorphisms (SNPs) within the BARD1 gene showing the most significant enrichment in neuroblastoma patients, and also discovered pathogenic (P) or likely pathogenic (LP) rare germline loss-of-function variants in this gene. The functional implications of these findings remain poorly understood.

Objective

To define the functional relevance of BARD1 germline variation in children with neuroblastoma.

Design

We correlated BARD1 genotype with BARD1 expression in normal and tumor cells and the cellular burden of DNA damage in tumors. To validate the functional consequences of rare germline P-LP BARD1 variants, we generated isogenic cellular models harboring heterozygous BARD1 loss-of-function (LOF) variants and conducted multiple complementary assays to measure the efficiency of DNA repair.

Setting

(N/A).

Participants

(N/A).

Interventions/Exposures

(N/A).

Main Outcomes and Measures

BARD1 expression, efficiency of DNA repair, and genome-wide burden of DNA damage in neuroblastoma tumors and cellular models harboring disease-associated BARD1 germline variants.

Results

Both common and rare neuroblastoma associated BARD1 germline variants were significantly associated with lower levels of BARD1 mRNA and an increased burden of DNA damage. Using neuroblastoma cellular models engineered to harbor disease-associated heterozygous BARD1 LOF variants, we functionally validated this association with inefficient DNA repair. These BARD1 LOF variant isogenic models exhibited reduced efficiency in repairing Cas9-induced DNA damage, ineffective RAD51 focus formation at DNA doublestrand break sites, and enhanced sensitivity to cisplatin and poly-ADP ribose polymerase (PARP) inhibition.

Conclusions and Relevance

Considering that at least 1 in 10 children diagnosed with cancer carry a predicted pathogenic mutation in a cancer predisposition gene, it is critically important to understand their functional relevance. Here, 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, and these findings may also extend to other cancers harboring germline variants in genes essential for DNA damage repair.

Key Points

Question: How do neuroblastoma patient BRCA1-associated RING domain 1 ( BARD1 ) germline variants impact DNA repair? Findings: Neuroblastoma-associated germline BARD1 variants disrupt DNA repair fidelity. Common risk variants correlate with decreased BARD1 expression and increased DNA double-strand breaks in neuroblastoma tumors and rare heterozygous loss-of-function variants induce BARD1 haploinsufficiency, resulting in defective DNA repair and genomic instability in neuroblastoma cellular models. Meaning: Germline variation in BARD1 contributes to neuroblastoma pathogenesis via dysregulation of critical cellular DNA repair functions, with implications for neuroblastoma treatment, risk stratification, and cancer predisposition.

The children's brain tumor network (CBTN) - Accelerating research in pediatric central nervous system tumors through collaboration and open science

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.

Upfront Biology-Guided Therapy in Diffuse Intrinsic Pontine Glioma: Therapeutic, Molecular, and Biomarker Outcomes from PNOC003

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.

Development of GPC2-directed chimeric antigen receptors using mRNA for pediatric brain tumors

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.

Abstract 2725: Molecular mechanisms of aberrant splicing in pediatric high grade gliomas

Pediatric brain cancer is the number one cause of disease-related death in children. While much of the somatic mutation variation in these tumors have been profiled, there have been few studies exploring their transcriptional splicing programs. Here, we characterize aberrant alternative splicing across pediatric brain tumors. We identified three splicing-based clusters that spanned across seven broad histological tumor classifications. Patients in clusters 1 and 2 had significantly lower overall survival (p &lt; 0.001) than those in cluster 3. High-grade gliomas (HGGs) displayed the greatest splicing heterogeneity, spanning all three clusters. We focused on midline HGGs (n = 84), identifying 19,275 recurrent and significant (20% change from control, P &lt; 0.05, FDR &lt; 0.05) aberrant splicing events in 8,587 genes compared to non-diseased brainstem controls. We also identified several splicing factor genes, including SRSF11, that were differentially expressed (p &lt; 0.05), suggesting a potential tumor-specific novel mechanism of splicing dysregulation. Strikingly, 27% (5,157) of these splicing events resulted in either a gain or loss of a functional site in 3,294 genes: disulfide bonding (n = 391), localization signaling (n = 216), modifications, such as phosphorylation sites (n = 1,670), and other, such as ion-binding sites (n = 4,245). We prioritized splice variants affecting targetable kinases and found that mRNAs encoding CDC-like kinase 1 (CLK1), a known modulator of the SR family of master splicing regulators, exhibit significantly increased exon 4 inclusion. This leads to a gain of two known phosphorylation sites in CLK1, increased protein expression, and hyper-phosphorylation of SR proteins, inducing cooperative binding of SR proteins to RNA. To assess the impact of this CLK1 event, we compared splicing in tumors with high exon 4 inclusion (n = 5) to those with low exon 4 inclusion (n = 5) and discovered 862 differential splicing events (DSE). Samples with high CLK1 exon 4 inclusion had significantly lower percent spliced in events, suggesting that increased activity of CLK1 resulting from the gain of the exon 4 promotes overall exon skipping. Additionally, more than 80% of the DSE (n = 695), including those in known brain tumor oncogenes such as CDC27 and E2H2, were recurrently mis-spliced in midline HGG tumors, suggesting that CLK1 splicing dysregulation leads to a global aberrant splicing phenotype in HGGs. Finally, we determined a subset of events significantly associated with prolonged survival, suggesting that stratification of patients based on splicing phenotypes can be utilized as a prognostic indicator. In summary, we demonstrate aberrant splicing in pediatric HGGs as an alternative mechanism driving tumorigenesis and future work will focus on urgently needed therapeutic targeting of these pathways toward our overall goal of improving survival outcomes in children with brain tumors.

Citation Format: Ammar S. Naqvi, Brian Ennis, Run Jin, Karina L. Conkrite, Komal S. Rathi, Adam A. Kraya, Poonam Sonawane, Phillip B. Storm, Adam C. Resnick, Jo Lynne Rokita. Molecular mechanisms of aberrant splicing in pediatric high grade gliomas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2725.

HGG-54. CLK1 aberrant splicing in pediatric high-grade gliomas disrupts key oncogenic transcriptional programs

While much of the somatic coding variation underlying the oncogenic transformation of pediatric high-grade gliomas (HGGs) has been profiled, transcriptional splicing programs of these tumors remain under-explored. Here, we characterize aberrant alternative splicing in pediatric midline HGGs (n = 84). We identified 19,275 recurrent and significant (20% change from control, P &lt; 0.05, FDR &lt; 0.05) aberrant splicing events in 8,587 genes compared to non-diseased brainstem controls. Of those, 27% (n = 5,157) resulted in either a gain or loss of a known protein functional site within 3,294 genes. We prioritized splice variants affecting targetable kinases and found that mRNAs encoding CDC-like kinase 1 (CLK1), a known modulator of master splicing regulators, exhibit significantly increased exon 4 inclusion in midline HGGs. This leads to a gain of two known phosphorylation sites in CLK1, increased CLK1 protein expression and hyper-phosphorylation of Serine-rich splicing factors. To assess the impact of this event, we performed differential splicing and expression analyses, comparing tumors with the highest (n= 5) and lowest (n = 5) exon 4 inclusion. We discovered 3,037 genes to be differentially up-regulated in high exon 4 inclusion tumors with an enrichment of cancer-related pathways, including DNA repair, mitotic spindle, myogenesis and EMT. We next integrated these gene signatures with protein-protein interaction networks of kinase and transcription factors and show that increased CLK1 exon 4 inclusion disrupts critical regulatory networks, such as those involving FOXM1, which is implicated in cell cycle and proliferation processes. In summary, we describe aberrant splicing in pediatric HGGs as an additional mechanism that could drive tumorigenesis. Future work will focus on molecular validation and therapeutic targeting of CLK1 in available HGG models. Characterizing tumor-specific splicing variation has the potential to open new therapeutic strategies and understand mechanisms of treatment resistance in children with central nervous system tumors.

Abstract 3024: Targeting the RNA binding protein LIN28B in Group 3 medulloblastoma decreases proliferation and promotes apoptosis

Medulloblastoma (MB) is the most common pediatric malignant brain tumor and is currently divided into WNT, SHH, Group 3 and Group 4 subtypes. Even with multimodal chemotherapy, radiotherapy and surgery, many children with Group 3 MBs do not survive. While the molecular aberrations underlying WNT- and SHH-driven MBs are relatively well understood, the oncogenic drivers that lead to Group 3/4 MBs are poorly defined, limiting therapeutic progress. In addition to genetic mutations and alterations, cancers display dysregulated transcription and translation. RNA-binding proteins (RBPs) play key roles in both transcription and translation, and a subset of RBPs are differentially expressed in many different cancers. Indeed, we have previously demonstrated an oncogenic role for the RBP LIN28B in neuroblastoma and it is known to be upregulated in Wilms tumor, hepatoblastoma, germ cell tumors, leukemia among others. LIN28B is a key regulator of let-7 family miRNAs, which in turn inhibit LIN28B and other oncogenes. We hypothesize that LIN28B plays an important role in Group 3 MB and that a better understanding of LIN28B and LIN28B-driven networks will reveal novel therapeutic vulnerabilities. In support of our hypothesis we find that among the four subtypes, LIN28B levels are highest in Group 3 MB, and that overexpression is associated with significantly worse survival. Down-regulation of LIN28B results in significant reduction in cell proliferation by CellTiter-Glo and increased apoptosis by Caspase-Glo (as well as induction of cleaved PARP on immunoblots). In contrast overexpression of LIN28B increases Group 3 cell proliferation and tumor sphere formation. In addition we find that PDZ-binding kinase (PBK) a downstream target of LIN28B is downregulated when LIN28B is depleted. PBK knock down also leads to decreased proliferation of Group 3 MB cells. Finally, in order to robustly define the signaling networks downstream from LIN28B that are involved in Group 3 MB metastasis, we have performed who transcriptome RNA-seq profiling of two group 3 cell lines following LIN28B depletion and plan to interrogate a subset of these based on expression change and functional relevance to LIN28B-mediated Group 3 MB metastasis. This work will help define the role for LIN28B in Group 3 MB aggressiveness and pave the way for similar studies in other cancers.

Citation Format: Shubin W. Shahab, Jo Lynne Rokita, Kyle Juraschka, Sachin Kumar, Michael Taylor, Robert W. Schnepp, Tobey J. MacDonald, Anna M. Kenney. Targeting the RNA binding protein LIN28B in Group 3 medulloblastoma decreases proliferation and promotes apoptosis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 3024.

OMIC-14. OPENPBTA: AN OPEN PEDIATRIC BRAIN TUMOR ATLAS

Pediatric brain tumors comprise a heterogeneous molecular and histological landscape that challenges most current precision-medicine approaches. While recent large-scale efforts to molecularly characterize distinct histological entities have dramatically advanced the field’s capacity to classify and further define molecular subtypes, developing therapeutic and less toxic molecularly-defined clinical approaches remains a challenge. To define new approaches to meet these challenges and advance scalable, shared biospecimen- and data-resources for pediatric brain tumors, the Children’s Brain Tumor Network and Pacific Pediatric Neuro-Oncology Consortium, in partnership with the Alex’s Lemonade Stand Foundation Childhood Cancer Data Lab, launched OpenPBTA, a global open science Pediatric Brain Tumor Atlas initiative to comprehensively define the molecular landscape of pediatric brain tumors. The initiative contains multi-modal analyses of research- and clinical-trial based DNA and RNA sequences from nearly 1,000 subjects (with 1,256 tumors) along with their longitudinal clinical data. The OpenPBTA’s open science framework for analysis tests the capacity of crowd-sourced collaborative architectures to advance more rapid, iterative and integrated discovery of the underlying mechanisms of disease across pediatric brain and spinal cord tumors. Since the launch of the project, OpenPBTA has collaboratively created reproducible workflows for integrated consensus SNV, CNV, and fusion calling, enabled RNA-Seq-based classification of medulloblastoma subtypes, and more than 25 additional DNA- and RNA-based analyses. The open-science platform and associated datasets and processed results provide a continuously updated, global view of the integrated cross-disease molecular landscape of pediatric brain tumors. Such biospecimen- and clinically-linked scalable data resources provide unprecedented collaborative opportunities for precision-based, personalized therapeutic discovery and drug development with the upcoming further integration of proteomic sample data (N >300) and drug response datasets, additionally diversifying the multimodal discovery potential of crowd-sourced approaches for accelerated impact for children with brain tumors.

Macrophages in SHH subgroup medulloblastoma display dynamic heterogeneity that varies with treatment modality

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.

GAS7 Deficiency Promotes Metastasis in MYCN-Driven Neuroblastoma

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.

Integrated Proteogenomic Characterization across Major Histological Types of Pediatric Brain Cancer

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.

annoFuse: an R Package to annotate, prioritize, and interactively explore putative oncogenic RNA fusions

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.

EPID-14. GABRIELLA MILLER KIDS FIRST DATA RESOURCE CENTER: COLLABORATIVE PLATFORMS FOR ACCELERATING RESEARCH IN PEDIATRIC CANCERS & STRUCTURAL BIRTH DEFECTS

Since launching to the public in September 2018, the Gabriella Miller Kids First Data Resource Center (DRC) has made an increasing number of pediatric genomic studies available to the research community. Currently, 1.3 PBs of genomic and clinical data drawn from 12,000 participants are available across a variety of pediatric cancer and structural birth defect studies. The DRC has architected a secure, cloud-based platform with over 1,300 users that supports the ability of researchers to not only find, access, and reuse data, but also integrate, collaborate, and analyze data quickly at scale. Users can use integrations with platforms such as Cavatica for bioinformatics workflows and PedcBioPortal for cancer genomic visualizations. Additionally, a set of framework services, powered by Gen3, provide a foundation for interoperability with other large-scale data sources, platforms, and a growing ecosystem of analysis and visualization applications. These integrations allow users to search across both TARGET and Kids First clinical data in one location while allowing data governance to be maintained by the original approvers. The new “explore data” feature allows users to search across all studies in order to identify virtual cohorts. Within the portal, these cohorts can be saved and shared with collaborators for iterative refinement and analysis. With appropriate approvals, the associated genomic data can be accessed and analyzed seamlessly in Cavatica or other platforms with interoperable framework services. Additionally, gene searching capabilities will be available in 2020. Data is free to download and cloud credits are available for analysis support.

DIPG-31. MOLECULAR MECHANISMS AND FUNCTIONAL IMPACT OF ABERRANT SPLICING IN DIFFUSE MIDLINE GLIOMAS

Fewer than 1% of children diagnosed with diffuse-midline glioma (DMG) survive for more than 5 years, because no effective therapies exist for these patients. Here, we sought to identify and characterize mechanisms of aberrant splicing (AS) in primary DMG tumors. We observed transcriptome-wide AS (9,805 differential splicing variations in 4,734 genes), and identified a DMG-specific splicing signature, that included known cancer genes. We hypothesize that AS of cancer genes play a role in DMG tumor formation. Assessing whether splicing factor dysregulation impacted known cancer transcripts, we discovered several splicing factors, including SRRM4, SRRM3 and RBFOX3 to be down-regulated in DMG. Additionally, we found an enrichment of binding motifs for these proteins within flanking regions of these mis-spliced exons. We also observed recurrent significant exon inclusion in tumor suppressor SMARCA4, an integral member of the SWI/SNF family of proteins involved in chromatin remodeling. Further, we identified AS in 16 of the 27 members of the SWI/SNF complex, including increased skipping of exon 7 in DPF2, representing a complete mRNA transcript switch in DMG. Since SRRM4, SRRM3 and RBFOX3 are known regulators for neural-specific microexons, we focused on microexon splicing changes, hypothesizing that these regulators may be driving microexon mis-splicing in these tumors. We identified 245 known microexons lost or gained in DMG. Moreover, a quarter of which were observed in known cancer genes, with the most frequent splice event causing gain of a clathrin-binding site in the tumor suppressor BIN1 with a concurrent loss of an out-of-frame microexon in the oncogene BAK1, presumably activating it. Altogether, our results suggest that aberrant splicing may be an alternative mechanism driving DMG tumorigenesis and we are currently molecularly validating a subset of these events with the overall goal of identifying novel therapeutic targets for DMG tumors.

A transcriptome-based classifier to determine molecular subtypes in medulloblastoma

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.

Abstract 445: Integrated proteogenomic characterization across seven histological types of pediatric brain tumors

We performed a comprehensive proteogenomic analysis across seven major types of childhood brain tumors for a deeper understanding of their functional biology. Whole genome seq, RNAseq, quantitative proteomic and phosphoproteomic profiling were performed on 219 fresh frozen tumor samples representing the histologic diagnoses of: low grade astrocytoma (93), ependymoma (32), high grade astrocytoma (26), medulloblastoma (22), ganglioglioma (18), craniopharyngioma (16) and atypical teratoid rhabdoid tumor (12). Characterization of the tumor microenvironment through multi-omics based deconvolution analyses revealed 5 distinct tumor clusters associated with different populations of infiltrating immune cells: Cold-medullo, Cold-mixed, Epithelial, Neuronal and Hot. The two cold-tumor clusters have the lowest immune cell infiltration, one characterized by the enrichment of medulloblastoma tumors; while the other is a mixture of ependymoma, ATRT, HGG and medulloblastoma. The Epithelial group, on the other hand, was enriched in craniopharyngioma samples, an epithelium derived tumor. Interestingly, the RNA levels of PD-1 and CTLA4 were significantly upregulated in this Epithelial group, confirming that craniopharingioma could potentially benefit from anti PD-1 and/or CTLA-4 therapies as previously reported. LLG and ganglioglioma were allocated into two groups of Neuronal and Hot, the former characterized by the presence of neuronal cells, and the latter by the presence of macrophages, microglia, and dendritic cells. Adenosine producers (e.g., ENTPD1 and NT5E), which act as inhibitory immune checkpoint molecular, showed up-regulation in the Hot cluster based on both RNAseq and proteome data, suggesting patients in this group might benefit from adenosine reducing treatments. Among LGG tumors, there is a significant difference between microglial and macrophage polarization across BRAF statuses: BRAF-fusion promoted more pro-regenerative (immune suppressive) microglia than pro-inflammatory microglia, while BRAF-V600E promoted more pro-regenerative macrophages than pro-inflammatory macrophages, implying different immunosuppressive mechanisms in the BRAF-V600E and fusion tumors. This study reports the first large-scale deep comprehensive proteogenomic analysis crossing traditional histologic boundaries to uncover foundational pediatric brain tumor biology relating to tumor microenvironment. The incorporation of the proteomic and phosphoproteomic dimension into this large-scale multi-omic study adds functional insight that helps drive translational efforts.

Citation Format: Francesca Petralia, Nicole Tignor, Boris Reva, Pichai Raman, Shrabanti Chowdhury, Dmitry Rykunov, Azra Krek, Weiping Ma, Jiayi Ji, Xiaoyu Song, Yuankun Zhu, Jo Lynne Rokita, Antonio Colaprico, Anna Calinawan, Jeffrey R. Whiteaker, Richard G. Ivey, Zeynep Gumus, Selim Kalayci, Gonzalo L. Garcia, Seungyeul Yoo, Lizabeth Katsnelson, Ying Wang, Jacob J. Kennedy, Uliana J. Voytovich, Lei Zhao, Felipe Leprevost, Hui-Yin Chang, Krutika S. Gaonkar, Elizabeth M. Appert, Ximena Cuellar, Jena Lilly, Jun Zhu, Eric E. Schadt, Medhi Mesri, Emily Boja, Tara Hiltka, Henry Rodriguez, Li Ding, Antonio Iavarone, Maciej Wiznerowicz, Alexey I. Nesvizhskii, David Fenyo, Steven Gygi, Amanda Paulovich, Adam C. Resnick, Phillip B. Storm, Brian Rood, Pei Wang, Children's Brain Tumor Tissue Consortium and Clinical Proteomic Tumor Analysis Consortium. Integrated proteogenomic characterization across seven histological types of pediatric brain tumors [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 445.

Telomere Maintenance Mechanisms Define Clinical Outcome in High-Risk Neuroblastoma

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.

Limited antitumor activity of combined BET and MEK inhibition in neuroblastoma

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.

LIN28B promotes neuroblastoma metastasis and regulates PDZ binding kinase

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.

Epigenomic profiling of neuroblastoma cell lines

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.

CAMKV Is a Candidate Immunotherapeutic Target in MYCN Amplified 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.

ASCL1 is a MYCN- and LMO1-dependent member of the adrenergic neuroblastoma core regulatory circuitry

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.

Combined innate and adaptive immunotherapy overcomes resistance of immunologically cold syngeneic murine neuroblastoma to checkpoint inhibition

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.

Abstract LB-B04: A multi-omic surfaceome study identifies DLK1 as an epigenetically regulated protein and immunotherapeutic target in neuroblastoma

Background: Neuroblastoma (NB) is an embryonal tumor of the sympathetic nervous system that accounts for 12% of childhood cancer deaths. While the introduction of GD2 immunotherapy provides an improvement in time to progression, the therapy is toxic and impact on overall survival is minimal, supporting an urgent need for novel immunotherapies. To date, the cell surface landscape (surfaceome) of NB remains undefined, hindering the identification of immunotherapeutic targets. Methods: To identify NB surfaceome proteins, we performed plasma membrane protein extraction using sucrose gradient ultracentrifugation coupled to mass spectrometry (nLC-MS/MS) in NB cell lines (n=12) and patient derived xenografts (PDX; n=10). These data were integrated with existing RNA-sequencing (NB=153; Normal=7859) and H3K27ac chromatin immunoprecipitation (ChIP)-sequencing data (from overlapping NB cell lines) to evaluate extracellular proteins differentially expressed in NB compared to normal tissues. Candidate targets were validated by immunohistochemistry on NB tumor and normal tissue microarrays (TMAs), flow cytometry and immunofluorescence. In-vitro functional studies were performed following genetic manipulation of candidate targets to assess cell proliferation, differentiation and viability. Finally, we tested ADCT-701 (a DLK1-directed antibody drug conjugate [ADC] with a pyrrolobenzodiazepine [PBD] warhead) in eight PDX models (study ongoing, total of 12 models initiated) with varying levels of DLK1 expression. At enrollment, two mice were each treated with a single dose of saline or 1mg/kg of B12-PL1601 (non-targeting PBD-conjugated ADC) or 1mg/kg ADCT-701 and mice were evaluated for 100 days or until tumor reached 2.0cm3. Results: We yielded on average 66% (range:60-68%) membrane protein enrichment with high reproducibility between biological replicates (80%; range:78-84%) and identified 4826 unique membrane proteins. Our approach confirmed known cell surface proteins in development as immunotherapeutic targets in NB (ALK, GPC2, NCAM1, DLL3 and CD276). Here, we prioritized DLK1 for further evaluation due to it being the only candidate with expression directly associated with a super enhancer element (P=6.09X10-5). RNA-sequencing and tissue microarray analysis of NB and normal tissues showed DLK1 to be overexpressed in a large subset of high-risk NB with minimal expression in normal tissues, excepting adrenal medulla and pituitary. Flow cytometry and immunofluorescence confirmed cell surface expression of DLK1 in a panel of NB cell lines. Genetic depletion of DLK1 using shRNA resulted in neurite outgrowth (P=7.26X10-5) and terminal differentiation. Full proteome analysis of DLK1 knockdown and control cell lines using MS showed regulation of proteins that control outgrowth of neurites (P=3.37X10-3) and development of neurons (P=3.76X10-3). To date, ADCT-701 treatment resulted in maintained complete response (N=2), complete response (N=3) and stable disease (N=1) in models with high DLK1 expression, while those with low/no expression showed disease progression (N=2). Conclusion: DLK1 is an epigenetically regulated immunotherapeutic target in neuroblastoma. ADCT-701 shows potent activity in preclinical models of NB and should be prioritized for clinical development.

Citation Format: Amber K. Weiner, Alexander B. Radaoui, Matthew Tsang, Dan Martinez, Simone Sidoli, Karina L. Conkrite, Alberto Delaidelli, Jo Lynne Rokita, Maria V. Lane, Zalman Vaksman, Komal S. Rathi, Pichai Raman, Jennifer Pogoriler, Tricia Bhatti, Bruce Pawel, Beverly Teicher, Stephen W. Erickson, Poul Sorensen, Yael P. Mosse, Kateryna Krytska, Francesca Zammarchi, Patrick H. van Berkel, Malcolm A. Smith, Benjamin A. Garcia, John M. Maris, Sharon J. Diskin. A multi-omic surfaceome study identifies DLK1 as an epigenetically regulated protein and immunotherapeutic target in neuroblastoma [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference on Molecular Targets and Cancer Therapeutics; 2019 Oct 26-30; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2019;18(12 Suppl):Abstract nr LB-B04. doi:10.1158/1535-7163.TARG-19-LB-B04

PDTM-16. PEDIATRIC HIGH GRADE GLIOMA RESOURCES FROM THE CHILDREN’S BRAIN TUMOR TISSUE CONSORTIUM (CBTTC) AND PEDIATRIC BRAIN TUMOR ATLAS (PBTA)

BACKGROUND

Pediatric high grade glioma (pHGG) remains a fatal disease. Access to richly annotated biospecimens and patient derived tumor models will accelerate pHGG research and support translation of research discoveries. This work describes the pediatric high grade glioma set of the Children’s Brain Tumor Tissue Consortium (CBTTC) from the first release of the Pediatric Brain Tumor Atlas (PBTA).

METHODS

pHGG tumors with associated clinical data and imaging were prospectively collected through the CBTTC and analyzed as the Pediatric Brain Tumor Atlas (PBTA) with processed genomic data deposited into PedcBioPortal for broad access and visualization. Matched tumor was cultured to create high grade glioma cell lines analyzed by targeted and WGS and RNA-seq. A tissue microarray (TMA) of primary pHGG tumors was also created.

RESULTS

The pHGG set includes 87 collection events (73 patients, 60% at diagnosis, median age of 9 yrs, 55% female, 46% hemispheric). Operative reports, pathology reports and histology images are available for nearly all cases. Pre- and post-operative MRI images and reports are also available for a subset. Tumor WGS/RNAseq is available for 70 subjects. Analysis of somatic mutations and copy number alterations of known glioma genes were of expected distribution (36% H3.3, 47% TP53, 24% ATRX and 7% BRAFV600E variants). In our panel of pHGG, six patients (8 tumors) harbored germline mismatch repair mutations with tumor hyper-mutation. A pHGG TMA (n=77), includes 36 patient tumors with matched sequencing. At least one established glioma cell line was generated from 23 patients (32%). Unique reagents include those derived from a H3.3 G34R glioma and from tumors with mismatch repair deficiency.

CONCLUSION

The CBTTC and PBTA have created an openly available integrated resource of over 2,000 tumors, including a rich set of pHGG primary tumors, corresponding cell lines and archival fixed tissue to advance translational research for pHGG.

Genomic Profiling of Childhood Tumor Patient-Derived Xenograft Models to Enable Rational Clinical Trial Design

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.

YAP1 Mediates Resistance to MEK1/2 Inhibition in Neuroblastomas with Hyperactivated RAS Signaling

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₁-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.

The Pediatric Cell Atlas: Defining the Growth Phase of Human Development at Single-Cell Resolution

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.