Childhood neuroblastoma has a remarkable variability in outcome. Age at diagnosis is one of the most important prognostic factors, with children less than 1 year old having favorable outcomes. Here we study single-cell and single-nuclei transcriptomes of neuroblastoma with different clinical risk groups and stages, including healthy adrenal gland. We compare tumor cell populations with embryonic mouse sympatho-adrenal derivatives, and post-natal human adrenal gland. We provide evidence that low and high-risk neuroblastoma have different cell identities, representing two disease entities. Low-risk neuroblastoma presents a transcriptome that resembles sympatho- and chromaffin cells, whereas malignant cells enriched in high-risk neuroblastoma resembles a subtype of TRKB+ cholinergic progenitor population identified in human post-natal gland. Analyses of these populations reveal different gene expression programs for worst and better survival in correlation with age at diagnosis. Our findings reveal two cellular identities and a composition of human neuroblastoma tumors reflecting clinical heterogeneity and outcome.

Neuroblastoma is an aggressive pediatric cancer with a high rate of metastasis to the BM. Despite intensive treatments including high-dose chemotherapy, the overall survival rate for children with metastatic neuroblastoma remains dismal. Understanding the cellular and molecular mechanisms of the metastatic tumor microenvironment is crucial for developing new therapies and improving clinical outcomes. Here, we used single-cell RNA-Seq to characterize immune and tumor cell alterations in neuroblastoma BM metastases by comparative analysis with patients without metastases. Our results reveal remodeling of the immune cell populations and reprogramming of gene expression profiles in the metastatic niche. In particular, within the BM metastatic niche, we observed the enrichment of immune cells, including tumor-associated neutrophils, macrophages, and exhausted T cells, as well as an increased number of Tregs and a decreased number of B cells. Furthermore, we highlighted cell communication between tumor cells and immune cell populations, and we identified prognostic markers in malignant cells that are associated with worse clinical outcomes in 3 independent neuroblastoma cohorts. Our results provide insight into the cellular, compositional, and transcriptional shifts underlying neuroblastoma BM metastases that contribute to the development of new therapeutic strategies.

Neuroblastoma is a peripheral nervous system tumor that almost exclusively occurs in young children. Although intensified treatment modalities have led to increased patient survival, the prognosis for patients with high-risk disease is still around 50%, signifying neuroblastoma as a leading cause of cancer-related deaths in children. Neuroblastoma is an embryonal tumor and is shaped by its origin from cells within the neural crest. Hence, neuroblastoma usually presents with a low mutational burden and is, in the majority of cases, driven by epigenetically deregulated transcription networks. The recent development of Omic techniques has given us detailed knowledge of neuroblastoma evolution, heterogeneity, and plasticity, as well as intra- and intercellular molecular communication networks within the neuroblastoma microenvironment. Here, we discuss the potential of these recent discoveries with emphasis on new treatment modalities, including immunotherapies which hold promise for better future treatment regimens.

Deregulation of the MYC family of transcription factors c-MYC (encoded by MYC), MYCN, and MYCL is prevalent in most human cancers, with an impact on tumor initiation and progression, as well as response to therapy. In neuroblastoma (NB), amplification of the MYCN oncogene and over-expression of MYC characterize approximately 40% and 10% of all high-risk NB cases, respectively. However, the mechanism and stage of neural crest development in which MYCN and c-MYC contribute to the onset and/or progression of NB are not yet fully understood. Here, we hypothesized that subtle differences in the expression of MYCN and/or c-MYC targets could more accurately stratify NB patients in different risk groups rather than using the expression of either MYC gene alone. We employed an integrative approach using the transcriptome of 498 NB patients from the SEQC cohort and previously defined c-MYC and MYCN target genes to model a multigene transcriptional risk score. Our findings demonstrate that defined sets of c-MYC and MYCN targets with significant prognostic value, effectively stratify NB patients into different groups with varying overall survival probabilities. In particular, patients exhibiting a high-risk signature score present unfavorable clinical parameters, including increased clinical risk, higher INSS stage, MYCN amplification, and disease progression. Notably, target genes with prognostic value differ between c-MYC and MYCN, exhibiting distinct expression patterns in the developing sympathoadrenal system. Genes associated with poor outcomes are mainly found in sympathoblasts rather than in chromaffin cells during the sympathoadrenal development.

To generate energy efficiently, the cell is uniquely challenged to co-ordinate the abundance of electron transport chain protein subunits expressed from both nuclear and mitochondrial genomes. How an effective stoichiometry of this many constituent subunits is co-ordinated post-transcriptionally remains poorly understood. Here we show that Cerox1 , an unusually abundant cytoplasmic long noncoding RNA (lncRNA), modulates the levels of mitochondrial complex I subunit transcripts in a manner that requires binding to microRNA-488-3p. Increased abundance of Cerox1 cooperatively elevates complex I subunit protein abundance and enzymatic activity, decreases reactive oxygen species production, and protects against the complex I inhibitor rotenone. Cerox1 function is conserved across placental mammals: human and mouse orthologues effectively modulate complex I enzymatic activity in mouse and human cells, respectively. Cerox1 is the first lncRNA demonstrated, to our knowledge, to regulate mitochondrial oxidative phosphorylation (OXPHOS) and, with miR-488-3p, represent novel targets for the modulation of complex I activity.