Integrating Bulk-seq and Single-cell-seq Reveals Estrogen and MAPK Pathways Associating with Neuroblastoma Outcome

INTRODUCTION: Neuroblastoma is the most common extracranial solid tumor in children. Patients with high-risk neuroblastoma have a 5-year survival rate less than 50% after extensive treatment. Signaling pathways control cell fate decisions that dictate the behavior of tumor cells. The deregulation of signaling pathways is etiological in cancer cells. Thus, we speculated that the pathway activity of neuroblastoma contains more prognostic information and therapeutic targets. METHODS: Using a footprint-based method, we calculated the activity of fourteen pathways in neuroblastoma. Through stepwise Cox regression analyses, we established a three-gene prognostic signature whose predictive performance was evaluated by external validation. Combining a single-cell sequencing dataset, the most active pathways in high-risk neuroblastoma were found. RESULTS: We found that several pathway activities were correlated with neuroblastoma outcomes. We built a three-gene model comprising DLK1, FLT3, and NTRK1, which exhibited superior internal and external performances. We created a nomogram that combines clinical characteristics to aid in the selection and visualization of high-risk neuroblastoma patients. Furthermore, by integrating a single-cell sequencing dataset, we found that estrogen and MAPK were the most active pathways in high-risk neuroblastoma. CONCLUSION: Our findings suggest that pathway-related therapies may hold promise for the treatment of high-risk neuroblastoma.