整合单细胞转录组学(single?cell transcriptomics)与孟德尔随机化(Mendelian Randomization, MR)鉴定RAC1作为系统性硬化病周细胞(pericyte)功能障碍因果性代谢驱动因子

《MEDIATORS OF INFLAMMATION》：Integrating Single-Cell Transcriptomics and Mendelian Randomization to Identify RAC1 as a Causal Metabolic Driver of Pericyte Dysfunction in Systemic Sclerosis

　　

背景：系统性硬化病(Systemic sclerosis, SSc)是一种以血管损伤及进行性纤维化为特征的自身免疫性疾病。尽管微血管损伤是始发事件且周细胞(pericyte)被认为是肌成纤维细胞的主要来源，SSc微环境中周细胞的精确表型异质性及驱动其病理性转化的遗传机制仍不清楚。方法：研究人员通过整合SSc患者单细胞RNA测序(scRNA?seq)数据与全基因组关联研究(Genome?Wide Association Study, GWAS)数据，采用双向孟德尔随机化(Mendelian Randomization, MR)分析系统探讨周细胞功能障碍的细胞与遗传学基础。利用拟时序轨迹(pseudotime trajectory)分析重建发育谱系，开展细胞?细胞互作(cell–cell communication)及代谢通路分析以阐明潜在机制，并使用外部批量RNA测序(bulk RNA?seq)数据集进行多组学验证。结果：单细胞分析揭示周细胞亚群存在显著异质性，特异地鉴定出祖细胞样周细胞(progenitor?like pericytes)明显扩增，且该亚群位于向纤维化表型分化轨迹的根部。双向MR分析确定RAC1是SSc显著的因果危险因素(OR = 2.0756, p = 0.0046)。机制上，RAC1阳性周细胞通过MIF?(CD74+CD44)信号轴增强与巨噬细胞的促炎互作，且这些活化周细胞呈现特征性代谢重编程——核黄素(riboflavin)及硫胺素(thiamine)代谢上调以满足生物能量需求。转录组验证进一步确认RAC1在SSc组织中异常过表达。结论：本研究确立RAC1介导的祖细胞样周细胞活化与SSc发病机制间的机制联系。RAC1作为因果驱动因子促进病理性周细胞转化，并通过代谢重编程与免疫招募调控促炎微环境，为SSc提供了新的治疗靶点。

本文对发表于《MEDIATORS OF INFLAMMATION》的研究论文《整合单细胞转录组学(single?cell transcriptomics)与孟德尔随机化(Mendelian Randomization, MR)鉴定RAC1作为系统性硬化病周细胞(pericyte)功能障碍因果性代谢驱动因子》进行解读总结。

研究背景

系统性硬化病(Systemic sclerosis, SSc)是一种异质性自身免疫性风湿病，以广泛血管异常、免疫系统失调及皮肤与内脏进行性纤维化为特征，其特定死亡率在结缔组织病中最高，主要死因为肺纤维化与肺动脉高压。现有证据表明微血管损伤是SSc发病的最早始动事件，包绕毛细血管内皮细胞的周细胞(pericyte)可通过周细胞?肌成纤维细胞转化(pericyte?to?myofibroblast transition, PMT)成为肌成纤维细胞的重要来源，参与血管稀疏及胶原沉积。然而SSc微环境中周细胞的具体表型异质性及驱动其活化的上游遗传因果机制尚不明晰，既往观察性研究无法排除混杂因素与反向因果，亟需结合高分辨率单细胞图谱与因果推断遗传学分型加以阐明。

研究人员整合公开SSc患者scRNA?seq数据集(GSE138669，12例SSc与10例健康对照)与FinnGen SSc GWAS汇总统计数据、eQTLGen全血表达数量性状位点(expression Quantitative Trait Loci, eQTL)数据，通过Seurat进行周细胞亚群分群注释、Slingshot拟时序分析、CellChat细胞互作分析、双向两样本孟德尔随机化(以IVW法为主，辅以MR?Egger与加权中位数法)判定周细胞标志基因表达对SSc风险的因果方向，并以独立bulk RNA?seq数据集(GSE130955)行CIBERSORTx去卷积与相关分析验证，探讨RAC1在周细胞活化、炎症互作及代谢重编程中的作用。

---

主要技术方法概述

研究人员获取GEO数据库SSc患者scRNA?seq数据(GSE138669，12例SSc/10例对照)与bulk RNA?seq验证数据(GSE130955)，使用Seurat(v4.4.0)质控滤除低质量细胞及双细胞(DoubletFinder)，Harmony算法去除批次效应后UMAP降维Louvain聚类，依经典标记基因(RGS5、PDGFRB等)注释周细胞并独立再分群。以Slingshot包(根设定为祖细胞样周细胞)与CytoTRACE行拟时序及干性评分，GeneSwitches识别沿拟时序开关基因。CellChat推断祖细胞样周细胞与其他细胞分泌信号配体?受体互作。双向两样本MR分析以eQTLGen全血eQTL为工具变量(IV，p<5×10^?8，LD r²<0.001，F>10)检验scRNA?seq筛选候选基因(含RAC1)表达对SSc(FinnGen ID: finn?b?M13_SYSTSLCE)的因果效应，反向MR排除反向因果，Cochran's Q与MR?Egger截距检验异质性与水平多效性。外部bulk RNA?seq做Wilcoxon秩和检验差异表达、CIBERSORTx估算祖细胞样周细胞丰度并与RAC1表达及改良Rodnan皮肤评分(modified Rodnan Skin Score, mRSS)做Pearson相关，CMap检索可逆转RAC1⁺促纤维化特征FDA批准药物。

---

研究结果

3.1. Single?Cell Transcriptional Landscape and Cellular Heterogeneity Identification

经严格质控与Harmony批次校正，UMAP将SSc及对照皮肤细胞分为27个亚簇并注释为11种主要谱系(角质形成细胞KRT1/5、成纤维细胞COL1A1/DCN、T细胞CD3D/E、周细胞RGS5/PDGFRB、巨噬细胞LYZ/CD14等)，明确组织微环境细胞组成异质性。

3.2. Subpopulation Identification and Functional Characterization of Pericytes

从全局景观中提取周细胞再行Harmony整合与UMAP聚类得到14个周细胞亚簇，依标记基因注释为5种功能表型：静息周细胞(resting pericytes)、祖细胞样周细胞(progenitor?like pericytes)、纤维化相关周细胞(fibrosis?associated pericytes)、肌成纤维细胞样周细胞(myofibroblast?like pericytes)及活化周细胞(activated pericytes)，SSc组中纤维化相关与活化亚型比例呈升高趋势。

3.3. Expansion and Functional Trajectory of Progenitor?Like Pericytes in SSc

比例分析显示SSc组祖细胞样周细胞与纤维化相关周细胞显著扩增，静息周细胞减少；CytoTRACE证实祖细胞样周细胞干性评分最高。Slingshot拟时序示祖细胞样周细胞居轨迹根部，沿分化方向趋向纤维化相关及活化周细胞。CellChat分析发现祖细胞样周细胞与血管内皮细胞(VEGFA?VEGFR、ANGPTL?ITGA5/ITGB1)、巨噬细胞/树突状细胞(MIF?(CD74+CD44/CXCR4)、MDK?NCL/SDC1)存在强互作。

3.4. Causal Association Between Progenitor?Like Pericyte Marker Genes and SSc

以周细胞标志基因为暴露的双向两样本MR显示，RAC1高表达增加SSc患病风险(IVW: OR = 2.0756, 95% CI: 1.2530–3.4383, 名义p = 0.0046；加权中位数法OR = 2.0886, p = 0.0072)，敏感性分析无显著异质性与水平多效性(Egger截距p = 0.613，Cochran's Q p = 0.452)，个体SNP效应一致；反向MR示SSc对RAC1表达无显著因果(p = 0.270)，排除反向因果。LRRFIP1(OR = 0.5457, p = 0.0052)、CDC37(OR = 0.2806, p = 0.0028)、PSMB9为保护性因子。

3.5. Intercellular Crosstalk and Metabolic Reprogramming of RAC1⁺Pericytes

按RAC1表达分层，RAC1⁺免疫调节(immunomodulatory, IM)周细胞与巨噬细胞(CCL2?ACKR1、CXCL2?ACKR1、MIF?(CD74+CD44))、内皮细胞(CXCL8?ACKR1)互作增强，处微环境信号中枢。GeneSwitches沿拟时序示FOSL1、NR4A1、NFATC2随RAC1上调。通路富集示RAC1⁺IM周细胞显著上调核黄素(riboflavin)代谢、硫胺素(thiamine)代谢及视黄醇(retinol)代谢转录特征，静息周细胞富集生物素(biotin)代谢，提示RAC1驱动活化需特异性代谢重编程提供生物能量与生物合成支持。

3.6. Transcriptional Validation and Expression Profiling of Key MR?Identified Genes

scRNA?seq点图示RAC1广泛表达于周细胞、成纤维细胞及巨噬细胞，PSMB9主要富集免疫细胞。Feature plot确认RAC1在各周细胞表型尤其活化及纤维化相关簇高表达。SSc组RAC1表达较对照显著升高(Wilcoxon p = 3.7×10^?5)，LRRFIP1、CDC37显著下调。bulk队列CIBERSORTx示RAC1高表达亚组祖细胞样周细胞估算丰度更高(p = 0.038)；bulk RAC1表达与临床改良Rodnan皮肤评分(mRSS)正相关(R = 0.31, p = 0.023)。CMap筛选出可逆转RAC1⁺促纤维化特征FDA批准化合物(如利福平rifampicin、马赛替尼masitinib)。

---

讨论与结论总结

讨论部分指出SSc微血管损伤触发氧化应激与炎症介质释放，致驻留周细胞病理性扩增与PMT，是血管损伤向不可逆纤维化转化的关键桥梁。本研究通过scRNA?seq联合双向MR突破既往观察性局限，首次在SSc中鉴定祖细胞样周细胞扩增及RAC1为该群细胞促纤维化转化的因果遗传驱动因子。RAC1编码RAS超家族小GTP酶，调控肌动蛋白 cytoskeleton及NADPH氧化酶(NOX)产活性氧(Reactive Oxygen Species, ROS)，其过表达促进间充质细胞迁移、增殖及表型转化；RAC1⁺周细胞不仅经PMT贡献肌成纤维细胞池，还通过MIF?(CD74+CD44)轴招募活化巨噬细胞、分泌CCL2/CXCL8塑造促炎微环境，并伴随核黄素与硫胺素代谢上调满足高生物能量需求，构成"代谢?免疫?纤维化"协同致病轴。LRRFIP1与CDC37的保护性作用提示SSc受促纤维化与免疫调节失衡共同调控。外部bulk验证与MR因果推断方向一致，支持RAC1变异通过改变自身表达水平驱动下游周细胞病理表型。研究局限性含代谢重编程结论源于转录组通路富集，缺直接代谢组验证，需在动物模型确认RAC1?MIF?CD74轴靶向干预抗纤维化效果。CMap提示利福平及马赛替尼等已有药物具重定位潜力。

结论(翻译)

研究人员通过交叉比对scRNA?seq谱与双向MR数据，解析了周细胞功能失调关联SSc的因果架构与共同病理通路。研究在周细胞群体中表征出一独特祖细胞样亚群，该亚群在SSc条件下显著增殖并主动向纤维化细胞命运定向分化。机制探索揭示RAC1是驱动此病理性转化的因果危险因素，促进以增强免疫细胞招募及特异性代谢重编程为特征的促炎状态。本研究填补了RAC1依赖的周细胞功能紊乱与SSc发生发展间分子空白，所揭示的RAC1信号轴为针对系统性硬化病开发创新性靶向治疗范式奠定理论基础。