上海交通大学学报(医学版) ›› 2022, Vol. 42 ›› Issue (1): 51-62.doi: 10.3969/j.issn.1674-8115.2022.01.008

• 论著 · 基础研究 • 上一篇    下一篇

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

王建茹(), 彭广操, 朱明军()   

  1. 河南中医药大学第一附属医院心内科,郑州 450099
  • 收稿日期:2021-09-27 出版日期:2022-01-28 发布日期:2022-02-18
  • 通讯作者: 朱明军 E-mail:wangjianru1988@126.com;zhumingjun317@163.com
  • 作者简介:王建茹(1988—),男,主治医师,博士;电子信箱:wangjianru1988@126.com
  • 基金资助:
    国家自然科学基金(82004311);河南省科技攻关项目(202102310492);河南省中医药科学研究专项课题(2019JDZX2013);河南省博士后科研项目(202001045);中原院士基金-中原名医

Screening potential hub genes associated with myocardial ischemia-reperfusion injury in mice based on GEO database and bioinformatics analysis

Jianru WANG(), Guangcao PENG, Mingjun ZHU()   

  1. Department of Cardiology, The First Affiliated Hospital of Henan University of Chinese Medicine, Zhengzhou 450099, China
  • Received:2021-09-27 Online:2022-01-28 Published:2022-02-18
  • Contact: Mingjun ZHU E-mail:wangjianru1988@126.com;zhumingjun317@163.com
  • Supported by:
    National Natural Science Foundation of China(82004311);Key Research and Development and Promotion Special Project of Henan Province (Scientific Problem Tackling) Project(202102310492);Traditional Chinese Medicine Scientific Research Special Project of Henan Province(2019JDZX2013);Postdoctoral Research Project of Henan Province(202001045);Zhongyuan Academician Foundation-Famous Doctor of Zhongyuan

摘要:

目的·基于基因表达数据库(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 β 1,Sirpb1)、低亲和力免疫球蛋白γ Fc区受体Ⅱb(low affinity immunoglobulin gamma Fc region receptor Ⅱb,Fcgr2b白细胞表面抗原Cd53(leukocyte surface antigen CD53,Cd53花生四烯酸5-脂氧合酶激活蛋白(arachidonate 5-lipoxygenase 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(complement C1q subcomponent subunit B,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的分子机制和治疗靶点提供新的思路和切入点。

关键词: 心肌缺血再灌注损伤, 生物信息学分析, 差异表达基因, 关键基因, 基因表达数据库

Abstract:

Objective·To screen the potential hub genes associated with myocardial ischemia-reperfusion injury (MIRI) in mice by bioinformatics analysis based on gene expression omnibus (GEO) database.

Methods·The mouse MIRI data sets GSE61592, GSE83472 and GSE160516 were obtained from GEO database. The differentially expressed genes (DEGs) in each data set were screened by limma package, and then robust DEGs were screened by robust sorting integration (RRA) method. The protein-protein interaction (PPI) network of robust DEGs was constructed, and the submodules and hub genes in the PPI network were screened. The clusterProfiler package was used to analyze the robust DEGs, the most important submodule genes and hub genes by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Eighteen male C57BL/6 J mice aged 6?8 weeks were randomly divided into sham group and MIRI group, 9 mice each group. The MIRI model was constructed by left anterior descending coronary artery ligation for 30 min ischemia and 24 h reperfusion, and the mRNA expression of hub genes was detected by reverse transcription-quantitative polymerase chain reaction (RT-qPCR).

Results·RRA method identified 294 robust DEGs in three data sets. In PPI network, a total of 14 sub-modules were screened, of which module 1 was the most important and 17 hub genes were found. GO and KEGG analysis showed that the robust DEGs, module 1 genes, and the hub genes were mainly involved in regulating the migration of inflammatory cells, the activity of chemokines and cytokines and their receptors, Toll-like receptors and other biological function and signaling pathways. RT-qPCR results showed that compared with the sham group, the expressions of chemokine (C-C motif) ligand 4 (Ccl4), Ccl6, Ccl7, chemokine (C-X-C motif) receptor 4 (Cxcr4), chemokine (C-C motif) receptor 2 (Ccr2), signal-regulatory protein β1 (Sirpb1), low affinity immunoglobulin gamma Fc region receptor Ⅱb (Fcgr2b), leukocyte surface antigen CD53 (Cd53), arachidonate 5-lipoxygenase activating protein (Alox5ap), myeloid differentiation primary response gene 88 (Myd88), macrophage scavenger receptor 1 (Msr1), matrix metallopeptidase 14 (Mmp14), triggering receptor expressed on myeloid cells 2 (Trem2) and leupaxin (Lpxn) were up-regulated in the myocardium of the MIRI group, but there was no difference inlow affinity immunoglobulin gamma Fc region receptor Ⅲ (Fcgr3), complement C1q subcomponent subunit B (C1qb) and a disintegrin and metalloproteinase domain-containing protein 8 (Adam8). By reviewing the literatures, Trem2, Lpxn, Cd53, Alox5ap, Sirpb1 and Fcgr2b were not reported to participate in MIRI.

Conclusion·This study has unearthed 6 potential hub genes for MIRI in mice, and the results can provide new ideas and entry points for further exploring the molecular mechanism and therapeutic targets of MIRI.

Key words: myocardium ischemia-reperfusion injury (MIRI), bioinformatics analysis, differentially expressed gene (DEG), hub gene, gene expression omnibus

中图分类号: