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

Genome-Wide Association Studies (GWAS) have been decisive in elucidating the genetic predisposition of neuroblastoma (NB). The majority of genetic variants identified in GWAS are found in non-coding regions, suggesting that they can be causative of pathogenic dysregulations of gene expression. Nonetheless, pinpointing the potential causal genes within implicated genetic loci remains a major challenge. In this study, we integrated NB GWAS and expression Quantitative Trait Loci (eQTL) data from adrenal gland to identify candidate genes impacting NB susceptibility. We found that ZMYM1, CBL, GSKIP and WDR81 expression was dysregulated by NB predisposing variants. We further investigated the functional role of the identified genes through computational analysis of RNA sequencing (RNA-seq) data from single-cell and whole-tissue samples of NB, neural crest, and adrenal gland tissues, as well as through in vitro differentiation assays in NB cell cultures. Our results indicate that dysregulation of ZMYM1, CBL, GSKIP, WDR81 may lead to malignant transformation by affecting early and late stages of normal program of neuronal differentiation. Our findings enhance the understanding of how specific genes contribute to NB pathogenesis by highlighting their influence on neuronal differentiation and emphasizing the impact of genetic risk variants on the regulation of genes involved in critical biological processes., (© 2024. The Author(s).)

Immunotherapeutic targeting of cell surface proteins is an increasingly effective cancer therapy. However, given the limited number of current targets, the identification of new surface proteins, particularly those with biological importance, is critical. Here, we uncover delta-like non-canonical Notch ligand 1 (DLK1) as a cell surface protein with limited normal tissue expression and high expression in multiple refractory adult metastatic cancers including small cell lung cancer (SCLC) and adrenocortical carcinoma (ACC), a rare cancer with few effective therapies. In ACC, ADCT-701, a DLK1 targeting antibody-drug conjugate (ADC), shows potent in vitro activity among established cell lines and a new cohort of patient-derived organoids as well as robust in vivo anti-tumor responses in cell line-derived and patient-derived xenografts. However, ADCT-701 efficacy is overall limited in ACC due to high expression and activity of the drug efflux protein ABCB1 (MDR1, P-glycoprotein). In contrast, ADCT-701 is extremely potent and induces complete responses in DLK1 + ACC and SCLC in vivo models with low or no ABCB1 expression. Genetic deletion of DLK1 in ACC dramatically downregulates ABCB1 and increases ADC payload and chemotherapy sensitivity through NOTCH1-mediated adrenocortical de-differentiation. Single cell RNA-seq of ACC metastatic tumors reveals significantly decreased adrenocortical differentiation in DLK low or negative cells compared to DLK1 positive cells. This works identifies DLK1 as a novel immunotherapeutic target that regulates tumor cell plasticity and chemoresistance in ACC. Our data support targeting DLK1 with an ADC in ACC and neuroendocrine neoplasms in an active first-in-human phase I clinical trial (NCT06041516)., Competing Interests: Nitin Roper and Jaydira Del Rivero have received research funding from ADC Therapeutics for this study. The other authors have no competing interests to report.

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.

Cancer remains a leading cause of mortality globally. Recent improvements in survival have been facilitated by the development of less toxic immunotherapies; however, identifying targets for immunotherapies remains a challenge in the field. To address this challenge, we developed IMMUNOTAR, a computational tool that systematically prioritizes and identifies candidate immunotherapeutic targets. IMMUNOTAR integrates user-provided RNA-sequencing or proteomics data with quantitative features extracted from publicly available databases based on predefined optimal immunotherapeutic target criteria and quantitatively prioritizes potential surface protein targets. We demonstrate the utility and flexibility of IMMUNOTAR using three distinct datasets, validating its effectiveness in identifying both known and new potential immunotherapeutic targets within the analyzed cancer phenotypes. Overall, IMMUNOTAR enables the compilation of data from multiple sources into a unified platform, allowing users to simultaneously evaluate surface proteins across diverse criteria. By streamlining target identification, IMMUNOTAR empowers researchers to efficiently allocate resources and accelerate immunotherapy development., Competing Interests: Competing Interests: None

Temporal regulation of super-enhancer (SE) driven transcription factors (TFs) underlies normal developmental programs. Neuroblastoma (NB) arises from an inability of sympathoadrenal progenitors to exit a self-renewal program and terminally differentiate. To identify SEs driving TF regulators, we use all-trans retinoic acid (ATRA) to induce NB growth arrest and differentiation. Time-course H3K27ac ChIP-seq and RNA-seq reveal ATRA coordinated SE waves. SEs that decrease with ATRA link to stem cell development (MYCN, GATA3, SOX11). CRISPR-Cas9 and siRNA verify SOX11 dependency, in vitro and in vivo. Silencing the SOX11 SE using dCAS9-KRAB decreases SOX11 mRNA and inhibits cell growth. Other TFs activate in sequential waves at 2, 4 and 8 days of ATRA treatment that regulate neural development (GATA2 and SOX4). Silencing the gained SOX4 SE using dCAS9-KRAB decreases SOX4 expression and attenuates ATRA-induced differentiation genes. Our study identifies oncogenic lineage drivers of NB self-renewal and TFs critical for implementing a differentiation program., (© 2024. The Author(s).)

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: Ewing sarcoma is the second most common bone sarcoma in children, with 1 case per 1.5 million in the United States. Although the survival rate of patients diagnosed with localized disease is approximately 70%, this decreases to approximately 30% for patients with metastatic disease and only approximately 10% for treatment-refractory disease, which have not changed for decades. Therefore, new therapeutic strategies are urgently needed for metastatic and refractory Ewing sarcoma., Experimental Design: This study analyzed 19 unique Ewing sarcoma patient- or cell line-derived xenografts (from 14 primary and 5 metastatic specimens) using proteomics to identify surface proteins for potential immunotherapeutic targeting. Plasma membranes were enriched using density gradient ultracentrifugation and compared with a reference standard of 12 immortalized non-Ewing sarcoma cell lines prepared in a similar manner. In parallel, global proteome analysis was carried out on each model to complement the surfaceome data. All models were analyzed by Tandem Mass Tags-based mass spectrometry to quantify identified proteins., Results: The surfaceome and global proteome analyses identified 1,131 and 1,030 annotated surface proteins, respectively. Among surface proteins identified, both approaches identified known Ewing sarcoma-associated proteins, including IL1RAP, CD99, STEAP1, and ADGRG2, and many new cell surface targets, including ENPP1 and CDH11. Robust staining of ENPP1 was demonstrated in Ewing sarcoma tumors compared with other childhood sarcomas and normal tissues., Conclusions: Our comprehensive proteomic characterization of the Ewing sarcoma surfaceome provides a rich resource of surface-expressed proteins in Ewing sarcoma. This dataset provides the preclinical justification for exploration of targets such as ENPP1 for potential immunotherapeutic application in Ewing sarcoma. See related commentary by Bailey, p. 934., (©2023 The Authors; Published by the American Association for Cancer Research.)

Background: Neuroblastoma is an embryonal cancer of the developing sympathetic nervous system. The genetic contribution of rare pathogenic or likely pathogenic germline variants in patients without a family history remains unclear., Methods: Germline DNA sequencing was performed on 786 neuroblastoma patients. The frequency of rare cancer predisposition gene pathogenic or likely pathogenic variants in patients was compared with 2 cancer-free control cohorts. Matched tumor DNA sequencing was evaluated for second hits, and germline DNA array data from 5585 neuroblastoma patients and 23 505 cancer-free control children were analyzed to identify rare germline copy number variants. Patients with germline pathogenic or likely pathogenic variants were compared with those without to test for association with clinical characteristics, tumor features, and survival., Results: We observed 116 pathogenic or likely pathogenic variants involving 13.9% (109 of 786) of neuroblastoma patients, representing a statistically significant excess burden compared with cancer-free participants (odds ratio [OR] = 1.60, 95% confidence interval [CI] = 1.27 to 2.00). BARD1 harbored the most statistically significant enrichment of pathogenic or likely pathogenic variants (OR = 32.30, 95% CI = 6.44 to 310.35). Rare germline copy number variants disrupting BARD1 were identified in patients but absent in cancer-free participants (OR = 29.47, 95% CI = 1.52 to 570.70). Patients harboring a germline pathogenic or likely pathogenic variant had a worse overall survival compared with those without (P = 8.6 x 10-3)., Conclusions: BARD1 is an important neuroblastoma predisposition gene harboring both common and rare germline pathogenic or likely pathogenic variations. The presence of any germline pathogenic or likely pathogenic variant in a cancer predisposition gene was independently predictive of worse overall survival. As centers move toward paired tumor-normal sequencing at diagnosis, efforts should be made to centralize data and provide an infrastructure to support cooperative longitudinal prospective studies of germline pathogenic variation., (© The Author(s) 2023. Published by Oxford University Press. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

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

Long noncoding RNAs (lncRNA) play an important role in gene regulation and contribute to tumorigenesis. While pan-cancer studies of lncRNA expression have been performed for adult malignancies, the lncRNA landscape across pediatric cancers remains largely uncharted. Here, we curated RNA sequencing data for 1,044 pediatric leukemia and extracranial solid tumors and integrated paired tumor whole genome sequencing and epigenetic data in relevant cell line models to explore lncRNA expression, regulation, and association with cancer. A total of 2,657 lncRNAs were robustly expressed across six pediatric cancers, including 1,142 exhibiting histotype-elevated expression. DNA copy number alterations contributed to lncRNA dysregulation at a proportion comparable to protein coding genes. Application of a multidimensional framework to identify and prioritize lncRNAs impacting gene networks revealed that lncRNAs dysregulated in pediatric cancer are associated with proliferation, metabolism, and DNA damage hallmarks. Analysis of upstream regulation via cell type-specific transcription factors further implicated distinct histotype-elevated and developmental lncRNAs. Integration of these analyses prioritized lncRNAs for experimental validation, and silencing of TBX2-AS1, the top-prioritized neuroblastoma-specific lncRNA, resulted in significant growth inhibition of neuroblastoma cells, confirming the computational predictions. Taken together, these data provide a comprehensive characterization of lncRNA regulation and function in pediatric cancers and pave the way for future mechanistic studies., Significance: Comprehensive characterization of lncRNAs in pediatric cancer leads to the identification of highly expressed lncRNAs across childhood cancers, annotation of lncRNAs showing histotype-specific elevated expression, and prediction of lncRNA gene regulatory networks., (©2023 American Association for Cancer Research.)

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

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like mesenchymal cell state and a more differentiated sympathetic adrenergic cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional cofactor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism, G/T, that affects a GATA motif in the first intron of LMO1. Here, we showed that WT zebrafish with the GATA genotype developed adrenergic neuroblastoma, while knock-in of the protective TATA allele at this locus reduced the penetrance of MYCN-driven tumors, which were restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrated that TATA/TATA tumors also exhibited a mesenchymal cell state and were low risk at diagnosis. Thus, conversion of the regulatory GATA to a TATA allele in the first intron of LMO1 reduced the neuroblastoma-initiation rate by preventing formation of the adrenergic cell state. This mechanism was conserved over 400 million years of evolution, separating zebrafish and humans.

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

Neuroblastoma is a childhood cancer that originates in the developing sympathetic nervous system. We previously reported a crucial role of mitochondrial DNA haplogroups in the pathology of neuroblastoma. To pinpoint mitochondrial DNA variants associated with neuroblastoma risk, we applied a mitochondrial genome imputation pipeline to the single nucleotide polymorphisms array data of 2 pediatric cohorts containing a total of 2404 neuroblastoma patients and 9310 cancer-free controls. All statistical tests were 2-sided. The single nucleotide variant, rs2853493, was statistically significantly associated with neuroblastoma risk in the discovery cohort (odds ratio = 0.62, 95% confidence interval = 0.53 to 0.72, P < .001) and further confirmed in the replication cohort (odds ratio = 0.75, 95% confidence interval = 0.62 to 0.90, P = .002). Further, expression quantitative trait loci analysis indicated genotypes of rs2853493 were associated with expression levels of MT-CYB gene expression in neuroblastoma cells, suggesting rs2853493 may confer risk to neuroblastoma via regulating the expression level of its nearby genes., (© The Author(s) 2022. Published by Oxford University Press. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Background: Neuroblastoma is rarer in African American (AA) children compared with American children of European descent. AA children affected with neuroblastoma, however, more frequently develop the high-risk form of the disease., Methods: We have genotyped an AA cohort of 629 neuroblastoma cases (254 high-risk) and 2,990 controls to investigate genetic susceptibility to neuroblastoma in AAs., Results: We confirmed the known neuroblastoma susceptibility gene BARD1 at genome-wide significance in the subset of high-risk cases. We also estimated local admixture across the autosomal genome in the AA cases and controls and detected a signal at 4q31.22 where cases show an increase in European ancestry. A region at 17p13.1 showed increased African ancestry in the subgroup of high-risk cases with respect to intermediate- and low-risk cases. Using results from our published European American (EA) genome-wide association study (GWAS), we found that a polygenic score that included all independent SNPs showed a highly significant association (P value = 1.8 × 10-73) and explained 19% of disease risk variance in an independent EA cohort. In contrast, the best fit polygenic score (P value = 3.2 × 10-11) in AAs included only 22 independent SNPs with association P value < 2.75 × 10-6 in the EA GWAS, and explained 2% of neuroblastoma risk variance. The significance of the polygenic score dropped rapidly with inclusion of additional SNPs., Conclusions: These findings suggest that several common variants contribute to risk of neuroblastoma in an ancestry-specific fashion., Impact: This work supports the need for GWAS to be performed in populations of all races and ethnicities., (©2022 American Association for Cancer Research.)

Motivation: Despite widespread prevalence of somatic structural variations (SVs) across most tumor types, understanding of their molecular implications often remains poor. SVs are extremely heterogeneous in size and complexity, hindering the interpretation of their pathogenic role. Tools integrating large SV datasets across platforms are required to fully characterize the cancer's somatic landscape., Results: svpluscnv R package is a swiss army knife for the integration and interpretation of orthogonal datasets including copy number variant segmentation profiles and sequencing-based structural variant calls. The package implements analysis and visualization tools to evaluate chromosomal instability and ploidy, identify genes harboring recurrent SVs and detects complex rearrangements such as chromothripsis and chromoplexia. Further, it allows systematic identification of hot-spot shattered genomic regions, showing reproducibility across alternative detection methods and datasets., Availability and Implementation: https://github.com/ccbiolab/svpluscnv., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Ornithine decarboxylase (ODC1), a critical regulatory enzyme in polyamine biosynthesis, is a direct transcriptional target of MYCN, amplification of which is a powerful marker of aggressive neuroblastoma. A single nucleotide polymorphism (SNP), G316A, within the first intron of ODC1 , results in genotypes wildtype GG, and variants AG/AA. CRISPR-cas9 technology was used to investigate the effects of AG clones from wildtype MYCN -amplified SK-N-BE(2)-C cells and the effect of the SNP on MYCN binding, and promoter activity was investigated using EMSA and luciferase assays. AG clones exhibited decreased ODC1 expression, growth rates, and histone acetylation and increased sensitivity to ODC1 inhibition. MYCN was a stronger transcriptional regulator of the ODC1 promoter containing the G allele, and preferentially bound the G allele over the A. Two neuroblastoma cohorts were used to investigate the clinical impact of the SNP. In the study cohort, the minor AA genotype was associated with improved survival, while poor prognosis was associated with the GG genotype and AG/GG genotypes in MYCN -amplified and non-amplified patients, respectively. These effects were lost in the GWAS cohort. We have demonstrated that the ODC1 G316A polymorphism has functional significance in neuroblastoma and is subject to allele-specific regulation by the MYCN oncoprotein.

Background: Neuroblastoma is a childhood malignancy that arises from the developing sympathetic nervous system. Although mitochondrial dysfunctions have been implicated in the pathophysiology of neuroblastoma, the role of mitochondrial DNA (mtDNA) has not been extensively investigated., Methods: A total of 2404 Caucasian children diagnosed with neuroblastoma and 9310 ancestry-matched controls were recruited at the Children's Hospital of Philadelphia. The mtDNA haplogroups were identified from SNP array data of two independent cohorts. We conducted a case-control study to explore potential associations of mtDNA haplogroups with the susceptibility of neuroblastoma. The genetic effect of neuroblastoma was measured by odds ratios (ORs) of mitochondrial haplogroups. All tests were two-sided., Results: Haplogroup K was statistically significantly associated with reduced risk of neuroblastoma in the discovery cohort consisting of 1474 cases and 5699 controls (OR = 0.72, 95% confidence interval [CI] = 0.57 to 0.90; P = 4.8 × 10-3). The association was replicated in an independent cohort (OR = 0.69, 95% CI = 0.53 to 0.92; P = .01) of 930 cases and 3611 controls. Pooled analysis was performed by combining the two data sets. The association remained highly statistically significant after correction for multiple testing (OR = 0.71, 95% CI = 0.59 to 0.84, P = 1.96 × 10-4, Pcorrected = .002). Further analysis focusing on neuroblastoma subtypes indicated haplogroup K was more associated with high-risk neuroblastoma (OR = 0.57, 95% CI = 0.43 to 0.76; P = 1.46 × 10-4) than low-risk and intermediate-risk neuroblastoma., Conclusions: Haplogroup K is an independent genetic factor associated with reduced risk of developing neuroblastoma in European descents. These findings provide new insights into the genetic basis of neuroblastoma, implicating mitochondrial DNA encoded proteins in the etiology of neuroblastoma., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Background: Ependymoma is a histologically defined central nervous system tumor most commonly occurring in childhood. Population-level incidence differences by race/ethnicity are observed, with individuals of European ancestry at highest risk. We aimed to determine whether extent of European genetic ancestry is associated with ependymoma risk in US populations., Methods: In a multi-ethnic study of Californian children (327 cases, 1970 controls), we estimated the proportions of European, African, and Native American ancestry among recently admixed Hispanic and African American subjects and estimated European admixture among non-Hispanic white subjects using genome-wide data. We tested whether genome-wide ancestry differences were associated with ependymoma risk and performed admixture mapping to identify associations with local ancestry. We also evaluated race/ethnicity-stratified ependymoma incidence data from the Central Brain Tumor Registry of the United States (CBTRUS)., Results: CBTRUS data revealed that African American and Native American children have 33% and 36%, respectively, reduced incidence of ependymoma compared with non-Hispanic whites. In genetic analyses, a 20% increase in European ancestry was associated with a 1.31-fold higher odds of ependymoma among self-reported Hispanics and African Americans (95% CI: 1.08-1.59, Pmeta = 6.7 × 10-3). Additionally, eastern European ancestral substructure was associated with increased ependymoma risk in non-Hispanic whites (P = 0.030) and in Hispanics (P = 0.043). Admixture mapping revealed a peak at 20p13 associated with increased local European ancestry, and targeted fine-mapping identified a lead variant at rs6039499 near RSPO4 (odds ratio = 1.99; 95% CI: 1.45-2.73; P = 2.2 × 10-5) but which was not validated in an independent set of posterior fossa type A patients., Conclusions: Interethnic differences in ependymoma risk are recapitulated in the genomic ancestry of ependymoma patients, implicating regions to target in future association studies., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Neuro-Oncology. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Ependymoma is the third most common brain tumor in children, with well-described molecular characterization but poorly understood underlying germline risk factors. To investigate whether genetic predisposition to longer telomere length influences ependymoma risk, we utilized case-control data from three studies: a population-based pediatric and adolescent ependymoma case-control sample from California (153 cases, 696 controls), a hospital-based pediatric posterior fossa type A (EPN-PF-A) ependymoma case-control study from Toronto's Hospital for Sick Children and the Children's Hospital of Philadelphia (83 cases, 332 controls), and a multicenter adult-onset ependymoma case-control dataset nested within the Glioma International Case-Control Consortium (GICC) (103 cases, 3287 controls). In the California case-control sample, a polygenic score for longer telomere length was significantly associated with increased risk of ependymoma diagnosed at ages 12-19 (P = 4.0 × 10 -3 ), but not with ependymoma in children under 12 years of age (P = 0.94). Mendelian randomization supported this observation, identifying a significant association between genetic predisposition to longer telomere length and increased risk of adolescent-onset ependymoma (OR PRS  = 1.67; 95% CI 1.18-2.37; P = 3.97 × 10 -3 ) and adult-onset ependymoma (P MR-Egger  = 0.042), but not with risk of ependymoma diagnosed before age 12 (OR = 1.12; 95% CI 0.94-1.34; P = 0.21), nor with EPN-PF-A (P MR-Egger  = 0.59). These findings complement emerging literature suggesting that augmented telomere maintenance is important in ependymoma pathogenesis and progression, and that longer telomere length is a risk factor for diverse nervous system malignancies.

Neuroblastoma is a malignancy of the developing sympathetic nervous system that accounts for 12% of childhood cancer deaths. Like many childhood cancers, neuroblastoma shows a relative paucity of somatic single-nucleotide variants (SNVs) and small insertions and deletions (indels) compared to adult cancers. Here, we assessed the contribution of somatic structural variation (SV) in neuroblastoma using a combination of whole-genome sequencing (WGS) of tumor-normal pairs ( n = 135) and single-nucleotide polymorphism (SNP) genotyping of primary tumors ( n = 914). Our study design allowed for orthogonal validation and replication across platforms. SV frequency, type, and localization varied significantly among high-risk tumors. MYCN nonamplified high-risk tumors harbored an increased SV burden overall, including a significant excess of tandem duplication events across the genome. Genes disrupted by SV breakpoints were enriched in neuronal lineages and associated with phenotypes such as autism spectrum disorder (ASD). The postsynaptic adapter protein-coding gene, SHANK2 , located on Chromosome 11q13, was disrupted by SVs in 14% of MYCN nonamplified high-risk tumors based on WGS and 10% in the SNP array cohort. Expression of SHANK2 was low across human-derived neuroblastoma cell lines and high-risk neuroblastoma tumors. Forced expression of SHANK2 in neuroblastoma cells resulted in significant growth inhibition ( P = 2.6 × 10 -2 to 3.4 × 10 -5 ) and accelerated neuronal differentiation following treatment with all- trans retinoic acid ( P = 3.1 × 10 -13 to 2.4 × 10 -30 ). These data further define the complex landscape of somatic structural variation in neuroblastoma and suggest that events leading to deregulation of neurodevelopmental processes, such as inactivation of SHANK2 , are key mediators of tumorigenesis in this childhood cancer., (© 2020 Lopez et al.; Published by Cold Spring Harbor Laboratory Press.)

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

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.

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

Background: Neuroblastoma is a biologically and clinically heterogeneous disease. Based on recent studies demonstrating an association between the primary tumor site, prognosis, and commonly measured tumor biological features, we hypothesized that neuroblastomas arising in different sites would show distinct genomic features reflective of the developmental biology of the sympathicoadrenal nervous system., Methods: We first compared genomic and epigenomic data of primary diagnostic neuroblastomas originating in the adrenal gland (n = 646) compared to thoracic sympathetic ganglia (n = 118). We also evaluated association of common germline variation with these primary sites in 1027 European-American neuroblastoma patients., Results: We observed higher rates of MYCN amplification, chromosome 1q gain, and chromosome 11q deletion among adrenal tumors, which were highly predictive of functional RNA signatures. Surprisingly, thoracic neuroblastomas were more likely to harbor ALK driver mutations than adrenal cases among all cases (odds ratio = 1.89, 95% confidence interval = 1.04 to 3.43), and among cases without MYCN amplification (odds ratio = 2.86, 95% confidence interval = 1.48 to 5.49). Common germline single nucleotide polymorphisms (SNPs) in BARD1 (previously associated with high-risk neuroblastoma) were found to be strongly associated with predisposition for origin at adrenal, rather than thoracic, sites., Conclusions: Neuroblastomas arising in the adrenal gland are more likely to harbor structural DNA aberrations including MYCN amplification, whereas thoracic tumors show defects in mitotic checkpoints resulting in hyperdiploidy. Despite the general association of ALK mutations with high-risk disease, thoracic tumors are more likely to harbor gain-of-function ALK aberrations. Site of origin is likely reflective of stage of sympathetic nervous system development when malignant transformation occurs and is a surrogate for underlying tumor biology., (© The Author(s) 2019. Published by Oxford University Press. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

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

Background: There is evidence vitamin A plays a role in neuroblastoma. Not only is 13-cis-retinoic acid used as maintenance therapy for high-risk cases, but prenatal vitamin intake use may decrease neuroblastoma risk. We hypothesized that single nucleotide polymorphisms (SNPs) in vitamin A-related genes are may be associated with neuroblastoma risk and potentially be modified by vitamin A intake., Methods: The Neuroblastoma Epidemiology in North America (NENA) study recruited 563 case-parent sets through the Children's Oncology Group's Childhood Cancer Research Network. We ascertained dietary nutrient intake through questionnaires and genotyped 463 SNPs in vitamin A-related genes from saliva DNA. Offspring and maternal log-additive risk ratios (RR) and stratum-specific RR for gene-environment interaction were estimated with a log-linear model. We avoided false positives due to multiple testing by using the false discovery rate (FDR)., Results: When all neuroblastoma cases were considered together, no offspring variants met the significance criteria (FDR Q-value < 0.2). One maternal SNP (rs12442054) was associated with decreased risk of neuroblastoma (RR: 0.61; 95% Confidence Interval (CI): 0.47-0.79, Q = 0.076). When the cases were categorized according to prognostic risk category and age at onset, nine offspring SNPs were significantly associated with intermediate-risk neuroblastoma. Maternal rs6776706 was associated with (RR: 0.49; 95% CI: 0.33-0.72, Q = 0.161) high-risk neuroblastoma and maternal rs11103603 (RR: 0.60; 95% CI: 0.45-0.79, Q = 0.127) was associated with neuroblastoma aged <1 year. For gene-environment interaction, maternal rs729147 was associated with decreased risk of neuroblastoma among mothers with vitamin A consumption above the recommendation., Conclusions: Although there is biologic plausibility for the role of vitamin A in neuroblastoma, we found weak evidence of a relationship between vitamin A related genes and neuroblastoma., (Copyright © 2019 Elsevier Ltd. 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.)

Here we present a free interactive web tool to process and visualize proteomics data sets with a single click. GiaPronto can process all proteomics quantification methods, i.e. label-free, SILAC and isobaric labeling, and analyze post-translational modifications (PTMs). The software performs normalization and statistics, assists determination of regulated proteins, biomarkers and Gene Ontology (GO) enrichment, and provides high resolution images and tables for further data analysis. We foresee that GiaPronto will become the most rapid and simple tool for assessing data quality and determining most relevant features of proteomic data sets. GiaPronto is available at giapronto.diskinlab.org., (© 2018 by The American Society for Biochemistry and Molecular Biology, Inc.)

Neuroblastoma, a malignancy of the developing peripheral nervous system that affects infants and young children, is a complex genetic disease. Over the past two decades, significant progress has been made toward understanding the genetic determinants that predispose to this often lethal childhood cancer. Approximately 1-2% of neuroblastomas are inherited in an autosomal dominant fashion and a combination of co-morbidity and linkage studies has led to the identification of germline mutations in PHOX2B and ALK as the major genetic contributors to this familial neuroblastoma subset. The genetic basis of "sporadic" neuroblastoma is being studied through a large genome-wide association study (GWAS). These efforts have led to the discovery of many common susceptibility alleles, each with modest effect size, associated with the development and progression of sporadic neuroblastoma. More recently, next-generation sequencing efforts have expanded the list of potential neuroblastoma-predisposing mutations to include rare germline variants with a predicted larger effect size. The evolving characterization of neuroblastoma's genetic basis has led to a deeper understanding of the molecular events driving tumorigenesis, more precise risk stratification and prognostics and novel therapeutic strategies. This review details the contemporary understanding of neuroblastoma's genetic predisposition, including recent advances and discusses ongoing efforts to address gaps in our knowledge regarding this malignancy's complex genetic underpinnings.

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.