整合微阵列数据与单细胞测序分析以探索与心力衰竭中巨噬细胞浸润相关的关键基因

Integration of Microarray Data and Single-Cell Sequencing Analysis to Explore Key Genes Associated with Macrophage Infiltration in Heart Failure.

作者信息Jin Rao, Xuefu Wang, Zhinong Wang
PMID39717663
发布时间2024-12-19
DOI10.2147/JIR.S475633

实验完整度

结合公开转录组数据挖掘、四种机器学习建模、无监督聚类、WGCNA、细胞通讯分析,并通过TAC小鼠模型进行qPCR和Western blot体内验证,证据层级多样。

主要模型

GSE222144单细胞数据集(正常和HF样本) GSE57338批量RNA-seq数据集(313例有或无HF个体) GSE26887外部验证数据集 TAC诱导的C57BL/6小鼠HF模型

重点核对

小鼠TAC模型:雄性C57BL/6,27G针头缩窄主动脉弓,术后4周 SERPINA3在HF小鼠心肌中表达下调(qPCR和Western blot验证) 诊断模型基于SERPINA3、GPAT3、ANPEP、FCER1G,AUC评估预测能力 GSE57338数据集用于差异表达分析(调整p<0.05,|logFC|>0.5) WGCNA软阈值β=6,排除异常样本GSM1379815和GSM1380018

摘要

背景:心脏巨噬细胞是一种具有高度可塑性和适应性的异质群体,其在心力衰竭(HF)中的作用机制仍不完全清楚。方法:我们利用单细胞和批量RNA测序数据揭示非心肌细胞的异质性,并评估每个亚群的免疫反应性。此外,我们采用四种集成的机器学习算法来识别具有诊断价值的巨噬细胞相关基因,并进行了体内验证。为了评估HF中的免疫浸润特征,我们使用了CIBERSORT和单样本基因集富集分析(ssGSEA)。应用无监督共识聚类算法来识别巨噬细胞相关的HF亚型。此外,采用scMetabolism探讨巨噬细胞亚型的特定代谢模式。最后,使用CellChat研究已识别亚型之间的细胞间相互作用。结果:HF中巨噬细胞的免疫反应性评分高于其他细胞类型。巨噬细胞簇的GSEA显示白细胞介导的免疫过程、抗原加工和呈递显著富集。机器学习结果的交集显示SERPINA3、GPAT3、ANPEP和FCER1G可作为特征基因,并形成具有良好预测能力的诊断模型。无监督共识聚类算法揭示了巨噬细胞的免疫和代谢亚型。巨噬细胞亚群的代谢异质性可导致巨噬细胞极化为不同类型,这可能与糖酵解和线粒体氧化磷酸化之间的代谢重编程有关。细胞通讯揭示巨噬细胞与中性粒细胞形成相互作用网络,以支持彼此的功能和维持。复杂的传出和传入信号与心肌纤维化密切相关。结论:SERPINA3、GPAT3、ANPEP和FCER1G可能作为免疫治疗靶点和核心生物标志物。巨噬细胞的免疫和代谢异质性可能为探索HF的潜在机制和新型免疫疗法提供更精确的方向。

实验结论

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

研究问题
心力衰竭中巨噬细胞浸润相关的关键基因及其免疫代谢异质性机制是什么?
核心机制
巨噬细胞亚群之间的代谢重编程(如糖酵解与氧化磷酸化)驱动巨噬细胞极化,影响HF的免疫炎症状态;ANXA1-FPR1/FPR2信号介导巨噬细胞-中性粒细胞通讯,与心肌纤维化相关。
主要证据
单细胞测序显示HF中巨噬细胞比例增加且免疫评分最高;机器学习鉴定出SERPINA3、GPAT3、ANPEP、FCER1G四个诊断基因,并在外部数据集验证;TAC小鼠模型中SERPINA3表达下调;无监督聚类区分出免疫型(Cluster B)和代谢型(Cluster A)两种HF亚型。
研究意义
SERPINA3、GPAT3、ANPEP和FCER1G可作为HF的潜在生物标志物和免疫治疗靶点,巨噬细胞的免疫代谢异质性为HF的精准治疗提供新方向。

研究路径

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

1

单细胞转录组图谱构建

揭示HF中非心肌细胞的异质性,特别是巨噬细胞的组成和状态。

对GSE222144的scRNA-seq数据进行质量控制(排除低质量细胞),使用Harmony整合去批次,PCA降维,基于marker基因注释细胞类型,共得到9种细胞系。

2

免疫活性评分与细胞通讯分析

评估不同细胞类型的免疫活性,并探究巨噬细胞与其他细胞的相互作用。

使用AUCell计算免疫相关基因集的AUC评分,比较各细胞类型的免疫活性;利用CellChat分析细胞间配体-受体相互作用,重点关注巨噬细胞-中性粒细胞通讯。

3

巨噬细胞相关诊断模型的构建与验证

筛选具有诊断价值的巨噬细胞相关基因,并构建预测模型。

对GSE57338进行差异表达分析,与巨噬细胞差异基因取交集;采用LightGBM、XGBoost、RF和SVM四种机器学习算法筛选特征基因,取交集得到SERPINA3、GPAT3、ANPEP、FCER1G;构建nomogram,通过ROC曲线评估诊断性能,并在GSE26887中验证。

4

TAC小鼠HF模型建立与关键基因验证

验证SERPINA3在HF中的表达变化。

对雄性C57BL/6小鼠进行TAC手术(27G针头),术后4周通过超声心动图评估心功能,HE和Masson染色观察组织病理,qPCR和Western blot检测SERPINA3表达。

5

无监督聚类识别HF亚型

基于巨噬细胞相关基因表达模式区分HF的不同分子亚型。

使用ConsensusClusterPlus对177个HF样本基于27个巨噬细胞相关基因进行共识聚类,确定最佳聚类数K=2,并比较两个亚型在通路富集(GSVA)和免疫浸润(ssGSEA)上的差异。

6

WGCNA识别代谢相关hub基因

识别与脂肪酸代谢和氧化磷酸化相关的关键基因模块。

对177个HF样本进行WGCNA,选择软阈值β=6,构建共表达网络,鉴定与脂肪代谢和氧化磷酸化显著相关的模块(MEgreen),并与差异表达基因取交集,构建PPI网络,使用cytoHubba筛选hub基因。

7

巨噬细胞亚型代谢景观分析

在单细胞分辨率下揭示巨噬细胞亚型的代谢异质性。

对巨噬细胞进行更高分辨率的聚类,识别出四种亚型(Mac_M1, Mac_M2, Mac_Fib, Mac_TREM2),使用scMetabolism和KEGG代谢基因集评估各亚型的代谢通路活性。

研究方法

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

产品清单

实验环节名称品牌货号
便携式小动物麻醉机Shanghai YuyanACS
异氟烷Shanghai Yuyan--
SERPINA3抗体Biossbs-0094R
预染蛋白markerThermo Scientific26616

关键环节

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

环节核对要点
单细胞数据分析
质控阈值(基因数<300,线粒体基因>25%),Harmony整合参数,细胞类型注释marker基因
阅读提示:Materials and Methods - scRNA-Seq Data Processing
诊断模型构建
GSE57338差异表达阈值(adjusted p<0.05,|logFC|>0.5),机器学习算法具体参数,ROC曲线评估方法
阅读提示:Materials and Methods - Establishment and Validation of a Diagnostic Model
TAC小鼠模型
小鼠品系(雄性C57BL/6)、年龄、体重、手术细节(27G针头,6-0缝线)、麻醉方式、术后护理、超声心动图测量指标、样本量(n=3)
阅读提示:Materials and Methods - Construction and Evaluation of Mouse Models of Heart Failure
qPCR验证
引物序列、内参基因(Actin)、RNA提取方法、逆转录和PCR程序
阅读提示:Materials and Methods - Experimental Validation of Macrophage-Related Characterized Genes
Western blot验证
抗体浓度、蛋白上样量、抗体孵育条件、显影方法
阅读提示:Materials and Methods - Experimental Validation of Macrophage-Related Characterized Genes
WGCNA
软阈值β=6,异常样本排除(GSM1379815, GSM1380018),模块识别方法
阅读提示:Materials and Methods - Identification of Metabolism-Related Hub Genes in Heart Failure Based on WGCNA