Identification and verification of International Neuroblastoma Staging System (INSS) stage-related genes as potential biomarkers for neuroblastoma prognostic models

Can Qi ¹ Ziwei Zhao ¹ Yanwei Qi ¹ Yun Zhou ¹ Fang Yue ¹ Huizhong Niu ¹ Guochen Duan ² Zhiyong Zhong ¹ Le Wang ^1*

• ¹ Hebei Provincial Children's Hospital, Shijiazhuang, Hebei Province, China
• ² Hebei General Hospital, Shijiazhuang, Hebei Province, China

Background: Neuroblastoma (NB), one of the most common malignant extracranial solid tumours in children, is highly invasive and lethal with limited treatment efficacy. This study aimed to establish a prognostic model of advanced-stage NB.Methods: Differentially expressed genes were screened and validated using two training datasets and one validation dataset from the Therapeutically Applicable Research to Generate Effective Treatments and Gene Expression Omnibus databases. Protein-protein interactions networks were developed using the with MCode plug-in, and the top three key clusters were used to produce candidate genes. We performed gene set enrichment analysis (GSEA), gene ontology analysis (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, immune infiltration, and drug sensitivity analysis to further understand the functions of these candidate genes. Kaplan-Meier (K-M) and receiver operating characteristic (ROC) were used to check their prognosis value. Real-time quantitative polymerase chain reaction (qPCR), western blot (WB) and immunohistochemistry (IHC) were employed to verify the mRNA and protein levels in clinical sample.Results: A total 699 differentially expressed genes were identified, including 294 upregulated and 405 downregulated genes. CNR1, PRKACB, CDKN3, and PCLAF were found to significantly affected the overall survival and event-free survival of neuroblastoma patients and were positively correlated with the INSS advanced stages. The functional analysis of these four genes revealed their cancer-promoting effects and correlations with immune-inflammatory, cell cycle and p53 signalling pathways. After stratifying patients using the established model containing the above four genes, significantly different patterns were observed in terms of infiltrating immune cell proportion, drug sensitivity, and the expression of immune checkpoints. Finally, both of the mRNA and protein expression verification assays demonstrated that the CDKN3 and PCLAF were upregulated, while the PRKACB was downregulated in advanced-stage neuroblastoma tissue samples.Conclusions: The model containing CNR1, PRKACB, CDKN3 and PCLAF can serve as a new prognostic biomarker for predicting the prognosis of patients with neuroblastoma. Findings on immune infiltration and immune checkpoints provide novel insights for the immunotherapy of neuroblastoma.