Genomic Amplifications and Distal 6q Loss: Novel Markers for Poor Survival in High-risk Neuroblastoma Patients

Overexpression of H2AFX gene in neuroblastoma is associated with worse prognosis

BACKGROUND

Neuroblastoma (NB) is the most common solid tumour in childhood, and rises in the sympathetic nervous system. Here, we addressed the in silico analysis of the association between the expression of H2AFX gene involved in DNA damage response, and the survival of a cohort of 786 NB patients.

METHODS

In silico gene expression was retrieved from the publicly available dataset summarised by Cangelosi et al., including 13,696 gene expression profiles of 786 NB tumours at onset of disease. The prognostic value of H2AFX (H2A histone family member X) gene expression for event-free survival (EFS) and overall survival (OS) was evaluated by Kaplan-Meier and Cox regression analysis. The main results were validated on another openly accessible in silico database (NRC-283) containing 13,489 gene expressions in 283 NB patients. The expression of H2AFX protein was then tested by immunofluorescence on 48 primary NB samples of different tumour stages. H2AFX activity as an oncogene has been further validated in vitro by silencing the molecule in two NB cell lines, characterised by MYCN amplified or not, and performing cell growth and migration assays.

RESULTS

A strong inverse association between H2AFX expression and patients' survival was observed and confirmed by immunofluorescence results on primary NB tissue sections. Cox regression analysis also disclosed H2AFX as an independent predictor of EFS and OS. The gene-silencing experiments strongly suggested an oncogenic role for H2AFX on NB cells, regardless of MYCN amplification.

CONCLUSIONS

H2AFX is a prognostic marker for unfavourable NB and could be considered a target for therapeutic interventions.

MOXD1 is a lineage-specific gene and a tumor suppressor in neuroblastoma

Neuroblastoma is a childhood developmental cancer; however, its embryonic origins remain poorly understood. Moreover, in-depth studies of early tumor-driving events are limited because of the lack of appropriate models. Herein, we analyzed RNA sequencing data obtained from human neuroblastoma samples and found that loss of expression of trunk neural crest-enriched gene MOXD1 associates with advanced disease and worse outcome. Further, by using single-cell RNA sequencing data of human neuroblastoma cells and fetal adrenal glands and creating in vivo models of zebrafish, chick, and mouse, we show that MOXD1 is a determinate of tumor development. In addition, we found that MOXD1 expression is highly conserved and restricted to mesenchymal neuroblastoma cells and Schwann cell precursors during healthy development. Our findings identify MOXD1 as a lineage-restricted tumor-suppressor gene in neuroblastoma, potentiating further stratification of these tumors and development of novel therapeutic interventions.

Identification of an Ultra-High-Risk Subgroup of Neuroblastoma Patients within the High-Risk Cohort Using a Computed Tomography-Based Radiomics Approach

RATIONALE AND OBJECTIVES

To identify ultra-high-risk (UHR) neuroblastoma patients who experienced disease-related mortality within 18 months of diagnosis within the high-risk cohort using computed tomography (CT)-based radiomics analysis.

MATERIALS AND METHODS

A retrospective analysis was conducted on 105 high-risk neuroblastoma patients, divided into a training set (n = 74) and a test set (n = 31). Radiomics features were extracted and selected from arterial phase CT images, and an optimal radiomics signature was established using the support vector machine algorithm. Evaluation metrics, including area under the curve (AUC) and 95% confidence interval (CI), were calculated. Furthermore, the fit and clinical benefit of the signature, along with its correlation with overall survival (OS), were analyzed.

RESULTS

The optimal radiomics signature comprised 11 features. In the training set, AUC and accuracy were 0.911 (95% CI: 0.840-0.982) and 0.892, respectively. In the test set, AUC and accuracy were 0.828 (95% CI: 0.669-0.987) and 0.839, respectively. There was no significant difference between predicted probability and actual probability, and the signature demonstrated net benefit. The concordance index of this signature for predicting OS was 0.743 (95% CI: 0.672-0.814) in the training set and 0.688 (95% CI: 0.566-0.810) in the test set. Moreover, the signature achieved AUC values of 0.832, 0.863, and 0.721 for 1-year, 2-year, and 3-year OS in the training set, and 0.870, 0.836, and 0.638 in the test set for the respective time periods.

CONCLUSION

The utilization of CT-based radiomics signature to identify an UHR subgroup of neuroblastoma patients within the high-risk cohort can help aid in predicting early disease progression.

Passenger Gene Coamplifications Create Collateral Therapeutic Vulnerabilities in Cancer

DNA amplifications in cancer do not only harbor oncogenes. We sought to determine whether passenger coamplifications could create collateral therapeutic vulnerabilities. Through an analysis of >3,000 cancer genomes followed by the interrogation of CRISPR-Cas9 loss-of-function screens across >700 cancer cell lines, we determined that passenger coamplifications are accompanied by distinct dependency profiles. In a proof-of-principle study, we demonstrate that the coamplification of the bona fide passenger gene DEAD-Box Helicase 1 (DDX1) creates an increased dependency on the mTOR pathway. Interaction proteomics identified tricarboxylic acid (TCA) cycle components as previously unrecognized DDX1 interaction partners. Live-cell metabolomics highlighted that this interaction could impair TCA activity, which in turn resulted in enhanced mTORC1 activity. Consequently, genetic and pharmacologic disruption of mTORC1 resulted in pronounced cell death in vitro and in vivo. Thus, structurally linked coamplification of a passenger gene and an oncogene can result in collateral vulnerabilities.

SIGNIFICANCE

We demonstrate that coamplification of passenger genes, which were largely neglected in cancer biology in the past, can create distinct cancer dependencies. Because passenger coamplifications are frequent in cancer, this principle has the potential to expand target discovery in oncology. This article is featured in Selected Articles from This Issue, p. 384.

Prognostic and predictive biomarkers in paediatric solid tumours

Characterisation of histological, immunohistochemical and molecular prognostic and predictive biomarkers has contributed significantly to precision medicine and better outcomes in the management of paediatric solid tumours. Prognostic biomarkers allow predictions to be made regarding a tumour's aggressiveness and clinical course, whereas predictive biomarkers help determine responses to a specific treatment. This review summarises prognostic biomarkers currently used in the more common paediatric solid tumours, with a brief commentary on the most relevant less common predictive biomarkers. MYCN amplification is the most important genetic alteration in neuroblastoma prognosis, and the histological classification devised by Shimada in 1999 is still used in routine diagnosis. Moreover, a new subgrouping of unfavourable histology neuroblastoma enables immunohistochemical characterisation of tumours with markedly different genetic features and prognosis. The predominant histology and commonly observed cytogenetic abnormalities are recognised outcome predictors in Wilms tumour. Evaluation for anaplasia, which is tightly associated with TP53 gene mutations and poor outcomes, is central in both the International Society of Paediatric Oncology and the Children's Oncology Group approaches to disease classification. Characterisation of distinct genotype-phenotype subclasses and critical mutations has expanded overall understanding of hepatoblastoma outcomes. The C1 subclass hepatoblastoma and CTNNB1 mutations are associated with good prognosis. In contrast, the C2 subclass, NFE2L2 mutations, TERT promoter mutations and high expression of oncofetal proteins and stem cell markers are associated with poor outcomes. Risk stratification in sarcomas is highly variable depending on the entity. The prognosis of rhabdomyosarcoma, for example, primarily depends on histological and molecular characteristics. Advances in our understanding of clinically significant biomarkers will translate into more precise diagnoses, improved risk stratification and more effective and less toxic treatment in this challenging group of patients.

Integration of clinical characteristics and molecular signatures of the tumor microenvironment to predict the prognosis of neuroblastoma

This study aimed to analyze the clinical characteristics, cell types, and molecular characteristics of the tumor microenvironment to better predict the prognosis of neuroblastoma (NB). The gene expression data and corresponding clinical information of 498 NB patients were obtained from the Gene Expression Omnibus (GEO: GSE62564) and ArrayExpress (accession: E-MTAB-8248). The relative cell abundances were estimated using single-sample gene set enrichment analysis (ssGSEA) with the R gene set variation analysis (GSVA) package. We performed Cox regression analyses to identify marker genes indicating cell subsets and combined these with prognostically relevant clinical factors to develop a new prognostic model. Data from the E-MTAB-8248 cohort verified the predictive accuracy of the prognostic model. Single-cell RNA-seq data were analyzed by using the R Seurat package. Multivariate survival analysis for each gene, using clinical characteristics as cofactors, identified 34 prognostic genes that showed a significant correlation with both event-free survival (EFS) and overall survival (OS) (log-rank test, P value < 0.05). The pathway enrichment analysis revealed that these prognostic genes were highly enriched in the marker genes of NB cells with mesenchymal features and protein translation. Ultimately, USP39, RPL8, IL1RAPL1, MAST4, CSRP2, ATP5E, International Neuroblastoma Staging System (INSS) stage, age, and MYCN status were selected to build an optimized Cox model for NB risk stratification. These samples were divided into two groups using the median of the risk score as a cutoff. The prognosis of samples in the poor prognosis group (PP) was significantly worse than that of samples in the good prognosis group (GP) (log-rank test, P value < 0.0001, median EFS: 640.5 vs. 2247 days, median OS: 1279.5 vs. 2519 days). The risk model was also regarded as a prognostic indicator independent of MYCN status, age, and stage. Finally, through scRNA-seq data, we found that as an important prognostic marker, USP39 might participate in the regulation of RNA splicing in NB. Our study established a multivariate Cox model based on gene signatures and clinical characteristics to better predict the prognosis of NB and revealed that mesenchymal signature genes of NB cells, especially USP39, were more abundant in patients with a poor prognosis than in those with a good prognosis. KEY MESSAGES: Our study established a multivariate Cox model based on gene signatures and clinical characteristics to better predict the prognosis of NB and revealed that mesenchymal signature genes of NB cells, especially USP39, were more abundant in patients with a poor prognosis than in those with a good prognosis. USP39, RPL8, IL1RAPL1, MAST4, CSRP2, ATP5E, International Neuroblastoma Staging System (INSS) stage, age, and MYCN status were selected to build an optimized Cox model for NB risk stratification. These samples were divided into two groups using the median of the risk score as a cutoff. The prognosis of samples in the poor prognosis group (PP) was significantly worse than that of samples in the good prognosis group (GP). Finally, through scRNA-seq data, we found that as an important prognostic marker, USP39 might participate in the regulation of RNA splicing in NB.

Use of Optical Genome Mapping to Detect Structural Variants in Neuroblastoma

BACKGROUND

Neuroblastoma is the most common extracranial solid tumour in children, accounting for 15% of paediatric cancer deaths. Multiple genetic abnormalities have been identified as prognostically significant in neuroblastoma patients. Optical genome mapping (OGM) is a novel cytogenetic technique used to detect structural variants, which has not previously been tested in neuroblastoma. We used OGM to identify copy number and structural variants (SVs) in neuroblastoma which may have been missed by standard cytogenetic techniques.

METHODS

Five neuroblastoma cell lines (SH-SY5Y, NBLW, GI-ME-N, NB1691 and SK-N-BE2(C)) and two neuroblastoma tumours were analysed using OGM with the Bionano Saphyr® instrument. The results were analysed using Bionano Access software and compared to previous genetic analyses including G-band karyotyping, FISH (fluorescent in situ hybridisation), single-nucleotide polymorphism (SNP) array and RNA fusion panels for cell lines, and SNP arrays and whole genome sequencing (WGS) for tumours.

RESULTS

OGM detected copy number abnormalities found using previous methods and provided estimates for absolute copy numbers of amplified genes. OGM identified novel SVs, including fusion genes in two cell lines of potential clinical significance.

CONCLUSIONS

OGM can reliably detect clinically significant structural and copy number variations in a single test. OGM may prove to be more time- and cost-effective than current standard cytogenetic techniques for neuroblastoma.

Neuroblastoma Patients' Outcome and Chromosomal Instability

Chromosomal instability (CIN) induces a high rate of losses or gains of whole chromosomes or parts of chromosomes. It is a hallmark of most human cancers and one of the causes of aneuploidy and intra-tumor heterogeneity. The present study aimed to evaluate the potential prognostic role of CIN in NB patients at diagnosis. We performed array comparative genomic hybridization analyses on 451 primary NB patients at the onset of the disease. To assess global chromosomal instability with high precision, we focused on the total number of DNA breakpoints of gains or losses of chromosome arms. For each tumor, an array-CGH-based breakpoint instability index (BPI) was assigned which defined the total number of chromosomal breakpoints per genome. This approach allowed us to quantify CIN related to whole genome disruption in all NB cases analyzed. We found differences in chromosomal breakages among the NB clinical risk groups. High BPI values are negatively associated with survival of NB patients. This association remains significant when correcting for stage, age, and MYCN status in the Cox model. Stratified analysis confirms the prognostic effect of BPI index in low-risk NB patients with non-amplified MYCN and with segmental chromosome aberrations.

Functionalized GD2 Electrochemical Immunosensor to Diagnose Minimum Residual Disease of Bone Marrow in Neuroblastoma Effectively

Neuroblastoma (NB) is known as the "king of childhood tumors" due to its highly metastatic, recurrence-prone, and difficult-to-treat characteristics. International Neuroblastoma Risk Grading Group (INRG) has recommended GD2, a disialoganglioside expressed on neuroectodermal tumor cells, as the target for detecting minimal residual disease in bone marrow metastases of high-risk neuroblastoma in children. Therefore, accurately identifying GD2-positive cells is crucial for diagnosing children with high-risk NB. Here, we designed a graphene/AuNP/GD2 Ab-functionalized electrochemical biosensor for GD2 detection. A three-electrode system was processed using a screen-printed technique with a working electrode of indium tin oxide, a counter electrode of carbon, and a reference electrode of silver/silver chloride. Graphene/AuNPs were modified on the indium tin oxide electrode using chronoamperometric scans, and then, the GD2 antibody was modified on the biosensor by electrostatic adsorption to achieve sensitive and specific detection of GD2-positive cells in bone marrow fluid. The results showed that a graphene/AuNP/GD2 Ab-functionalized electrochemical biosensor achieved GD2-positive cell detection in the range of 102 cells/mL~105 cells/mL by differential pulse voltammetry. Bone marrow fluid samples from 12 children with high-risk NB were retained for testing on our biosensor and showed 100% compliance with the clinical application of the gold-standard immunocytochemical staining technique for detecting GD2-positive cells qualitatively. The GD2-based electrochemical assay can accurately detect children with high-risk NB, providing a rapidly quantitative basis for clinical diagnosis and treatment.

Development and validation of a nomogram based on common biochemical indicators for survival prediction of children with high-risk neuroblastoma: A valuable tool for resource-limited hospitals

BACKGROUND

Despite multiple attempts have been made to develop risk stratification within high-risk neuroblastoma (NB) patients (age of diagnosis ≥ 18 month-old with metastatic NB), the definition of "ultra high-risk NB" is still lack of consensus, and indicators for identifying this subgroup are still unclear. This study aimed to develop a nomogram based on easy-to-obtain blood-derived biofactors for identifying ultra high-risk NB patients with highest risk of death within 3 or 5 years.

METHODS

One hundred sixty-seven NB patients who treated at Sun Yat-sen University Cancer Center between 2015 and 2023 were recruited and clustered randomly into training and validation cohorts (116 and 51 cases, respectively). Univariate and multivariate Cox analysis were performed in training set to screen independent prognostic indicators for constructing nomogram model of predicting 1-, 3- and 5-year overall survival (OS). The discrimination power of the nomogram in training and validation sets were assessed by concordance index (C-index) and calibration plot. Based on the risk score obtained from nomogram model, the prognostic accuracy of 1-, 3- and 5-year OS rates in training and validation cohorts were further evaluated using the area under receiver operating characteristic (ROC) curves (AUC).

RESULTS

Through univariate and multivariate Cox analysis, independent prognostic indicators, including serum lactate dehydrogenase (LDH) and albumin (ALB), were identified in training set, and used to establish a nomogram model. The model showed good discrimination power with C-index in training cohort being 0.706 (95%CI: 0.633-0.788). According to the cut-point calculated based on the established nomogram, patients with a nomogram score > 34 points could be stratified to ultra high-risk NB subgroup, and this subgroup had poorer OS than those in non-ultra one (p < 0.001). AUC values of ROC curves for 3- and 5-year OS rates in the training set were 0.758 and 0.756, respectively. Moreover, based on the cut-point score (34 points) developed in training set, The model also showed good discrimination power with C-index of 0.773 (95%CI: 0.664-0.897) and powerful prognostic accuracy of AUC for 3- and 5-year OS rates being 0.825 and 0.826, respectively, in validation cohort.

CONCLUSIONS

We developed a simple-to-use nomogram based on common laboratory indicators to identify the subgroup of ultra high-risk NB before treatment, providing these children even from developing countries or regions access to intensified multimodal treatments earlier and thus improving their long-term outcome.

A MYCN-independent mechanism mediating secretome reprogramming and metastasis in MYCN-amplified neuroblastoma

MYCN amplification (MNA) is a defining feature of high-risk neuroblastoma (NB) and predicts poor prognosis. However, whether genes within or in close proximity to the MYCN amplicon also contribute to MNA+ NB remains poorly understood. Here, we identify that GREB1, a transcription factor encoding gene neighboring the MYCN locus, is frequently coexpressed with MYCN and promotes cell survival in MNA+ NB. GREB1 controls gene expression independently of MYCN, among which we uncover myosin 1B (MYO1B) as being highly expressed in MNA+ NB and, using a chick chorioallantoic membrane (CAM) model, as a crucial regulator of invasion and metastasis. Global secretome and proteome profiling further delineates MYO1B in regulating secretome reprogramming in MNA+ NB cells, and the cytokine MIF as an important pro-invasive and pro-metastatic mediator of MYO1B activity. Together, we have identified a putative GREB1-MYO1B-MIF axis as an unconventional mechanism promoting aggressive behavior in MNA+ NB and independently of MYCN.

Multiple Genes with Potential Tumor Suppressive Activity Are Present on Chromosome 10q Loss in Neuroblastoma and Are Associated with Poor Prognosis

Neuroblastoma (NB) is a tumor affecting the peripheral sympathetic nervous system that substantially contributes to childhood cancer mortality. Despite recent advances in understanding the complexity of NB, the mechanisms determining its progression are still largely unknown. Some recurrent segmental chromosome aberrations (SCA) have been associated with poor survival. However, the prognostic role of most SCA has not yet been investigated. We examined a cohort of 260 NB primary tumors at disease onset for the loss of chromosome 10q, by array-comparative genomic hybridization (a-CGH) and Single Nucleotide Polymorphism (SNP) array and we found that 26 showed 10q loss, while the others 234 displayed different SCA. We observed a lower event-free survival for NB patients displaying 10q loss compared to patients with tumors carrying other SCA. Furthermore, analyzing the region of 10q loss, we identified a cluster of 75 deleted genes associated with poorer outcome. Low expression of six of these genes, above all CCSER2, was significantly correlated to worse survival using in silico data from 786 NB patients. These potential tumor suppressor genes can be partly responsible for the poor prognosis of NB patients with 10q loss.

SOX11 regulates SWI/SNF complex components as member of the adrenergic neuroblastoma core regulatory circuitry

The pediatric extra-cranial tumor neuroblastoma displays a low mutational burden while recurrent copy number alterations are present in most high-risk cases. Here, we identify SOX11 as a dependency transcription factor in adrenergic neuroblastoma based on recurrent chromosome 2p focal gains and amplifications, specific expression in the normal sympatho-adrenal lineage and adrenergic neuroblastoma, regulation by multiple adrenergic specific (super-)enhancers and strong dependency on high SOX11 expression in adrenergic neuroblastomas. SOX11 regulated direct targets include genes implicated in epigenetic control, cytoskeleton and neurodevelopment. Most notably, SOX11 controls chromatin regulatory complexes, including 10 SWI/SNF core components among which SMARCC1, SMARCA4/BRG1 and ARID1A. Additionally, the histone deacetylase HDAC2, PRC1 complex component CBX2, chromatin-modifying enzyme KDM1A/LSD1 and pioneer factor c-MYB are regulated by SOX11. Finally, SOX11 is identified as a core transcription factor of the core regulatory circuitry (CRC) in adrenergic high-risk neuroblastoma with a potential role as epigenetic master regulator upstream of the CRC.

Mitochondria-associated gene expression perturbation predicts clinical outcomes and shows potential for targeted therapy in neuroblastoma

Background

Mitochondria have long been considered a potential target in cancer therapy because malignant cells are known for their altered energy production. However, there is a lack of comprehensive research on the involvement of mitochondria-associated proteins (MAPs) in neuroblastoma (NB), and their potential as therapeutic targets is yet to be fully explored.

Methods

MAP genes were defined based on the protein-coding genes with mitochondrial localization. The mRNA expression patterns and dynamics of MAP genes associated with NB were investigated by integrating publicly available transcriptional profiles at the cellular and tissue levels. Multivariate Cox regression analysis was conducted to reveal the association of MAP genes with the overall survival (OS) and clinical subgroups of NB patients. The single-cell RNA-seq dataset and gene dependency screening datasets were analyzed to reveal the therapeutic potential of targeting MAP genes.

Results

We compiled a total of 1,712 MAP genes. We found the global and cell type-specific mRNA expression changes of the MAP genes associated with NB status and survival. Our analyses revealed a group of MAP gene signatures independent of MYCN-amplification status associated with NB outcome. We provided computational evidence with selected MAP genes showing good performance in predicting long-term prognosis. By analyzing gene dependency of the MAP genes in NB cell lines and ex vivo human primary T cells, we demonstrated the therapeutic potential of targeting several MAP genes in NB tumors.

Conclusions

Collectively, our study provides evidence for the MAP genes as extended candidates in NB tumor stratification and staging, prognostic prediction, and targeted drug development.

Loss of ALK hotspot mutations in relapsed neuroblastoma

ALK is the most commonly mutated oncogene in neuroblastoma with increased mutation frequency reported at relapse. Here we report the loss of an ALK mutation in two patients at relapse and a paired neuroblastoma cell line at relapse. ALK detection methods including Sanger sequencing, targeted next-generation sequencing and a new ALK Agena MassARRAY technique were used to detect common hotspot ALK variants in tumors at diagnosis and relapse from two high-risk neuroblastoma patients. Copy number analysis including single nucleotide polymorphism array and array comparative genomic hybridization confirmed adequate tumor cell content in DNA used for mutation testing. Case 1 presented with an ALK F1174L mutation at diagnosis with a variant allele frequency (VAF) ranging between 23.5% and 28.5%, but the mutation was undetectable at relapse. Case 2 presented with an ALK R1257Q mutation at diagnosis (VAF = 39%-47.4%) which decreased to <0.01% at relapse. Segmental chromosomal aberrations were maintained between diagnosis and relapse confirming sufficient tumor cell content for mutation detection. The diagnostic SKNBE1n cell line harbors an ALK F1174S mutation, which was lost in the relapsed SKNBE2c cell line. To our knowledge, these are the first reported cases of loss of ALK mutations at relapse in neuroblastoma in the absence of ALK inhibitor therapy, reflecting intra-tumoral spatial and temporal heterogeneity. As ALK inhibitors are increasingly used in the treatment of refractory/relapsed neuroblastoma, our study highlights the importance of confirming whether an ALK mutation detected at diagnosis is still present in clones leading to relapse.

Multi-omics integration reveals a six-malignant cell maker gene signature for predicting prognosis in high-risk neuroblastoma

Background

Neuroblastoma is the most common extracranial solid tumor of childhood, arising from the sympathetic nervous system. High-risk neuroblastoma (HRNB) remains a major therapeutic challenge with low survival rates despite the intensification of therapy. This study aimed to develop a malignant-cell marker gene signature (MMGS) that might serve as a prognostic indicator in HRNB patients.

Methods

Multi-omics datasets, including mRNA expression (single-cell and bulk), DNA methylation, and clinical information of HRNB patients, were used to identify prognostic malignant cell marker genes. MMGS was established by univariate Cox analysis, LASSO, and stepwise multivariable Cox regression analysis. Kaplan–Meier (KM) curve and time-dependent receiver operating characteristic curve (tROC) were used to evaluate the prognostic value and performance of MMGS, respectively. MMGS further verified its reliability and accuracy in the independent validation set. Finally, the characteristics of functional enrichment, tumor immune features, and inflammatory activity between different MMGS risk groups were also investigated.

Results

We constructed a prognostic model consisting of six malignant cell maker genes (MAPT, C1QTNF4, MEG3, NPW, RAMP1, and CDT1), which stratified patients into ultra-high-risk (UHR) and common-high-risk (CHR) group. Patients in the UHR group had significantly worse overall survival (OS) than those in the CHR group. MMGS was verified as an independent predictor for the OS of HRNB patients. The area under the curve (AUC) values of MMGS at 1-, 3-, and 5-year were 0.78, 0.693, and 0.618, respectively. Notably, functional enrichment, tumor immune features, and inflammatory activity analyses preliminarily indicated that the poor prognosis in the UHR group might result from the dysregulation of the metabolic process and immunosuppressive microenvironment.

Conclusion

This study established a novel six-malignant cell maker gene prognostic model that can be used to predict the prognosis of HRNB patients, which may provide new insight for the treatment and personalized monitoring of HRNB patients.

Chromosome instability in neuroblastoma: A pathway to aggressive disease

For over 100-years, genomic instability has been investigated as a central player in the pathogenesis of human cancer. Conceptually, genomic instability includes an array of alterations from small deletions/insertions to whole chromosome alterations, referred to as chromosome instability. Chromosome instability has a paradoxical impact in cancer. In most instances, the introduction of chromosome instability has a negative impact on cellular fitness whereas in cancer it is usually associated with a worse prognosis. One exception is the case of neuroblastoma, the most common solid tumor outside of the brain in children. Neuroblastoma tumors have two distinct patterns of genome instability: whole-chromosome aneuploidy, which is associated with a better prognosis, or segmental chromosomal alterations, which is a potent negative prognostic factor. Through a computational screen, we found that low levels of the de- ubiquitinating enzyme USP24 have a highly significant negative impact on survival in neuroblastoma. At the molecular level, USP24 loss leads to destabilization of the microtubule assembly factor CRMP2 - producing mitotic errors and leading to chromosome missegregation and whole-chromosome aneuploidy. This apparent paradox may be reconciled through a model in which whole chromosome aneuploidy leads to the subsequent development of segmental chromosome alterations. Here we review the mechanisms behind chromosome instability and the evidence for the progressive development of segmental alterations from existing numerical aneuploidy in support of a multi-step model of neuroblastoma progression.

E2F3 gene expression is a potential negative prognostic marker for localised and MYCN not-amplified neuroblastoma: Results of in silico analysis of 786 samples

BACKGROUND

Neuroblastoma (NB) is an enigmatic childhood malignancy characterised by a wide range of clinical behaviour. Many potential oncogenes for NB have recently been identified. Among them, E2 transcription factor 3 (E2F3) expression was associated with a poor survival in 134 stage 4S patients, but evidence for other stage groups remains poorly investigated.

METHODS

We have analysed the expression of E2F3 gene from a database of 786 NB samples. Overall and event-free survivals (EFS) were assessed by the Kaplan-Meier method, splitting the data on the median and tertile expression values. The Cox model was applied to control for the confounding by stage, age and MYCN amplification. Validation was performed by an in silico analysis of an independent cohort of 283 NB patients. Furthermore, an immunofluorescence analysis on 48 formalin-fixed, paraffin-embedded NB specimens was also performed.

RESULTS

E2F3 overexpression was associated with a poor survival (EFS = 84%, 95% CI: 79%-95%, for low expression levels; EFS = 62%, 95% CI: 56%-68% for middle levels; EFS = 30%, 95% CI: 24%-36%, for high levels, p < .001). This association was confirmed in multivariable analysis and was more evident in patients with MYCN not-amplified and localised stages. Immunofluorescence results and the validation on an independent cohort of NB primary samples confirmed these findings.

CONCLUSIONS

E2F3 is a new potential prognostic marker in NB with favourable characteristics at diagnosis. Further studies are needed to elucidate the potential role of E2F3 in NB oncogenesis and progression, in order to identify new targets for therapeutic interventions.

Tumoral heterogeneity in neuroblastoma

Neuroblastoma is a solid, neuroendocrine tumor with divergent clinical behavior ranging from asymptomatic to fatal. The diverse clinical presentations of neuroblastoma are directly linked to the high intra- and inter-tumoral heterogeneity it presents. This heterogeneity is strongly associated with therapeutic resistance and continuous relapses, often leading to fatal outcomes. The development of successful risk assessment and tailored treatment strategies lies in evaluating the extent of heterogeneity via the accurate genetic and epigenetic profiling of distinct cell subpopulations present in the tumor. Recent studies have focused on understanding the molecular mechanisms that drive tumoral heterogeneity in pursuing better therapeutic and diagnostic approaches. This review describes the cellular, genetic, and epigenetic aspects of neuroblastoma heterogeneity. In addition, we summarize the recent findings on three crucial factors that can lead to heterogeneity in solid tumors: the inherent diversity of the progenitor cells, the presence of cancer stem cells, and the influence of the tumor microenvironment.

Amplification of CDK4 and MDM2: a detailed study of a high-risk neuroblastoma subgroup

In neuroblastoma, MYCN amplification and 11q-deletion are important, although incomplete, markers of high-risk disease. It is therefore relevant to characterize additional alterations that can function as prognostic and/or predictive markers. Using SNP-microarrays, a group of neuroblastoma patients showing amplification of one or multiple 12q loci was identified. Two loci containing CDK4 and MDM2 were commonly co-amplified, although amplification of either locus in the absence of the other was observed. Pharmacological inhibition of CDK4/6 with ribociclib or abemaciclib decreased proliferation in a broad set of neuroblastoma cell lines, including CDK4/MDM2-amplified, whereas MDM2 inhibition by Nutlin-3a was only effective in p53^wild-type cells. Combined CDK4/MDM2 targeting had an additive effect in p53^wild-type cell lines, while no or negative additive effect was observed in p53^mutated cells. Most 12q-amplified primary tumors were of abdominal origin, including those of intrarenal origin initially suspected of being Wilms' tumor. An atypical metastatic pattern was also observed with low degree of bone marrow involvement, favoring other sites such as the lungs. Here we present detailed biological data of an aggressive neuroblastoma subgroup hallmarked by 12q amplification and atypical clinical presentation for which our in vitro studies indicate that CDK4 and/or MDM2 inhibition also could be beneficial.

Cellular senescence in neuroblastoma

Neuroblastoma is a tumour that arises from the sympathoadrenal lineage occurring predominantly in children younger than five years. About half of the patients are diagnosed with high-risk tumours and undergo intensive multi-modal therapy. The success rate of current treatments for high-risk neuroblastoma is disappointingly low and survivors suffer from multiple therapy-related long-term side effects. Most chemotherapeutics drive cancer cells towards cell death or senescence. Senescence has long been considered to represent a terminal non-proliferative state and therefore an effective barrier against tumorigenesis. This dogma, however, has been challenged by recent observations that infer a much more dynamic and reversible nature for this process, which may have implications for the efficacy of therapy-induced senescence-oriented treatment strategies. Neuroblastoma cells in a dormant, senescent-like state may escape therapy, whilst their senescence-associated secretome may promote inflammation and invasiveness, potentially fostering relapse. Conversely, due to its distinct molecular identity, senescence may also represent an opportunity for the development of novel (combination) therapies. However, the limited knowledge on the molecular dynamics and diversity of senescence signatures demands appropriate models to study this process in detail. This review summarises the molecular knowledge about cellular senescence in neuroblastoma and investigates current and future options towards therapeutic exploration.

Advancing therapy for neuroblastoma

Neuroblastomas are tumours of sympathetic origin, with a heterogeneous clinical course ranging from localized or spontaneously regressing to widely metastatic disease. Neuroblastomas recapitulate many of the features of sympathoadrenal development, which have been directly targeted to improve the survival outcomes in patients with high-risk disease. Over the past few decades, improvements in the 5-year survival of patients with metastatic neuroblastomas, from <20% to >50%, have resulted from clinical trials incorporating high-dose chemotherapy with autologous stem cell transplantation, differentiating agents and immunotherapy with anti-GD2 monoclonal antibodies. The next generation of trials are designed to improve the initial response rates in patients with high-risk neuroblastomas via the addition of immunotherapies, targeted therapies (such as ALK inhibitors) and radiopharmaceuticals to standard induction regimens. Other trials are focused on testing precision medicine strategies for patients with relapsed and/or refractory disease, enhancing the antitumour immune response and improving the effectiveness of maintenance regimens, in order to prolong disease remission. In this Review, we describe advances in delineating the pathogenesis of neuroblastoma and in identifying the drivers of high-risk disease. We then discuss how this knowledge has informed improvements in risk stratification, risk-adapted therapy and the development of novel therapies.

miR-99b-5p, miR-380-3p, and miR-485-3p are novel chemosensitizing miRNAs in high-risk neuroblastoma

Neuroblastoma is a deadly childhood cancer arising in the developing sympathetic nervous system. High-risk patients are currently treated with intensive chemotherapy, which is curative in only 50% of children and leaves some surviving patients with life-long side effects. microRNAs (miRNAs) are critical regulators of neural crest development and are deregulated during neuroblastoma tumorigenesis, making miRNA-based drugs an attractive therapeutic avenue. A functional screen of >1,200 miRNA mimics was conducted in neuroblastoma cell lines to discover miRNAs that sensitized cells to low doses (30% inhibitory concentration [IC30]) of doxorubicin and vincristine chemotherapy used in the treatment of the disease. Three miRNAs, miR-99b-5p, miR-380-3p, and miR-485-3p, had potent chemosensitizing activity with doxorubicin in multiple models of high-risk neuroblastoma. These miRNAs underwent genomic loss in a subset of neuroblastoma patients, and low expression predicted poor survival outcome. In vitro functional assays revealed each of these miRNAs enhanced the anti-proliferative and pro-apoptotic effects of doxorubicin. We used RNA sequencing (RNA-seq) to show that miR-99b-5p represses neuroblastoma dependency genes LIN28B and PHOX2B both in vitro and in patient-derived xenograft (PDX) tumors. Luciferase reporter assays demonstrate that PHOX2B is a direct target of miR-99b-5p. We anticipate that restoring the function of the tumor-suppressive miRNAs discovered here may be a valuable therapeutic strategy for the treatment of neuroblastoma patients.

Single-cell landscape analysis reveals distinct regression trajectories and novel prognostic biomarkers in primary neuroblastoma

Neuroblastoma (NB), which is the most common pediatric extracranial solid tumor, varies widely in its clinical presentation and outcome. NB has a unique ability to spontaneously differentiate and regress, suggesting a potential direction for therapeutic intervention. However, the underlying mechanisms of regression remain largely unknown, and more reliable prognostic biomarkers are needed for predicting trajectories for NB. We performed scRNA-seq analysis on 17 NB clinical samples and three peritumoral adrenal tissues. Primary NB displayed varied cell constitution, even among tumors of the same pathological subtype. Copy number variation patterns suggested that neuroendocrine cells represent the malignant cell type. Based on the differential expression of sets of related marker genes, a subgroup of neuroendocrine cells was identified and projected to differentiate into a subcluster of benign fibroblasts with highly expressed CCL2 and ZFP36, supporting a progressive pathway of spontaneous NB regression. We also identified prognostic markers (STMN2, TUBA1A, PAGE5, and ETV1) by evaluating intra-tumoral heterogeneity. Lastly, we determined that ITGB1 in M2-like macrophages was associated with favorable prognosis and may serve as a potential diagnostic marker and therapeutic target. In conclusion, our findings reveal novel mechanisms underlying regression and potential prognostic markers and therapeutic targets of NB.

Single-cell transcriptomics of neuroblastoma identifies chemoresistance-associated genes and pathways

High-Risk neuroblastoma (NB) survival rate is still <50%, despite treatments being more and more aggressive. The biggest hurdle liable to cancer therapy failure is the drug resistance by tumor cells that is likely due to the intra-tumor heterogeneity (ITH). To investigate the link between ITH and therapy resistance in NB, we performed a single cell RNA sequencing (scRNAseq) of etoposide and cisplatin resistant NB and their parental cells. Our analysis showed a clear separation of resistant and parental cells for both conditions by identifying 8 distinct tumor clusters in etoposide-resistant/parental and 7 in cisplatin-resistant/parental cells. We discovered that drug resistance can affect NB cell identities; highlighting the bi-directional ability of adrenergic-to-mesenchymal transition of NB cells. The biological processes driving the identified resistant cell subpopulations revealed genes such as (BARD1, BRCA1, PARP1, HISTH1 axis, members of RPL family), suggesting a potential drug resistance due to the acquisition of DNA repair mechanisms and to the modification of the drug targets. Deconvolution analysis of bulk RNAseq data from 498 tumors with cell subpopulation signatures showed that the transcriptional heterogeneity of our cellular models reflected the ITH of NB tumors and allowed the identification of clusters associated with worse/better survival.

Our study demonstrates the distinct cell populations characterized by genes involved in different biological processes can have a role in NB drug treatment failure. These findings evidence the importance of ITH in NB drug resistance studies and the chance that scRNA-seq analysis offers in the identification of genes and pathways liable for drug resistance.

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.

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.

High Grade of Amplification of Six Regions on Chromosome 2p in a Neuroblastoma Patient with Very Poor Outcome: The Putative New Oncogene TSSC1

Simple Summary

Here, a case of neuroblastoma (NB) carrying a high-grade amplification of six loci besides MYCN is described. Since the patient had a very poor outcome, we postulated that these DNA co-amplifications might have a synergistic effect in increasing NB cell proliferation. In order to verify this hypothesis, we analyzed in silico the impact of high expression of the genes located within the amplifications on the NB patients’ outcome using a large dataset integrating three different platforms. These analyses disclosed that high expression of the TSSC1 gene was the most significantly associated with reduced overall survival of NB patients, suggesting that it may have a potential prognostic role in NB in both MYCN amplified and MYCN not amplified tumors. Further studies on TSSC1 interactions and functioning could lead to possible focused therapies for high-risk NB patients.

Abstract

We observed a case of high-risk neuroblastoma (NB) carried by a 28-month-old girl, displaying metastatic disease and a rapid decline of clinical conditions. By array-CGH analysis of the tumor tissue and of the metastatic bone marrow aspirate cells, we found a high-grade amplification of six regions besides MYCN on bands 2p25.3–p24.3. The genes involved in these amplifications were MYT1L, TSSC1, CMPK2, RSAD2, RNF144A, GREB1, NTSR2, LPIN1, NBAS, and the two intergenic non-protein coding RNAs LOC730811 and LOC339788. We investigated if these DNA co-amplifications may have an effect on enhancing tumor aggressiveness. We evaluated the association between the high expression of the amplified genes and NB patient’s outcome using the integration of gene expression data of 786 NB samples profiled with different public platforms from patients with at least five-year follow-up. NB patients with high expression of the TSSC1 gene were associated with a reduced survival rate. Immunofluorescence staining on primary tumor tissues confirmed that the TSSC1 protein expression was high in the relapsed or dead stage 4 cases, but it was generally low in NB patients in complete remission. TSSC1 appears as a putative new oncogene in NB.

Integrated CGH/WES Analyses Advance Understanding of Aggressive Neuroblastoma Evolution: A Case Study

Neuroblastoma (NB) is the most common extra-cranial malignancy in preschool children. To portray the genetic landscape of an overly aggressive NB leading to a rapid clinical progression of the disease, tumor DNA collected pre- and post-treatment has been analyzed. Array comparative genomic hybridization (aCGH), whole-exome sequencing (WES), and pharmacogenetics approaches, respectively, have identified relevant copy number alterations (CNAs), single nucleotide variants (SNVs), and polymorphisms (SNPs) that were then combined into an integrated analysis. Spontaneously formed 3D tumoroids obtained from the recurrent mass have also been characterized. The results prove the power of combining CNAs, SNVs, and SNPs analyses to assess clonal evolution during the disease progression by evidencing multiple clones at disease onset and dynamic genomic alterations during therapy administration. The proposed molecular and cytogenetic integrated analysis empowers the disease follow-up and the prediction of tumor recurrence.

MEIS2 Is an Adrenergic Core Regulatory Transcription Factor Involved in Early Initiation of TH-MYCN-Driven Neuroblastoma Formation

Simple Summary

Neuroblastoma is a pediatric tumor originating from the sympathetic nervous system responsible for 10–15% of all childhood cancer deaths. Half of all neuroblastoma patients present with high-risk disease, of which nearly 50% relapse and die of their disease. In addition, standard therapies cause serious lifelong side effects and increased risk for secondary tumors. Further research is crucial to better understand the molecular basis of neuroblastomas and to identify novel druggable targets. Neuroblastoma tumorigenesis has to this end been modeled in both mice and zebrafish. Here, we present a detailed dissection of the gene expression patterns that underlie tumor formation in the murine TH-MYCN-driven neuroblastoma model. We identified key factors that are putatively important for neuroblastoma tumor initiation versus tumor progression, pinpointed crucial regulators of the observed expression patterns during neuroblastoma development and scrutinized which factors could be innovative and vulnerable nodes for therapeutic intervention.

Abstract

Roughly half of all high-risk neuroblastoma patients present with MYCN amplification. The molecular consequences of MYCN overexpression in this aggressive pediatric tumor have been studied for decades, but thus far, our understanding of the early initiating steps of MYCN-driven tumor formation is still enigmatic. We performed a detailed transcriptome landscaping during murine TH-MYCN-driven neuroblastoma tumor formation at different time points. The neuroblastoma dependency factor MEIS2, together with ASCL1, was identified as a candidate tumor-initiating factor and shown to be a novel core regulatory circuit member in adrenergic neuroblastomas. Of further interest, we found a KEOPS complex member (gm6890), implicated in homologous double-strand break repair and telomere maintenance, to be strongly upregulated during tumor formation, as well as the checkpoint adaptor Claspin (CLSPN) and three chromosome 17q loci CBX2, GJC1 and LIMD2. Finally, cross-species master regulator analysis identified FOXM1, together with additional hubs controlling transcriptome profiles of MYCN-driven neuroblastoma. In conclusion, time-resolved transcriptome analysis of early hyperplastic lesions and full-blown MYCN-driven neuroblastomas yielded novel components implicated in both tumor initiation and maintenance, providing putative novel drug targets for MYCN-driven 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.

Engineering large-scale chromosomal deletions by CRISPR-Cas9

Large-scale chromosomal deletions are a prevalent and defining feature of cancer. A high degree of tumor-type and subtype specific recurrencies suggest a selective oncogenic advantage. However, due to their large size it has been difficult to pinpoint the oncogenic drivers that confer this advantage. Suitable functional genomics approaches to study the oncogenic driving capacity of large-scale deletions are limited. Here, we present an effective technique to engineer large-scale deletions by CRISPR-Cas9 and create isogenic cell line models. We simultaneously induce double-strand breaks (DSBs) at two ends of a chromosomal arm and select the cells that have lost the intermittent region. Using this technique, we induced large-scale deletions on chromosome 11q (65 Mb) and chromosome 6q (53 Mb) in neuroblastoma cell lines. A high frequency of successful deletions (up to 30% of selected clones) and increased colony forming capacity in the 11q deleted lines suggest an oncogenic advantage of these deletions. Such isogenic models enable further research on the role of large-scale deletions in tumor development and growth, and their possible therapeutic potential.

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.

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.

USP24 Is a Cancer-Associated Ubiquitin Hydrolase, Novel Tumor Suppressor, and Chromosome Instability Gene Deleted in Neuroblastoma

Deubiquitinating enzymes are increasingly recognized to play important roles in cancer, with many acting as oncogenes or tumor suppressors. In this study, we employed a bioinformatics approach to screen for enzymes from this family involved in cancer and found USP24 as a potent predictor of poor outcomes in neuroblastoma, an aggressive childhood cancer. USP24 resides in a region commonly deleted in neuroblastoma, yet was independently associated with poor outcomes in this disease. Deletion of Usp24 in a murine model resulted in degradation of collapsin response mediator protein 2 (CRMP2), a regulator of axon growth, guidance, and neuronal polarity. Cells lacking USP24 had significant increases in spindle defects, chromosome missegregation, and aneuploidy, phenotypes that were rescued by the restoration of CRMP2. USP24 prevented aneuploidy by maintaining spindle-associated CRMP2, which is required for mitotic accuracy. Our findings further indicate that USP24 is a tumor suppressor that may play an important role in the pathogenesis of neuroblastoma. SIGNIFICANCE: This study identifies the chromosome instability gene USP24 as frequently deleted in neuroblastoma and provides important insight into the pathogenesis of this aggressive childhood cancer.

A nomogram of clinical and biologic factors to predict survival in children newly diagnosed with high-risk neuroblastoma: An International Neuroblastoma Risk Group project

BACKGROUND

Long-term outcome remains poor for children with high-risk neuroblastoma (five-year overall survival [OS] ∼50%). Our objectives were to (a) identify prognostic biomarkers and apply them in a nomogram to identify the subgroup of ultra-high-risk patients at highest risk of disease progression/death, for whom novel frontline therapy is urgently needed; and (b) validate the nomogram in an independent cohort.

METHODS

A total of 1820 high-risk patients (≥18 months old with metastatic neuroblastoma), diagnosed 1998-2015, from the International Neuroblastoma Risk Groups (INRG) Data Commons were analyzed in a retrospective cohort study. Using multivariable Cox regression of OS from diagnosis, a nomogram was created from prognostic biomarkers to predict three-year OS. External validation was performed using the SIOPEN HR-NBL1 trial cohort (n = 521), evidenced by receiver operating characteristic curves.

RESULTS

The nomogram, including MYCN status (P < 0.0001), lactate dehydrogenase (LDH) (P = 0.0007), and presence of bone marrow metastases (P = 0.004), had robust performance and was validated. Applying the nomogram at diagnosis (a) gives prognosis of an individual patient and (b) identifies patients predicted to have poor outcome (three-year OS was 30% ± 5% for patients with a nomogram score of > 82 points; 58% ± 1% for those ≤82 points). Median follow-up time was 5.5 years (range, 0-14.1).

CONCLUSIONS

In high-risk neuroblastoma, a novel, publicly available nomogram using prognostic biomarkers (MYCN status, LDH, presence of bone marrow metastases; https://neuroblastoma.shinyapps.io/High-Risk-Neuroblastoma-Nomogram/) has the flexibility to apply a clinically suitable and context-specific cutoff to identify patients at highest risk of death. This will facilitate testing urgently needed new frontline treatment options to improve outcome for these children.

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.

Ezrin interacts with the tumor suppressor CHL1 and promotes neuronal differentiation of human neuroblastoma

In a previous study, we demonstrated that CHL1, the neuronal cell adhesion molecule close homolog of L1, acts as a tumor suppressor in human neuroblastoma (NB), a still highly lethal childhood malignancy, influencing its differentiation and proliferation degree. Here we found that ezrin, one of the ERM (ezrin, radixin, moesin) proteins involved in cytoskeleton organization, strongly interacts with CHL1. The low expression of EZRIN, as well as the low expression of CHL1 and of the neuronal differentiation marker MAP2, correlates with poor outcome in NB patients. Knock-down of ezrin in HTLA-230 cell line induces neurite retraction, enhances cell proliferation and migration, and triggers anchorage-independent growth, with effects very similar to those already obtained by CHL1 silencing. Furthermore, lack of ezrin inhibits the expression of MAP2 and of the oncosuppressor molecule p53, whereas it enhances MAPK activation, all typical features of tumor aggressiveness. As already described, CHL1 overexpression in IMR-32 cell line provokes an opposite trend, but the co-silencing of ezrin reduces these effects, confirming the hypothesis that CHL1 acts in close connection with ezrin. Overall, our data show that ezrin reinforces the differentiating and oncosuppressive functions of CHL1, identifying this ERM protein as a new targetable molecule for NB therapy.

Association of image-defined risk factors with clinical features, histopathology, and outcomes in neuroblastoma

BACKGROUND

Clinical, molecular, and histopathologic features guide treatment for neuroblastoma, but obtaining tumor tissue may cause complications and is subject to sampling error due to tumor heterogeneity. We hypothesized that image-defined risk factors (IDRFs) would reflect molecular features, histopathology, and clinical outcomes in neuroblastoma.

METHODS

We performed a retrospective cohort study of 76 patients with neuroblastoma or ganglioneuroblastoma. Diagnostic CT scans were reviewed for 20 IDRFs, which were consolidated into five IDRF groups (involvement of multiple body compartments, vascular encasement, tumor infiltration of adjacent organs/structures, airway compression, or intraspinal extension). IDRF groups were analyzed for association with clinical, molecular, and histopathologic features of neuroblastoma.

RESULTS

Patients with more IDRF groups had a higher risk of surgical complications (OR = 3.1, p = 0.001). Tumor vascular encasement was associated with increased risk of surgical complications (OR = 5.40, p = 0.009) and increased risk of undifferentiated/poorly differentiated histologic grade (OR = 11.11, p = 0.013). Tumor infiltration of adjacent organs and structures was associated with decreased survival (HR = 8.90, p = 0.007), MYCN amplification (OR = 9.91, p = 0.001), high MKI (OR = 6.20, p = 0.003), and increased risk of International Neuroblastoma Staging System stage 4 disease (OR = 8.96, p < 0.001).

CONCLUSIONS

The presence of IDRFs at diagnosis was associated with high-risk clinical, molecular, and histopathologic features of neuroblastoma. The IDRF group tumor infiltration into adjacent organs and structures was associated with decreased survival. Collectively, these findings may assist surgical planning and medical management for neuroblastoma patients.

Enhancer hijacking determines extrachromosomal circular MYCN amplicon architecture in neuroblastoma

MYCN amplification drives one in six cases of neuroblastoma. The supernumerary gene copies are commonly found on highly rearranged, extrachromosomal circular DNA (ecDNA). The exact amplicon structure has not been described thus far and the functional relevance of its rearrangements is unknown. Here, we analyze the MYCN amplicon structure using short-read and Nanopore sequencing and its chromatin landscape using ChIP-seq, ATAC-seq and Hi-C. This reveals two distinct classes of amplicons which explain the regulatory requirements for MYCN overexpression. The first class always co-amplifies a proximal enhancer driven by the noradrenergic core regulatory circuit (CRC). The second class of MYCN amplicons is characterized by high structural complexity, lacks key local enhancers, and instead contains distal chromosomal fragments harboring CRC-driven enhancers. Thus, ectopic enhancer hijacking can compensate for the loss of local gene regulatory elements and explains a large component of the structural diversity observed in MYCN amplification.

MYCN amplification is common in neuroblastomas. Here the authors analyse the MYCN amplicon structure and its epigenetic regulation by integrating short- and longread genomic and epigenomic data and find two classes of MYCN amplicons in neuroblastomas, one driven by local enhancers and the other by hijacking of distal regulatory elements.

Imbalance between genomic gain and loss identifies high-risk neuroblastoma patients with worse outcomes

Survival in high-risk neuroblastoma (HR-NB) patients remains poor despite multimodal treatment. We aimed to identify HR-NB patients with worse outcomes by analyzing the genomic instability derived from segmental chromosomal aberrations. We calculated 3 genomic instability indexes for primary tumor SNP array profiles from 127 HR-NB patients: (1) Copy number aberration burden (%gains_length+%losses_length), (2) copy number load (CNL) (%gains_length-%losses_length) and (3) net genomic load (NGL) (%gains_amount-%losses_amount). Tumors were classified according to positive or negative CNL and NGL genomic subtypes. The impact of the genomic instability indexes on overall survival (OS) was assessed with Cox regression. We identified 38% of HR-NB patients with poor 5-year OS. A negative CNL genomic background was related to poor prognosis in patients ≥18 months showing tumors with homogeneous MYCN amplification (9.5% survival probability, P < 0.05) and patients with non-MYCN amplified NB (18.8% survival probability related to >2.4% CNL, P < 0.01). A positive CNL genomic background was associated with worse outcome in patients with heterogeneous MYCN amplification (22.5% survival probability, P < 0.05). We conclude that characterizing a tumor genomic background according to predominance of genome gained or lost contributes toward improved outcome prediction and brings greater insight into the tumor biology of HR-NB patients.

Distribution of segmental chromosomal alterations in neuroblastoma

BACKGROUND

Neuroblastoma (NB) is a heterogeneous tumor with extremely diverse prognosis according to clinical and genetic factors such as specific combinations of chromosomal imbalances.

METHODS

Molecular karyotyping data from a national neuroblastic tumor database of 155 NB samples were analyzed and related to clinical data.

RESULTS

Segmental chromosomal alterations (SCA) were detected in 102 NB, whereas 45 only displayed numerical alterations. Incidence of SCA was higher in stage M (92%) and MYCN amplified (MNA) NB (96%). Presence of SCA was associated with older age, especially 1q gain and 3p deletion. 96% of the deaths were observed in the SCA group and 85% of the relapsed NB contained SCA. The alteration most commonly associated with a higher number of other segmental rearrangements was 11q deletion, followed by 4p deletion. Whole-chromosome 19 gain was associated with lower stages, absence of SCA and better outcome.

CONCLUSIONS

SCA are not randomly distributed and are concentrated on recurrent chromosomes. The most frequently affected chromosomes identify prognostic factors in specific risk groups. SCA are associated with older age and MNA. We have identified a small subset of patients with better outcome that share whole-chromosome 19 numeric gain, suggesting its use as a prognostic biomarker in NB.

Analysis and identification of novel biomarkers involved in neuroblastoma via integrated bioinformatics

Neuroblastoma (NB) is the most common extracranial solid tumor in children. Under various treatments, some patients still have a poor prognosis. Hence, it is necessary to find new valid targets for NB therapy. In this study, a comprehensive bioinformatic analysis was used to identify differentially expressed genes (DEGs) between NB and control cells, and to select hub genes associated with NB. GSE66586 and GSE78061 datasets were downloaded from the Gene Expression Omnibus (GEO) database and DEGs were selected. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were applied to the selected DEGs. The STRING database and Cytoscape software were used to construct protein-protein interaction (PPI) networks and perform modular analysis of the DEGs. The R2 database was used for prognostic analysis. We identified a total of 238 DEGs from two microarray databases. GO enrichment analysis shows that these DEGs are mainly concentrated in the regulation of cell growth, cell migration, cell fate determination, and cell maturation. KEGG pathway analysis showed that these DEGs are mainly involved in focal adhesion, the TNF signaling pathway, cancer-related pathways, and signaling pathways regulating stem cell pluripotency. We identified the 15 most closely related DEGs from the PPI network, and performed R2 database prognostic analysis to select five hub genes - CTGF, EDN1, GATA2, LOX, and SERPINE1. This study distinguished hub genes and related signaling pathways that can potentially serve as diagnostic indicators and therapeutic biomarkers for NB, thereby improving understanding of the molecular mechanisms involved in NB.

Kinome multigenic panel identified novel druggable EPHB4-V871I somatic variant in high-risk neuroblastoma

Neuroblastoma (NB) is the most common extracranial neoplasm in children. The overall outcome for high‐risk NB patients is still unacceptable, therefore, it is critical to deeply understand molecular mechanisms associated with NB, which in turn can be utilized for developing drugs towards the treatment of NB. Protein kinases (TKs) play an essential role in the regulation of cell survival and proliferation. Different kinases, such as anaplastic lymphoma kinase (ALK), Aurora kinase, RET receptor tyrosine kinase, are potential therapeutic targets in various cancers, including NB. We analysed a cohort of 45 high‐risk NB patients and 9 NB cell lines by a targeted—(t)NGS custom gene panel (genes codifying for the kinase domains of 90 TKs). We identified somatic variants in four TK genes (ALK, EPHB4, LMTK3 and EPHB6) in NB patients and we functionally characterized an interesting somatic variant, V871I, in EPHB4 gene. EPHB4 plays a crucial role in cardiovascular development and regulates vascularization in cancer‐promoting angiogenesis, tumour growth and metastasis. Several EPHB4 mutations have previously been identified in solid and haematological tumour specimens but EPHB4 mutations were not described until now in NB. Interestingly, a re‐analysis of public CGH‐array showed that the EPHB4 gain is associated with advanced diseases in NB. We further demonstrated that higher EPHB4 expression is correlated to stage 4 of NB and with poor overall survival. Additionally, we also revealed that the EPHB4‐V871I accounts for increased proliferation, migration and invasion properties in two NB cell lines by acting on VEGF, c‐RAF and CDK4 target genes and by increasing the phosphorylation of ERK1‐2 pathway. The use of two EPHB4 inhibitors, JI‐101 and NVP‐BHG712, was able to rescue the phenotype driven by the variant. Our study suggested that EPHB4 is a promising therapeutic target in high‐risk NB.