基于生信分析对强直性脊柱炎关键基因及药物靶点的预测

目的:基于强直性脊柱炎(AS)基因芯片数据,应用生物信息学分析方法识别免疫细胞亚型相关的关键基因及药物靶点预测。方法:GEO数据库下载AS基因表达数据,筛选差异表达基因(DEGs);运用R软件包进行基因本体论(GO)和京都基因和基因组百科全书(KEGG)富集分析;应用CIBERSORT反卷积算法和加权基因共表达网络分析(WGCNA)建立免疫细胞亚型与基因表达的关联性,结合蛋白质互作(PPI)网络筛选的枢纽基因进而确定AS关键基因,药物优先指数(Pi)数据库预测关键基因药物靶点指标。结果:共筛选出126个DEGs,通路富集分析显示DEGs主要富集于淋巴及自然杀伤(NK)细胞免疫调节和毒性信号通路。WGCNA确定了3个模块与AS的中性粒细胞、CD8^(+)T细胞和激活型NK细胞亚型相关性较强:棕色模块与中性粒细胞、绿色模块与中性粒细胞、蓝绿色模块与CD8^(+)T细胞、蓝绿色模块与激活型NK细胞(cor=0.83、0.68、0.52、0.28,均P<0.05)。PPI与模块基因结合筛选出5个关键基因,分别为CXCR1、IKZF1、RUNX3、ID2和ITGB3,基因Pi分析IKZF1和ITGB3在AS中排名较高。结论:CXCR1、IKZF1、RUNX3、ID2和ITGB3在AS中发挥重要作用,IKZF1和ITGB3可作为潜在药物治疗靶点。