摘要:目的利用生物信息学分析筛选扩张型心肌病(DCM)发生的相关核心基因并分析其免疫细胞浸润情况。方法从GEO数据库中下载基因芯片数据集,利用多种算法分析差异表达基因,进行功能富集,构建加权共表达网络及蛋白互作网络,鉴定核心基因,挖掘免疫细胞的浸润情况,并利用外部数据集验证核心基因的诊断效能。结果共筛选出35个差异基因和包含114个基因的核心模块。其中,F13A1、VSIG4、CD163、RNASE2及LYVE1被鉴定为核心基因。DCM组织与健康心肌中的记忆B细胞、浆细胞、Tregs细胞、活化NK细胞、M2型巨噬细胞及中性粒细胞的含量具有统计学意义。结论筛选出的5个核心基因和6种免疫细胞可能参与了DCM的发生发展。
扩张型心肌病(DCM)是一种涉及遗传等多种因素的原发性心脏病,其特点是心室腔扩张和心肌收缩能力受损,是最常见的心肌病之一[1]。目前加权基因共表达网络分析(WGCNA)是一种新兴的生物信息学技术[2],该技术通过分析基因表达谱,构建相应模块,将模块与感兴趣的临床信息关联起来。评估组织中免疫细胞浸润情况的CIBERSORT算法已成为免疫学领域常用的技术方法[3],但目前在DCM组织中暂未使用。本研究通过分析GEO数据库中相应的数据集,探讨与DCM发生相关的核心基因和免疫细胞浸润情况。
1、资料和方法
1.1 数据集来源及预处理
从GEO基因表达综合数据库[4]下载基因芯片数据集(GSE79962、GSE3585和GSE42955)。将3组数据集的DCM样本以SVA程序包进行批次矫正后合并,用于后续免疫细胞浸润分析。GSE79962数据集包括11个健康对照样本及7个DCM样本;GSE3585数据集包括5个健康对照样本及7个DCM样本;GSE4295数据集包括5个健康对照样本及12个DCM样本。所有数据集以Affy程序包依据各自平台文件,进行背景消除,去除缺失、低表达、无对应的基因探针,最终共注释到21632个基因。
1.2 筛选差异基因及加权共表达网络的构建
采用Limma程序包在GSE79962数据集中,从21632个基因中进行差异表达基因的筛选,阈值设置为|log2foldchange(FC)|≥1.5,矫正后P<0.05。采用WGCNA程序包[2]对GSE79962数据集中所有纳入样本及基因进行聚类,以合适的剪切线剔除离群样本。基因的相关矩阵由基因间的相关系数组成,根据适当的软阈值(本研究软阈值为8)将邻接矩阵转换为拓扑重叠矩阵,并将相似的基因放入同一个模块(本研究每个模块基因最小为50,切割高度为0.25)中。将临床信息与基因模块相结合,以分析基因显著性(genesignificance,GS)和模块成员(modulemembership,MM)。
1.3 DCM发生核心基因的鉴定
将核心模块中的基因与差异表达基因进行重叠,鉴定出共同基因。将共同基因输入在线STRING数据库(https://string-db.org)[5]预测蛋白质水平的互作关系,然后输入Cytoscape软件,进行PPI网络的构建。根据最大集团中心(maximalcliquecentrality,MCC)的方法选择核心枢纽基因。将上述核心基因在外部基因集中进行验证,确定核心基因在DCM组织中的表达异常,并采用ROC曲线下面积(AUC)值对核心基因的诊断效能进行判定。
1.4 功能富集分析
利用R软件中的相关程序包(Clusterprofiler、GOplot以及ggplot2)对来自核心模块基因进行基因本体论(GO)和京都基因与基因组百科全书(KEGG)分析。
1.5 免疫细胞浸润鉴定
CIBERSORT是一种使用547个基因表达值来估算组织中22个免疫细胞含量的算法,可用于肿瘤、非肿瘤组织、血液样本[6]。计算DCM组织中22个免疫细胞的比例,并与正常心肌进行对比。
1.6 统计学方法
应用R软件v3.6.3版本对转录组数据进行统计学分析。其中Fisher精确检验计算GO及KEGG富集条目值,Pearson法进行相关性分析,Wilcox检验分析免疫细胞含量差异。P<0.05为差异具有统计学意义。
2、结果
2.1 差异表达基因的鉴定结果
在GSE79962数据集中鉴定出35个差异表达基因(图1)。其中相对于正常心肌组织的显著上调基因共计21个,下调基因共14个。
2.2 加权共表达网络的构建和核心模块的识别
在GSE79962数据集中共富集出5个基因模块(图2A),蓝色模块与DCM发生的相关性最高(图2B),相关系数为0.91,该模块共包含114个基因。
2.3 DCM发生核心基因的鉴定与验证
核心模块中的114个基因与35个差异表达基因有11个交集基因(图3A),构建PPI网络(图3B),筛选出5个核心基因F13A1、VSIG4、CD163、RNASE2及LYVE1(图3C)。筛选数据集GSE79962中,核心基因在DCM心肌组织中显著低表达,外部数据集GSE42955和GSE3585同样显示核心基因在DCM心肌组织中表达下调(P<0.05),且核心基因均具有较好的诊断效能(AUC>0.7;图3D)。
2.4 基因富集分析
核心模块中的114个基因主要富集在急性炎症反应、炎症反应调节等生物学过程,并且与吞噬体、补体、凝血级联等通路有关(图4)。
2.5 免疫细胞浸润结果
DCM组织中浆细胞、Tregs细胞、中性粒细胞较健康心肌组织的含量更高(P<0.05),而记忆B细胞和M2型巨噬细胞的含量低于健康心肌组织(P<0.05;图5)。
3、讨论
DCM是全球猝死和心力衰竭的主要原因之一[6],然而,导致DCM的病因尚不完全清楚,可能包括感染、非感染性炎症、中毒、内分泌代谢紊乱、遗传、外伤等原因[7]。本文GO和KEGG富集分析结果显示,与DCM发生有关的基因模块显著富集于吞噬体、补体系统等相关通路。这提示原发性DCM发生发展过程中,免疫也参与了主要环节,相关文献[8]也揭示了使用免疫吸附方法可改善DCM的左心室功能。CIBERSORT算法结果显示DCM组织中浆细胞、Tregs细胞、活化NK细胞和中性粒细胞表现出较健康心肌组织更高的浸润模式,而记忆B细胞和M2型巨噬细胞浸润的程度较低。在一项163例DCM患者心肌组织的免疫组化研究中[9],多因素回归分析中提示CD163与胶原面积显著相关,提示M2型巨噬细胞与胶原形成之间的关联,并表明巨噬细胞向M2型分化与CD163表达可能与DCM中的心室重构有关。RNASE2基因作为一种RNA结合蛋白相关基因已参与包括胃癌[10]、肾透明细胞癌[11]的预后鉴定,该基因在所有入组的自身免疫性疾病患者中均较健康患者下调[12]。相关的生物信息学文章将该基因鉴定为M2型巨噬细胞共表达因子,并可参与调节肿瘤的免疫微环境[13]。与VSIG4基因类似,F13A1与IL-1Ra、IL-10和MMP9被鉴定为动脉粥样硬化中NOR1的潜在靶基因,并在人类替代性巨噬细胞极化时被诱导,并刺激M2表型标志物的表达[14]。
心肌组织微环境中的免疫细胞是免疫疗法的核心,在DCM病理过程中,由单核细胞/巨噬细胞、中性粒细胞和嗜酸性粒细胞介导的先天反应与适应性反应发挥最终效应。另外Th17是导致DCM过程的主要驱动因素[15]。本文结果仅CD163基因与活化的肥大细胞呈正相关性,可能由于CD163基因表达巨噬细胞特异性蛋白,提示在DCM免疫微环境中,被鉴定出的核心基因CD163与肥大细胞的增殖互为因果。
综上所述,本研究基于生物信息学方法,针对DCM患者心肌组织的测序集,进行靶向基因及免疫微环境的探究。但该研究有仍一些局限性。①必须要对核心基因进行体外、体内实验验证,来进一步解释DCM的潜在机制;②虽然CIBERSORT算法作为经典的bulk-RNA反卷积工具,但也需要进一步的实验来明确相关免疫细胞浸润的丰度。③目前临床上将DCM主要分为原发性和继发性,本文纳入的测序数据虽仅关注原发性DCM,但原发和继发并不相互排斥。
参考文献:
[7]中华医学会心血管病学分会,中国心肌炎心肌病协作组.中国扩张型心肌病诊断和治疗指南[J].临床心血管病杂志,2018,34(5):421-434.
文章来源:张清泉,吴春宇,范勐慷,潘海燕,潘闽.扩张型心肌病核心基因及其免疫浸润生物信息学分析[J].中南医学科学杂志,2022,50(02):202-205.
分享:
先天性心脏病(CHD)是由于胚胎期心血管发育障碍引起患儿出生后心脏大血管的形态、结构及功能异常的一类心脏畸形疾病的总称。根据《中国心血管健康与疾病报告2020概要》,2018年中国城市居民CHD死亡率为0.84/10万,农村居民为1.02/10万,农村地区高于城市地区。2019年全国724家开展心脏外科手术医院共进行CHD手术81 246例。
2024-04-26心动过速诱发心肌病(TIC)定义为心率异常正常化后完全或部分可逆的高和/或不规则心室率引起的心功能障碍。局灶性房性心动过速(AT)是TIC的众所周知的原因。尽管局灶性AT通常是良性的,它偶尔会表现出不间断的性质,并可能导致T-CM。在一系列接受导管消融的局灶性AT患者中,10%的患者在导管消融后左心功能得到改善。
2024-04-26目前心血管疾病已经成为全球的主要致死疾病,不仅危害公众健康,而且对个人及国家都造成了巨大的负担。快速性心律失常包括室上性心动过速、室性心动过速、心房颤动(以下简称房颤)、心室颤动等,是心血管疾病中重要的一组疾病,对患者生命健康造成极大威胁。
2024-04-24心力衰竭(heart failure,HF)[1]是多种病因导致的心脏功能损害的终末阶段,作为一个严重的、快速增长的公共卫生问题,全球约超过3 370万人患病,并且因心力衰竭症状,如呼吸困难、疲劳、运动耐量下降等,导致患者生活质量严重下降或病死率增高。根据2020中国心力衰竭医疗质量控制报告表明,心力衰竭患者平均年龄在67岁,男性占比60.8%,心衰住院患者的病死率为2.8%[2]。
2024-04-15心血管疾病是全球发病率和死亡率的主要因素,也是许多慢性病的共同终点,而心肌肥厚是心力衰竭的一个主要病因。心肌肥厚有生理性和病理性两种类型。心肌肥大最初发展为对生理和病理刺激的适应性反应,但病理性心肌肥大通常进展为心力衰竭。每种形式的肥大都受到不同细胞信号通路的调节。
2024-04-12抽取2020年5月—2021年12月本院收治的60例急性心肌梗死后心律失常患者,均经心电图、心脏彩超、心肌酶等检查确诊,患者及家属签署知情同意书,资料齐全,本研究经伦理委员会批准。排除合并其他重大疾病,精神、意识障碍,交流、沟通困难,拒绝参与研究者。随机将其分为观察组与参照组,每组各30例。
2024-04-08心脏骤停(cardiac arrest, CA)是受到广泛关注的公共卫生问题,对社会以及家庭的影响巨大。在欧洲CA的发生率接近每年每千人2次[1],在美国每年约有35万名成年人出现院外心脏骤停(out-of-hospital cardiac arrest, OHCA)[2],而我国调查显示每年发生心源性猝死的人数在54万左右。
2024-04-07多柔比星(doxorubicin,DOX)又称阿霉素,是一种具有抗肿瘤活性的蒽环类抗生素,抗癌谱较广,能抑制肿瘤细胞RNA和DNA的合成,可用于恶性淋巴瘤、乳腺癌、肺癌、膀胱癌等的临床治疗[1]。然而临床实践显示,DOX具有心脏毒性、骨髓造血功能障碍等毒副作用,且以心脏毒性最为显著,严重时会引发心力衰竭(heart failure,HF)。
2024-04-02缺血性心脏病仍是我国的主要死亡原因之一,其具有较高的疾病负担。研究显示,在无阻塞性冠状动脉狭窄的心绞痛患者中,约50%是由冠状动脉微血管功能异常所致,与心血管不良结局密切相关。
2024-03-23据统计,到2050年,全球范围内超过65岁的老年人口将超过15亿,占总人口的16%,80岁以上的人数将增加2倍。人口老龄化已经逐渐成为一个全球性的热点话题,随之而来的是衰老相关疾病发生率的增加,给社会医疗体系带来了巨大的负担。心脏在40岁之后就会开始衰老,随着年龄增加,心肌细胞发生变性和坏死,导致心脏功能下降。
2024-03-19人气:14048
人气:13229
人气:12802
人气:11630
人气:11440
我要评论
期刊名称:心电图杂志(电子版)
期刊人气:1693
主管单位:中央军委后勤保障部卫生局
主办单位:人民军医出版社
出版地方:北京
专业分类:医学
国际刊号:2095-4220
国内刊号:11-9318/R
创刊时间:2012年
发行周期:季刊
期刊开本:16开
见刊时间:7-9个月
影响因子:0.000
影响因子:0.582
影响因子:1.464
影响因子:0.000
影响因子:0.623
400-069-1609
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!