动脉粥样硬化肿瘤学视角:利用机器学习识别动脉粥样硬化诊断的核心基因

Athero-oncology perspective: identifying hub genes for atherosclerosis diagnosis using machine learning.

作者信息Liyan Zhao, Xuzhen Lv, Wen Chen, Xinru Li, Jie Zhou, Qi Ai, Qinhui Tuo
PMID41262249
发布时间2025-11-04
DOI10.3389/fimmu.2025.1616096

实验完整度

中

包含生物信息学分析、机器学习建模、单细胞RNA-seq分析以及免疫荧光验证,但缺乏功能实验和机制验证。

主要模型

GSE100927动脉粥样硬化组织样本 GSE159677人颈动脉斑块样本 SMC谱系追踪小鼠模型(GSE155514) ApoE-/-小鼠动脉粥样硬化模型

重点核对

诊断模型基于GSE100927训练集和GSE43292、GSE28829验证集 单细胞数据质控标准:基因数500-5000,线粒体比例<5% 小鼠模型:ApoE-/-小鼠,高脂饮食诱导12周 免疫荧光验证:人颈动脉斑块和ApoE-/-小鼠主动脉,CD68共定位

摘要

背景:平滑肌细胞(SMCs)向替代表型的转化是动脉粥样硬化发病机制中的关键过程。近期研究揭示了动脉粥样硬化与癌症之间的肿瘤学相似性,如SMCs中的DNA损伤和致癌通路激活,但确切的分子机制仍不清楚。本研究利用生物信息学整合癌症基因集,识别与动脉粥样硬化相关的关键核心基因,并探讨其免疫分子机制。方法:分析来自基因表达综合数据库(GEO)的数据集,使用Limma和WGCNA识别差异表达基因(DEGs)和模块基因。采用机器学习算法(SVM-RFE、LASSO回归和随机森林)识别用于早期动脉粥样硬化诊断的癌症相关核心基因。构建并验证了诊断模型。利用单细胞RNA测序数据的UMAP图分析核心基因的表达模式,特别关注SMC谱系追踪小鼠模型中的巨噬细胞样SMCs。在人和小鼠实验中验证了生物标志物的表达。结果:识别出四个癌症相关核心基因(CRGs):干扰素调节因子7(IRF7)、Formin同源结构域蛋白1(FHOD1)、肿瘤坏死因子(TNF)和锌指SWIM结构域蛋白3(ZSWIM3)。使用IRF7、FHOD1和TNF构建的诊断列线图在训练和验证数据集中均显示出高准确性和可靠性。免疫微环境分析显示动脉粥样硬化组和对照组之间存在显著差异。Spearman相关分析强调了核心基因与免疫细胞浸润之间的关联。单细胞RNA测序识别出不同的SMC衍生细胞簇和表型转变,在人类颈动脉斑块和小鼠模型中,可能来源于SMCs的巨噬细胞中IRF7和FHOD1的表达增加。结论:本研究整合癌症基因集识别动脉粥样硬化的关键核心基因,强调其与癌症的相似性。基于IRF7、FHOD1和TNF的诊断列线图为早期诊断提供了可靠工具,而对SMC表型转换和免疫微环境调节的见解提供了潜在的治疗靶点。

实验结论

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

研究问题
动脉粥样硬化与癌症是否存在共享的分子机制,并能否利用癌症相关基因识别动脉粥样硬化的诊断核心基因?
核心机制
平滑肌细胞向巨噬细胞样表型转化,伴随IRF7和FHOD1上调,可能涉及癌症相关通路如TNFα信号、炎症反应和干扰素应答。
主要证据
基于GEO数据集,通过WGCNA和机器学习识别出核心基因,单细胞RNA-seq显示IRF7和FHOD1在SMC源性巨噬细胞中上调,免疫荧光证实其在人和小鼠斑块巨噬细胞中高表达。
研究意义
提供“动脉粥样硬化肿瘤学”视角,为动脉粥样硬化早期诊断提供潜在生物标志物,并为靶向SMC表型转换的治疗策略提供新思路。

研究路径

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

1

差异表达基因和模块基因鉴定

识别动脉粥样硬化相关差异表达基因和共表达模块

从GEO获取GSE100927等数据集,使用limma进行差异表达分析,使用WGCNA构建共表达网络并识别临床相关模块。

2

筛选癌症相关基因

确定与癌症基因重叠的动脉粥样硬化相关基因

将WGCNA模块基因与AS差异表达基因取交集,再与癌症驱动基因(NCG 7.1数据库)取交集,获得癌症相关基因。

3

机器学习筛选核心基因

利用多种机器学习算法缩小关键基因范围

对214个CRGs应用SVM-RFE、LASSO和随机森林算法,取交集识别核心基因。

4

构建并验证诊断模型

评估核心基因对动脉粥样硬化的诊断价值

用TNF、IRF7、FHOD1构建列线图,在训练集和验证集中用ROC、PR曲线和DCA评估性能。

5

免疫浸润分析

评估核心基因与免疫细胞浸润的关系

使用ssGSEA对GSE100927样本进行28种免疫细胞浸润分析,并与CIBERSORT结果比较,计算核心基因与免疫细胞的相关性。

6

单细胞RNA-seq分析

解析核心基因在动脉粥样硬化中的细胞类型定位和动态变化

分析GSE159677(人颈动脉)和GSE155514(小鼠SMC谱系追踪)数据,UMAP聚类,分析核心基因在细胞簇中的表达。

7

免疫荧光验证

在人和小鼠动脉粥样硬化组织中验证IRF7和FHOD1表达

对人颈动脉斑块和ApoE-/-小鼠主动脉切片进行免疫荧光染色,检测IRF7和FHOD1与CD68共定位,并定量阳性细胞比例。

研究方法

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

产品清单

实验环节名称品牌货号
IRF7抗体Proteintech22392-1-AP
FHOD1抗体AbmartPC13158s
CD68抗体ElabscienceE-AB-22013
CoraLite488偶联山羊抗兔IgGAbcamab150081
Alexa Fluor 594偶联山羊抗小鼠IgGAbcamab150120
含DAPI的Fluoroshield封片剂SigmaF6057
尼康A1R HD25共聚焦显微镜Nikon--
GEOquery包----
limma包----
ggplot2包----
pheatmap包----
clusterProfiler包----
GOplot包----
glmnet包----
rms包----
pROC包----
GSVA包----
Seurat包----
GraphPad Prism V8.4.3GraphPad--

关键环节

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

环节核对要点
差异表达分析
数据集GSE100927样本量:AS 69个(实际67个),对照35个;筛选阈值:调整P<0.05,|log2FC|>0.5
阅读提示:Methods 2.2和Results 3.1
WGCNA
软阈值β=14,R²>0.9;minModuleSize=30,mergeCutHeight=0.25;样本聚类后保留67个AS和35个对照
阅读提示:Methods 2.4和Results 3.2
机器学习筛选
SVM-RFE准确率0.981;LASSO使用1-SE标准,10折交叉验证;随机森林500棵树,重要性>0.25;最终交集得到4个hub基因
阅读提示:Methods 2.5和Results 3.4
诊断模型
列线图使用TNF、IRF7、FHOD1;训练集GSE100927,验证集GSE43292和GSE28829;AUC值:训练集0.995,GSE43292 0.819,GSE28829 0.889
阅读提示:Methods 2.6和Results 3.5
免疫浸润分析
ssGSEA分析28种免疫细胞;CIBERSORT验证;Wilcoxon检验和Spearman相关
阅读提示:Methods 2.7和Results 3.7
单细胞RNA-seq分析
GSE159677质控阈值:基因数500-5000,线粒体比例<5%;GSE155514质控阈值:基因数200-4000,线粒体比例<5%;细胞聚类数为11(GSE159677)和5(GSE155514)
阅读提示:Methods 2.8和Results 3.8-3.9
免疫荧光验证
人颈动脉斑块和邻近血管(n=4);小鼠:ApoE-/-,高脂饮食(77.5%标准饲料、20%猪油、2%胆固醇、0.5%胆酸钠)12周,n=5;抗体浓度:IRF7和FHOD1 1:50,CD68 1:50,二抗1:1000
阅读提示:Methods 2.9-2.11和Results 3.10