Search

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

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

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

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.

The genetic aetiology of sporadic neuroblastoma is still largely unknown. We have identified diverse neuroblastoma susceptibility loci by genomewide association studies (GWASs); however, additional SNPs that likely contribute to neuroblastoma susceptibility prompted this investigation for identification of additional variants that are likely hidden among signals discarded by the multiple testing corrections used in the analysis of genomewide data. There is evidence suggesting the CDKN1B, coding for the cycle inhibitor p27Kip1, is involved in neuroblastoma. We thus assess whether genetic variants of CDKN1B are associated with neuroblastoma. We imputed all possible genotypes across CDKN1B locus on a discovery case series of 2101 neuroblastoma patients and 4202 genetically matched controls of European ancestry. The most significantly associated rs34330 was analysed in an independent Italian cohort of 311 cases and 709 controls. In vitro functional analysis was carried out in HEK293T and in neuroblastoma cell line SHEP-2, both transfected with pGL3-CDKN1B-CC or pGL3-CDKN1B-TT constructs. We identified an association of the rs34330 T allele (-79C/T) with the neuroblastoma risk (P combined = 0.002; OR = 1.17). The risk allele (T) of this single nucleotide polymorphism led to a lower transcription rate in cells transfected with a luciferase reporter driven by the polymorphic p27Kip1 promoter (P < 0.05). Three independent sets of neuroblastoma tumours carrying -79TT genotype showed a tendency towards lower CDKN1B mRNA levels. Our study shows that a functional variant, associated with a reduced CDKN1B gene transcription, influences neuroblastoma susceptibility., (© 2017 The Authors. Journal of Cellular and Molecular Medicine published by John Wiley & Sons Ltd and Foundation for Cellular and Molecular Medicine.)

To investigate genetic predispositions for MYCN-amplified neuroblastoma, we performed a meta-analysis of three genome-wide association studies totaling 615 MYCN-amplified high-risk neuroblastoma cases and 1869 MYCN-nonamplified non-high-risk neuroblastoma cases as controls using a fixed-effects model with inverse variance weighting. All statistical tests were two-sided. We identified a novel locus at 3p21.31 indexed by the single nucleotide polymorphism (SNP) rs80059929 (odds ratio [OR] = 2.95, 95% confidence interval [CI] = 2.17 to 4.02, Pmeta = 6.47 × 10-12) associated with MYCN-amplified neuroblastoma, which was replicated in 127 MYCN-amplified cases and 254 non-high-risk controls (OR = 2.30, 95% CI = 1.12 to 4.69, Preplication = .02). To confirm this signal is exclusive to MYCN-amplified tumors, we performed a second meta-analysis comparing 728 MYCN-nonamplified high-risk patients to identical controls. rs80059929 was not statistically significant in MYCN-nonamplified high-risk patients (OR = 1.24, 95% CI = 0.90 to 1.71, Pmeta = .19). SNP rs80059929 is within intron 16 in the KIF15 gene. Additionally, the previously reported LMO1 neuroblastoma risk locus was statistically significant only in patients with MYCN-nonamplified high-risk tumors (OR = 0.63, 95% CI = 0.53 to 0.75, Pmeta = 1.51 × 10-8; Pmeta = .95). Our results indicate that common genetic variation predisposes to different neuroblastoma genotypes, including the likelihood of somatic MYCN-amplification., (© The Author 2017. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

MYCN amplification and 11q deletion are two inversely correlated prognostic factors of poor outcome in neuroblastoma. Here we identify common variants at 11q22.2 within MMP20 that associate with neuroblastoma cases harboring 11q deletion (rs10895322), using GWAS in 113 European-American cases and 5109 ancestry-matched controls. The association is replicated in 44 independent cases and 1902 controls. Our study yields novel insights into the genetic underpinnings of neuroblastoma, demonstrating that the inherited common variants reported contribute to the origin of intra-tumor genetic heterogeneity in neuroblastoma.Chromosomal abnormalities such as 11q deletion are associated with poor prognosis in neuroblastoma. Here, the authors perform a genome-wide association study and identify an association between a variant within a Matrix metalloproteinase (MMP) gene member, MMP20, and 11q-deletion subtype neuroblastoma.

Neuroblastoma is a cancer of the developing sympathetic nervous system that most commonly presents in young children and accounts for approximately 12% of pediatric oncology deaths. Here, we report on a genome-wide association study (GWAS) in a discovery cohort or 2,101 cases and 4,202 controls of European ancestry. We identify two new association signals at 3q25 and 4p16 that replicated robustly in multiple independent cohorts comprising 1,163 cases and 4,396 controls (3q25: rs6441201 combined P = 1.2x10-11, Odds Ratio 1.23, 95% CI:1.16-1.31; 4p16: rs3796727 combined P = 1.26x10-12, Odds Ratio 1.30, 95% CI: 1.21-1.40). The 4p16 signal maps within the carboxypeptidase Z (CPZ) gene. The 3q25 signal resides within the arginine/serine-rich coiled-coil 1 (RSRC1) gene and upstream of the myeloid leukemia factor 1 (MLF1) gene. Increased expression of MLF1 was observed in neuroblastoma cells homozygous for the rs6441201 risk allele (P = 0.02), and significant growth inhibition was observed upon depletion of MLF1 (P < 0.0001) in neuroblastoma cells. Taken together, we show that common DNA variants within CPZ at 4p16 and upstream of MLF1 at 3q25 influence neuroblastoma susceptibility and MLF1 likely plays an important role in neuroblastoma tumorigenesis.

Neuroblastoma is a paediatric malignancy that typically arises in early childhood, and is derived from the developing sympathetic nervous system. Clinical phenotypes range from localized tumours with excellent outcomes to widely metastatic disease in which long-term survival is approximately 40% despite intensive therapy. A previous genome-wide association study identified common polymorphisms at the LMO1 gene locus that are highly associated with neuroblastoma susceptibility and oncogenic addiction to LMO1 in the tumour cells. Here we investigate the causal DNA variant at this locus and the mechanism by which it leads to neuroblastoma tumorigenesis. We first imputed all possible genotypes across the LMO1 locus and then mapped highly associated single nucleotide polymorphism (SNPs) to areas of chromatin accessibility, evolutionary conservation and transcription factor binding sites. We show that SNP rs2168101 G>T is the most highly associated variant (combined P = 7.47 × 10(-29), odds ratio 0.65, 95% confidence interval 0.60-0.70), and resides in a super-enhancer defined by extensive acetylation of histone H3 lysine 27 within the first intron of LMO1. The ancestral G allele that is associated with tumour formation resides in a conserved GATA transcription factor binding motif. We show that the newly evolved protective TATA allele is associated with decreased total LMO1 expression (P = 0.028) in neuroblastoma primary tumours, and ablates GATA3 binding (P < 0.0001). We demonstrate allelic imbalance favouring the G-containing strand in tumours heterozygous for this SNP, as demonstrated both by RNA sequencing (P < 0.0001) and reporter assays (P = 0.002). These findings indicate that a recently evolved polymorphism within a super-enhancer element in the first intron of LMO1 influences neuroblastoma susceptibility through differential GATA transcription factor binding and direct modulation of LMO1 expression in cis, and this leads to an oncogenic dependency in tumour cells.

The MYC oncogene encodes MYC, a transcription factor that binds the genome through sites termed E-boxes (5'-CACGTG-3'), which are identical to the binding sites of the heterodimeric CLOCK-BMAL1 master circadian transcription factor. Hence, we hypothesized that ectopic MYC expression perturbs the clock by deregulating E-box-driven components of the circadian network in cancer cells. We report here that deregulated expression of MYC or N-MYC disrupts the molecular clock in vitro by directly inducing REV-ERBα to dampen expression and oscillation of BMAL1, and this could be rescued by knockdown of REV-ERB. REV-ERBα expression predicts poor clinical outcome for N-MYC-driven human neuroblastomas that have diminished BMAL1 expression, and re-expression of ectopic BMAL1 in neuroblastoma cell lines suppresses their clonogenicity. Further, ectopic MYC profoundly alters oscillation of glucose metabolism and perturbs glutaminolysis. Our results demonstrate an unsuspected link between oncogenic transformation and circadian and metabolic dysrhythmia, which we surmise to be advantageous for cancer., (Copyright © 2015 Elsevier Inc. All rights reserved.)

A more complete understanding of aberrant oncogenic signaling in neuroblastoma, a malignancy of the developing sympathetic nervous system, is paramount to improving patient outcomes. Recently, we identified LIN28B as an oncogenic driver in high-risk neuroblastoma. Here, we identify the oncogene RAN as a LIN28B target and show regional gain of chromosome 12q24 as an additional somatic alteration resulting in increased RAN expression. We show that LIN28B influences RAN expression by promoting RAN Binding Protein 2 expression and by directly binding RAN mRNA. Further, we demonstrate a convergence of LIN28B and RAN signaling on Aurora kinase A activity. Collectively, these findings demonstrate that LIN28B-RAN-AURKA signaling drives neuroblastoma oncogenesis, suggesting that this pathway may be amenable to therapeutic targeting., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Chromosome 6p22 was identified recently as a neuroblastoma susceptibility locus, but its mechanistic contributions to tumorigenesis are as yet undefined. Here we report that the most highly significant single-nucleotide polymorphism (SNP) associations reside within CASC15, a long noncoding RNA that we define as a tumor suppressor at 6p22. Low-level expression of a short CASC15 isoform (CASC15-S) associated highly with advanced neuroblastoma and poor patient survival. In human neuroblastoma cells, attenuating CASC15-S increased cellular growth and migratory capacity. Gene expression analysis revealed downregulation of neuroblastoma-specific markers in cells with attenuated CASC15-S, with concomitant increases in cell adhesion and extracellular matrix transcripts. Altogether, our results point to CASC15-S as a mediator of neural growth and differentiation, which impacts neuroblastoma initiation and progression., (©2015 American Association for Cancer Research.)