基于双硫死亡的心肌缺血再灌注损伤相关风险基因筛选及分析

目的:筛选基于双硫死亡的心肌缺血再灌注损伤(myocardial ischemia reperfusion injury,MIRI)相关风险基因及探索其可能的作用途径。方法:用R语言limma分析并筛选MIRI相关风险基因24 h组与Control组之间的差异表达基因。基于GO和KEGG数据库利用R语言cluterProfiler包对其进行功能富集分析。根据双硫死亡基因集以及GSE160516的表达数据,利用R语言GSVA包中的ssgsea法对其进行富集分析,并利用Venn图获得心肌缺血再灌注中双硫死亡相关基因的表达。用R语言psych包中的corr.test计算双硫死亡相关基因与凋亡因子、线粒体、铁死亡及炎症等相关基因的相关性,并绘制热图。利用R语言Mfuzz包来进行不同时间点心肌缺血再灌注转录组数据的表达模式聚类分析并利用Venn图获取时序性差异表达的双硫死亡基因。结果:本研究共筛选出17个差异表达双硫死亡相关基因,经过双硫死亡相关基因与炎症、凋亡、铁死亡、线粒体基因相关性分析,以及双硫死亡相关基因与心肌缺血再灌注时序性分析的差异表达基因筛选,最终得到双硫死亡相关基因Flna、Myl6、Tln1在MIRI不同时间的特异性表达。结论:Flna、Myl6、Tln1在MIRI中可能起着关键作用,其基于双硫死亡同时与线粒体相关基因密切相关,可作为MIRI患者风险因素的筛选指标。