类风湿关节炎中阶段特异性的ROMO1：通过整合GWAS、转录组、单细胞和空间图谱对MIF通路及HLA-DR/IL2RA轴的预测性免疫洞察

《Autoimmunity》：Stage-specific ROMO1 in rheumatoid arthritis: predictive immune insights into the MIF pathway and HLA-DR/IL2RA axis via integrated GWAS, transcriptomic, single-cell, and spatial profiling

　　

新证据将活性氧调节因子1（reactive oxygen species modulator 1，ROMO1）——一个关键的线粒体ROS（mitochondrial reactive oxygen species，mROS）调节因子——与类风湿关节炎（rheumatoid arthritis，RA）发病机制联系起来。然而，由于其具体功能存在矛盾证据，其确切机制仍不清楚。研究人员采用四级整合框架，结合多组学数据与文献支持的机制推断。在遗传水平，进行孟德尔随机化（Mendelian randomization，MR）以探索ROMO1、IL2RA、HLA-DR、MIF与RA风险之间的潜在因果关系，随后进行差异表达分析与基于机器学习的特征选择以识别关键mROS基因，并评估ROMO1在RA进展中的时序表达动态。在细胞和组织水平，整合单细胞RNA测序与空间转录组学，绘制ROMO1相关免疫细胞与通路的细胞类型特异性表达及滑膜定位。最后，将多组学发现置于文献支持的机制推断中。（1）MR结果符合ROMO1对RA的潜在保护作用（OR=0.52），及其对危险因素IL2RA（OR=0.46）与HLA-DR（OR=0.40）的潜在调控；而IL2RA（OR=1.42）、HLA-DR（OR=1.88）与MIF（OR=1.17）与RA风险正相关；ROMO1被确定为具有阶段特异性动态的关键诊断预测因子：早期下调、晚期/缓解期上调。（2）单细胞RNA测序显示ROMO1在CD14+ HLA-DR+ CD74+单核细胞与CD4+ IL2RA+ T细胞中细胞特异性表达；细胞通讯分析提示这些细胞可能参与MIF通路调控；空间转录组学发现ROMO1相关细胞定位于滑膜病理区域，MIF通路变化与RA进展相关。（3）文献支持的机制推断提示，ROMO1可能通过调节mROS水平促进抗炎M2巨噬细胞极化，从而在理论上减轻全身炎症与RA多器官衰退。该整合多组学研究支持ROMO1作为RA中阶段依赖性生物标志物候选与潜在免疫调节因子。

研究背景方面，类风湿关节炎（rheumatoid arthritis，RA）是最常见自身免疫性炎症关节炎，全球疾病负担持续上升，约1760万人受累且患病率较1990年增加14.1%，预计2050年前病例数继续增长。RA主要累及关节，导致慢性疼痛、僵硬与肿胀，个体功能严重受损。尽管治疗选择多样，因复杂自身免疫机制与高疾病负担仍不可治愈，需长期管理。早期干预可延缓关节破坏与残疾，但缺乏可靠早期生物标志物制约及时诊断与精准治疗。近期线粒体功能障碍尤其是线粒体活性氧（mitochondrial reactive oxygen species，mROS）失调在自身免疫病中的作用受关注；活性氧调节因子1（reactive oxygen species modulator 1，ROMO1）作为线粒体内膜跨膜蛋白与mROS关键调节因子，在肿瘤中研究多，在RA中机制未明。既往跨病遗传分析提示ROMO1为自身免疫病共享基因，可调控mROS信号并影响HLA-DR、IL2RA与MIF通路；但近期研究认为ROMO1上调/激活，而研究人员前期RA数据却显示下调，表达模式与生物学意义不一致尚未系统解决，其在RA中经MIF–HLA-DR/IL2RA轴的致病作用也不清，且尚无从DNA到RNA再到细胞、组织与机制推断多层级的综合证据。因此研究人员开展本研究，论文发表于《Autoimmunity》。

主要关键技术方法方面，研究人员使用GEO微阵列数据、MSigDB的mROS基因集做差异表达与机器学习（SVM、LASSO、随机 forest）特征筛选；用GSE159117与GSE244515单细胞数据注释免疫细胞；用GSE195598空间转录组学分析TNF-Tg小鼠滑膜；用胶原诱导关节炎（CIA）小鼠时间序列微阵列做阶段动态；以TwoSampleMR做ROMO1、IL2RA、HLA-DR、MIF与RA的MR；用CellChat做细胞通讯、Monocle3做轨迹、WGCNA做人鼠共表达；并以PubMed/Web of Science/EMBASE文献做机制推断验证。

研究结果部分，Genetic level：通过MR发现ROMO1上调具RA保护效应（IVW OR=0.52），并可能抑制HLA-DR（OR=0.40）与IL2RA（OR=0.46）；IL2RA、HLA-DR、MIF均正向关联RA。差异分析与机器学习从mROS相关基因中锁定ATF2、ROMO1、DHFR、FXN、CHCHD2、MAPK9、SOD2、SLC25A14，随机森林中ROMO1预测价值最高；CIA时序显示ROMO1早期降、晚期/缓解升，H2-Ab1与之拮抗。Cellular level：单细胞聚类得7类免疫细胞，ROMO1在RA中于CD8⁺ T、CD4⁺ T、B、NK、CD14⁺ HLA-DR⁺ CD74⁺单核、CD4⁺ IL2RA⁺ T中显著下调；细胞通讯中MIF–CD74–CD44–CXCR4配体受体对贡献最高，单核细胞为主要信号源。Tissue level：TNF-Tg小鼠滑膜空间分析见MIF相关基因各期共定位，CD74⁺⁺与CD74⁺ H2-Ab1⁺群体随阶段变化；拟时轨迹从WT到Adv2渐复杂并形成反馈环；Mif、Cd44、Cd74、H2-Ab1早期高、晚期降。Cross-species analysis：人鼠中MIF通路分与HLA-DR同源基因强正相关，ROMO1与RA表型负相关，HLA-DR/MIF组分正协表达，IL2RA–MIF–HLA-DR网络跨物种保守。Mechanistic inference level：文献与已有研究提示ROMO1过表达促BMDM向M2极化，转基因小鼠中减轻多器官衰退与慢性炎症、降外周单核/中性粒与循环IL-6，支持ROMO1经mROS–M2极化发挥免疫调节。

讨论部分总结：研究人员指出，相较既往相关分析，本文以多组学加跨物种建立遗传→细胞→组织→机制推断四级证据链，提出ROMO1–ROS–MIF信号级联，填补线粒体功能障碍如何具体调节自身免疫反应的机制空白；HLA-DR4等提升RA易感、IL2RA关联关节破坏已被知，本研究新增ROMO1→IL2RA与ROMO1→HLA-DR因果链假设；阶段特异性表达解释RA异质性与分期精准治疗依据；药物–基因互作与甲状腺自身免疫富集提示ROMO1靶点或超RA。局限包括MR用放宽p值与100 kb窗口、CIA每时点n=3、缺原代Western/qPCR/流式与转基因动物验证，机制仍为计算与相关性推断。

结论部分翻译：本研究中，通过综合多组学调查，研究人员系统刻画了ROMO1在RA中的动态表达模式，并探索其参与MIF信号通路的推定作用；进一步识别其与HLA-DR和IL2RA的潜在协同关联。这些发现为涉及线粒体ROS调节因子的自身免疫病可能机制提供了计算性、假设生成性洞察，并提示ROMO1可作为RA中值得进一步研究的候选靶点。此外，结果也指向ROMO1在其他自身免疫病中的潜在临床相关性，或可启示未来自身免疫病管理的探索策略。