A mechanistic classification of clinical phenotypes in neuroblastoma

Suppressed miR-128-3p combined with TERT overexpression predicts dismal outcomes for neuroblastoma

BACKGROUND

Molecular and clinical diversity of neuroblastomas is notorious. The activating TERT rearrangements have been associated with dismal prognosis. Suppression of miR-128-3p may complement and enhance the adverse effects of TERT overexpression.

OBJECTIVE

The study aimed at evaluation of prognostic significance of the miR-128-3p/TERT expression in patients with primary neuroblastoma.

METHODS

RNA samples isolated from fresh-frozen tumor specimens (n= 103) were reverse transcribed for evaluation of miR-128-3p and TERT expression by qPCR. The normalized expression levels were tested for correlations with the event-free survival (EFS). ROC-analysis was used to establish threshold expression levels (TLs) for the possible best prediction of the outcomes. The median follow-up was 57 months.

RESULTS

Both TERT overexpression and miR-128-3p downregulation were independently associated with superior rates of adverse events (p= 0.027, TL =-2.32 log10 and p= 0.080, TL =-1.33 log10, respectively). The MYCN single-copy patients were stratified into groups based on the character of alterations in expression of the studied transcripts. Five-year EFS in the groups of patients with elevated TERT/normal miR-128-3p expression and normal TERT/reduced miR-128-3p expression were 0.74 ± 0.08 and 0.60 ± 0.16, respectively. The patients with elevated TERT/reduced miR-128-3p expression had the worst outcomes, with 5-year EFS of 0.40 ± 0.16 compared with 0.91 ± 0.06 for the patients with unaltered levels of both transcripts (p< 0.001). Cumulative incidence of relapse/progression for the groups constituted 0.23 ± 0.08, 0.40 ± 0.16, 0.60 ± 0.16 and 0.09 ± 0.06, respectively. Moreover, the loss of miR-128-3p was qualified as independent adverse predictor which outperformed the conventional clinical and genetic risk factors in the multivariate Cox regression model of EFS.

CONCLUSIONS

Combined expression levels of miR-128-3p and TERT represent a novel prognostic biomarker for neuroblastoma.

Feasibility of whole genome and transcriptome profiling in pediatric and young adult cancers

The utility of cancer whole genome and transcriptome sequencing (cWGTS) in oncology is increasingly recognized. However, implementation of cWGTS is challenged by the need to deliver results within clinically relevant timeframes, concerns about assay sensitivity, reporting and prioritization of findings. In a prospective research study we develop a workflow that reports comprehensive cWGTS results in 9 days. Comparison of cWGTS to diagnostic panel assays demonstrates the potential of cWGTS to capture all clinically reported mutations with comparable sensitivity in a single workflow. Benchmarking identifies a minimum of 80× as optimal depth for clinical WGS sequencing. Integration of germline, somatic DNA and RNA-seq data enable data-driven variant prioritization and reporting, with oncogenic findings reported in 54% more patients than standard of care. These results establish key technical considerations for the implementation of cWGTS as an integrated test in clinical oncology.

Cancer whole-genome and transcriptome sequencing (cWGTS) has been challenging to implement in clinical settings. Here, the authors develop a workflow to deliver robust cWGTS analyses and reports within clinically-relevant timeframes for paediatric, adolescent and young adult solid tumour patients.

Environmental cues from neural crest derivatives act as metastatic triggers in an embryonic neuroblastoma model

Embryonic malignant transformation is concomitant to organogenesis, often affecting multipotent and migratory progenitors. While lineage relationships between malignant cells and their physiological counterparts are extensively investigated, the contribution of exogenous embryonic signals is not fully known. Neuroblastoma (NB) is a childhood malignancy of the peripheral nervous system arising from the embryonic trunk neural crest (NC) and characterized by heterogeneous and interconvertible tumor cell identities. Here, using experimental models mimicking the embryonic context coupled to proteomic and transcriptomic analyses, we show that signals released by embryonic sympathetic ganglia, including Olfactomedin-1, induce NB cells to shift from a noradrenergic to mesenchymal identity, and to activate a gene program promoting NB metastatic onset and dissemination. From this gene program, we extract a core signature specifically shared by metastatic cancers with NC origin. This reveals non-cell autonomous embryonic contributions regulating the plasticity of NB identities and setting pro-dissemination gene programs common to NC-derived cancers.

Diagnosis and Management of Opsoclonus-Myoclonus-Ataxia Syndrome in Children: An International Perspective

BACKGROUND AND OBJECTIVES

Opsoclonus-myoclonus-ataxia syndrome (OMAS) is a rare disorder of the nervous system that classically presents with a combination of characteristic eye movement disorder and myoclonus, in addition to ataxia, irritability, and sleep disturbance. There is good evidence that OMAS is an immune-mediated condition that may be paraneoplastic in the context of neuroblastoma. This syndrome may be associated with long-term cognitive impairment, yet it remains unclear how this is influenced by disease course and treatment. Treatment is largely predicated on immune suppression, but there is limited evidence to indicate an optimal regimen.

METHODS

Following an international multiprofessional workshop in 2004, a body of clinicians and scientists comprising the International OMS Study group continued to meet biennially in a joint professionals and family workshop focusing on pediatric OMAS. Seventeen years after publication of the first report, a writing group was convened to provide a clinical update on the definitions and clinical presentation of OMAS, biomarkers and the role of investigations in a child presenting with OMAS, treatment and management strategies including identification and support of long-term sequelae.

RESULTS

The clinical criteria for diagnosis were reviewed, with a proposed approach to laboratory and radiologic investigation of a child presenting with possible OMAS. The evidence for an upfront vs escalating treatment regimen was reviewed, and a treatment algorithm proposed to recognize both these approaches. Importantly, recommendations on monitoring of immunotherapy response and longer-term follow-up based on an expert consensus are provided.

DISCUSSION

OMAS is a rare neurologic condition that can be associated with poor cognitive outcomes. This report proposes an approach to investigation and treatment of children presenting with OMAS, based on expert international opinion recognizing the limited data available.

The role of ncRNAs in neuroblastoma: mechanisms, biomarkers and therapeutic targets

Neuroblastoma (NB) is a malignant tumor in young children that originates from the neural crest of the sympathetic nervous system. Generally, NB occurs in the adrenal glands, but it can also affect the nerve tissues of the neck, chest, abdomen, and pelvis. Understanding the pathophysiology of NB and developing novel therapeutic approaches are critical. Noncoding RNAs (ncRNAs) are associated with crucial aspects of pathology, metastasis and drug resistance in NB. Here, we summarized the pretranscriptional, transcriptional and posttranscriptional regulatory mechanisms of ncRNAs involved in NB, especially focusing on regulatory pathways. Furthermore, ncRNAs with the potential to serve as biomarkers for risk stratification, drug resistance and therapeutic targets are also discussed, highlighting the clinical application of ncRNAs in NB.

The mutational spectrum of ATRX aberrations in neuroblastoma and the associated patient and tumor characteristics

Neuroblastoma is the most common extracranial solid tumor in children. The chromatin remodeler ATRX is frequently mutated in high‐risk patients with a poor prognosis. Although many studies have reported ATRX aberrations and the associated clinical characteristics in neuroblastoma, a comprehensive overview is currently lacking. In this study, we extensively characterize the mutational spectrum of ATRX aberrations in neuroblastoma tumors reported in previous studies and present an overview of patient and tumor characteristics. We collected the data of a total of 127 neuroblastoma patients and three cell lines with ATRX aberrations originating from 20 papers. We subdivide the ATRX aberrations into nonsense, missense, and multiexon deletions (MEDs) and show that 68% of them are MEDs. Of these MEDs, 75% are predicted to be in‐frame. Furthermore, we identify a missense mutational hotspot region in the helicase domain. We also confirm that all three ATRX mutation types are more often identified in patients diagnosed at an older age, but still approximately 40% of the patients are aged 5 years or younger at diagnosis. Surprisingly, we found that 11q deletions are enriched in neuroblastomas with ATRX deletions compared to a reference cohort, but not in neuroblastomas with ATRX point mutations. Taken together, our data emphasizes a distinct ATRX mutation spectrum in neuroblastoma, which should be considered when studying molecular phenotypes and therapeutic strategies.

EIF4EBP1 is transcriptionally upregulated by MYCN and associates with poor prognosis in neuroblastoma

Neuroblastoma (NB) accounts for 15% of cancer-related deaths in childhood despite considerable therapeutic improvements. While several risk factors, including MYCN amplification and alterations in RAS and p53 pathway genes, have been defined in NB, the clinical outcome is very variable and difficult to predict. Since genes of the mechanistic target of rapamycin (mTOR) pathway are upregulated in MYCN-amplified NB, we aimed to define the predictive value of the mTOR substrate-encoding gene eukaryotic translation initiation factor 4E-binding protein 1 (EIF4EBP1) expression in NB patients. Using publicly available data sets, we found that EIF4EBP1 mRNA expression is positively correlated with MYCN expression and elevated in stage 4 and high-risk NB patients. In addition, high EIF4EBP1 mRNA expression is associated with reduced overall and event-free survival in the entire group of NB patients in three cohorts, as well as in stage 4 and high-risk patients. This was confirmed by monitoring the clinical value of 4EBP1 protein expression, which revealed that high levels of 4EBP1 are significantly associated with prognostically unfavorable NB histology. Finally, functional analyses revealed that EIF4EBP1 expression is transcriptionally controlled by MYCN binding to the EIF4EBP1 promoter in NB cells. Our data highlight that EIF4EBP1 is a direct transcriptional target of MYCN whose high expression is associated with poor prognosis in NB patients. Therefore, EIF4EBP1 may serve to better stratify patients with NB.

Small extracellular vesicles induce resistance to anti-GD2 immunotherapy unveiling tipifarnib as an adjunct to neuroblastoma immunotherapy

Background

Anti-GD2 monoclonal antibody immunotherapy has significantly improved the overall survival rate for high-risk neuroblastoma patients. However, 40% of patients fail to respond or develop resistance to treatment, and the molecular mechanisms by which this occurs remain poorly understood. Tumor-derived small extracellular vesicles (sEVs) have emerged as critical regulators in modulating the response to immunotherapy. In this study, we investigated the role of neuroblastoma-derived sEVs in promoting resistance to the anti-GD2 monoclonal antibody dinutuximab. Moreover, to determine whether pharmacologic inhibition of sEV secretion sensitizes tumors to dinutuximab treatment, we combined dinutuximab with tipifarnib, a farnesyltransferase inhibitor that inhibits sEV secretion.

Methods

We investigated the role of neuroblastoma-derived sEVs in modulating the response to dinutuximab by utilizing the syngeneic 9464D-GD2 mouse model. The effect of neuroblastoma-derived sEVs in modulating the tumor microenvironment (TME) and host immune system were evaluated by RNA-sequencing and flow cytometry. Importantly, we used this mouse model to investigate the efficacy of tipifarnib in sensitizing neuroblastoma tumors to dinutuximab. The effect of tipifarnib on both the TME and host immune system were assessed by flow cytometry.

Results

We demonstrated that neuroblastoma-derived sEVs significantly attenuated the efficacy of dinutuximab in vivo and modulated tumor immune cell infiltration upon dinutuximab treatment to create an immunosuppressive TME that contains more tumor-associated macrophages and fewer tumor-infiltrating NK cells. In addition, we demonstrated that neuroblastoma-derived sEVs suppress splenic NK cell maturation in vivo and dinutuximab-induced NK cell-mediated antibody-dependent cellular cytotoxicity in vitro. Importantly, tipifarnib drastically enhanced the efficacy of dinutuximab-mediated inhibition of tumor growth and prevented the immunosuppressive effects of neuroblastoma-derived sEVs in vivo.

Conclusions

These preclinical findings uncover a novel mechanism by which neuroblastoma-derived sEVs modulate the immune system to promote resistance to dinutuximab and suggest that tipifarnib-mediated inhibition of sEV secretion may serve as a viable treatment strategy to enhance the antitumor efficacy of anti-GD2 immunotherapy in high-risk neuroblastoma patients.

Mechanisms of Senescence-Related NKG2D Ligands Release and Immune Escape Induced by Chemotherapy in Neuroblastoma Cells

Chemotherapy-induced senescence promotes immunocyte aggregation in the tumor microenvironment by upregulating the surface expression of activating ligands in cancer cells. However, these senescent tumor cells cannot be completely cleared and can induce tumor recurrence. Previous studiesshowed that soluble natural killer (NK) group 2D (NKG2D) ligands impair the recognition of multiple immune cells. In this study, we established an in vitro senescence model using neuroblastoma cells subjected to low-dose Chemotherapeutic drug doxorubicin or the Aurora A inhibitor MLN8237. The results showed that different neuroblastoma cell lines showed increased secretion of the NKG2D ligand MHC class I polypeptide-related sequence A/B (MICA/B) following proteolysis after treatment, with MICA/B subsequently recruited to exosomes to downregulate NKG2D expression in NK cells. Interestingly, disintegrin and metalloproteinase domain-containing 10 (ADAM10) was upregulated in senescent tumor cells, and combined treatment with the ADAM10 inhibitor GI254023X and chemotherapeutic drugs inhibited MICA/B secretion and enhanced recognition and killing by NK cells. Additionally, we found that expression of the long noncoding RNA MALAT1 was significantly increased in senescent neuroblastoma cells, and that MALAT1 served as a sponge for microRNA (miR)-92a-3p to counteract miR-92a-3p-mediated repression of ADAM10 levels. Furthermore, administration of a MALAT1 inhibitor or an miR-92a-3p mimic reduced the MICA/B shedding and enhanced recognition and killing by NK cells. These results confirmed that low-dose chemotherapy induces senescence in neuroblastoma cells, and that senescent tumor cells promote the shedding of the NKG2D ligand MICA/B through the MALAT1/miR-92a/ADAM10 axis, thereby contributing to the formation of a suppressive immune microenvironment and promoting immune escape.

Recent advances in the developmental origin of neuroblastoma: an overview

Neuroblastoma (NB) is a pediatric tumor that originates from neural crest-derived cells undergoing a defective differentiation due to genomic and epigenetic impairments. Therefore, NB may arise at any final site reached by migrating neural crest cells (NCCs) and their progeny, preferentially in the adrenal medulla or in the para-spinal ganglia.

NB shows a remarkable genetic heterogeneity including several chromosome/gene alterations and deregulated expression of key oncogenes that drive tumor initiation and promote disease progression.

NB substantially contributes to childhood cancer mortality, with a survival rate of only 40% for high-risk patients suffering chemo-resistant relapse. Hence, NB remains a challenge in pediatric oncology and the need of designing new therapies targeted to specific genetic/epigenetic alterations become imperative to improve the outcome of high-risk NB patients with refractory disease or chemo-resistant relapse.

In this review, we give a broad overview of the latest advances that have unraveled the developmental origin of NB and its complex epigenetic landscape.

Single-cell RNA sequencing with spatial transcriptomics and lineage tracing have identified the NCC progeny involved in normal development and in NB oncogenesis, revealing that adrenal NB cells transcriptionally resemble immature neuroblasts or their closest progenitors. The comparison of adrenal NB cells from patients classified into risk subgroups with normal sympatho-adrenal cells has highlighted that tumor phenotype severity correlates with neuroblast differentiation grade.

Transcriptional profiling of NB tumors has identified two cell identities that represent divergent differentiation states, i.e. undifferentiated mesenchymal (MES) and committed adrenergic (ADRN), able to interconvert by epigenetic reprogramming and to confer intra-tumoral heterogeneity and high plasticity to NB.

Chromatin immunoprecipitation sequencing has disclosed the existence of two super-enhancers and their associated transcription factor networks underlying MES and ADRN identities and controlling NB gene expression programs.

The discovery of NB-specific regulatory circuitries driving oncogenic transformation and maintaining the malignant state opens new perspectives on the design of innovative therapies targeted to the genetic and epigenetic determinants of NB. Remodeling the disrupted regulatory networks from a dysregulated expression, which blocks differentiation and enhances proliferation, toward a controlled expression that prompts the most differentiated state may represent a promising therapeutic strategy for NB.

Somatic mutations enriched in cis-regulatory elements affect genes involved in embryonic development and immune system response in neuroblastoma

Noncoding cis-regulatory variants have gained interest as cancer drivers, yet progress in understanding their significance is hindered by the numerous challenges and limitations of variant prioritization. To overcome these limitations, we focused on active cis-regulatory elements (aCRE) in order to design a customized panel for the deep sequencing of 56 neuroblastoma tumor and normal DNA sample pairs. In order to search for driver mutations, aCREs were defined by reanalysis of H3K27ac ChiP-seq peaks in 25 neuroblastoma cell lines. These regulatory genomic regions were tested for an excess of somatic mutations and assessed for statistical significance using a global approach that accounted for chromatin accessibility and replication timing. Additional validation was provided by whole genome sequence analysis of 151 neuroblastomas. Analysis of Hi-C data determined the presence of candidate target genes interacting with mutated regions. An excess of somatic mutations in aCREs of diverse genes were identified, including IPO7, HAND2, and ARID3A. CRISPR-Cas9 editing was utilized to assess the functional consequences of mutations in the IPO7 aCRE. Patients with noncoding mutations in aCREs showed inferior overall and event-free survival independent of age at diagnosis, stage, risk stratification, and MYCN status. Expression of aCRE-interacting genes correlated strongly with negative prognostic markers and low survival rates. Moreover, a convergence between the biological functions of aCRE target genes and transcription factors with mutated binding motifs was associated with embryonic development and immune system response. Overall, this strategy enabled the identification of somatic mutations in regulatory elements that collectively promote neuroblastoma tumorigenesis.

Dynamic Activity of Histone H3-Specific Chaperone Complexes in Oncogenesis

Canonical histone H3.1 and variant H3.3 deposit at different sites of the chromatin via distinct histone chaperones. Histone H3.1 relies on chaperone CAF-1 to mediate replication-dependent nucleosome assembly during S-phase, while H3.3 variant is regulated and incorporated into the chromatin in a replication-independent manner through HIRA and DAXX/ATRX. Current literature suggests that dysregulated expression of histone chaperones may be implicated in tumor progression. Notably, ectopic expression of CAF-1 can promote a switch between canonical H3.1 and H3 variants in the chromatin, impair the chromatic state, lead to chromosome instability, and impact gene transcription, potentially contributing to carcinogenesis. This review focuses on the chaperone proteins of H3.1 and H3.3, including structure, regulation, as well as their oncogenic and tumor suppressive functions in tumorigenesis.

Urinary 3-Methoxytyramine Is a Biomarker for MYC Activity in Patients With Neuroblastoma

Elevated urinary 3-methoxytyramine (3MT) level at diagnosis was recently put forward as independent risk factor for poor prognosis in neuroblastoma. Here, we investigated the biologic basis underlying the putative association between elevated 3MT levels and poor prognosis.

Elevated urinary 3-methoxytyramine (3MT) at diagnosis is an independent risk factor for poor event-free survival and poor overall survival in neuroblastoma. Our analysis including gene-expression data and urinary 3MT showed that increased urinary 3MT corresponds with tumor MYC activity.

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.