Metastatic burden defines clinically and biologically distinct subgroups of stage 4 high-risk neuroblastoma

This study aimed to identify the prognostic subgroups of stage 4 high-risk neuroblastoma based on metastatic burden and explore their distinct clinical and genomic features. Patients aged ≥18 months with stage 4 and metaiodobenzyl guanidine-avid neuroblastoma were enrolled. One hundred and thirty eligible patients were treated under the tandem high-dose chemotherapy scheme. Prognostic significance of metastatic burden measured by the modified Curie score was analyzed using a competing risk approach, and the optimal cut-point was determined. Metastasis-specific subgroups (cut-point: 26) were compared using clinicopathological variables, and differential gene expression analysis and gene set variation analysis (GSVA) were performed using RNA sequencing (RNA-seq). Metastatic burden at diagnosis showed a progressive association with relapse/progression. After applying the cut-point, patients with high metastatic burden showed >3-fold higher risk of relapse/progression than those with low metastatic burden. Moreover, patients with high metastatic burden showed smaller primary tumors and higher biochemical marker levels than those with low metastatic burden. In the genomic analysis, 51 genes were found to be differentially expressed based on the set criteria. GSVA revealed 55 gene sets, which significantly distinguished patients with high metastatic burden from those with low metastatic burden at a false discovery rate <0.25. The results indicated the prognostic significance of metastatic burden in stage 4 high-risk neuroblastoma, and we identified the distinct clinicopathological and genomic features based on metastatic burden. This study may aid in the better understanding and risk-stratification of stage 4 high-risk neuroblastoma patients.