《International Journal of Genomics》：Arg1 and Tgfb1 Identify a Partially Shared Macrophage Polarization Program in Postinfarction Cardiac Repair: A Single-Cell Transcriptomic Analysis
编辑推荐：
背景：梗死后心脏修复由巨噬细胞协调，这些巨噬细胞经历从促炎M1到修复性M2状态的顺序极化。这些巨噬细胞亚群调节梗死心肌中的炎症消退、血管生成和纤维化重塑。乳腺癌中的肿瘤相关巨噬细胞（TAMs）与心脏M2巨噬细胞共享转录特征，为跨疾病比较共享免疫程序提供了机会。
背景：梗死后心脏修复由巨噬细胞协调，这些巨噬细胞经历从促炎M1到修复性M2状态的顺序极化。这些巨噬细胞亚群调节梗死心肌中的炎症消退、血管生成和纤维化重塑。乳腺癌中的肿瘤相关巨噬细胞（TAMs）与心脏M2巨噬细胞共享转录特征，为跨疾病比较共享免疫程序提供了机会。理解跨病理背景下修复性巨噬细胞共有的分子程序，可能识别出未来治疗研究中的候选调控节点。
方法：单细胞RNA测序数据来自GEO：GSE136088（小鼠AMI后心脏巨噬细胞）、GSE176078和GSE167036（人类原发性乳腺癌）。分析包括质量控制、UMAP/t-SNE降维、基于图的聚类、扩散伪时间轨迹重建、配体-受体通讯分析、主成分分析和跨数据集相关性分析。所有分析均在Python中使用scanpy进行，并在R（biomaRt，14,892个一对一小鼠-人类直系同源基因）中进行跨物种直系同源映射，后续分析在Python（scanpy）中进行。
结果：分析QC后数据集（GSE136088：约14,900个细胞；GSE176078：100,064个细胞；GSE167036：49,141个细胞）。UMAP和t-SNE在所有三个数据集中解析出多个巨噬细胞和免疫亚群。扩散伪时间分析重建了一个具有M1/M2命运分叉的连续激活轨迹。在心脏巨噬细胞中鉴定出一个7基因共表达模块（Spp1、Mif、Tgfb1、Arg1、Il10、Tnf和Mrc1），并在乳腺癌巨噬细胞群体中显示出部分保守的表达模式（GSE176078：Pearson r = 0.72；GSE167036：检测到7个基因中的6个，ARG1未检测到）。配体-受体通讯分析确定MIF-ACKR3、SPP1-CD44和TGFB1-TGFBR2为心脏和肿瘤微环境中介导巨噬细胞-基质通讯的共享配体-受体对。心脏巨噬细胞与乳腺癌巨噬细胞群体之间7基因模块的跨数据集Pearson相关性在GSE176078中达到r = 0.72，而在GSE167036中为r = ?0.07（无显著差异）。
结论：这项单细胞转录组分析表征了梗死后心脏修复中的巨噬细胞极化程序，并评估了其在乳腺癌中的部分保守性。鉴定的7基因模块与GSE176078表现出中等跨数据集相关性（r = 0.72），但在GSE167036中未得到复制，该数据集中未检测到M2极化巨噬细胞。MIF-ACKR3信号轴代表了一个计算推导的假设，用于未来实验验证。这些发现突出了跨疾病状态下巨噬细胞极化的共享特征和背景依赖性差异，强调了在配对心脏和肿瘤模型中进行功能研究的必要性。
**论文解读：梗死后心脏修复与乳腺癌中巨噬细胞极化程序的单细胞转录组比较分析**
**一、研究背景与问题**
急性心肌梗死（AMI）后，心脏修复过程由巨噬细胞精确调控，其经历从促炎M1状态向修复性M2状态的顺序极化。M1巨噬细胞负责清除坏死碎片并放大危险信号，而M2巨噬细胞则促进炎症消退、血管生成和瘢痕巩固。若这一免疫转换失败（如M1持续存在或M2失调），将导致不良心脏重塑和心力衰竭。然而，在单细胞分辨率下，控制M1到M2转换的分子回路及特定免疫基因程序的功能角色仍不完全清楚。此外，缺乏能够忠实标记功能不同免疫状态并预测其旁分泌信号后果的验证分子标志物，是理解心脏巨噬细胞免疫生物学的主要障碍。乳腺癌中的肿瘤相关巨噬细胞（TAMs）被报道与心脏M2巨噬细胞具有转录相似性，但跨疾病比较的分子证据尚不充分。因此，研究人员旨在通过单细胞转录组分析，鉴定心脏修复中巨噬细胞极化的核心分子程序，并利用乳腺癌TAM数据集进行跨疾病验证，以识别共享的免疫调控节点。
**二、研究内容与结论**
研究人员使用单细胞RNA测序（scRNA-seq）数据，分析了小鼠心肌梗死后心脏巨噬细胞及两种人类乳腺癌TAM数据集。通过降维、聚类、伪时间轨迹重建、配体-受体通讯分析及跨数据集相关性分析，鉴定出一个7基因共表达模块（Spp1、Mif、Tgfb1、Arg1、Il10、Tnf、Mrc1），该模块在心脏M2巨噬细胞中富集，并在GSE176078乳腺癌TAM中部分保守（Pearson r = 0.72），但在GSE167036中未复制（r = ?0.07）。同时，MIF-ACKR3、SPP1-CD44和TGFB1-TGFBR2被确定为共享的配体-受体轴，介导巨噬细胞与基质细胞的通讯。该研究揭示了跨疾病背景下巨噬细胞极化的共享特征与背景依赖性差异，并提出MIF-Ackr3轴作为候选调控节点。论文发表在《International Journal of Genomics》。
**三、主要关键技术与方法**
1. **单细胞RNA测序数据获取与质量控制**：数据来自GEO数据库，包括小鼠心脏scRNA-seq（GSE136088，14天左前降支结扎后）和人类乳腺癌scRNA-seq（GSE176078，26例原发性乳腺癌；GSE167036，8例未经治疗患者）。质量控制阈值：GSE136088中nFeature_RNA 500–5000，线粒体比例<20%；GSE176078/GSE167036中UMI≥400，线粒体比例<25%。
2. **降维、聚类与细胞注释**：使用UMAP和t-SNE进行可视化，基于图的聚类（分辨率0.5），通过典型标志物（如Cd68、Nos2、Arg1等）注释免疫细胞类型。
3. **伪时间轨迹与细胞通讯分析**：扩散伪时间分析重建巨噬细胞激活顺序，根设为最早病理状态；配体-受体相互作用分析计算每个细胞类型中配体与受体平均表达值的乘积作为通讯分数，显著交互定义为p<0.05。
4. **跨物种和跨数据集比较**：使用biomaRt（Ensembl v105）映射14,892个一对一小鼠-人类直系同源基因；Pearson相关分析评估7基因模块在心脏与乳腺癌巨噬细胞群体间的相关性，并采用置换检验（100次迭代）进行统计验证。
**四、研究结果**
**结果1：质量控制、批次效应验证与高变异免疫基因鉴定**
严格质量控制后，各数据集细胞数分别为：GSE136088（8555个MI_day14细胞、6337个Sham细胞）、GSE176078（100,064个细胞）、GSE167036（49,141个细胞）。心脏与乳腺癌数据集中高变异基因（HVG）重叠分析显示，共享HVG富集了巨噬细胞极化调控因子，包括Cd68、Arg1、Nos2、Spp1、Mif和Tgfb1。三种独立参考数据库的直系同源映射敏感性分析得到一致的重叠计数（619–641个基因），证实共享免疫特征并非单一映射方法的人为产物。
**结果2：UMAP和t-SNE可视化解析出不同的巨噬细胞免疫亚群**
在心脏数据集（GSE136088）中，UMAP解析出四种主要细胞群体：巨噬细胞、成纤维细胞、内皮细胞和单核细胞。Leiden亚聚类进一步区分出M1样（Nos2
+、Il1b
+、Tnf
+）和M2样（Arg1
+、Mrc1
+、Il10
+）巨噬细胞亚群。在乳腺癌数据集中，UMAP界定出免疫和基质群体，包括巨噬细胞、T细胞、B细胞、NK细胞、上皮细胞和成纤维细胞。t-SNE可视化提供了互补验证。心脏和乳腺癌微环境中观察到的平行免疫亚群结构（特别是M2/TAM-M2簇定位），支持跨疾病转录组比较的概念基础。
**结果3：伪时间轨迹重建梗死后心脏修复中的顺序免疫状态转换**
伪时间轨迹分析重建了心脏巨噬细胞在AMI后修复中的免疫激活动态。主轨迹跨越9个有序状态：状态1–3（M0/单核细胞前体，高Cd68/Itgam，低效应标志物）；状态4–6（M1促炎效应细胞，峰值Nos2、Il1b、Tnf表达）；状态7–9（M2调节性巨噬细胞，上升Arg1、Mrc1、Il10、Tgfb1表达）。在伪时间约40%处发现关键免疫命运分叉，细胞分向M1或M2分支。分叉点基因表达分析鉴定出Mif、Tgfb1和Spp1（7基因模块成员）在分叉节点富集，提示这些基因参与M1/M2决策点的命运决定。
**结果4：心脏与乳腺癌数据集中标志基因表达分析**
在GSE136088中，10个经典巨噬细胞标志物的点图证实了亚群特异性表达模式：泛巨噬细胞标志物（Cd68、Itgam）广泛表达，M1相关基因（Nos2、Il1b、Tnf）在特定亚群富集，M2相关基因（Arg1、Mrc1、Il10、Tgfb1、Spp1、Mif）在M2样亚群富集。跨疾病比较显示，GSE176078中M2相关基因（ARG1、MRC1、IL10、TGFB1、SPP1、MIF）在巨噬细胞群体中富集，M1相关基因（NOS2、IL1B、TNF）标记不同亚群。GSE167036中观察到部分重叠模式，7个模块基因中6个检测到，ARG1缺失。置换检验显示，心脏与GSE176078数据集间7基因模块相关性趋势高于随机预期（p = 0.069）。
**结果5：热图和基因集富集分析（GSEA）支持7基因免疫共表达模块**
热图分析显示，在GSE136088中，11个经典巨噬细胞极化基因在心脏免疫亚群中表现出清晰分界的亚群特异性表达模式。在乳腺癌TAM数据集中，观察到结构类似的表达架构：GSE176078中巨噬细胞群体在Leiden亚簇中富集NOS2、IL1B、TNF（M1相关）以及ARG1、MRC1、IL10、TGFB1、SPP1、MIF（M2相关）。GSE167036中检测到6个模块基因（ARG1缺失）。GSEA证实巨噬细胞群体中M2极化、细胞因子信号转导和免疫调节基因集富集。
**结果6：主成分分析（PCA）揭示共享免疫转录架构**
PCA分析显示，共享模块基因（Spp1、Arg1、Mif、Tgfb1、Il10、Tnf）在心脏和乳腺癌数据集的PC1中均贡献最大。跨疾病将巨噬细胞转录组投影到共享主成分空间，显示心脏与乳腺癌巨噬细胞群体部分重叠，表明M2程序存在转录趋同，同时反映物种和组织背景差异。
**结果7：配体-受体通讯分析揭示共享免疫通讯枢纽**
在GSE136088中，巨噬细胞与成纤维细胞间检测到Mif-Ackr3、Spp1-CD44和Tgfb1-Tgfbr2显著信号。在GSE176078中，对应的人类配体-受体对（MIF-ACKR3、SPP1-CD44、TGFB1-TGFBR2）显示出类似的巨噬细胞-癌症相关成纤维细胞（CAF）通讯模式。跨疾病比较确认MIF-ACKR3、SPP1-CD44和TGFB1-TGFBR2为心脏和肿瘤微环境中共有的主要通讯界面。GSE167036中因缺乏M2极化巨噬细胞，限制了巨噬细胞-CAF轴的直接比较。
**结果8：置换验证的MIF-Ackr3免疫信号轴及独立复制**
MIF-Ackr3信号轴被鉴定为心脏和乳腺癌微环境中共有的显著调控界面。心脏巨噬细胞与GSE176078 TAM群体间7基因模块的跨疾病共表达相关性呈现中等正相关（Pearson r = 0.72，置换p = 0.069），而GSE167036中相关性不显著（r = ?0.07，p = 0.905），与该数据集中M2极化巨噬细胞和ARG1表达缺失一致。MIF表达在心脏和乳腺癌巨噬细胞群体中均检测到，而ACKR3表达在成纤维细胞和CAF中富集，支持MIF-ACKR3信号轴在配体-受体通讯分析中的推断。
**五、讨论与结论总结**
讨论部分指出，该研究核心发现是鉴定了一个7基因共表达模块，其在心脏M2巨噬细胞中富集，并在乳腺癌TAM中部分保守。每个模块基因（Arg1、Mrc1、Il10、Tgfb1、Spp1、Mif、Tnf）在巨噬细胞生物学中均有明确功能，共同支持抗炎、促重塑表型。跨疾病比较的优势在于利用乳腺癌TAM丰富的功能文献，但不足之处包括小鼠心脏数据与人类乳腺癌数据的跨物种差异，以及GSE167036中复制失败凸显数据集特异性巨噬细胞组成的重要性。MIF-Ackr3轴被提出作为候选调控节点，但所有计算发现均需实验验证。此外，讨论强调了治疗悖论：共享M2程序在心脏修复中有利，在乳腺癌中却可能促进免疫抑制，因此靶向该模块需考虑组织特异性递送或时间窗口策略。
**结论**：这项单细胞转录组分析表征了梗死后心脏修复与乳腺癌之间共享的巨噬细胞极化程序。鉴定出一个7基因共表达模块（Spp1/Mif/Tgfb1/Arg1/Il10/Tnf/Mrc1）作为M2极化巨噬细胞的转录特征，该模块与GSE176078表现出中等跨数据集保守性（r = 0.72），但在GSE167036中未复制（r = ?0.07）。MIF-ACKR3、SPP1-CD44和TGFB1-TGFBR2配体-受体轴被确定为心脏和肿瘤微环境中巨噬细胞-基质通讯网络的共享特征。所有发现均为计算推导，作为未来实验验证的假设。靶向一个在一种疾病背景中有益而在另一种中有害的程序的治疗悖论被明确承认，并提出了组织靶向、时间特异性干预策略作为未来研究方向。
