缺氧驱动的脓毒症巨噬细胞异质性：单细胞转录组揭示免疫重编程与治疗靶向

《Frontiers in Immunology》：Hypoxia-driven macrophage heterogeneity in sepsis: single-cell transcriptomic insights into immune reprogramming and therapeutic targeting

　　

脓毒症以深度免疫失调和代谢重塑为特征，然而缺氧如何塑造巨噬细胞状态多样性仍未完全明确。本研究整合早期与晚期脓毒症及健康对照的单细胞转录组图谱，构建了一个以缺氧为中心的免疫图谱。单细胞景观识别出23个簇共103,851个细胞，缺氧评分显示巨噬细胞表现出最强的缺氧活性，尤其在脓毒症早期。巨噬细胞被进一步分为低、中、高缺氧状态，揭示了一个连续转录梯度而非离散亚群。高缺氧巨噬细胞表现出CytoTRACE定义的增强可塑性，占据拟时间中间至晚期状态，并显示炎症与代谢通路激活增强。细胞间通讯分析将这些巨噬细胞定位为双向信号枢纽，具有放大的MIF、CCL、ANNEXIN和IL1相关相互作用。整合高维加权基因共表达网络分析（hdWGCNA）与机器学习分析识别出缺氧相关枢纽基因，其中TIPARP作为与巨噬细胞缺氧重塑相关的关键调控因子脱颖而出。TIPARP在高缺氧巨噬细胞中优先富集，并经实验验证对缺氧和炎症刺激有响应，其敲低可增强HIF1A和VEGF的诱导。药物敏感性预测进一步提示缺氧适应巨噬细胞具有状态特异性脆弱性。综上，这些发现将缺氧定义为脓毒症中巨噬细胞异质性的核心决定因素，并提名TIPARP作为免疫代谢干预的候选靶点。

论文解读

研究背景与问题

脓毒症是由感染引起的致命综合征，以宿主免疫反应失衡和微环境显著重塑为特征。尽管诊断和治疗手段不断进步，脓毒症的高死亡率仍未得到根本改善。越来越多的证据表明，免疫细胞的功能异质性和动态演变是疾病进程和预后的关键决定因素，但传统群体水平转录组研究无法阐明复杂免疫细胞网络内部的精细调控机制，限制了对脓毒症潜在发病机制的深入理解。单细胞RNA测序（scRNA-seq）技术的出现为揭示不同疾病阶段细胞类型特异性转录特征及其动态变化提供了重要机遇。在炎症微环境中，缺氧被认为是驱动免疫功能重编程的关键因素之一，不仅调控能量代谢，还通过多种信号通路影响炎症反应的激活、分化和放大。然而，在脓毒症高度动态的病理过程中，缺氧在不同免疫细胞中的作用模式及其对细胞命运决定的影响仍缺乏系统性研究。巨噬细胞作为关键效应细胞，其在炎症反应、抗原呈递和组织修复中的核心作用已被广泛认可，但缺氧驱动的巨噬细胞功能状态层级特征及分子调控网络尚未完全阐明。基于此，本研究旨在整合多队列单细胞转录组数据，系统构建脓毒症细胞图谱，揭示缺氧驱动的巨噬细胞异质性及其在脓毒症免疫微环境中的核心作用。

研究内容与结论

研究人员整合了来自GEO数据库的两个脓毒症单细胞数据集（GSE167363和GSE175453）及其健康对照，经批次效应校正后共纳入103,851个细胞，通过无监督聚类识别出23个细胞簇，涵盖NK细胞、B细胞、巨噬细胞、中性粒细胞、内皮细胞和上皮细胞等主要细胞类型。利用396个缺氧相关基因，通过五种评分方法（AUCell、UCell、singscore、ssGSEA和AddModuleScore）标准化后构建复合缺氧评分，发现巨噬细胞在脓毒症中表现出显著升高的缺氧活性，尤其在早期脓毒症中最为突出。依据复合评分将巨噬细胞分为低缺氧（LHM）、中缺氧（IHM）和高缺氧（HHM）三组，揭示了一个连续转录梯度而非离散亚群。高缺氧巨噬细胞表现出增强的CytoTRACE可塑性，占据拟时间中间至晚期状态，并显示炎症和代谢通路激活增强。细胞间通讯分析表明高缺氧巨噬细胞作为双向信号枢纽，具有放大的MIF、CCL、ANNEXIN和IL1相关相互作用。通过hdWGCNA和多种机器学习算法筛选出27个共同枢纽基因，其中TIPARP被确定为核心调控因子。实验验证表明，TIPARP在缺氧和炎症刺激下表达上调，其敲低可进一步增强HIF1A和VEGF的诱导，提示TIPARP在缺氧反应中发挥负反馈调控功能。药物敏感性预测进一步提示高缺氧巨噬细胞具有状态特异性治疗脆弱性。该研究将缺氧定义为脓毒症中巨噬细胞异质性的核心决定因素，为免疫代谢干预提供了新靶点。

关键技术方法

研究人员主要运用了以下关键技术：整合多队列单细胞转录组数据分析，包括来自GSE167363（脓毒症发病0-6小时）和GSE175453（脓毒症发病14-21天）两个数据集及其健康对照组，经Seurat锚点整合流程校正批次效应；利用AUCell、UCell、singscore、ssGSEA和AddModuleScore五种算法构建复合缺氧评分系统；采用Monocle3进行拟时间轨迹分析；使用CellChat进行细胞间通讯分析；应用高维加权基因共表达网络分析（hdWGCNA）识别与高缺氧巨噬细胞表型相关的共表达模块；结合LASSO、随机森林（RF）、决策树（DT）、梯度提升机（GBM）和ABESS五种机器学习算法筛选枢纽基因；通过mlr3框架评估八种机器学习模型的预测性能；利用RAW264.7巨噬细胞进行Western blot和qRT-PCR实验验证。

研究结果

脓毒症单细胞景观：通过整合单细胞测序数据，研究人员评估了正常样本与脓毒症样本的差异，共纳入103,851个细胞。UMAP图显示细胞聚为23个簇，利用已知标记基因注释出NK细胞、B细胞、巨噬细胞、中性粒细胞、内皮细胞和上皮细胞等细胞类型。细胞组成比例分析揭示了样本特异性免疫微环境的潜在差异，标记基因热图和特征图进一步验证了细胞类型注释的准确性。

评估各细胞类型缺氧活性：通过多种缺氧评分方法，研究人员发现脓毒症样本的缺氧活性显著高于健康样本。巨噬细胞在脓毒症中表现出相对较高的缺氧相关基因表达，五种评分方法结果一致。早期脓毒症样本的缺氧活性显著高于晚期脓毒症和健康样本，UMAP可视化进一步证实了这一空间异质性分布。

鉴定脓毒症巨噬细胞缺氧活性：基于缺氧评分分布，巨噬细胞被分为LHM、IHM和HHM三组，反映缺氧活性梯度。UMAP显示HHM细胞主要分布于中右区域，LHM集中于左侧和外周区域，IHM广泛分布于两者之间，呈现连续转变。CytoTRACE分析显示缺氧评分与细胞发育潜能呈正相关（Pearson r = 0.55），高缺氧细胞具有更高转录复杂性和代谢/炎症激活状态。单细胞轨迹分析显示HHM主要分布于轨迹中后期，LHM聚集于起始端，形成三个主要分支，反映缺氧状态对巨噬细胞分化轨迹的显著影响。

脓毒症各细胞间通讯分析：细胞间通讯网络显示HHM、单核细胞和CD4+ T细胞具有最大节点，HHM细胞与Cycling、cDC、Mast和pDC细胞通讯最频繁。HHM的信号输出和输入强度均显著高于LHM，在MIF、GALECTIN、RESISTIN、ANNEXIN、CCL、VISFATIN和IL1等多条信号通路中表现出强大的输出能力。HHM巨噬细胞在细胞通讯网络中兼具信号发送和接收枢纽功能，主要通过ANXA1–FPR1/FPR2、CCL3–CCR1、IL1B–IL1R2等通路放大炎症信号。GSVA分析显示HHM组缺氧和TNFA_SIGNALING_VIA_NFKB通路显著上调。

hdWGCNA鉴定与HHM相关的共表达基因模块：研究人员构建了无标度共表达网络，选择最优软阈值为6，识别出八个基因共表达模块，其中pink、red、brown、blue和yellow五个模块在HHM中高度激活。保留模块特征基因连接度（kME）≥0.30的基因作为模块相关枢纽基因，并构建了相应的蛋白质-蛋白质相互作用网络。

利用机器学习算法进一步筛选枢纽HRGs：HHM与LHM差异表达分析鉴定出417个上调基因，与hdWGCNA五个模块中492个基因取交集获得222个候选缺氧相关基因（HRGs）。Pearson相关分析识别出192个与缺氧活性显著正相关的基因。疾病本体论（DO）富集分析显示与肝炎显著相关，基因本体论（GO）分析显示与中性粒细胞活化和中性粒细胞介导免疫密切相关。通过LASSO（169个基因）、随机森林（192个）、决策树（188个）、GBM（146个）和ABESS（50个）五种机器学习算法交叉分析，最终确定27个共同枢纽基因。

机器学习基准模型评估：AUC诊断模型显示候选基因具有中等偏高的判别能力，LDHA、ETS2、FKBP5、ENO1和CD63排名靠前。训练集和测试集的PR曲线和ROC曲线均显示良好判别效能，决策曲线分析（DCA）证明模型具有临床实用性。八种机器学习模型中，ranger（随机森林）模型具有最高中位AUC，表现出最佳预测性能和稳健性。SHAP分析显示TIPARP、ACSL1、PLIN2、FKBP5、ETS2和B3GNT5对模型输出贡献最大。

HHM关键基因验证与药物敏感性分析：27个关键基因在HC、Early和Late三组中的表达水平被评估。TIPARP被选为深入研究基因，其在高缺氧巨噬细胞中表达显著升高，且在巨噬细胞中表达水平高于其他细胞类型。轨迹分析显示TIPARP高表达细胞主要分布于轨迹中间及部分分支区域。代谢通路分析显示HHM组在氨基酸代谢、脂质代谢、糖代谢和外源物质代谢通路中富集水平更高。药物敏感性预测显示WH-4025、BMS-777607和TWS-119位于高差异敏感性药物区域，其中BMS-777607具有最广泛的覆盖范围，三种候选药物的bcscore均呈现HHM > IHM > LHM趋势。

TIPARP与缺氧的实验验证：研究人员使用LPS刺激RAW264.7巨噬细胞建立细胞模型，Western blot和qRT-PCR分析显示LPS处理显著增加HIF1A和TIPARP的蛋白及mRNA表达。在缺氧（2% O₂）条件下，TIPARP蛋白随时间逐渐增加，伴随HIF1A和VEGF的平行诱导。炎症刺激以浓度依赖方式增加TIPARP、HIF1A和VEGF表达。通过shRNA敲低TIPARP后，在缺氧条件下HIF1A蛋白积累显著增加，BNIP3表达上调，HIF1A和VEGF的缺氧诱导表达进一步增强。这表明TIPARP虽与HIF-1α同步激活，但其后续效应可抑制HIF1A通路的过度激活，支持TIPARP在细胞缺氧反应中发挥负反馈调控功能。

讨论与结论

讨论部分指出，脓毒症日益被理解为一种动态的免疫-代谢紊乱，炎症、组织缺氧和细胞代谢重塑并行演变。缺氧并非灌注受损的被动后果，而是通过HIF依赖性转录程序主动重塑先天免疫功能。本研究的核心发现是巨噬细胞在脓毒症免疫群体中表现出最显著的缺氧相关转录活性，这与巨噬细胞作为高度可塑性哨兵细胞的公认作用一致。传统M1/M2极化框架不足以描述脓毒症中巨噬细胞的极化状态，巨噬细胞占据由炎症配体、代谢底物和组织微环境线索塑造的连续激活状态。将巨噬细胞分为低、中、高缺氧状态，将缺氧定位为巨噬细胞异质性的组织轴，是概念性进展。机制上，高缺氧巨噬细胞表型与炎症放大紧密耦合，HIF1A对髓系炎症反应至关重要，支持缺氧和炎症信号形成前馈回路。高缺氧巨噬细胞可能作为脓毒症免疫微环境中的通讯枢纽，通过MIF–CD74/CD44、CCL家族趋化因子和IL1相关信号协调白细胞募集和炎症传播。基因调控层面，LDHA和ENO1的鉴定强化了糖酵解在缺氧适应巨噬细胞中的核心作用。TIPARP/PARP7作为单ADP-核糖基转移酶，参与缺氧反应、先天免疫调控和转录反馈控制，是引人注目的候选调控因子。药物敏感性预测提示高缺氧巨噬细胞可能存在可靶向的脆弱性。

研究结论部分总结：本研究支持缺氧作为脓毒症中巨噬细胞异质性的核心决定因素这一模型。高缺氧巨噬细胞代表一种可塑、代谢重编程且通讯活跃的状态，可能放大炎症损伤，同时暴露状态特异性治疗脆弱性。通过整合单细胞转录组学、轨迹推断、细胞间通讯、共表达网络分析、机器学习和实验验证，该工作为理解缺氧驱动的巨噬细胞重编程提供了框架，并突出TIPARP作为未来机制和治疗研究的候选调控因子。