A mechanistic classification of clinical phenotypes in neuroblastoma

TWIST1 expression is associated with high-risk neuroblastoma and promotes primary and metastatic tumor growth

The embryonic transcription factors TWIST1/2 are frequently overexpressed in cancer, acting as multifunctional oncogenes. Here we investigate their role in neuroblastoma (NB), a heterogeneous childhood malignancy ranging from spontaneous regression to dismal outcomes despite multimodal therapy. We first reveal the association of TWIST1 expression with poor survival and metastasis in primary NB, while TWIST2 correlates with good prognosis. Secondly, suppression of TWIST1 by CRISPR/Cas9 results in a reduction of tumor growth and metastasis colonization in immunocompromised mice. Moreover, TWIST1 knockout tumors display a less aggressive cellular morphology and a reduced disruption of the extracellular matrix (ECM) reticulin network. Additionally, we identify a TWIST1-mediated transcriptional program associated with dismal outcome in NB and involved in the control of pathways mainly linked to the signaling, migration, adhesion, the organization of the ECM, and the tumor cells versus tumor stroma crosstalk. Taken together, our findings confirm TWIST1 as promising therapeutic target in NB.

Retrospective Analysis of INRG Clinical and Genomic Factors for 605 Neuroblastomas in Japan: A Report from the Japan Children’s Cancer Group Neuroblastoma Committee (JCCG-JNBSG)

Neuroblastomas (NBs) exhibit broad and divergent clinical behaviors and tumor risk classification at diagnosis is crucial for the selection of an appropriate therapeutic strategy for each patient. The present study aimed to validate the clinical relevance of International Neuroblastoma Risk Group (INRG) prognostic and genomic markers in a Japanese NB cohort using a retrospective analysis. Follow-up data based on 30 common INRG queries in 605 NB cases diagnosed in Japan between 1990 and 2014 were collected and the genome signature of each tumor sample was integrated. As previously indicated, age, tumor stage, MYCN, DNA ploidy, the adrenals as the primary tumor site, serum ferritin and lactate dehydrogenase (LDH) levels, segmental chromosome aberrations, and the number of chromosome breakpoints (BP) correlated with lower survival rates, while the thorax as the primary tumor site and numerical chromosome aberrations correlated with a favorable prognosis. In the patient group with stage 4, MYCN non-amplified tumors (n = 225), one of the challenging subsets for risk stratification, age ≥ 18 months, LDH ≥ 1400 U/L, and BP ≥ 7 correlated with lower overall and event-free survival rates (p < 0.05). The genome subgroup GG-P2s (partial chromosome gain/loss type with 1p/11q losses and 17q gain, n = 30) was strongly associated with a lower overall survival rate (5-year survival rate: 34%, p < 0.05). Therefore, the combination of the tumor genomic pattern (GG-P2s and BP ≥ 7) with age at diagnosis and LDH will be a promising predictor for MYCN-non-amplified high-risk NBs in patient subsets, in accordance with previous findings from the INRG project.

Integrated analysis of the functions and prognostic values of RNA-binding proteins in neuroblastoma

Background

Neuroblastoma (NB) is the most common solid tumor in children. NB treatment has made significant progress; however, given the high degree of heterogeneity, basic research findings and their clinical application to NB still face challenges. Herein, we identify novel prognostic models for NB.

Methods

We obtained RNA expression data of NB and normal nervous tissue from TARGET and GTEx databases and determined the differential expression patterns of RNA binding protein (RBP) genes between normal and cancerous tissues. Lasso regression and Cox regression analyses identified the five most important differentially expressed genes and were used to construct a new prognostic model. The function and prognostic value of these RBPs were systematically studied and the predictive accuracy verified in an independent dataset.

Results

In total, 348 differentially expressed RBPs were identified. Of these, 166 were up-regulated and 182 down-regulated RBPs. Two hubs RBPs (CPEB3 and CTU1) were identified as prognostic-related genes and were chosen to build the prognostic risk score models. Multivariate Cox analysis was performed on genes from univariate Cox regression and Lasso regression analysis using proportional hazards regression model. A five gene prognostic model: Risk score = (-0.60901*expCPEB3)+(0.851637*expCTU1) was built. Based on this model, the overall survival of patients in the high-risk subgroup was lower (P = 2.152e-04). The area under the curve (AUC) of the receiver-operator characteristic curve of the prognostic model was 0.720 in the TARGET cohort. There were significant differences in the survival rate of patients in the high and low-risk subgroups in the validation data set GSE85047 (P = 0.1237e-08), with the AUC 0.730. The risk model was also regarded as an independent predictor of prognosis (HR = 1.535, 95% CI = 1.368–1.722, P = 2.69E-13).

Conclusions

This study identified a potential risk model for prognosis in NB using Cox regression analysis. RNA binding proteins (CPEB3 and CTU1) can be used as molecular markers of NB.

From DNA Copy Number Gains and Tumor Dependencies to Novel Therapeutic Targets for High-Risk Neuroblastoma

Neuroblastoma is a pediatric tumor arising from the sympatho-adrenal lineage and a worldwide leading cause of childhood cancer-related deaths. About half of high-risk patients die from the disease while survivors suffer from multiple therapy-related side-effects. While neuroblastomas present with a low mutational burden, focal and large segmental DNA copy number aberrations are highly recurrent and associated with poor survival. It can be assumed that the affected chromosomal regions contain critical genes implicated in neuroblastoma biology and behavior. More specifically, evidence has emerged that several of these genes are implicated in tumor dependencies thus potentially providing novel therapeutic entry points. In this review, we briefly review the current status of recurrent DNA copy number aberrations in neuroblastoma and provide an overview of the genes affected by these genomic variants for which a direct role in neuroblastoma has been established. Several of these genes are implicated in networks that positively regulate MYCN expression or stability as well as cell cycle control and apoptosis. Finally, we summarize alternative approaches to identify and prioritize candidate copy-number driven dependency genes for neuroblastoma offering novel therapeutic opportunities.

Variant profiling of colorectal adenomas from three patients of two families with MSH3-related adenomatous polyposis

The spectrum of somatic genetic variation in colorectal adenomas caused by biallelic pathogenic germline variants in the MSH3 gene, was comprehensively analysed to characterise mutational signatures and identify potential driver genes and pathways of MSH3-related tumourigenesis. Three patients from two families with MSH3-associated polyposis were included. Whole exome sequencing of nine adenomas and matched normal tissue was performed. The amount of somatic variants in the MSH3-deficient adenomas and the pattern of single nucleotide variants (SNVs) was similar to sporadic adenomas, whereas the fraction of small insertions/deletions (indels) (21-42% of all small variants) was significantly higher. Interestingly, pathogenic somatic APC variants were found in all but one adenoma. The vast majority (12/13) of these were di-, tetra-, or penta-base pair (bp) deletions. The fraction of APC indels was significantly higher than that reported in patients with familial adenomatous polyposis (FAP) (p < 0.01) or in sporadic adenomas (p < 0.0001). In MSH3-deficient adenomas, the occurrence of APC indels in a repetitive sequence context was significantly higher than in FAP patients (p < 0.01). In addition, the MSH3-deficient adenomas harboured one to five (recurrent) somatic variants in 13 established or candidate driver genes for early colorectal carcinogenesis, including ACVR2A and ARID genes. Our data suggest that MSH3-related colorectal carcinogenesis seems to follow the classical APC-driven pathway. In line with the specific function of MSH3 in the mismatch repair (MMR) system, we identified a characteristic APC mutational pattern in MSH3-deficient adenomas, and confirmed further driver genes for colorectal tumourigenesis.

Spatial and temporal intratumour heterogeneity has potential consequences for single biopsy-based neuroblastoma treatment decisions

Intratumour heterogeneity is a major cause of treatment failure in cancer. We present in-depth analyses combining transcriptomic and genomic profiling with ultra-deep targeted sequencing of multiregional biopsies in 10 patients with neuroblastoma, a devastating childhood tumour. We observe high spatial and temporal heterogeneity in somatic mutations and somatic copy-number alterations which are reflected on the transcriptomic level. Mutations in some druggable target genes including ALK and FGFR1 are heterogeneous at diagnosis and/or relapse, raising the issue whether current target prioritization and molecular risk stratification procedures in single biopsies are sufficiently reliable for therapy decisions. The genetic heterogeneity in gene mutations and chromosome aberrations observed in deep analyses from patient courses suggest clonal evolution before treatment and under treatment pressure, and support early emergence of metastatic clones and ongoing chromosomal instability during disease evolution. We report continuous clonal evolution on mutational and copy number levels in neuroblastoma, and detail its implications for therapy selection, risk stratification and therapy resistance.

Reciprocal impacts of telomerase activity and ADRN/MES differentiation state in neuroblastoma tumor biology

Telomere maintenance and tumor cell differentiation have been separately implicated in neuroblastoma malignancy. Their mechanistic connection is unclear. We analyzed neuroblastoma cell lines and morphologic subclones representing the adrenergic (ADRN) and mesenchymal (MES) differentiation states and uncovered sharp differences in their telomere protein and telomerase activity levels. Pharmacologic conversion of ADRN into MES cells elicited consistent and robust changes in the expression of telomere-related proteins. Conversely, stringent down-regulation of telomerase activity triggers the differentiation of ADRN into MES cells, which was reversible upon telomerase up-regulation. Interestingly, the MES differentiation state is associated with elevated levels of innate immunity factors, including key components of the DNA-sensing pathway. Accordingly, MES but not ADRN cells can mount a robust response to viral infections in vitro. A gene expression signature based on telomere and cell lineage-related factors can cluster neuroblastoma tumor samples into predominantly ADRN or MES-like groups, with distinct clinical outcomes. Our findings establish a strong mechanistic connection between telomere and differentiation and suggest that manipulating telomeres may suppress malignancy not only by limiting the tumor growth potential but also by inducing tumor cell differentiation and altering its immunogenicity.

Predicting amplification of MYCN using CpG methylation biomarkers in neuroblastoma

Background: Neuroblastoma is the most common extracranial solid tumor in childhood. Amplification of MYCN in neuroblastoma is a predictor of poor prognosis. Materials and methods: DNA methylation data from the TARGET data matrix were stratified into MYCN amplified and non-amplified groups. Differential methylation analysis, clustering, recursive feature elimination (RFE), machine learning (ML), Cox regression analysis and Kaplan-Meier estimates were performed. Results and Conclusion: 663 CpGs were differentially methylated between the two groups. A total of 25 CpGs were selected by RFE for clustering and ML, and a 100% clustering accuracy was obtained. ML validation on three external datasets produced high accuracy scores of 100%, 97% and 93%. Eight survival-associated CpGs were also identified. Therapeutic interventions may need to be targeted to patient subgroups.

Neuroblastoma signalling models unveil combination therapies targeting feedback-mediated resistance

Very high risk neuroblastoma is characterised by increased MAPK signalling, and targeting MAPK signalling is a promising therapeutic strategy. We used a deeply characterised panel of neuroblastoma cell lines and found that the sensitivity to MEK inhibitors varied drastically between these cell lines. By generating quantitative perturbation data and mathematical modelling, we determined potential resistance mechanisms. We found that negative feedbacks within MAPK signalling and via the IGF receptor mediate re-activation of MAPK signalling upon treatment in resistant cell lines. By using cell-line specific models, we predict that combinations of MEK inhibitors with RAF or IGFR inhibitors can overcome resistance, and tested these predictions experimentally. In addition, phospho-proteomic profiling confirmed the cell-specific feedback effects and synergy of MEK and IGFR targeted treatment. Our study shows that a quantitative understanding of signalling and feedback mechanisms facilitated by models can help to develop and optimise therapeutic strategies. Our findings should be considered for the planning of future clinical trials introducing MEKi in the treatment of neuroblastoma.

NTRK1/TrkA Signaling in Neuroblastoma Cells Induces Nuclear Reorganization and Intra-Nuclear Aggregation of Lamin A/C

Simple Summary

Neuroblastoma (NB) accounts for 15% of all cancer-related deaths of children. While the amplification of the Myc-N proto-oncogene (MYCN) is a major driver of aggressive NB, the expression of the neurotrophin receptor, NTRK1/TrkA, has been shown to be associated with an excellent outcome. MYCN downregulates NTRK1 expression, but it is unknown if the molecular effects of NTRK1 signaling also affect MYCN-induced networks. The aim of this study was to decipher NTRK1 signaling using an unbiased proteome and phosphoproteome approach. To this end, we realized inducible ectopic NTRK1 expression in a NB cell line with MYCN amplification and analyzed the proteomic changes upon NTRK1 activation in a time-dependent manner. In line with the phenotypes observed, NTRK1 activation induced markers of neuronal differentiation and cell cycle arrest. Most prominently, NTRK1 upregulated the expression and phosphorylation of the nuclear lamina component Lamin A/C. Moreover, NTRK1 signaling also induced the aggregation of LMNA within nucleic foci, which accompanies differentiation in other cell types.

Abstract

(1) Background: Neuroblastomas (NBs) are the most common extracranial solid tumors of children. The amplification of the Myc-N proto-oncogene (MYCN) is a major driver of NB aggressiveness, while high expression of the neurotrophin receptor NTRK1/TrkA is associated with mild disease courses. The molecular effects of NTRK1 signaling in MYCN-amplified NB, however, are still poorly understood and require elucidation. (2) Methods: Inducible NTRK1 expression was realized in four NB cell lines with (IMR5, NGP) or without MYCN amplification (SKNAS, SH-SY5Y). Proteome and phosphoproteome dynamics upon NTRK1 activation by its ligand, NGF, were analyzed in a time-dependent manner in IMR5 cells. Target validation by immunofluorescence staining and automated image processing was performed using the three other NB cell lines. (3) Results: In total, 230 proteins and 134 single phosphorylated class I phosphosites were found to be significantly regulated upon NTRK1 activation. Among known NTRK1 targets, Stathmin and the neurosecretory protein VGF were recovered. Additionally, we observed the upregulation and phosphorylation of Lamin A/C (LMNA) that accumulated inside nuclear foci. (4) Conclusions: We provide a comprehensive picture of NTRK1-induced proteome and phosphoproteome dynamics. The phosphorylation of LMNA within nucleic aggregates was identified as a prominent feature of NTRK1 signaling independent of the MYCN status of NB cells.

Intra-Tumour Genetic Heterogeneity and Prognosis in High-Risk Neuroblastoma

Spatial ITH is defined by genomic and biological variations within a tumour acquired by tumour cell evolution under diverse microenvironments, and its role in NB patient prognosis is understudied. In this work, we applied pangenomic techniques to detect chromosomal aberrations in at least two different areas of each tumour and/or in simultaneously obtained solid and liquid biopsies, detecting ITH in the genomic profile of almost 40% of HR-NB. ITH was better detected when comparing one or more tumour pieces and liquid biopsy (50%) than between different tumour pieces (21%). Interestingly, we found that patients with ITH analysed by pangenomic techniques had a significantly better survival rate that those with non-heterogeneous tumours, especially in cases without MYCN amplification. Moreover, all patients in the studied cohort with high ITH (defined as 50% or more genomic aberration differences between areas of a tumour or simultaneously obtained samples) survived after 48 months. These results clearly support analysing at least two solid tumour areas (separately or mixed) and liquid samples to provide more accurate genomic diagnosis, prognosis and therapy options in HR-NB.

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.

Revised Neuroblastoma Risk Classification System: A Report From the Children's Oncology Group

PURPOSE

Treatment planning for children with neuroblastoma requires accurate assessment of prognosis. The most recent Children's Oncology Group (COG) risk classification system used tumor stage as defined by the International Neuroblastoma Staging System. Here, we validate a revised classifier using the International Neuroblastoma Risk Group Staging System (INRGSS) and incorporate segmental chromosome aberrations (SCA) as an additional genomic biomarker.

METHODS

Newly diagnosed patients enrolled on the COG neuroblastoma biology study ANBL00B1 between 2007 and 2017 with known age, International Neuroblastoma Staging System, and INRGSS stage were identified (N = 4,832). Tumor MYCN status, ploidy, SCA status (1p and 11q), and International Neuroblastoma Pathology Classification histology were determined centrally. Survival analyses were performed for combinations of prognostic factors used in COG risk classification according to the prior version 1, and to validate a revised algorithm (version 2).

RESULTS

Most patients with locoregional tumors had excellent outcomes except for those with image-defined risk factors (INRGSS L2) with MYCN amplification (5-year event-free survival and overall survival: 76.3% ± 5.8% and 79.9% ± 5.5%, respectively) or patients age ≥ 18 months with L2 MYCN nonamplified tumors with unfavorable International Neuroblastoma Pathology Classification histology (72.7% ± 5.4% and 82.4% ± 4.6%), which includes the majority of L2 patients with SCA. For patients with stage M (metastatic) and MS (metastatic, special) disease, genomic biomarkers affected risk group assignment for those < 12 months (MYCN) or 12-18 months (MYCN, histology, ploidy, and SCA) of age. In a retrospective analysis of patient outcome, the 5-year event-free survival and overall survival using COG version 1 were low-risk: 89.4% ± 1.1% and 97.9% ± 0.5%; intermediate-risk: 86.1% ± 1.3% and 94.9% ± 0.8%; high-risk: 50.8% ± 1.4% and 61.9% ± 1.3%; and using COG version 2 were low-risk: 90.7% ± 1.1% and 97.9% ± 0.5%; intermediate-risk: 85.1% ± 1.4% and 95.8% ± 0.8%; high-risk: 51.2% ± 1.4% and 62.5% ± 1.3%, respectively.

CONCLUSION

A revised 2021 COG neuroblastoma risk classifier (version 2) that uses the INRGSS and incorporates SCAs has been adopted to prospectively define COG clinical trial eligibility and treatment assignment.

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

Simple Summary

Preclinical analyses identified APR-246 as a potent treatment option for neuroblastoma. However, a specific mode of action, sufficient biomarkers and promising combination partners are still missing. Here, we analyze the susceptibilities of different entities and relate them to gene expression profiles and previously described biomarkers. We propose a gene signature, consisting of 13 genes, as a novel predictive biomarker. Furthermore, we provide evidence that APR-246 directly targets metabolic weaknesses in neuroblastoma cell lines, thus hampering ROS detoxification. This makes APR-246 suitable to be combined with ROS-inducing HDAC inhibitors, a treatment combination that has not been described for neuroblastoma thus far.

Abstract

APR-246 (Eprenetapopt/PRIMA-1^Met) 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.

B cell infiltration is highly associated with prognosis and an immune-infiltrated tumor microenvironment in neuroblastoma

Aim

Neuroblastoma is the most common extracranial solid tumor in children. Recent advances in immunotherapy Approaches, including in neuroblastoma, have shown the important role of the immune system in mounting an effective anti-tumor response. In this study, we aimed to provide a comprehensive investigation of immune cell infiltration in neuroblastoma utilizing a large number of gene expression datasets.

Methods

We inferred immune cell infiltration using an established immune inference method and evaluated the association between immune cell abundance and patient prognosis as well as common chromosomal abnormalities found in neuroblastoma. In addition, we evaluated co-infiltration patterns among distinct immune cell types.

Results

The infiltration of naïve B cells, NK cells, and CD8+ T cells was associated with improved patient prognosis. Naïve B cells were the most consistent indicator of prognosis and associated with an active immune tumor microenvironment. Patients with high B cell infiltration showed high co-infiltration of other immune cell types and the enrichment of immune-related pathways. The presence of high B cell infiltration was associated with both recurrence-free and overall survival, even after adjusting for clinical variables.

Conclusion

In this study, we have provided a comprehensive evaluation of immune cell infiltration in neuroblastoma using gene expression data. We propose an important role for B cells in the neuroblastoma tumor microenvironment and suggest that B cells can be used as a prognostic biomarker to predict recurrence-free and overall survival independently of currently utilized prognostic variables.

Genetic Alterations and Resectability Predict Outcome in Patients with Neuroblastoma Assigned to High-Risk Solely by MYCN Amplification

Simple Summary

Currently, patients with high-risk neuroblastoma are uniformly treated with maximum therapy. This study investigated a high-risk subgroup characterized by the presence of the amplified MYCN oncogene in the tumor regardless of the stage. In contrast to the corresponding high-risk subgroup consisting of patients with metastases and age at diagnosis over 18 months, the investigated subgroup had generally a superior survival chance. However, the detection of mutations of specific genes in the tumor tissue (RAS and p53 pathway including ALK) had a strong, negative impact. These genes should be therefore also investigated in the future. Complete surgical removal of the primary tumor proved to be beneficial for high-risk neuroblastoma patients assigned to the high-risk category solely by MYCN amplification.

Abstract

Background: To identify variables predicting outcome in neuroblastoma patients assigned to the high-risk group solely by the presence of MYCN oncogene amplification (MNA). Methods: Clinical characteristics, genomic information, and outcome of 190 patients solely assigned to high-risk neuroblastoma by MNA were analyzed and compared to 205 patients with stage 4 neuroblastoma aged ≥18 months with MNA (control group). Results: Event-free survival (EFS) and overall survival (OS) at 10 years were 47% (95%-CI 39–54%) and 56% (95%-CI 49–63%), respectively, which was significantly better than EFS and OS of the control group (EFS 25%, 95%-CI 18–31%, p < 0.001; OS 32% 95%-CI 25–39%, p < 0.001). The presence of RAS-/p53-pathway gene alterations was associated with impaired 10-year EFS and OS (19% vs. 55%, and 19% vs. 67%, respectively; both p < 0.001). In time-dependent multivariable analyses, alterations of RAS-/p53-pathway genes and the extent of the best primary tumor resection were the only independent prognostic variables for OS (p < 0.001 and p = 0.011, respectively). Conclusions: Neuroblastoma patients attributed to high risk solely by MYCN amplification have generally a more favorable outcome. Mutations of genes of the RAS and/or p53 pathways and incomplete resection are the main risk factors predicting poor outcome.

ALT neuroblastoma chemoresistance due to telomere dysfunction-induced ATM activation is reversible with ATM inhibitor AZD0156

Cancers overcome replicative immortality by activating either telomerase or an alternative lengthening of telomeres (ALT) mechanism. ALT occurs in ~25% of high-risk neuroblastomas, and progression in patients with ALT neuroblastoma during or after front-line therapy is frequent and often fatal. Temozolomide + irinotecan is commonly used as salvage therapy for neuroblastoma. Patient-derived cell lines and xenografts established from patients with relapsed ALT neuroblastoma demonstrated de novo resistance to temozolomide + irinotecan [SN-38 in vitro, P < 0.05; in vivo mouse event-free survival (EFS), P < 0.0001] vs. telomerase-positive neuroblastomas. We observed that ALT neuroblastoma cells manifested constitutive ataxia-telangiectasia mutated (ATM) activation due to spontaneous telomere dysfunction which was not observed in telomerase-positive neuroblastoma cells. We demonstrated that induction of telomere dysfunction resulted in ATM activation that, in turn, conferred resistance to temozolomide + SN-38 (4.2-fold change in IC50, P < 0.001). ATM knockdown (shRNA) or inhibition using a clinical-stage small-molecule inhibitor (AZD0156) reversed resistance to temozolomide + irinotecan in ALT neuroblastoma cell lines in vitro (P < 0.001) and in four ALT xenografts in vivo (EFS, P < 0.0001). AZD0156 showed modest to no enhancement of temozolomide + irinotecan activity in telomerase-positive neuroblastoma cell lines and xenografts. Ataxia telangiectasia and Rad3 related (ATR) inhibition using AZD6738 did not enhance temozolomide + SN-38 activity in ALT neuroblastoma cells. Thus, ALT neuroblastoma chemotherapy resistance occurs via ATM activation and is reversible with ATM inhibitor AZD0156. Combining AZD0156 with temozolomide + irinotecan warrants clinical testing for neuroblastoma.

Clinical and molecular characterization of patients with stage 4(M) neuroblastoma aged less than 18 months without MYCN amplification

INTRODUCTION

The survival of children with stage 4(M) neuroblastoma without MYCN amplification and below the age of 18 months is considered better than the still dismal outcome of older high-risk neuroblastoma patients. This study analyzes the impact of clinical and molecular characteristics on the long-term outcome.

PATIENTS AND METHODS

Clinical presentation, survival, and recurrence patterns of patients enrolled onto trials NB90, NB97, and NB2004 were retrospectively analyzed. Gene expression signatures based on RNA microarrays (TH10) were investigated if tumor material was available.

RESULTS

Between 1990 and 2015, 177 patients with stage 4(M) MYCN nonamplified neuroblastoma aged less than 18 months at diagnosis were eligible. After a median follow-up of 9.7 years (IQR 5.0, 13.4), the proportions of 10-year event-free survival (EFS) and overall survival (OS) were 73% (95% confidence interval [CI] 67-79%) and 86% (95% CI 80-92%), respectively. Of the 27 neuroblastoma recurrences, 44% occurred in more than one site. Four additional patients presented histologically mature ganglioneuroma at recurrence. Six patients developed a secondary malignancy. The secondary 5-year EFS and OS of the 27 patients with neuroblastoma recurrence were 44% and 59%, respectively. TH10 gene expression signature was not prognostically predictive in the investigated subcohort.

CONCLUSION

The outcome of patients with stage 4(M) neuroblastoma aged less than 18 months is favorable when treated with high-risk or otherwise intensive therapy. The development of secondary malignancies and the potential of maturation to ganglioneuroma call for a controlled stepwise reduction of treatment intensity.

How Do Telomere Abnormalities Regulate the Biology of Neuroblastoma?

Telomere maintenance plays important roles in genome stability and cell proliferation. Tumor cells acquire replicative immortality by activating a telomere-maintenance mechanism (TMM), either telomerase, a reverse transcriptase, or the alternative lengthening of telomeres (ALT) mechanism. Recent advances in the genetic and molecular characterization of TMM revealed that telomerase activation and ALT define distinct neuroblastoma (NB) subgroups with adverse outcomes, and represent promising therapeutic targets in high-risk neuroblastoma (HRNB), an aggressive childhood solid tumor that accounts for 15% of all pediatric-cancer deaths. Patients with HRNB frequently present with widely metastatic disease, with tumors harboring recurrent genetic aberrations (MYCN amplification, TERT rearrangements, and ATRX mutations), which are mutually exclusive and capable of promoting TMM. This review provides recent insights into our understanding of TMM in NB tumors, and highlights emerging therapeutic strategies as potential treatments for telomerase- and ALT-positive tumors.

Identification of Prognostic Genes in Neuroblastoma in Children by Weighted Gene Coexpression Network Analysis

Background

Neuroblastoma is a malignant neuroendocrine tumor from the sympathetic nervous system, the most common extracranial tumor in children. Identifying potential prognostic markers of neuroblastoma can provide clues for early diagnosis, recurrence, and treatment.

Methods

RNA sequence data and clinical features of 147 neuroblastomas were obtained from the TARGET (Therapeutically Applicable Research to Generate Effective Treatments project) database. Application weighted gene coexpression network analysis (WGCNA) was used to construct a free-scale gene coexpression network, to study the interrelationship between its potential modules and clinical features, and to identify hub genes in the module. We performed Lasso regression and Cox regression analyses to identify the three most important genes and develop a new prognostic model. Data from the GSE85047 cohort verified the predictive accuracy of the prognostic model.

Results

14 coexpression modules were constructed using WGCNA. Brown coexpression modules were found to be significantly associated with disease survival status. Multivariate Cox analysis was performed on genes from univariate Cox regression and Lasso regression analyses using the Cox proportional hazards regression model. Finally, we constructed a three-gene prognostic model: risk score = (0.003812659∗CKB) + (−0.152376975∗expDST) + (0.032032815∗expDUT). The prognosis of samples in the high-risk group was significantly poorer than that of samples in the low-risk group (P=1.225e − 06). The risk model was also regarded as an independent predictor of prognosis (HR = 1.632; 95% CI = 1.391–1.934; P < 0.001).

Conclusion

Our study constructed a neuroblastoma coexpressing gene module and identified a prognostic potential risk model for prognosis in neuroblastoma.

Neuroblastoma Risk Assessment and Treatment Stratification with Hybrid Capture-Based Panel Sequencing

For many years, the risk-based therapy stratification of children with neuroblastoma has relied on clinical and molecular covariates. In recent years, genome analysis has revealed further alterations defining risk, tumor biology, and therapeutic targets. The implementation of a robust and scalable method for analyzing traditional and new molecular markers in routine diagnostics is an urgent clinical need. Here, we investigated targeted panel sequencing as a diagnostic approach to analyze all relevant genomic neuroblastoma risk markers in one assay. Our "neuroblastoma hybrid capture sequencing panel" (NB-HCSP) assay employs a technology for the high-coverage sequencing (>1000×) of 55 selected genes and neuroblastoma-relevant genomic regions, which allows for the detection of single nucleotide changes, structural rearrangements, and copy number alterations. We validated our assay by analyzing 15 neuroblastoma cell lines and a cohort of 20 neuroblastomas, for which reference routine diagnostic data and genome sequencing data were available. We observed a high concordance for risk markers identified by the NB-HSCP assay, clinical routine diagnostics, and genome sequencing. Subsequently, we demonstrated clinical applicability of the NB-HCSP assay by analyzing routine clinical samples. We conclude that the NB-HCSP assay may be implemented into routine diagnostics as a single assay that covers all essential covariates for initial neuroblastoma classification, extended risk stratification, and targeted therapy selection.

Discovery of Spatial Peptide Signatures for Neuroblastoma Risk Assessment by MALDI Mass Spectrometry Imaging

Simple Summary

The childhood tumor, neuroblastoma, has a broad clinical presentation. Risk assessment at diagnosis is particularly difficult in molecularly heterogeneous high-risk cases. Here we investigate the potential of imaging mass spectrometry to directly detect intratumor heterogeneity on the protein level in tissue sections. We show that this approach can produce discriminatory peptide signatures separating high- from low- and intermediate-risk tumors, identify 8 proteins aassociated with these signatures and validate two marker proteins using tissue immunostaining that have promise for further basic and translational research in neuroblastoma. We provide proof-of-concept that mass spectrometry-based technology could assist early risk assessment in neuroblastoma and provide insights into peptide signature-based detection of intratumor heterogeneity.

Abstract

Risk classification plays a crucial role in clinical management and therapy decisions in children with neuroblastoma. Risk assessment is currently based on patient criteria and molecular factors in single tumor biopsies at diagnosis. Growing evidence of extensive neuroblastoma intratumor heterogeneity drives the need for novel diagnostics to assess molecular profiles more comprehensively in spatial resolution to better predict risk for tumor progression and therapy resistance. We present a pilot study investigating the feasibility and potential of matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) to identify spatial peptide heterogeneity in neuroblastoma tissues of divergent current risk classification: high versus low/intermediate risk. Univariate (receiver operating characteristic analysis) and multivariate (segmentation, principal component analysis) statistical strategies identified spatially discriminative risk-associated MALDI-based peptide signatures. The AHNAK nucleoprotein and collapsin response mediator protein 1 (CRMP1) were identified as proteins associated with these peptide signatures, and their differential expression in the neuroblastomas of divergent risk was immunohistochemically validated. This proof-of-concept study demonstrates that MALDI-MSI combined with univariate and multivariate analysis strategies can identify spatially discriminative risk-associated peptide signatures in neuroblastoma tissues. These results suggest a promising new analytical strategy improving risk classification and providing new biological insights into neuroblastoma intratumor heterogeneity.

Targeting Telomere Biology in Acute Lymphoblastic Leukemia

Increased cell proliferation is a hallmark of acute lymphoblastic leukemia (ALL), and genetic alterations driving clonal proliferation have been identified as prognostic factors. To evaluate replicative history and its potential prognostic value, we determined telomere length (TL) in lymphoblasts, B-, and T-lymphocytes, and measured telomerase activity (TA) in leukocytes of patients with ALL. In addition, we evaluated the potential to suppress the in vitro growth of B-ALL cells by the telomerase inhibitor imetelstat. We found a significantly lower TL in lymphoblasts (4.3 kb in pediatric and 2.3 kb in adult patients with ALL) compared to B- and T-lymphocytes (8.0 kb and 8.2 kb in pediatric, and 6.4 kb and 5.5 kb in adult patients with ALL). TA in leukocytes was 3.2 TA/C for pediatric and 0.7 TA/C for adult patients. Notably, patients with high-risk pediatric ALL had a significantly higher TA of 6.6 TA/C compared to non-high-risk patients with 2.2 TA/C. The inhibition of telomerase with imetelstat ex vivo led to significant dose-dependent apoptosis of B-ALL cells. These results suggest that TL reflects clonal expansion and indicate that elevated TA correlates with high-risk pediatric ALL. In addition, telomerase inhibition induces apoptosis of B-ALL cells cultured in vitro. TL and TA might complement established markers for the identification of patients with high-risk ALL. Moreover, TA seems to be an effective therapeutic target; hence, telomerase inhibitors, such as imetelstat, may augment standard ALL treatment.

Frequency and Prognostic Impact of ALK Amplifications and Mutations in the European Neuroblastoma Study Group (SIOPEN) High-Risk Neuroblastoma Trial (HR-NBL1)

PURPOSE

In neuroblastoma (NB), the ALK receptor tyrosine kinase can be constitutively activated through activating point mutations or genomic amplification. We studied ALK genetic alterations in high-risk (HR) patients on the HR-NBL1/SIOPEN trial to determine their frequency, correlation with clinical parameters, and prognostic impact.

MATERIALS AND METHODS

Diagnostic tumor samples were available from 1,092 HR-NBL1/SIOPEN patients to determine ALK amplification status (n = 330), ALK mutational profile (n = 191), or both (n = 571).

RESULTS

Genomic ALK amplification (ALKa) was detected in 4.5% of cases (41 out of 901), all except one with MYCN amplification (MNA). ALKa was associated with a significantly poorer overall survival (OS) (5-year OS: ALKa [n = 41] 28% [95% CI, 15 to 42]; no-ALKa [n = 860] 51% [95% CI, 47 to 54], [P < .001]), particularly in cases with metastatic disease. ALK mutations (ALKm) were detected at a clonal level (> 20% mutated allele fraction) in 10% of cases (76 out of 762) and at a subclonal level (mutated allele fraction 0.1%-20%) in 3.9% of patients (30 out of 762), with a strong correlation between the presence of ALKm and MNA (P < .001). Among 571 cases with known ALKa and ALKm status, a statistically significant difference in OS was observed between cases with ALKa or clonal ALKm versus subclonal ALKm or no ALK alterations (5-year OS: ALKa [n = 19], 26% [95% CI, 10 to 47], clonal ALKm [n = 65] 33% [95% CI, 21 to 44], subclonal ALKm (n = 22) 48% [95% CI, 26 to 67], and no alteration [n = 465], 51% [95% CI, 46 to 55], respectively; P = .001). Importantly, in a multivariate model, involvement of more than one metastatic compartment (hazard ratio [HR], 2.87; P < .001), ALKa (HR, 2.38; P = .004), and clonal ALKm (HR, 1.77; P = .001) were independent predictors of poor outcome.

CONCLUSION

Genetic alterations of ALK (clonal mutations and amplifications) in HR-NB are independent predictors of poorer survival. These data provide a rationale for integration of ALK inhibitors in upfront treatment of HR-NB with ALK alterations.

FOXR2 Stabilizes MYCN Protein and Identifies Non-MYCN-Amplified Neuroblastoma Patients With Unfavorable Outcome

Clinical outcomes of patients with neuroblastoma range from spontaneous tumor regression to fatality. Hence, understanding the mechanisms that cause tumor progression is crucial for the treatment of patients. In this study, we show that FOXR2 activation identifies a subset of neuroblastoma tumors with unfavorable outcome and we investigate the mechanism how FOXR2 relates to poor outcome in patients.

Precision medicine in pediatric solid cancers

Despite huge advances in the diagnosis and treatment of pediatric cancers over the past several decades, it remains one of the leading causes of death during childhood in developed countries. The development of new targeted treatments for these diseases has been hampered by two major factors. First, the extremely heterogeneous nature of the types of tumors encountered in this age group, and their fundamental differences from common adult carcinomas, has made it hard to truly get a handle on the complexities of the underlying biology driving tumor growth. Second, a reluctance of the pharmaceutical industry to develop products or trials for this population due to the relatively small size of the 'market', and a too-easy mechanism of obtaining waivers for pediatric development of adult oncology drugs based on disease type rather than mechanism of action, led to significant difficulties in getting access to new drugs. Thankfully, the field has now started to change, both scientifically and from a regulatory perspective, in order to address some of these challenges. In this review, we will examine some of the recent insights into molecular features which make pediatric tumors so unique and how these might represent therapeutic targets; highlight ongoing international initiatives for providing comprehensive, personalized genomic profiling of childhood tumors in a clinically-relevant timeframe, and look briefly at where the field of pediatric precision oncology may be heading in future.

A review of the biological and clinical implications of RAS-MAPK pathway alterations in neuroblastoma

Neuroblastoma is the most common extra-cranial solid tumor in children, representing approximately 8% of all malignant childhood tumors and 15% of pediatric cancer-related deaths. Recent sequencing and transcriptomics studies have demonstrated the RAS-MAPK pathway’s contribution to the development and progression of neuroblastoma. This review compiles up-to-date evidence of this pathway’s involvement in neuroblastoma. We discuss the RAS-MAPK pathway’s general functioning, the clinical implications of its deregulation in neuroblastoma, and current promising therapeutics targeting proteins involved in signaling.

The synergy of BET inhibitors with aurora A kinase inhibitors in MYCN-amplified neuroblastoma is heightened with functional TP53

Amplification of MYCN is a poor prognostic feature in neuroblastoma (NBL) indicating aggressive disease. We and others have shown BET bromodomain inhibitors (BETi) target MYCN indirectly by downregulating its transcription. Here we sought to identify agents that synergize with BETi and to identify biomarkers of resistance. We previously performed a viability screen of ∼1,900 oncology-focused compounds combined with BET bromodomain inhibitors against MYCN-amplified NBL cell lines. Reanalysis of our screening results prominently identified inhibitors of aurora kinase A (AURKAi) to be highly synergistic with BETi. We confirmed the anti-proliferative effects of several BETi+AURKAi combinations in MYCN-amplified NBL cell lines. Compared to single agents, these combinations cooperated to decrease levels of N-myc. We treated both TP53-wild type and mutant, MYCN-amplified cell lines with the BETi JQ1 and the AURKAi Alisertib. The combination had improved efficacy in the TP53-WT context, notably driving apoptosis in both genetic backgrounds. JQ1+Alisertib combination treatment of a MYCN-amplified, TP53-null or TP53-restored genetically engineered mouse model of NBL prolonged survival better than either single agent. This was most profound with TP53 restored, with marked tumor shrinkage and apoptosis induction in response to combination JQ1+Alisertib. BETi+AURKAi in MYCN-amplified NBL, particularly in the context of functional TP53, provided anti-tumor benefits in preclinical models. This combination should be studied more closely in a pediatric clinical trial.

Neuromedin U and neurotensin may promote the development of the tumour microenvironment in neuroblastoma

Stage 4S neuroblastoma, as defined by the International Neuroblastoma Staging System committee (INSS), is known to regress spontaneously and have a more favourable outcome compared with stage 4 tumours. Comparing the molecular differences between these two stages may provide insights into the progression of neuroblastoma. Our study aimed to explore the molecular differences in the tumour microenvironment (TME) between INSS stage 4S and stage 4 tumours to provide an insight into the mechanisms underlying the biological processes of neuroblastoma. We downloaded the datasets GSE120572 and GSE73517 from the GEO database and pre-processed them using the limma package. CIBERSORT deconvolution agorithm was applied to analyse the differences in 22 infiltrating immune leukocyte subsets between the two stages. We used gene ontology (GO) enrichment analysis to determine the biological process (BP) annotation of differentially expressed genes (DEGs) using the online WebGestalt tool. Hub genes were determined in the STRING database and Cytoscape, and the expression of these genes was verified in the Oncomine database. Then these critical genes were performed survival analysis in TARGET database. We further validated the hub genes using a transwell assay and wound healing assay to detect the function of the genes in the neuroblastoma cell line SK-N-BE(2). GO analysis revealed that the 216 DEGs between stage 4S and stage 4 were enriched in aggressive biological processes. Neuromedin U (NMU) and neurotensin (NTS), which were significantly associated with patients’ overall survival rate, were verified to be elevated in stage 4, and to promote the proliferation and invasion of the SK-N-BE(2) cell. Tumour infiltrating leukocyte analysis showed a high infiltration of regulatory T cells and type 2 tumour-associated macrophages in stage 4 but not in stage 4S. Results of gene co-expression correlation, and the results of previous studies, suggest that NMU and NTS may play certain roles in modulating TME, thus facilitating the progression of neuroblastoma.

Molecular Genetics in Neuroblastoma Prognosis

In recent years, much research has been carried out to identify the biological and genetic characteristics of the neuroblastoma (NB) tumor in order to precisely define the prognostic subgroups for improving treatment stratification. This review will describe the major genetic features and the recent scientific advances, focusing on their impact on diagnosis, prognosis, and therapeutic solutions in NB clinical management.

Acceleration or Brakes: Which Is Rational for Cell Cycle-Targeting Neuroblastoma Therapy?

Unrestrained proliferation is a common feature of malignant neoplasms. Targeting the cell cycle is a therapeutic strategy to prevent unlimited cell division. Recently developed rationales for these selective inhibitors can be subdivided into two categories with antithetical functionality. One applies a “brake” to the cell cycle to halt cell proliferation, such as with inhibitors of cell cycle kinases. The other “accelerates” the cell cycle to initiate replication/mitotic catastrophe, such as with inhibitors of cell cycle checkpoint kinases. The fate of cell cycle progression or arrest is tightly regulated by the presence of tolerable or excessive DNA damage, respectively. This suggests that there is compatibility between inhibitors of DNA repair kinases, such as PARP inhibitors, and inhibitors of cell cycle checkpoint kinases. In the present review, we explore alterations to the cell cycle that are concomitant with altered DNA damage repair machinery in unfavorable neuroblastomas, with respect to their unique genomic and molecular features. We highlight the vulnerabilities of these alterations that are attributable to the features of each. Based on the assessment, we offer possible therapeutic approaches for personalized medicine, which are seemingly antithetical, but both are promising strategies for targeting the altered cell cycle in unfavorable neuroblastomas.

Single-cell transcriptomic analyses provide insights into the developmental origins of neuroblastoma

Neuroblastoma is a pediatric tumor of the developing sympathetic nervous system. However, the cellular origin of neuroblastoma has yet to be defined. Here we studied the single-cell transcriptomes of neuroblastomas and normal human developing adrenal glands at various stages of embryonic and fetal development. We defined normal differentiation trajectories from Schwann cell precursors over intermediate states to neuroblasts or chromaffin cells and showed that neuroblastomas transcriptionally resemble normal fetal adrenal neuroblasts. Importantly, neuroblastomas with varying clinical phenotypes matched different temporal states along normal neuroblast differentiation trajectories, with the degree of differentiation corresponding to clinical prognosis. Our work highlights the roles of oncogenic MYCN and loss of TFAP2B in blocking differentiation and may provide the basis for designing therapeutic interventions to overcome differentiation blocks.

Druggable epigenetic suppression of interferon-induced chemokine expression linked to MYCN amplification in neuroblastoma

BACKGROUND

Amplification of the MYCN oncogene is a molecular hallmark of aggressive neuroblastoma (NB), a childhood cancer of the sympathetic nervous system. There is evidence that MYCN promotes a non-inflamed and T-cell infiltration-poor ('cold') tumor microenvironment (TME) by suppressing interferon signaling. This may explain, at least in part, why patients with NB seem to have little benefit from single-agent immune checkpoint blockade (ICB) therapy. Targeting MYCN or its effectors could be a strategy to convert a cold TME into a 'hot' (inflamed) TME and improve the efficacy of ICB therapy.

METHODS

NB transcriptome analyses were used to identify epigenetic drivers of a T-cell infiltration-poor TME. Biological and molecular responses of NB cells to epigenetic drugs and interferon (IFN)-γ exposure were assessed by proliferation assays, immunoblotting, ELISA, qRT-PCR, RNA-seq and ChIP-qPCR as well as co-culture assays with T cells.

RESULTS

We identified H3K9 euchromatic histone-lysine methyltransferases EHMT2 and EHMT1, also known as G9a and GLP, as epigenetic effectors of the MYCN-driven malignant phenotype and repressors of IFN-γ transcriptional responses in NB cells. EHMT inhibitors enhanced IFN-γ-induced expression of the Th1-type chemokines CXCL9 and CXCL10, key factors of T-cell recruitment into the TME. In MYCN-amplified NB cells, co-inhibition of EZH2 (enhancer of zeste homologue 2), a H3K27 histone methyltransferase cooperating with EHMTs, was needed for strong transcriptional responses to IFN-γ, in line with histone mark changes at CXCL9 and CXCL10 chemokine gene loci. EHMT and EZH2 inhibitor response gene signatures from NB cells were established as surrogate measures and revealed high EHMT and EZH2 activity in MYCN-amplified high-risk NBs with a cold immune phenotype.

CONCLUSION

Our results delineate a strategy for targeted epigenetic immunomodulation of high-risk NBs, whereby EHMT inhibitors alone or in combination with EZH2 inhibitors (in particular, MYCN-amplified NBs) could promote a T-cell-infiltrated TME via enhanced Th1-type chemokine expression.

An Overview of Long Non-Coding (lnc)RNAs in Neuroblastoma

Neuroblastoma (NB) is a heterogeneous developmental tumor occurring in childhood, which arises from the embryonic sympathoadrenal cells of the neural crest. Although the recent progress that has been done on this tumor, the mechanisms involved in NB are still partially unknown. Despite some genetic aberrations having been identified, the sporadic cases represent the majority. Due to its wide heterogeneity in clinical behavior and etiology, NB represents a challenge in terms of prevention and treatment. Since a definitive therapy is lacking so far, there is an urgent necessity to unveil the molecular mechanisms behind NB onset and progression to develop new therapeutic approaches. Long non-coding RNAs (lncRNAs) are a group of RNAs longer than 200 nucleotides. Whether lncRNAs are destined to become a protein or not, they exert multiple biological functions such as regulating gene expression and functions. In recent decades, different research has highlighted the possible role of lncRNAs in the pathogenesis of many diseases, including cancer. Moreover, lncRNAs may represent potential markers or targets for diagnosis and treatment of diseases. This mini-review aimed to briefly summarize the most recent findings on the involvement of some lncRNAs in NB disease by focusing on their mechanisms of action and possible role in unveiling NB onset and progression.

Exosomes from Plasma of Neuroblastoma Patients Contain Doublestranded DNA Reflecting the Mutational Status of Parental Tumor Cells

Neuroblastoma (NB) is an aggressive infancy tumor, leading cause of death among preschool age diseases. Here we focused on characterization of exosomal DNA (exo-DNA) isolated from plasma cell-derived exosomes of neuroblastoma patients, and its potential use for detection of somatic mutations present in the parental tumor cells. Exosomes are small extracellular membrane vesicles secreted by most cells, playing an important role in intercellular communications. Using an enzymatic method, we provided evidence for the presence of double-stranded DNA in the NB exosomes. Moreover, by whole exome sequencing, we demonstrated that NB exo-DNA represents the entire exome and that it carries tumor-specific genetic mutations, including those occurring on known oncogenes and tumor suppressor genes in neuroblastoma (ALK, CHD5, SHANK2, PHOX2B, TERT, FGFR1, and BRAF). NB exo-DNA can be useful to identify variants responsible for acquired resistance, such as mutations of ALK, TP53, and RAS/MAPK genes that appear in relapsed patients. The possibility to isolate and to enrich NB derived exosomes from plasma using surface markers, and the quick and easy extraction of exo-DNA, gives this methodology a translational potential in the clinic. Exo-DNA can be an attractive non-invasive biomarker for NB molecular diagnostic, especially when tissue biopsy cannot be easily available.

Subcellular Distribution of p53 by the p53-Responsive lncRNA NBAT1 Determines Chemotherapeutic Response in Neuroblastoma

Neuroblastoma has a low mutation rate for the p53 gene. Alternative ways of p53 inactivation have been proposed in neuroblastoma, such as abnormal cytoplasmic accumulation of wild-type p53. However, mechanisms leading to p53 inactivation via cytoplasmic accumulation are not well investigated. Here we show that the neuroblastoma risk-associated locus 6p22.3-derived tumor suppressor NBAT1 is a p53-responsive lncRNA that regulates p53 subcellular levels. Low expression of NBAT1 provided resistance to genotoxic drugs by promoting p53 accumulation in cytoplasm and loss from mitochondrial and nuclear compartments. Depletion of NBAT1 altered CRM1 function and contributed to the loss of p53-dependent nuclear gene expression during genotoxic drug treatment. CRM1 inhibition rescued p53-dependent nuclear functions and sensitized NBAT1-depleted cells to genotoxic drugs. Combined inhibition of CRM1 and MDM2 was even more effective in sensitizing aggressive neuroblastoma cells with p53 cytoplasmic accumulation. Thus, our mechanistic studies uncover an NBAT1-dependent CRM1/MDM2-based potential combination therapy for patients with high-risk neuroblastoma. SIGNIFICANCE: This study shows how a p53-responsive lncRNA mediates chemotherapeutic response by modulating nuclear p53 pathways and identifies a potential treatment strategy for patients with high-risk neuroblastoma.

Oncogenic N-Ras Mitigates Oxidative Stress-Induced Apoptosis of Hematopoietic Stem Cells

Leukemic relapse is believed to be driven by transformed hematopoietic stem cells (HSC) that harbor oncogenic mutations or have lost tumor suppressor function. Recent comprehensive sequencing studies have shown that mutations predicted to activate Ras signaling are highly prevalent in hematologic malignancies and, notably, in refractory and relapsed cases. To better understand what drives this clinical phenomenon, we expressed oncogenic NrasG12D within the hematopoietic system in mice and interrogated its effects on HSC survival. N-RasG12D conferred a survival benefit to HSCs and progenitors following metabolic and genotoxic stress. This effect was limited to HSCs and early progenitors and was independent of autophagy and cell proliferation. N-RasG12D-mediated HSC survival was not affected by inhibition of canonical Ras effectors such as MEK and PI3K. However, inhibition of the noncanonical Ras effector pathway protein kinase C (PKC) ameliorated the protective effects of N-RasG12D. Mechanistically, N-RasG12D lowered levels of reactive oxygen species (ROS), which correlated with reduced mitochondrial membrane potential and ATP levels. Inhibition of PKC restored the levels of ROS to that of control HSCs and abrogated the protective effects granted by N-RasG12D. Thus, N-RasG12D activation within HSCs promotes cell survival through the mitigation of ROS, and targeting this mechanism may represent a viable strategy to induce apoptosis during malignant transformation of HSCs. SIGNIFICANCE: Targeting oncogenic N-Ras-mediated reduction of ROS in hematopoietic stem cells through inhibition of the noncanonical Ras effector PKC may serve as a novel strategy for treatment of leukemia and other Ras-mutated cancers.

Alternative lengthening of telomeres in childhood neuroblastoma from genome to proteome

Telomere maintenance by telomerase activation or alternative lengthening of telomeres (ALT) is a major determinant of poor outcome in neuroblastoma. Here, we screen for ALT in primary and relapsed neuroblastomas (n = 760) and characterize its features using multi-omics profiling. ALT-positive tumors are molecularly distinct from other neuroblastoma subtypes and enriched in a population-based clinical sequencing study cohort for relapsed cases. They display reduced ATRX/DAXX complex abundance, due to either ATRX mutations (55%) or low protein expression. The heterochromatic histone mark H3K9me3 recognized by ATRX is enriched at the telomeres of ALT-positive tumors. Notably, we find a high frequency of telomeric repeat loci with a neuroblastoma ALT-specific hotspot on chr1q42.2 and loss of the adjacent chromosomal segment forming a neo-telomere. ALT-positive neuroblastomas proliferate slowly, which is reflected by a protracted clinical course of disease. Nevertheless, children with an ALT-positive neuroblastoma have dismal outcome.

PCAT: an integrated portal for genomic and preclinical testing data of pediatric cancer patient-derived xenograft models

Although cancer is the leading cause of disease-related mortality in children, the relative rarity of pediatric cancers poses a significant challenge for developing novel therapeutics to further improve prognosis. Patient-derived xenograft (PDX) models, which are usually developed from high-risk tumors, are a useful platform to study molecular driver events, identify biomarkers and prioritize therapeutic agents. Here, we develop PDX for Childhood Cancer Therapeutics (PCAT), a new integrated portal for pediatric cancer PDX models. Distinct from previously reported PDX portals, PCAT is focused on pediatric cancer models and provides intuitive interfaces for querying and data mining. The current release comprises 324 models and their associated clinical and genomic data, including gene expression, mutation and copy number alteration. Importantly, PCAT curates preclinical testing results for 68 models and 79 therapeutic agents manually collected from individual agent testing studies published since 2008. To facilitate comparisons of patterns between patient tumors and PDX models, PCAT curates clinical and molecular data of patient tumors from the TARGET project. In addition, PCAT provides access to gene fusions identified in nearly 1000 TARGET samples. PCAT was built using R-shiny and MySQL. The portal can be accessed at http://pcat.zhenglab.info or http://www.pedtranscriptome.org.

MYCN-Amplified Neuroblastoma Is Addicted to Iron and Vulnerable to Inhibition of the System Xc-/Glutathione Axis

MYCN is amplified in 20% to 25% of neuroblastoma, and MYCN-amplified neuroblastoma contributes to a large percent of pediatric cancer-related deaths. Therapy improvements for this subtype of cancer are a high priority. Here we uncover a MYCN-dependent therapeutic vulnerability in neuroblastoma. Namely, amplified MYCN rewires the cell through expression of key receptors, ultimately enhancing iron influx through increased expression of the iron import transferrin receptor 1. Accumulating iron causes reactive oxygen species (ROS) production, and MYCN-amplified neuroblastomas show enhanced reliance on the system Xc- cystine/glutamate antiporter for ROS detoxification through increased transcription of this receptor. This dependence creates a marked vulnerability to targeting the system Xc-/glutathione (GSH) pathway with ferroptosis inducers. This reliance can be exploited through therapy with FDA-approved rheumatoid arthritis drugs sulfasalazine (SAS) and auranofin: in MYCN-amplified, patient-derived xenograft models, both therapies blocked growth and induced ferroptosis. SAS and auranofin activity was largely mitigated by the ferroptosis inhibitor ferrostatin-1, antioxidants like N-acetyl-L-cysteine, or by the iron scavenger deferoxamine (DFO). DFO reduced auranofin-induced ROS, further linking increased iron capture in MYCN-amplified neuroblastoma to a therapeutic vulnerability to ROS-inducing drugs. These data uncover an oncogene vulnerability to ferroptosis caused by increased iron accumulation and subsequent reliance on the system Xc-/GSH pathway. SIGNIFICANCE: This study shows how MYCN increases intracellular iron levels and subsequent GSH pathway activity and demonstrates the antitumor activity of FDA-approved SAS and auranofin in patient-derived xenograft models of MYCN-amplified neuroblastoma.

Association of high-risk neuroblastoma classification based on expression profiles with differentiation and metabolism

Neuroblastoma, the most common extracranial solid malignancy among children, originates from undifferentiated neural crest cells (NCC). Despite recent intensified treatment, high-risk patients still have a high mortality rate. To explore a new therapeutic strategy, we performed an integrated genomic and transcriptomic analysis of 30 high-risk neuroblastoma cases. Based on the expression profiling of RNA sequencing, neuroblastoma was classified into Mesenchymal (MES; n = 5) and Noradrenergic (ADRN; n = 25) clusters, as previously reported in the super-enhancer landscape. The expression patterns in MES-cluster cases were similar to normal adrenal glands, with enrichment in secretion-related pathways, suggesting chromaffin cell-like features built from NCC-derived Schwann cell precursors (SCPs). In contrast, neuron-related pathways were enriched in the ADRN-cluster, indicating sympathoblast features reported to originate from NCC but not via SCPs. Thus, MES- and ADRN-clusters were assumed to be corresponding to differentiation pathways through SCP and sympathoblast, respectively. ADRN-cluster cases were further classified into MYCN- and ATRX-clusters, characterized by genetic alterations, MYCN amplifications and ATRX alterations, respectively. MYCN-cluster cases showed high expression of ALDH18A1, encoding P5CS related to proline production. As reported in other cancers, this might cause reprogramming of proline metabolism leading to tumor specific proline vulnerability candidate for a target therapy of metabolic pathway. In ATRX-cluster, SLC18A2 (VMAT2), an enzyme known to prevent cell toxicity due to the oxidation of dopamine, was highly expressed and VMAT2 inhibitor (GZ-793A) represented significant attenuation of cell growth in NB-69 cell line (high SLC18A2 expression, no MYCN amplification) but not in IMR-32 cell line (MYCN amplification). In addition, the correlation of VMAT2 expression with metaiodobenzylguanidine (MIBG) avidity suggested a combination of VMAT2 inhibitor and MIBG radiation for a novel potential therapeutic strategy in ATRX-cluster cases. Thus, targeting the characteristics of unique neuroblastomas may prospectively improve prognosis.

Integrated analysis of telomerase enzymatic activity unravels an association with cancer stemness and proliferation

Active telomerase is essential for stem cells and most cancers to maintain telomeres. The enzymatic activity of telomerase is related but not equivalent to the expression of TERT, the catalytic subunit of the complex. Here we show that telomerase enzymatic activity can be robustly estimated from the expression of a 13-gene signature. We demonstrate the validity of the expression-based approach, named EXTEND, using cell lines, cancer samples, and non-neoplastic samples. When applied to over 9,000 tumors and single cells, we find a strong correlation between telomerase activity and cancer stemness. This correlation is largely driven by a small population of proliferating cancer cells that exhibits both high telomerase activity and cancer stemness. This study establishes a computational framework for quantifying telomerase enzymatic activity and provides new insights into the relationships among telomerase, cancer proliferation, and stemness.

Neuroblastoma and the epigenome

Neuroblastoma (NB) is a pediatric cancer of the sympathetic nervous system and one of the most common solid tumors in infancy. Amplification of MYCN, copy number alterations, numerical and segmental chromosomal aberrations, mutations, and rearrangements on a handful of genes, such as ALK, ATRX, TP53, RAS/MAPK pathway genes, and TERT, are attributed as underlying causes that give rise to NB. However, the heterogeneous nature of the disease-along with the relative paucity of recurrent somatic mutations-reinforces the need to understand the interplay of genetic factors and epigenetic alterations in the context of NB. Epigenetic mechanisms tightly control gene expression, embryogenesis, imprinting, chromosomal stability, and tumorigenesis, thereby playing a pivotal role in physio- and pathological settings. The main epigenetic alterations include aberrant DNA methylation, disrupted patterns of posttranslational histone modifications, alterations in chromatin composition and/or architecture, and aberrant expression of non-coding RNAs. DNA methylation and demethylation are mediated by DNA methyltransferases (DNMTs) and ten-eleven translocation (TET) proteins, respectively, while histone modifications are coordinated by histone acetyltransferases and deacetylases (HATs, HDACs), and histone methyltransferases and demethylases (HMTs, HDMs). This article focuses predominately on the crosstalk between the epigenome and NB, and the implications it has on disease diagnosis and treatment.

Super enhancers define regulatory subtypes and cell identity in neuroblastoma

Half of the children diagnosed with neuroblastoma (NB) have high-risk disease, disproportionately contributing to overall childhood cancer-related deaths. In addition to recurrent gene mutations, there is increasing evidence supporting the role of epigenetic deregulation in disease pathogenesis. Yet, comprehensive cis-regulatory network descriptions from NB are lacking. Here, using genome-wide H3K27ac profiles across 60 NBs, covering the different clinical and molecular subtypes, we identified four major super-enhancer-driven epigenetic subtypes and their underlying master regulatory networks. Three of these subtypes recapitulated known clinical groups; namely, MYCN-amplified, MYCN non-amplified high-risk and MYCN non-amplified low-risk NBs. The fourth subtype, exhibiting mesenchymal characteristics, shared cellular identity with multipotent Schwann cell precursors, was induced by RAS activation and was enriched in relapsed disease. Notably, CCND1, an essential gene in NB, was regulated by both mesenchymal and adrenergic regulatory networks converging on distinct super-enhancer modules. Overall, this study reveals subtype-specific super-enhancer regulation in NBs.

Whole-genome sequencing of recurrent neuroblastoma reveals somatic mutations that affect key players in cancer progression and telomere maintenance

Neuroblastoma is the most common and deadly childhood tumor. Relapsed or refractory neuroblastoma has a very poor prognosis despite recent treatment advances. To investigate genomic alterations associated with relapse and therapy resistance, whole-genome sequencing was performed on diagnostic and relapsed lesions together with constitutional DNA from seven children. Sequencing of relapsed tumors indicates somatic alterations in diverse genes, including those involved in RAS-MAPK signaling, promoting cell cycle progression or function in telomere maintenance and immortalization. Among recurrent alterations, CCND1-gain, TERT-rearrangements, and point mutations in POLR2A, CDK5RAP, and MUC16 were shown in ≥ 2 individuals. Our cohort contained examples of converging genomic alterations in primary-relapse tumor pairs, indicating dependencies related to specific genetic lesions. We also detected rare genetic germline variants in DNA repair genes (e.g., BARD1, BRCA2, CHEK2, and WRN) that might cooperate with somatically acquired variants in these patients with highly aggressive recurrent neuroblastoma. Our data indicate the importance of monitoring recurrent neuroblastoma through sequential genomic characterization and that new therapeutic approaches combining the targeting of MAPK signaling, cell cycle progression, and telomere activity are required for this challenging patient group.

Telomere-related Markers for Cancer

Telomeres are structurally nucleoprotein complexes at termini of linear chromosomes and essential to chromosome stability/integrity. In normal human cells, telomere length erodes progressively with each round of cell divisions, which serves as an important barrier to uncontrolled proliferation and malignant transformation. In sharp contrast, telomere maintenance is a key feature of human malignant cells and required for their infinite proliferation and maintenance of other cancer hallmarks as well. Thus, a telomere-based anti-cancer strategy has long been suggested. However, clinically efficient and specific drugs targeting cancer telomere-maintenance have still been in their infancy thus far. To achieve this goal, it is highly necessary to elucidate how exactly cancer cells maintain functional telomeres. In the last two decades, numerous studies have provided profound mechanistic insights, and the identified mechanisms include the aberrant activation of telomerase or the alternative lengthening of telomere pathway responsible for telomere elongation, dysregulation and mutation of telomere-associated factors, and other telomere homeostasis-related signaling nodes. In the present review, these various strategies employed by malignant cells to regulate their telomere length, structure and function have been summarized, and potential implications of these findings in the rational development of telomere-based cancer therapy and other clinical applications for precision oncology have been discussed.

Basic research in childhood cancer: Progress and future directions in China

Childhood cancer is a leading cause of death in children. Some childhood cancers have a particularly high mortality rate. Following the World Health Organization's emphasis on child health, most governments worldwide have taken measures to facilitate childhood cancer research. Thus, the scientific community is showing increasing interest in this area. Chinese government has prominence in building a system for the diagnosis and treatment of childhood cancer, thereby promoting the development of childhood cancer research. This review summarizes the research progress, challenges, and perspectives in childhood cancer, and the increasing contributions of National Natural Science Foundation of China (NSFC) in the past decade (2008-2018).

Refining immunotherapeutic approaches to high-risk neuroblastoma based on tumor genomic profiles

In this issue, Coronado et al. attempt to improve our understanding of the factors affecting the response to immunotherapy in a large subset of high‐risk neuroblastoma with hemizygous deletion of chromosome 11q. By using several computational approaches, the authors study potential transcriptional and post‐transcriptional pathways that may affect the response to immunotherapy and further be leveraged therapeutically in a biomarker‐directed fashion.

Combination of tumor necrosis factor-α and epidermal growth factor induces the adrenergic-to-mesenchymal transdifferentiation in SH-SY5Y neuroblastoma cells

Neuroblastoma, a type of cancer that is common in children, is composed of two genetically clonal but epigenetically distinct cell types: mesenchymal (MES) and adrenergic (ADRN) types, controlled by super-enhancer-associated lineage-specific transcription factor networks. Mesenchymal-type cells are more migratory, resistant to chemotherapy, and prevalent in relapse tumors. Importantly, both cell types spontaneously transdifferentiate into one another, and this interconversion can be induced by genetic manipulations. However, the mechanisms of their spontaneous transdifferentiation and extracellular factors inducing this phenomenon have not yet been elucidated. Using a unique approach involving gene set enrichment analysis, we selected six ADRN and 10 MES candidate factors, possibly inducing ADRN and MES phenotypes, respectively. Treatment with a combination of 10 MES factors clearly induced the MES gene expression profile in ADRN-type SH-SY5Y neuroblastoma cells. Considering the effects on gene expression profile, migration ability, and chemoresistance, a combination of tumor necrosis factor alpha (TNF-α) and epidermal growth factor (EGF) was sufficient to synergistically induce the ADRN-to-MES transdifferentiation in SH-SY5Y cells. In addition, human neuroblastoma cohort analysis revealed that the expression of TNF and EGF receptors was strongly associated with MES gene expression signatures, supporting their important roles in transdifferentiation in vivo. Collectively, we propose a mechanism of neuroblastoma transdifferentiation induced by extracellular growth factors, which can be controlled in clinical situations, providing a new therapeutic possibility.