A mechanistic classification of clinical phenotypes in neuroblastoma

NBAtlas: A harmonized single-cell transcriptomic reference atlas of human neuroblastoma tumors

Neuroblastoma, a rare embryonic tumor arising from neural crest development, is responsible for 15% of pediatric cancer-related deaths. Recently, several single-cell transcriptome studies were performed on neuroblastoma patient samples to investigate the cell of origin and tumor heterogeneity. However, these individual studies involved a small number of tumors and cells, limiting the conclusions that could be drawn. To overcome this limitation, we integrated seven single-cell or single-nucleus datasets into a harmonized cell atlas covering 362,991 cells across 61 patients. We use this atlas to decipher the transcriptional landscape of neuroblastoma at single-cell resolution, revealing associations between transcriptomic profiles and clinical outcomes within the tumor compartment. In addition, we characterize the complex immune-cell landscape and uncover considerable heterogeneity among tumor-associated macrophages. Finally, we showcase the utility of our atlas as a resource by expanding it with additional data and using it as a reference for data-driven cell-type annotation.

R-loops' m6A modification and its roles in cancers

R-loops are three-stranded nucleic acid structures composed of an RNA-DNA hybrid and a displaced DNA strand. They are widespread and play crucial roles in regulating gene expression, DNA replication, and DNA and histone modifications. However, their regulatory mechanisms remain unclear. As R-loop detection technology advances, changes in R-loop levels have been observed in cancer models, often associated with transcription-replication conflicts and genomic instability. N6-methyladenosine (m6A) is an RNA epigenetic modification that regulates gene expression by affecting RNA localization, splicing, translation, and degradation. Upon reviewing the literature, we found that R-loops with m6A modifications are implicated in tumor development and progression. This article summarizes the molecular mechanisms and detection methods of R-loops and m6A modifications in gene regulation, and reviews recent research on m6A-modified R-loops in oncology. Our goal is to provide new insights into the origins of genomic instability in cancer and potential strategies for targeted therapy.

TERT Expression and Clinical Outcome in Pulmonary Carcinoids

PURPOSE

The clinical course of pulmonary carcinoids ranges from indolent to fatal disease, suggesting that specific molecular alterations drive progression toward the fully malignant state. A similar spectrum of clinical phenotypes occurs in pediatric neuroblastoma, in which activation of telomerase reverse transcriptase (TERT) is decisive in determining the course of disease. We therefore investigated whether TERT expression defines the clinical fate of patients with pulmonary carcinoid.

METHODS

TERT expression was examined by RNA sequencing in a test cohort and a validation cohort of pulmonary carcinoids (n = 88 and n = 105, respectively). A natural TERT expression cutoff was determined in the test cohort on the basis of the distribution of TERT expression, and its prognostic value was assessed by Kaplan-Meier survival estimates and multivariable analyses. Telomerase activity was validated by telomere repeat amplification protocol assay.

RESULTS

Similar to neuroblastoma, TERT expression exhibited a bimodal distribution in pulmonary carcinoids, separating tumors into TERT-high and TERT-low subgroups. A natural TERT cutoff discriminated unfavorable from favorable clinical courses with high accuracy both in the test cohort (5-year overall survival [OS], 0.547 ± 0.132 v 1.0; P < .001) and the validation cohort (5-year OS, 0.788 ± 0.063 v 0.913 ± 0.048; P < .001). In line with these findings, telomerase activity was largely absent in TERT-low tumors, whereas it was readily detectable in TERT-high carcinoids. In multivariable analysis considering TERT expression, histology (typical v atypical carcinoid), and stage (≤IIA v ≥IIB), high TERT expression was an independent prognostic marker for poor survival, with a hazard ratio of 5.243 (95% CI, 1.943 to 14.148; P = .001).

CONCLUSION

Our data demonstrate that high TERT expression defines clinically aggressive pulmonary carcinoids with fatal outcome, similar to neuroblastoma, indicating that activation of TERT may be a defining feature of lethal cancers.

Pharmacologically Targeting Ferroptosis and Cuproptosis in Neuroblastoma

Neuroblastoma is a deadly pediatric cancer that originates from the neural crest and frequently develops in the abdomen or adrenal gland. Although multiple approaches, including chemotherapy, radiotherapy, targeted therapy, and immunotherapy, are recommended for treating neuroblastoma, the tumor will eventually develop resistance, leading to treatment failure and cancer relapse. Therefore, a firm understanding of the molecular mechanisms underlying therapeutic resistance is vital for the development of new effective therapies. Recent research suggests that cancer-specific modifications to multiple subtypes of nonapoptotic regulated cell death (RCD), such as ferroptosis and cuproptosis, contribute to therapeutic resistance in neuroblastoma. Targeting these specific types of RCD may be viable novel targets for future drug discovery in the treatment of neuroblastoma. In this review, we summarize the core mechanisms by which the inability to properly execute ferroptosis and cuproptosis can enhance the pathogenesis of neuroblastoma. Therefore, we focus on emerging therapeutic compounds that can induce ferroptosis or cuproptosis, delineating their beneficial pharmacodynamic effects in neuroblastoma treatment. Cumulatively, we suggest that the pharmacological stimulation of ferroptosis and ferroptosis may be a novel and therapeutically viable strategy to target neuroblastoma.

Whole-Exome Sequencing Reveals Novel Candidate Driver Mutations and Potential Druggable Mutations in Patients with High-Risk Neuroblastoma

Neuroblastoma is the most prevalent solid tumor in early childhood, with a 5-year overall survival rate of 40-60% in high-risk cases. Therefore, the identification of novel biomarkers for the diagnosis, prognosis, and therapy of neuroblastoma is crucial for improving the clinical outcomes of these patients. In this study, we conducted the whole-exome sequencing of 48 freshly frozen tumor samples obtained from the Biobank. Somatic variants were identified and selected using a bioinformatics analysis pipeline. The mutational signatures were determined using the Mutalisk online tool. Cancer driver genes and druggable mutations were predicted using the Cancer Genome Interpreter. The most common mutational signature was single base substitution 5. MUC4, MUC16, and FLG were identified as the most frequently mutated genes. Using the Cancer Genome Interpreter, we identified five recurrent cancer driver mutations spanning MUC16, MUC4, ALK, and CTNND1, with the latter being novel and containing a missense mutation, R439C. We also identified 11 putative actionable mutations including NF1 Q1798*, Q2616*, and S636X, ALK F1174L and R1275Q, SETD2 P10L and Q1829E, BRCA1 R612S, NOTCH1 D1670V, ATR S1372L, and FGFR1 N577K. Our findings provide a comprehensive overview of the novel information relevant to the underlying molecular pathogenesis and therapeutic targets of neuroblastoma.

KLF7 promotes neuroblastoma differentiation through the GTPase signaling pathway by upregulating neuroblast differentiation-associated protein AHNAKs and glycerophosphodiesterase GDPD5

The arrest of neural crest-derived sympathoadrenal neuroblast differentiation contributes to neuroblastoma formation, and overriding this blocked differentiation is a clear strategy for treating high-risk neuroblastoma. A better understanding of neuroblast or neuroblastoma differentiation is essential for developing new therapeutic approaches. It has been proposed that Krueppel-like factor 7 (KLF7) is a neuroblastoma super-enhancer-associated transcription factor gene. Moreover, KLF7 was found to be intensely active in postmitotic neuroblasts of the developing nervous system during embryogenesis. However, the role of KLF7 in the differentiation of neuroblast or neuroblastoma is unknown. Here, we find a strong association between high KLF7 expression and favorable clinical outcomes in neuroblastoma. KLF7 induces differentiation of neuroblastoma cells independently of the retinoic acid (RA) pathway and acts cooperatively with RA to induce neuroblastoma differentiation. KLF7 alters the GTPase activity and multiple differentiation-related genes by binding directly to the promoters of neuroblast differentiation-associated protein (AHNAK and AHNAK2) and glycerophosphodiester phosphodiesterase domain-containing protein 5 (GDPD5) and regulating their expression. Furthermore, we also observe that silencing KLF7 in neuroblastoma cells promotes the adrenergic-to-mesenchymal transition accompanied by changes in enhancer-mediated gene expression. Our results reveal that KLF7 is an inducer of neuroblast or neuroblastoma differentiation with prognostic significance and potential therapeutic value.

Higher tumor mutational burden is associated with inferior outcomes among pediatric patients with neuroblastoma

INTRODUCTION

Neuroblastoma is a pediatric malignancy with heterogeneous clinical outcomes. Our aim was to identify prognostic genetic markers for patients with neuroblastoma, who were treated with the Taiwan Pediatric Oncology Group (TPOG) neuroblastoma N2002 protocol, to improve risk stratification and inform treatment.

METHODS

Our analysis was based on 53 primary neuroblastoma specimens, diagnosed pre-chemotherapy, and 11 paired tumor relapse specimens. Deep sequencing of 113 target genes was performed using a custom panel. Multiplex ligation-dependent probe amplification was performed to identify clinical outcomes related to copy-number variations.

RESULTS

We identified 128 variations associated with survival, with the number of variations being higher in the relapse than that in the diagnostic specimen (p = .03). The risk of event and mortality was higher among patients with a tumor mutational burden ≥10 than that in patients with a lower burden (p < .0001). Multivariate analysis identified tumor mutational burden, MYCN amplification, and chromosome 3p deletion as significant prognostic factors, independent of age at diagnosis, sex, and tumor stage. The 5-year event-free survival and overall survival rate was lower among patients with high tumor burden than in patients with low tumor burden. Furthermore, there was no survival of patients with an ALK F1147L variation at 5 years after diagnosis.

CONCLUSIONS

Genome sequencing to determine the tumor mutational burden and ALK variations can improve the risk classification of neuroblastoma and inform treatment.

Altered methylation of imprinted genes in neuroblastoma: implications for prognostic refinement

BACKGROUND

Neuroblastoma (NB) is a complex disease, and the current understanding of NB biology is limited. Deregulation in genomic imprinting is a common event in malignancy. Since imprinted genes play crucial roles in early fetal growth and development, their role in NB pathogenesis could be suggested.

METHODS

We examined alterations in DNA methylation patterns of 369 NB tumours at 49 imprinted differentially methylated regions (DMRs) and assessed its association with overall survival probabilities and selected clinical and genomic features of the tumours. In addition, an integrated analysis of DNA methylation and allele-specific copy number alterations (CNAs) was performed, to understand the correlation between the two molecular events.

RESULTS

Several imprinted regions with aberrant methylation patterns in NB were identified. Regions that underwent loss of methylation in > 30% of NB samples were DMRs annotated to the genes NDN, SNRPN, IGF2, MAGEL2 and HTR5A and regions with gain of methylation were NNAT, RB1 and GPR1. Methylation alterations at six of the 49 imprinted DMRs were statistically significantly associated with reduced overall survival: MIR886, RB1, NNAT/BLCAP, MAGEL2, MKRN3 and INPP5F. RB1, NNAT/BLCAP and MKRN3 were further able to stratify low-risk NB tumours i.e. tumours that lacked MYCN amplification and 11q deletion into risk groups. Methylation alterations at NNAT/BLCAP, MAGEL2 and MIR886 predicted risk independently of MYCN amplification or 11q deletion and age at diagnosis. Investigation of the allele-specific CNAs demonstrated that the imprinted regions that displayed most alterations in NB tumours harbor true epigenetic changes and are not result of the underlying CNAs.

CONCLUSIONS

Aberrant methylation in imprinted regions is frequently occurring in NB tumours and several of these regions have independent prognostic value. Thus, these could serve as potentially important clinical epigenetic markers to identify individuals with adverse prognosis. Incorporation of methylation status of these regions together with the established risk predictors may further refine the prognostication of NB patients.

Defining Neuroblastoma: from origin to precision medicine

Neuroblastoma (NB), an heterogenous pediatric tumor of the sympathetic nervous system, is the most common and deadly extracranial solid malignancy diagnosed in infants. Numerous efforts have been invested in understanding its origin and in development of novel curative targeted therapies. Here, we summarize the recent advances in the identification of the cell of origin and the genetic alterations occurring during development contributing to NB. We discuss current treatment regimens, present and future directions for identification of novel therapeutic metabolic targets, differentiation agents, as well as personalized combinatory therapies as potential approaches for improving survival and quality of life of children with NB.

Long-term follow-up of patients with intermediate-risk neuroblastoma treated with response- and biology-based therapy: A report from the Children's Oncology Group study ANBL0531

BACKGROUND

We previously reported excellent three-year overall survival (OS) for patients with newly diagnosed intermediate-risk neuroblastoma treated with a biology- and response-based algorithm on the Children's Oncology Group study ANBL0531. We now present the long-term follow-up results.

METHODS

All patients who met the age, stage, and tumor biology criteria for intermediate-risk neuroblastoma were eligible. Treatment was based on prognostic biomarkers and overall response. Event-free survival (EFS) and OS were estimated by the Kaplan-Meier method.

RESULTS

The 10-year EFS and OS for the entire study cohort (n = 404) were 82.0% (95% confidence interval (CI), 77.2%-86.9%) and 94.7% (95% CI, 91.8%-97.5%), respectively. International Neuroblastoma Staging System stage 4 patients (n = 133) had inferior OS compared with non-stage 4 patients (n = 271; 10-year OS: 90.8% [95% CI, 84.5%-97.0%] vs 96.6% [95% CI, 93.9%-99.4%], p = .02). Infants with stage 4 tumors with ≥1 unfavorable biological feature (n = 47) had inferior EFS compared with those with favorable biology (n = 61; 10-year EFS: 66.8% [95% CI, 50.4%-83.3%] vs 86.9% [95% CI, 76.0%-97.8%], p = .02); OS did not differ (10-year OS: 84.4% [95% CI, 71.8%-97.0%] vs 95.0% [95% CI, 87.7%-100.0%], p = .08). Inferior EFS but not OS was observed among patients with tumors with (n = 26) versus without (n = 314) 11q loss of heterozygosity (10-year EFS: 68.4% [95% CI, 44.5%-92.2%] vs 83.9% [95% CI, 78.7%-89.2%], p = .03; 10-year OS: 88.0% [95% CI, 72.0%-100.0%] vs 95.7% [95% CI, 92.8%-98.6%], p = .09).

CONCLUSIONS

The ANBL0531 trial treatment algorithm resulted in excellent long-term survival. More effective treatments are needed for subsets of patients with unfavorable biology tumors.

Single-cell RNA-seq reveals the transcriptional program underlying tumor progression and metastasis in neuroblastoma

Neuroblastoma (NB) is one of the most common childhood malignancies. Sixty percent of patients present with widely disseminated clinical signs at diagnosis and exhibit poor outcomes. However, the molecular mechanisms triggering NB metastasis remain largely uncharacterized. In this study, we generated a transcriptomic atlas of 15 447 NB cells from eight NB samples, including paired samples of primary tumors and bone marrow metastases. We used time-resolved analysis to chart the evolutionary trajectory of NB cells from the primary tumor to the metastases in the same patient and identified a common 'starter' subpopulation that initiates tumor development and metastasis. The 'starter' population exhibited high expression levels of multiple cell cycle-related genes, indicating the important role of cell cycle upregulation in NB tumor progression. In addition, our evolutionary trajectory analysis demonstrated the involvement of partial epithelial-to-mesenchymal transition (p-EMT) along the metastatic route from the primary site to the bone marrow. Our study provides insights into the program driving NB metastasis and presents a signature of metastasis-initiating cells as an independent prognostic indicator and potential therapeutic target to inhibit the initiation of NB metastasis.

Preclinical spheroid models identify BMX as a therapeutic target for metastatic MYCN nonamplified neuroblastoma

The development of targeted therapies offers new hope for patients affected by incurable cancer. However, multiple challenges persist, notably in controlling tumor cell plasticity in patients with refractory and metastatic illness. Neuroblastoma (NB) is an aggressive pediatric malignancy originating from defective differentiation of neural crest-derived progenitors with oncogenic activity due to genetic and epigenetic alterations and remains a clinical challenge for high-risk patients. To identify critical genes driving NB aggressiveness, we performed combined chromatin and transcriptome analyses on matched patient-derived xenografts (PDXs), spheroids, and differentiated adherent cultures derived from metastatic MYCN nonamplified tumors. Bone marrow kinase on chromosome X (BMX) was identified among the most differentially regulated genes in PDXs and spheroids versus adherent models. BMX expression correlated with high tumor stage and poor patient survival and was crucial to the maintenance of the self-renewal and tumorigenic potential of NB spheroids. Moreover, BMX expression positively correlated with the mesenchymal NB cell phenotype, previously associated with increased chemoresistance. Finally, BMX inhibitors readily reversed this cellular state, increased the sensitivity of NB spheroids toward chemotherapy, and partially reduced tumor growth in a preclinical NB model. Altogether, our study identifies BMX as a promising innovative therapeutic target for patients with high-risk MYCN nonamplified NB.

Antitumor Activity of a Pyrrolobenzodiazepine Antibody-Drug Conjugate Targeting LGR5 in Preclinical Models of Neuroblastoma

Neuroblastoma (NB) is a cancer of the peripheral nervous system found in children under 15 years of age. It is the most frequently diagnosed cancer during infancy, accounting for ~12% of all cancer-related deaths in children. Leucine-rich repeat-containing G-protein-coupled receptor 5 (LGR5) is a membrane receptor that is associated with the primary tumor formation and metastasis of cancers in the gastrointestinal system. Remarkably, high levels of LGR5 are found in NB tumor cells, and high LGR5 expression is strongly correlated with poor survival. Antibody-drug conjugates (ADCs) are monoclonal antibodies that are covalently linked to cell-killing cytotoxins to deliver the payloads into cancer cells. We generated an ADC with an anti-LGR5 antibody and pyrrolobenzodiazepine (PBD) dimer-based payload SG3199 using a chemoenzymatic conjugation method. The resulting anti-LGR5 ADC was able to inhibit the growth of NB cells expressing LGR5 with high potency and specificity. Importantly, the ADC was able to completely inhibit the growth of NB xenograft tumors in vivo at a clinically relevant dose for the PBD class of ADCs. The findings support the potential of targeting LGR5 using the PBD class of payload for the treatment of high-risk NBs.

Loss of Chromosome Y in Neuroblastoma Is Associated With High-Risk Disease, 11q-Deletion, and Telomere Maintenance

Neuroblastoma (NB) is a heterogeneous childhood cancer with a slightly higher incidence in boys than girls, with the reason for this gender disparity unknown. Given the growing evidence for the involvement of loss of the Y chromosome (LoY) in male diseases including cancer, we investigated Y chromosome status in NB. Male NB tumor samples from a Swedish cohort, analyzed using Cytoscan HD SNP-microarray, were selected. Seventy NB tumors were analyzed for aneuploidy of the Y chromosome, and these data were correlated with other genetic, biological, and clinical parameters. LoY was found in 21% of the male NB tumors and it was almost exclusively found in those with high-risk genomic profiles. Furthermore, LoY was associated with increased age at diagnosis and enriched in tumors with 11q-deletion and activated telomere maintenance mechanisms. In contrast, tumors with an MYCN-amplified genomic profile retained their Y chromosome. The understanding of LoY in cancer is limited, making it difficult to conclude whether LoY is a driving event in NB or function of increased genomic instability. Gene expression analysis of Y chromosome genes in male NB tumors showed low expression of certain genes correlating with worse overall survival. KDM5D, encoding a histone demethylase stands out as an interesting candidate for further studies. LoY has been shown to impact the epigenomic layer of autosomal loci in nonreproductive tissues, and KDM5D has been reported as downregulated and/or associated with poor survival in different malignancies. Further studies are needed to explore the mechanisms and functional consequences of LoY in NB.

Preclinical Therapeutic Efficacy of RAF/MEK/ERK and IGF1R/AKT/mTOR Inhibition in Neuroblastoma

Activating mutations in the RAS/MAPK pathway are observed in relapsed neuroblastoma. Preclinical studies indicate that these tumors have an increased sensitivity to inhibitors of the RAS/MAPK pathway, such as MEK inhibitors. MEK inhibitors do not induce durable responses as single agents, indicating a need to identify synergistic combinations of targeted agents to provide therapeutic benefit. We previously showed preclinical therapeutic synergy between a MEK inhibitor, trametinib, and a monoclonal antibody specific for IGF1R, ganitumab in RAS-mutated rhabdomyosarcoma. Neuroblastoma cells, like rhabdomyosarcoma cells, are sensitive to the inhibition of the RAS/MAPK and IGF1R/AKT/mTOR pathways. We hypothesized that the combination of trametinib and ganitumab would be effective in RAS-mutated neuroblastoma. In this study, trametinib and ganitumab synergistically suppressed neuroblastoma cell proliferation and induced apoptosis in cell culture. We also observed a delay in tumor initiation and prolongation of survival in heterotopic and orthotopic xenograft models treated with trametinib and ganitumab. However, the growth of both primary and metastatic tumors was observed in animals receiving the combination of trametinib and ganitumab. Therefore, more preclinical work is necessary before testing this combination in patients with relapsed or refractory RAS-mutated neuroblastoma.

SAP30, an oncogenic driver of progression, poor survival, and drug resistance in neuroblastoma

Neuroblastoma is the most devastating extracranial solid malignancy in children. Despite an intense treatment regimen, the prognosis for high-risk neuroblastoma patients remains poor, with less than 40% survival. So far, MYCN amplification status is considered the most prognostic factor but corresponds to only ∼25% of neuroblastoma patients. Therefore, it is essential to identify a better prognosis and therapy response marker in neuroblastoma patients. We applied robust bioinformatic data mining tools, such as weighted gene co-expression network analysis, cisTarget, and single-cell regulatory network inference and clustering on two neuroblastoma patient datasets. We found Sin3A-associated protein 30 (SAP30), a driver transcription factor positively associated with high-risk, progression, stage 4, and poor survival in neuroblastoma patient cohorts. Tumors of high-risk neuroblastoma patients and relapse-specific patient-derived xenografts showed higher SAP30 levels. The advanced pharmacogenomic analysis and CRISPR-Cas9 screens indicated that SAP30 essentiality is associated with cisplatin resistance and further showed higher levels in cisplatin-resistant patient-derived xenograft tumor cell lines. Silencing of SAP30 induced cell death in vitro and led to a reduced tumor burden and size in vivo. Altogether, these results indicate that SAP30 is a better prognostic and cisplatin-resistance marker and thus a potential drug target in high-risk neuroblastoma.

Differential RNA Expression Between Metastatic and Primary Neuroblastoma Cells

INTRODUCTION

Neuroblastoma (NB) is the most common extra-cranial malignancy in children. Poor survival in high-risk NB is attributed to recurrent metastatic disease. To better study metastatic disease, we used a novel mouse model to investigate differential gene expression between primary tumor cells and metastatic cells. We hypothesized that metastatic NB cells have a different gene expression profile from primary tumor cells and cultured cells.

METHODS

Using three human NB cell lines (NGP, CHLA255, and SH-SY5Y), orthotopic xenografts were established in immunodeficient nod/scid gamma mice via subcapsular renal injection. Mice were sacrificed and NB cells were isolated from the primary tumor and from sites of metastasis (bone marrow, liver). RNA sequencing, gene set analysis, and pathway analysis were performed to identify differentially expressed genes and molecular pathways in the metastatic cells compared to primary tumor cells.

RESULTS

There were 266 differentially expressed genes in metastatic tumor cells (bone marrow and liver combined) compared to primary tumor cells. The top upregulated gene was KCNK1 and the top downregulated genes were PDE7B and NEBL. Top upregulated pathways in the metastatic cells were involved in ion transport, cell signaling, and cell proliferation. Top downregulated pathways were involved in DNA synthesis, transcription, and cellular metabolism.

CONCLUSIONS

In metastatic NB cells, our study identified the upregulation of biologic processes involved in cell cycle regulation, cell proliferation, migration, and invasion. Ongoing studies aim to validate downstream translation of these genomic alterations, as well as target these pathways to more effectively suppress and inhibit recurrent metastatic disease in NB.

Proteomic insights into paediatric cancer: Unravelling molecular signatures and therapeutic opportunities

Survival rates in some paediatric cancers have improved greatly over recent decades, in part due to the identification of diagnostic, prognostic and predictive molecular signatures, and the development of risk-directed therapies. However, other paediatric cancers have proved difficult to treat, and there is an urgent need to identify novel biomarkers that reveal therapeutic opportunities. The proteome is the total set of expressed proteins present in a cell or tissue at a point in time, and is vastly more dynamic than the genome. Proteomics holds significant promise for cancer research, as proteins are ultimately responsible for cellular phenotype and are the target of most anticancer drugs. Here, we review the discoveries, opportunities and challenges of proteomic analyses in paediatric cancer, with a focus on mass spectrometry (MS)-based approaches. Accelerating incorporation of proteomics into paediatric precision medicine has the potential to improve survival and quality of life for children with cancer.

Developing targeted therapies for neuroblastoma by dissecting the effects of metabolic reprogramming on tumor microenvironments and progression

Rationale: Synergic reprogramming of metabolic dominates neuroblastoma (NB) progression. It is of great clinical implications to develop an individualized risk prognostication approach with stratification-guided therapeutic options for NB based on elucidating molecular mechanisms of metabolic reprogramming. Methods: With a machine learning-based multi-step program, the synergic mechanisms of metabolic reprogramming-driven malignant progression of NB were elucidated at single-cell and metabolite flux dimensions. Subsequently, a promising metabolic reprogramming-associated prognostic signature (MPS) and individualized therapeutic approaches based on MPS-stratification were developed and further validated independently using pre-clinical models. Results: MPS-identified MPS-I NB showed significantly higher activity of metabolic reprogramming than MPS-II counterparts. MPS demonstrated improved accuracy compared to current clinical characteristics [AUC: 0.915 vs. 0.657 (MYCN), 0.713 (INSS-stage), and 0.808 (INRG-stratification)] in predicting prognosis. AZD7762 and etoposide were identified as potent therapeutics against MPS-I and II NB, respectively. Subsequent biological tests revealed AZD7762 substantially inhibited growth, migration, and invasion of MPS-I NB cells, more effectively than that of MPS-II cells. Conversely, etoposide had better therapeutic effects on MPS-II NB cells. More encouragingly, AZD7762 and etoposide significantly inhibited in-vivo subcutaneous tumorigenesis, proliferation, and pulmonary metastasis in MPS-I and MPS-II samples, respectively; thereby prolonging survival of tumor-bearing mice. Mechanistically, AZD7762 and etoposide-induced apoptosis of the MPS-I and MPS-II cells, respectively, through mitochondria-dependent pathways; and MPS-I NB resisted etoposide-induced apoptosis by addiction of glutamate metabolism and acetyl coenzyme A. MPS-I NB progression was fueled by multiple metabolic reprogramming-driven factors including multidrug resistance, immunosuppressive and tumor-promoting inflammatory microenvironments. Immunologically, MPS-I NB suppressed immune cells via MIF and THBS signaling pathways. Metabolically, the malignant proliferation of MPS-I NB cells was remarkably supported by reprogrammed glutamate metabolism, tricarboxylic acid cycle, urea cycle, etc. Furthermore, MPS-I NB cells manifested a distinct tumor-promoting developmental lineage and self-communication patterns, as evidenced by enhanced oncogenic signaling pathways activated with development and self-communications. Conclusions: This study provides deep insights into the molecular mechanisms underlying metabolic reprogramming-mediated malignant progression of NB. It also sheds light on developing targeted medications guided by the novel precise risk prognostication approaches, which could contribute to a significantly improved therapeutic strategy for NB.

EGCG Disrupts the LIN28B/Let-7 Interaction and Reduces Neuroblastoma Aggressiveness

Neuroblastoma (NB) is the most commonly diagnosed extracranial solid tumor in children, accounting for 15% of all childhood cancer deaths. Although the 5-year survival rate of patients with a high-risk disease has increased in recent decades, NB remains a challenge in pediatric oncology, and the identification of novel potential therapeutic targets and agents is an urgent clinical need. The RNA-binding protein LIN28B has been identified as an oncogene in NB and is associated with a poor prognosis. Given that LIN28B acts by negatively regulating the biogenesis of the tumor suppressor let-7 miRNAs, we reasoned that selective interference with the LIN28B/let-7 miRNA interaction would increase let-7 miRNA levels, ultimately leading to reduced NB aggressiveness. Here, we selected (-)-epigallocatechin 3-gallate (EGCG) out of 4959 molecules screened as the molecule with the best inhibitory activity on LIN28B/let-7 miRNA interaction and showed that treatment with PLC/PLGA-PEG nanoparticles containing EGCG (EGCG-NPs) led to an increase in mature let-7 miRNAs and a consequent inhibition of NB cell growth. In addition, EGCG-NP pretreatment reduced the tumorigenic potential of NB cells in vivo. These experiments suggest that the LIN28B/let-7 miRNA axis is a good therapeutic target in NB and that EGCG, which can interfere with this interaction, deserves further preclinical evaluation.

Targeting NRAS via miR-1304-5p or farnesyltransferase inhibition confers sensitivity to ALK inhibitors in ALK-mutant neuroblastoma

Targeting Anaplastic lymphoma kinase (ALK) is a promising therapeutic strategy for aberrant ALK-expressing malignancies including neuroblastoma, but resistance to ALK tyrosine kinase inhibitors (ALK TKI) is a distinct possibility necessitating drug combination therapeutic approaches. Using high-throughput, genome-wide CRISPR-Cas9 knockout screens, we identify miR-1304-5p loss as a desensitizer to ALK TKIs in aberrant ALK-expressing neuroblastoma; inhibition of miR-1304-5p decreases, while mimics of this miRNA increase the sensitivity of neuroblastoma cells to ALK TKIs. We show that miR-1304-5p targets NRAS, decreasing cell viability via induction of apoptosis. It follows that the farnesyltransferase inhibitor (FTI) lonafarnib in addition to ALK TKIs act synergistically in neuroblastoma, inducing apoptosis in vitro. In particular, on combined treatment of neuroblastoma patient derived xenografts with an FTI and an ALK TKI complete regression of tumour growth is observed although tumours rapidly regrow on cessation of therapy. Overall, our data suggests that combined use of ALK TKIs and FTIs, constitutes a therapeutic approach to treat high risk neuroblastoma although prolonged therapy is likely required to prevent relapse.

ALKBH5 gene polymorphisms and risk of neuroblastoma in Chinese children from Jiangsu Province

Background

Neuroblastoma is one of the most common extracranial malignant solid tumors in children. AlkB homolog 5 (ALKBH5) is an RNA N6-methyladenosine (m6A) demethylase that plays a critical role in tumorigenesis and development. We assessed the association between single nucleotide polymorphisms (SNPs) in ALKBH5 and the risk of neuroblastoma in a case-control study including 402 patients and 473 non-cancer controls.

Methods

Genotyping was determined by the TaqMan method. The association between ALKBH5 polymorphisms (rs1378602 and rs8400) and the risk of neuroblastoma was evaluated using the odds ratio (OR) and 95% confidence interval (CI).

Results

We found no strong association of ALKBH5 rs1378602 and rs8400 with neuroblastoma risk. Further stratification analysis by age, sex, primary site, and clinical stage showed that the rs1378602 AG/AA genotype was associated with a lower risk of neuroblastoma in males (adjusted OR = 0.58, 95% CI = 0.35-0.97, p = 0.036) and children with retroperitoneal neuroblastoma (adjusted OR = 0.58, 95% CI = 0.34-0.98, p = 0.040).

Conclusions

ALKBH5 SNPs do not seem to be associated with neuroblastoma risk. More studies are required to confirm this negative result and reveal the relationship between gene polymorphisms of the m6A modifier ALKBH5 and neuroblastoma.

METTL3 drives telomere targeting of TERRA lncRNA through m6A-dependent R-loop formation: a therapeutic target for ALT-positive neuroblastoma

Telomerase-negative tumors maintain telomere length by alternative lengthening of telomeres (ALT), but the underlying mechanism behind ALT remains poorly understood. A proportion of aggressive neuroblastoma (NB), particularly relapsed tumors, are positive for ALT (ALT+), suggesting that a better dissection of the ALT mechanism could lead to novel therapeutic opportunities. TERRA, a long non-coding RNA (lncRNA) derived from telomere ends, localizes to telomeres in a R-loop-dependent manner and plays a crucial role in telomere maintenance. Here we present evidence that RNA modification at the N6 position of internal adenosine (m6A) in TERRA by the methyltransferase METTL3 is essential for telomere maintenance in ALT+ cells, and the loss of TERRA m6A/METTL3 results in telomere damage. We observed that m6A modification is abundant in R-loop enriched TERRA, and the m6A-mediated recruitment of hnRNPA2B1 to TERRA is critical for R-loop formation. Our findings suggest that m6A drives telomere targeting of TERRA via R-loops, and this m6A-mediated R-loop formation could be a widespread mechanism employed by other chromatin-interacting lncRNAs. Furthermore, treatment of ALT+ NB cells with a METTL3 inhibitor resulted in compromised telomere targeting of TERRA and accumulation of DNA damage at telomeres, indicating that METTL3 inhibition may represent a therapeutic approach for ALT+ NB.

Asymmetric Adrenals: Sexual Dimorphism of Adrenal Tumors

CONTEXT

Sexual dimorphism has direct consequences on the incidence and survival of cancer. Early and accurate diagnosis is crucial to improve prognosis.

OBJECTIVE

This work aimed to characterized the influence of sex and adrenal asymmetry on the emergence of adrenal tumors.

METHODS

We conducted a multicenter, observational study involving 8037 patients with adrenal tumors, including adrenocortical carcinoma (ACC), aldosterone-producing adenoma (APA), cortisol-secreting adrenocortical adenomas (CSAs), non-aldosterone-producing adrenal cortical adenoma (NAPACA), pheochromocytoma (PCC), and neuroblastoma (NB), and investigated tumor lateralization according to sex. Human adrenal tissues (n = 20) were analyzed with a multiomics approach that allows determination of gene expression, catecholamine, and steroid contents in a single sample. In addition, we performed a literature review of computed tomography and magnetic resonance imaging-based studies examining adrenal gland size.

RESULTS

ACC (n = 1858); CSA (n = 68), NAPACA (n = 2174), and PCC (n = 1824) were more common in females than in males (female-to-male ratio: 1.1:1-3.8:1), whereas NBs (n = 2320) and APAs (n = 228) were less prevalent in females (0.8:1). ACC, APA, CSA, NAPACA, and NB occurred more frequently in the left than in the right adrenal (left-to-right ratio: 1.1:1-1.8:1), whereas PCC arose more often in the right than in the left adrenal (0.8:1). In both sexes, the left adrenal was larger than the right adrenal; females have smaller adrenals than males.

CONCLUSION

Adrenal asymmetry in both sexes may be related to the pathogenesis of adrenal tumors and should be considered during the diagnosis of these tumors.

Telomere Maintenance Mechanisms in a Cohort of High-Risk Neuroblastoma Tumors and Its Relation to Genomic Variants in the TERT and ATRX Genes

Tumor cells are hallmarked by their capacity to undergo unlimited cell divisions, commonly accomplished either by mechanisms that activate TERT or through the alternative lengthening of telomeres pathway. Neuroblastoma is a heterogeneous pediatric cancer, and the aim of this study was to characterize telomere maintenance mechanisms in a high-risk neuroblastoma cohort. All tumor samples were profiled with SNP microarrays and, when material was available, subjected to whole genome sequencing (WGS). Telomere length was estimated from WGS data, samples were assayed for the ALT biomarker c-circles, and selected samples were subjected to methylation array analysis. Samples with ATRX aberration in this study were positive for c-circles, whereas samples with either MYCN amplification or TERT re-arrangement were negative for c-circles. Both ATRX aberrations and TERT re-arrangement were enriched in 11q-deleted samples. An association between older age at diagnosis and 1q-deletion was found in the ALT-positive group. TERT was frequently placed in juxtaposition to a previously established gene in neuroblastoma tumorigenesis or cancer in general. Given the importance of high-risk neuroblastoma, means for mitigating active telomere maintenance must be therapeutically explored.

Rapid detection of telomerase expression of neuroblastoma in paraffin-embedded tissue: combination of in situ hybridisation and quantitative PCR

Neuroblastoma is a heterogeneous paediatric malignant tumour. Telomere maintenance mechanism (TMM) by telomerase activation or alternative lengthening of telomeres (ALT) is a hallmark of high-risk neuroblastoma. However, the prior assays for telomerase, such as TERT expression by RNA sequencing or microarrays, may not be easy to perform in many histopathology laboratories in hospitals. The aims of this study are to assess the utility of ultrasensitive single-cell RNA in situ hybridisation (RNAscope), immunohistochemistry, and RT-qPCR on formalin-fixed, paraffin-embedded tumour samples as diagnostic tools for detecting TERT expression in neuroblastoma. In this study, we detected MYCN amplification in 22 of 222 cases (10%), TERT rearrangements in 18 of 220 cases (8%), and ALT activation in 39 of 222 cases (18%) using fluorescence in situ hybridisation (FISH). By RNA in situ hybridisation, 36 of 210 (17%) pretreatment neuroblastomas were found to have TERT overexpression, which was significantly associated with the high-risk group (33/78, 42%), TERT rearrangements (16/18, 89%), and MYCN amplification (13/22, 59%). None of the tumours with ALT showed TERT staining. In our study, 19 of the 55 MYCN non-amplified high-risk neuroblastomas displayed TERT mRNA expression, including 13 of the 14 TERT rearrangements, none of the 30 ALT-positive cases, and a significant proportion (6/11, 55%) that did not have the aforementioned genomic anomalies. RT-qPCR results correlated well with RNAscope levels (Spearman's rho=0.621, p<0.001, n=94). In conclusion, TERT RNA in situ hybridisation and RT-qPCR are suitable methods to evaluate TERT expression in neuroblastoma. The combination of detection of the genomic alterations and TERT mRNA expression is a powerful strategy for TMM activation detection, which can categorise neuroblastomas into multiple clinical subgroups for risk stratification in routine histopathology practice.

Risk factors and novel predictive models for metastatic neuroblastoma in children

INTRODUCTION

Neuroblastoma (NB) with distant metastasis (DM) is a high-risk condition with a poor prognosis. Early identify the risk and prognostic differences of DM in children, which is helpful for the development of clinical diagnosis and treatment.

METHODS

The study cohort included patients with NB in surveillance, epidemiological, and final outcome databases between 2010 and 2018. To identify the risk and prognostic factors for DM, both univariate and multivariate logistic and Cox regression analyses were conducted. In addition, we created and verified three online clinical prediction models. Finally, we assess the performance of the proposed predictive model.

RESULTS

Among the 1224 children with NB included in the study, 599 developed DM. Primary site is the most important factor affecting metastasis and prognosis. The training and validation groups of the diagnostic nomograms had area under curves (AUC) values of 0.872 and 0.824, respectively. In addition, in the training group, the AUC values at 12, 36, and 60 months were 0.68, 0.71, and 0.75 for the OS nomogram and 0.70, 0.72, and 0.75 for the CSS nomogram. In the validation group, the AUC values at 12, 36, and 60 months were 0.68, 0.72, and 0.70 for the OS nomogram and 0.67, 0.71, and 0.69 for the CSS nomogram, respectively. Calibration curve and decision curve analyses revealed good performance of the nomogram.

CONCLUSIONS

The nomogram developed in this study could appropriately predict DM and assess its prognosis in patients with NB.

Optimal combination of MYCN differential gene and cellular senescence gene predicts adverse outcomes in patients with neuroblastoma

Introduction

Neuroblastoma (NB) is a common extracranial tumor in children and is highly heterogeneous. The factors influencing the prognosis of NB are not simple.

Methods

To investigate the effect of cell senescence on the prognosis of NB and tumor immune microenvironment, 498 samples of NB patients and 307 cellular senescence-related genes were used to construct a prediction signature.

Results

A signature based on six optimal candidate genes (TP53, IL-7, PDGFRA, S100B, DLL3, and TP63) was successfully constructed and proved to have good prognostic ability. Through verification, the signature had more advantages than the gene expression level alone in evaluating prognosis was found. Further T cell phenotype analysis displayed that exhausted phenotype PD-1 and senescence-related phenotype CD244 were highly expressed in CD8+ T cell in MYCN-amplified group with higher risk-score.

Conclusion

A signature constructed the six MYCN-amplified differential genes and aging-related genes can be used to predict the prognosis of NB better than using each high-risk gene individually and to evaluate immunosuppressed and aging tumor microenvironment.

Translational practice of fluorescence in situ hybridisation to identify neuroblastic tumours with TERT rearrangements

Recently, telomerase reverse transcriptase (TERT) gene rearrangements have been identified in neuroblastoma (NB), the typical pathological type of neuroblastic tumours (NTs); however, the prevalence of TERT rearrangements in other types of NT remains unknown. This study aimed to develop a practical method for detecting TERT defects and to evaluate the clinical relevance of TERT rearrangements as a biomarker for NT prognosis. A TERT break-apart probe for fluorescence in situ hybridisation (FISH) was designed, optimised, and applied to assess the genomic status of TERT in Chinese children with NTs at the Beijing Children's Hospital from 2016 to 2019. Clinical, histological, and genetic characteristics of TERT-rearranged NTs were further addressed. Genomic TERT rearrangements could be effectively detected by FISH and were mutually exclusive with MYCN amplification. TERT rearrangements were identified in 6.0% (38/633) of NTs overall, but 12.4% (31/250) in high-risk patients. TERT rearrangements identified a subtype of aggressive NTs with the characteristics of Stage 3/4, high-risk category, over 18 months old, and presenting all histological subtypes of NB and ganglioneuroblastoma nodular. Moreover, TERT rearrangements were significantly associated with elevated TERT expression levels and decreased survival chances. Multivariable analysis confirmed that it was an independent prognostic marker for NTs. FISH is an easily applicable method for evaluating TERT defects, which define a subgroup of NTs with unfavourable prognosis. TERT rearrangements would contribute to characterising NT molecular signatures in clinical practice.

[Peripheral neuroblastic tumors in childhood]

Peripheral neuroblastic tumors represent the fourth-largest group of malignant tumors in childhood. The majority of these tumors are neuroblastomas, which can be classified into undifferentiated, poorly differentiated, and differentiating subtypes. In addition, peripheral neuroblastic tumors include ganglioneuroblastoma, a composite tumor composed of Schwannian cell stroma and neuroblasts as well as benign ganglioneuroma. In this overview, histopathological diagnostic criteria and grading systems, as well as common molecular alterations that are of prognostic and therapeutic significance, are discussed.

A multiscale orchestrated computational framework to reveal emergent phenomena in neuroblastoma

Neuroblastoma is a complex and aggressive type of cancer that affects children. Current treatments involve a combination of surgery, chemotherapy, radiotherapy, and stem cell transplantation. However, treatment outcomes vary due to the heterogeneous nature of the disease. Computational models have been used to analyse data, simulate biological processes, and predict disease progression and treatment outcomes. While continuum cancer models capture the overall behaviour of tumours, and agent-based models represent the complex behaviour of individual cells, multiscale models represent interactions at different organisational levels, providing a more comprehensive understanding of the system. In 2018, the PRIMAGE consortium was formed to build a cloud-based decision support system for neuroblastoma, including a multi-scale model for patient-specific simulations of disease progression. In this work we have developed this multi-scale model that includes data such as patient's tumour geometry, cellularity, vascularization, genetics and type of chemotherapy treatment, and integrated it into an online platform that runs the simulations on a high-performance computation cluster using Onedata and Kubernetes technologies. This infrastructure will allow clinicians to optimise treatment regimens and reduce the number of costly and time-consuming clinical trials. This manuscript outlines the challenging framework's model architecture, data workflow, hypothesis, and resources employed in its development.

Chromosome 1p36 candidate gene ZNF436 predicts the prognosis of neuroblastoma: a bioinformatic analysis

BACKGROUND

Genetic 1p deletion is reported in 30% of all neuroblastoma and is associated with the unfavorable prognosis of neuroblastoma. The expressions and prognosis of 1p candidate genes in neuroblastoma are unclear.

METHODS

Public neuroblastoma cohorts were obtained for secondary analysis. The prognosis of 1p candidate genes in neuroblastoma was determined using Kaplan-Meier and cox regression analysis. The prediction of the nomogram model was determined using timeROC.

RESULTS

First, we confirmed the bad prognosis of 1p deletion in neuroblastoma. Moreover, zinc finger protein 436 (ZNF436) located at 1p36 region was down-regulated in 1p deleted neuroblastoma and higher ZNF436 expression was associated with the longer event free survival and overall survival of neuroblastoma. The expression levels of ZNF436 were lower in neuroblastoma patients with MYCN amplification or age at diagnosis ≥ 18months, or with stage 4 neuroblastoma. ZNF436 had robust predictive values of MYCN amplification and overall survival of neuroblastoma. Furthermore, the prognostic significance of ZNF436 in neuroblastoma was independent of MYCN amplification and age of diagnosis. Combinations of ZNF436 with MYCN amplification or age of diagnosis achieved better prognosis. At last, we constructed a nomogram risk model based on age, MYCN amplification and ZNF436. The nomogram model could predict the overall survival of neuroblastoma with high specificity and sensitivity.

CONCLUSIONS

Chromosome 1p36 candidate gene ZNF436 was a prognostic maker of neuroblastoma.

Neuroblastoma Patients' Outcome and Chromosomal Instability

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

DUBing Primary Tumors of the Central Nervous System: Regulatory Roles of Deubiquitinases

The ubiquitin proteasome system (UPS) utilizes an orchestrated enzymatic cascade of E1, E2, and E3 ligases to add single or multiple ubiquitin-like molecules as post-translational modification (PTM) to proteins. Ubiquitination can alter protein functions and/or mark ubiquitinated proteins for proteasomal degradation but deubiquitinases (DUBs) can reverse protein ubiquitination. While the importance of DUBs as regulatory factors in the UPS is undisputed, many questions remain on DUB selectivity for protein targeting, their mechanism of action, and the impact of DUBs on the regulation of diverse biological processes. Furthermore, little is known about the expression and role of DUBs in tumors of the human central nervous system (CNS). In this comprehensive review, we have used publicly available transcriptional datasets to determine the gene expression profiles of 99 deubiquitinases (DUBs) from five major DUB families in seven primary pediatric and adult CNS tumor entities. Our analysis identified selected DUBs as potential new functional players and biomarkers with prognostic value in specific subtypes of primary CNS tumors. Collectively, our analysis highlights an emerging role for DUBs in regulating CNS tumor cell biology and offers a rationale for future therapeutic targeting of DUBs in CNS tumors.

MYCN Amplifications and Metabolic Rewiring in Neuroblastoma

Cancer is a disease caused by (epi)genomic and gene expression abnormalities and characterized by metabolic phenotypes that are substantially different from the normal phenotypes of the tissues of origin. Metabolic reprogramming is one of the key features of tumors, including those established in the human nervous system. In this work, we emphasize a well-known cancerous genomic alteration: the amplification of MYCN and its downstream effects in neuroblastoma phenotype evolution. Herein, we extend our previous computational biology investigations by conducting an integrative workflow applied to published genomics datasets and comprehensively assess the impact of MYCN amplification in the upregulation of metabolism-related transcription factor (TF)-encoding genes in neuroblastoma cells. The results obtained first emphasized overexpressed TFs, and subsequently those committed in metabolic cellular processes, as validated by gene ontology analyses (GOs) and literature curation. Several genes encoding for those TFs were investigated at the mechanistic and regulatory levels by conducting further omics-based computational biology assessments applied on published ChIP-seq datasets retrieved from MYCN-amplified- and MYCN-enforced-overexpression within in vivo systems of study. Hence, we approached the mechanistic interrelationship between amplified MYCN and overexpression of metabolism-related TFs in neuroblastoma and showed that many are direct targets of MYCN in an amplification-inducible fashion. These results illuminate how MYCN executes its regulatory underpinnings on metabolic processes in neuroblastoma.

Children's Oncology Group's 2023 blueprint for research: Neuroblastoma

Neuroblastoma is the most common extra-cranial solid tumor in children and is known for its clinical heterogeneity. A greater understanding of the biology of this disease has led to both improved risk stratification and new approaches to therapy. Outcomes for children with low and intermediate risk disease are excellent overall, and efforts to decrease therapy for such patients have been largely successful. Although survival has improved over time for patients with high-risk disease and treatments evaluated in the relapse setting are now being moved into earlier phases of treatment, much work remains to improve survival and decrease therapy-related toxicities. Studies of highly annotated biobanked samples continue to lead to important insights regarding neuroblastoma biology. Such studies, along with correlative biology studies incorporated into therapeutic trials, are expected to continue to provide insights that lead to new and more effective therapies. A focus on translational science is accompanied by an emphasis on new agent development, optimized risk stratification, and international collaboration to address questions relevant to molecularly defined subsets of patients. In addition, the COG Neuroblastoma Committee is committed to addressing the patient/family experience, mitigating late effects of therapy, and studying social determinants of health in patients with neuroblastoma.

Two opposing gene expression patterns within ATRX aberrant neuroblastoma

Neuroblastoma is the most common extracranial solid tumor in children. A subgroup of high-risk patients is characterized by aberrations in the chromatin remodeller ATRX that is encoded by 35 exons. In contrast to other pediatric cancer where ATRX point mutations are most frequent, multi-exon deletions (MEDs) are the most frequent type of ATRX aberrations in neuroblastoma. 75% of these MEDs are predicted to produce in-frame fusion proteins, suggesting a potential gain-of-function effect compared to nonsense mutations. For neuroblastoma there are only a few patient-derived ATRX aberrant models. Therefore, we created isogenic ATRX aberrant models using CRISPR-Cas9 in several neuroblastoma cell lines and one tumoroid and performed total RNA-sequencing on these and the patient-derived models. Gene set enrichment analysis (GSEA) showed decreased expression of genes related to both ribosome biogenesis and several metabolic processes in our isogenic ATRX exon 2-10 MED model systems, the patient-derived MED models and in tumor data containing two patients with an ATRX exon 2-10 MED. In sharp contrast, these same processes showed an increased expression in our isogenic ATRX knock-out and exon 2-13 MED models. Our validations confirmed a role of ATRX in the regulation of ribosome homeostasis. The two distinct molecular expression patterns within ATRX aberrant neuroblastomas that we identified imply that there might be a need for distinct treatment regimens.

A Comparative Study of Neuroendocrine Heterogeneity in Small Cell Lung Cancer and Neuroblastoma

Lineage plasticity has long been documented in both small cell lung cancer (SCLC) and neuroblastoma, two clinically distinct neuroendocrine (NE) cancers. In this study, we quantified the NE features of cancer as NE scores and performed a systematic comparison of SCLC and neuroblastoma. We found neuroblastoma and SCLC cell lines have highly similar molecular profiles and shared therapeutic sensitivity. In addition, NE heterogeneity was observed at both the inter- and intra-cell line levels. Surprisingly, we did not find a significant association between NE scores and overall survival in SCLC or neuroblastoma. We described many shared and unique NE score-associated features between SCLC and neuroblastoma, including dysregulation of Myc oncogenes, alterations in protein expression, metabolism, drug resistance, and selective gene dependencies.

IMPLICATIONS

Our work establishes a reference for molecular changes and vulnerabilities associated with NE to non-NE transdifferentiation through mutual validation of SCLC and neuroblastoma samples.

Effective methods for bulk RNA-seq deconvolution using scnRNA-seq transcriptomes

BACKGROUND

RNA profiling technologies at single-cell resolutions, including single-cell and single-nuclei RNA sequencing (scRNA-seq and snRNA-seq, scnRNA-seq for short), can help characterize the composition of tissues and reveal cells that influence key functions in both healthy and disease tissues. However, the use of these technologies is operationally challenging because of high costs and stringent sample-collection requirements. Computational deconvolution methods that infer the composition of bulk-profiled samples using scnRNA-seq-characterized cell types can broaden scnRNA-seq applications, but their effectiveness remains controversial.

RESULTS

We produced the first systematic evaluation of deconvolution methods on datasets with either known or scnRNA-seq-estimated compositions. Our analyses revealed biases that are common to scnRNA-seq 10X Genomics assays and illustrated the importance of accurate and properly controlled data preprocessing and method selection and optimization. Moreover, our results suggested that concurrent RNA-seq and scnRNA-seq profiles can help improve the accuracy of both scnRNA-seq preprocessing and the deconvolution methods that employ them. Indeed, our proposed method, Single-cell RNA Quantity Informed Deconvolution (SQUID), which combines RNA-seq transformation and dampened weighted least-squares deconvolution approaches, consistently outperformed other methods in predicting the composition of cell mixtures and tissue samples.

CONCLUSIONS

We showed that analysis of concurrent RNA-seq and scnRNA-seq profiles with SQUID can produce accurate cell-type abundance estimates and that this accuracy improvement was necessary for identifying outcomes-predictive cancer cell subclones in pediatric acute myeloid leukemia and neuroblastoma datasets. These results suggest that deconvolution accuracy improvements are vital to enabling its applications in the life sciences.

Identification of the molecular subtypes and construction of risk models in neuroblastoma

The heterogeneity of neuroblastoma directly affects the prognosis of patients. Individualization of patient treatment to improve prognosis is a clinical challenge at this stage and the aim of this study is to characterize different patient populations. To achieve this, immune-related cell cycle genes, identified in the GSE45547 dataset using WGCNA, were used to classify cases from multiple datasets (GSE45547, GSE49710, GSE73517, GES120559, E-MTAB-8248, and TARGET) into subgroups by consensus clustering. ESTIMATES, CIBERSORT and ssGSEA were used to assess the immune status of the patients. And a 7-gene risk model was constructed based on differentially expressed genes between subtypes using randomForestSRC and LASSO. Enrichment analysis was used to demonstrate the biological characteristics between different groups. Key genes were screened using randomForest to construct neural network and validated. Finally, drug sensitivity was assessed in the GSCA and CellMiner databases. We classified the 1811 patients into two subtypes based on immune-related cell cycle genes. The two subtypes (Cluster1 and Cluster2) exhibited distinct clinical features, immune levels, chromosomal instability and prognosis. The same significant differences were demonstrated between the high-risk and low-risk groups. Through our analysis, we identified neuroblastoma subtypes with unique characteristics and established risk models which will improve our understanding of neuroblastoma heterogeneity.

Defining the landscape of circular RNAs in neuroblastoma unveils a global suppressive function of MYCN

Circular RNAs (circRNAs) are a regulatory RNA class. While cancer-driving functions have been identified for single circRNAs, how they modulate gene expression in cancer is not well understood. We investigate circRNA expression in the pediatric malignancy, neuroblastoma, through deep whole-transcriptome sequencing in 104 primary neuroblastomas covering all risk groups. We demonstrate that MYCN amplification, which defines a subset of high-risk cases, causes globally suppressed circRNA biogenesis directly dependent on the DHX9 RNA helicase. We detect similar mechanisms in shaping circRNA expression in the pediatric cancer medulloblastoma implying a general MYCN effect. Comparisons to other cancers identify 25 circRNAs that are specifically upregulated in neuroblastoma, including circARID1A. Transcribed from the ARID1A tumor suppressor gene, circARID1A promotes cell growth and survival, mediated by direct interaction with the KHSRP RNA-binding protein. Our study highlights the importance of MYCN regulating circRNAs in cancer and identifies molecular mechanisms, which explain their contribution to neuroblastoma pathogenesis.

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.

Preclinical Models of Neuroblastoma-Current Status and Perspectives

Preclinical in vitro and in vivo models remain indispensable tools in cancer research. These classic models, including two- and three-dimensional cell culture techniques and animal models, are crucial for basic and translational studies. However, each model has its own limitations and typically does not fully recapitulate the course of the human disease. Therefore, there is an urgent need for the development of novel, advanced systems that can allow for efficient evaluation of the mechanisms underlying cancer development and progression, more accurately reflect the disease pathophysiology and complexity, and effectively inform therapeutic decisions for patients. Preclinical models are especially important for rare cancers, such as neuroblastoma, where the availability of patient-derived specimens that could be used for potential therapy evaluation and screening is limited. Neuroblastoma modeling is further complicated by the disease heterogeneity. In this review, we present the current status of preclinical models for neuroblastoma research, discuss their development and characteristics emphasizing strengths and limitations, and describe the necessity of the development of novel, more advanced and clinically relevant approaches.

Survival of Patients With Neuroblastoma After Assignment to Reduced Therapy Because of the 12- to 18-Month Change in Age Cutoff in Children's Oncology Group Risk Stratification

PURPOSE

In 2006, Children's Oncology Group (COG) reclassified subgroups of toddlers diagnosed with neuroblastoma from high-risk to intermediate-risk, when the age cutoff for high-risk assignment was raised from 365 days (12 months) to 547 days (18 months). The primary aim of this retrospective study was to determine if excellent outcome was maintained after assigned reduction of therapy.

PATIENTS AND METHODS

Children <3 years old at diagnosis, enrolled on a COG biology study from 1990 to 2018, were eligible (n = 9,189). Assigned therapy was reduced for two cohorts of interest on the basis of the age cutoff change: 365-546 days old with International Neuroblastoma Staging System (INSS) stage 4, MYCN not amplified (MYCN-NA), favorable International Neuroblastoma Pathology Classification (INPC), hyperdiploid tumors (12-18mo/Stage4/FavBiology), and 365-546 days old with INSS stage 3, MYCN-NA, and unfavorable INPC tumors (12-18mo/Stage3/MYCN-NA/Unfav). Log-rank tests compared event-free survival (EFS) and overall survival (OS) curves.

RESULTS

For 12-18mo/Stage4/FavBiology, 5-year EFS/OS (± SE) before (≤2006; n = 40) versus after (>2006; n = 55) assigned reduction in therapy was similar: 89% ± 5.1%/89% ± 5.1% versus 87% ± 4.6%/94% ± 3.2% (P = .7; P = .4, respectively). For 12-18mo/Stage3/MYCN-NA/Unfav, the 5-year EFS and OS were both 100%, before (n = 6) and after (n = 4) 2006. The 12-18mo/Stage4/FavBiology plus 12-18mo/Stage3/MYCN-NA/Unfav classified as high-risk ≤2006 had an EFS/OS of 91% ± 4.4%/91% ± 4.5% versus 38% ± 1.3%/43% ± 1.3% for all other high-risk patients <3 years old (P < .0001; P < .0001, respectively). The 12-18mo/Stage4/FavBiology plus 12-18mo/Stage3/MYCN-NA/Unfav classified as intermediate-risk >2006 had an EFS/OS of 88% ± 4.3%/95% ± 2.9% versus 88% ± 0.9%/95% ± 0.6% for all other intermediate-risk patients <3 years old (P = .87; P = .85, respectively).

CONCLUSION

Excellent outcome was maintained among subsets of toddlers with neuroblastoma assigned to reduced treatment after reclassification of risk group from high to intermediate on the basis of new age cutoffs. Importantly, as documented in prior trials, intermediate-risk therapy is not associated with the degree of acute toxicity and late effects commonly observed with high-risk regimens.

HIF and MYC signaling in adrenal neoplasms of the neural crest: implications for pediatrics

Pediatric neural crest-derived adrenal neoplasms include neuroblastoma and pheochromocytoma. Both entities are associated with a high degree of clinical heterogeneity, varying from spontaneous regression to malignant disease with poor outcome. Increased expression and stabilization of HIF2α appears to contribute to a more aggressive and undifferentiated phenotype in both adrenal neoplasms, whereas MYCN amplification is a valuable prognostic marker in neuroblastoma. The present review focuses on HIF- and MYC signaling in both neoplasms and discusses the interaction of associated pathways during neural crest and adrenal development as well as potential consequences on tumorigenesis. Emerging single-cell methods together with epigenetic and transcriptomic analyses provide further insights into the importance of a tight regulation of HIF and MYC signaling pathways during adrenal development and tumorigenesis. In this context, increased attention to HIF-MYC/MAX interactions may also provide new therapeutic options for these pediatric adrenal neoplasms.

Clonal evolution during metastatic spread in high-risk neuroblastoma

Patients with high-risk neuroblastoma generally present with widely metastatic disease and often relapse despite intensive therapy. As most studies to date focused on diagnosis-relapse pairs, our understanding of the genetic and clonal dynamics of metastatic spread and disease progression remain limited. Here, using genomic profiling of 470 sequential and spatially separated samples from 283 patients, we characterize subtype-specific genetic evolutionary trajectories from diagnosis through progression and end-stage metastatic disease. Clonal tracing timed disease initiation to embryogenesis. Continuous acquisition of structural variants at disease-defining loci (MYCN, TERT, MDM2-CDK4) followed by convergent evolution of mutations targeting shared pathways emerged as the predominant feature of progression. At diagnosis metastatic clones were already established at distant sites where they could stay dormant, only to cause relapses years later and spread via metastasis-to-metastasis and polyclonal seeding after therapy.

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.

MYCN Amplification Is Associated with Reduced Expression of Genes Encoding γ-Secretase Complex and NOTCH Signaling Components in Neuroblastoma

Amplification of the MYCN oncogene is found in ~20% of neuroblastoma (NB) cases and correlates with high-risk disease and poor prognosis. Despite the plethora of studies describing the role of MYCN in NB, the exact molecular mechanisms underlying MYCN's contribution to high-risk disease are not completely understood. Herein, we implemented an integrative approach combining publicly available RNA-Seq and MYCN ChIP-Seq datasets derived from human NB cell lines to define biological processes directly regulated by MYCN in NB. Our approach revealed that MYCN-amplified NB cell lines, when compared to non-MYCN-amplified cell lines, are characterized by reduced expression of genes involved in NOTCH receptor processing, axoneme assembly, and membrane protein proteolysis. More specifically, we found genes encoding members of the γ-secretase complex, which is known for its ability to liberate several intracellular signaling molecules from membrane-bound proteins such as NOTCH receptors, to be down-regulated in MYCN-amplified NB cell lines. Analysis of MYCN ChIP-Seq data revealed an enrichment of MYCN binding at the transcription start sites of genes encoding γ-secretase complex subunits. Notably, using publicly available gene expression data from NB primary tumors, we revealed that the expression of γ-secretase subunits encoding genes and other components of the NOTCH signaling pathway was also reduced in MYCN-amplified tumors and correlated with worse overall survival in NB patients. Genetic or pharmacological depletion of MYCN in NB cell lines induced the expression of γ-secretase genes and NOTCH-target genes. Chemical inhibition of γ-secretase activity dampened the expression of NOTCH-target genes upon MYCN depletion in NB cells. In conclusion, this study defines a set of MYCN-regulated pathways that are specific to MYCN-amplified NB tumors, and it suggests a novel role for MYCN in the suppression of genes of the γ-secretase complex, with an impact on the NOTCH-target gene expression in MYCN-amplified NB.

TAF1D promotes proliferation by transcriptionally activating G2/M phase-related genes in MYCN-amplified neuroblastoma

High-risk neuroblastoma (HR-NB) is an aggressive childhood cancer that responds poorly to currently available therapies and is associated with only about a 50% 5-year survival rate. MYCN amplification is a critical driver of these aggressive tumors, but so far there have not been any approved treatments to effectively treat HR-NB by targeting MYCN or its downstream effectors. Thus, the identification of novel molecular targets and therapeutic strategies to treat children diagnosed with HR-NB represents an urgent unmet medical need. Here, we conducted a targeted siRNA screening and identified TATA box-binding protein-associated factor RNA polymerase I subunit D, TAF1D, as a critical regulator of the cell cycle and proliferation in HR-NB cells. Analysis of three independent primary NB cohorts determined that high TAF1D expression correlated with MYCN-amplified, high-risk disease and poor clinical outcomes. TAF1D knockdown more robustly inhibited cell proliferation in MYCN-amplified NB cells compared with MYCN-non-amplified NB cells, as well as suppressed colony formation and inhibited tumor growth in a xenograft mouse model of MYCN-amplified NB. RNA-seq analysis revealed that TAF1D knockdown downregulates the expression of genes associated with the G2/M transition, including the master cell-cycle regulator, cell-cycle-dependent kinase 1 (CDK1), resulting in cell-cycle arrest at G2/M. Our findings demonstrate that TAF1D is a key oncogenic regulator of MYCN-amplified HR-NB and suggest that therapeutic targeting of TAF1D may be a viable strategy to treat HR-NB patients by blocking cell-cycle progression and the proliferation of tumor cells.

Whole-exome sequencing reveals genetic variants in low-risk and high-risk neuroblastoma

This study aimed to investigate the genetic aberrations in neuroblastoma (NB) by comparing high and low-risk NB patients by whole-exome sequencing (WES) and to reveal the heterogeneity and association between somatic variants and clinical features. Seven NB patients with available clinical data were included in the study (4 in the low-risk group and 3 in the high-risk group). WES was performed and somatic variants associated with NB genes in the COSMIC database were selected through bioinformatics pipeline analysis. Variants were determined using the Integrative Genomics Viewer (IGV). Some gene variations were found in both groups, including variations in oncogene and tumor suppressor genes. In general, candidate gene variations were associated with chromatin remodeling complexes, the RAS pathway, cell proliferation, and DNA repair mechanism. Some variations in CSF1R, MSH6, PTPN11, SOX9, RET, TSC1, and DNMT1 genes were detected only in high-risk patients, while EP300, TET2, MYCN, PRDM1, and ARID2 gene variations were detected only in low-risk patients. When high-risk gene variants were compared with the cBioportal cancer genomic database, two common gene variants (ARID1A and NCOR2) were identified. However, when low-risk gene variants were compared with the cBioportal cancer genomic database, no common genes were found. GO/KEGG enrichment analysis was performed to find relevant biological processes and molecular pathways related to gene variants, which will help to decipher the molecular mechanisms of NB tumorigenesis and the phenotypic differences between high-risk and low-risk patients.