MYCN Amplification Is Associated with Repressed Cellular Immunity in Neuroblastoma: An In Silico Immunological Analysis of TARGET Database

Natural killer cells in neuroblastoma: immunological insights and therapeutic perspectives

Natural killer (NK) cells have multifaceted roles within the complex tumor milieu. They are pivotal components of innate immunity and shape the dynamic landscape of tumor-immune cell interactions, and thus can be leveraged for use in therapeutic interventions. NK-based immunotherapies have had remarkable success in hematological malignancies, but these therapies are met with many challenges in solid tumors, including neuroblastoma (NB), a childhood tumor arising from the sympathetic nervous system. With a focus on NB, this review outlines the mechanisms employed by NK cells to recognize and eliminate malignant cells, delving into the dynamic relationship between ligand-receptor interactions, cytokines, and other molecules that facilitate the cross talk between NK and NB cells. We discuss the immunomodulatory functions of NK cells and the mechanisms that contribute to loss of this immunosurveillance in NB, with a focus on how this dynamic has been utilized in recent immunotherapy advancements for NB.

Survival Distinctions for Cases Representing Immunologically Cold Tumors via Intrinsic Disorder Assessments for Blood-Sourced TRB Variable Regions

T cell receptor beta (TRB) sequences were recovered from the Cancer Genome Atlas Uveal Melanoma blood exome files. Intrinsic disorder scores for amino acid (AA) sequences of the entire TRB variable region were obtained and evaluated as potentially representative of overall survival (OS) distinctions, i.e., for cases representing the upper and lower 50th percentiles for intrinsic disorder scores. Analyses using four intrinsic disorder assessment tools indicated that a lower intrinsic disorder of the blood-sourced TRB variable regions, including continuous AA sequences of the V-gene segment, the complementarity-determining region-3, and the J-gene segment, was associated with a better OS probability (with log-rank p-values ranging from 0.002 to 0.014). We further determined that intrinsic disorder assessments could be used for OS stratification for a second, immunologically cold cancer: MYCN amplified neuroblastoma. Thus, intrinsic disorder assessments of blood-sourced, full TRB variable regions may provide a novel patient stratification approach for patients with immunologically cold cancers.

Immune evasion: An imperative and consequence of MYC deregulation

MYC has been implicated in the pathogenesis of a wide range of human tumors and has been described for many years as a transcription factor that regulates genes with pleiotropic functions to promote tumorigenic growth. However, despite extensive efforts to identify specific target genes of MYC that alone could be responsible for promoting tumorigenesis, the field is yet to reach a consensus whether this is the crucial function of MYC. Recent work shifts the view on MYC's function from being a gene-specific transcription factor to an essential stress resilience factor. In highly proliferating cells, MYC preserves cell integrity by promoting DNA repair at core promoters, protecting stalled replication forks, and/or preventing transcription-replication conflicts. Furthermore, an increasing body of evidence demonstrates that MYC not only promotes tumorigenesis by driving cell-autonomous growth, but also enables tumors to evade the host's immune system. In this review, we summarize our current understanding of how MYC impairs antitumor immunity and why this function is evolutionarily hard-wired to the biology of the MYC protein family. We show why the cell-autonomous and immune evasive functions of MYC are mutually dependent and discuss ways to target MYC proteins in cancer therapy.

Current Knowledge and Perspectives of Immunotherapies for Neuroblastoma

Neuroblastoma (NBL) cells highly express disialoganglioside GD2, which is restricted and weakly expressed in selected healthy cells, making it a desirable target of immunotherapy. Over the past two decades, application of dinutuximab, an anti-GD2 monoclonal antibody (mAb), has been one of the few new therapies to substantially improve outcomes to current levels. Given the persistent challenge of relapse and therapeutic resistance, there is an urgent need for new effective and tolerable treatment options for high-risk NBL. Recent breakthroughs in immune checkpoint inhibitor (ICI) therapeutics have not translated into high-risk NBL, like many other major pediatric solid tumors. Given the suppressed tumor microenvironment (TME), single ICIs like anti-CTLA4 and anti-PD1 have not demonstrated significant antitumor response rates. Meanwhile, emerging studies are reporting novel advancements in GD2-based therapies, targeted therapies, nanomedicines, and other immunotherapies such as adoptive transfer of natural killer (NK) cells and chimeric antigen receptors (CARs), and these hold interesting promise for the future of high-risk NBL patient care. Herein, we summarize the current state of the art in NBL therapeutic options and highlight the unique challenges posed by NBL that have limited the successful adoption of immune-modifying therapies. Through this review, we aim to direct the field's attention to opportunities that may benefit from a combination immunotherapy strategy.

N-MYC impairs innate immune signaling in high-grade serous ovarian carcinoma

High-grade serous ovarian cancer (HGSC) is a challenging disease, especially for patients with immunologically "cold" tumors devoid of tumor-infiltrating lymphocytes (TILs). We found that HGSC exhibits among the highest levels of MYCN expression and transcriptional signature across human cancers, which is strongly linked to diminished features of antitumor immunity. N-MYC repressed basal and induced IFN type I signaling in HGSC cell lines, leading to decreased chemokine expression and T cell chemoattraction. N-MYC inhibited the induction of IFN type I by suppressing tumor cell-intrinsic STING signaling via reduced STING oligomerization, and by blunting RIG-I-like receptor signaling through inhibition of MAVS aggregation and localization in the mitochondria. Single-cell RNA sequencing of human clinical HGSC samples revealed a strong negative association between cancer cell-intrinsic MYCN transcriptional program and type I IFN signaling. Thus, N-MYC inhibits tumor cell-intrinsic innate immune signaling in HGSC, making it a compelling target for immunotherapy of cold tumors.

CKLF instigates a "cold" microenvironment to promote MYCN-mediated tumor aggressiveness

Solid tumors, especially those with aberrant MYCN activation, often harbor an immunosuppressive microenvironment to fuel malignant growth and trigger treatment resistance. Despite this knowledge, there are no effective strategies to tackle this problem. We found that chemokine-like factor (CKLF) is highly expressed by various solid tumor cells and transcriptionally up-regulated by MYCN. Using the MYCN-driven high-risk neuroblastoma as a model system, we demonstrated that as early as the premalignant stage, tumor cells secrete CKLF to attract CCR4-expressing CD4+ cells, inducing immunosuppression and tumor aggression. Genetic depletion of CD4+ T regulatory cells abolishes the immunorestrictive and protumorigenic effects of CKLF. Our work supports that disrupting CKLF-mediated cross-talk between tumor and CD4+ suppressor cells represents a promising immunotherapeutic approach to battling MYCN-driven tumors.

HDAC6-MYCN-CXCL3 axis mediates allergic inflammation and is necessary for allergic inflammation-promoted cellular interactions

Histone deacetylase 6 (HDAC6) has been shown to play an important role in allergic inflammation. This study hypothesized that novel downstream targets of HDAC6 would mediate allergic inflammation. Experiments employing HDAC6 knock out C57BL/6 mice showed that HDAC6 mediated passive cutaneous anaphylaxis (PCA) and passive systemic anaphylaxis (PSA). Antigen stimulation increased expression of N-myc (MYCN) and CXCL3 in an HDAC6-dependent manner in the bone marrow-derived mast cells. MYCN and CXCL3 were necessary for both PCA and PSA. The role of early growth response 3 (EGR3) in the regulation of HDAC6 expression has been reported. ChIP assays showed EGR3 as a direct regulator of MYCN. miR-34a-5p was predicted to be a negative regulator of MYCN. Luciferase activity assays showed miR-34a-5p as a direct regulator of MYCN. miR-34a-5p mimic negatively regulated PCA and PSA. MYCN decreased miR-34a-5p expression in antigen-stimulated rat basophilic leukemia cells (RBL2H3). MYCN was shown to bind to the promoter sequence of CXCL3. In an IgE-independent manner, recombinant CXCL3 protein increased expression of HDAC6, MYCN, and β-hexosaminidase activity in RBL2H3 cells. Mouse recombinant CXCL3 protein enhanced the angiogenic potential of the culture medium of RBL2H3. CXCL3 was necessary for the enhanced angiogenic potential of the culture medium of antigen-stimulated RBL2H3. The culture medium of RBL2H3 was able to induce M2 macrophage polarization in a CXCL3-dependent manner. Recombinant CXCL3 protein also increased the expression of markers of M2 macrophage. Thus, the identification of the novel role of HDAC6-MYCN-CXCL3 axis can help better understand the pathogenesis of anaphylaxis.

Divergent HLA variations and heterogeneous expression but recurrent HLA loss-of- heterozygosity and common HLA-B and TAP transcriptional silencing across advanced pediatric solid cancers

The human leukocyte antigen (HLA) system is a major factor controlling cancer immunosurveillance and response to immunotherapy, yet its status in pediatric cancers remains fragmentary. We determined high-confidence HLA genotypes in 576 children, adolescents and young adults with recurrent/refractory solid tumors from the MOSCATO-01 and MAPPYACTS trials, using normal and tumor whole exome and RNA sequencing data and benchmarked algorithms. There was no evidence for narrowed HLA allelic diversity but discordant homozygosity and allele frequencies across tumor types and subtypes, such as in embryonal and alveolar rhabdomyosarcoma, neuroblastoma MYCN and 11q subtypes, and high-grade glioma, and several alleles may represent protective or susceptibility factors to specific pediatric solid cancers. There was a paucity of somatic mutations in HLA and antigen processing and presentation (APP) genes in most tumors, except in cases with mismatch repair deficiency or genetic instability. The prevalence of loss-of-heterozygosity (LOH) ranged from 5.9 to 7.7% in HLA class I and 8.0 to 16.7% in HLA class II genes, but was widely increased in osteosarcoma and glioblastoma (~15-25%), and for DRB1-DQA1-DQB1 in Ewing sarcoma (~23-28%) and low-grade glioma (~33-50%). HLA class I and HLA-DR antigen expression was assessed in 194 tumors and 44 patient-derived xenografts (PDXs) by immunochemistry, and class I and APP transcript levels quantified in PDXs by RT-qPCR. We confirmed that HLA class I antigen expression is heterogeneous in advanced pediatric solid tumors, with class I loss commonly associated with the transcriptional downregulation of HLA-B and transporter associated with antigen processing (TAP) genes, whereas class II antigen expression is scarce on tumor cells and occurs on immune infiltrating cells. Patients with tumors expressing sufficient HLA class I and TAP levels such as some glioma, osteosarcoma, Ewing sarcoma and non-rhabdomyosarcoma soft-tissue sarcoma cases may more likely benefit from T cell-based approaches, whereas strategies to upregulate HLA expression, to expand the immunopeptidome, and to target TAP-independent epitopes or possibly LOH might provide novel therapeutic opportunities in others. The consequences of HLA class II expression by immune cells remain to be established. Immunogenetic profiling should be implemented in routine to inform immunotherapy trials for precision medicine of pediatric cancers.

Targeting the myeloid microenvironment in neuroblastoma

Myeloid cells (granulocytes and monocytes/macrophages) play an important role in neuroblastoma. By inducing a complex immunosuppressive network, myeloid cells pose a challenge for the adaptive immune system to eliminate tumor cells, especially in high-risk neuroblastoma. This review first summarizes the pro- and anti-tumorigenic functions of myeloid cells, including granulocytes, monocytes, macrophages, and myeloid-derived suppressor cells (MDSC) during the development and progression of neuroblastoma. Secondly, we discuss how myeloid cells are engaged in the current treatment regimen and explore novel strategies to target these cells in neuroblastoma. These strategies include: (1) engaging myeloid cells as effector cells, (2) ablating myeloid cells or blocking the recruitment of myeloid cells to the tumor microenvironment and (3) reprogramming myeloid cells. Here we describe that despite their immunosuppressive traits, tumor-associated myeloid cells can still be engaged as effector cells, which is clear in anti-GD2 immunotherapy. However, their full potential is not yet reached, and myeloid cell engagement can be enhanced, for example by targeting the CD47/SIRPα axis. Though depletion of myeloid cells or blocking myeloid cell infiltration has been proven effective, this strategy also depletes possible effector cells for immunotherapy from the tumor microenvironment. Therefore, reprogramming of suppressive myeloid cells might be the optimal strategy, which reverses immunosuppressive traits, preserves myeloid cells as effectors of immunotherapy, and subsequently reactivates tumor-infiltrating T cells.

Correlations between contrast-enhanced CT-measured extracellular volume fraction, histopathological features, and MYCN amplification status in abdominal neuroblastoma: a retrospective study

PURPOSE

To retrospectively investigate the correlations between contrast-enhanced CT (CECT)-measured extracellular volume fraction (fECV) and histopathological features, as well as MYCN amplification status, in abdominal neuroblastoma.

MATERIALS AND METHODS

One hundred and forty-one patients with abdominal neuroblastoma who underwent CECT scanning were retrospectively enrolled. Calculation of fECV involved the measurement of CT values within regions of interest located within the neuroblastoma and aorta on both non-contrast-enhanced CT and equilibrium CECT. The correlations between fECV and various factors, including pathological subtype, mitosis karyorrhexis index (MKI), Shimada classification, MYCN amplification status, International Neuroblastoma Risk Group (INRG) stage, and risk group were analyzed using either the Mann-Whitney U test or Kruskal-Wallis test.

RESULTS

Neuroblastoma and ganglioneuroblastoma exhibited fECV values of 0.349 (0.252, 0.424) and 0.438 (0.327, 0.508), respectively, indicating a statistically significant difference (Z = 2.200, P = 0.028). Additionally, the fECV decreased significantly in neuroblastoma with high MKI (H = 8.314, P = 0.016) or unfavorable histology (Z = 3.880, P < 0.001), as well as in those with MYCN amplification (Z = 5.486, P < 0.001). Notably, a significant variation in fECV was observed among different INRG stages (H = 16.881, P <0.001) and risk groups (H = 29.014, P < 0.001).

CONCLUSION

CECT-derived fECV is associated with histopathological features, MYCN amplification status, INRG stage, and risk stratification of abdominal neuroblastoma, reflecting a potential correlation between the extracellular matrix and the biological behavior of neuroblastoma.

Neuroblastoma and Glioblastoma Cases With Amplified Oncogenes Have Reduced Numbers of Tumor-Resident Adaptive Immune Receptor Recombinations

PURPOSE

In certain cancers, oncogene amplification is correlated with an immunologically cold or noninflamed, tumor immune microenvironment (TIME) and a worse prognosis, for example, in the case of MYCN-amplified neuroblastoma (NBL). However, for other cancer types, the relationship between oncogene amplification and immune response is more complicated or unresolved. One such cancer is glioblastoma multiforme (GBM), in which the epidermal growth factor receptor (EGFR) oncogene is commonly amplified. Unlike MYCN-amplified NBL, EGFR-amplified GBM has not been shown to correlate with a distinct survival probability.

METHODS

Given this contrasting state for NBL and GBM, we sought to apply a genomics approach to evaluating the immune response for cases with gene amplification.

RESULTS

Our results confirmed and added further specificity to the cold TIME of MYCN-amplified NBL. Moreover, we demonstrated a novel state of immunologically cold EGFR-amplified GBM tumors.

CONCLUSION

This approach to using copy number variation and immune receptor recombination read recovery levels to assess gene amplification and TIME, respectively, may be particularly efficient for the rapid evaluation of many other cancer types.

Single-cell transcriptomics and epigenomics unravel the role of monocytes in neuroblastoma bone marrow metastasis

Metastasis is the major cause of cancer-related deaths. Neuroblastoma (NB), a childhood tumor has been molecularly defined at the primary cancer site, however, the bone marrow (BM) as the metastatic niche of NB is poorly characterized. Here we perform single-cell transcriptomic and epigenomic profiling of BM aspirates from 11 subjects spanning three major NB subtypes and compare these to five age-matched and metastasis-free BM, followed by in-depth single cell analyses of tissue diversity and cell-cell interactions, as well as functional validation. We show that cellular plasticity of NB tumor cells is conserved upon metastasis and tumor cell type composition is NB subtype-dependent. NB cells signal to the BM microenvironment, rewiring via macrophage mgration inhibitory factor and midkine signaling specifically monocytes, which exhibit M1 and M2 features, are marked by activation of pro- and anti-inflammatory programs, and express tumor-promoting factors, reminiscent of tumor-associated macrophages. The interactions and pathways characterized in our study provide the basis for therapeutic approaches that target tumor-to-microenvironment interactions.

Biological Insight and Recent Advancement in the Treatment of Neuroblastoma

One of the most frequent solid tumors in children is neuroblastoma, which has a variety of clinical behaviors that are mostly influenced by the biology of the tumor. Unique characteristics of neuroblastoma includes its early age of onset, its propensity for spontaneous tumor regression in newborns, and its high prevalence of metastatic disease at diagnosis in individuals older than 1 year of age. Immunotherapeutic techniques have been added to the previously enlisted chemotherapeutic treatments as therapeutic choices. A groundbreaking new treatment for hematological malignancies is adoptive cell therapy, specifically chimeric antigen receptor (CAR) T cell therapy. However, due to the immunosuppressive nature of the tumor microenvironment (TME) of neuroblastoma tumor, this treatment approach faces difficulties. Numerous tumor-associated genes and antigens, including the MYCN proto-oncogene (MYCN) and disialoganglioside (GD2) surface antigen, have been found by the molecular analysis of neuroblastoma cells. The MYCN gene and GD2 are two of the most useful immunotherapy findings for neuroblastoma. The tumor cells devise numerous methods to evade immune identification or modify the activity of immune cells. In addition to addressing the difficulties and potential advancements of immunotherapies for neuroblastoma, this review attempts to identify important immunological actors and biological pathways involved in the dynamic interaction between the TME and immune system.

Macrophage infiltration promotes regrowth in MYCN-amplified neuroblastoma after chemotherapy

Despite aggressive treatment, the 5-year event-free survival rate for children with high-risk neuroblastoma is <50%. While most high-risk neuroblastoma patients initially respond to treatment, often with complete clinical remission, many eventually relapse with therapy-resistant tumors. Novel therapeutic alternatives that prevent the recurrence of therapy-resistant tumors are urgently needed. To understand the adaptation of neuroblastoma under therapy, we analyzed the transcriptomic landscape in 46 clinical tumor samples collected before (PRE) or after (POST) treatment from 22 neuroblastoma patients. RNA sequencing revealed that many of the top-upregulated biological processes in POST MYCN amplified (MNA⁺) tumors compared to PRE MNA⁺ tumors were immune-related, and there was a significant increase in numerous genes associated with macrophages. The infiltration of macrophages was corroborated by immunohistochemistry and spatial digital protein profiling. Moreover, POST MNA⁺ tumor cells were more immunogenic compared to PRE MNA⁺ tumor cells. To find support for the macrophage-induced outgrowth of certain subpopulations of immunogenic tumor cells following treatment, we examined the genetic landscape in multiple clinical PRE and POST tumor samples from nine neuroblastoma patients revealing a significant correlation between an increased amount of copy number aberrations (CNA) and macrophage infiltration in POST MNA⁺ tumor samples. Using an in vivo neuroblastoma patient-derived xenograft (PDX) chemotherapy model, we further show that inhibition of macrophage recruitment with anti-CSF1R treatment prevents the regrowth of MNA⁺ tumors following chemotherapy. Taken together, our work supports a therapeutic strategy for fighting the relapse of MNA⁺ neuroblastoma by targeting the immune microenvironment.

The immune landscape of solid pediatric tumors

Background

Large immunogenomic analyses have demonstrated the prognostic role of the functional orientation of the tumor microenvironment in adult solid tumors, this variable has been poorly explored in the pediatric counterpart.

Methods

We performed a systematic analysis of public RNAseq data (TARGET) for five pediatric tumor types (408 patients): Wilms tumor (WLM), neuroblastoma (NBL), osteosarcoma (OS), clear cell sarcoma of the kidney (CCSK) and rhabdoid tumor of the kidney (RT). We assessed the performance of the Immunologic Constant of Rejection (ICR), which captures an active Th1/cytotoxic response. We also performed gene set enrichment analysis (ssGSEA) and clustered more than 100 well characterized immune traits to define immune subtypes and compared their outcome.

Results

A higher ICR score was associated with better survival in OS and high risk NBL without MYCN amplification but with poorer survival in WLM. Clustering of immune traits revealed the same five principal modules previously described in adult tumors (TCGA). These modules divided pediatric patients into six immune subtypes (S1-S6) with distinct survival outcomes. The S2 cluster showed the best overall survival, characterized by low enrichment of the wound healing signature, high Th1, and low Th2 infiltration, while the reverse was observed in S4. Upregulation of the WNT/Beta-catenin pathway was associated with unfavorable outcomes and decreased T-cell infiltration in OS.

Conclusions

We demonstrated that extracranial pediatric tumors could be classified according to their immune disposition, unveiling similarities with adults’ tumors. Immunological parameters might be explored to refine diagnostic and prognostic biomarkers and to identify potential immune-responsive tumors.

Opsoclonus-myoclonus syndrome associated with neuroblastoma: Insights into antitumor immunity

Opsoclonus-myoclonus syndrome (OMS) is a rare neurological disorder. Half of these cases occur in children with neuroblastoma. Neuroblastoma patients with OMS usually have better oncological outcomes than those without OMS even after stratification by tumor stage and age, indicating that factors mediating OMS may also inhibit tumor cell proliferation. Although the mechanisms underlying OMS remain undefined, the cytokines and lymphocytes alterations in the cerebrospinal fluid support the concept that it is a pattern of neuroinflammation due to an autoimmune effect. The presence of lymphoid follicles consisting of follicular dendritic cells, CD20⁺ B lymphocytes, CD3⁺ T lymphocytes, and CD68⁺ macrophages in the tumor microenvironment in OMS-associated neuroblastoma support the autoimmune nature of this disorder. This review focuses on the clinical and genetic features of OMS-associated neuroblastoma, and we update readers on immune features of neuroblastoma with or without OMS to gain insights into antitumor immunity as it relates to tumor biology and prognosis.

Engineered tumor cell-derived vaccines against cancer: The art of combating poison with poison

Tumor vaccination is a promising approach for tumor immunotherapy because it presents high specificity and few side effects. However, tumor vaccines that contain only a single tumor antigen can allow immune system evasion by tumor variants. Tumor antigens are complex and heterogeneous, and identifying a single antigen that is uniformly expressed by tumor cells is challenging. Whole tumor cells can produce comprehensive antigens that trigger extensive tumor-specific immune responses. Therefore, tumor cells are an ideal source of antigens for tumor vaccines. A better understanding of tumor cell-derived vaccines and their characteristics, along with the development of new technologies for antigen delivery, can help improve vaccine design. In this review, we summarize the recent advances in tumor cell-derived vaccines in cancer immunotherapy and highlight the different types of engineered approaches, mechanisms, administration methods, and future perspectives. We discuss tumor cell-derived vaccines, including whole tumor cell components, extracellular vesicles, and cell membrane-encapsulated nanoparticles. Tumor cell-derived vaccines contain multiple tumor antigens and can induce extensive and potent tumor immune responses. However, they should be engineered to overcome limitations such as insufficient immunogenicity and weak targeting. The genetic and chemical engineering of tumor cell-derived vaccines can greatly enhance their targeting, intelligence, and functionality, thereby realizing stronger tumor immunotherapy effects. Further advances in materials science, biomedicine, and oncology can facilitate the clinical translation of tumor cell-derived vaccines.

Impaired Antitumor Immune Response in MYCN-amplified Neuroblastoma Is Associated with Lack of CCL2 Secretion and Poor Dendritic Cell Recruitment

In neuroblastoma, MYCN amplification is associated with sparse immune infiltrate and poor prognosis. Dendritic cells (DC) are crucial immune sentinels but their involvement in neuroblastoma pathogenesis is poorly understood. We observed that the migration of monocytes, myeloid and plasmacytoid DC induced by MYCN-nonamplified neuroblastoma supernatants was abrogated by the addition of anti-CCL2 antibodies, demonstrating the involvement of the CCR2/CCL2 axis in their recruitment by these tumors. Using public RNA sequencing and microarray datasets, we describe lower level of expression of CCL2 in MYCN-amplified neuroblastoma tumors, and we propose a working model for T-cell recruitment in neuroblastoma tumors in which CCL2 produced by neuroblastoma cells initiates the recruitment of monocytes, myeloid and plasmacytoid DCs. Among these cells, the CD1c⁺ subset may recruit T cells by means of CCL19/CCL22 secretion. In vitro, supernatants from DCs cocultured with neuroblastoma cell lines and activated contain CCL22 and CCL19, and are chemotactic for both CD4⁺ and CD8⁺ T cells. We also looked at immunomodulation induced by neuroblastoma cell lines, and found MYCN-nonamplified neuroblastoma cell lines were able to create a microenvironment where DC activation is enhanced. Overall, our findings highlight a major role for CCL2/CCR2 axis in monocytes, myeloid and plasmacytoid cells recruitment toward MYCN-nonamplified neuroblastoma, allowing further immune cell recruitment, and show that these tumors present a microenvironment that can favor DC responses.

Significance

In MYCN-nonamplified neuroblastoma, CCL2 produced by neuroblastoma cells induces the recruitment of antigen-presenting cells (DCs and monocytes/macrophages), allowing infiltration by T cells, in link with CCL19 and CCL22 production, hence favoring immune responses.

The immune cell atlas of human neuroblastoma

Understanding the complete immune cell composition of human neuroblastoma (NB) is crucial for the development of immunotherapeutics. Here, we perform single-cell RNA sequencing (scRNA-seq) on 19 human NB samples coupled with multiplex immunohistochemistry, survival analysis, and comparison with normal fetal adrenal gland data. We provide a comprehensive immune cell landscape and characterize cell-state changes from normal tissue to NB. Our analysis reveals 27 immune cell subtypes, including distinct subpopulations of myeloid, NK, B, and T cells. Several different cell types demonstrate a survival benefit. In contrast to adult cancers and previous NB studies, we show an increase in inflammatory monocyte cell state when contrasting normal and tumor tissue, while no differences in cytotoxicity and exhaustion score for T cells, nor in Treg activity, are observed. Our receptor-ligand interaction analysis reveals a highly complex interactive network of the NB microenvironment from which we highlight several interactions that we suggest for future therapeutic studies.

Macrophage-mediated anti-tumor immunity against high-risk neuroblastoma

Neuroblastoma is the most common extracranial childhood solid tumor. The majority of high-risk neuroblastoma is resistant/refractory to the current high intensity therapy. Neuroblastoma lacks classical HLA Class I expression and exhibits low mutation burden, allowing neuroblastoma cells to evade CD8+ T cell-mediated immunity. Neuroblastoma cells do not express PD-L1, and tumor-associated macrophages are the predominant PD-L1+ cells in the tumor. In this study, we performed gene expression profiling and survival analyses on large neuroblastoma datasets to address the prognostic effect of PD-L1 gene expression and the possible involvement of the SLAMF7 pathway in the anti-neuroblastoma immunity. High-level expression of PD-L1 was found significantly associated with better outcome of high-risk neuroblastoma patients; two populations of PD-1+ PD-L1+ macrophages could be present in high-risk tumors with PD-1/PD-L1 ratios, ≈1 and >1. Patients with the PD-1/PD-L1 ratio >1 tumor showed inferior survival. High-level co-expression of SLAMF7 and SH2D1B was significantly associated with better survival of the high-risk neuroblastoma patients. Together, this study supports the hypothesis that macrophages are important effector cells in the anti-high-risk neuroblastoma immunity, that PD-1 blockade therapy can be beneficial to the high-risk neuroblastoma subset with the PD-1/PD-L1 expression ratio >1, and that SLAMF7 is a new therapeutic target of high-risk neuroblastoma.

Neuroblastoma Heterogeneity, Plasticity, and Emerging Therapies

Purpose of Review

The evolving information of the initiation, tumor cell heterogeneity, and plasticity of childhood neuroblastoma has opened up new perspectives for developing therapies based on detailed knowledge of the disease.

Recent Findings

The cellular origin of neuroblastoma has begun to unravel and there have been several reports on tumor cell heterogeneity based on transcriptional core regulatory circuitries that have given us important information on the biology of neuroblastoma as a developmental disease. This together with new insight of the tumor microenvironment which acts as a support for neuroblastoma growth has given us the prospect for designing better treatment approaches for patients with high-risk neuroblastoma. Here, we discuss these new discoveries and highlight some emerging therapeutic options.

Summary

Neuroblastoma is a disease with multiple facets. Detailed biological and molecular knowledge on neuroblastoma initiation, heterogeneity, and the communications between cells in the tumor microenvironment holds promise for better therapies.

The Thermal Dose of Photothermal Therapy Generates Differential Immunogenicity in Human Neuroblastoma Cells

Photothermal therapy (PTT) is an effective method for tumor eradication and has been successfully combined with immunotherapy. However, besides its cytotoxic effects, little is known about the effect of the PTT thermal dose on the immunogenicity of treated tumor cells. Therefore, we administered a range of thermal doses using Prussian blue nanoparticle-based photothermal therapy (PBNP-PTT) and assessed their effects on tumor cell death and concomitant immunogenicity correlates in two human neuroblastoma cell lines: SH-SY5Y (MYCN-non-amplified) and LAN-1 (MYCN-amplified). PBNP-PTT generated thermal dose-dependent tumor cell killing and immunogenic cell death (ICD) in both tumor lines in vitro. However, the effect of the thermal dose on ICD and the expression of costimulatory molecules, immune checkpoint molecules, major histocompatibility complexes, an NK cell-activating ligand, and a neuroblastoma-associated antigen were significantly more pronounced in SH-SY5Y cells compared with LAN-1 cells, consistent with the high-risk phenotype of LAN-1 cells. In functional co-culture studies in vitro, T cells exhibited significantly higher cytotoxicity toward SH-SY5Y cells relative to LAN-1 cells at equivalent thermal doses. This preliminary report suggests the importance of moving past the traditional focus of using PTT solely for tumor eradication to one that considers the immunogenic effects of PTT thermal dose to facilitate its success in cancer immunotherapy.

Epigenetic state determines inflammatory sensing in neuroblastoma

Immunotherapy has revolutionized cancer treatment, but many cancers are not impacted by currently available immunotherapeutic strategies. Here, we investigated inflammatory signaling pathways in neuroblastoma, a classically "cold" pediatric cancer. By testing the functional response of a panel of 20 diverse neuroblastoma cell lines to three different inflammatory stimuli, we found that all cell lines have intact interferon signaling, and all but one lack functional cytosolic DNA sensing via cGAS-STING. However, double-stranded RNA (dsRNA) sensing via Toll-like receptor 3 (TLR3) was heterogeneous, as was signaling through other dsRNA sensors and TLRs more broadly. Seven cell lines showed robust response to dsRNA, six of which are in the mesenchymal epigenetic state, while all unresponsive cell lines are in the adrenergic state. Genetically switching adrenergic cell lines toward the mesenchymal state fully restored responsiveness. In responsive cells, dsRNA sensing results in the secretion of proinflammatory cytokines, enrichment of inflammatory transcriptomic signatures, and increased tumor killing by T cells in vitro. Using single-cell RNA sequencing data, we show that human neuroblastoma cells with stronger mesenchymal signatures have a higher basal inflammatory state, demonstrating intratumoral heterogeneity in inflammatory signaling that has significant implications for immunotherapeutic strategies in this aggressive childhood cancer.

Sensitizing tumors to anti-PD-1 therapy by promoting NK and CD8+ T cells via pharmacological activation of FOXO3

BACKGROUND

Stimulating antitumor immunity by blocking programmed death-1 (PD-1) or its ligand (programmed death-ligand 1 (PD-L1) is a promising antitumor therapy. However, numerous patients respond poorly to PD-1/PD-L1 blockade. Unresponsiveness to immune-checkpoint blockade (ICB) can cast significant challenges to the therapeutic options for patients with hard-to-treat tumors. There is an unmet clinical need to establish new therapeutic approaches for mitigating ICB unresponsiveness in patients. In this study, we investigated the efficacy and role of low-dose antineoplastic agent SN-38 or metformin in sensitizing unresponsive tumors to respond to ICB therapy.

METHODS

We assessed the significant pathological relationships between PD-L1 and FOXO3 expression and between PD-L1 and c-Myc or STAT3 expression in patients with various tumors. We determined the efficacy of low-dose SN-38 or metformin in sensitizing unresponsive tumors to respond to anti-PD-1 therapy in a syngeneic tumor system. We deciphered novel therapeutic mechanisms underlying the SN-38 and anti-PD-1 therapy-mediated engagement of natural killer (NK) or CD8+ T cells to infiltrate tumors and boost antitumor immunity.

RESULTS

We showed that PD-L1 protein level was inversely associated with FOXO3 protein level in patients with ovarian, breast, and hepatocellular tumors. Low-dose SN-38 or metformin abrogated PD-L1 protein expression, promoted FOXO3 protein level, and significantly increased the animal survival rate in syngeneic mouse tumor models. SN-38 or metformin sensitized unresponsive tumors responding to anti-PD-1 therapy by engaging NK or CD8+ T cells to infiltrate the tumor microenvironment (TME) and secret interferon-γ and granzyme B to kill tumors. SN-38 suppressed the levels of c-Myc and STAT3 proteins, which controlled PD-L1 expression. FOXO3 was essential for SN38-mediated PD-L1 suppression. The expression of PD-L1 was compellingly linked to that of c-Myc or STAT3 in patients with the indicated tumors.

CONCLUSION

We show that SN-38 or metformin can boost antitumor immunity in the TME by inhibiting c-Myc and STAT3 through FOXO3 activation. These results may provide novel insight into ameliorating patient response to overarching immunotherapy for tumors.

MYC regulates metabolism through vesicular transfer of glycolytic kinases

Amplification of the proto-oncogene MYCN is a key molecular aberration in high-risk neuroblastoma and predictive of poor outcome in this childhood malignancy. We investigated the role of MYCN in regulating the protein cargo of extracellular vesicles (EVs) secreted by tumour cells that can be internalized by recipient cells with functional consequences. Using a switchable MYCN system coupled to mass spectrometry analysis, we found that MYCN regulates distinct sets of proteins in the EVs secreted by neuroblastoma cells. EVs produced by MYCN-expressing cells or isolated from neuroblastoma patients induced the Warburg effect, proliferation and c-MYC expression in target cells. Mechanistically, we linked the cancer-promoting activity of EVs to the glycolytic kinase pyruvate kinase M2 (PKM2) that was enriched in EVs secreted by MYC-expressing neuroblastoma cells. Importantly, the glycolytic enzymes PKM2 and hexokinase II were detected in the EVs circulating in the bloodstream of neuroblastoma patients, but not in those of non-cancer children. We conclude that MYC-activated cancers might spread oncogenic signals to remote body locations through EVs.

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.

An Engineered Prussian Blue Nanoparticles-based Nanoimmunotherapy Elicits Robust and Persistent Immunological Memory in a TH-MYCN Neuroblastoma Model

A combination therapy using Prussian blue nanoparticles (PBNP) as photothermal therapy (PTT) agents coated with CpG oligodeoxynucleotides, an immunologic adjuvant, as a nanoimmunotherapy (CpG-PBNP-PTT) for neuroblastoma (NB) is described. NB driven by MYCN amplification confers high risk and correlates with a dismal prognosis, accounting for the majority of NB-related mortality. The efficacy of the CpG-PBNP-PTT nanoimmunotherapy in a clinically relevant, TH-MYCN murine NB model (9464D) overexpressing MYCN is tested. When administered to 9464D NB cells in vitro, CpG-PBNP-PTT triggers thermal dose-dependent immunogenic cell death and tumor cell priming for immune recognition in vitro, measured by the expression of specific costimulatory and antigen-presenting molecules. In vivo, intratumorally administered CpG-PBNP-PTT generates complete tumor regression and significantly higher long-term survival compared to controls. Furthermore, CpG-PBNP-PTT-treated mice reject tumor rechallenge. Ex vivo studies confirm these therapeutic responses result from the generation of robust T cell-mediated immunological memory. Consequently, in a synchronous 9464D tumor model, CpG-PBNP-PTT induces complete tumor regression on the treated flank and significantly slows tumor progression on the untreated flank, improving animal survival. These findings demonstrate that localized administration of the CpG-PBNP-PTT nanoimmunotherapy drives potent systemic T cell responses in solid tumors such as NB and therefore has therapeutic implications for NB.

GD2-directed bispecific trifunctional antibody outperforms dinutuximab beta in a murine model for aggressive metastasized neuroblastoma

Background

Neuroblastoma is the most common extracranial solid tumor of childhood. Patients with high-risk disease undergo extremely aggressive therapy and nonetheless have cure rates below 50%. Treatment with the ch14.18 monoclonal antibody (dinutuximab beta), directed against the GD2 disialoganglioside, improved 5-year event-free survival in high-risk patients when administered in postconsolidation therapy and was recently implemented in standard therapy. Relapse still occurred in 57% of these patients, necessitating new therapeutic options. Bispecific trifunctional antibodies (trAbs) are IgG-like molecules directed against T cells and cancer surface antigens, redirecting T cells (via their CD3 specificity) and accessory immune cells (via their functioning Fc-fragment) toward tumor cells. We sought proof-of-concept for GD2/CD3-directed trAb efficacy against neuroblastoma.

Methods

We used two GD2-specific trAbs differing only in their CD3-binding specificity: EKTOMUN (GD2/human CD3) and SUREK (GD2/mouse Cd3). This allowed trAb evaluation in human and murine experimental settings. Tumor-blind trAb and the ch14.18 antibody were used as controls. A coculture model of human peripheral blood mononuclear cells (PBMCs) and neuroblastoma cell lines was established to evaluate trAb antitumor efficacy by assessing expression of T-cell surface markers for activation, proinflammatory cytokine release and cytotoxicity assays. Characteristics of tumor-infiltrating T cells and response of neuroblastoma metastases to SUREK treatment were investigated in a syngeneic immunocompetent neuroblastoma mouse model mimicking minimal residual disease.

Results

We show that EKTOMUN treatment caused effector cell activation and release of proinflammatory cytokines in coculture with neuroblastoma cell lines. Furthermore, EKTOMUN mediated GD2-dependent cytotoxic effects in human neuroblastoma cell lines in coculture with PBMCs, irrespective of the level of target antigen expression. This effect was dependent on the presence of accessory immune cells. Treatment with SUREK reduced the intratumor Cd4/Cd8 ratio and activated tumor infiltrating T cells in vivo. In a minimal residual disease model for neuroblastoma, we demonstrated that single-agent treatment with SUREK strongly reduced or eliminated neuroblastoma metastases in vivo. SUREK as well as EKTOMUN demonstrated superior tumor control compared with the anti-GD2 antibody, ch14.18.

Conclusions

Here we provide proof-of-concept for EKTOMUN preclinical efficacy against neuroblastoma, presenting this bispecific trAb as a promising new agent to fight neuroblastoma.

Mitochondrial DNA Haplogroups and Susceptibility to Neuroblastoma

BACKGROUND

Neuroblastoma is a childhood malignancy that arises from the developing sympathetic nervous system. Although mitochondrial dysfunctions have been implicated in the pathophysiology of neuroblastoma, the role of mitochondrial DNA (mtDNA) has not been extensively investigated.

METHODS

A total of 2404 Caucasian children diagnosed with neuroblastoma and 9310 ancestry-matched controls were recruited at the Children's Hospital of Philadelphia. The mtDNA haplogroups were identified from SNP array data of two independent cohorts. We conducted a case-control study to explore potential associations of mtDNA haplogroups with the susceptibility of neuroblastoma. The genetic effect of neuroblastoma was measured by odds ratios (ORs) of mitochondrial haplogroups. All tests were two-sided.

RESULTS

Haplogroup K was statistically significantly associated with reduced risk of neuroblastoma in the discovery cohort consisting of 1474 cases and 5699 controls (OR = 0.72, 95% confidence interval [CI] = 0.57 to 0.90; P = 4.8 × 10-3). The association was replicated in an independent cohort (OR = 0.69, 95% CI = 0.53 to 0.92; P = .01) of 930 cases and 3611 controls. Pooled analysis was performed by combining the two data sets. The association remained highly statistically significant after correction for multiple testing (OR = 0.71, 95% CI = 0.59 to 0.84, P = 1.96 × 10-4, Pcorrected = .002). Further analysis focusing on neuroblastoma subtypes indicated haplogroup K was more associated with high-risk neuroblastoma (OR = 0.57, 95% CI = 0.43 to 0.76; P = 1.46 × 10-4) than low-risk and intermediate-risk neuroblastoma.

CONCLUSIONS

Haplogroup K is an independent genetic factor associated with reduced risk of developing neuroblastoma in European descents. These findings provide new insights into the genetic basis of neuroblastoma, implicating mitochondrial DNA encoded proteins in the etiology of neuroblastoma.

Clinical Relevance of CD4 Cytotoxic T Cells in High-Risk Neuroblastoma

Neuroblastoma is the most common extracranial childhood solid tumor. The majority of high-risk neuroblastoma is resistant/refractory to the current high intensity therapy, and the survival of these patients remains poor for the last three decades. To effectively treat these extremely unfavorable neuroblastomas, innovative immunotherapy approaches would be the most promising. In this article, we discuss the identity of tumor-infiltrating effector cells and immunosuppressive cells in high-risk neuroblastoma. Neuroblastoma is unique in that it expresses little or no classical HLA Class I and II. In contrast, high-risk neuroblastomas express the stress-responsive non-classical Class I, HLA-E molecule. HLA-E is the ligand of activating receptors NKG2C/E that are expressed on memory NK cells, CD8+T cells and CD4 CTLs. By examining a comprehensive RNA-seq gene expression dataset, we detected relatively high levels of CD4 expression in high-risk neuroblastoma tissues. The majority of CD4+ cells were CD3+, and thus they were likely tumor-associated CD4+T cells. In addition, high-level of both CD4 and NKG2C/E expression was associated with prolonged survival of the high-risk neuroblastoma patients, but CD8 levels were not, further suggesting that the CD4+ NKG2C/E+ T cells or CD4 CTL conferred cytotoxicity against the neuroblastoma cells. However, this T cell mediated- "protective effect" declined over time, in part due to the progressive formation of immunosuppressive tumor microenvironment. These observations suggest that to improve survival of high-risk neuroblastoma patients, it is essential to gain insights into how to enhance CD4 CTL cytotoxicity and control the immunosuppressive tumor microenvironment during the course of the disease.

Artificial Tumor Microenvironments in Neuroblastoma

Simple Summary

Children with high-risk neuroblastoma have limited therapeutic options poor survival rates. The neuroblastoma tumor microenvironment contributes the lack of response to many interventions so innovative methods are needed to study the effects of the tumor microenvironment on new therapies. In this manuscript, we review the current literature related to the components of the tumor microenvironment and to the use of three-dimensional printing as modality to study cancer. This review highlights the potential for using three-dimensional printing to create an artificial tumor microenvironment in the presence of neuroblastoma to provide improved preclinical testing of novel therapies.

Abstract

In the quest to advance neuroblastoma therapeutics, there is a need to have a deeper understanding of the tumor microenvironment (TME). From extracellular matrix proteins to tumor associated macrophages, the TME is a robust and diverse network functioning in symbiosis with the solid tumor. Herein, we review the major components of the TME including the extracellular matrix, cytokines, immune cells, and vasculature that support a more aggressive neuroblastoma phenotype and encumber current therapeutic interventions. Contemporary treatments for neuroblastoma are the result of traditional two-dimensional culture studies and in vivo models that have been translated to clinical trials. These pre-clinical studies are costly, time consuming, and neglect the study of cofounding factors such as the contributions of the TME. Three-dimensional (3D) bioprinting has become a novel approach to studying adult cancers and is just now incorporating portions of the TME and advancing to study pediatric solid. We review the methods of 3D bioprinting, how researchers have included TME pieces into the prints, and highlight present studies using neuroblastoma. Ultimately, incorporating the elements of the TME that affect neuroblastoma responses to therapy will improve the development of innovative and novel treatments. The use of 3D bioprinting to achieve this aim will prove useful in developing optimal therapies for children with neuroblastoma.

MYC oncogene is associated with suppression of tumor immunity and targeting Myc induces tumor cell immunogenicity for therapeutic whole cell vaccination

Background

MYC oncogene is deregulated in 70% of all human cancers and is associated with multiple oncogenic functions including immunosuppression in the tumor microenvironment. The role of MYC in the immune microenvironment of neuroblastoma and melanoma is investigated and the effect of targeting Myc on immunogenicity of cancer cells is evaluated.

Methods

Immune cell infiltrates and immunogenic pathway signatures in the context of MYCN amplification were analyzed in human neuroblastoma tumors and in metastatic melanoma. Dose response and cell susceptibility to MYC inhibitors (I-BET726 and JQ1) were determined in mouse cell lines. The influence of downregulating Myc in tumor cells was characterized by immunogenic pathway signatures and functional assays. Myc-suppressed tumor cells were used as whole cell vaccines in preclinical neuroblastoma and melanoma models.

Results

Analysis of immune phenotype in human neuroblastoma and melanoma tumors revealed that MYCN or c-MYC amplified tumors respectively are associated with suppressed immune cell infiltrates and functional pathways. Targeting Myc in cancer cells with I-BET726 and JQ1 results in cell cycle arrest and induces cell immunogenicity. Combining vaccination of Myc-inhibited tumor cells with checkpoint inhibition induced robust antitumor immunity and resulted in therapeutic cancer vaccine therapy in mouse neuroblastoma tumors. Despite vigorous antitumor immunity in the mouse melanoma model, upregulation of immunosuppressive pathways enabled tumor escape.

Conclusions

This study demonstrates that the Myc oncogene is an appropriate target for inducing tumor cell immunogenicity and suggests that Myc-suppressed whole tumor cells combined with checkpoint therapy could be used for formulating a personalized therapeutic tumor vaccine.

Prognostic Signature of Immune Genes and Immune-Related LncRNAs in Neuroblastoma: A Study Based on GEO and TARGET Datasets

Background

The prognostic value of immune-related genes and lncRNAs in neuroblastoma has not been elucidated, especially in subgroups with different outcomes. This study aimed to explore immune-related prognostic signatures.

Materials and Methods

Immune-related prognostic genes and lncRNAs were identified by univariate Cox regression analysis in the training set. The top 20 C-index genes and 17 immune-related lncRNAs were included in prognostic model construction, and random forest and the Least Absolute Shrinkage and Selection Operator (LASSO) regression algorithms were employed to select features. The risk score model was constructed and assessed using the Kaplan-Meier plot and the receiver operating characteristic curve. Functional enrichment analysis of the immune-related lncRNAs was conducted using the STRING database.

Results

In GSE49710, five immune genes (CDK4, PIK3R1, THRA, MAP2K2, and ULBP2) were included in the risk score five genes (RS5_G) signature, and eleven immune-related lncRNAs (LINC00260, FAM13A1OS, AGPAT4-IT1, DUBR, MIAT, TSC22D1-AS1, DANCR, MIR137HG, ERC2-IT1, LINC01184, LINC00667) were brought into risk score LncRNAs (RS_Lnc) signature. Patients were divided into high/low-risk score groups by the median. Overall survival and event/progression-free survival time were shortened in patients with high scores, both in training and validation cohorts. The same results were found in subgroups. In grouping ability assessment, the area under the curves (AUCs) in distinguishing different groups ranged from 0.737 to 0.94, better in discriminating MYCN status and high risk in training cohort (higher than 0.9). Multivariate Cox analysis demonstrated that RS5_G and RS_Lnc were the independent risk factors for overall and event/progression-free survival (all p-values <0.001). Correlation analysis showed that RS5_G and RS_Lnc were negatively associated with aDC, CD8+ T cells, but positively correlated with Th2 cells. Functional enrichment analyzes demonstrated that immune-related lncRNAs are mainly enriched in cancer-related pathways and immune-related pathways.

Conclusion

We identified the immune-related prognostic signature RS5_G and RS_Lnc. The predicting and grouping ability is close to being even better than those reported in other studies, especially in subgroups. This study provided prognostic signatures that may help clinicians to choose optimal treatment strategies and showed a new insight for NB treatment. These results need further biological experiments and clinical validation.

An elderly-onset neuroblastoma concomitant with minimal change nephrotic syndrome: the first autopsy case report

BACKGROUND

Neuroblastoma is a well-known embryonal cancer; however, adult-onset neuroblastomas are rare. The systemic symptoms are related to catecholamine excretion or intraabdominal mass effects. Only two cases of neuroblastoma with nephrotic syndrome have previously been reported. We herein present the first autopsy case of neuroblastoma in an older individual associated with minimal change nephrotic syndrome.

CASE PRESENTATION

A 63-year-old man was admitted to our hospital for investigation of general fatigue. His renal function was normal and his urine was negative for protein. A computed tomography scan showed a renal tumor and intraabdominal lymph node swelling. Approximately 4 months after admission, he suddenly developed acute renal failure and severe proteinuria, and hemodialysis was instituted. A computed tomography scan revealed an increase in the size of the renal tumor and lymph nodes. He died 1 month later and an autopsy was performed. The tumor exhibited diffuse proliferation of tumor cells with scant cytoplasm, namely small blue cell tumor with rosette formation. As a result of immunohistochemical study, a neuroblastoma was diagnosed. Despite the patient's severe renal failure, most glomeruli showed no remarkable changes. The tubular epithelium exhibited detachment and vacuolation. Electron microscopic study of the glomeruli showed diffuse effacement of the foot processes. These features indicate a diagnosis of minimal change nephrotic syndrome with acute tubular injury.

CONCLUSIONS

Minimal change nephrotic syndrome is the most common renal manifestation associated with lymphoproliferative malignancies. We here present an extremely rare case of adult-onset neuroblastoma with minimal change nephrotic syndrome.

Low Expression of IL-15 and NKT in Tumor Microenvironment Predicts Poor Outcome of MYCN-Non-Amplified Neuroblastoma

Immune tumor microenvironment (TME) in neuroblastoma (NBL) contributes to tumor behavior and treatment response. T cells and natural killer (NK) cells have been shown to play important roles in the neuroblastoma TME. However, few reports address the clinical relevance of natural killer T cells (NKTs) and interleukin-15 (IL-15), one of the crucial cytokines controlling the activation and expansion of NK/NKT cells, in NBL. In this study, we examined NKT immunoscores and IL-15 expression in both MYCN-amplified and MYCN-non-amplified NBL to correlate with clinical outcomes such as event-free survival (EFS) and overall survival (OS). From Gene Expression Omnibus (GEO) datasets GSE45480 (n = 643) and GSE49711 (n = 493), we found that NKT immunoscore and IL-15 expression were both significantly lower in MYCN-amplified NBL, and similar results were observed using our clinical NBL samples (n = 53). Moreover, NBL patients (GEO dataset GSE49711 and our clinical samples) with both lower NKT immunoscore and IL-15 expression exhibited decreased EFS and OS regardless of MYCN gene amplification status. Multivariate analysis further showed that the combination of low NKT immunoscore and low IL-15 expression level was an independent prognostic factor for poor EFS and OS in our NBL patients. These findings provide the rationale for the development of strategy to incorporate IL-15 and NKT cell therapy into the treatment regimen for neuroblastoma.

Intratumoral immunosuppression profiles in 11q-deleted neuroblastomas provide new potential therapeutic targets

High‐risk neuroblastoma (NB) patients with 11q deletion frequently undergo late but consecutive relapse cycles with fatal outcome. To date, no actionable targets to improve current multimodal treatment have been identified. We analyzed immune microenvironment and genetic profiles of high‐risk NB correlating with 11q immune status. We show in two independent cohorts that 11q‐deleted NB exhibits various immune inhibitory mechanisms, including increased CD4+ resting T cells and M2 macrophages, higher expression of programmed death‐ligand 1, interleukin‐10, transforming growth factor‐beta‐1, and indoleamine 2,3‐dioxygenase 1 (P < 0.05), and also higher chromosomal breakages (P ≤ 0.02) and hemizygosity of immunosuppressive miRNAs than MYCN‐amplified and other 11q‐nondeleted high‐risk NB. We also analyzed benefits of maintenance treatment in 83 high‐risk stage M NB patients focusing on 11q status, either with standard anti‐GD2 immunotherapy (n = 50) or previous retinoic acid‐based therapy alone (n = 33). Immunotherapy associated with higher EFS (50 vs. 30, P = 0.028) and OS (72 vs. 52, P = 0.047) at 3 years in the overall population. Despite benefits from standard anti‐GD2 immunotherapy in high‐risk NB patients, those with 11q deletion still face poor outcome. This NB subgroup displays intratumoral immune suppression profiles, revealing a potential therapeutic strategy with combination immunotherapy to circumvent this immune checkpoint blockade.

Combined immune checkpoint blockade increases CD8+CD28+PD-1+ effector T cells and provides a therapeutic strategy for patients with neuroblastoma

Immune checkpoint therapy has resulted in minimal clinical response in many pediatric cancers. We sought to understand the influence of immune checkpoint inhibition using anti-PD-1 and anti-CTLA-4 antibodies individually, in combination, and after chemotherapy on immune responses in minimal and established murine neuroblastoma models. We also sought to understand the role of the tumor microenvironment (TME) and PD-L1 expression and their alteration post-chemotherapy in our models and human tissues. PD-L1 expression was enriched in human tumor-associated macrophages and up-regulated after chemotherapy. In a murine minimal disease model, single and dual immune checkpoint blockade promoted tumor rejection, improved survival, and established immune memory with long-term anti-tumor immunity against re-challenge. In an established tumor model, only dual immune checkpoint blockade showed efficacy. Interestingly, dual immune checkpoint therapy distinctly influenced adaptive and innate immune responses, with significant increase in CD8⁺CD28⁺PD-1⁺ T cells and inflammatory macrophages (CD11b^hiCD11c⁻F4/80⁺Ly6C^hi) in tumor-draining lymph nodes. Adding chemotherapy before immunotherapy provided significant survival benefit for mice with established tumors receiving anti-PD-1 or dual immune checkpoint blockade. Our findings demonstrate anti-PD-1 and anti-CTLA-4 therapy induces a novel subset of effector T cells, and support administration of induction chemotherapy immediately prior to immune checkpoint blockade in children with high-risk neuroblastoma.

The immune landscape of neuroblastoma: Challenges and opportunities for novel therapeutic strategies in pediatric oncology

Immunotherapy holds great promise for the treatment of pediatric cancers. In neuroblastoma, the recent implementation of anti-GD2 antibody Dinutuximab into the standard of care has improved patient outcomes substantially. However, 5-year survival rates are still below 50% in patients with high-risk neuroblastoma, which has sparked investigations into novel immunotherapeutic approaches. T cell-engaging therapies such as immune checkpoint blockade, antibody-mediated therapy and adoptive T cell therapy have proven remarkably successful in a range of adult cancers but still meet challenges in pediatric oncology. In neuroblastoma, their limited success may be due to several factors. Neuroblastoma displays low immunogenicity due to its low mutational load and lack of MHC-I expression. Tumour infiltration by T and NK cells is especially low in high-risk neuroblastoma and is prognostic for survival. Only a small fraction of tumour-infiltrating lymphocytes shows tumour reactivity. Moreover, neuroblastoma tumours employ a variety of immune evasion strategies, including expression of immune checkpoint molecules, induction of immunosuppressive myeloid and stromal cells, as well as secretion of immunoregulatory mediators, which reduce infiltration and reactivity of immune cells. Overcoming these challenges will be key to the successful implementation of novel immunotherapeutic interventions. Combining different immunotherapies, as well as personalised strategies, may be promising approaches. We will discuss the composition, function and prognostic value of tumour-infiltrating lymphocytes (TIL) in neuroblastoma, reflect on challenges for immunotherapy, including a lack of TIL reactivity and tumour immune evasion strategies, and highlight opportunities for immunotherapy and future perspectives with regard to state-of-the-art developments in the tumour immunology space.

The Tumor Microenvironment in Neuroblastoma: New Players, New Mechanisms of Interaction and New Perspectives

The contribution of the tumor microenvironment (TME) to cancer progression has been well recognized in recent decades. As cancer therapeutic strategies are increasingly precise and include immunotherapies, knowledge of the nature and function of the TME in a tumor becomes essential. Our understanding of the TME in neuroblastoma (NB), the second most common solid tumor in children, has significantly progressed from an initial focus on its Schwannian component to a better awareness of its complex nature, which includes not only immune but also non-immune cells such as cancer-associated fibroblasts (CAFs), the contribution of which to inflammation and interaction with tumor-associated macrophages (TAMs) is now recognized. Recent studies on the TME landscape of NB tumors also suggest significant differences between MYCN-amplified (MYCN-A) and non-amplified (MYCN-NA) tumors, in their content in stromal and inflammatory cells and their immunosuppressive activity. Extracellular vesicles (EVs) released by cells in the TME and microRNAs (miRs) present in their cargo could play important roles in the communication between NB cells and the TME. This review article discusses these new aspects of the TME in NB and the impact that information on the TME landscape in NB will have in the design of precise, biomarker-integrated clinical trials.

Significance of hematopoietic surface antigen CD34 in neuroblastoma prognosis and the genetic landscape of CD34-expressing neuroblastoma CSCs

High-risk neuroblastoma (HR-NB) is branded with hematogenous metastasis, relapses, and dismal long-term survival. Intensification of consolidation therapy with tandem/triple autologous stem cell (SC) rescue (with bone marrow [BM]/peripheral blood [PB] CD34⁺ selection) after myeloablative chemotherapy has improved long-term survival. However, the benefit is limited by the indication of NB cells in CD34⁺ PBSCs, CD34 expression in NB cells, and the risk of reinfusing NB cancer stem cells (NB CSCs) that could lead to post-transplant relapse. We investigated the association of CD34 surface expression (92 patients) with NB evolution/clinical outcomes. CD34 gene-level status in NB was assessed through RNA-Seq data mining (18 cohorts, n, 3324). Genetic landscape of CD34-expressing NB CSCs (CD133⁺CD34⁺) was compared with CD34^- CSCs (CD133⁺CD34^-). RNA-seq data revealed equivocal association patterns of CD34 expression with patient survival. Our immunohistochemistry data revealed definite, but rare (mean, 0.73%; range 0.00-7.87%; median, 0.20%) CD34 positivity in NB. CD34⁺ significantly associated with MYCN amplification (p, 0.003), advanced disease stage (p, 0.016), and progressive disease (PD, p < 0.0009) after clinical therapy. A general high-is-worse tendency was observed in patients with relapsed disease. High CD34⁺ correlated with poor survival in patients with N-MYC-amplified HR-NB. Gene expression analysis of CD34⁺-NB CSCs identified significant up (4631) and downmodulation (4678) of genes compared with NB CSCs that lack CD34. IPA recognized the modulation of crucial signaling elements (EMT, stemness maintenance, differentiation, inflammation, clonal expansion, drug resistance, metastasis) that orchestrate NB disease evolution in CD34⁺ CSCs compared with CD34^- CSCs. While the function of CD34 in NB evolution requires further in-depth investigation, careful consideration should be exercised for autologous stem cell rescue with CD34⁺ selection in NB patients. Graphical abstract.

MYCN drives chemoresistance in small cell lung cancer while USP7 inhibition can restore chemosensitivity

Small cell lung cancer (SCLC) is an aggressive neuroendocrine cancer characterized by initial chemosensitivity followed by emergence of chemoresistant disease. To study roles for MYCN amplification in SCLC progression and chemoresistance, we developed a genetically engineered mouse model of MYCN-overexpressing SCLC. In treatment-naïve mice, MYCN overexpression promoted cell cycle progression, suppressed infiltration of cytotoxic T cells, and accelerated SCLC. MYCN overexpression also suppressed response to cisplatin-etoposide chemotherapy, with similar findings made upon MYCL overexpression. We extended these data to genetically perturb chemosensitive patient-derived xenograft (PDX) models of SCLC. In chemosensitive PDX models, overexpression of either MYCN or MYCL also conferred a switch to chemoresistance. To identify therapeutic strategies for MYCN-overexpressing SCLC, we performed a genome-scale CRISPR-Cas9 sgRNA screen. We identified the deubiquitinase USP7 as a MYCN-associated synthetic vulnerability. Pharmacological inhibition of USP7 resensitized chemoresistant MYCN-overexpressing PDX models to chemotherapy in vivo. Our findings show that MYCN overexpression drives SCLC chemoresistance and provide a therapeutic strategy to restore chemosensitivity.

Microsatellite instability status differentially associates with intratumoral immune microenvironment in human cancers

Based on clinical outcomes in colorectal cancer, high microsatellite instability (MSI-H) has recently been approved by the Food and Drug Administration (FDA) as a genetic test to select patients for immunotherapy targeting PD-1 and/or CTLA-4 without limitation to cancer type. However, it is unclear whether the MSI-H would broadly alter the tumor microenvironment to confer the therapeutic response of different cancer types to immunotherapy. To fill in this gap, we performed an in silico analysis of tumor immunity among different MSI statuses in five cancer types. We found that consistent with clinical responses to immunotherapy, MSI-H and non-MSI-H samples from colorectal cancer (COAD-READ) exhibited distinct infiltration levels and immune phenotypes. Surprisingly, the immunological difference between MSI-H and non-MSI-H samples was diminished in stomach adenocarcinoma and esophageal carcinoma (STAD-ESCA) and completely disappeared in uterine corpus endometrial carcinoma (UCEC). Regardless of cancer types, the abundance of tumor-infiltrating immune cells, rather than MSI status, strongly associated with the clinical outcome. Since preexisting antitumor immune response in the tumor (hot cancer) is accepted as a prerequisite to the therapeutic response to anti-PD-1/CTLA-4 immunotherapy, our data demonstrate that the impact of MSI varied on immune contexture will lead to the further evaluation of predictive immunotherapy responsiveness based on the universal biomarker of MSI status.

5-aza-2'-Deoxycytidine Induces a RIG-I-Related Innate Immune Response by Modulating Mitochondria Stress in Neuroblastoma

Background: Neuroblastoma (NB) is one of the most common malignant solid tumors to occur in children, characterized by a wide range of genetic and epigenetic aberrations. We studied whether modifications of the latter with a 5-aza-2′-deoxycytidine (decitabine, Dac) DNA methyltransferase inhibitor can provide a therapeutic advantage in NB. Methods: NB cells with or without MYCN amplification were treated with Dac. We used flow cytometry to measure cell apoptosis and death and mitochondrial reactive oxygen species (mtROS), microarray to analyze gene expression profile and bisulfite pyrosequencing to determine the methylation level of the DDX58/RIG-I promoter. Western blot was used to detect markers related to innate immune response and apoptotic signaling, while immunofluorescent imaging was used to determine dsRNA. We generated mtDNA depleted ρ⁰ cells using long-term exposure to low-dose ethidium bromide. Results: Dac preferentially induced a RIG-I-predominant innate immune response and cell apoptosis in SK-N-AS NB cells, significantly reduced the methylation level of the DDX58/RIG-I promoter and increased dsRNA accumulation in the cytosol. Dac down regulated mitochondrial genes related to redox homeostasis, but augmented mtROS production. ρ⁰ cells demonstrated a blunted response in innate immune response and apoptotic cell death, as well as greatly diminished dsRNA. The response of NB cells to CDDP and poly(I:C) was potentiated by Dac in association with increased mtROS, which was blunted in ρ⁰ cells. Conclusions: This study indicates that Dac effectively induces a RIG-I-related innate immune response and apoptotic signaling primarily in SK-N-AS NB cells by hypomethylating DDX58/RIG-I promoter, elevated mtROS and increased dsRNA. Dac can potentiate the cytotoxic effects of CDDP and poly(I:C) in NB cells.

Exploration of the molecular characteristics of the tumor-immune interaction and the development of an individualized immune prognostic signature for neuroblastoma

Neuroblastoma (NBL) exists in a complex tumor-immune microenvironment. Immune cell infiltration and tumor-immune molecules play a critical role in tumor development and significantly impact the prognosis of patients. However, the molecular characteristics describing the NBL-immune interaction and their prognostic potential have yet to be investigated systematically. We first employed multiple machine learning algorithms, such as Gene Sets Enrichment Analysis and cell type identification by estimating relative subsets of RNA transcripts, to identify immunophenotypes and immunological characteristics in NBL patient data from public databases and then investigated the prognostic potential and regulatory networks of identified immune-related genes involved in the NBL-immune interaction. The immunity signature combining nine immunity genes was confirmed as more effective for individual risk stratification and survival outcome prediction in NBL patients than common clinical characteristics (area under the curve [AUC] = 0.819, C-index = 0.718, p < .001). A mechanistic exploration revealed the regulatory network of molecules involved in the NBL-immune interaction. These immune molecules were also discovered to possess a significant correlation with plasma cell infiltration, MYCN status, and the level of chemokines and macrophage-related molecules (p < .001). A nomogram was constructed based on the immune signature and clinical characteristics, which showed high potential for prognosis prediction (AUC = 0.856, C-index = 0.755, p < .001). We systematically elucidated the complex regulatory mechanisms and characteristics of the molecules involved in the NBL-immune interaction and their prognostic potential, which may have important implications for further understanding the molecular mechanism of the NBL-immune interaction and identifying high-risk NBL patients to guide clinical treatment.

Formation of Human Neuroblastoma in Mouse-Human Neural Crest Chimeras

Neuroblastoma (NB), derived from the neural crest (NC), is the most common pediatric extracranial solid tumor. Here, we establish a platform that allows the study of human NBs in mouse-human NC chimeras. Chimeric mice were produced by injecting human NC cells carrying NB relevant oncogenes in utero into gastrulating mouse embryos. The mice developed tumors composed of a heterogenous cell population that resembled that seen in primary NBs of patients but were significantly different from homogeneous tumors formed in xenotransplantation models. The human tumors emerged in immunocompetent hosts and were extensively infiltrated by mouse cytotoxic T cells, reflecting a vigorous host anti-tumor immune response. However, the tumors blunted the immune response by inducing infiltration of regulatory T cells and expression of immune-suppressive molecules similar to escape mechanisms seen in human cancer patients. Thus, this experimental platform allows the study of human tumor initiation, progression, manifestation, and tumor-immune-system interactions in an animal model system.

Monitoring Immune Responses in Neuroblastoma Patients during Therapy

Neuroblastoma (NBL) is the most common extracranial solid tumor in childhood. Despite intense treatment, children with this high-risk disease have a poor prognosis. Immunotherapy showed a significant improvement in event-free survival in high-risk NBL patients receiving chimeric anti-GD2 in combination with cytokines and isotretinoin after myeloablative consolidation therapy. However, response to immunotherapy varies widely, and often therapy is stopped due to severe toxicities. Objective markers that help to predict which patients will respond or develop toxicity to a certain treatment are lacking. Immunotherapy guided via immune monitoring protocols will help to identify responders as early as possible, to decipher the immune response at play, and to adjust or develop new treatment strategies. In this review, we summarize recent studies investigating frequency and phenotype of immune cells in NBL patients prior and during current treatment protocols and highlight how these findings are related to clinical outcome. In addition, we discuss potential targets to improve immunogenicity and strategies that may help to improve therapy efficacy. We conclude that immune monitoring during therapy of NBL patients is essential to identify predictive biomarkers to guide patients towards effective treatment, with limited toxicities and optimal quality of life.

Cellular components in tumor microenvironment of neuroblastoma and the prognostic value

Background

Tumor microenvironment (TME) contributes to tumor development, progression, and treatment response. In this study, we detailed the cell composition of the TME in neuroblastoma (NB) and constructed a cell risk score model to predict the prognosis of NB.

Methods

xCell score was calculated through transcriptomic data from the datasets GSE49711 and GSE45480 based on the xCell algorithm. The random forest method was employed to select important features and the coefficient was obtained via multivariate cox regression analysis to construct a prognostic model, and the performance was validated in another two independent datasets, GSE16476 and TARGET-NBL.

Results

We found that both immune and non-immune cells varies significantly in different prognostic groups, and were correlated with survival time. The proposed prognostic cell risk score (pCRS) model we constructed can be an independent prognostic indicator for overall survival (OS) and event-free survival (EFS) (training: OS, HR 1.579, EFS, HR 1.563; validation: OS, HR 1.665, 3.848, EFS, HR 2.203, all p-values < 0.01) and only independent prognostic factor in International Neuroblastoma Risk Group high risk patients (HR 1.339, 3.631; p-value 1.76e–2, 3.71e–5), rather than MYCN amplification. Besides, pCRS model showed good performance in grouping, in discriminating MYCN status, the area under the curve (AUC) was 0.889, 0.933, and 0.861 in GSE49711, GSE45480, and GSE16476, respectively. In separating high risk groups, the AUC was 0.904 in GSE49711.

Conclusion

This study details the cellular components in the TME of NB through gene expression data, the proposed pCRS model might provide a basis for treatment selection of high risk patients or targeting cellular components of TME in NB.

MYCN gene polymorphisms and Wilms tumor susceptibility in Chinese children

BACKGROUND

Wilms tumor, derived from embryonic cells, accounts for a large proportion of pediatric renal tumors. MYCN encoded by MYCN proto-oncogene, a member of the MYC family, is a BHLH transcription factor. It plays a critical role in tumorigenesis and predicts poor clinical outcomes in various types of cancer. However, the role of MYCN remained unclarified in Wilms tumor. In this study, we investigated the association between MYCN gene polymorphisms and Wilms tumor susceptibility.

METHODS

Four MYCN gene polymorphisms (rs57961569 G > A, rs9653226 T > C, rs13034994 A > G, and rs60226897 G > A) were genotyped in 183 cases and 603 controls. Adjusted odds ratios (AORs) and 95% confidence intervals (CIs) were calculated to evaluate the association between MYCN gene polymorphisms and Wilms tumor susceptibility.

RESULTS

Overall, no significant association was found for any of the four MYCN gene polymorphisms. Interestingly, in the stratification analysis, the rs57961569 was found to be associated with decreased Wilms tumor susceptibility in the children older than 18 months (AOR = 0.65, 95% CI = 0.42-1.00, P = .050). Moreover, older children carrying 2-4 risk genotypes were at increased risk of Wilms tumor (OR = 1.55, 95% CI = 1.001-2.40, P = .0497). Haplotype GCAA was shown to significantly increased Wilms tumor risk (AOR = 2.40, 95% CI = 1.12-5.14, P = .024).

CONCLUSION

Our study demonstrated that these MYCN gene polymorphisms might be low penetrant variants in Wilms tumor.

Interleukin-15 Enhances Anti-GD2 Antibody-Mediated Cytotoxicity in an Orthotopic PDX Model of Neuroblastoma

PURPOSE

Immunotherapy with IL2, GM-CSF, and an anti-disialoganglioside (GD2) antibody significantly increases event-free survival in children with high-risk neuroblastoma. However, therapy failure in one third of these patients and IL2-related toxicities pose a major challenge. We compared the immunoadjuvant effects of IL15 with those of IL2 for enhancing antibody-dependent cell-mediated cytotoxicity (ADCC) in neuroblastoma.

EXPERIMENTAL DESIGN

We tested ADCC against neuroblastoma patient-derived xenografts (PDX) in vitro and in vivo and examined the functional and migratory properties of NK cells activated with IL2 and IL15.

RESULTS

In cell culture, IL15-activated NK cells induced higher ADCC against two GD⁺ neuroblastoma PDXs than did IL2-activated NK cells (P < 0.001). This effect was dose-dependent (P < 0.001) and was maintained across several effector-to-tumor ratios. As compared with IL2, IL15 also improved chemotaxis of NK cells, leading to higher numbers of tumorsphere-infiltrating NK cells in vitro (P = 0.002). In an orthotopic PDX model, animals receiving chemoimmunotherapy with an anti-GD2 antibody, GM-CSF, and a soluble IL15/IL15Rα complex had greater tumor regression than did those receiving chemotherapy alone (P = 0.012) or combined with anti-GD2 antibody and GM-CSF with (P = 0.016) or without IL2 (P = 0.035). This was most likely due to lower numbers of immature tumor-infiltrating NK cells (DX5⁺CD27⁺) after IL15/IL15Rα administration (P = 0.029) and transcriptional upregulation of Gzmd.

CONCLUSIONS

The substitution of IL15 for IL2 leads to significant tumor regression in vitro and in vivo and supports clinical testing of IL15 for immunotherapy in pediatric neuroblastoma.

Integrative analysis identifies lincRNAs up- and downstream of neuroblastoma driver genes

Long intergenic non-coding RNAs (lincRNAs) are emerging as integral components of signaling pathways in various cancer types. In neuroblastoma, only a handful of lincRNAs are known as upstream regulators or downstream effectors of oncogenes. Here, we exploit RNA sequencing data of primary neuroblastoma tumors, neuroblast precursor cells, neuroblastoma cell lines and various cellular perturbation model systems to define the neuroblastoma lincRNome and map lincRNAs up- and downstream of neuroblastoma driver genes MYCN, ALK and PHOX2B. Each of these driver genes controls the expression of a particular subset of lincRNAs, several of which are associated with poor survival and are differentially expressed in neuroblastoma tumors compared to neuroblasts. By integrating RNA sequencing data from both primary tumor tissue and cancer cell lines, we demonstrate that several of these lincRNAs are expressed in stromal cells. Deconvolution of primary tumor gene expression data revealed a strong association between stromal cell composition and driver gene status, resulting in differential expression of these lincRNAs. We also explored lincRNAs that putatively act upstream of neuroblastoma driver genes, either as presumed modulators of driver gene activity, or as modulators of effectors regulating driver gene expression. This analysis revealed strong associations between the neuroblastoma lincRNAs MIAT and MEG3 and MYCN and PHOX2B activity or expression. Together, our results provide a comprehensive catalogue of the neuroblastoma lincRNome, highlighting lincRNAs up- and downstream of key neuroblastoma driver genes. This catalogue forms a solid basis for further functional validation of candidate neuroblastoma lincRNAs.