鞘脂代谢重编程在类风湿关节炎中医证候辨证中的整合转录组学-代谢组学研究

《Frontiers in Immunology》：Integrative transcriptomics–metabolomics reveals sphingolipid metabolic reprogramming in traditional Chinese medicine syndrome differentiation of rheumatoid arthritis

　　

摘要专业翻译 目的： 探讨鞘脂代谢在类风湿关节炎（Rheumatoid Arthritis，RA）湿热证与肝肾亏虚证鉴别中的作用，识别代谢生物标志物，并整合转录组数据以阐明其分子机制。 方法： 纳入40例RA患者（肝肾亏虚证，n=20；湿热证，n=20），采用超高效液相色谱-串联质谱法（ultra-performance liquid chromatography-tandem mass spectrometry，UPLC-MS/MS）进行血清非靶向代谢组学分析。构建最小绝对收缩和选择算子（least absolute shrinkage and selection operator，LASSO）逻辑回归诊断模型，并通过受试者工作特征（receiver operating characteristic，ROC）曲线和决策曲线分析（decision curve analysis，DCA）评估其性能。从基因表达综合数据库（Gene Expression Omnibus，GEO）下载三个独立转录组数据集（GSE55235：10例RA滑膜、10例骨关节炎滑膜和10例正常滑膜；GSE93272：RA全血与健康对照全血；GSE17755：RA外周血单个核细胞与健康对照外周血单个核细胞），分析鞘脂代谢相关基因的差异表达及组织特异性调控。采用标记基因集平均表达评分法估计免疫细胞浸润，并利用Spearman相关分析构建代谢物-临床指标网络。 结果： 湿热证组鞘氨醇激酶1（sphingosine kinase 1，SPHK1）活性（以1-磷酸鞘氨醇/C17-二氢鞘氨醇比值表示）显著高于肝肾亏虚证组（P = 3.45×10-6）。C17-二氢鞘氨醇与1-磷酸鞘氨醇（sphingosine-1-phosphate，S1P）呈显著负相关（r = -0.555）。单独以C17-二氢鞘氨醇鉴别肝肾亏虚证与湿热证，曲线下面积（area under the curve，AUC）达0.945，准确率为95%。1-磷酸鞘氨醇裂解酶1（sphingosine-1-phosphate lyase 1，SGPL1）在RA滑膜中较正常滑膜显著上调（log2 fold change = 1.124，P = 2.17×10-4），在外周血中也显著上调但效应量较小（log2 fold change = 0.134，P = 8.44×10-4），呈现组织特异性幅度差异调控。SPHK1与中性粒细胞浸润呈负相关趋势。C8-神经酰胺-1-磷酸与C反应蛋白（C-reactive protein，CRP）呈负相关（r = -0.47），与平均红细胞血红蛋白（mean corpuscular hemoglobin，MCH）呈正相关（r = 0.51）。 结论： C17-二氢鞘氨醇是区分肝肾亏虚证与湿热证的核心代谢生物标志物，单次检测即可实现精准证候分类。SGPL1的组织特异性幅度差异调控及SPHK1-中性粒细胞相关趋势，为鞘脂代谢参与RA中医证候辨证提供了新见解。

论文解读

一、研究背景与科学问题

类风湿关节炎（Rheumatoid Arthritis，RA）是一种以慢性滑膜炎症和进行性关节破坏为特征的系统性自身免疫病，全球患病率约为0.5%–1%。尽管改善病情抗风湿药和生物制剂已显著改善RA患者预后，仍有30%–40%的患者对现有治疗应答不佳，提示该病存在显著的分子异质性，亟需更精准的个体化治疗策略。中医药治疗RA历史悠久，辨证论治是个体化治疗的核心。湿热证和肝肾亏虚证是RA临床最常见的两种中医证型：湿热证多见于疾病活动期，表现为关节红、肿、热、痛及炎症指标升高；肝肾亏虚证多见于疾病后期，表现为关节隐痛、腰膝酸软及贫血倾向。既往非靶向代谢组学研究已发现鞘脂代谢是这两种证型间显著扰动的通路。鞘脂代谢是维持细胞膜结构和信号转导的重要网络，其核心信号分子1-磷酸鞘氨醇（sphingosine-1-phosphate，S1P）通过五种S1P受体亚型参与免疫细胞迁移和炎性细胞因子释放。S1P由鞘氨醇激酶1/2（SPHK1/2）合成，并被S1P裂解酶（SGPL1）不可逆降解。已有研究证实S1P信号通路在RA滑膜炎症中发挥关键作用，SPHK1可促进成纤维样滑膜细胞释放炎性细胞因子。然而，氧化应激与鞘脂代谢之间的串扰在RA证候辨证中的作用仍不清楚。本研究旨在整合血清代谢组学与多数据集转录组学，系统评估鞘脂代谢物对RA证候分类的诊断价值，并表征鞘脂代谢基因的组织特异性调控，为RA中医证候的客观化诊断提供分子基础。

二、研究设计与主要技术方法

研究人员开展了一项整合多组学的研究。研究分为三部分：第一部分基于既往研究的代谢组学数据，对40例RA患者（肝肾亏虚证n=20，湿热证n=20）进行血清非靶向代谢组学分析，采用超高效液相色谱-四极杆飞行时间质谱（UPLC-QTOF/MS）检测，通过偏最小二乘判别分析和t检验筛选差异代谢物，并基于KEGG数据库进行代谢通路富集分析。构建LASSO逻辑回归诊断模型，采用ROC曲线和决策曲线分析评估诊断性能。第二部分从GEO数据库下载三个独立转录组数据集（GSE55235：RA滑膜、骨关节炎滑膜和正常滑膜各10例；GSE93272：RA全血与健康对照全血；GSE17755：RA外周血单个核细胞与健康对照外周血单个核细胞），基于KEGG通路定义29个核心鞘脂代谢基因，利用GEO2R进行差异表达分析。第三部分为整合分析，包括代谢通路映射、组织特异性调控分析、免疫细胞浸润评估（采用标记基因集平均表达评分法估计六种免疫细胞相对丰度）及代谢物-临床指标Spearman相关网络构建。

三、研究结果

3.1 研究人群基线特征

40例RA患者按中医证型分层后，肝肾亏虚证（GS）与湿热证（HD）两组在年龄、性别、BMI、病程及常规炎症指标（CRP、ESR）方面均衡可比，用药情况（甲氨蝶呤、糖皮质激素、来氟米特、羟氯喹）亦无显著差异，排除了用药作为主要混杂因素的可能。但两组间胆碱酯酶、尿酸和单核细胞百分比存在显著差异。经严格的FDR校正后，鞘脂代谢物与40项临床指标间无显著相关，表明鞘脂代谢特征独立于常规临床参数，支持C17-二氢鞘氨醇作为证候特异性生物标志物而非系统性炎症或器官功能障碍的替代指标。

3.2 鞘脂代谢物对中医证候分类的诊断性能

研究发现C17-二氢鞘氨醇在GS组显著升高（P = 2.32×10^-4），而S1P(d18:1)在HD组显著升高（P = 4.19×10^-4），SPHK1相关代谢比值（S1P/C17-二氢鞘氨醇峰面积比）在HD组显著高于GS组（中位数4.66 vs. 0.11，P = 3.45×10^-6），提示GS以C17-二氢鞘氨醇蓄积伴SPHK1活性受抑为特征，HD则与S1P产生增强和SPHK1活性升高相关。Spearman相关分析显示C17-二氢鞘氨醇与S1P呈显著负相关（r = -0.555，P < 0.001），且在逐步校正CRP、ESR、胱抑素C及甲氨蝶呤和糖皮质激素使用后，该负相关仍稳定存在（完全校正模型偏相关r = -0.505，P = 0.0032），表明这种反向模式独立于系统性炎症、肾功能状态和临床用药。ROC曲线分析显示，C17-二氢鞘氨醇单独诊断AUC达0.945，优于S1P(d18:1)（AUC = 0.809），与C17/S1P比值（AUC = 0.894）相当（DeLong检验Z = 0.892，P = 0.372）。联合C17-二氢鞘氨醇和C17/S1P比值的逻辑回归模型正确分类20/20例GS患者（100%特异性）和18/20例HD患者（90%敏感性），总体准确率95%。多因素逻辑回归校正CRP、ESR和甲氨蝶呤使用后，C17-二氢鞘氨醇仍与中医证型显著相关（OR = 0.9999，95% CI：0.9998–0.99997，P = 0.018），确立其作为独立代谢生物标志物的价值。

3.3 决策曲线分析

决策曲线分析显示，C17-二氢鞘氨醇在0–1.0全阈值范围内（0.85–1.00）保持正净获益，显著优于“全部按HD治疗”和“均不治疗”两种极端策略。随着阈值概率增加，“全部治疗”策略的净获益迅速降至零以下，而C17-二氢鞘氨醇始终维持稳定的正净获益，表明无论临床医生要求高或低的诊断确定性，检测C17-二氢鞘氨醇辅助中医证候辨证均具有临床净获益和广泛实用价值。

3.4 RA滑膜中鞘脂代谢基因的差异表达

基于GSE55235数据集（RA滑膜n=10、OA滑膜n=10、正常滑膜n=10），在鞘脂代谢基因panel中，18个基因在三种两两比较中差异表达。RA与OA比较中，ASAH1（log₂FC = -1.635，adj.P.Val = 0.041）和CERS4（log₂FC = -1.582，adj.P.Val = 0.017）显著下调，UGCG（log₂FC = 1.000，adj.P.Val = 0.035）和SGPL1（log₂FC = 0.636，adj.P.Val = 0.017）显著上调，SPHK1无明显变化（log₂FC = 0.053，P = 0.818）。RA与正常比较中，SGPL1变化最显著（log₂FC = 1.124，adj.P.Val = 2.63×10^-6），伴随SPTLC2、CERK和SMPD2同步上调（均log₂FC > 0.67，P < 0.015），SPHK1和SPHK2呈非显著下降趋势。值得注意的是，OA与正常比较中SGPL1呈现截然相反的模式，显著下调（log₂FC = -0.488，adj.P.Val = 0.048），CERK亦下调（log₂FC = -0.657，adj.P.Val = 0.024）。整合热图分析显示SGPL1具有RA特异的双向调控特征，CERS4和ASAH1在RA与OA比较中选择性下调，SPTLC2上调在RA和OA中均存在，提示滑膜病理中鞘脂从头合成普遍增强。SPHK1和SPHK2在所有比较中均无显著差异表达，提示RA滑膜中S1P合成主要受翻译后水平或底物可用性调控而非转录调控。

3.5 多数据集验证鞘脂代谢基因表达

为验证上述发现的跨组织一致性，纳入GSE93272（全血）和GSE17755（外周血单个核细胞，PBMCs）两个数据集。SGPL1在外周血中的表达模式与滑膜部分一致：在GSE93272中显著上调（log₂FC = +0.134，P = 8.44×10^-4），但在GSE17755中无差异表达，且效应量小于滑膜，提示转录调控幅度具有组织背景依赖性。SPHK1在两个外周血数据集中均无显著差异表达（GSE93272：log₂FC = +0.084，P = 0.151；GSE17755：log₂FC = +0.018，P = 0.767）。SPHK2仅在GSE93272中显著上调（log₂FC = +0.176，P = 6.09×10^-6），在GSE17755中无变化。CERS4在GSE93272中无显著差异表达（log₂FC = +0.091，P = 0.189），在GSE17755中未检出。ASAH1在GSE93272中显著下调（log₂FC = -0.441，P = 1.40×10^-4），方向与其在RA与OA滑膜比较中的下调一致，在GSE17755中无显著变化。UGCG在GSE93272中显著下调（log₂FC = -0.384，P = 8.65×10^-6），在GSE17755中未检出。跨数据集整合热图显示，SGPL1在所有比较中效应量最大，是跨数据集验证中最稳健的差异表达基因。

3.6 鞘脂代谢通路的整合转录组-代谢组分析

将GSE55235三比较转录组数据映射至KEGG鞘脂代谢通路并与既往血清代谢组学发现整合。通路分析显示RA滑膜中鞘脂代谢呈双向调控特征：从头合成端SPTLC2显著上调（log₂FC = 1.31），提示上游合成通路激活；神经酰胺代谢分支中SMPD2（log₂FC = 0.68）、ASAH1（log₂FC = 0.42）和CERK（log₂FC = 0.68）均上调，分别促进鞘磷脂水解、神经酰胺降解和C1P产生；S1P代谢轴中，尽管SPHK1（log₂FC = -0.64）和SPHK2（log₂FC = -0.38）呈下降趋势，但催化S1P降解的SGPL1显著上调（log₂FC = 1.12），SGPP1亦上调（log₂FC = 0.49），提示RA滑膜中S1P降解增强。既往代谢组学数据显示HD组血清S1P显著高于GS组，与滑膜中降解酶上调形成对比，提示关节可能通过SGPL1上调加速S1P清除以限制局部炎症。组织特异性分析进一步揭示SGPL1在滑膜和GSE93272全血中均显著上调但外周血效应量小得多，揭示了鞘脂代谢基因调控在RA不同组织区室间的幅度差异模式。

3.7 免疫细胞浸润分析及与鞘脂基因的相关性

基于GSE55235数据集，采用标记基因集平均表达评分法估计正常（n=10）、OA（n=10）和RA（n=10）滑膜中六种免疫细胞特征评分。所有六种免疫细胞特征评分在RA滑膜中均显著高于正常和OA滑膜（均P < 0.05）。CD8⁺ T细胞、CD4⁺ T细胞、B细胞和NK细胞在RA中较正常滑膜高约3–4倍，巨噬细胞和中性粒细胞亦显著升高。对RA和OA合并滑膜样本（n=20）中六个鞘脂相关基因与免疫细胞特征评分进行Spearman相关分析，经Benjamini-Hochberg FDR校正后：SGPL1与CD8⁺ T细胞（r = 0.737，q = 0.009）、CD4⁺ T细胞（r = 0.720，q = 0.009）、B细胞（r = 0.653，q = 0.017）、巨噬细胞（r = 0.657，q = 0.017）和NK细胞（r = 0.592，q = 0.031）显著正相关。CERS4与B细胞（r = -0.679，q = 0.012）和NK细胞（r = -0.599，q = 0.027）负相关。UGCG与CD8⁺ T细胞（r = 0.630，q = 0.022）和CD4⁺ T细胞（r = 0.555，q = 0.050）正相关。SPHK1与中性粒细胞呈负相关趋势（r = -0.400，q = 0.161），但FDR校正后未达显著。多数涉及SPHK2和ASAH1的相关性较弱（|r| < 0.40，q > 0.1）。

3.8 谷胱甘肽-鞘脂代谢串扰机制

基于谷胱甘肽代谢（GS与HD间变化最显著的通路）和鞘脂代谢（HD与GS间差异通路）的双重扰动，以及转录组数据中SGPL1的差异表达，研究人员提出氧化应激驱动的谷胱甘肽-鞘脂代谢串扰模型：RA关节局部产生的过量活性氧通过两条途径加剧炎症——一方面消耗谷胱甘肽并削弱GPX4介导的铁死亡抑制；另一方面促进神经酰胺生成并驱动鞘脂代谢流向S1P合成。S1P进而激活NF-κB信号驱动促炎细胞因子分泌，形成氧化应激驱动鞘脂重塑、鞘脂重塑又放大炎症的正反馈循环。该模型为HD患者氧化应激严重且炎症水平高的临床特征提供了代谢机制解释，也为靶向谷胱甘肽-鞘脂代谢轴的治療策略提供了理论基础。

四、讨论与结论

讨论部分指出，本研究首次系统证明了鞘脂代谢重编程在RA中医证候辨证中的关键作用。C17-二氢鞘氨醇单独即可高效区分GS与HD（AUC = 0.945，准确率95%），且与CRP、ESR无显著相关，表明该代谢生物标志物反映独立于传统炎症指标的代谢特征。转录水平上，SGPL1是表现最突出的基因——在RA滑膜和外周血中均显著上调但血液中效应量较小。SGPL1编码S1P裂解酶，负责S1P不可逆降解。研究人员推测RA滑膜中SGPL1升高代表对慢性局部S1P信号的反馈性代偿反应，但这种代偿可能不足：上游S1P合成可能超过降解能力，且翻译后机制可能限制SGPL1蛋白功能。SPHK1在滑膜和外周血数据集中均无显著差异表达，SPHK2仅在GSE93272全血中显著上调，提示RA中S1P合成主要受翻译后或底物可用性调控。CERS4和ASAH1在RA与OA比较中选择性下调，可能作为区分RA与OA的转录生物标志物。免疫微环境中SPHK1和CERS4与中性粒细胞呈负相关趋势，SGPL1与CD8⁺和CD4⁺ T细胞呈中度正相关，提示参与适应性免疫细胞浸润。研究局限性包括：数据来自不同队列缺乏样本水平匹配、样本量小且单中心、SPHK1活性为间接计算、免疫浸润为算法估计、串扰模型需实验验证、所有机制推断均为相关性而非因果关系。

结论部分：本研究揭示了鞘脂代谢重编程在RA中医证候辨证中的关键作用。C17-二氢鞘氨醇是区分肝肾亏虚证与湿热证的核心代谢生物标志物，单次检测即可实现精准证候分类。SGPL1的组织特异性幅度差异调控和SPHK1-中性粒细胞相关趋势为理解RA证候辨证的分子机制提供了新视角，谷胱甘肽-鞘脂代谢串扰模型为RA精准诊疗提供了潜在靶点。