The TERT Promoter is Polycomb-Repressed in Neuroblastoma Cells with Long Telomeres.

Acquiring a telomere maintenance mechanism is a hallmark of high-risk neuroblastoma and commonly occurs by expressing telomerase (TERT). Telomerase-negative neuroblastoma has long telomeres and utilizes the telomerase-independent alternative lengthening of telomeres (ALT) mechanism. Conversely, no discernable telomere maintenance mechanism is detected in a fraction of neuroblastoma with long telomeres. Here, we show, unlike most cancers, DNA of the TERT promoter is broadly hypomethylated in neuroblastoma. In telomerase-positive neuroblastoma cells, the hypomethylated DNA promoter is approximately 1.5 kb. The TERT locus shows active chromatin marks with low enrichment for the repressive mark, H3K27me3. MYCN, a commonly amplified oncogene in neuroblstoma, binds to the promoter and induces TERT expression. Strikingly, in neuroblastoma with long telomeres, the hypomethylated region spans the entire TERT locus, including multiple nearby genes with enrichment for the repressive H3K27me3 chromatin mark. Furthermore, subtelomeric regions showed enrichment of repressive chromatin marks in neuroblastomas with long telomeres relative to those with short telomeres. These repressive marks were even more evident at the genic loci, suggesting a telomere position effect (TPE). Inhibiting H3K27 methylation by three different EZH2 inhibitors induced the expression of TERT in cell lines with long telomeres and H3K27me3 marks in the promoter region. EZH2 inhibition facilitated MYCN binding to the TERT promoter in neuroblastoma cells with long telomeres. Taken together, these data suggest that epigenetic regulation of TERT expression differs in neuroblastoma depending on the telomere maintenance status, and H3K27 methylation is important in repressing TERT expression in neuroblastoma with long telomeres. SIGNIFICANCE: The epigenetic landscape of the TERT locus is unique in neuroblastoma. The DNA at the TERT locus, unlike other cancer cells and similar to normal cells, are hypomethylated in telomerase-positive neuroblastoma cells. The TERT locus is repressed by polycomb repressive complex-2 complex in neuroblastoma cells that have long telomeres and do not express TERT. Long telomeres in neuroblastoma cells are also associated with repressive chromatin states at the chromosomal termini, suggesting TPE.

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.

KAP1 stabilizes MYCN mRNA and promotes neuroblastoma tumorigenicity by protecting the RNA m6A reader YTHDC1 protein degradation.

BACKGROUND: Neuroblastoma (NB) patients with amplified MYCN often face a grim prognosis and are resistant to existing therapies, yet MYCN protein is considered undruggable. KAP1 (also named TRIM28) plays a crucial role in multiple biological activities. This study aimed to investigate the relationship between KAP1 and MYCN in NB. METHODS: Transcriptome analyses and luciferase reporter assay identified that KAP1 was a downstream target of MYCN. The effects of KAP1 on cancer cell proliferation and colony formation were explored using the loss-of-function assays in vitro and in vivo. RNA stability detection was used to examine the influence of KAP1 on MYCN expression. The mechanisms of KAP1 to maintain MYCN mRNA stabilization were mainly investigated by mass spectrum, immunoprecipitation, RIP-qPCR, and western blotting. In addition, a xenograft mouse model was used to reveal the antitumor effect of STM2457 on NB. RESULTS: Here we identified KAP1 as a critical regulator of MYCN mRNA stability by protecting the RNA N6-methyladenosine (m6A) reader YTHDC1 protein degradation. KAP1 was highly expressed in clinical MYCN-amplified NB and was upregulated by MYCN. Reciprocally, KAP1 knockdown reduced MYCN mRNA stability and inhibited MYCN-amplified NB progression. Mechanistically, KAP1 regulated the stability of MYCN mRNA in an m6A-dependent manner. KAP1 formed a complex with YTHDC1 and RNA m6A writer METTL3 to regulate m6A-modified MYCN mRNA stability. KAP1 depletion decreased YTHDC1 protein stability and promoted MYCN mRNA degradation. Inhibiting MYCN mRNA m6A modification synergized with chemotherapy to restrain tumor progression in MYCN-amplified NB. CONCLUSIONS: Our research demonstrates that KAP1, transcriptionally activated by MYCN, forms a complex with YTHDC1 and METTL3, which in turn maintain the stabilization of MYCN mRNA in an m6A-dependent manner. Targeting m6A modification by STM2457, a small-molecule inhibitor of METTL3, could downregulate MYCN expression and attenuate tumor proliferation. This finding provides a new alternative putative therapeutic strategy for MYCN-amplified NB.

The MYCN oncoprotein is an RNA-binding accessory factor of the nuclear exosome targeting complex.

The MYCN oncoprotein binds active promoters in a heterodimer with its partner protein MAX. MYCN also interacts with the nuclear exosome, a 3'-5' exoribonuclease complex, suggesting a function in RNA metabolism. Here, we show that MYCN forms stable high-molecular-weight complexes with the exosome and multiple RNA-binding proteins. MYCN binds RNA in vitro and in cells via a conserved sequence termed MYCBoxI. In cells, MYCN associates with thousands of intronic transcripts together with the ZCCHC8 subunit of the nuclear exosome targeting complex and enhances their processing. Perturbing exosome function results in global re-localization of MYCN from promoters to intronic RNAs. On chromatin, MYCN is then replaced by the MNT(MXD6) repressor protein, inhibiting MYCN-dependent transcription. RNA-binding-deficient alleles show that RNA-binding limits MYCN's ability to activate cell growth-related genes but is required for MYCN's ability to promote progression through S phase and enhance the stress resilience of neuroblastoma cells.

The KAT module of the SAGA complex maintains the oncogenic gene expression program in MYCN-amplified neuroblastoma.

Pediatric cancers are frequently driven by genomic alterations that result in aberrant transcription factor activity. Here, we used functional genomic screens to identify multiple genes within the transcriptional coactivator Spt-Ada-Gcn5-acetyltransferase (SAGA) complex as selective dependencies for MYCN-amplified neuroblastoma, a disease of dysregulated development driven by an aberrant oncogenic transcriptional program. We characterized the DNA recruitment sites of the SAGA complex in neuroblastoma and the consequences of loss of SAGA complex lysine acetyltransferase (KAT) activity on histone acetylation and gene expression. We demonstrate that loss of SAGA complex KAT activity is associated with reduced MYCN binding on chromatin, suppression of MYC/MYCN gene expression programs, and impaired cell cycle progression. Further, we showed that the SAGA complex is pharmacologically targetable in vitro and in vivo with a KAT2A/KAT2B proteolysis targeting chimeric. Our findings expand our understanding of the histone-modifying complexes that maintain the oncogenic transcriptional state in this disease and suggest therapeutic potential for inhibitors of SAGA KAT activity in MYCN-amplified neuroblastoma.

In vivo cisplatin-resistant neuroblastoma metastatic model reveals tumour necrosis factor receptor superfamily member 4 (TNFRSF4) as an independent prognostic factor of survival in neuroblastoma.

Neuroblastoma is the most common solid extracranial tumour in children. Despite major advances in available therapies, children with drug-resistant and/or recurrent neuroblastoma have a dismal outlook with 5-year survival rates of less than 20%. Therefore, tackling relapsed tumour biology by developing and characterising clinically relevant models is a priority in finding targetable vulnerability in neuroblastoma. Using matched cisplatin-sensitive KellyLuc and resistant KellyCis83Luc cell lines, we developed a cisplatin-resistant metastatic MYCN-amplified neuroblastoma model. The average number of metastases per mouse was significantly higher in the KellyCis83Luc group than in the KellyLuc group. The vast majority of sites were confirmed as having lymph node metastasis. Their stiffness characteristics of lymph node metastasis values were within the range reported for the patient samples. Targeted transcriptomic profiling of immuno-oncology genes identified tumour necrosis factor receptor superfamily member 4 (TNFRSF4) as a significantly dysregulated MYCN-independent gene. Importantly, differential TNFRSF4 expression was identified in tumour cells rather than lymphocytes. Low TNFRSF4 expression correlated with poor prognostic indicators in neuroblastoma, such as age at diagnosis, stage, and risk stratification and significantly associated with reduced probability of both event-free and overall survival in neuroblastoma. Therefore, TNFRSF4 Low expression is an independent prognostic factor of survival in neuroblastoma.

Lineage specific transcription factor waves reprogram neuroblastoma from self-renewal to differentiation.

Temporal regulation of super-enhancer (SE) driven transcription factors (TFs) underlies normal developmental programs. Neuroblastoma (NB) arises from an inability of sympathoadrenal progenitors to exit a self-renewal program and terminally differentiate. To identify SEs driving TF regulators, we use all-trans retinoic acid (ATRA) to induce NB growth arrest and differentiation. Time-course H3K27ac ChIP-seq and RNA-seq reveal ATRA coordinated SE waves. SEs that decrease with ATRA link to stem cell development (MYCN, GATA3, SOX11). CRISPR-Cas9 and siRNA verify SOX11 dependency, in vitro and in vivo. Silencing the SOX11 SE using dCAS9-KRAB decreases SOX11 mRNA and inhibits cell growth. Other TFs activate in sequential waves at 2, 4 and 8 days of ATRA treatment that regulate neural development (GATA2 and SOX4). Silencing the gained SOX4 SE using dCAS9-KRAB decreases SOX4 expression and attenuates ATRA-induced differentiation genes. Our study identifies oncogenic lineage drivers of NB self-renewal and TFs critical for implementing a differentiation program.

Development and validation of a novel nomogram for predicting overall survival patients with neuroblastoma.

PURPOSE: The aim of this study was to develop a nomogram specially for predicting overall survival (OS) for Chinese patients with neuroblastoma (NB). METHODS: Patients with pathologically confirmed NB who were newly diagnosed and received treatments at our hospital from October 2013 to October 2021 were retrospectively reviewed. The nomogram for OS were built based on Cox regression analysis. The validation of the prognostic model was evaluated by concordance index (C-index), calibration curves, and decision curve analyses (DCAs). RESULTS: A total of 254 patients with NB were included in this study. They were randomly divided into a training cohort (n = 178) and a validation cohort (n = 76) at a ratio of 7:3. Multivariate analyses revealed that prognostic variables significantly related to the OS were age at diagnosis, bone metastasis, hepatic metastasis, INSS stage, MYCN status and DNA ploidy. The nomogram was constructed based on above 6 factors. The C-index values of the nomogram for predicting 3-year and 5-year OS were 0.926 and 0.964, respectively. The calibration curves of the nomogram showed good consistency between nomogram prediction and actual survival. The DCAs showed great clinical usefulness of the nomograms. Furthermore, patients with low-risk identified by our nomogram had much higher OS than those with high-risk (p < 0.001). CONCLUSION: The nomogram we constructed exhibited good predictive performance and could be used to assist clinicians in their decision-making process.

Analysis of the relationship between MYCN gene and serum NSE and urinary VMA levels and neuroblastoma pathological features and prognosis.

The purpose of this study was to explore the relationship between the MYCN gene, serum neuron-specific enolase (NSE), urinary vanillylmandelic acid (VMA) levels, and neuroblastoma pathological features and prognosis. Ninety-four children with neuroblastoma treated in the hospital were selected to compare the differences in MYCN gene amplification, serum NSE, and urinary VMA levels in children with different clinicopathological features and prognoses. The proportion of children with MYCN gene copy number ≥10 in INSS stage 3-4 was higher than that of children with INSS stage 1-2 (P < 0.05); the proportion of children with MYCN gene copy number ≥10 in high-risk children in the COG risk stratification was higher than that of children with intermediate and low risk (P < 0.05); the serum NSE of children aged >12 months higher than that of children aged ≤12 months (P < 0.05); serum NSE of children with tumors >500 cm3 higher than that of children with tumors ≤500 cm3 (P < 0.05); serum NSE and urinary VMA of children with INSS staging of stages 3-4 were higher than that of children with stages 1 to 2 (P < 0.05); serum NSE and urinary VMA in children with lymph node metastasis were higher than that of children without lymph node metastasis (P < 0.05); serum NSE of children with MYCN gene copy number ≥10 was higher than that of children without lymph node metastasis (P < 0.05); the proportion of children with MYCN gene copy number ≥10 who died, and the percentages of serum NSE and urinary VMA were higher than those of the surviving children (P < 0.05). MYCN gene amplification and serum NSE and urinary VMA levels were related to the age, tumor size, INSS stage, COG stage, lymph node metastasis, and prognosis of the children with neuroblastoma.

The MYCN 5' UTR as a therapeutic target in neuroblastoma.

Tumor MYCN amplification is seen in high-risk neuroblastoma, yet direct targeting of this oncogenic transcription factor has been challenging. Here, we take advantage of the dependence of MYCN-amplified neuroblastoma cells on increased protein synthesis to inhibit the activity of eukaryotic translation initiation factor 4A1 (eIF4A1) using an amidino-rocaglate, CMLD012824. Consistent with the role of this RNA helicase in resolving structural barriers in 5' untranslated regions (UTRs), CMLD012824 increased eIF4A1 affinity for polypurine-rich 5' UTRs, including that of the MYCN and associated transcripts with critical roles in cell proliferation. CMLD012824-mediated clamping of eIF4A1 spanned the full lengths of mRNAs, while translational inhibition was mediated through 5' UTR binding in a cap-dependent and -independent manner. Finally, CMLD012824 led to growth inhibition in MYCN-amplified neuroblastoma models without generalized toxicity. Our studies highlight the key role of eIF4A1 in MYCN-amplified neuroblastoma and demonstrate the therapeutic potential of disrupting its function.

MYCN drives oncogenesis by cooperating with the histone methyltransferase G9a and the WDR5 adaptor to orchestrate global gene transcription.

MYCN activates canonical MYC targets involved in ribosome biogenesis, protein synthesis, and represses neuronal differentiation genes to drive oncogenesis in neuroblastoma (NB). How MYCN orchestrates global gene expression remains incompletely understood. Our study finds that MYCN binds promoters to up-regulate canonical MYC targets but binds to both enhancers and promoters to repress differentiation genes. MYCN binding also increases H3K4me3 and H3K27ac on canonical MYC target promoters and decreases H3K27ac on neuronal differentiation gene enhancers and promoters. WDR5 facilitates MYCN promoter binding to activate canonical MYC target genes, whereas MYCN recruits G9a to enhancers to repress neuronal differentiation genes. Targeting both MYCN's active and repressive transcriptional activities using both WDR5 and G9a inhibitors synergistically suppresses NB growth. We demonstrate that MYCN cooperates with WDR5 and G9a to orchestrate global gene transcription. The targeting of both these cofactors is a novel therapeutic strategy to indirectly target the oncogenic activity of MYCN.

Withaferin A alters the expression of microRNAs 146a-5p and 34a-5p and associated hub genes in MDA-MB-231 cells.

Triple-negative breast cancer (TNBC) is a highly metastatic subtype of breast cancer. Due to the absence of obvious therapeutic targets, microRNAs (miRNAs) provide possible hope to treat TNBC. Withaferin A (WA), a steroidal lactone, possesses potential anticancer activity with lesser side effects. The present study identifies hub genes (CDKN3, TRAF6, CCND1, JAK1, MET, AXIN2, JAG1, VEGFA, BRCA1, E2F3, WNT1, CDK6, KRAS, MYB, MYCN, TGFßR2, NOTCH1, SIRT1, MYCN, NOTCH2, WNT3A) from the list of predicted targets of the differentially expressed miRNAs (DEMs) in WA-treated MDA-MB-231 cells using in silico protein-protein interaction network analysis. CCND1, CDK6, and TRAF6 hub genes were predicted as targets of miR-34a-5p and miR-146a-5p, respectively. The study found the lower expression of miR-34a-5p and miR-146a-5p in MDA-MB-231 cells, and further, it was observed that WA treatment effectively restored the lost expression of miR-34a-5p and miR-146a-5p in MDA-MB-231 cells. An anti-correlation expression pattern was found among the miR-34a-5p and miR-146a-5p and the respective target hub genes in WA-treated TNBC cells. In conclusion, WA might exert anti-cancer effect in TNBC cells by inducing miR-34a-5p and miR-146a-5p expressions and decreasing CCND1, CDK6, and TARF6 target hub genes in TNBC cells.

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

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

Ependymomas of the spinal region in adults: Clinical and pathological features and MYCN expression levels in spinal ependymomas and myxopapillary ependymomas.

BACKGROUND: Ependymomas (EPNs) of the spinal region are a heterogeneous group of tumors that account for 17.6 % in adults. Four types have been recognized: subependymoma, spinal ependymoma (Sp-EPN), myxopapillary ependymoma (MPE), and Sp-EPN-MYCN amplified, each with distinct histopathological and molecular features. METHODS: This study investigated the clinical and pathological characteristics and MYCN expression levels of 35 Sp-EPN and MPE cases diagnosed at a tertiary university hospital over a decade-long period. RESULTS: Twenty-five cases were Sp-EPN and 10 cases were MPE, and were graded as WHO grade 2, except for 1 Sp-EPN case with grade 3 features. The most common symptoms were lower back pain and difficulty in walking. Radiology showed different tumor sizes and locations along the spinal cord, with MPEs exclusively in the lumbosacral region. Surgery was the main treatment, and gross total resection was achieved in all cases except for one. Immunohistochemistry showed low Ki-67 proliferation indices in all cases, and no MYCN expression. During follow-up, 3 (8.6 %) cases recurred and/or metastasized and 5 cases (14.3 %) died. No significant difference was found in disease-free survival or overall survival between Sp-EPN and MPE cases. However, 3 cases with grade 2 histology demonstrated recurrence and/or metastasis, despite the lack of MYCN expression. CONCLUSION: Our results underscore the multifactorial nature of tumor aggressiveness in EPNs of the spinal region. This study enhances our knowledge of the clinical and pathological features of Sp-EPNs and MPEs and highlights the need for better diagnostic and prognostic markers in these rare tumors.

SESN1 functions as a new tumor suppressor gene via Toll-like receptor signaling pathway in neuroblastoma.

AIMS: Neuroblastoma (NB) is the most common extracranial solid tumor in children, with a 5-year survival rate of <50% in high-risk patients. MYCN amplification is an important factor that influences the survival rate of high-risk patients. Our results indicated MYCN regulates the expression of SESN1. Therefore, this study aimed to investigate the role and mechanisms of SESN1 in NB. METHODS: siRNAs or overexpression plasmids were used to change MYCN, SESN1, or MyD88's expression. The role of SESN1 in NB cell proliferation, migration, and invasion was elucidated. Xenograft mice models were built to evaluate SESN1's effect in vivo. The correlation between SESN1 expression and clinicopathological data of patients with NB was analyzed. RNA-Seq was done to explore SESN1's downstream targets. RESULTS: SESN1 was regulated by MYCN in NB cells. Knockdown SESN1 promoted NB cell proliferation, cell migration, and cell invasion, and overexpressing SESN1 had opposite functions. Knockdown SESN1 promoted tumor growth and shortened tumor-bearing mice survival time. Low expression of SESN1 had a positive correlation with poor prognosis in patients with NB. RNA-Seq showed that Toll-like receptor (TLR) signaling pathway, and PD-L1 expression and PD-1 checkpoint pathway in cancer were potential downstream targets of SESN1. Knockdown MyD88 or TLRs inhibitor HCQ reversed the effect of knockdown SESN1 in NB cells. High expression of SESN1 was significantly associated with a higher immune score and indicated an active immune microenvironment for patients with NB. CONCLUSIONS: SESN1 functions as a new tumor suppressor gene via TLR signaling pathway in NB.

Histone lysine demethylase 4 family proteins maintain the transcriptional program and adrenergic cellular state of MYCN-amplified neuroblastoma.

Neuroblastoma with MYCN amplification (MNA) is a high-risk disease that has a poor survival rate. Neuroblastoma displays cellular heterogeneity, including more differentiated (adrenergic) and more primitive (mesenchymal) cellular states. Here, we demonstrate that MYCN oncoprotein promotes a cellular state switch in mesenchymal cells to an adrenergic state, accompanied by induction of histone lysine demethylase 4 family members (KDM4A-C) that act in concert to control the expression of MYCN and adrenergic core regulatory circulatory (CRC) transcription factors. Pharmacologic inhibition of KDM4 blocks expression of MYCN and the adrenergic CRC transcriptome with genome-wide induction of transcriptionally repressive H3K9me3, resulting in potent anticancer activity against neuroblastomas with MNA by inducing neuroblastic differentiation and apoptosis. Furthermore, a short-term KDM4 inhibition in combination with conventional, cytotoxic chemotherapy results in complete tumor responses of xenografts with MNA. Thus, KDM4 blockade may serve as a transformative strategy to target the adrenergic CRC dependencies in MNA neuroblastomas.

ALK upregulates POSTN and WNT signaling to drive neuroblastoma.

Neuroblastoma is the most common extracranial solid tumor of childhood. While MYCN and mutant anaplastic lymphoma kinase (ALKF1174L) cooperate in tumorigenesis, how ALK contributes to tumor formation remains unclear. Here, we used a human stem cell-based model of neuroblastoma. Mis-expression of ALKF1174L and MYCN resulted in shorter latency compared to MYCN alone. MYCN tumors resembled adrenergic, while ALK/MYCN tumors resembled mesenchymal, neuroblastoma. Transcriptomic analysis revealed enrichment in focal adhesion signaling, particularly the extracellular matrix genes POSTN and FN1 in ALK/MYCN tumors. Patients with ALK-mutant tumors similarly demonstrated elevated levels of POSTN and FN1. Knockdown of POSTN, but not FN1, delayed adhesion and suppressed proliferation of ALK/MYCN tumors. Furthermore, loss of POSTN reduced ALK-dependent activation of WNT signaling. Reciprocally, inhibition of the WNT pathway reduced expression of POSTN and growth of ALK/MYCN tumor cells. Thus, ALK drives neuroblastoma in part through a feedforward loop between POSTN and WNT signaling.

Identification of MYCN non-amplified neuroblastoma subgroups points towards molecular signatures for precision prognosis and therapy stratification.

BACKGROUND: Despite the extensive study of MYCN-amplified neuroblastomas, there is a significant unmet clinical need in MYCN non-amplified cases. In particular, the extent of heterogeneity within the MYCN non-amplified population is unknown. METHODS: A total of 1566 samples from 16 datasets were identified in Gene Expression Omnibus (GEO) and ArrayExpress. Characterisation of the subtypes was analysed by ConsensusClusterPlus. Independent predictors for subgrouping were constructed from the single sample predictor based on the multiclassPairs package. Findings were verified using immunohistochemistry and CIBERSORTx analysis. RESULTS: We demonstrate that MYCN non-amplified neuroblastomas are heterogeneous and can be classified into 3 subgroups based on their transcriptional signatures. Within these groups, subgroup_2 has the worst prognosis and this group shows a 'MYCN' signature that is potentially induced by the overexpression of Aurora Kinase A (AURKA); whilst subgroup_3 is characterised by an 'inflamed' gene signature. The clinical implications of this subtype classification are significant, as each subtype demonstrates a unique prognosis and vulnerability to investigational therapies. A total of 420 genes were identified as independent subgroup predictors with average balanced accuracy of 0.93 and 0.84 for train and test datasets, respectively. CONCLUSION: We propose that transcriptional subtyping may enhance precision prognosis and therapy stratification for patients with MYCN non-amplified neuroblastomas.