加权基因共表达网络分析识别强迫症中的关键基因

目的:利用加权基因共表达网络分析(WGCNA)及套索算法(LASSO)筛选可能影响强迫症发生发展的生物学通路及关键基因,为进一步发现强迫症的分子标志物,确立新的诊断及治疗提供理论依据.方法:从GEO数据库中下载强迫症基因芯片数据GSE78104,包括30例强迫症患者的外周血样本和30例正常人的外周血样本.用WGCNA构建一个基因共表达网络,选择出与强迫症发生高度相关的枢纽模块,分析枢纽模块的生物学功能.通过构建LASSO回归模型筛选与强迫症发生相关的关键基因并检验模型的准确性,最后对关键基因进行差异表达分析.结果:分析发现GSE78104的13模块中,magenta模块与强迫症最为相关(r=0.32,P=0.01).将magen-ta作为枢纽模块进一步进行GO和KEGG分析,表明magenta模块主要涉及中性粒细胞活化、调控白细胞凋亡、白细胞介素-1介导的信号通路和TNF信号通路等功能和信号通路(P<0.05).采用LASSO回归模型筛选得出WDFY3、ZCCHC6、HAL3个关键基因并发现模型准确性较高(AUC=0.78),最后发现3个关键基因在强迫症中均为高表达.结论:WDFY3、ZCCHC6、HAL 3个关键基因可能参与强迫症的发生发展过程.