Neuroblastoma arises in early fetal development and its evolutionary duration predicts outcome.

Neuroblastoma, the most frequent solid tumor in infants, shows very diverse outcomes from spontaneous regression to fatal disease. When these different tumors originate and how they evolve are not known. Here we quantify the somatic evolution of neuroblastoma by deep whole-genome sequencing, molecular clock analysis and population-genetic modeling in a comprehensive cohort covering all subtypes. We find that tumors across the entire clinical spectrum begin to develop via aberrant mitoses as early as the first trimester of pregnancy. Neuroblastomas with favorable prognosis expand clonally after short evolution, whereas aggressive neuroblastomas show prolonged evolution during which they acquire telomere maintenance mechanisms. The initial aneuploidization events condition subsequent evolution, with aggressive neuroblastoma exhibiting early genomic instability. We find in the discovery cohort (n = 100), and validate in an independent cohort (n = 86), that the duration of evolution is an accurate predictor of outcome. Thus, insight into neuroblastoma evolution may prospectively guide treatment decisions.

Inherited rare variants in homologous recombination and neurodevelopmental genes are associated with increased risk of neuroblastoma.

BACKGROUND: Neuroblastoma (NB) is the most common solid extracranial paediatric tumour. Genome-wide association studies have driven the discovery of common risk variants, but no large study has investigated the contribution of rare variants to NB susceptibility. Here, we conducted a whole-exome sequencing (WES) of 664 NB cases and 822 controls and used independent validation datasets to identify genes with rare risk variants and involved pathways. METHODS: WES was performed at 50× depth and variants were jointly called in cases and controls. We developed two models to identify mutations with high clinical impact (P/LP model) and to discover less penetrant risk mutations affecting non-canonical cancer pathways (RPV model). We performed a gene-level collapsing test using Firth's logistic regression in 242 selected cancer predisposition genes (CPGs) and a gene-sets burden analysis of biologically-informed pathways. FINDINGS: Twelve percent of patients carried P/LP variants in CPGs and showed a significant enrichment (P = 2.3 × 10-4) compared to controls (6%). We identified P/LP variants in 45 CPGs enriched in homologous recombination (HR) pathway. The most P/LP enriched genes in NB were BRCA1, ALK and RAD51C. Additionally, we found higher RPV burden in gene-sets of neuron differentiation, neural tube development and synapse assembly, and in gene-sets associated with neurodevelopmental disorders (NDD). INTERPRETATION: The high fraction of NB patients with P/LP variants indicates the need of genetic counselling. Furthermore, inherited rare variants predispose to NB development by affecting mechanisms related to HR and neurodevelopmental processes, and demonstrate that NDD genes are altered in NB at the germline level. FUNDING: Associazione Italiana per la Ricerca sul Cancro, Fondazione Italiana per la Lotta al Neuroblastoma, Associazione Oncologia Pediatrica e Neuroblastoma, Regione Campania, Associazione Giulio Adelfio onlus, and Italian Health Ministry.

MYCN mediates cysteine addiction and sensitizes neuroblastoma to ferroptosis.

Aberrant expression of MYC transcription factor family members predicts poor clinical outcome in many human cancers. Oncogenic MYC profoundly alters metabolism and mediates an antioxidant response to maintain redox balance. Here we show that MYCN induces massive lipid peroxidation on depletion of cysteine, the rate-limiting amino acid for glutathione (GSH) biosynthesis, and sensitizes cells to ferroptosis, an oxidative, non-apoptotic and iron-dependent type of cell death. The high cysteine demand of MYCN-amplified childhood neuroblastoma is met by uptake and transsulfuration. When uptake is limited, cysteine usage for protein synthesis is maintained at the expense of GSH triggering ferroptosis and potentially contributing to spontaneous tumor regression in low-risk neuroblastomas. Pharmacological inhibition of both cystine uptake and transsulfuration combined with GPX4 inactivation resulted in tumor remission in an orthotopic MYCN-amplified neuroblastoma model. These findings provide a proof of concept of combining multiple ferroptosis targets as a promising therapeutic strategy for aggressive MYCN-amplified tumors.

Functional Therapeutic Target Validation Using Pediatric Zebrafish Xenograft Models.

The survival rate among children with relapsed tumors remains poor, due to tumor heterogeneity, lack of directly actionable tumor drivers and multidrug resistance. Novel personalized medicine approaches tailored to each tumor are urgently needed to improve cancer treatment. Current pediatric precision oncology platforms, such as the INFORM (INdividualized Therapy FOr Relapsed Malignancies in Childhood) study, reveal that molecular profiling of tumor tissue identifies targets associated with clinical benefit in a subgroup of patients only and should be complemented with functional drug testing. In such an approach, patient-derived tumor cells are exposed to a library of approved oncological drugs in a physiological setting, e.g., in the form of animal avatars injected with patient tumor cells. We used molecularly fully characterized tumor samples from the INFORM study to compare drug screen results of individual patient-derived cell models in functional assays: (i) patient-derived spheroid cultures within a few days after tumor dissociation; (ii) tumor cells reisolated from the corresponding mouse PDX; (iii) corresponding long-term organoid-like cultures and (iv) drug evaluation with the corresponding zebrafish PDX (zPDX) model. Each model had its advantage and complemented the others for drug hit and drug combination selection. Our results provide evidence that in vivo zPDX drug screening is a promising add-on to current functional drug screening in precision medicine platforms.

ALT-FISH quantifies alternative lengthening of telomeres activity by imaging of single-stranded repeats.

Alternative lengthening of telomeres (ALT) occurs in ∼10% of cancer entities. However, little is known about the heterogeneity of ALT activity since robust ALT detection assays with high-throughput in situ readouts are lacking. Here, we introduce ALT-FISH, a method to quantitate ALT activity in single cells from the accumulation of single-stranded telomeric DNA and RNA. It involves a one-step fluorescent in situ hybridization approach followed by fluorescence microscopy imaging. Our method reliably identified ALT in cancer cell lines from different tumor entities and was validated in three established models of ALT induction and suppression. Furthermore, we successfully applied ALT-FISH to spatially resolve ALT activity in primary tissue sections from leiomyosarcoma and neuroblastoma tumors. Thus, our assay provides insights into the heterogeneity of ALT tumors and is suited for high-throughput applications, which will facilitate screening for ALT-specific drugs.

Neuroblastoma-Telomere maintenance, deregulated signaling transduction and beyond.

The childhood malignancy neuroblastoma belongs to the group of embryonal tumors and originates from progenitor cells of the sympathoadrenal lineage. Treatment options for children with high-risk and relapsed disease are still very limited. In recent years, an ever-growing molecular diversity was identified using (epi)-genetic profiling of neuroblastoma tumors, indicating that molecularly targeted therapies could be a promising therapeutic option. In this review article, we summarize the various molecular subtypes and genetic events associated with neuroblastoma and describe recent advances in targeted therapies. We lay a strong emphasis on the importance of telomere maintenance mechanisms for understanding tumor progression and risk classification of neuroblastoma.

iTReX: Interactive exploration of mono- and combination therapy dose response profiling data.

High throughput screening methods, measuring the sensitivity and resistance of tumor cells to drug treatments have been rapidly evolving. Not only do these screens allow correlating response profiles to tumor genomic features for developing novel predictors of treatment response, but they can also add evidence for therapy decision making in precision oncology. Recent analysis methods developed for either assessing single agents or combination drug efficacies enable quantification of dose-response curves with restricted symmetric fit settings. Here, we introduce iTReX, a user-friendly and interactive Shiny/R application, for both the analysis of mono- and combination therapy responses. The application features an extended version of the drug sensitivity score (DSS) based on the integral of an advanced five-parameter dose-response curve model and a differential DSS for combination therapy profiling. Additionally, iTReX includes modules that visualize drug target interaction networks and support the detection of matches between top therapy hits and the sample omics features to enable the identification of druggable targets and biomarkers. iTReX enables the analysis of various quantitative drug or therapy response readouts (e.g. luminescence, fluorescence microscopy) and multiple treatment strategies (drug treatments, radiation). Using iTReX we validate a cost-effective drug combination screening approach and reveal the application's ability to identify potential sample-specific biomarkers based on drug target interaction networks. The iTReX web application is accessible at https://itrex.kitz-heidelberg.de.

Combining APR-246 and HDAC-Inhibitors: A Novel Targeted Treatment Option for Neuroblastoma.

APR-246 (Eprenetapopt/PRIMA-1Met) is a very potent anti-cancer drug in clinical trials and was initially developed as a p53 refolding agent. As an alternative mode of action, the elevation of reactive oxygen species (ROS) has been proposed. Through an in silico analysis, we investigated the responses of approximately 800 cancer cell lines (50 entities; Cancer Therapeutics Response Portal, CTRP) to APR-246 treatment. In particular, neuroblastoma, lymphoma and acute lymphocytic leukemia cells were highly responsive. With gene expression data from the Cancer Cell Line Encyclopedia (CCLE; n = 883) and patient samples (n = 1643) from the INFORM registry study, we confirmed that these entities express low levels of SLC7A11, a previously described predictive biomarker for APR-246 responsiveness. Combining the CTRP drug response data with the respective CCLE gene expression profiles, we defined a novel gene signature, predicting the effectiveness of APR-246 treatment with a sensitivity of 90% and a specificity of 94%. We confirmed the predicted APR-246 sensitivity in 8/10 cell lines and in ex vivo cultures of patient samples. Moreover, the combination of ROS detoxification-impeding APR-246 with approved HDAC-inhibitors, known to elevate ROS, substantially increased APR-246 sensitivity in cell cultures and in vivo in two zebrafish neuroblastoma xenograft models. These data provide evidence that APR-246, in combination with HDAC-inhibitors, displays a novel potent targeted treatment option for neuroblastoma patients.