2型糖尿病中线粒体相关关键基因的鉴定及知母-黄柏药对作用机制的阐明:生物信息学、机器学习与实验验证的整合方法

Identification of mitochondria-related key genes in type 2 diabetes mellitus and elucidation of the Zhimu-Huangbai herb Pair's mechanism: an integrated approach of bioinformatics, machine learning, and experimental validation.

作者信息Sufang Zhang, Minghe Gu, Baoting Chen, Yiming Liu, Aihua Lin
PMID41909133
发布时间2026-03-13
DOI10.3389/fcell.2026.1763178

实验完整度

包含生物信息学分析、机器学习建模、分子对接及体外细胞实验等多层级验证,并通过Western blot及AMPK抑制剂验证机制。

主要模型

INS-1细胞高糖高脂(HGHL)损伤模型 GSE76894、GSE25724、GSE38642数据集 人胰岛β细胞及胰岛组织样本

重点核对

GSE76894作为训练集,GSE25724和GSE38642合并作为验证集 HGHL模型:33.3 mmol/L葡萄糖+600 μmol/L棕榈酸处理24小时 ZB剂量:50、100、200 μg/mL;二甲双胍1 mmol/L 关键蛋白:BCAT2、CASP8、EPHX2、UCP2、p-AMPK/AMPK、SIRT1、PGC-1α AMPK抑制剂Compound C(10 μM)预处理1小时

摘要

背景:线粒体损伤通过损害细胞能量代谢和胰岛素敏感性,在2型糖尿病(T2DM)发病机制中发挥关键作用。知母-黄柏药对(ZB)是由知母和黄柏组成的经典中药配方,对T2DM有疗效,但其分子机制尚不清楚。本研究旨在鉴定2型糖尿病的关键线粒体相关基因及ZB的潜在机制。方法:从GEO数据库获取T2DM基因表达数据集(GSE76894、GSE25724和GSE38642)。通过筛选多个数据库(包括TCMSP和HERB)鉴定ZB药对和T2DM的交集靶点。从人类线粒体相关数据库获取线粒体功能相关基因。采用WGCNA鉴定差异表达基因,然后与生物活性成分-靶基因和线粒体相关基因取交集,构建PPI网络。随后进行GO和KEGG富集分析。应用四种机器学习算法——SVM-RFE、RF、GLM和XGB——筛选特征基因并建立诊断模型。此外,分析特征靶点与免疫细胞浸润的相关性,进行单基因GSEA,并进行分子对接以研究特征靶点与ZB生物活性成分之间的相互作用。实验验证方面,将INS-1细胞分为六组:对照组、模型组、二甲双胍组以及ZB低、中、高剂量组。评估细胞活力、凋亡、ROS水平、线粒体膜电位以及线粒体形态和功能。进行Western blot分析以评估线粒体相关基因(BCAT2、CASP8、EPHX2和UCP2)和AMPK–SIRT1–PGC-1α信号通路组分的表达。结果:共鉴定出8个与ZB治疗T2DM相关的线粒体相关差异表达基因。GO分析显示在多个生物学过程中富集,包括对营养水平的反应;细胞组分,如孔复合物;分子功能,包括有毒物质结合。KEGG通路分析表明在凋亡、p53信号通路和坏死性凋亡等通路中显著富集。通过机器学习算法筛选出三个关键基因——BCAT2、CASP8和EPHX2,构建的T2DM诊断模型均显示出大于0.7的曲线下面积(AUC)值,表明具有令人满意的判别性能。免疫浸润分析显示,所有三个关键基因均与免疫细胞群体显著相关。分子对接结果表明,三个关键基因与其对应的ZB生物活性化合物表现出强结合亲和力(≤−5.0 kcal/mol),但CASP8-烟酰胺组合除外。实验验证显示,ZB显著增强了高糖高脂条件下INS-1细胞的活力,抑制了凋亡,减少了细胞内ROS生成,并改善了线粒体膜电位、线粒体形态和呼吸功能。同时,UCP2和BCAT2的蛋白表达水平显著上调,而CASP8和EPHX2显著下调。此外,ZB处理上调了p-AMPK/AMPK比值以及SIRT1和PGC-1α的表达。结论:以BCAT2、CASP8和EPHX2为特征的诊断模型为T2DM的诊断和治疗提供了新的见解。ZB的治疗机制涉及调节线粒体相关基因(BCAT2、CASP8、EPHX2、UCP2)和激活AMPK-SIRT1-PGC-1α通路,从而改善线粒体形态和功能,减少氧化损伤,增强能量代谢。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
T2DM中哪些线粒体相关基因是关键靶点,以及知母-黄柏药对通过何种机制改善T2DM?
核心机制
ZB通过调节线粒体相关基因(BCAT2、CASP8、EPHX2、UCP2)表达,并激活AMPK-SIRT1-PGC-1α信号通路,改善线粒体形态和功能,减少氧化损伤,增强能量代谢。
主要证据
基于GEO数据集与机器学习筛选出BCAT2、CASP8、EPHX2三个诊断基因;体外实验证实ZB在高糖高脂诱导的INS-1细胞中增强细胞活力、抑制凋亡、降低ROS、恢复线粒体膜电位、形态和呼吸功能,并调控相关蛋白表达。
研究意义
为T2DM的诊断和治疗提供新见解,为知母-黄柏药对治疗T2DM的机制提供实验依据,并支持中医“滋阴清热”理论的现代解释。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

T2DM差异表达基因与模块鉴定

鉴定T2DM相关差异表达基因及关键共表达模块

对GSE76894数据集用limma进行差异表达分析(P<0.05, |log2FC|>0.263),并用WGCNA构建共表达网络,确定与T2DM显著相关的模块(turquoise和brown)。

2

核心靶点筛选与富集分析

筛选ZB治疗T2DM的线粒体相关核心靶点并分析其功能

取DEGs、ZB靶点、WGCNA模块基因和线粒体基因的交集,获得8个靶点,构建PPI网络,并进行GO和KEGG富集分析。

3

机器学习筛选特征基因

筛选具有诊断价值的线粒体相关特征基因

基于8个靶点,使用SVM-RFE、RF、GLM和XGB四种算法建立预测模型,评估模型性能并选择重要变量。

4

特征基因表达验证与诊断模型评估

验证候选基因在独立数据集中的表达差异并评估诊断价值

在训练集和验证集中比较候选基因在T2DM与对照组间的表达,并进行ROC分析。

5

免疫浸润与GSEA分析

分析特征基因与免疫细胞浸润的关系及其相关信号通路

使用ssGSEA评估免疫细胞浸润,并分析特征基因与免疫细胞的相关性;使用单基因GSEA分析特征基因富集的通路。

6

分子对接预测活性成分与靶点结合

预测ZB活性成分与特征基因编码蛋白的结合亲和力

使用AutoDock Vina对3个靶蛋白和20个ZB化合物进行分子对接,计算结合自由能。

7

体外实验验证ZB的保护作用

验证ZB对高糖高脂诱导INS-1细胞损伤的保护作用

建立HGHL诱导的INS-1细胞损伤模型,评估细胞活力(CCK-8)、凋亡(流式)、ROS(DCFH-DA)、线粒体膜电位(JC-1)、线粒体形态(MitoTracker Deep Red)和呼吸功能(Seahorse)。

8

Western blot验证蛋白表达与通路机制

验证ZB对线粒体相关蛋白及AMPK-SIRT1-PGC-1α通路的影响

通过Western blot检测BCAT2、CASP8、EPHX2、UCP2、AMPK、p-AMPK、SIRT1、PGC-1α表达,并使用AMPK抑制剂Compound C验证通路依赖性。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
RPMI 1640培养基----
胎牛血清----
磷酸盐缓冲液Wuhan Pricella Biotechnology Co., Ltd.--
0.25%胰蛋白酶Wuhan Pricella Biotechnology Co., Ltd.--
D-葡萄糖溶液Wuhan Pricella Biotechnology Co., Ltd.--
PA-BSA高脂细胞添加剂Xi'an Kunchuang Technology Development Co., Ltd.--
CCK-8试剂盒Shanghai Beyotime Biotechnology Co., Ltd.--
ROS检测试剂盒Shanghai Beyotime Biotechnology Co., Ltd.--
MitoTracker深红荧光探针Shanghai Beyotime Biotechnology Co., Ltd.--
JC-1线粒体膜电位检测试剂盒Beijing Solarbio Science & Technology Co., Ltd.--
Seahorse XF细胞线粒体压力测试试剂盒Guangzhou Huanghe Instrument Technology Co., Ltd.--
BCAT2抗体Abcam--
CASP8抗体Abcam--
EPHX2抗体Abcam--
UCP2抗体Abcam--
AMPK抗体Abcam--
p-AMPK抗体Abcam--
SIRT1抗体Abcam--
PGC-1α抗体Abcam--
化合物CMCE--
知母干燥根茎Guangdong Kangmei Pharmaceutical Co., Ltd.--
黄柏干燥树皮Guangdong Kangmei Pharmaceutical Co., Ltd.--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
数据来源与处理
数据集GSE76894、GSE25724、GSE38642的样本类型和数量;差异表达阈值;WGCNA软阈值与模块筛选标准
阅读提示:Methods 2.1.1, 2.1.5, 2.1.6
机器学习和特征筛选
四种机器学习算法的具体参数;特征选择方法;验证集批次效应校正
阅读提示:Methods 2.1.9, 2.1.10
细胞培养与HGHL模型
INS-1细胞培养条件;HGHL浓度(33.3 mmol/L葡萄糖+600 μmol/L棕榈酸)和处理时间(24h)
阅读提示:Methods 2.2.3, Results 3.11
药物处理
ZB提取物浓度(50, 100, 200 μg/mL)和二甲双胍浓度(1 mmol/L);处理时间24h;细胞分组
阅读提示:Methods 2.2.5
功能检测
CCK-8、流式凋亡、ROS、JC-1、MitoTracker、Seahorse实验的具体操作条件(如探针浓度、孵育时间、仪器设置)
阅读提示:Methods 2.2.6-2.2.10
蛋白表达验证
Western blot一抗稀释比、二抗浓度、蛋白上样量、曝光条件;Compound C浓度(10 μM)及预处理时间
阅读提示:Methods 2.2.11