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.
Keywords: Biomarkers; DNA methylation; Neuroblastoma; STEM analysis; WGCNA.
© 2022. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.