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.

Identification of an Ultra-High-Risk Subgroup of Neuroblastoma Patients within the High-Risk Cohort Using a Computed Tomography-Based Radiomics Approach

RATIONALE AND OBJECTIVES

To identify ultra-high-risk (UHR) neuroblastoma patients who experienced disease-related mortality within 18 months of diagnosis within the high-risk cohort using computed tomography (CT)-based radiomics analysis.

MATERIALS AND METHODS

A retrospective analysis was conducted on 105 high-risk neuroblastoma patients, divided into a training set (n = 74) and a test set (n = 31). Radiomics features were extracted and selected from arterial phase CT images, and an optimal radiomics signature was established using the support vector machine algorithm. Evaluation metrics, including area under the curve (AUC) and 95% confidence interval (CI), were calculated. Furthermore, the fit and clinical benefit of the signature, along with its correlation with overall survival (OS), were analyzed.

RESULTS

The optimal radiomics signature comprised 11 features. In the training set, AUC and accuracy were 0.911 (95% CI: 0.840-0.982) and 0.892, respectively. In the test set, AUC and accuracy were 0.828 (95% CI: 0.669-0.987) and 0.839, respectively. There was no significant difference between predicted probability and actual probability, and the signature demonstrated net benefit. The concordance index of this signature for predicting OS was 0.743 (95% CI: 0.672-0.814) in the training set and 0.688 (95% CI: 0.566-0.810) in the test set. Moreover, the signature achieved AUC values of 0.832, 0.863, and 0.721 for 1-year, 2-year, and 3-year OS in the training set, and 0.870, 0.836, and 0.638 in the test set for the respective time periods.

CONCLUSION

The utilization of CT-based radiomics signature to identify an UHR subgroup of neuroblastoma patients within the high-risk cohort can help aid in predicting early disease progression.

FGFR1 is a potential therapeutic target in neuroblastoma

BACKGROUND

FGFR1 regulates cell-cell adhesion and extracellular matrix architecture and acts as oncogene in several cancers. Potential cancer driver mutations of FGFR1 occur in neuroblastoma (NB), a neural crest-derived pediatric tumor arising in sympathetic nervous system, but so far they have not been studied experimentally. We investigated the driver-oncogene role of FGFR1 and the implication of N546K mutation in therapy-resistance in NB cells.

METHODS

Public datasets were used to predict the correlation of FGFR1 expression with NB clinical outcomes. Whole genome sequencing data of 19 paired diagnostic and relapse NB samples were used to find somatic mutations. In NB cell lines, silencing by short hairpin RNA and transient overexpression of FGFR1 were performed to evaluate the effect of the identified mutation by cell growth, invasion and cologenicity assays. HEK293, SHSY5Y and SKNBE2 were selected to investigate subcellular wild-type and mutated protein localization. FGFR1 inhibitor (AZD4547), alone or in combination with PI3K inhibitor (GDC0941), was used to rescue malignant phenotypes induced by overexpression of FGFR1 wild-type and mutated protein.

RESULTS

High FGFR1 expression correlated with low relapse-free survival in two independent NB gene expression datasets. In addition, we found the somatic mutation N546K, the most recurrent point mutation of FGFR1 in all cancers and already reported in NB, in one out of 19 matched primary and recurrent tumors. Loss of FGFR1 function attenuated invasion and cologenicity in NB cells, whereas FGFR1 overexpression enhanced oncogenicity. The overexpression of FGFR1^N546K protein showed a higher nuclear localization compared to wild-type protein and increased cellular invasion and cologenicity. Moreover, N546K mutation caused the failure in response to treatment with FGFR1 inhibitor by activation of ERK, STAT3 and AKT pathways. The combination of FGFR1 and PI3K pathway inhibitors was effective in reducing the invasive and colonigenic ability of cells overexpressing FGFR1 mutated protein.

CONCLUSIONS

FGFR1 is an actionable driver oncogene in NB and a promising therapy may consist in targeting FGFR1 mutations in patients with therapy-resistant NB.

MYCN amplification plus 1p36 loss of heterozygosity predicts ultra high risk in bone marrow metastatic neuroblastoma

Background

This study aimed to better understand the prognostic effect of multiple genetic markers and identify more subpopulations at ultra high risk of poor outcome in bone marrow (BM) metastatic neuroblastoma (NB).

Methods

We screened the MYCN, 1p36 and 11q23 loss of heterozygosity (LOH) statuses of 154 patients by interphase fluorescence in situ hybridization of BM cells. The clinical characteristics of patients with the three markers and their associations with prognosis were analysed.

Results

MYCN amplification and LOH at 1p36 and 11q23 were identified in 16.2%, 33.1% and 30.5% of patients, respectively. There were strong associations between MYCN amplification and 1p36 LOH as well as 11q23 LOH. Both MYCN amplification and 1p36 LOH were strongly associated with high levels of lactate dehydrogenase (LDH) and neuron‐specific enolase, more than 3 metastatic organs, and more events. 11q23 LOH occurred mainly in patients older than 18 months, and those who had high LDH levels. In univariate analysis, patients with MYCN amplification had poorer prognosis than those without. Patients with 1p36 LOH had a 3‐year event‐free survival (EFS) and overall survival lower than those without. 11q23 LOH was associated with poorer EFS only for patients without MYCN amplification. In a multivariate model, MYCN amplification was independently associated with decreased EFS in all cohorts. 11q23 LOH was an independent prognostic factor for patients without MYCN amplification, whereas 1p36 LOH was not an independent marker regardless of MYCN amplification. Compared with all cohorts, patients with both MYCN amplification and 1p36 LOH had the worst outcome and clinical features.

Conclusions

Patients with both MYCN amplification and 1p36LOH had the worst survival rate, indicating an ultra high‐risk group. Our results may be applied in clinical practice for accurate risk stratification in future studies.

Somatic mutations enriched in cis-regulatory elements affect genes involved in embryonic development and immune system response in neuroblastoma

Noncoding cis-regulatory variants have gained interest as cancer drivers, yet progress in understanding their significance is hindered by the numerous challenges and limitations of variant prioritization. To overcome these limitations, we focused on active cis-regulatory elements (aCRE) in order to design a customized panel for the deep sequencing of 56 neuroblastoma tumor and normal DNA sample pairs. In order to search for driver mutations, aCREs were defined by reanalysis of H3K27ac ChiP-seq peaks in 25 neuroblastoma cell lines. These regulatory genomic regions were tested for an excess of somatic mutations and assessed for statistical significance using a global approach that accounted for chromatin accessibility and replication timing. Additional validation was provided by whole genome sequence analysis of 151 neuroblastomas. Analysis of Hi-C data determined the presence of candidate target genes interacting with mutated regions. An excess of somatic mutations in aCREs of diverse genes were identified, including IPO7, HAND2, and ARID3A. CRISPR-Cas9 editing was utilized to assess the functional consequences of mutations in the IPO7 aCRE. Patients with noncoding mutations in aCREs showed inferior overall and event-free survival independent of age at diagnosis, stage, risk stratification, and MYCN status. Expression of aCRE-interacting genes correlated strongly with negative prognostic markers and low survival rates. Moreover, a convergence between the biological functions of aCRE target genes and transcription factors with mutated binding motifs was associated with embryonic development and immune system response. Overall, this strategy enabled the identification of somatic mutations in regulatory elements that collectively promote neuroblastoma tumorigenesis.

Identification of novel prognostic markers of survival time in high-risk neuroblastoma using gene expression profiles

Neuroblastoma is the most common extracranial solid tumor in childhood. Patients in high-risk group often have poor outcomes with low survival rates despite several treatment options. This study aimed to identify a genetic signature from gene expression profiles that can serve as prognostic indicators of survival time in patients of high-risk neuroblastoma, and that could be potential therapeutic targets. RNA-seq count data was downloaded from UCSC Xena browser and samples grouped into Short Survival (SS) and Long Survival (LS) groups. Differential gene expression (DGE) analysis, enrichment analyses, regulatory network analysis and machine learning (ML) prediction of survival group were performed. Forty differentially expressed genes (DEGs) were identified including genes involved in molecular function activities essential for tumor proliferation. DEGs used as features for prediction of survival groups included EVX2, NHLH2, PRSS12, POU6F2, HOXD10, MAPK15, RTL1, LGR5, CYP17A1, OR10AB1P, MYH14, LRRTM3, GRIN3A, HS3ST5, CRYAB and NXPH3. An accuracy score of 82% was obtained by the ML classification models. SMIM28 was revealed to possibly have a role in tumor proliferation and aggressiveness. Our results indicate that these DEGs can serve as prognostic indicators of survival in high-risk neuroblastoma patients and will assist clinicians in making better therapeutic and patient management decisions.

Genetic predisposition and chromosome instability in neuroblastoma

Neuroblastoma (NB) is a pediatric tumor of embryonic origin. About 1-2% of all NBs are familial cases, and genetic predisposition is suspected for the remaining cases. During the last decade, genome-wide association studies (GWAS) and high-throughput sequencing approaches have been used to identify associations among common and rare genetic variants and NB risk. Substantial data has been produced by large patient cohorts that implicate various genes in NB tumorigenesis, such as CASC15, BARD1, CHEK2, LMO1, LIN28B, AXIN2, BRCA1, TP53, SMARCA4, and CDK1NB. NB, as well as other pediatric cancers, has few recurrent mutations but several copy number variations (CNVs). Almost all NBs show both numerical and structural CNVs. The proportion between numerical and structural CNVs differs between localized and metastatic tumors, with a greater prevalence of structural CNVs in metastatic NB. This genomic chaos frequently identified in NBs suggests that chromosome instability (CIN) could be one of the major actors in NB oncogenesis. Interestingly, many NB-predisposing variants occur in genes involved in the control of genome stability, mitosis, and normal chromosome separation. Here, we discuss the relationship between genetic predisposition and CIN in NB.

Gene Expression Signature of Acquired Chemoresistance in Neuroblastoma Cells

Drug resistance of childhood cancer neuroblastoma is a serious clinical problem. Patients with relapsed disease have a poor prognosis despite intense treatment. In the present study, we aimed to identify chemoresistance gene expression signatures in vincristine resistant neuroblastoma cells. We found that vincristine-resistant neuroblastoma cells formed larger clones and survived under reduced serum conditions as compared with non-resistant parental cells. To identify the possible mechanisms underlying vincristine resistance in neuroblastoma cells, we investigated the expression profiles of genes known to be involved in cancer drug resistance. This specific gene expression patterns could predict the behavior of a tumor in response to chemotherapy and for predicting the prognosis of high-risk neuroblastoma patients. Our signature could help chemoresistant neuroblastoma patients in avoiding useless and harmful chemotherapy cycles.

Through Predictive Personalized Medicine

Neuroblastoma (NBM) is a deadly form of solid tumor mostly observed in the pediatric age. Although survival rates largely differ depending on host factors and tumor-related features, treatment for clinically aggressive forms of NBM remains challenging. Scientific advances are paving the way to improved and safer therapeutic protocols, and immunotherapy is quickly rising as a promising treatment that is potentially safer and complementary to traditionally adopted surgical procedures, chemotherapy and radiotherapy. Improving therapeutic outcomes requires new approaches to be explored and validated. In-silico predictive models based on analysis of a plethora of data have been proposed by Lombardo et al. as an innovative tool for more efficacious immunotherapy against NBM. In particular, knowledge gained on intracellular signaling pathways linked to the development of NBM was used to predict how the different phenotypes could be modulated to respond to anti-programmed cell death-ligand-1 (PD-L1)/programmed cell death-1 (PD-1) immunotherapy. Prediction or forecasting are important targets of artificial intelligence and machine learning. Hopefully, similar systems could provide a reliable opportunity for a more targeted approach in the near future.

Kinome multigenic panel identified novel druggable EPHB4-V871I somatic variant in high-risk neuroblastoma

Neuroblastoma (NB) is the most common extracranial neoplasm in children. The overall outcome for high‐risk NB patients is still unacceptable, therefore, it is critical to deeply understand molecular mechanisms associated with NB, which in turn can be utilized for developing drugs towards the treatment of NB. Protein kinases (TKs) play an essential role in the regulation of cell survival and proliferation. Different kinases, such as anaplastic lymphoma kinase (ALK), Aurora kinase, RET receptor tyrosine kinase, are potential therapeutic targets in various cancers, including NB. We analysed a cohort of 45 high‐risk NB patients and 9 NB cell lines by a targeted—(t)NGS custom gene panel (genes codifying for the kinase domains of 90 TKs). We identified somatic variants in four TK genes (ALK, EPHB4, LMTK3 and EPHB6) in NB patients and we functionally characterized an interesting somatic variant, V871I, in EPHB4 gene. EPHB4 plays a crucial role in cardiovascular development and regulates vascularization in cancer‐promoting angiogenesis, tumour growth and metastasis. Several EPHB4 mutations have previously been identified in solid and haematological tumour specimens but EPHB4 mutations were not described until now in NB. Interestingly, a re‐analysis of public CGH‐array showed that the EPHB4 gain is associated with advanced diseases in NB. We further demonstrated that higher EPHB4 expression is correlated to stage 4 of NB and with poor overall survival. Additionally, we also revealed that the EPHB4‐V871I accounts for increased proliferation, migration and invasion properties in two NB cell lines by acting on VEGF, c‐RAF and CDK4 target genes and by increasing the phosphorylation of ERK1‐2 pathway. The use of two EPHB4 inhibitors, JI‐101 and NVP‐BHG712, was able to rescue the phenotype driven by the variant. Our study suggested that EPHB4 is a promising therapeutic target in high‐risk NB.

Metabolomic profiling in neuroblastoma

BACKGROUND AND OBJECTIVES

Previous studies on several cancer types show that metabolomics provides a potentially useful noninvasive screening approach for outcome prediction and accurate response to treatment assessment. Neuroblastoma (NB) accounts for at least 15% of cancer-related deaths in children. Although current risk-based treatment approaches in NB have resulted in improved outcome, survival for high-risk patients remains poor. This study aims to evaluate the use of metabolomics for improving patients' risk-group stratification and outcome prediction in NB.

DESIGN AND METHODS

Plasma samples from 110 patients with NB were collected at diagnosis prior to starting therapy and at the end of treatment if available. Metabolomic analysis of samples was carried out by ultra-performance liquid chromatography-time of flight mass spectrometry (UPLC-MS).

RESULTS

The metabolomic analysis was able to identify different plasma metabolic profiles in high-risk and low-risk NB patients at diagnosis. The metabolic model correctly classified 16 high-risk and 15 low-risk samples in an external validation set providing 84.2% sensitivity (60.4-96.6, 95% CI) and 93.7% specificity (69.8-99.8, 95% CI). Metabolomic profiling could also discriminate high-risk patients with active disease from those in remission. Notably, a plasma metabolomic signature at diagnosis identified a subset of high-risk NB patients who progressed during treatment.

CONCLUSIONS

To the best of our knowledge, this is the largest NB study investigating the prognostic power of plasma metabolomics. Our results support the potential of metabolomic profiling for improving NB risk-group stratification and outcome prediction. Additional validating studies with a large cohort are needed.

Gene signatures associated with genomic aberrations predict prognosis in neuroblastoma

BACKGROUND

Neuroblastoma (NB) is a heterogeneous disease with respect to genomic abnormalities and clinical behaviors. Despite recent advances in our understanding of the association between the genetic aberrations and clinical features, it remains one of the major challenges to predict prognosis and stratify patients for determining personalized therapy in this disease. The aim of this study was to develop an effective prognosis prediction model for NB patients.

METHODS

We integrated diverse computational analyses to define gene signatures that reflect MYCN activity and chromosomal aberrations including deletion of chromosome 1p (Chr1p_del) and chromosome 11q (Chr11q_del) as well as chromosome 11q whole loss (Chr11q_wls). We evaluated the prognostic and predictive values of these signatures in seven NB gene expression datasets (the number of samples ranges from 94 to 498, with a total of 2120) generated from both RNA sequencing and microarray platforms.

RESULTS

MYCN signature was a more effective prognostic marker than MYCN amplification status and MYCN expression. Similarly, the Chr1p_del score was more prognostic than Chr1p status. The activity scores of MYCN, Chr1p_del and Chr11q_del were associated with poor prognosis, while the Chr11q_wls score was linked to good outcome. We integrated the activity scores of MYCN, Chr1p_del, Chr11q_del, and Chr11q_wls and clinical variables into an integrative prognostic model, which displayed significant performance over the clinical variables or each genomic aberration alone.

CONCLUSIONS

Our integrative gene signature model shows a significantly improved forecast performance with prognostic and predictive information, and thereby can be served as a biomarker to stratify NB patients for prognosis evaluation and surveillance programs.

Association of PARP1 polymorphisms with response to chemotherapy in patients with high-risk neuroblastoma

The genetic aetiology and the molecular mechanisms that characterize high‐risk neuroblastoma are still little understood. The majority of high‐risk neuroblastoma patients do not take advantage of current induction therapy. So far, one of the main reasons liable for cancer therapeutic failure is the acquisition of resistance to cytotoxic anticancer drugs, because of the DNA repair system of tumour cells. PARP1 is one of the main DNA damage sensors involved in the DNA repair system and genomic stability. We observed that high PARP1 mRNA level is associated with unfavourable prognosis in 3 public gene expression NB patients’ datasets and in 20 neuroblastomas analysed by qRT‐PCR. Among 4983 SNPs in PARP1, we selected two potential functional SNPs. We investigated the association of rs907187, in PARP1 promoter, and rs2048426 in non‐coding region with response chemotherapy in 121 Italian patients with high‐risk NB. Results showed that minor G allele of rs907187 associated with induction response of patients (P = .02) and with decrease PARP1 mRNA levels in NB cell line (P = .003). Furthermore, rs907187 was predicted to alter the binding site of E2F1 transcription factor. Specifically, allele G had low binding affinity with E2F1 whose expression positively correlates with PARP1 expression and associated with poor prognosis of patients with NB. By contrast, we did not find genetic association for the SNP rs2048426. These data reveal rs907187 as a novel potential risk variant associated with the failure of induction therapy for high‐risk NB.

CCR5 and CXCL12 allelic variants: Possible association with childhood neuroblastoma susceptibility?

Neuroblastoma (NB) is a heterogeneous and particularly malignant childhood neoplasm in its higher stages, prone to form metastasis in selected organs and for which there is still no efficient treatment available beyond surgery. Evidence indicates that chemokines and their receptors present involvement as mediators of neuroinflammation and have a neurophysiological role. In the present study, we aimed to verify if CCR5 (rs333) and CXCL12 (rs1801157) allelic variants were associated with NB. For CCR5 (rs333) D32 carriers (OR: 5.96, IC: 2.21-16.06) and for CXCL12 genotype 3'A/3'A (OR:26.18, IC:6.15-111.4) there were statistically significant differences as well to allelic frequency (OR:4.20, IC: 2.19-8.03). Although no correlation was verified regarding prognostic parameters for both CCR5 and CXCL12 polymorphic variants, these polymorphisms may be associated with NB susceptibility which deserve attention for future investigations.

Transcription Factors Involved in Tumorigenesis Are Over-Represented in Mutated Active DNA-Binding Sites in Neuroblastoma

The contribution of coding mutations to oncogenesis has been largely clarified, whereas little is known about somatic mutations in noncoding DNA and their role in driving tumors remains controversial. Here, we used an alternative approach to interpret the functional significance of noncoding somatic mutations in promoting tumorigenesis. Noncoding somatic mutations of 151 neuroblastomas were integrated with ENCODE data to locate somatic mutations in regulatory elements specifically active in neuroblastoma cells, nonspecifically active in neuroblastoma cells, and nonactive. Within these types of elements, transcription factors (TF) were identified whose binding sites were enriched or depleted in mutations. For these TFs, a gene expression signature was built to assess their implication in neuroblastoma. DNA- and RNA-sequencing data were integrated to assess the effects of those mutations on mRNA levels. The pathogenicity of mutations was significantly higher in transcription factor binding site (TFBS) of regulatory elements specifically active in neuroblastoma cells, as compared with the others. Within these elements, there were 18 over-represented TFs involved mainly in cell-cycle phase transitions and 15 under-represented TFs primarily regulating cell differentiation. A gene expression signature based on over-represented TFs correlated with poor survival and unfavorable prognostic markers. Moreover, recurrent mutations in TFBS of over-represented TFs such as EZH2 affected MCF2L and ADP-ribosylhydrolase like 1 expression, among the others. We propose a novel approach to study the involvement of regulatory variants in neuroblastoma that could be extended to other cancers and provide further evidence that alterations of gene expression may have relevant effects in neuroblastoma development. SIGNIFICANCE: These findings propose a novel approach to study regulatory variants in neuroblastoma and suggest that noncoding somatic mutations have relevant implications in neuroblastoma development.

The challenge of defining "ultra-high-risk" neuroblastoma

Given the biological and clinical heterogeneity of neuroblastoma, risk stratification is vital to determining appropriate treatment. Historically, most patients with high-risk neuroblastoma (HR-NBL) have been treated uniformly without further stratification. Attempts have been made to identify factors that can be used to risk stratify these patients and to characterize an "ultra-high-risk" (UHR) subpopulation with particularly poor outcome. However, among published data, there is a lack of consensus in the definition of the UHR population and heterogeneity in the endpoints and statistical methods used. This review summarizes our current understanding of stratification of HR-NBL and discusses the complex issues in defining UHR neuroblastoma.

MAPT (Tau) expression is a biomarker for an increased rate of survival in pediatric neuroblastoma

Although the impact of MAPT (Tau) expression has been well documented for neuronal cells in the context of tauopathies and neurodegenerative diseases, the impact and role of Tau expression in cancer, and specifically cancers of neuronal origin, is in its infancy. To determine the correlation between MAPT expression and survival in pediatric neuroblastoma, MAPT gene expression for samples from the TARGET pediatric neuroblastoma dataset was assessed. Initial analyses indicated that increased MAPT expression correlated with increased overall survival in neuroblastoma but not in ovarian cancer. Expression of apoptosis- and proliferation-effector genes in the neuroblastoma samples was consistent with the MAPT related survival result. Furthermore, we determined that higher neuroblastoma expression of APP also associated with neurodegeneration, correlated with better neuroblastoma survival rates. In sum, Gene expression associated with neuronal degenerative diseases was associated with a better neuroblastoma outcome. Abbreviations: ALS: Amyotrophic Lateral Sclerosis; APP: Amyloid Precursor Protein gene; CASP3: Caspase 3 gene; CASP9: Caspase 9 gene; H2AFX: H2A histone family, member X gene; HIST1H2AL: Histone H2A type 1 gene; HIST1H2BK: Histone H2B type 1-K gene; HIST1H3J: Histone H3J gene; HIST1H4B: Histone H4B gene; HIST2H2BE: Histone H2B type 2-E gene; HUGO: human genome organization; KM: Kaplan-Meier survival curve; MAPT: Tau gene; OV: Ovarian cancer; SNCA: alpha-syneculin gene; TARDBP: Transactive response DNA binding protein 43 kDa; TCGA: the cancer genome atlas.

3-Methoxytyramine: An independent prognostic biomarker that associates with high-risk disease and poor clinical outcome in neuroblastoma patients

INTRODUCTION

Prognosis of neuroblastoma patients is very diverse, indicating the need for more accurate prognostic parameters. The excretion of catecholamine metabolites by most neuroblastomas is used for diagnostic purposes, but their correlation with prognosis has hardly been investigated. Therefore, we performed an in-depth analysis of a panel of elevated urinary catecholamine metabolites at diagnosis and their correlation with prognosis.

PATIENTS AND METHODS

Retrospective study of eight urinary catecholamine metabolites in a test (n = 96) and validation (n = 205) cohort of patients with neuroblastoma (all stages) at diagnosis.

RESULTS

Multivariate analyses, including risk factors such as stage and MYCN amplification, revealed that 3-methoxytyramine (3MT) was an independent risk factor for event-free survival (EFS) and overall survival (OS). Furthermore, only 3MT appeared to be an independent risk factor for both EFS and OS in high-risk patients, which was independent of modern high-risk therapy and immunotherapy. Among high-risk patients, those with elevated 3MT and older than 18 months had an extremely poor prognosis compared to patients with non-elevated 3MT and younger than 18 months (5-year EFS of 14.3% ± 4% and 66.7% ± 18%, respectively, p = 0.001; 5-year OS of 21.8% ± 5% and 87.5% ± 12%, respectively, p < 0.001).

CONCLUSIONS

Elevated 3MT at diagnosis was associated with high-risk disease and poor prognosis. For high-risk patients, elevated 3MT at diagnosis was the only significant risk factor for EFS and OS. 3MT was also able to identify subgroups of high-risk patients with favourable and extremely poor prognosis.