基于生物信息学分析、机器学习与实验验证揭示与二硫化物应激性死亡相关的子宫内膜异位症发病机制、生物标志物及潜在治疗药物

Unraveling pathogenesis, biomarkers and potential therapeutic agents for endometriosis associated with disulfidptosis based on bioinformatics analysis, machine learning and experiment validation.

作者信息Xiaoxuan Zhao, Yang Zhao, Yuanyuan Zhang, Qingnan Fan, Huanxiao Ke, Xiaowei Chen, Linxi Jin, Hongying Tang, Yuepeng Jiang, Jing Ma
PMID39061076
发布时间2024-07-26
DOI10.1186/s13036-024-00437-0

实验完整度

包含GEO转录组数据分析、多种机器学习建模、独立数据集验证,以及EMs小鼠模型的体内药效和基因表达验证,具有多个独立证据层级。

主要模型

GEO转录组数据集(GSE25628、GSE51981、GSE141549、GSE7305、GSE58178) EMs小鼠模型(C57BL/6小鼠,子宫内膜腹腔接种)

重点核对

机器学习模型构建时数据集按80%和20%分为训练集和测试集,采用十折交叉验证 EMs小鼠模型采用C57BL/6小鼠,每组8只,共24只,分为假手术组、EMs组和EMs+维A酸组 维A酸给药剂量为10 mg/kg,假手术组和EMs组给予等量生理盐水,建模14天后处死 qRT-PCR检测IL-1β和TNF-α的mRNA表达水平 免疫荧光染色检测CD31、VEGF、MYH10、PDLIM1和CD2AP的表达

摘要

背景:子宫内膜异位症(EMs)是一种发病机制尚未明确的疑难疾病。二硫化物应激性死亡(disulfidptosis)是一种新发现的由二硫键应激引起的程序性细胞死亡形式,有望治疗多种疾病。然而,二硫化物应激性死亡相关基因(DRGs)在EMs中的潜在作用仍不明确。本研究旨在基于生物信息学分析、机器学习和动物实验,从二硫化物应激性死亡的角度深入探索参与EMs的关键二硫化物应激性死亡基因,并探寻新的诊断标志物和候选治疗化合物。结果:对EMs在位和异位子宫内膜组织的关键模块基因和差异表达基因(DEGs)的富集分析提示,EMs与二硫化物应激性死亡密切相关。随后,我们在在位和异位子宫内膜组织中分别获得了20个和16个二硫化物应激性死亡相关DEGs。蛋白质-蛋白质相互作用(PPI)网络揭示了基因之间的复杂相互作用,并分别在在位和异位子宫内膜组织中筛选出9个和10个枢纽基因。此外,免疫浸润分析揭示了与健康对照相比,在位和异位子宫内膜组织中的免疫细胞、人类白细胞抗原(HLA)基因集和免疫检查点存在显著差异。此外,上述枢纽基因与EMs的免疫微环境密切相关。进一步地,我们应用四种机器学习算法在在位和异位子宫内膜组织中筛选特征基因,包括二元逻辑回归(BLR)、最小绝对收缩和选择算子(LASSO)、支持向量机-递归特征消除(SVM-RFE)和极端梯度提升(XGBoost)。模型训练和超参数调整在80%的数据上采用十折交叉验证方法进行,并在测试集中进行测试,通过六个指标(灵敏度、特异度、阳性预测值、阴性预测值、准确性和曲线下面积)确定了这些模型出色的诊断性能。基于机器学习最终确定了7个在位特征基因(ACTB, GYS1, IQGAP1, MYH10, NUBPL, SLC7A11, TLN1)和5个异位特征基因(CAPZB, CD2AP, MYH10, OXSM, PDLIM1)。独立验证数据集也显示了特征基因(IQGAP1, SLC7A11, CD2AP, MYH10, PDLIM1)在预测EMs方面的高准确性。此外,我们基于异位特征基因筛选出12种针对EMs的特异性化合物,并在EMs小鼠模型的异位病灶中进一步验证了维A酸对特征基因的药理作用。结论:本研究基于生物信息学分析、机器学习和动物实验证实了二硫化物应激性死亡与EMs之间的密切关联。进一步研究二硫化物应激性死亡在EMs中的生物学机制,有望为寻找潜在的诊断生物标志物和革命性的EMs治疗方法带来新的进展。

实验结论

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

研究问题
二硫化物应激性死亡相关基因在子宫内膜异位症发病机制中的作用,以及它们作为诊断标志物和治疗靶点的潜力。
核心机制
二硫化物应激性死亡相关基因通过调控肌动蛋白细胞骨架、硫代谢、免疫炎症和血管生成等通路参与EMs的发病机制。
主要证据
基于GEO数据集,通过WGCNA、差异表达分析、PPI网络、免疫浸润分析和多种机器学习算法,筛选出特征基因,并在独立数据集和小鼠模型中验证。
研究意义
为EMs的诊断和治疗提供了新的生物标志物和潜在药物靶点,特别是维A酸可能通过调节二硫化物应激性死亡相关基因发挥治疗作用。

研究路径

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

1

数据收集与预处理

获取子宫内膜异位症相关转录组数据并进行标准化处理。

从GEO数据库筛选符合条件的EMs数据集,使用R软件affy包进行RMA背景校正、log2转换、分位数归一化和中位数算法汇总。

2

WGCNA共表达网络构建

识别与在位和异位子宫内膜组织显著相关的基因模块。

使用WGCNA包构建基因共表达网络,通过样本聚类去除离群样本,建立加权邻接矩阵,识别基因模块并计算模块与表型的相关性。

3

差异表达基因分析

鉴定在位和异位子宫内膜组织中的差异表达基因。

使用limma包对训练集中的EMs样本和正常样本分别进行差异基因分析,筛选显著差异表达基因。

4

富集分析

探究关键模块基因和差异表达基因的生物学功能。

使用Metascape网站对WGCNA关键模块基因和DEGs进行GO和KEGG通路富集分析。

5

二硫化物应激性死亡相关DEGs的筛选与PPI网络分析

确定与二硫化物应激性死亡相关的DEGs并分析其相互作用。

将DRGs与DEGs取交集获得二硫化物应激性死亡相关DEGs,使用STRING数据库进行PPI网络分析,并用Cytohubba筛选枢纽基因。

6

免疫浸润景观分析

评估EMs在位和异位组织中的免疫细胞浸润情况及其与二硫化物应激性死亡相关DEGs的关联。

使用CIBERSORT算法计算22种免疫细胞的比例,比较EMs组和对照组之间的差异,并进行相关性分析。

7

机器学习诊断模型构建

筛选具有诊断价值的二硫化物应激性死亡相关特征基因并构建诊断模型。

使用BLR、LASSO、SVM-RFE和XGBoost四种算法筛选特征基因,并用十个指标评估模型性能,构建列线图预测模型。

8

候选化合物筛选

筛选针对异位特征基因的潜在治疗药物。

使用CTD、Drugbank和DGIdb数据库筛选靶向异位特征基因的化合物。

9

EMs小鼠模型验证

验证维A酸对特征基因表达和EMs病变的影响。

建立EMs小鼠模型,给予维A酸(10 mg/kg)或生理盐水处理,通过qRT-PCR和免疫荧光染色检测相关指标。

研究方法

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

产品清单

实验环节名称品牌货号
TRIzol试剂----
HiFiScript cDNA合成试剂盒CWBIOCW2569M
SYBR Green Pro Taq HS预混qPCR试剂盒AGBIOAG11701
CD31抗体signalway antibody38365
VEGF抗体signalway antibody29301
MYH10抗体abcamab230823
PDLIM1抗体signalway antibody56823
CD2AP抗体CST2135
山羊抗兔IgG H&L (Alexa Fluor® 488)abcamab150077
DAPI(4',6-二脒基-2-苯基吲哚)----
数字病理切片(荧光)扫描分析仪OLYMPUSVS120-S6-W
Image J软件National Institutes of Health--
Graphpad Prism 9GraphPad Software--
维A酸----

关键环节

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

环节核对要点
数据收集与预处理
GEO数据集的选择标准、样本分组(在位/异位EMs vs 正常对照)、数据归一化方法
阅读提示:Materials and methods: Data collection, pre-processing and normalization; Table 1
WGCNA共表达网络分析
软阈值选择、模块最小基因数、关键模块的识别标准(相关系数和P值)
阅读提示:Materials and methods: Construction of the co-expression network and hub module identification by WGCNA; Results: Identification of key modules via WGCNA
差异表达基因分析
差异基因筛选阈值(P值和log2FC)、使用的R包(limma)
阅读提示:Materials and methods: Identification of DEGs; Results: DEGs and functional enrichment analysis
二硫化物应激性死亡相关DEGs的筛选与PPI分析
二硫化物应激性死亡相关基因列表、筛选阈值、PPI网络构建参数、枢纽基因筛选算法
阅读提示:Materials and methods: Screening for disulfidptosis-related DEGs, PPI analysis and correlation analysis; Results: Identification of disulfidptosis-related DEGs
机器学习诊断模型构建
特征基因筛选算法(BLR, LASSO, SVM-RFE, XGBoost)、训练/测试集划分比例、交叉验证方法、模型性能评估指标
阅读提示:Materials and methods: Construction of disulfidptosis-related diagnosis model based on machine learning; Table 2
EMs小鼠模型构建与药物处理
小鼠品系、周龄、样本量(每组8只)、分组(假手术/EMs/EMs+维A酸)、给药途径和剂量(维A酸10 mg/kg)、建模时间(14天)
阅读提示:Materials and methods: Animals and administration schedule; Results: Murine model experiment validation
qRT-PCR检测
引物序列(TNF-α、IL-1β)、内参基因、数据分析方法(2^-△△CT)
阅读提示:Materials and methods: Quantitative real-time polymerase chain reaction (qRT-PCR)
免疫荧光染色检测
抗体浓度(1:200)、组织切片厚度(10µm)、固定条件(4% PFA 30分钟)、荧光强度定量方法(Image J)
阅读提示:Materials and methods: Immunofluorescence staining; Results: Murine model experiment validation
统计分析
统计检验方法(t检验、ANOVA、Welch检验)、方差齐性检验、显著性阈值(P<0.05)
阅读提示:Materials and methods: Statistical analysis