Neuroblastoma (NB) has the highest incidence of all extracranial solid tumors in children and is highly lethal. This study aims to establish a prognostic model of NB with MYCN-related genes. We determined the gene expression profiles of 900 NB samples from the UCSC database and four Gene Expression Omnibus (GEO) data sets, and performed a comprehensive bioinformatics analysis and clinical sample verification. After univariate Cox regression, least absolute shrinkage and selection operator (Lasso), and multivariate Cox regression analyses, four (AKR1C1, CHD5, PDE4DIP, and PRKACB) genes were finally selected and used to construct a risk score prognostic model. In the UCSC data set, the high-risk group exhibited a significantly worse prognosis than the low-risk group. In addition, the nomogram, which includes prognostic markers and clinical factors, demonstrates high prognostic value. Finally, the differential expression of the four genes in the model was verified by quantitative real-time PCR in clinical tissues. These findings of MYCN-related genes provide a new and reliable prognostic model for NB related to MYCN.
Keywords: MYCN; bioinformatics analysis; neuroblastoma; prognostic model; qRT-PCR.