Identification of hub genes and their correlation with infiltration of immune cells in MYCN positive neuroblastoma based on WGCNA and LASSO algorithm

The prognosis of MYCN positive NB is poor, and there is no targeted drug for N-myc at present. This study aims to screen out hub genes closely related to MYCN, analyze the relationship between hub genes and NB microenvironment, and provide basis for molecular targeted therapy of MYCN positive NB.

ZNF695, CHEK1 and C15ORF42 are highly expressed in MYCN positive NB, and their expression levels are negatively correlated with the prognosis of children with NB. The infiltration levels of Activated CD4 T cell and Type2 T helper cell increased in the microenvironment of MYCN positive NB and were significantly positively correlated with the expression levels of the three hub genes. The results of this study provide that ZNF695, CHEK1 and C15ORF42 may be potential prognostic markers and immunotherapy targets for MYCN positive NB.

We combined the microarray data of GSE45547 (n=649) and GSE49710 (n=498), screened the DEGs between MYCN positive (n=185) and MYCN negative NB (n=951), performed WGCNA, Lasso regression and Roc analyses on the merged matrix, and obtained the hub genes related to MYCN in the training group. We performed ssGSEA on the experimental group to calculate the infiltration level of 28 kinds of immune cells in each sample, compared the differences of immune cell infiltration between MYCN positive and MYCN negative group. The influences of hub genes on the distribution of each immune cell were also analyzed by ssGSEA. The expression differences of the three hub genes were verified in the E-MTAB-8248 cohort (n=223), and the correlation between hub genes and prognosis of NB was calculated by Kaplan-Meier method in GSE62564 (n=498) and the validation group. We also verified the expression differences of hub genes by qRT-PCR in SK-N-BE(2), SKNDZ, Kelly and SH-SY5Y cell lines.

Here were 880 DEGs including 420 upregulated and 460 downregulated genes in MYCN positive NB in the training group. Overlap of the DEGs and WGCNA networks identified four shared genes, namely, ZNF695, CHEK1, C15ORF42 and EXO1, as candidate hub genes in MYCN positive NB. Three core genes, ZNF695, CHEK1 and C15ORF42, were finally identified by Lasso regression and Roc analyses. ZNF695, CHEK1 and C15ORF42 were highly expressed in MYCN positive NB tissues and cell lines. These three genes were closely related to the prognosis of children with NB. Except that Activated CD4 T cell and Type2 T helper cell increased, the infiltration levels of the other 26 cells decreased significantly in MYCN positive NB tissues. The infiltration levels of Type2 T helper cell and Activated CD4 T cell were also significantly positively correlated with the expression levels of the three hub genes.