Whole-exome sequencing reveals genetic variants in low-risk and high-risk neuroblastoma

This study aimed to investigate the genetic aberrations in neuroblastoma (NB) by comparing high and low-risk NB patients by whole-exome sequencing (WES) and to reveal the heterogeneity and association between somatic variants and clinical features. Seven NB patients with available clinical data were included in the study (4 in the low-risk group and 3 in the high-risk group). WES was performed and somatic variants associated with NB genes in the COSMIC database were selected through bioinformatics pipeline analysis. Variants were determined using the Integrative Genomics Viewer (IGV). Some gene variations were found in both groups, including variations in oncogene and tumor suppressor genes. In general, candidate gene variations were associated with chromatin remodeling complexes, the RAS pathway, cell proliferation, and DNA repair mechanism. Some variations in CSF1R, MSH6, PTPN11, SOX9, RET, TSC1, and DNMT1 genes were detected only in high-risk patients, while EP300, TET2, MYCN, PRDM1, and ARID2 gene variations were detected only in low-risk patients. When high-risk gene variants were compared with the cBioportal cancer genomic database, two common gene variants (ARID1A and NCOR2) were identified. However, when low-risk gene variants were compared with the cBioportal cancer genomic database, no common genes were found. GO/KEGG enrichment analysis was performed to find relevant biological processes and molecular pathways related to gene variants, which will help to decipher the molecular mechanisms of NB tumorigenesis and the phenotypic differences between high-risk and low-risk patients.