Tumor to normal single-cell mRNA comparisons reveal a pan-neuroblastoma cancer cell

Developmental origins shape the paediatric cancer genome

In the past two decades, technological advances have brought unprecedented insights into the paediatric cancer genome revealing characteristics distinct from those of adult cancer. Originating from developing tissues, paediatric cancers generally have low mutation burden and are driven by variants that disrupt the transcriptional activity, chromatin state, non-coding cis-regulatory regions and other biological functions. Within each tumour, there are multiple populations of cells with varying states, and the lineages of some can be tracked to their fetal origins. Genome-wide genetic screening has identified vulnerabilities associated with both the cell of origin and transcription deregulation in paediatric cancer, which have become a valuable resource for designing new therapeutic approaches including those for small molecules, immunotherapy and targeted protein degradation. In this Review, we present recent findings on these facets of paediatric cancer from a pan-cancer perspective and provide an outlook on future investigations.

Small-molecule inhibition of the METTL3/METTL14 complex suppresses neuroblastoma tumor growth and promotes differentiation

The N6-methyladenosine (m6A) RNA modification is an important regulator of gene expression. m6A is deposited by a methyltransferase complex that includes methyltransferase-like 3 (METTL3) and methyltransferase-like 14 (METTL14). High levels of METTL3/METTL14 drive the growth of many types of adult cancer, and METTL3/METTL14 inhibitors are emerging as new anticancer agents. However, little is known about the m6A epitranscriptome or the role of the METTL3/METTL14 complex in neuroblastoma, a common pediatric cancer. Here, we show that METTL3 knockdown or pharmacologic inhibition with the small molecule STM2457 leads to reduced neuroblastoma cell proliferation and increased differentiation. These changes in neuroblastoma phenotype are associated with decreased m6A deposition on transcripts involved in nervous system development and neuronal differentiation, with increased stability of target mRNAs. In preclinical studies, STM2457 treatment suppresses the growth of neuroblastoma tumors in vivo. Together, these results support the potential of METTL3/METTL14 complex inhibition as a therapeutic strategy against neuroblastoma.

Transcriptional repression of the oncofetal LIN28B gene by the transcription factor SOX6

The identification of regulatory networks contributing to fetal/adult gene expression switches is a major challenge in developmental biology and key to understand the aberrant proliferation of cancer cells, which often reactivate fetal oncogenes. One key example is represented by the developmental gene LIN28B, whose aberrant reactivation in adult tissues promotes tumor initiation and progression. Despite the prominent role of LIN28B in development and cancer, the mechanisms of its transcriptional regulation are largely unknown. Here, by using quantitative RT-PCR and single cell RNA sequencing data, we show that in erythropoiesis the expression of the transcription factor SOX6 matched a sharp decline of LIN28B mRNA during human embryo/fetal to adult globin switching. SOX6 overexpression repressed LIN28B not only in a panel of fetal-like erythroid cells (K562, HEL and HUDEP1; ≈92% p < 0.0001, 54% p = 0.0009 and ≈60% p < 0.0001 reduction, respectively), but also in hepatoblastoma HepG2 and neuroblastoma SH-SY5H cells (≈99% p < 0.0001 and ≈59% p < 0.0001 reduction, respectively). SOX6-mediated repression caused downregulation of the LIN28B/Let-7 targets, including MYC and IGF2BP1, and rapidly blocks cell proliferation. Mechanistically, Lin28B repression is accompanied by SOX6 physical binding within its locus, suggesting a direct mechanism of LIN28B downregulation that might contribute to the fetal/adult erythropoietic transition and restrict cancer proliferation.

A human neural crest model reveals the developmental impact of neuroblastoma-associated chromosomal aberrations

Early childhood tumours arise from transformed embryonic cells, which often carry large copy number alterations (CNA). However, it remains unclear how CNAs contribute to embryonic tumourigenesis due to a lack of suitable models. Here we employ female human embryonic stem cell (hESC) differentiation and single-cell transcriptome and epigenome analysis to assess the effects of chromosome 17q/1q gains, which are prevalent in the embryonal tumour neuroblastoma (NB). We show that CNAs impair the specification of trunk neural crest (NC) cells and their sympathoadrenal derivatives, the putative cells-of-origin of NB. This effect is exacerbated upon overexpression of MYCN, whose amplification co-occurs with CNAs in NB. Moreover, CNAs potentiate the pro-tumourigenic effects of MYCN and mutant NC cells resemble NB cells in tumours. These changes correlate with a stepwise aberration of developmental transcription factor networks. Together, our results sketch a mechanistic framework for the CNA-driven initiation of embryonal tumours.

Targeting c-Myc transactivation by LMNA inhibits tRNA processing essential for malate-aspartate shuttle and tumour progression

BACKGROUND

A series of studies have demonstrated the emerging involvement of transfer RNA (tRNA) processing during the progression of tumours. Nevertheless, the roles and regulating mechanisms of tRNA processing genes in neuroblastoma (NB), the prevalent malignant tumour outside the brain in children, are yet unknown.

METHODS

Analysis of multi-omics results was conducted to identify crucial regulators of downstream tRNA processing genes. Co-immunoprecipitation and mass spectrometry methods were utilised to measure interaction between proteins. The impact of transcriptional regulators on expression of downstream genes was measured by dual-luciferase reporter, chromatin immunoprecipitation, western blotting and real-time quantitative reverse transcription-polymerase chain reaction (RT-PCR) methods. Studies have been conducted to reveal impact and mechanisms of transcriptional regulators on biological processes of NB. Survival differences were analysed using the log-rank test.

RESULTS

c-Myc was identified as a transcription factor driving tRNA processing gene expression and subsequent malate-aspartate shuttle (MAS) in NB cells. Mechanistically, c-Myc directly promoted the expression of glutamyl-prolyl-tRNA synthetase (EPRS) and leucyl-tRNA synthetase (LARS), resulting in translational up-regulation of glutamic-oxaloacetic transaminase 1 (GOT1) as well as malate dehydrogenase 1 (MDH1) via inhibiting general control nonrepressed 2 or activating mechanistic target of rapamycin signalling. Meanwhile, lamin A (LMNA) inhibited c-Myc transactivation via physical interaction, leading to suppression of MAS, aerobic glycolysis, tumourigenesis and aggressiveness. Pre-clinically, lobeline was discovered as a LMNA-binding compound to facilitate its interaction with c-Myc, which inhibited aminoacyl-tRNA synthetase expression, MAS and tumour progression of NB, as well as growth of organoid derived from c-Myc knock-in mice. Low levels of LMNA or elevated expression of c-Myc, EPRS, LARS, GOT1 or MDH1 were linked to a worse outcome and a shorter survival time of clinical NB patients.

CONCLUSIONS

These results suggest that targeting c-Myc transactivation by LMNA inhibits tRNA processing essential for MAS and tumour progression.

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models. Here, we generated single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We developed an unsupervised machine learning approach ('automatic consensus nonnegative matrix factorization' (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirmed a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly however, this weak-mesenchymal-like program was maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 hours, suggesting an uncharacterized therapy-escape mechanism. Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

RBM39 degrader invigorates natural killer cells to eradicate neuroblastoma despite cancer cell plasticity

The cellular plasticity of neuroblastoma is defined by a mixture of two major cell states, adrenergic (ADRN) and mesenchymal (MES), which may contribute to therapy resistance. However, how neuroblastoma cells switch cellular states during therapy remains largely unknown and how to eradicate neuroblastoma regardless of their cell states is a clinical challenge. To better understand the lineage switch of neuroblastoma in chemoresistance, we comprehensively defined the transcriptomic and epigenetic map of ADRN and MES types of neuroblastomas using human and murine models treated with indisulam, a selective RBM39 degrader. We showed that cancer cells not only undergo a bidirectional switch between ADRN and MES states, but also acquire additional cellular states, reminiscent of the developmental pliancy of neural crest cells. The lineage alterations are coupled with epigenetic reprogramming and dependency switch of lineage-specific transcription factors, epigenetic modifiers and targetable kinases. Through targeting RNA splicing, indisulam induces an inflammatory tumor microenvironment and enhances anticancer activity of natural killer cells. The combination of indisulam with anti-GD2 immunotherapy results in a durable, complete response in high-risk transgenic neuroblastoma models, providing an innovative, rational therapeutic approach to eradicate tumor cells regardless of their potential to switch cell states.

Conditional c-MYC activation in catecholaminergic cells drives distinct neuroendocrine tumors: neuroblastoma vs somatostatinoma

The MYC proto-oncogenes (c-MYC, MYCN , MYCL ) are among the most deregulated oncogenic drivers in human malignancies including high-risk neuroblastoma, 50% of which are MYCN -amplified. Genetically engineered mouse models (GEMMs) based on the MYCN transgene have greatly expanded the understanding of neuroblastoma biology and are powerful tools for testing new therapies. However, a lack of c-MYC-driven GEMMs has hampered the ability to better understand mechanisms of neuroblastoma oncogenesis and therapy development given that c-MYC is also an important driver of many high-risk neuroblastomas. In this study, we report two transgenic murine neuroendocrine models driven by conditional c-MYC induction in tyrosine hydroxylase (Th) and dopamine β-hydroxylase (Dbh)-expressing cells. c-MYC induction in Th-expressing cells leads to a preponderance of Pdx1 + somatostatinomas, a type of pancreatic neuroendocrine tumor (PNET), resembling human somatostatinoma with highly expressed gene signatures of δ cells and potassium channels. In contrast, c-MYC induction in Dbh-expressing cells leads to onset of neuroblastomas, showing a better transforming capacity than MYCN in a comparable C57BL/6 genetic background. The c-MYC murine neuroblastoma tumors recapitulate the pathologic and genetic features of human neuroblastoma, express GD2, and respond to anti-GD2 immunotherapy. This model also responds to DFMO, an FDA-approved inhibitor targeting ODC1, which is a known MYC transcriptional target. Thus, establishing c-MYC-overexpressing GEMMs resulted in different but related tumor types depending on the targeted cell and provide useful tools for testing immunotherapies and targeted therapies for these diseases.

A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

Differentiated neuroblastoma cells remain epigenetically poised for de-differentiation to an immature state

Neuroblastoma is the most common extracranial solid tumor of childhood and accounts for a significant share of childhood cancer deaths. Prior studies utilizing RNA sequencing of bulk tumor populations showed two predominant cell states characterized by high and low expression of neuronal genes. Although cells respond to treatment by altering their gene expression, it is unclear whether this reflects shifting balances of distinct subpopulations or plasticity of individual cells. Using mouse and human neuroblastoma cell lines lacking MYCN amplification, we show that the antigen CD49b (also known as ITGA2) distinguishes these subpopulations. CD49b expression marked proliferative cells with an immature gene expression program, whereas CD49b-negative cells expressed differentiated neuronal marker genes and were non-cycling. Sorted populations spontaneously switched between CD49b expression states in culture, and CD49b-negative cells could generate rapidly growing, CD49b-positive tumors in mice. Although treatment with the chemotherapy drug doxorubicin selectively killed CD49b-positive cells in culture, the CD49b-positive population recovered when treatment was withdrawn. We profiled histone 3 (H3) lysine 27 acetylation (H3K27ac) to identify enhancers and super enhancers that were specifically active in each population and found that CD49b-negative cells maintained the priming H3 lysine 4 methylation (H3K4me1) mark at elements that were active in cells with high expression of CD49b. Improper maintenance of primed enhancer elements might thus underlie cellular plasticity in neuroblastoma, representing potential therapeutic targets for this lethal tumor.

Pan-cancer landscape of epigenetic factor expression predicts tumor outcome

Oncogenic pathways that drive cancer progression reflect both genetic changes and epigenetic regulation. Here we stratified primary tumors from each of 24 TCGA adult cancer types based on the gene expression patterns of epigenetic factors (epifactors). The tumors for five cancer types (ACC, KIRC, LGG, LIHC, and LUAD) separated into two robust clusters that were better than grade or epithelial-to-mesenchymal transition in predicting clinical outcomes. The majority of epifactors that drove the clustering were also individually prognostic. A pan-cancer machine learning model deploying epifactor expression data for these five cancer types successfully separated the patients into poor and better outcome groups. Single-cell analysis of adult and pediatric tumors revealed that expression patterns associated with poor or worse outcomes were present in individual cells within tumors. Our study provides an epigenetic map of cancer types and lays a foundation for discovering pan-cancer targetable epifactors.

Ewing Sarcoma Single-cell Transcriptome Analysis Reveals Functionally Impaired Antigen-presenting Cells

Novel therapeutic strategies are urgently needed for patients with high-risk Ewing sarcoma and for the reduction of severe side effects for all patients. Immunotherapy may fill this need, but its successful application has been hampered by a lack of knowledge on the composition and function of the Ewing sarcoma immune microenvironment. Here, we explore the immune microenvironment of Ewing sarcoma, by single-cell RNA sequencing of 18 Ewing sarcoma primary tissue samples. Ewing sarcoma is infiltrated by natural killer, T, and B cells, dendritic cells, and immunosuppressive macrophages. Ewing sarcoma-associated T cells show various degrees of dysfunction. The antigen-presenting cells found in Ewing sarcoma lack costimulatory gene expression, implying functional impairment. Interaction analysis reveals a clear role for Ewing sarcoma tumor cells in turning the Ewing sarcoma immune microenvironment into an immunosuppressive niche. These results provide novel insights into the functional state of immune cells in the Ewing sarcoma tumor microenvironment and suggest mechanisms by which Ewing sarcoma tumor cells interact with, and shape, the immune microenvironment.

SIGNIFICANCE

This study is the first presenting a detailed analysis of the Ewing sarcoma microenvironment using single-cell RNA sequencing. We provide novel insight into the functional state of immune cells and suggests mechanisms by which Ewing tumor cells interact with, and shape, their immune microenvironment. These insights provide help in understanding the failures and successes of immunotherapy in Ewing sarcoma and may guide novel targeted (immuno) therapeutic approaches.

Two bullets in the gun: combining immunotherapy with chemotherapy to defeat neuroblastoma by targeting adrenergic-mesenchymal plasticity

Neuroblastoma (NB) is a childhood tumor that originates in the peripheral sympathetic nervous system and is responsible for 15% of cancer-related deaths in the pediatric population. Despite intensive multimodal treatment, many patients with high-risk NB relapse and develop a therapy-resistant tumor. One of the phenomena related to therapeutic resistance is intratumor heterogeneity resulting from the adaptation of tumor cells in response to different selective environmental pressures. The transcriptional and epigenetic profiling of NB tissue has recently revealed the existence of two distinct cellular identities in the NB, termed adrenergic (ADRN) and mesenchymal (MES), which can spontaneously interconvert through epigenetic regulation. This phenomenon, known as tumor plasticity, has a major impact on cancer pathogenesis. The aim of this review is to describe the peculiarities of these two cell states, and how their plasticity affects the response to current therapeutic treatments, with special focus on the immunogenic potential of MES cells. Furthermore, we will discuss the opportunity to combine immunotherapy with chemotherapy to counteract NB phenotypic interconversion.

Pathway to Independence: the future of developmental biology

In 2022, Development launched its Pathway to Independence (PI) Programme, aimed at supporting postdocs as they transition to their first independent position. We selected eight talented researchers as the first cohort of PI Fellows. In this article, each of our Fellows provides their perspective on the future of their field. Together, they paint an exciting picture of the current state of and open questions in developmental biology.

Target Genes of c-MYC and MYCN with Prognostic Power in Neuroblastoma Exhibit Different Expressions during Sympathoadrenal Development

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.

Therapeutic targeting of FUBP3 phase separation by GATA2-AS1 inhibits malate-aspartate shuttle and neuroblastoma progression via modulating SUZ12 activity

Malate-aspartate shuttle (MAS) is essential for maintaining glycolysis and energy metabolism in tumors, while its regulatory mechanisms in neuroblastoma (NB), the commonest extracranial malignancy during childhood, still remain to be elucidated. Herein, by analyzing multi-omics data, GATA binding protein 2 (GATA2) and its antisense RNA 1 (GATA2-AS1) were identified to suppress MAS during NB progression. Mechanistic studies revealed that GATA2 inhibited the transcription of glutamic-oxaloacetic transaminase 2 (GOT2) and malate dehydrogenase 2 (MDH2). As a long non-coding RNA destabilized by RNA binding motif protein 15-mediated N6-methyladenosine methylation, GATA2-AS1 bound with far upstream element binding protein 3 (FUBP3) to repress its liquid-liquid phase separation and interaction with suppressor of zest 12 (SUZ12), resulting in decrease of SUZ12 activity and epigenetic up-regulation of GATA2 and other tumor suppressors. Rescue experiments revealed that GATA2-AS1 inhibited MAS and NB progression via repressing interaction between FUBP3 and SUZ12. Pre-clinically, administration of lentivirus carrying GATA2-AS1 suppressed MAS, aerobic glycolysis, and aggressive behaviors of NB xenografts. Notably, low GATA2-AS1 or GATA2 expression and high FUBP3, SUZ12, GOT2 or MDH2 levels were linked with unfavorable outcome of NB patients. These findings suggest that GATA2-AS1 inhibits FUBP3 phase separation to repress MAS and NB progression via modulating SUZ12 activity.

An integrated single-cell analysis of human adrenal cortex development

The adrenal glands synthesize and release essential steroid hormones such as cortisol and aldosterone, but many aspects of human adrenal gland development are not well understood. Here, we combined single-cell and bulk RNA sequencing, spatial transcriptomics, IHC, and micro-focus computed tomography to investigate key aspects of adrenal development in the first 20 weeks of gestation. We demonstrate rapid adrenal growth and vascularization, with more cell division in the outer definitive zone (DZ). Steroidogenic pathways favored androgen synthesis in the central fetal zone, but DZ capacity to synthesize cortisol and aldosterone developed with time. Core transcriptional regulators were identified, with localized expression of HOPX (also known as Hop homeobox/homeobox-only protein) in the DZ. Potential ligand-receptor interactions between mesenchyme and adrenal cortex were seen (e.g., RSPO3/LGR4). Growth-promoting imprinted genes were enriched in the developing cortex (e.g., IGF2, PEG3). These findings reveal aspects of human adrenal development and have clinical implications for understanding primary adrenal insufficiency and related postnatal adrenal disorders, such as adrenal tumor development, steroid disorders, and neonatal stress.

Transcriptional signatures associated with persisting CD19 CAR-T cells in children with leukemia

In the context of relapsed and refractory childhood pre-B cell acute lymphoblastic leukemia (R/R B-ALL), CD19-targeting chimeric antigen receptor (CAR)-T cells often induce durable remissions, which requires the persistence of CAR-T cells. In this study, we systematically analyzed CD19 CAR-T cells of 10 children with R/R B-ALL enrolled in the CARPALL trial via high-throughput single-cell gene expression and T cell receptor sequencing of infusion products and serial blood and bone marrow samples up to 5 years after infusion. We show that long-lived CAR-T cells developed a CD4/CD8 double-negative phenotype with an exhausted-like memory state and distinct transcriptional signature. This persistence signature was dominant among circulating CAR-T cells in all children with a long-lived treatment response for which sequencing data were sufficient (4/4, 100%). The signature was also present across T cell subsets and clonotypes, indicating that persisting CAR-T cells converge transcriptionally. This persistence signature was also detected in two adult patients with chronic lymphocytic leukemia with decade-long remissions who received a different CD19 CAR-T cell product. Examination of single T cell transcriptomes from a wide range of healthy and diseased tissues across children and adults indicated that the persistence signature may be specific to long-lived CAR-T cells. These findings raise the possibility that a universal transcriptional signature of clinically effective, persistent CD19 CAR-T cells exists.

Neural crest-related NXPH1/α-NRXN signaling opposes neuroblastoma malignancy by inhibiting organotropic metastasis

Neuroblastoma is a pediatric cancer that can present as low- or high-risk tumors (LR-NBs and HR-NBs), the latter group showing poor prognosis due to metastasis and strong resistance to current therapy. Whether LR-NBs and HR-NBs differ in the way they exploit the transcriptional program underlying their neural crest, sympatho-adrenal origin remains unclear. Here, we identified the transcriptional signature distinguishing LR-NBs from HR-NBs, which consists mainly of genes that belong to the core sympatho-adrenal developmental program and are associated with favorable patient prognosis and with diminished disease progression. Gain- and loss-of-function experiments revealed that the top candidate gene of this signature, Neurexophilin-1 (NXPH1), has a dual impact on NB cell behavior in vivo: whereas NXPH1 and its receptor α-NRXN1 promote NB tumor growth by stimulating cell proliferation, they conversely inhibit organotropic colonization and metastasis. As suggested by RNA-seq analyses, these effects might result from the ability of NXPH1/α-NRXN signalling to restrain the conversion of NB cells from an adrenergic state to a mesenchymal one. Our findings thus uncover a transcriptional module of the sympatho-adrenal program that opposes neuroblastoma malignancy by impeding metastasis, and pinpoint NXPH1/α-NRXN signaling as a promising target to treat HR-NBs.

Single-cell transcriptomics uncovers cellular architecture and developmental trajectories in hepatoblastoma

BACKGROUND AND AIMS

Hepatoblastoma (HB) is the predominant type of childhood liver cancer. Treatment options for the clinically advanced HB remain limited. We aimed to dissect the cellular and molecular basis underlying HB oncogenesis and heterogeneity at the single-cell level, which could facilitate a better understanding of HB at both the biological and clinical levels.

APPROACH AND RESULTS

Single-cell transcriptome profiling of tumor and paired distal liver tissue samples from five patients with HB was performed. Deconvolution analysis was used for integrating the single-cell transcriptomic profiles with the bulk transcriptomes of our HB cohort of post-neoadjuvant chemotherapy tumor samples. A single-cell transcriptomic landscape of early human liver parenchymal development was established for exploring the cellular root and hierarchy of HB oncogenesis. As a result, seven distinct tumor cell subpopulations were annotated, and an effective HB subtyping method was established based on their compositions. A HB tumor cell hierarchy was further revealed to not only fit with the classical cancer stem cell (CSC) model but also mirror the early human liver parenchymal development. Moreover, FACT inhibition, which could disrupt the oncogenic positive feedback loop between MYC and SSRP1 in HB, was identified as a promising epigenetic-targeted therapeutic strategy against the CSC-like HB1-Pro-like1 subpopulation and its related high-risk "Pro-like1" subtype of HB.

CONCLUSIONS

Our findings illustrate the cellular architecture and developmental trajectories of HB via integrative bulk and single-cell transcriptome analyses, thus establishing a resourceful framework for the development of targeted diagnostics and therapeutics in the future.

Identification of a novel gene signature for neuroblastoma differentiation using a Boolean implication network

Although induction of differentiation represents an effective strategy for neuroblastoma treatment, the mechanisms underlying neuroblastoma differentiation are poorly understood. We generated a computational model of neuroblastoma differentiation consisting of interconnected gene clusters identified based on symmetric and asymmetric gene expression relationships. We identified a differentiation signature consisting of series of gene clusters comprised of 1251 independent genes that predicted neuroblastoma differentiation in independent datasets and in neuroblastoma cell lines treated with agents known to induce differentiation. This differentiation signature was associated with patient outcomes in multiple independent patient cohorts and validated the role of MYCN expression as a marker of neuroblastoma differentiation. Our results further identified novel genes associated with MYCN via asymmetric Boolean implication relationships that would not have been identified using symmetric computational approaches and that were associated with both neuroblastoma differentiation and patient outcomes. Our differentiation signature included a cluster of genes involved in intracellular signaling and growth factor receptor trafficking pathways that is strongly associated with neuroblastoma differentiation, and we validated the associations of UBE4B, a gene within this cluster, with neuroblastoma cell and tumor differentiation. Our findings demonstrate that Boolean network analyses of symmetric and asymmetric gene expression relationships can identify novel genes and pathways relevant for neuroblastoma tumor differentiation that could represent potential therapeutic targets.

IGF2BP1 induces neuroblastoma via a druggable feedforward loop with MYCN promoting 17q oncogene expression

BACKGROUND

Neuroblastoma is the most common solid tumor in infants accounting for approximately 15% of all cancer-related deaths. Over 50% of high-risk neuroblastoma relapse, emphasizing the need of novel drug targets and therapeutic strategies. In neuroblastoma, chromosomal gains at chromosome 17q, including IGF2BP1, and MYCN amplification at chromosome 2p are associated with adverse outcome. Recent, pre-clinical evidence indicates the feasibility of direct and indirect targeting of IGF2BP1 and MYCN in cancer treatment.

METHODS

Candidate oncogenes on 17q were identified by profiling the transcriptomic/genomic landscape of 100 human neuroblastoma samples and public gene essentiality data. Molecular mechanisms and gene expression profiles underlying the oncogenic and therapeutic target potential of the 17q oncogene IGF2BP1 and its cross-talk with MYCN were characterized and validated in human neuroblastoma cells, xenografts and PDX as well as novel IGF2BP1/MYCN transgene mouse models.

RESULTS

We reveal a novel, druggable feedforward loop of IGF2BP1 (17q) and MYCN (2p) in high-risk neuroblastoma. This promotes 2p/17q chromosomal gains and unleashes an oncogene storm resulting in fostered expression of 17q oncogenes like BIRC5 (survivin). Conditional, sympatho-adrenal transgene expression of IGF2BP1 induces neuroblastoma at a 100% incidence. IGF2BP1-driven malignancies are reminiscent to human high-risk neuroblastoma, including 2p/17q-syntenic chromosomal gains and upregulation of Mycn, Birc5, as well as key neuroblastoma circuit factors like Phox2b. Co-expression of IGF2BP1/MYCN reduces disease latency and survival probability by fostering oncogene expression. Combined inhibition of IGF2BP1 by BTYNB, MYCN by BRD inhibitors or BIRC5 by YM-155 is beneficial in vitro and, for BTYNB, also.

CONCLUSION

We reveal a novel, druggable neuroblastoma oncogene circuit settling on strong, transcriptional/post-transcriptional synergy of MYCN and IGF2BP1. MYCN/IGF2BP1 feedforward regulation promotes an oncogene storm harboring high therapeutic potential for combined, targeted inhibition of IGF2BP1, MYCN expression and MYCN/IGF2BP1-effectors like BIRC5.

Reversible transitions between noradrenergic and mesenchymal tumor identities define cell plasticity in neuroblastoma

Noradrenergic and mesenchymal identities have been characterized in neuroblastoma cell lines according to their epigenetic landscapes and core regulatory circuitries. However, their relationship and relative contribution in patient tumors remain poorly defined. We now document spontaneous and reversible plasticity between the two identities, associated with epigenetic reprogramming, in several neuroblastoma models. Interestingly, xenografts with cells from each identity eventually harbor a noradrenergic phenotype suggesting that the microenvironment provides a powerful pressure towards this phenotype. Accordingly, such a noradrenergic cell identity is systematically observed in single-cell RNA-seq of 18 tumor biopsies and 15 PDX models. Yet, a subpopulation of these noradrenergic tumor cells presents with mesenchymal features that are shared with plasticity models, indicating that the plasticity described in these models has relevance in neuroblastoma patients. This work therefore emphasizes that intrinsic plasticity properties of neuroblastoma cells are dependent upon external cues of the environment to drive cell identity.

Molecular Evidence for Olfactory Neuroblastoma as a Tumor of Malignant Globose Basal Cells

Olfactory neuroblastoma (ONB, esthesioneuroblastoma) is a sinonasal cancer with an underdeveloped diagnostic toolkit, and is the subject of many incidents of tumor misclassification throughout the literature. Despite its name, connections between the cancer and normal cells of the olfactory epithelium have not been systematically explored and markers of olfactory epithelial cell types are not deployed in clinical practice. Here, we utilize an integrated human-mouse single-cell atlas of the nasal mucosa, including the olfactory epithelium, to identify transcriptomic programs that link ONB to a specific population of stem/progenitor cells known as olfactory epithelial globose basal cells (GBCs). Expression of a GBC transcription factor NEUROD1 distinguishes both low- and high-grade ONB from sinonasal undifferentiated carcinoma, a potential histologic mimic with a distinctly unfavorable prognosis. Furthermore, we identify a reproducible subpopulation of highly proliferative ONB cells expressing the GBC stemness marker EZH2, suggesting that EZH2 inhibition may play a role in the targeted treatment of ONB. Finally, we study the cellular states comprising ONB parenchyma using single-cell transcriptomics and identify evidence of a conserved GBC transcriptional regulatory circuit that governs divergent neuronal-versus-sustentacular differentiation. These results link ONB to a specific cell type for the first time and identify conserved developmental pathways within ONB that inform diagnostic, prognostic, and mechanistic investigation.

Neuroblastoma arises in early fetal development and its evolutionary duration predicts outcome

Neuroblastoma, the most frequent solid tumor in infants, shows very diverse outcomes from spontaneous regression to fatal disease. When these different tumors originate and how they evolve are not known. Here we quantify the somatic evolution of neuroblastoma by deep whole-genome sequencing, molecular clock analysis and population-genetic modeling in a comprehensive cohort covering all subtypes. We find that tumors across the entire clinical spectrum begin to develop via aberrant mitoses as early as the first trimester of pregnancy. Neuroblastomas with favorable prognosis expand clonally after short evolution, whereas aggressive neuroblastomas show prolonged evolution during which they acquire telomere maintenance mechanisms. The initial aneuploidization events condition subsequent evolution, with aggressive neuroblastoma exhibiting early genomic instability. We find in the discovery cohort (n = 100), and validate in an independent cohort (n = 86), that the duration of evolution is an accurate predictor of outcome. Thus, insight into neuroblastoma evolution may prospectively guide treatment decisions.

Adaptive risk-aware sharable and individual subspace learning for cancer survival analysis with multi-modality data

Biomedical multi-modality data (also named multi-omics data) refer to data that span different types and derive from multiple sources in clinical practices (e.g. gene sequences, proteomics and histopathological images), which can provide comprehensive perspectives for cancers and generally improve the performance of survival models. However, the performance improvement of multi-modality survival models may be hindered by two key issues as follows: (1) how to learn and fuse modality-sharable and modality-individual representations from multi-modality data; (2) how to explore the potential risk-aware characteristics in each risk subgroup, which is beneficial to risk stratification and prognosis evaluation. Additionally, learning-based survival models generally refer to numerous hyper-parameters, which requires time-consuming parameter setting and might result in a suboptimal solution. In this paper, we propose an adaptive risk-aware sharable and individual subspace learning method for cancer survival analysis. The proposed method jointly learns sharable and individual subspaces from multi-modality data, whereas two auxiliary terms (i.e. intra-modality complementarity and inter-modality incoherence) are developed to preserve the complementary and distinctive properties of each modality. Moreover, it equips with a grouping co-expression constraint for obtaining risk-aware representation and preserving local consistency. Furthermore, an adaptive-weighted strategy is employed to efficiently estimate crucial parameters during the training stage. Experimental results on three public datasets demonstrate the superiority of our proposed model.

The potential mechanisms of neuroblastoma in children based on bioinformatics big data

BACKGROUND

In recent years, miRNAs have become a research hotspot, which is related to the occurrence and development of a variety of malignant tumors, but there are few studies in neuroblastoma. In this study, the differentially expressed microRNAs (miRNAs) in neuroblastoma were identified and analyzed using bioinformatics, and their biological functions and related signaling pathways were examined.

METHODS

The neuroblastoma miRNA chip GSE121513 was obtained from the Gene Expression Omnibus (GEO) database and the data of 95 neuroblastoma samples and normal fetal adrenal neuroblastoma samples were analyzed to screen the differential miRNAs. The target genes of the differentially expressed miRNAs were predicted using |log fold change (FC)| ≥4. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses were performed to construct a protein-protein interaction network and identify the core target genes.

RESULTS

A total of 91 differentially expressed miRNAs were identified (P<0.05, |logFC| ≥1), including 52 upregulated and 39 downregulated miRNAs. The target genes of the differential miRNAs (P<0.05, |logFC| ≥4) were pretested, and 602 target genes were obtained. Functional analysis showed that these genes were mainly located in the extracellular matrix region of proteins, and were involved in the negative regulation of cytoplasmic translation, mRNA 3'-untranslated region (UTR) binding, and binding to nucleic acid to inhibit the activity of translation factors. They were also involved in RNA degradation, adhesion pathways, and the phosphatidylinositol-3-kinase (PI3K)-Akt signaling pathway. Ten key target genes were identified via protein interaction network screening.

CONCLUSIONS

The differential miRNAs may be related to the occurrence of neuroblastoma were screened.

Defining the origin, evolution, and immune composition of SDH-deficient renal cell carcinoma

Succinate dehydrogenase (SDH)-deficient renal cell carcinoma represents a rare subtype of hereditary kidney cancer. Clinical diagnosis can be challenging and there is little evidence to guide systemic therapeutic options. We performed genomic profiling of a cohort of tumors through the analysis of whole genomes, transcriptomes, as well as flow cytometry and immunohistochemistry in order to gain a deeper understanding of their molecular biology. We find neutral evolution after early tumor activation with a lack of secondary driver events. We show that these tumors have epithelial derivation, possibly from the macula densa, a specialized paracrine cell of the renal juxtaglomerular apparatus. They subsequently develop into immune excluded tumors. We provide transcriptomic and protein expression evidence of a highly specific tumor marker, PAPPA2. These translational findings have implications for the diagnosis and treatment for this rare tumor subtype.

Clinically relevant treatment of PDX models reveals patterns of neuroblastoma chemoresistance

Chemotherapy resistance and relapses are common in high-risk neuroblastoma (NB). Here, we developed a clinically relevant in vivo treatment protocol mimicking the first-line five-chemotherapy treatment regimen of high-risk NB and applied this protocol to mice with MYCN-amplified NB patient-derived xenografts (PDXs). Genomic and transcriptomic analyses were used to reveal NB chemoresistance mechanisms. Intrinsic resistance was associated with high genetic diversity and an embryonic phenotype. Relapsed NB with acquired resistance showed a decreased adrenergic phenotype and an enhanced immature mesenchymal-like phenotype, resembling multipotent Schwann cell precursors. NBs with a favorable treatment response presented a lineage-committed adrenergic phenotype similar to normal neuroblasts. Novel integrated phenotypic gene signatures reflected treatment response and patient prognosis. NB organoids established from relapsed PDX tumors retained drug resistance, tumorigenicity, and transcriptional cell states. This work sheds light on the mechanisms of NB chemotherapy response and emphasizes the importance of transcriptional cell states in chemoresistance.

Tumoral heterogeneity in neuroblastoma

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

Neuroblastoma: When differentiation goes awry

Neuroblastoma is a leading cause of cancer-related death in children. Accumulated data suggest that differentiation arrest of the neural-crest-derived sympathoadrenal lineage contributes to neuroblastoma formation. The developmental arrest of these cell types explains many biological features of the disease, including its cellular heterogeneity, mutational spectrum, spontaneous regression, and response to drugs that induce tumor cell differentiation. In this review, we provide evidence that supports the notion that arrested neural-crest-derived progenitor cells give rise to neuroblastoma and discuss how this concept could be exploited for clinical management of the disease.

Precise identification of cancer cells from allelic imbalances in single cell transcriptomes

A fundamental step of tumour single cell mRNA analysis is separating cancer and non-cancer cells. We show that the common approach to separation, using shifts in average expression, can lead to erroneous biological conclusions. By contrast, allelic imbalances representing copy number changes directly detect the cancer genotype and accurately separate cancer from non-cancer cells. Our findings provide a definitive approach to identifying cancer cells from single cell mRNA sequencing data.

The identification of cancer cells from single cell transcriptomes can be improved by detecting allelic imbalances due to copy number changes.

Connecting telomere maintenance and regulation to the developmental origin and differentiation states of neuroblastoma tumor cells

A cardinal feature that distinguishes clinically high-risk neuroblastoma from low-risk tumors is telomere maintenance. Specifically, neuroblastoma tumors with either active telomerase or alternative lengthening of telomeres exhibit aggressive growth characteristics that lead to poor outcomes, whereas tumors without telomere maintenance can be managed with observation or minimal treatment. Even though the need for cancer cells to maintain telomere DNA-in order to sustain cell proliferation-is well established, recent studies suggest that the neural crest origin of neuroblastoma may enforce unique relationships between telomeres and tumor malignancy. Specifically in neuroblastoma, telomere structure and telomerase activity are correlated with the adrenergic/mesenchymal differentiation states, and manipulating telomerase activity can trigger tumor cell differentiation. Both findings may reflect features of normal neural crest development. This review summarizes recent advances in the characterization of telomere structure and telomere maintenance mechanisms in neuroblastoma and discusses the findings in the context of relevant literature on telomeres during embryonic and neural development. Understanding the canonical and non-canonical roles of telomere maintenance in neuroblastoma could reveal vulnerabilities for telomere-directed therapies with potential applications to other pediatric malignancies.

Extensive patient-to-patient single nucleus transcriptome heterogeneity in pheochromocytomas and paragangliomas

Pheochromocytoma, neuroendocrine tumor, single cell RNA-sequencing, transcriptome, heterogeneity, SDHB, RET, paraganglinoma; Pheochromocytomas (PC) and paragangliomas (PG) are rare neuroendocrine tumors with varied genetic makeup and are associated with high cardiovascular morbidity and a variable risk of malignancy. The source of the transcriptional heterogeneity of the disease and the underlying biological processes that determine the outcome of PCPG remain largely unclear. We focused on PCPG tumors with germline SDHB and RET mutations, which represent distinct prognostic groups with worse or better prognoses, respectively. We applied single-nuclei RNA sequencing (snRNA-seq) to tissue samples from 11 patients and found high patient-to-patient transcriptome heterogeneity in neuroendocrine tumor cells. The tumor microenvironment also showed heterogeneous profiles, mainly contributed by macrophages of the immune cell clusters and Schwann cells of the stroma. By performing non-negative matrix factorization, we identified common transcriptional programs active in RET and SDHB, as well as distinct modules, including neuronal development, hormone synthesis and secretion, and DNA replication. Similarities between the transcriptomes of the tumor cells and those of the chromaffin- and precursor cell types suggests different developmental stages at which PC and PG tumors appear to be arrested.

Linking neural crest development to neuroblastoma pathology

Although rare, childhood (paediatric) cancers are a major cause of death in young children. Unlike many adult cancers, paediatric cancers, such as neuroblastoma (NB), are developmental diseases that rarely show genetic predispositions. NB is the most common extracranial solid tumour in children, accounting for ∼15% of paediatric cancer deaths. This heterogeneous cancer arises from undifferentiated neural crest-derived progenitor cells. As neural crest cells are multipotent and migratory, they are often considered the embryonic paradigm of cancer stem cells. However, very little is known about the events that trigger tumour initiation and progression. Here, we discuss recent insights into sympathoadrenal lineage specification, as well as genetic factors associated with NB. With this in mind, we consider the molecular underpinnings of NB in the context of developmental trajectories of the neural crest lineage. This allows us to compare distinct subtypes of the disease and gene-function interactions during sensitive phases of neural crest development.

Expression and activation of nuclear hormone receptors result in neuronal differentiation and favorable prognosis in neuroblastoma

BACKGROUND

Neuroblastoma (NB), a childhood tumor derived from the sympathetic nervous system, presents with heterogeneous clinical behavior. While some tumors regress spontaneously without medical intervention, others are resistant to therapy, associated with an aggressive phenotype. MYCN-amplification, frequently occurring in high-risk NB, is correlated with an undifferentiated phenotype and poor prognosis. Differentiation induction has been proposed as a therapeutic approach for high-risk NB. We have previously shown that MYCN maintains an undifferentiated state via regulation of the miR-17 ~ 92 microRNA cluster, repressing the nuclear hormone receptors (NHRs) estrogen receptor alpha (ERα) and the glucocorticoid receptor (GR).

METHODS

Cell viability was determined by WST-1. Expression of differentiation markers was analyzed by Western blot, RT-qPCR, and immunofluorescence analysis. Metabolic phenotypes were studied using Agilent Extracellular Flux Analyzer, and accumulation of lipid droplets by Nile Red staining. Expression of angiogenesis, proliferation, and neuronal differentiation markers, and tumor sections were assessed by immunohistochemistry. Gene expression from NB patient as well as adrenal gland cohorts were analyzed using GraphPad Prism software (v.8) and GSEA (v4.0.3), while pseudo-time progression on post-natal adrenal gland cells from single-nuclei transcriptome data was computed using scVelo.

RESULTS

Here, we show that simultaneous activation of GR and ERα potentiated induction of neuronal differentiation, reduced NB cell viability in vitro, and decreased tumor burden in vivo. This was accompanied by a metabolic reprogramming manifested by changes in the glycolytic and mitochondrial functions and in lipid droplet accumulation. Activation of the retinoic acid receptor alpha (RARα) with all-trans retinoic acid (ATRA) further enhanced the differentiated phenotype as well as the metabolic switch. Single-cell nuclei transcriptome analysis of human adrenal glands indicated a sequential expression of ERα, GR, and RARα during development from progenitor to differentiated chromaffin cells. Further, in silico analysis revealed that patients with higher combined expression of GR, ERα, and RARα mRNA levels had elevated expression of neuronal differentiation markers and a favorable outcome.

CONCLUSION

Together, our findings suggest that combination therapy involving activation of several NHRs could be a promising pharmacological approach for differentiation treatment of NB patients.

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.

•

27 different immune subtypes form a complex neuroblastoma microenvironment

•

Inflammatory monocytes correlate with neuroblastoma disease progression

•

Active NK cells correlate with improved survival in human neuroblastoma

•

T cell exhaustion and Treg activity demonstrate no difference between risk groups

Cellular plasticity in the neural crest and cancer

In vertebrates, neural crest cells (NCCs) are a multipotent embryonic population generating both neural/neuronal and mesenchymal derivatives, and thus the neural crest (NC) is often referred to as the fourth germ layer. NC development is a dynamic process, where NCCs possess substantial plasticity in transcriptional and epigenomic profiles. Recent technical advances in single-cell and low-input sequencing have empowered fine-resolution characterisation of NC development. In this review, we summarise the latest models underlying NC-plasticity acquirement and cell-fate restriction, outline the connections between NC plasticity and NC-derived cancer and envision the new opportunities in studying NC plasticity and its link to cancer.

The need for a standard for informed consent for collection of human fetal material

The ISSCR has developed the Informed Consent Standards for Human Fetal Tissue Donation and Research to promote uniformity and transparency in tissue donation and collection. This standard is designed to assist those working with and overseeing the regulation of such tissue and reassure the wider community and public.

BET and CDK Inhibition Reveal Differences in the Proliferation Control of Sympathetic Ganglion Neuroblasts and Adrenal Chromaffin Cells

Neuroblastoma arising from the adrenal differ from ganglionic neuroblastoma both genetically and clinically, with adrenal tumors being associated with a more severe prognosis. The different tumor properties may be linked to specific tumor founder cells in adrenal and sympathetic ganglia. To address this question, we first set up cultures of mouse sympathetic neuroblasts and adrenal chromaffin cells. These cultures were then treated with various proliferation inhibitors to identify lineage-specific responses. We show that neuroblast and chromaffin cell proliferation was affected by WNT, ALK, IGF1, and PRC2/EZH2 signaling inhibitors to a similar extent. However, differential effects were observed in response to bromodomain and extraterminal (BET) protein inhibitors (JQ1, GSK1324726A) and to the CDK-7 inhibitor THZ1, with BET inhibitors preferentially affecting chromaffin cells, and THZ1 preferentially affecting neuroblasts. The differential dependence of chromaffin cells and neuroblasts on BET and CDK signaling may indicate different mechanisms during tumor initiation in sympathetic ganglia and adrenal.

Identifying tumor cells at the single-cell level using machine learning

Tumors are complex tissues of cancerous cells surrounded by a heterogeneous cellular microenvironment with which they interact. Single-cell sequencing enables molecular characterization of single cells within the tumor. However, cell annotation—the assignment of cell type or cell state to each sequenced cell—is a challenge, especially identifying tumor cells within single-cell or spatial sequencing experiments. Here, we propose ikarus, a machine learning pipeline aimed at distinguishing tumor cells from normal cells at the single-cell level. We test ikarus on multiple single-cell datasets, showing that it achieves high sensitivity and specificity in multiple experimental contexts.

Advancing therapy for neuroblastoma

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

TH-MYCN tumors, but not tumor-derived cell lines, are adrenergic lineage, GD2+, and responsive to anti-GD2 antibody therapy

Neuroblastoma is a commonly lethal solid tumor of childhood and intensive chemoradiotherapy treatment cures ~50% of children with high-risk disease. The addition of immunotherapy using dinutuximab, a monoclonal antibody directed against the GD2 disialoganglioside expressed on neuroblasts, improves survival when incorporated into front-line therapy and shows robust activity in regressing relapsed disease when combined with chemotherapy. Still, many children succumb to neuroblastoma progression despite receiving dinutuximab-based immunotherapy, and efforts to counteract the immune suppressive signals responsible are warranted. Animal models of human cancers provide useful platforms to study immunotherapies. TH-MYCN transgenic mice are immunocompetent and develop neuroblastomas at autochthonous sites due to enforced MYCN expression in developing neural crest tissues. However, GD2-directed immunotherapy in this model has been underutilized due to the prevailing notion that TH-MYCN neuroblasts express insufficient GD2 to be targeted. We demonstrate that neuroblasts in TH-MYCN-driven tumors express GD2 at levels comparable to human neuroblastomas but rapidly lose GD2 expression when explanted ex vivo to establish tumor cell lines. This occurs in association with a transition from an adrenergic to mesenchymal differentiation state. Importantly, not only is GD2 expression retained on tumors in situ, treatment with a murine anti-GD2 antibody, 14G2a, markedly extends survival in such mice, including durable complete responses. Tumors in 14G2a-treated mice have fewer macrophage and myeloid-derived suppressor cells in their tumor microenvironment. Our findings support the utility of this model to inform immunotherapy approaches for neuroblastoma and potential opportunities to investigate drivers of adrenergic to mesenchymal fate decisions.

Intrinsic transcriptional heterogeneity in neuroblastoma guides mechanistic and therapeutic insights

Cell state is controlled by master transcription factors (mTFs) that determine the cellular gene expression program. Cancer cells acquire dysregulated gene expression programs by mutational and non-mutational processes. Intratumoral heterogeneity can result from cells displaying distinct mTF-regulated cell states, which co-exist within the tumor. One archetypal tumor associated with transcriptionally regulated heterogeneity is high-risk neuroblastoma (NB). Patients with NB have poor overall survival despite intensive therapies, and relapsed patients are commonly refractory to treatment. The cellular populations that comprise NB are marked by different cohorts of mTFs and differential sensitivity to conventional therapies. Recent studies have highlighted mechanisms by which NB cells dynamically shift the cell state with treatment, revealing new opportunities to control the cellular response to treatment by manipulating cell-state-defining transcriptional programs. Here, we review recent advances in understanding transcriptionally defined cancer heterogeneity. We offer challenges to the field to encourage translation of basic science into clinical benefit.

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.

Single-cell transcriptomics reveals a distinct developmental state of KMT2A-rearranged infant B-cell acute lymphoblastic leukemia

KMT2A-rearranged infant ALL is an aggressive childhood leukemia with poor prognosis. Here, we investigated the developmental state of KMT2A-rearranged infant B-cell acute lymphoblastic leukemia (B-ALL) using bulk messenger RNA (mRNA) meta-analysis and examination of single lymphoblast transcriptomes against a developing bone marrow reference. KMT2A-rearranged infant B-ALL was uniquely dominated by an early lymphocyte precursor (ELP) state, whereas less adverse NUTM1-rearranged infant ALL demonstrated signals of later developing B cells, in line with most other childhood B-ALLs. We compared infant lymphoblasts with ELP cells and revealed that the cancer harbored hybrid myeloid-lymphoid features, including nonphysiological antigen combinations potentially targetable to achieve cancer specificity. We validated surface coexpression of exemplar combinations by flow cytometry. Through analysis of shared mutations in separate leukemias from a child with infant KMT2A-rearranged B-ALL relapsing as AML, we established that KMT2A rearrangement occurred in very early development, before hematopoietic specification, emphasizing that cell of origin cannot be inferred from the transcriptional state.

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.

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.