Identification of Key Genes Associated with Risk and Prognosis of Neuroblastoma

Neuroblastoma is a childhood malignancy with high morbidity and mortality. We identified key biomarkers associated with neuroblastoma risk and prognosis. The gene modules most associated with neuroblastoma risk were derived by WGCNA. Modular genes were intersected with differentially expressed genes between patients with high-risk (HR) and non-high-risk (NHR) to obtain risk genes, and enrichment analysis was performed. After incorporating risk genes into Cox regression analysis, LASSO algorithm, and K-M survival analysis, key genes were identified and introduced into four external datasets for validation. We performed short time-series expression miner analysis and single-sample genome enrichment analysis. Finally, we evaluated the difference in DNA methylation levels to identify meaningful methylation marks. We identified 5 key genes (ANO6, CPNE2, DST, PLXNC1, SCN3A) for neuroblastoma risk and prognosis, which correlated closely with known neuroblastoma biomarkers. All key genes showed a progressive downregulation trend with increasing risk levels of neuroblastoma. The immune infiltration of 14 immune cells was significantly different between HR-NB and NHR-NB, and most immune cells were negatively correlated with key genes. Furthermore, the expression of ANO6, CPNE2, DST, and PLXNC1 was modified by DNA methylation. This study identified 5 key genes for neuroblastoma risk and prognosis that were potential biomarkers.