《Nature Communications》：Large-scale molecular endotype discovery in synovial fluid reveals osteoarthritis as a single biological continuum
编辑推荐：
膝关节骨关节炎（Osteoarthritis, OA）累及人群广泛且临床表型异质性强，其分子发病机制尚不明确，学界长期争议OA是单一疾病还是由多条不同通路驱动的疾病集合。研究人员建立"滑液无偏蛋白质组学检测OA分子内型（Synovial fluid To de
膝关节骨关节炎（Osteoarthritis, OA）累及人群广泛且临床表型异质性强，其分子发病机制尚不明确，学界长期争议OA是单一疾病还是由多条不同通路驱动的疾病集合。研究人员建立"滑液无偏蛋白质组学检测OA分子内型（Synovial fluid To detect Endotypes by Unbiased Proteomics in OA, STEpUP OA）"联盟，对1361例膝OA患者滑液样本进行无偏蛋白质组学分析以搜寻分子内型（Molecular Endotype），并解析关键病理生物学通路。采用无监督聚类（Unsupervised Clustering）分析后，仅观测到代表生物学连续谱（Biological Continuum）的单一集群，主要由"上皮间质转化（Epithelial Mesenchymal Transition, EMT）"驱动，未检测到截然不同的分子内型。按临床表型（肥胖状态、生物学性别）分层后，"血管生成（Angiogenesis）""补体（Complement）"及"凝血（Coagulation）"通路显著富集；补体与凝血通路与炎症标志物C反应蛋白（C-reactive Protein, CRP）相关；与患者自我报告的膝痛关联性较弱。上述发现支持膝OA为生物学连续谱而非离散内型集合，鉴定出共有及表型富集的可靶向通路，为临床试验设计中受试者分层提供依据。
一、研究背景与目的
膝关节骨关节炎（Knee Osteoarthritis, OA）是常见关节退行性疾病，60岁以上成人患病率可达三分之一，目前尚无获批的结构改善药物（Disease-Modifying Osteoarthritis Drugs, DMOADs），症状治疗药物疗效有限且不良反应明显。OA临床表现高度异质，长期存在"OA是单一疾病还是多种亚型集合"的学术争议，既往提出的临床表型缺乏治疗反应预测及疾病进展验证价值，基于分子特征定义的分子内型（Molecular Endotype）被认为可能更好解释表型差异并指导分层。既往OA分子分型多使用外周血，而滑液（Synovial Fluid, SF）直接接触病变关节组织，较血液更能反映局部病理分子特征。本研究由STEpUP OA（Synovial fluid To detect Endotypes by Unbiased Proteomics in OA）联盟开展，旨在通过大样本膝OA滑液无偏蛋白质组学分析，检验膝OA是否存在可区分的分子内型，解析与结构及症状性疾病相关的蛋白特征与信号通路，相关成果发表于《Nature Communications》。
二、主要关键技术方法
研究人员整合17个队列共1361例确诊膝OA患者及36例非OA对照者的单份膝关节滑液样本，按发现集（N=708）与验证集（N=653）划分。滑液采用SomaScanTMv4.1适配体平台检测6596种人源蛋白质，经标准化及批次效应校正（ComBat校正离心状态）。无监督聚类采用PCA降维结合k-means及f(K)统计量判断最佳簇数，引入细胞内蛋白评分（Intracellular Protein Score, IPS）回归以排除胞内蛋白造成的方差干扰。蛋白与临床特征关联采用线性回归、logistic回归或有序比例优势模型，校正年龄、性别及队列，显著性采用Benjamini-Hochberg法（FDR<5%），复制定义为两数据集同方向padj≤0.05。通路富集分析采用MSigDB Hallmark基因集及fgsea算法，蛋白互作网络基于STRING数据库构建。部分分析按肥胖（BMI≥30 kg/m2）及性别分层，并与已发表OA软骨及滑膜RNA-seq数据进行比对。
三、研究结果
Endotype detection in OA SF（OA滑液中分子内型的检测）
在未进行IPS回归的数据中，f(K)指标提示存在2个聚类，但在IPS回归后聚类消失；UMAP可视化显示两簇界限模糊且可被连续IPS二分化解释，按性别、影像学严重度（KL分级）、血样污染分层后均仅在非IPS回归数据中见模糊双簇。表明未调整IPS时出现的分组实为同一生物学连续谱，非独立分子内型，且此连续分布不依赖疾病分期、性别或肉眼血染。
Synovial fluid protein associations with radiographic OA（滑液蛋白与影像学OA的关联）
以晚期影像学OA（KL 3-4级）vs非晚期（KL 0-2级）做关联分析，发现3815个蛋白显著相关，跨数据集复制688个，多数随病情进展上调。Hallmark通路富集显示上皮间质转化（Epithelial Mesenchymal Transition, EMT）、补体（Complement）、血管生成（Angiogenesis）在全部数据集中与晚期OA显著相关，调整血红蛋白A（血污染替代指标）后仍显著。EMT通路含TIMP1、TIMP3、MMP-2、TGFβ1、VEGFA等已知软骨基质重塑相关分子；与已发表OA软骨及滑膜RNA-seq比对显示EMT在两种组织均富集，补体与凝血仅见于滑膜。
Synovial fluid protein associations with advanced radiographic OA after stratification by BMI or biological sex（按BMI或性别分层后滑液蛋白与晚期影像学OA的关联）
按肥胖分层后，肥胖与非肥胖组分别与晚期OA关联的蛋白效应高度相关（r=0.72），无显著交互项；两组均强富集EMT，仅肥胖组额外显著富集补体与凝血通路。正式交互检验未发现BMI显著修饰任一蛋白-疾病关联。按性别分层后男女蛋白效应亦高度相关（r=0.69），均无显著性别交互；两性均富集EMT，仅男性额外显著富集血管生成与凝血通路（调整血红蛋白A后仍显著）。补体与凝血通路蛋白水平与血清CRP正相关（调整BMI后亦然）。
Synovial fluid protein associations with WOMAC pain in OA（OA滑液蛋白与WOMAC疼痛的关联）
以WOMAC膝痛亚分为结局，联合数据集发现797个显著相关蛋白，但Discovery与Replication集无复制、跨集相关性弱（r=0.36），调整影像学严重度后关联蛋白数减少。仅Noelin-2（NOE2）和Ecto-ADP-ribosyltransferase 3（NAR3）在合并集及调整队列后保持显著且与KL分级无关。神经生长因子（Nerve Growth Factor, NGF）与晚期影像学OA正相关但与WOMAC疼痛无显著关联。提示滑液蛋白与自我报告疼痛关联不稳定，可能受心理、镇痛药及疼痛评分变异性影响。
四、讨论与结论总结
研究人员通过迄今最大规模膝OA滑液无偏复制性横断面蛋白质组学分析发现，膝OA滑液未呈现独立可分的分子内型，而是表现为沿EMT等通路强度渐变的生物学连续谱（Biological Continuum），个体在同一通路上呈连续分布。滑液蛋白关联分析证实EMT通路稳健激活，反映关节组织（软骨及滑膜）活跃重塑；补体及凝血通路富集具表型特异性——肥胖者补体/凝血增强、男性血管生成/凝血增强，且与系统性炎症标志物CRP关联，提示可用于临床试验受试者富集分层。与结构改变相比，滑液蛋白与自我报告疼痛关联弱且不可跨队列复制，符合疼痛受多因素影响的特性。局限性包括SomaScan平台对滑液蛋白检测的新颖性及部分队列穿刺成功率记录不全。研究支持膝OA作为生物学连续谱而非离散内型疾病的观点，鉴定出共有及临床表型富集的可靶向通路（EMT、补体、凝血、血管生成），为OA精准分型及试验设计提供依据。STEpUP OA数据已开放外部研究申请。
原文结论译文：本研究呈现STEpUP OA联盟主要结果——迄今最大规模无偏、可重复、横断面的膝OA滑液蛋白质组学分析。研究人员阐明疾病中生物学通路的平衡及其随结构性与症状性严重度、患者因素（肥胖、性别）的变化规律。数据未提供截然可分分子内型的证据，提示OA为一生物学连续谱，个体沿给定生物学通路呈连续分布，该信息有助于为合适个体选择恰当治疗。滑液反映关节局部过程，与晚期影像学OA关联的SF蛋白通路分析显示EMT稳健激活（提示组织主动重塑），该特征在OA软骨及滑膜RNA-seq中亦见；补体与凝血仅见于滑膜富集，随患者分层（BMI、性别）及CRP水平变化。肥胖者补体/凝血显著、男性血管生成/凝血显著。与结构疾病相比疼痛相关蛋白关联未能跨队列复制，符合患者自报症状的变异性。上述发现支持膝OA为生物学连续谱，确定共有及表型富集的可靶向通路，为临床试验分层设计提供理论依据。
