Skip to main content

Metastasis pattern and prognosis in children with neuroblastoma



We aimed to investigate the different metastases and prognoses of neuroblastoma (NB) and determine the risk factors of metastasis.


Data of 1224 patients with NB were obtained from the Surveillance, Epidemiology and End Results database (2010–2018). Pearson’s chi-square test, Kaplan–Meier analysis, multivariable logistic regression and Cox regression analysis were used to determine the factors associated with prognosis.


The overall incidence of NB was an age-adjusted rate of 8.2 patients per 1,000,000 children. In total, 1224 patients were included in our study, with 599 patients (48.9%) exhibiting distant metastases. Compared to patients with non-metastatic NB, a greater proportion of patients with metastatic NB were under 1 year, male, had an adrenal primary site, unilateral tumour, a tumour size > 10 cm, neuroblastoma-not otherwise specified (NB-NOS), second malignant neoplasms and were more likely to choose radiotherapy and chemotherapy. Multivariate Cox regression showed that metastasis was an independent risk factor for overall survival (OS) and cancer-specific survival (CSS). The survival rate of non-metastatic patients with NB was better than those with metastasis (OS: hazard ratio (HR): 0.248, P < 0.001; CSS: HR: 0.267, P < 0.001). The bone and liver were the two most common isolated metastatic sites in NB. However, no statistical difference was observed in OS and CSS between the only bone metastasis group, only liver metastasis group and bone metastasis combined with liver metastasis group (all P > 0.05). Additionally, age at diagnosis > 1 year (odds ratio (OR): 3.295, P < 0 .001), grades III–IV (OR: 26.228, P < 0 .001) and 5–10 cm tumours (OR: 1.781, P < 0 .001) increased the risk of bone metastasis of NB. Moreover, no surgical treatment (OR: 2.441, P < 0 .001) increased the risk of liver metastasis of NB.


Metastatic NB has unique clinicopathological features, with the bone and liver as the most common single metastatic sites of NB. Therefore, more aggressive treatment is recommended for high-risk children with NB displaying distant metastases.


Neuroblastoma (NB) is a highly metastatic malignant tumour, with 70% of patients exhibiting metastasis at the time of diagnosis [1], along with poor prognosis and resistance to conventional chemotherapy. The most common sites of distant metastasis are reported to be the bone marrow and bone, followed by the liver and regional lymph nodes. The 5-year survival rate of low- and medium-risk children has been estimated to reach 75–98% but that of high-risk children is less than 50% [2, 3]. Relapse or progression at the metastatic site is a major cause of NB-related mortality [4].

Currently, internationally recognised treatment-based risk stratification of NB include age, histology, grade, MYCN, 11q aberration and DNA ploidy [5, 6]. The 5-year survival rate of patients with metastatic NB has increased from 20 to > 50% [7], largely owing to the combination of high-dose chemotherapy and autologous stem cell transplantation and differentiation agents and immunotherapy with anti-GD2 monoclonal antibody [8]. However, only a few large-scale clinical studies on the metastasis pattern of NB in children have been reported.

Therefore, this study aims to use the Surveillance, Epidemiology, and End Results (SEER) database to analyse the metastatic pattern of NB, screen for metastasis-related risk factors and determine the clinicopathological features and independent prognostic factors of metastatic patients.

Materials and methods

Patients and ethics

We used the SEER database and SEER-stat software (SEER*Stat 8.3.9) to identify and collect the data of patients aged ≤ 18 years with a confirmed NB diagnosis according to the International Classification of Diseases for Oncology (ICD-O-3) between 2010 and 2018. Data on the following variables were extracted: age at diagnosis, gender, race, year of diagnosis, primary site, tumour grade, laterality, tumour size, histology, chemotherapy, surgery, radiotherapy, second malignant neoplasms (SMN), metastasis, survival time and cause of death. SMN was defined as the appearance of other malignant tumours at least 2 months after NB diagnosis.

Tumour grades were defined as follows: (1) grade I, also called well-differentiated; (2) grade II, moderately differentiated; (3) grade III, poorly differentiated; and (4) grade IV, undifferentiated or anaplastic.

The metastasis of distant organs is defined in SEER as the state of metastasis of distant organs at the time of the first diagnosis of cancer. Among the distant organs, bone metastasis did not include bone marrow metastasis.

The inclusion criteria were as follows: (1) patients with ICD-O-3 Hist/behave = neuroblastoma (ICD-O-3 = 9500) or ganglioneuroblastoma (ICD-O-3 = 9490) and (2) patients with NB as the primary tumour. The exclusion criteria were as follows: (1) lack of follow-up data; (2) diagnosis not confirmed by pathological examinations; and (3) age > 18 years. The research flow is presented in Fig. 1.

Fig. 1
figure 1

Flow chart of the overall study design

All data in the SEER database are anonymised; therefore, ethical approval was not required for this study.

Statistical analysis

The age-standardised incidence of NB based on the 2000 US Standard Population was identified from the cancer registries of the SEER database. The annual percentage change (APC) of incidence was calculated using the weighted least squares method. Moreover, as the cause of death of seven children was unknown, CSS was excluded from the analysis. Kaplan–Meier was used for survival analysis and between-group differences were assessed using Pearson’s chi-square test and Fisher’s exact probability method. Univariate Cox analysis was performed to identify factors (P < 0.05) for inclusion in the multivariate Cox regression models. Furthermore, logistic regression was used to estimate the odds ratio (OR) of metastasis in patients with NB.

Data analyses were performed using SPSS 26.0 for Windows (Armonk, NY: IBM Corp) and R software (version 3.6.3), namely the ‘survival’ and ‘ggplot2’ packages. P < 0.05 was considered statistically significant.


The age-adjusted incidence of NB

A total of 1224 patients were included in our study. The median follow-up time was 42.5 months (inter-quartile range (IQR) 15.0–67.0 months). The overall age-adjusted invasive NB incidence rate for SEER from 2010 to 2018 was 8.2 per 1,000,000 children annually. Notably, the annual incidence rate showed no downtrend decline annually. Furthermore, the APC was 0.1% [95% confidence interval (CI), − 0.2 to 0.4%; P > 0.05] (Fig. 2).

Fig. 2
figure 2

The age-adjusted incidence of NB in the USA (2010–2018). Abbreviations: NB, neuroblastoma

Patient clinicopathological data according to metastases status

A total of 599 (48.9%) patients displayed distant metastases. Moreover, 1003 (81.9%) tumours were NB-not otherwise specified (NOS) and 221 (18.1%) were ganglioneuroblastoma. Table 1 shows the baseline characteristics of the patients. The two groups had significant differences in age, gender, primary site, laterality, tumour size, histology and SMN (all P ≤ 0.05). Furthermore, the metastatic group underwent more radiation therapy and chemotherapy and had lesser surgical interventions than the non-metastatic group. However, both groups were identical in race and tumour grade.

Table 1 Characteristics of patients with metastatic neuroblastoma

Prognostic factors affecting overall survival (OS) and cancer-specific survival (CSS)

Univariate analysis showed that age, primary site, laterality, tumour size, histology, chemotherapy, radiotherapy and metastasis were the risk factors for OS and CSS (all P < 0.05) (Table 2). Additionally, multivariate Cox regression analyses revealed that tumour size > 10 cm and distant metastasis correlated with poor OS (Table 3). Furthermore, tumour size ≥ 5 cm, lack of chemotherapy and distant metastasis correlated with poor CSS. The survival of patients with no metastases was better than those with metastases (OS: hazard ratio (HR): 0.248, P < 0.001; CSS: HR: 0.267, P < 0.001).

Table 2 Univariate Cox proportional hazard regression model of OS and CSS
Table 3 Multivariate Cox proportional hazard regression model of OS and CSS

Metastasis pattern

As shown in Table 4, NB has 15 possible metastatic forms, including 4 single metastases and 11 combinations of metastases. Among patients with single-organ metastasis, we found that the bone was the most common site of metastasis (42.6%), followed by the liver (10.9%), lung (1.3%) and brain (0.2%). For multiple site metastases, the combination of bone and liver metastases was estimated to be 9.5%, accounting for the majority of all cases. Additionally, 329 (54.9%), 107 (17.9%), 22 (3.7%) and 2 (0.3%) patients had one, two, three and four sites of metastasis of NB, respectively.

Table 4 Organ metastasis patterns in patients with neuroblastoma

Furthermore, we compared the differences in metastasis pattern between the ≤ 1-year group and > 1-year group, wherein the > 1-year group had a higher incidence of bone metastasis and a lower incidence of liver metastasis than the ≤ 1-year group (bone metastasis only: 27.6% vs 53.8%, P < 0.001; liver metastasis only: 22.6% vs 2.0%, P < 0.001). In terms of multiple organ metastasis, the incidence of bone + liver and liver + lung in the ≤ 1-year group was higher than that in > 1 year group (bone + liver: 13.2% vs 6.7%, P = 0.007; liver + lung: 1.6% vs 0.0%, P = 0.021).

Survival analysis

Kaplan–Meier analysis revealed that the patients without bone metastasis had a better 5-year OS and CSS than those with bone metastasis (OS: 87.2% vs 66.9%, P < 0.001; CSS: 89.1% vs 69.7%, P < 0.001, Fig. 3A, B). Similar results were also seen in patients without and with liver metastasis (OS: 81.3% vs 74%, P = 0.010; CSS: 83.4% vs 77.6%, P = 0.038, Fig. 3C, D). As bone and liver metastases were most commonly observed in NB, survival differences were observed between the different combinations in this study. As shown in Fig. 4A, B, the OS and CSS of the metastasis group were worse than the without metastasis group. Furthermore, no statistical difference was observed between the only bone metastasis group, only liver metastasis group and the combined bone and liver metastases group (all P > 0.05). As shown in Fig. 4C, D, survival rates decreased as the number of metastatic sites increased in patients with NB.

Fig. 3
figure 3

Kaplan-Meier curves for OS and CSS in NB patients with a single site of metastasis. A, B OS and CSS between patients with and without bone metastasis. C, D OS and CSS between patients with and without liver metastasis. Abbreviations: OS, overall survival; CSS, cancer-specific survival; NB, neuroblastoma

Fig. 4
figure 4

Kaplan–Meier curves for OS and CSS of patients with NB at different metastasis sites. A, B The survival difference among the different metastasis sites in patients with NB. C, D The survival difference among the different numbers of metastasis sites in patients with NB. Abbreviations: OS, overall survival; CSS, cancer-specific survival; NB, neuroblastoma

Risk factors for bone metastasis and liver metastasis

Multivariate logistic regression analysis revealed that age at diagnosis > 1 year (OR: 3.295, 95% CI: 2.335–4.648; P < 0 0.001), tumour grades III–IV (OR: 26.228, 95% CI: 6.972–98.668, P < 0 0.001) and tumour size 5–10 cm (OR: 1.781, 95% CI: 1.050–3.019, P < 0 0.001) increased bone metastasis risk in NB. However, a non-adrenal primary site and lack of chemotherapy or radiotherapy without liver, brain or lung metastasis decreased the risk of bone metastasis (Fig. 5A).

Fig. 5
figure 5

The risk factors for the different primary sites in the developing metastatic disease. A Bone metastasis. B Liver metastasis. Abbreviations: SMN, second malignant neoplasms; CI, confidence interval; OR, odds ratio; NOS, not otherwise specified

When screening the risk factors of hepatoma with liver metastasis, we observed that a lack of surgery (OR: 2.441, 95% CI: 1.545–3.855, P < 0 0.001) increased liver metastasis risk in NB; however, diagnosis > 1 year, non-adrenal primary site, ganglioneuroblastoma, lack of chemotherapy and lack of bone or lung metastasis decreased the risk of liver metastasis (Fig. 5B).


NB is a biologically and clinically heterogeneous group of tumours. It is also the most common paediatric extracranial solid tumour and is responsible for up to 15% of all paediatric oncologic deaths [9]. According to our results, the overall age-adjusted incidence rate of NB was 8.2 per 1,000,000 children annually, and the annual incidence rate showed no significant decrease, which could be highly correlated to the lack of specific clinical manifestations and detection methods in the early stage of the disease. A recent study revealed that nearly 80% of patients with metastatic NB show an ineffective response to current treatment methods [10]. This study compared the clinical features and prognosis of metastatic and non-metastatic NB and screened out high-risk factors for metastasis, thereby aiding early diagnosis. In addition to conventional surgery, chemotherapy and radiotherapy, for high-risk children with distant metastasis, combined tandem autologous transplantation or targeted immunotherapy can be considered.

In this study, 48.9% of patients with NB displayed distant metastasis at the time of diagnosis, which is consistent with previous reports [11]. Compared to patients with non-metastatic NB, a higher proportion of patients with metastatic NB with age > 1 year, male, the presence of an adrenal primary site, unilateral tumours, a tumour size > 10 cm, NB-NOS and SMN features were observed. Consistent with the previous report, our results also showed that primary adrenal NB is prone to metastasis [12, 13]. These differences in biological behaviour could be attributed to the fact that (1) primary adrenal neuroblastoma is more likely to carry structural DNA distortion, including MYCN amplification [13], and (2) the tumour microenvironment is different [12,13,14]. Accordingly, patients with metastatic NB often lose the opportunity of surgery and are more likely to choose radiotherapy and chemotherapy. However, we compared the OS and CSS of the non-metastasis group, the bone metastasis-only group, the liver metastasis-only group and the combined bone and liver metastasis group and found that the survival prognosis of children with distant metastasis was poor, indicating that conventional therapies such as chemotherapy, surgery and radiotherapy have limited efficacy in the treatment of metastatic NB.

Currently, new therapies for NB mainly include immunotherapy, targeted therapy, tandem-autologous transplantation or a combination of multiple methods. In immunotherapy, the combination of GD2 monoclonal antibody (mAb), chemotherapy regimen and autologous NK cells or the combination of anti-GD2 antibody and other antibodies or CART, such as CD47 mAb, have been reported to improve the tumour immunosuppressive microenvironment and enhance the therapeutic effect [15, 16]. Cheung et al. reported that the anti-GD2 antibody has good efficacy in high-risk NB, with an OS of 71% ± 7% at 5 years [17]. Regarding targeted therapy, currently, the majority of research is focused on the ALK, RAS-MAPK and N-Myc-Myc pathways. Mueller et al. reported that 22 (41.5%) of the 53 children treated with chemotherapy combined with anti-GD2 antibody had a shrunken tumour size or underwent complete remission and the 1-year OS rate reached up to 90% [18]. Compared with the previous single transplantation, tandem autologous stem cell transplant showed a better therapeutic effect in high-risk NB. Park et al. reported a clinical trial of 652 high-risk patients with NB, wherein the 3-year event-free survival (EFS) was 61.6% in the tandem transplant group and 48.4% in the single transplant group (P = 0.006) [19]. For high-risk patients with NB displaying distant metastases, the high-intensity multimodal treatment methods, such as tandem autologous transplantation and new targeted or immune treatment measures, are recommended for the improvement of their poor prognosis.

When the different metastasis patterns of NB were compared, we found that bone metastasis (42.6%) and liver metastasis (10.9%) were the most common single-organ metastasis, and the combination of bone and liver metastasis in multi-organ metastasis reached 9.5%. On exploring the risk factors of bone and liver metastases, we found that patients without lung metastasis had a lower risk of bone or liver metastasis. Furthermore, although only lung metastasis (1.3%) is uncommon, lung metastases were accelerated when tumours metastasised to other sites (6.6%). Therefore, it was speculated that once the tumour metastasises to one distant tissue, it could accelerate the metastasis in other tissues, which needs to be considered in a clinical setting. There are many conjectures about the mechanism of NB with distant metastasis: (1) NB originates from the sympathetic nervous system; hence, it might distribute along the sympathetic nervous system; (2) NB is an immune cold tumour with multiple immune escape mechanisms, such as Myc and N-Myc downregulating MHC-I to circumvent T cell-mediated cytotoxicity [20] or tumour-associated macrophages and myeloid suppressor cells that drive immunosuppression in NB [21]; (3) genetic variation and epigenetic changes: 1p and 11q deletion, 17q gain, MYCN amplification and the transcriptome up-regulation of N-Myc and down-regulation of EZH2 are closely related to NB metastasis [22,23,24]; (4) noncoding RNAs are involved in metastasis. For example, lncNB1 promotes the progression of NB by interacting with the ribosomal protein RPL35 [25]. Among them, gene mutation and metastasis of NB have received extensive attention. MYCN amplification was observed in approximately 20 to 25% of NB cases [26]. Moreover, the levels of MYCN correlated with metastatic behaviour. Notably, MYCN contributes to all facets of metastasis: adhesion, motility, invasion and degradation of surrounding matrices [27, 28]. Furthermore, ALK abnormalities were observed in approximately 10% of all sporadic NB [29]. The overexpression of ALK mutants or fusions results in increased migration and invasiveness in NB cells via the upregulation of the MAPK pathway targeting ETV5 [30]. Therefore, it is crucial to elucidate the biological properties of NB and actively explore and develop new and more effective therapies.

Next, we explored the risk factors of patients with bone metastases and liver metastases, respectively, and found that the adrenal primary site was a risk factor for liver and bone metastasis, which could be attributed to the different biological characteristics of NB at different origins, such as ALK mutation levels and MYCN amplification situation [12, 31]. Furthermore, the presence of intra-tumour heterogeneity and different evolutionary trajectories with discrete genotypes in different locations of the same tumour was also observed. Moreover, bone metastasis was more likely to be detected in patients with age > 1 year and was less likely to develop liver metastasis, which is consistent with previous studies [32]. Previous studies have demonstrated age-related changes in MYCN, TERT, PTPRD and Ras pathway alterations, suggesting that the sympathetic nervous system, the tissue of origin of NB, is susceptible to different oncogenic lesions at different times during development [33, 34]. When screening the risk factors, we found that children with NB undergoing surgery had a reduced risk of liver metastasis but had no effect on bone metastasis. Moreover, chemotherapy and radiation therapy were also observed to increase the risk of bone metastasis in NB. Mullassery et al. reported that radical surgery can improve the prognosis of metastatic NB [35]. However, some studies did not observe a survival benefit of surgical resection in patients with high-risk NB [36, 37]. Nonetheless, the reason and mechanism of the unsatisfactory therapeutic effect of chemotherapy and radiation therapy on patients with NB exhibiting bone metastasis require further exploration.

This study acknowledges its limitations. First, missing data in the sample inevitably introduces selection bias. Second, there is a lack of information on chemotherapy regimens, tumour recurrence, immunotherapy, targeted therapy, stem cell transplantation, INSS stage, monthly age and MYCN expansion. Third, owing to the difficulty in obtaining the data on other metastatic sites, only the information available on the liver, lung, bone and brain was used in this study. Finally, this is a retrospective study that requires further prospective, multicentre and large-scale studies to validate our findings.


The population-based analysis of 1224 patients with NB revealed the incidence of metastasis to be 48.9% of all NB cases. Moreover, metastasis was an independent risk factor for OS and CSS. Bone metastasis and liver metastasis were the common distant metastasis patterns of NB, with age > 1 year, grades III–IV and tumour size 5–10 cm increasing the risk of bone metastasis in NB. Additionally, surgical treatment reduced the risk of liver metastasis in NB. For high-risk NB with distant metastasis, a combination of tandem autologous stem cell transplantation, immunotherapy and targeted therapy is recommended to improve its poor prognosis.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author upon reasonable request.


  1. Ara T, DeClerck YA. Mechanisms of invasion and metastasis in human neuroblastoma. Cancer Metastasis Rev. 2006;25(4):645–57.

    Article  PubMed  Google Scholar 

  2. Mlakar V, JurkovicMlakar S, Lopez G, et al. 11q deletion in neuroblastoma: a review of biological and clinical implications. Mol Cancer. 2017;16(1):114.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Matthay KK, Maris JM, Schleiermacher G, et al. Neuroblastoma. Nat Rev Dis Primers. 2016;2:16078.

    Article  PubMed  Google Scholar 

  4. London WB, Bagatell R, Weigel BJ, et al. Historical time to disease progression and progression-free survival in patients with recurrent/refractory neuroblastoma treated in the modern era on Children’s Oncology Group early-phase trials. Cancer. 2017;123(24):4914–23.

    Article  CAS  PubMed  Google Scholar 

  5. Tolbert VP, Matthay KK. Neuroblastoma: clinical and biological approach to risk stratification and treatment. Cell Tissue Res. 2018;372(2):195–209.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Irwin MS, Naranjo A, Zhang FF, et al. Revised neuroblastoma risk classification system: a report from the Children’s Oncology Group. J Clin Oncol. 2021;39(29):3229–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Qiu B, Matthay KK. Advancing therapy for neuroblastoma. Nat Rev Clin Oncol. 2022;19(8):515–33.

    Article  CAS  PubMed  Google Scholar 

  8. Tas ML, Reedijk AMJ, Karim-Kos HE, et al. Neuroblastoma between 1990 and 2014 in the Netherlands: increased incidence and improved survival of high-risk neuroblastoma. Eur J Cancer. 2020;124:47–55.

    Article  CAS  PubMed  Google Scholar 

  9. Swift CC, Eklund MJ, Kraveka JM, et al. Updates in diagnosis, management, and treatment of neuroblastoma. Radiographics. 2018;38(2):566–80.

    Article  PubMed  Google Scholar 

  10. Morgenstern DA, Bagatell R, Cohn SL, et al. The challenge of defining “ultra-high-risk” neuroblastoma. Pediatr Blood Cancer. 2019;66(4):e27556.

    Article  PubMed  Google Scholar 

  11. Al-Battashi A, Al-Rahbi A, Al-Rawahi A, et al. Neuroblastoma among Omani children: clinical characteristics and survival outcome from a dedicated centre. Sultan Qaboos Univ Med J. 2021;21(4):578–84.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Oldridge DA, Truong B, Russ D, et al. Differences in genomic profiles and outcomes between thoracic and adrenal neuroblastoma. J Natl Cancer Inst. 2019;111(11):1192–201.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Vo KT, Matthay KK, Neuhaus J, et al. Clinical, biologic, and prognostic differences on the basis of primary tumor site in neuroblastoma: a report from the international neuroblastoma risk group project. J Clin Oncol. 2014;32(28):3169–76.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Borriello L, Seeger RC, Asgharzadeh S, et al. More than the genes, the tumor microenvironment in neuroblastoma. Cancer Lett. 2016;380(1):304–14.

    Article  CAS  PubMed  Google Scholar 

  15. Theruvath J, Menard M, Smith BAH, et al. Anti-GD2 synergizes with CD47 blockade to mediate tumor eradication. Nat Med. 2022;28(2):333–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Heczey A, Courtney AN, Montalbano A, et al. Anti-GD2 CAR-NKT cells in patients with relapsed or refractory neuroblastoma: an interim analysis. Nat Med. 2020;26(11):1686–90.

    Article  CAS  PubMed  Google Scholar 

  17. Cheung IY, Cheung NV, Modak S, et al. Survival impact of anti-GD2 antibody response in a phase II ganglioside vaccine trial among patients with high-risk neuroblastoma with prior disease progression. J Clin Oncol. 2021;39(3):215–26.

    Article  CAS  PubMed  Google Scholar 

  18. Mueller I, Ehlert K, Endres S, et al. Tolerability, response and outcome of high-risk neuroblastoma patients treated with long-term infusion of anti-GD(2) antibody ch14.18/CHO. mAbs. 2018;10(1):55–61.

    Article  CAS  PubMed  Google Scholar 

  19. Park JR, Kreissman SG, London WB, et al. Effect of tandem autologous stem cell transplant vs single transplant on event-free survival in patients with high-risk neuroblastoma: a randomized clinical trial. JAMA. 2019;322(8):746–55.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Layer JP, Kronmüller MT, Quast T, et al. Amplification of N-Myc is associated with a T-cell-poor microenvironment in metastatic neuroblastoma restraining interferon pathway activity and chemokine expression. Oncoimmunology. 2017;6(6):e1320626.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Miao L, Zhuo Z, Tang J, et al. FABP4 deactivates NF-κB-IL1α pathway by ubiquitinating ATPB in tumor-associated macrophages and promotes neuroblastoma progression. Clin Transl Med. 2021;11(4):e395.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Chen L, Alexe G, Dharia NV, et al. CRISPR-Cas9 screen reveals a MYCN-amplified neuroblastoma dependency on EZH2. J Clin Investig. 2018;128(1):446–62.

    Article  PubMed  Google Scholar 

  23. Carén H, Kryh H, Nethander M, et al. High-risk neuroblastoma tumors with 11q-deletion display a poor prognostic, chromosome instability phenotype with later onset. Proc Natl Acad Sci USA. 2010;107(9):4323–8.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Cetinkaya C, Martinsson T, Sandgren J, et al. Age dependence of tumor genetics in unfavorable neuroblastoma: arrayCGH profiles of 34 consecutive cases, using a Swedish 25-year neuroblastoma cohort for validation. BMC Cancer. 2013;13:231.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Liu PY, Tee AE, Milazzo G, et al. The long noncoding RNA lncNB1 promotes tumorigenesis by interacting with ribosomal protein RPL35. Nat Commun. 2019;10(1):5026.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Niemas-Teshiba R, Matsuno R, Wang LL, et al. MYC-family protein overexpression and prominent nucleolar formation represent prognostic indicators and potential therapeutic targets for aggressive high-MKI neuroblastomas: a report from the children’s oncology group. Oncotarget. 2018;9(5):6416–32.

    Article  PubMed  Google Scholar 

  27. Megison ML, Stewart JE, Nabers HC, et al. FAK inhibition decreases cell invasion, migration and metastasis in MYCN amplified neuroblastoma. Clin Exp Metas. 2013;30(5):555–68.

    Article  CAS  Google Scholar 

  28. Huang M, Weiss WA. Neuroblastoma and MYCN. Cold Spring Harb Perspect Med. 2013;3(10):a014415.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Mossé YP, Laudenslager M, Longo L, et al. Identification of ALK as a major familial neuroblastoma predisposition gene. Nature. 2008;455(7215):930–5.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Mus LM, Lambertz I, Claeys S, et al. The ETS transcription factor ETV5 is a target of activated ALK in neuroblastoma contributing to increased tumour aggressiveness. Sci Rep. 2020;10(1):218.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Fetahu IS, Taschner-Mandl S. Neuroblastoma and the epigenome. Cancer Metastasis Rev. 2021;40(1):173–89.

    Article  PubMed  PubMed Central  Google Scholar 

  32. He B, Mao J, Huang L. Clinical characteristics and survival outcomes in neuroblastoma with bone metastasis based on SEER database analysis. Front Oncol. 2021;11:677023.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Brady SW, Liu Y, Ma X, et al. Pan-neuroblastoma analysis reveals age- and signature-associated driver alterations. Nat Commun. 2020;11(1):5183.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Ognibene M, De Marco P, Parodi S, et al. Genomic analysis made it possible to identify gene-driver alterations covering the time window between diagnosis of neuroblastoma 4S and the progression to stage 4. Int J Mol Sci. 2022;23(12).

  35. Mullassery D, Farrelly P, Losty PD. Does aggressive surgical resection improve survival in advanced stage 3 and 4 neuroblastoma? A systematic review and meta-analysis. Pediatr Hematol Oncol. 2014;31(8):703–16.

    Article  PubMed  Google Scholar 

  36. von Allmen D, Davidoff AM, London WB, et al. Impact of extent of resection on local control and survival in patients from the COG A3973 study with high-risk neuroblastoma. J Clin Oncol. 2017;35(2):208–16.

    Article  Google Scholar 

  37. Simon T, Häberle B, Hero B, et al. Role of surgery in the treatment of patients with stage 4 neuroblastoma age 18 months or older at diagnosis. J Clin Oncol. 2013;31(6):752–8.

    Article  PubMed  Google Scholar 

Download references


The authors sincerely thank the Surveillance, Epidemiology and End Results (SEER) programme for providing the researchers with open resources.


This study was supported by the Fujian Province Natural Science Foundation (grant number 2022J0112).

Author information

Authors and Affiliations



L.S., Y.W.M., and L.Y.B. conceived and designed the study. H.S.H. collected the data. X.S.F., S.G.Y., and L.Y.B. analyzed the data. W.C.Y. supervised the study. L.S. and Y.W.M. wrote the manuscript. L.S. prepared figures 1-2 and L.Y.B. prepared figures 3-4. All authors read and approved the final manuscript prior to submission.

Corresponding author

Correspondence to Chengyi Wang.

Ethics declarations

Ethics approval and consent to participate

Since the SEER database is public, the approval of the ethics committee is not required. The SEER database has hidden the patient’s identity information, so the patient’s informed consent is not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, S., Yin, W., Lin, Y. et al. Metastasis pattern and prognosis in children with neuroblastoma. World J Surg Onc 21, 130 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Neuroblastoma
  • Distant metastases
  • Overall survival
  • Cancer-specific survival
  • SEER