Several neuroblastoma (NB) susceptibility loci have been identified within LINC00340, BARD1, LMO1, DUSP12, HSD17B12, DDX4, IL31RA, HACE1 and LIN28B by genome-wide association (GWA) studies including European American individuals. To validate and comprehensively evaluate the impact of the identified NB variants on disease risk and phenotype, we analyzed 16 single nucleotide polymorphisms (SNPs) in an Italian population (370 cases and 809 controls). We assessed their regulatory activity on gene expression in lymphoblastoid (LCLs) and NB cell lines. We evaluated the cumulative effect of the independent loci on NB risk and high-risk phenotype development in Italian and European American (1627 cases and 2575 controls) populations. All NB susceptibility genes replicated in the Italian dataset except for DDX4 and IL31RA, and the most significant SNP was rs6435862 in BARD1 (P = 8.4 × 10(-15)). BARD1 showed an additional and independent SNP association (rs7585356). This variant influenced BARD1 mRNA expression in LCLs and NB cell lines. No evidence of epistasis among the NB-associated variants was detected, whereas a cumulative effect of risk variants on NB risk (European Americans: P (trend) = 6.9 × 10(-30), Italians: P (trend) = 8.55 × 10(13)) and development of high-risk phenotype (European Americans: P (trend) = 6.9 × 10(-13), Italians: P (trend) = 2.2 × 10(-1)) was observed in a dose-dependent manner. These results provide further evidence that the risk loci identified in GWA studies contribute to NB susceptibility in distinct populations and strengthen the role of BARD1 as major genetic contributor to NB risk. This study shows that even in the absence of interaction the combination of several low-penetrance alleles has potential to distinguish subgroups of patients at different risks of developing NB.

The spectrum of somatic mutation of the most aggressive forms of neuroblastoma is not completely determined. We sought to identify potential cancer drivers in clinically aggressive neuroblastoma.Whole exome sequencing was conducted on 17 germline and tumor DNA samples from high-risk patients with adverse events within 36 months from diagnosis (HR-Event3) to identify somatic mutations and deep targeted sequencing of 134 genes selected from the initial screening in additional 48 germline and tumor pairs (62.5% HR-Event3 and high-risk patients), 17 HR-Event3 tumors and 17 human-derived neuroblastoma cell lines.We revealed 22 significantly mutated genes, many of which implicated in cancer progression. Fifteen genes (68.2%) were highly expressed in neuroblastoma supporting their involvement in the disease. CHD9, a cancer driver gene, was the most significantly altered (4.0% of cases) after ALK.Other genes (PTK2, NAV3, NAV1, FZD1 and ATRX), expressed in neuroblastoma and involved in cell invasion and migration were mutated at frequency ranged from 4% to 2%.Focal adhesion and regulation of actin cytoskeleton pathways, were frequently disrupted (14.1% of cases) thus suggesting potential novel therapeutic strategies to prevent disease progression.Notably BARD1, CHEK2 and AXIN2 were enriched in rare, potentially pathogenic, germline variants.In summary, whole exome and deep targeted sequencing identified novel cancer genes of clinically aggressive neuroblastoma. Our analyses show pathway-level implications of infrequently mutated genes in leading neuroblastoma progression.

Purpose: Activating ALK mutations are present in almost 10% of primary neuroblastomas and mark patients for treatment with small-molecule ALK inhibitors in clinical trials. However, recent studies have shown that multiple mechanisms drive resistance to these molecular therapies. We anticipated that detailed mapping of the oncogenic ALK-driven signaling in neuroblastoma can aid to identify potential fragile nodes as additional targets for combination therapies.Experimental Design: To achieve this goal, transcriptome profiling was performed in neuroblastoma cell lines with the ALK F1174L or ALK R1275Q hotspot mutations, ALK amplification, or wild-type ALK following pharmacologic inhibition of ALK using four different compounds. Next, we performed cross-species genomic analyses to identify commonly transcriptionally perturbed genes in MYCN/ALK F1174L double transgenic versus MYCN transgenic mouse tumors as compared with the mutant ALKdriven transcriptome in human neuroblastomas.Results: A 77-gene ALK signature was established and successfully validated in primary neuroblastoma samples, in a neuroblastoma cell line with ALK F1174L and ALK R1275Q regulable overexpression constructs and in other ALKomas. In addition to the previously established PI3K/AKT/mTOR, MAPK/ERK, and MYC/ MYCN signaling branches, we identified that mutant ALK drives a strong upregulation of MAPK negative feedback regulators and upregulates RET and RET-driven sympathetic neuronal markers of the cholinergic lineage.Conclusions: We provide important novel insights into the transcriptional consequences and the complexity of mutant ALK signaling in this aggressive pediatric tumor. The negative feedback loop of MAPK pathway inhibitors may affect novel ALK inhibition therapies, whereas mutant ALK induced RET signaling can offer novel opportunities for testing ALK-RET oriented molecular combination therapies.