Background: The molecular mechanisms driving the contrasting clinical courses of neuroblastoma have remained largely enigmatic to date. We therefore aimed to determine the genetic alterations underlying the distinct neuroblastoma phenotypes by massively parallel sequencing.
Methods: Mutation patterns of 416 diagnostic neuroblastoma specimens were determined by massively parallel sequencing (whole-genome sequencing, n=56; whole-exome sequencing, n=162; targeted sequencing, n=198).
Results: We discovered genomic rearrangements of the TERT locus in 12/56 neuroblastomas by whole-genome sequencing. TERT rearrangements affected exclusively high-risk tumors, in mutually exclusive fashion with MYCN amplification and ATRX mutations. In an extended cohort (n=217), TERT rearrangements were detected in 24% of high-risk neuroblastomas, and were associated with particularly poor patient outcome. We demonstrated that such rearrangements led to massive induction of TERT expression and telomerase activation caused by translocation of strong enhancer elements. Similarly, we observed upregulation of TERT and telomerase in MYCN-amplified tumors, while activation of the Alternative Lengthening of Telomeres pathway was detected in high-risk tumors lacking these alterations. In addition, we identified mutations in ALK and other RAS/MAPK pathway-acting genes as well as p53/MDM2 pathway genes in 36/402 (9.0%), 25/402 (6.2%) and 18/402 (4.5%) neuroblastomas, respectively. Such mutations affected predominantly high-risk tumors (51/227 cases, 22.5%) but occurred also in non-high-risk tumors (22/173 cases, 12.7%). Notably, survival of patients whose tumors harbored mutations in these cancer genes was significantly worse than those without such mutations in the entire cohort, and in both the high-risk and the non-high-risk subgroups (p<0.001 each).
Conclusion: Our data demonstrate that activation of telomere maintenance mechanisms is a hallmark of high-risk neuroblastoma, while low-risk tumors lack such mechanisms, which may explain their capacity of spontaneous regression. Furthermore, the clinical course of the disease is strongly influenced by the presence of mutations in RAS/MAPK and p53/MDM2 pathway genes at diagnosis. Together, our data strongly suggest that the clinical phenotype of neuroblastoma is driven by specific patterns of somatic mutations.