本研究通过整合生物信息学分析与实验验证,确定了2型糖尿病和动脉粥样硬化性心血管疾病共同的基因、生物学通路、免疫关联及潜在治疗靶点。
本研究通过整合生物信息学分析与实验验证,确定了2型糖尿病和动脉粥样硬化性心血管疾病共同的基因、生物学通路、免疫关联及潜在治疗靶点。
本研究探讨了2型糖尿病(T2DM)与动脉粥样硬化性心血管疾病(ASCVD)共病的潜在机制,并识别了可能的治疗靶点。从GSE78721和GSE12288数据集中提取了T2DM与ASCVD之间的共同差异表达基因(C-DEGs)。随后进行了基因本体(GO)和京都基因与基因组百科全书(KEGG)通路富集分析、蛋白质-蛋白质相互作用(PPI)网络构建、核心基因识别以及药物-基因相互作用数据库(DGIdb)分析。采用CIBERSORT方法分析了核心C-DEGs与免疫浸润细胞之间的关联。通过qRT-PCR和Western blot分析对核心C-DEGs的表达水平进行了定量检测。共鉴定出32个C-DEGs,其中包括20个上调基因和12个下调基因。C-DEGs主要富集于多个关键通路,包括病毒性心肌炎、致心律失常性右室心肌病、肥厚型心肌病和扩张型心肌病。PPI分析显示29个节点和39条边,从而在两个数据集中共同识别出八个核心C-DEGs(HSP90B1、PLAU、SLPI、TOP3A、NCF4、PRF1、TUBA1C和CS)。此外,部分核心C-DEGs(TOP3A、SLPI、NCF4、PRF1和PLAU)与免疫浸润细胞水平表现出显著相关性。针对这些核心C-DEGs的药物可能是治疗T2DM和ASCVD的潜在候选药物。同时,在T2DM合并ASCVD患者中,核心C-DEGs在mRNA和蛋白质水平的表达均得到了验证。综合生物信息学分析有助于筛选T2DM和ASCVD的候选治疗靶点、作用机制及相关药物,为这些疾病的分子治疗提供了新的见解。
糖尿病(DM)是一种病因复杂的慢性疾病,受多种遗传和环境因素的影响,目前全球约有4.63亿人患病,预计到2045年将上升至7亿人1,2。2型糖尿病(T2DM)是一种多因素疾病,其特征为胰岛素抵抗以及因β细胞功能衰竭导致的相对性胰岛素缺乏,进而引发高血糖和脂质代谢紊乱,占糖尿病病例的90%3,4,5。胰岛素抵抗是多种心血管代谢异常的关键前兆,这些异常共同促进动脉粥样硬化性心血管疾病(ASCVD)的发生,并可能影响其预后6,7。冠状动脉疾病(CAD)主要由动脉粥样硬化物质的进行性积聚引起,导致管腔狭窄和血流受阻8。动脉粥样硬化的发生由多种因素驱动,包括高胆固醇血症和糖尿病9,10。T2DM和CAD所带来的重大人类健康负担和经济成本,以及其长期管理中的挑战,凸显了制定有效预防策略的紧迫性。
合并2型糖尿病(T2DM)的冠心病(CAD)患者发生心脏事件的风险显著升高11。心血管疾病(CVD)是T2DM常见的并发症,也是T2DM患者死亡的首要原因,约有32.2%的T2DM患者受累于CVD12。LoDoCo2亚组研究证实,在慢性CAD人群中,合并T2DM与主要不良心血管事件风险增加相关13。T2DM患者发生动脉粥样硬化性心血管疾病(ASCVD)的风险升高14,15,而患有慢性CAD且合并糖尿病的个体后续发生心血管事件的风险显著更高16,凸显了对创新治疗方案存在的重要未满足临床需求。
越来越多的研究团队正在通过整合分析来自多个数据库的转录组数据,筛选2型糖尿病(T2DM)和动脉粥样硬化性心血管疾病(ASCVD)中新型差异表达基因(DEGs)。研究重点在于阐明这些疾病之间的关联,并识别有效的生物学功能。近期的生物信息学研究已探讨了T2DM与其他疾病的关系;例如,Song等人鉴定了T2DM与骨关节炎之间的共有差异表达基因17,Castillo-Velázquez等人预测了阿尔茨海默病与糖尿病(DM)之间的关联18,并通过生物信息学分析识别了T2DM与多囊卵巢综合征中的重要基因和通路19。然而,T2DM与ASCVD之间关联的分子机制仍不清楚,凸显出识别新的潜在治疗靶点以干预这两种疾病进展的重要性。
本研究旨在通过对GSE78721(T2DM)和GSE12288(ASCVD)数据集的转录组数据进行整合性生物信息学分析,系统性地识别T2DM与ASCVD之间共有的分子通路及潜在治疗靶点。选择这两个数据集是因为它们分别为每种疾病提供了匹配的患者和对照样本,并采用了可比较的微阵列平台,从而能够进行可靠的比较分析。本研究方法包括识别共同差异表达基因(C-DEGs)、进行功能富集和蛋白质相互作用分析、评估与免疫浸润的相关性以及预测相互作用的药物。此外,还进行了初步的实验验证。本研究为理解T2DM与ASCVD共病机制提供了新的见解和候选靶点。
本研究方案已通过宁夏医科大学伦理审查委员会(内部审查委员会编号:伦理[2023]-YCSKT-002)的审查与批准,并遵循《赫尔辛基宣言》中所述的原则。在患者参与研究之前,均已获得其签署的书面知情同意书。
2型糖尿病和动脉粥样硬化性心血管疾病表达数据的来源
T2DM 和 ASCVD 的基因表达数据分别来自 GEO 数据库中的 GSE78721 和 GSE12288 数据集。这两个数据集均使用了 [PrimeView] Affymetrix 人类基因表达芯片和 [HG-U133A] Affymetrix 人类基因组 U133A 芯片平台的微阵列数据。GSE78721 数据集包含来自 68 名 T2DM 患者和 62 个正常样本的表达数据(详见补充表 1),而 GSE12288 数据集包含来自 110 个 ASCVD 样本和 112 个正常对照样本的表达数据(详见补充表 2)。基因探针根据平台注释信息被转换为其对应的基因符号。
2型糖尿病和动脉粥样硬化性心血管疾病中差异表达基因的鉴定
从 GEO 数据库下载数据后,使用以下方法对 T2DM 和 ASCVD 的原始表达矩阵进行校正和标准化 "limma" 用于鉴定差异表达基因的软件包。P 值阈值为 < 0.05 被用作差异表达基因(DEG)筛选的阈值。差异表达分析结果分别利用 R 语言的 'heatmap' 和 'ggplot2' 软件包生成热图和火山图。筛选出 T2DM 与 ASCVD 之间的共有差异表达基因(C-DEGs)。 通过 维恩图
C-DEGs 的基因本体(GO)功能与京都基因与基因组百科全书(KEGG)通路富集分析
通过GO和KEGG通路富集分析确定了生物学功能和通路。使用R软件包org.Hs.eg.db将C-DEG名称转换为基因ID,并利用R软件包clusterProfiler、enrichplot、ggplot2和GOplot进行GO富集分析,以探索C-DEG的共同生物学功能。随后进行KEGG分析,以鉴定参与2型糖尿病(T2DM)与动脉粥样硬化性心血管疾病(ASCVD)共病分子机制的共有显著通路。GO和KEGG富集分析结果以条形图展示,图中包含基因数量、功能描述和P值,统计学显著性定义为P < 0.05。
C-DEGs 的蛋白质-蛋白质相互作用(PPI)网络分析
C-DEGs 和枢纽基因在本研究中被定义为共有基因。为评估这些共有基因之间的功能关联,使用 STRING 数据库(https://string-db.org)构建了蛋白质-蛋白质相互作用(PPI)网络。分析首先通过选择 "多种蛋白质" 在……之下 "搜索" 功能并输入C-DEG名称,同时指定物种为 "智人" C-DEGs的基本参数已配置用于展示PPI调控网络。使用Cytoscape软件(版本3.9.1)(https://apps.cytoscape.org/)对PPI网络进行可视化,以图形化方式呈现C-DEGs之间的相互作用。通过Cytoscape中的CytoHubba插件,基于Degree排序算法,在PPI网络中识别枢纽C-DEGs。
免疫细胞浸润分析
为了分析2型糖尿病(T2DM)组与正常对照组之间免疫细胞浸润的差异,采用CIBERSORT方法评估多种免疫浸润细胞类型的丰度。随后应用Wilcoxon检验算法比较T2DM组与对照组样本中这些浸润性免疫细胞的比例,以评估观察到的差异。数据通过小提琴图进行可视化展示。为了验证C-DEGs与免疫浸润细胞之间的相关性,使用R软件包limma、scales、tidyverse和ggplot2计算Spearman相关系数。该步骤量化了C-DEG表达水平与免疫细胞浸润程度之间的关系,结果以热图形式展示。统计学显著性以P值小于0.05为标准。
枢纽C-DEGs与潜在靶向药物的相互作用分析
药物-基因相互作用数据库(DGIdb)提供了已批准和潜在治疗药物的药物-基因相互作用及遗传特征的全面目录。在本研究中,将核心C-DEG列表输入该数据库,以获取核心C-DEG与其潜在靶向药物之间的相互作用信息,并利用Cytoscape软件进行可视化展示。
利用临床样本对核心基因进行实验验证
为实验验证核心基因,从临床受试者中采集外周血样本。研究队列包括实验组和对照组。实验组由15名被诊断为同时患有2型糖尿病(T2DM)和动脉粥样硬化性心血管疾病(ASCVD)的患者组成,代表共病状态。对照组包括10名年龄和性别相匹配、无糖尿病、心血管疾病(CVDs)或其他主要系统性疾病病史的健康志愿者。实验组所有患者均依据T2DM和ASCVD的既定临床指南确诊。
定量实时聚合酶链式反应(qRT-PCR)
使用TRIzol试剂提取总RNA,并利用逆转录系统试剂盒合成第一链cDNA。随后将cDNA稀释后用于SYBR Green PCR主混合液进行qRT-PCR,以GAPDH作为内参对照。相对表达倍数变化通过计算得出 通过 2−ΔΔCt 方法。分析所用引物序列如下:TUBA1C-上游:5′-GACCTCGTGTTGGACCGAAT-3′,下游:5′-CGAGGTGAACCCAGAACCAG-3′;NCF4-上游:5′-CTAAGTTTCAAAGCTGGAGATGTG-3′,下游:5′-GAATTCCTGAAAGGGCTCCTG-3′;PRF1-上游:5′-CCCAGTGGACACACAAAGGTT-3′,下游:5′-TCGTTGCGGATGCTACGAG-3′;SLPI-上游:5′-CCCTTCCTGGTGCTGCTT-3′,下游:5′-CCTCCTTGTTGGGTTTGG-3′;HSP90B1-上游:5′-CCGAAGAAGAACCTGAAGAG-3′,下游:5′-CATCTGTTCCCACATCCATT-3′;PLAU-上游:5′-GGGGAGATGAAGTTTGAGGTG-3′,下游:5′-GCAATGTCGTTGTGGTGAGC-3′;GAPDH-上游:5′-CAGGGCTGCTTTTAACTCTGG-3′,下游:5′-GGGTGGAATCATATTGGAACA-3′。
Western印迹检测
提取总蛋白后,用上样缓冲液稀释,并通过10% SDS-PAGE进行分离,随后转移至PVDF膜。用封闭液封闭膜后,使用以下特异性一抗进行孵育:抗TUBA1C(1:1000,ab222849)、抗HSP90B1(1:1000,ab238126)、抗NCF4(1:1000,14648-1-AP)、抗PRF1(1:1000,ab256453)、抗SLPI(1:1000,ab46763)、抗PLAU(1:1000,17968-1-AP)以及抗GAPDH(1:5000,10494-1-AP)。随后加入与辣根过氧化物酶(HRP)偶联的二抗(1:2000,ab6721),并使用增强型化学发光(ECL)底物显影蛋白条带。
统计学分析
采用 GraphPad Prism 软件进行统计分析,数据以均值 ± 标准差表示。组间差异通过学生t检验进行评估,当P值小于0.05时认为具有统计学显著性。
差异表达基因的鉴定
本研究利用了两个数据集 GSE78721 和 GSE12288 进行生物信息学分析。在 GSE78721 中,共鉴定出 2,266 个 T2DM 患者与正常样本之间的差异表达基因(DEGs),其中包括 1,209 个上调基因和 1,057 个下调基因(图1A)。在GSE12288中,共验证了ASCVD患者与对照组之间924个差异表达基因(DEGs),其中包括403个上调基因和521个下调基因(图1C)。在两个数据集中鉴定出的前50个差异表达基因(DEGs)以热图形式展示于 图1B,D通过整合这两个数据集,共鉴定出32个C-DEGs,其中包括20个上调基因和12个下调基因(图2A,B).
C-DEGs 的 GO 和 KEGG 富集分析
为了阐明与32个C-DEG相关的生物学功能及关键通路,进行了GO和KEGG分析。C-DEG中与生物过程(BP)、细胞组分(CC)和分子功能(MF)相关的前十项术语如图3A所示。生物过程主要涉及纤溶酶原激活的调控、纤溶作用的调控、心肌细胞分化的正向调控、蛋白质结合的负向调控以及心肌母细胞分化。细胞组分包括丝氨酸型内肽酶复合物、参与细胞-基质黏附的蛋白质复合物、凝缩染色体、特异性颗粒膜以及PML小体。分子功能涵盖细胞骨架的结构组分、超氧阴离子生成型NADPH氧化酶激活活性、分子内转移酶活性(磷酸转移酶)、丝氨酸型外肽酶活性以及IgG结合活性(图3A)。在KEGG分析方面,C-DEG主要富集于病毒性心肌炎、致心律失常性右室心肌病、补体与凝血级联反应、前列腺癌、肥厚型心肌病、扩张型心肌病、细胞凋亡、吞噬体以及内质网中的蛋白质加工等关键通路(图3B)。
蛋白质相互作用网络构建与核心基因筛选
在32个C-DEGs中,有三个基因在STRING数据库中未显示与其余C-DEGs存在已知相互作用,因此被排除在PPI网络构建之外。利用STRING平台对剩余的29个C-DEGs进行PPI分析,构建出包含29个节点和39条边的网络,并使用Cytoscape进行可视化(图4A随后,根据其连接度确定了八个核心C-DEG(HSP90B1、PLAU、SLPI、TOP3A、NCF4、PRF1、TUBA1C和CS)图4B).
免疫细胞浸润的差异及其与核心竞争性差异表达基因的相关性
本分析仅针对2型糖尿病数据集(GSE78721),旨在初步探讨T2DM背景下免疫微环境的变化及其与核心基因的关系。采用小提琴图可视化T2DM数据集与22种免疫细胞类型浸润水平之间的关系(图5A;详见补充表3)。免疫细胞浸润分析显示,与对照组相比,T2DM组中活化树突状细胞的比例更高(P = 0.049),而静息态树突状细胞的比例在T2DM组中显著降低(P = 0.023)。此外,利用热图可视化了8个核心C-DEG的表达水平与22种不同免疫细胞类型浸润水平之间的关系(图5B)。值得注意的是,TOP3A的表达与调节性T细胞(Tregs)的浸润水平呈强正相关(P < 0.001),SLPI的表达与γδ T细胞的浸润水平也呈强正相关(P < 0.001)。相反,NCF4的表达与活化NK细胞的浸润水平呈强负相关(P < 0.001),PRF1的表达与M2型巨噬细胞的浸润水平呈强负相关(P < 0.001),PLAU的表达与初始B细胞的浸润水平也呈强负相关(P < 0.001)。
枢纽C-DEG的候选药物靶点网络
在鉴定出的八个核心C-DEG中,CS和TOP3A未被DGIdb数据库列为已知药物靶点,因此被排除在后续候选药物-靶点网络分析之外。为筛选在2型糖尿病(T2DM)和动脉粥样硬化性心血管疾病(ASCVD)中具有潜在治疗作用的药物,将剩余的六个核心C-DEG提交至DGIdb数据库(详见补充表4)。随后,利用Cytoscape构建候选药物-靶点网络,以展示核心C-DEG与候选靶向药物之间的关联关系(图6)。该网络包含六个核心C-DEG、136种分子药物及137条相互作用边。具体而言,针对TUBA1C、NCF4、PRF1、SLPI、HSP90B1和PLAU分别鉴定出82、5、3、4、1和42种候选靶向药物,其中大多数被归类为抗肿瘤药物、抗炎药物和抗血栓药物。
枢纽C-DEGs表达的验证
由于CS和TOP3A在DGIdb数据库中未被列为候选药物靶点,且其生物学功能与T2DM和ASCVD研究背景下的治疗策略关联性较弱,因此本研究聚焦于其余六个枢纽C-DEG(HSP90B1、PLAU、SLPI、NCF4、PRF1和TUBA1C)进行实验验证。为进一步探讨这些枢纽C-DEG的作用,通过qRT-PCR和蛋白质印迹实验检测了不同组别中核心基因的表达水平(图7和图8)。与相应的健康对照组相比,T2DM合并ASCVD患者中TUBA1C、NCF4、SLPI、HSP90B1和PLAU的表达水平显著上调,而PRF1的表达则下调。这些结果提示,这些核心基因可能在T2DM合并ASCVD患者的疾病进展中发挥重要作用。
数据可用性:本研究中用于生物信息学分析的基因表达数据集 GSE78721 和 GSE12288 可从公共 NCBI GEO 数据库获取。 本研究实验验证部分生成的原始数据见补充文件 1。

图1: GSE78721和GSE12288数据集中差异表达基因(DEGs)的鉴定。(A,C)分别为GSE78721和GSE12288中差异表达基因的火山图。以|logFC|≥0.5且P < 0.05作为筛选阈值;上调基因标记为红色,下调基因标记为绿色。(B,D)分别为GSE78721和GSE12288中前50个差异表达基因的热图。请点击此处查看该图的放大版本。

图2:来自两个数据集的C-DEGs维恩图。(A)GSE78721和GSE12288数据集中T2DM与ASCVD上调C-DEGs的维恩图。(B)GSE78721和GSE12288数据集中T2DM与ASCVD下调C-DEGs的维恩图。 请点击此处查看此图的放大版本。

图 3:C-DEGs 的功能与通路富集分析。GO(A)和 KEGG(B)富集分析的条形图。纵坐标表示不同 GO 条目或 KEGG 通路的描述,横坐标表示富集的 C-DEGs 数量。请点击此处查看该图的放大版本。

图 4:C-DEG 的蛋白质-蛋白质相互作用(PPI)网络及枢纽基因筛选。(A)C-DEG 的 PPI 网络。红色和蓝色分别表示上调和下调的 C-DEG。(B)在 C-DEG 中鉴定了 8 个枢纽基因。节点代表一个 C-DEG,边代表 C-DEG 之间的相互作用。请点击此处查看该图的放大版本。

图5:免疫细胞浸润分析。(A)T2DM组与对照组之间22种免疫细胞类型的浸润差异。横轴:免疫细胞亚型;纵轴:免疫细胞的相对浸润水平。(B)8个核心C-DEG的表达与22种免疫细胞类型浸润水平的相关性热图。红色方块:正相关;蓝色方块:负相关。颜色越深,相关性越强。请点击此处查看该图的放大版本。

图6:六个枢纽C-DEG的候选药物靶点网络。红色节点表示枢纽C-DEG,绿色节点表示候选靶向药物,节点之间的连线表示相互作用对。请点击此处查看该图的放大版本。

图7:通过qRT-PCR验证T2DM合并ASCVD患者核心基因的表达。采用qRT-PCR检测HSP90B1(A)、PLAU(B)、SLPI(C)、NCF4(D)、PRF1(E)和TUBA1C(F)的mRNA表达水平。请点击此处查看该图的放大版本。

图8:通过蛋白质印迹法验证2型糖尿病合并动脉粥样硬化性心血管疾病患者核心基因的表达。采用蛋白质印迹实验检测了HSP90B1(A)、PLAU(B)、SLPI(C)、NCF4(D)、PRF1(E)和TUBA1C(F)的蛋白表达水平。请点击此处查看该图的放大版本。
补充表 1. 本研究中用于 2 型糖尿病转录组分析的 GSE78721 数据集基因表达矩阵。请点击此处下载该文件。
补充表 2. 本研究中用于 ASCVD 转录组分析的 GSE12288 数据集基因表达矩阵。请点击此处下载该文件。
补充表 3:GSE78721 数据集中的免疫细胞浸润分析。请点击此处下载该文件。
补充表 4:枢纽共同差异表达基因的候选药物-靶点相互作用。请点击此处下载该文件。
补充文件 1:本研究实验验证部分生成的原始数据。请点击此处下载该文件。
2型糖尿病(T2DM)和动脉粥样硬化性心血管疾病(ASCVD)是复杂的疾病,越来越多的证据表明两者之间存在显著关联。胰岛素抵抗是多种心血管代谢异常的关键前兆,包括T2DM,这些异常共同促进ASCVD的发生,并可能影响其预后6,7。此外,糖尿病(DM)是动脉粥样硬化进展的驱动因素20。T2DM患者发生冠状动脉疾病(CAD)的风险增加21,但关于共病潜在机制的问题仍有待阐明。本研究旨在探讨这两种疾病之间的关联及其机制,并通过生物信息学分析和实验验证识别潜在的生物标志物。
本研究共鉴定出32个T2DM与ASCVD之间的共有差异表达基因(C-DEGs),其中包括20个上调基因和12个下调基因。GO和KEGG富集分析显示,这些C-DEGs主要富集于与冠状动脉疾病(CAD)和T2DM密切相关的关键代谢通路22,23。在两种疾病中筛选出8个核心C-DEGs(HSP90B1、PLAU、SLPI、TOP3A、NCF4、PRF1、TUBA1C和CS)。值得注意的是,核心C-DEGs(TOP3A、SLPI、NCF4、PRF1和PLAU)与免疫细胞浸润水平显著相关。HSP90B1已被证实参与癌症进展24,也可能与T2DM存在关联25。PLAU已被证明在预防T2DM相关的过早血管老化中发挥作用26 ,并在动脉粥样硬化斑块进展中起关键作用,为早期预警生物标志物的发现及CAD进展机制的阐明提供了线索27。SLPI的表达水平随代谢功能障碍的进展而升高28,并已被鉴定为与糖耐量受损及T2DM相关的C-DEG29。NCF4与CAD有明确关联30,并参与调控NADPH氧化酶活性和炎性小体的激活31。由于PRF1在免疫监视和调控中的关键作用32,其具有诊断潜力,可用于区分心肌梗死与稳定型CAD33。此外,PRF1在肥胖相关胰岛素抵抗中发挥重要的免疫调节作用34。TUBA1C也与多种癌症相关,并与肿瘤微环境中的免疫浸润有关35,且已被鉴定为冠心病的易感基因36。然而,ASCVD与T2DM之间的关系仍需进一步深入研究。
KEGG 分析揭示了这些候选基因参与了诸如 "补体与凝血级联反应" 和 "内质网中的蛋白质加工" 在这些通路中对关键基因的实验验证,特别是对PLAU(一种参与纤溶和凝血的丝氨酸蛋白酶)和HSP90B1(内质网分子伴侣),证实了它们在我们的2型糖尿病模型中存在表达失调。这些发现将本研究置于特定的分子框架之中,提示内质网应激反应与溶栓通路的失衡可能构成了糖尿病病理机制的重要组成部分,超出了广泛认知的炎症过程。
大量研究已记录了糖尿病患者存在免疫反应功能障碍37。异常的免疫细胞活化在2型糖尿病(T2DM)的进展中起着关键作用38,39。免疫浸润分析显示,T2DM患者树突状细胞群体发生显著改变,表现为活化树突状细胞比例增加,而静息树突状细胞比例减少(图5A)。这表明T2DM中存在增强的抗原呈递和免疫活化状态,可能促进靶组织的慢性炎症性损伤。此外,特定枢纽基因与特定免疫细胞之间存在显著相关性(例如,NCF4与活化自然杀伤细胞,PRF1与M2型巨噬细胞)(图5B),提示代谢组织中的基因表达与局部免疫微环境组成之间存在潜在的机制联系,从而为免疫调节干预提供了新的靶点。
此外,进行了药物-靶点网络分析,以预测2型糖尿病(T2DM)和动脉粥样硬化性心血管疾病(ASCVD)的潜在治疗靶点、潜在作用机制及相关药物。该候选药物-靶点网络包含6个核心C-DEGs、136种分子药物以及137条相互作用边。分析结果显示,针对TUBA1C、NCF4、PRF1、SLPI、HSP90B1和PLAU的候选靶向药物分别为82、5、3、4、1和42种。这些药物主要归类为抗肿瘤药物、抗炎药物和抗血栓药物。值得注意的是,包括利妥昔单抗(Rituximab)、泼尼松(Prednisone)、抗氧化剂(Antioxidants)、维拉帕米(Verapamil)、抗生素(Antibiotics)、秋水仙碱(Colchicine)和肝素(Heparin)在内的部分药物已在T2DM15,40,41,42,43,44,45,46和冠状动脉疾病(CAD)15,47,48,49,50,51,52,53,54,55中显示出治疗作用,提示其在T2DM合并CAD患者中可能具有潜在治疗效果。此外,阿仑膦酸钠(Alendronate sodium)、骨化三醇(Calcitriol)、塞来昔布(Celecoxib)、马索普罗科(Masoprocol)、米非司酮(Mifepristone)和姜黄素(Curcumin)等多种药物已被证实对T2DM具有特定治疗作用56,57,58,59,60,61,62,63,这些药物也可能对ASCVD有益。同时,还对这些核心C-DEGs在T2DM和ASCVD中的表达情况进行了验证。
本研究存在若干局限性。首先,验证实验中的临床样本量不足,应进一步扩大。其次,尚需开展更深入的实验研究,以揭示2型糖尿病(T2DM)与动脉粥样硬化性心血管疾病(ASCVD)之间关联的有价值见解。未来的工作将包括利用更大的临床队列验证结果,并开展更全面的实验研究,以阐明这两种疾病之间关联的分子机制和病理通路。在实验验证部分,来自T2DM和ASCVD患者的样本被合并为一个单一的"实验组",并与健康对照组进行比较。尽管该设计有助于验证共病状态下所鉴定的核心C-DEG整体表达变化,但无法区分这些表达改变是特异性存在于T2DM、ASCVD,还是两者共有。后续研究应分别比较T2DM患者、ASCVD患者以及T2DM合并ASCVD患者的独立分组,以更准确地界定这些基因在单一疾病与共病状态中的特异性作用。此外,免疫细胞浸润分析仅在T2DM数据集(GSE78721)上进行,而未在ASCVD数据集(GSE12288)上平行开展。因此,所报道的枢纽C-DEG与免疫细胞浸润水平之间的相关性(例如TOP3A与调节性T细胞、SLPI与γδ T细胞)仅反映了它们在T2DM病理背景下的潜在免疫调节作用。另一项局限性在于,尚未对预测的药物-靶点相互作用进行功能验证。未来的研究应结合分子对接模拟以评估结合亲和力,并通过使用动物模型的in vivo实验来评估这些候选药物的治疗效果。
综上所述,本项整合生物信息学分析结合初步实验验证,鉴定了2型糖尿病与动脉粥样硬化性心血管疾病(ASCVD)之间多个共有的核心基因和通路,并提出了靶向这些基因的潜在药物候选物。这些发现为深入探究2型糖尿病与ASCVD之间复杂相互作用的机制及治疗策略,提供了新的假设和候选靶点。
作者声明无竞争利益。
作者衷心感谢宁夏自然科学基金项目(编号:2023AAC03498)和银川市科技创新项目(编号:2024SF005)提供的经费支持。作者同时感谢所有参与本研究的受试者和研究人员,正是他们的贡献使本研究得以顺利完成。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 抗GAPDH | Proteintech | 60004-1-Ig | |
| 抗HSP90B1 | Abcam | ab238126 | |
| 抗NCF4 | Proteintech | 14648-1-AP | |
| 抗PLAU | Proteintech | 17968-1-AP | |
| 抗PRF1 | Abcam | ab256453 | |
| 抗SLPI | Abcam | ab46763 | |
| 抗TUBA1C | Abcam | ab222849 | |
| CIBERSORT | 斯坦福大学 | https://cibersortx.stanford.edu | |
| clusterProfiler (R软件包) | Bioconductor | http://www.bioconductor.org/packages/release/bioc/html/clusterProfiler.html | |
| CytoHubba (Cytoscape插件) | Cytoscape联盟 | https://apps.cytoscape.org | |
| Cytoscape (版本 3.9.1) | Cytoscape联盟 | http://cytoscape.org/ | |
| DGIdb (药物-基因相互作用数据库) | 华盛顿大学 | http://www.dgidb.org/ | |
| 增强型化学发光(ECL)底物 | 赛默飞世尔科技 | 34080 | |
| enrichplot (R软件包) | Bioconductor | http://bioconductor.org/packages/release/bioc/html/enrichplot.html | |
| 凝胶成像系统 | 伯乐 | ChemiDox XRS | |
| GEO数据库 (数据集: GSE78721, GSE12288) | NCBI | https://www.ncbi.nlm.nih.gov/geo/ | |
| ggplot2 (R软件包) | CRAN | https://ggplot2.tidyverse.org | |
| GOplot (R软件包) | CRAN | https://cran.r-project.org/web/packages/GOplot/index.html | |
| GraphPad Prism(版本 9.1.2) | GraphPad软件公司 | https://www.graphpad.com/scientific-software/prism/ | |
| limma (R软件包) | Bioconductor | http://bioconductor.org/packages/release/bioc/html/limma.html | |
| 上样缓冲液 | 伯乐 | 1610737 | |
| org.Hs.eg.db (R软件包) | Bioconductor | https://bioconductor.org/packages/release/data/annotation/html/org.Hs.eg.db.htm | |
| PVDF膜 | 伯乐 | 1620177 | |
| R软件 (版本 4.0.2) | R基金会 | https://www.r-project.org/ | |
| 实时PCR系统 | 伯乐 | CFX-96 | |
| 逆转录系统试剂盒 | 宝生物工程 | D6110A | |
| scales (R软件包) | CRAN | https://CRAN.R-project.org/package = scales | |
| SDS-PAGE凝胶 | 伯乐 | 4561085 | |
| 二抗 | Abcam | ab6721 | |
| STRING数据库 | STRING联盟 | https://string-db.org/ | |
| SYBR Green PCR主混合液 | 诺唯赞生物科技 | 货号 #Q311-02 | |
| tidyverse (R软件包) | CRAN | https://cran.r-project.org/web/packages/tidyverse/index.html | |
| Trans-Blot转膜系统 | 伯乐 | 1704150 | |
| TRIzol试剂 | 赛默飞世尔科技 | #1559602 |