《Frontiers in Immunology》：Comorbidity risk characteristics of rheumatoid arthritis in the context of depression-associated lipid metabolism
编辑推荐：
摘要背景：类风湿关节炎（RA）与重性抑郁障碍（MDD）表现出显著共病性，并具有炎症及脂质代谢异常等共享病理机制。然而，连接这两种疾病的特异性分子特征与功能通路仍不明确。
方法：本研究整合了5个MDD队列和5个RA队列的公共转录组数据。在MDD队列中筛选脂质代
摘要背景：类风湿关节炎（RA）与重性抑郁障碍（MDD）表现出显著共病性，并具有炎症及脂质代谢异常等共享病理机制。然而，连接这两种疾病的特异性分子特征与功能通路仍不明确。
方法：本研究整合了5个MDD队列和5个RA队列的公共转录组数据。在MDD队列中筛选脂质代谢相关模块基因，并与RA中上调的差异表达基因进行交集分析，以构建候选基因集。采用最小绝对收缩和选择算子回归（LASSO）及机器学习方法〔极端梯度提升（XGBoost）、随机森林（Random Forest）〕进行特征筛选与模型构建。功能探索包括蛋白质-蛋白质相互作用（PPI）网络分析、免疫浸润分析以及多队列验证。为实验性检验所提出的脂质-免疫联系及潜在神经元-滑膜轴，研究人员建立了体外模型。海马HT-22神经元接受糖皮质激素与饱和脂肪酸联合刺激，以模拟MDD相关应激-脂毒性状态。随后，将这些神经元的条件培养基（CM）作用于人类成纤维样滑膜细胞（FLS）。通过实时定量聚合酶链反应（qPCR）、酶联免疫吸附测定（ELISA）及蛋白质印迹（Western blot, WB）评估FLS的炎症反应与功能变化。
结果：研究鉴定出一个由9个基因组成的脂质-免疫特征签名：ANXA3、IL18、CD59、TNFSF13B、BMX、WASF1、SLC8A1、COMMD8和NXT2。该特征签名在独立RA队列中表现出优异的判别性能〔曲线下面积（AUC）为0.86–1.00〕，但在MDD队列中的性能较弱（AUC为0.55–0.65）。TNFSF13B、IL18和CD59被一致识别为核心枢纽基因。该特征签名的表达与RA中的免疫细胞浸润显著相关。关键的是，体外验证表明，来源于脂质应激神经元的CM可直接激活FLS。其证据包括：（i） CD59 mRNA及蛋白水平均显著上调（经qPCR与WB证实）；（ii） 关键炎症介质IL-6和MMP3分泌增加（经ELISA证实）；以及（iii） FLS迁移能力增强。这些数据支持一种来自受压神经元至滑膜细胞的假定旁分泌通信轴，并提示脂质介导的炎症重编程可能参与其中。
结论：本研究鉴定出一个与MDD-RA共病风险相关的保守脂质-免疫基因特征签名，并提供了涵盖生物信息学预测、临床队列相关性及直接实验依据（qPCR、ELISA、WB）的多层级验证。该特征签名可能代表一部分MDD个体中被激活、且与RA病理重叠的免疫炎症状态。基于枢纽基因，研究提出了“脂质代谢—膜微环境—免疫炎症”轴假说。经实验证实的神经元-滑膜旁分泌轴以特定分子变化（CD59、IL-6、MMP3）为特征，为跨疾病病理过程提供了一个新颖且可检验的机制通路，也为理解精神-免疫共病提供了新的分子学见解与研究方向。
该文发表于《Frontiers in Immunology》，围绕类风湿关节炎（rheumatoid arthritis，RA）与重性抑郁障碍（major depressive disorder，MDD）的共病问题展开，试图从脂质代谢异常这一共同生物学背景中，解析两者之间潜在的分子连接机制。既往研究已显示，RA患者中抑郁发生率较高，而MDD人群未来发生RA的风险亦升高，提示两种疾病之间并非单向伴随关系，而可能具有双向关联基础。现有解释多集中于慢性炎症，但关于脂质代谢失衡如何参与精神疾病与自身免疫疾病之间的跨系统耦联，仍缺乏系统性的分子证据。尤其是，以往脂质组学或转录组学研究多局限于单病种分析，缺少跨疾病、跨队列的整合框架，也缺乏将生物信息学发现进一步延伸到功能实验验证的研究路径。因此，开展本研究的意义在于：在MDD相关脂质代谢异常中识别可能指向RA易感性的分子特征，并验证其是否能够支持一种“神经—滑膜”跨组织炎症通讯模式。
研究人员整合了5个MDD队列与5个RA队列的公共转录组数据，首先在MDD中构建脂质代谢相关的共表达模块，再与RA中上调差异表达基因取交集，筛选潜在共病候选基因；随后利用最小绝对收缩和选择算子回归（LASSO）、极端梯度提升（XGBoost）和随机森林（Random Forest）建立特征模型，并结合蛋白质-蛋白质相互作用（PPI）网络及免疫浸润分析进行功能解析；最后通过海马HT-22神经元应激/脂毒性模型及成纤维样滑膜细胞（FLS）条件培养基刺激实验，对脂质-免疫联系及神经元-滑膜旁分泌轴进行体外验证。研究最终鉴定出9基因脂质-免疫特征签名，并指出TNFSF13B、IL18与CD59为最稳定的核心枢纽。该签名在RA中具有很高判别能力，而在MDD中判别能力较弱，提示其更可能反映MDD部分个体中与RA病理重叠的免疫炎症状态，而非MDD普适性诊断标志。体外实验进一步证明，脂质应激神经元来源条件培养基可激活FLS，诱导CD59上调，并增强IL-6、MMP3分泌和细胞迁移。上述发现为“脂质代谢—膜微环境—免疫炎症”轴提供了支持，也提出了精神-免疫共病的一个可检验分子框架。
研究采用的主要关键技术方法包括：基于GEO数据库的10个独立转录组队列整合分析；使用单样本基因集富集分析（ssGSEA）计算脂质代谢评分，并通过加权基因共表达网络分析（WGCNA）识别MDD中与脂质代谢高度相关的模块；在合并并经ComBat校正批次效应的RA队列中进行差异表达分析；利用LASSO、XGBoost和Random Forest进行特征筛选、建模与解释；采用STRING、Cytoscape及拓扑学指标进行PPI网络与枢纽基因识别；在RA队列中应用CIBERSORT评估免疫细胞浸润；实验部分以小鼠海马HT-22细胞和人FLS构建体外旁分泌模型，并通过qPCR、ELISA、WB、划痕实验与CCK-8实验开展验证。
在研究结果部分，论文按多个层次展开了由计算筛选到实验验证的完整证据链。
3.1 Identification of lipid metabolism–related modules in MDD cohorts
研究人员先在MDD队列中通过ssGSEA获得每个样本的脂质代谢评分，再应用WGCNA识别与该评分最相关的共表达模块。在GSE98793中，turquoise模块与脂质代谢评分相关性最强；在GSE251778中，black模块相关性最高。由此说明，MDD内部存在与脂质代谢活性稳定相关的转录调控模块，为后续跨疾病整合提供了基础基因集。
3.2 PPI network and functional enrichment analysis of lipid metabolism–related genes
两个MDD队列关键模块取交集后得到1021个共享脂质模块基因。PPI网络显示这些基因间存在紧密连接，提示其具有较强的功能协同性。基因本体（GO）富集结果表明，这些基因不仅涉及脂质结合、离子结合、激酶活性和跨膜转运等分子功能，也显著富集于细胞因子应答、白细胞活化、炎症反应、细胞黏附、骨架重塑及物质运输等生物学过程。这一步结果说明，MDD中的脂质代谢相关程序并非孤立的代谢事件，而与免疫炎症调控深度耦联。
3.3 Differential expression analysis in RA combined cohorts and identification of candidate comorbidity genes
研究将GSE93272和GSE17755两个RA队列标准化后合并，并通过ComBat校正批次效应。差异表达分析共得到1434个差异表达基因，其中686个为上调基因。将RA上调基因与MDD脂质模块基因取交集后，得到59个候选共病相关基因。该结果表明，MDD脂质代谢相关转录程序中，确有一部分基因与RA炎症激活方向一致，构成潜在的跨疾病重叠分子层。
3.4 Model interpretation and network topology of the signature genes
在59个候选基因基础上，研究人员利用LASSO进一步压缩特征，最终得到9基因签名：ANXA3、IL18、CD59、TNFSF13B、BMX、WASF1、SLC8A1、COMMD8、NXT2。随后构建XGBoost与Random Forest模型，并通过SHAP分析解释各基因对RA分类的贡献。结果显示，IL18、TNFSF13B和CD59在不同模型中均具有较高重要性。与此同时，对59个候选基因建立PPI网络并开展拓扑中心性分析，也再次将TNFSF13B、IL18和CD59识别为核心节点，说明这3个基因不仅在统计建模中稳定，而且在网络结构中占据关键位置。
3.5 Selection and expression validation of the nine-gene signature
在表达层面，9个特征基因在RA患者中均显著上调；而在MDD队列中，这些基因在高、低脂质代谢评分组之间表现出差异性表达。该结果说明，这一签名兼具两方面特征：一方面在RA中呈现明确激活状态，另一方面其在MDD中与脂质代谢活性背景相关，支持其作为连接MDD脂质失衡与RA免疫病理的候选特征集合。
3.6 Predictive performance evaluation of the nine-gene signature
模型评估显示，9基因签名在RA发现队列和外部验证队列中均表现出优异判别能力，AUC达到0.86–1.00；但在MDD发现队列及外部验证队列中的AUC仅为0.55–0.65。论文据此指出，该签名并不适合作为一般意义上的MDD诊断工具，而更可能代表一类“RA样”的免疫炎症状态，此状态可能存在于部分MDD个体中，并与其未来RA共病风险增加相关。这一结果是全文最重要的解释支点之一，即“强于RA、弱于MDD”的表现模式本身，恰恰支持其共病风险标记而非单病种诊断标记的定位。
3.7 Independence of the lipid–immune signature from phenotypic BMI
研究进一步检查了体质指数（BMI）等肥胖表型对分子特征的潜在混杂影响。在可获得BMI信息的MDD队列中，BMI与脂质代谢评分、9基因签名评分及核心枢纽基因均未呈现稳定显著相关。作者据此认为，所识别的脂质-免疫轴并非主要由表型肥胖驱动。这一分析提高了研究结论对“脂质代谢异常并不等同于单纯肥胖相关炎症”的区分度。
3.8 Immune infiltration characteristics of the nine feature genes
在MDD队列中，9个特征基因与甘油磷脂代谢、鞘脂代谢、NOD样受体信号通路及自然杀伤细胞介导细胞毒性等多种脂质代谢和免疫相关通路显著相关。在RA队列中，这些基因与单核细胞、中性粒细胞、B细胞及多种T细胞亚群浸润水平存在显著相关性，并且两个RA队列中的相关模式较为一致。这说明该签名并非单纯转录标志，而与RA免疫微环境重塑具有实质关联。
3.9 Conditioned medium from lipid-stressed neurons activates synovial fibroblasts and upregulates hub genes
实验验证部分是本文的重要延伸。研究人员以糖皮质激素皮质酮（CORT）和棕榈酸（PA）联合刺激HT-22海马神经元，构建MDD相关应激-脂毒性模型。qPCR结果显示，FKBP5、SGK1、CHOP与GRP78等应激和内质网应激标志物显著升高，证明模型建立成功。随后，将不同处理组HT-22细胞的条件培养基作用于FLS。结果显示，FLS中CD59、IL18、TNFSF13B的mRNA表达受到调节，其中CD59上调最稳定，尤其在CORT+PA联合处理来源的条件培养基组最明显；而IL18和TNFSF13B在FLS中呈处理依赖性下调，提示核心基因在不同细胞类型中的调控方向可能存在差异。功能上，条件培养基可显著提高FLS分泌IL-6和MMP3，并增强细胞迁移能力；CCK-8结果显示细胞活性亦发生变化，其中联合处理组最为明显。进一步的WB证实，CD59蛋白在FLS中显著上调，且以联合处理组最强。该组结果直接支持“受压神经元—滑膜细胞”旁分泌通信轴的存在，并提示脂质相关神经元应激可诱导滑膜细胞炎症及侵袭表型激活。
在讨论部分，论文首先总结了其主要发现：借助跨疾病、多队列整合策略，研究人员从MDD脂质代谢相关模块出发，识别出RA高度相关的9基因脂质-免疫签名，并通过网络分析、免疫浸润分析和体外实验获得多层验证。作者认为，该签名最重要的意义不在于一般性疾病分类，而在于揭示一种可能存在于部分MDD个体中的免疫激活状态，这种状态与RA核心病理发生重叠。围绕TNFSF13B、IL18和CD59，论文提出“脂质代谢—膜微环境—免疫炎症”轴假说：脂质应激可能改变胆固醇和鞘脂组成，进而影响脂筏（lipid rafts，膜微区）中的受体聚集、补体调节和炎症信号放大，从而连接精神应激与滑膜免疫病理。讨论还指出，CD59在FLS中的稳定上调提示膜相关免疫调节可能是这一跨疾病轴线中更直接的执行环节，而IL18和TNFSF13B在FLS中的差异性变化则反映出明显的细胞类型特异性。此外，论文也强调目前证据仍以关联和体外验证为主，具体旁分泌介质、体内因果链条及反向“滑膜—神经元”通路仍待进一步研究。
论文结论部分可译为：
综上所述，本研究通过整合计算发现与实验验证，提出了一个以脂质代谢为中心、可供检验的MDD-RA共病机制假说。研究不仅筛选出新的分子枢纽，还为直接的神经-滑膜通信轴提供了初步功能学支持。经实验验证的神经元-滑膜旁分泌轴为跨疾病发病机制提供了新的理论框架。尽管这些分子通路的具体介质仍有待明确，但本研究为未来旨在阻断MDD-RA共病中脂质-免疫联系的转化研究提供了依据。
