基于GEO数据库和生物信息学分析筛选小鼠心肌缺血再灌注损伤相关的潜在枢纽基因

目的 ·基于基因表达数据库(gene expression omnibus,GEO)运用生物信息学分析筛选与小鼠心肌缺血再灌注损伤(myocardium ischemia-reperfusion injury,MIRI)相关的潜在枢纽(hub)基因.方法 ·从GEO数据库中获取小鼠MIRI数据集GSE61592、GSE83472和GSE160516.利用limma包筛选各数据集中差异表达基因(differentially expressed genes,DEGs),再用稳健排序整合(robust rank aggregation,RRA)方法筛选稳健DEGs.构建稳健DEGs的蛋白质-蛋白质相互作用(protein-protein interaction,PPI)网络,并筛选PPI网络中的子模块和hub基因,利用clusterProfiler包对稳健DEGs、最重要子模块基因和hub基因进行基因本体论(gene ontology,GO)和京都基因与基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,KEGG)通路富集分析.将18只6～8周龄的C57BL/6 J雄性小鼠随机分为假手术(sham)组和MIRI组,每组9只,通过结扎冠状动脉左前降支,缺血30 min再灌注24 h,来构建MIRI模型.采用反转录-定量聚合酶链反应(reverse transcription-quantitative polymerase chain reaction,RT-qPCR)检测hub基因mRNA的表达情况.结果 ·RRA法在3个数据集中共鉴定出294个稳健DEGs.在PPI网络中,共筛选14个子模块,其中模块1最重要;共发现17个关键基因.GO和KEGG分析显示,稳健DEGs、模块1中的基因和hub基因主要涉及调控炎性细胞的迁移、趋化因子和细胞因子及其受体活性、Toll样受体等生物学功能和通路.RT-qPCR结果显示,与sham组相比,MIRI组小鼠心肌中趋化因子(C-C基序)配体4[chemokine (C-C motif) ligand 4,Ccl4]、Ccl6、Ccl7、趋化因子(C-X-C基序)受体4[chemokine (C-X-C motif) receptor 4,Cxcr4]、趋化因子(C-C基序)受体2[chemokine (C-C motif) receptor 2,Ccr2]、信号调节蛋白β1(signal-regulatory protein β l,Sirpb1)、低亲和力免疫球蛋白γ Fc区受体Ⅱb(low affinity immunoglobulin gamma Fe region receptor Ⅱb,Fcgr2b)、白细胞表面抗原Cd53(leukocyte surface antigen CD53,Cd53)、花生四烯酸5-脂氧合酶激活蛋白(arachidonate 5-1ipoxygenase activating protein,Alox5ap)、髓样分化初级反应基因88(myeloid differentiation primary response gene 88,Myd88、巨噬细胞清道夫受体1(macrophage scavenger receptor 1,Msr1)、基质金属肽酶14 (matrix metallopeptidase 14,Mmp14)、髓样细胞上表达的触发受体2(triggering receptor expressed on myeloid cells 2,Trem2)、桩蛋白(leupaxin,Lpxn)的mRNA表达上调,而低亲和力免疫球蛋白γ Fc区受体Ⅲ(low affinity immunoglobulin gamma Fc region receptor Ⅲ,Fcgr3)、补体C1q亚组分亚单位B(complementC1q subcomponent subunitB,C1qb)、去整合素和含金属蛋白酶结构域蛋白(a disintegrin and metalloproteinase domain-containing protein 8,Adam8)的mRNA表达未见差异.回顾文献,17个hub基因中Trem2、Lpxn、Cd53、Alox5ap、Sirpb1、Fcgr2b这6个基因未见报道参与MIRI.结论 ·该研究挖掘出小鼠MIRI相关的6个潜在hub基因,可为进一步探讨MIRI的分子机制和治疗靶点提供新的思路和切入点.