Neuroblastoma, a tumour derived from the peripheral sympathetic nervous system, is one of the most frequent solid tumours in childhood. It usually occurs sporadically but familial cases are observed, with a subset of cases occurring in association with congenital malformations of the neural crest being linked to germline mutations of the PHOX2B gene. Here we conducted genome-wide comparative genomic hybridization analysis on a large series of neuroblastomas. Copy number increase at the locus encoding the anaplastic lymphoma kinase (ALK) tyrosine kinase receptor was observed recurrently. One particularly informative case presented a high-level gene amplification that was strictly limited to ALK, indicating that this gene may contribute on its own to neuroblastoma development. Through subsequent direct sequencing of cell lines and primary tumour DNAs we identified somatic mutations of the ALK kinase domain that mainly clustered in two hotspots. Germline mutations were observed in two neuroblastoma families, indicating that ALK is a neuroblastoma predisposition gene. Mutated ALK proteins were overexpressed, hyperphosphorylated and showed constitutive kinase activity. The knockdown of ALK expression in ALK-mutated cells, but also in cell lines overexpressing a wild-type ALK, led to a marked decrease of cell proliferation. Altogether, these data identify ALK as a critical player in neuroblastoma development that may hence represent a very attractive therapeutic target in this disease that is still frequently fatal with current treatments.

Neuroblastoma is a tumor of the peripheral sympathetic nervous system, derived from multipotent neural crest cells (NCCs). To define core regulatory circuitries (CRCs) controlling the gene expression program of neuroblastoma, we established and analyzed the neuroblastoma super-enhancer landscape. We discovered three types of identity in neuroblastoma cell lines: a sympathetic noradrenergic identity, defined by a CRC module including the PHOX2B, HAND2 and GATA3 transcription factors (TFs); an NCC-like identity, driven by a CRC module containing AP-1 TFs; and a mixed type, further deconvoluted at the single-cell level. Treatment of the mixed type with chemotherapeutic agents resulted in enrichment of NCC-like cells. The noradrenergic module was validated by ChIP-seq. Functional studies demonstrated dependency of neuroblastoma with noradrenergic identity on PHOX2B, evocative of lineage addiction. Most neuroblastoma primary tumors express TFs from the noradrenergic and NCC-like modules. Our data demonstrate a previously unknown aspect of tumor heterogeneity relevant for neuroblastoma treatment strategies.

Ewing sarcoma is a primary bone tumor initiated by EWSR1–ETS gene fusions. To identify secondary genetic lesions that contribute to tumor progression, we performed whole-genome sequencing of 112 Ewing sarcoma samples and matched germline DNA. Overall, Ewing sarcoma tumors had relatively few single-nucleotide variants, indels, structural variants and copy-number alterations. Apart from whole chromosome arm copy-number changes, the most common somatic mutations were detected in STAG2 (17%), CDKN2A (12%), TP53 (7%), EZH2, BCOR, and ZMYM3 (2.7% each). Strikingly, STAG2 mutations and CDKN2A deletions were mutually exclusive, as confirmed in Ewing sarcoma cell lines. In an expanded cohort of 299 patients with clinical data, we discovered that STAG2 and TP53 mutations are often concurrent and are associated with poor outcome. Finally, we detected subclonal STAG2 mutations in diagnostic tumors and expansion of STAG2 immuno-negative cells in relapsed tumors as compared with matched diagnostic samples.

The clinical behavior and molecular characteristics of alveolar cases without a fusion gene are indistinguishable from embryonal cases and significantly different from fusion-positive alveolar cases. This implies that fusion gene status irrespective of histology is a critical factor in risk stratification of RMS.

Gene signatures can improve current stratification of patients with RMS but will require complex assays to be developed and extensive validation before clinical application. A significant majority of their prognostic value was encapsulated by the fusion gene status. A continuous risk score derived from the combination of clinical parameters with the presence or absence of PAX3/FOXO1 represents a robust approach to improving current risk-adapted therapy for RMS.