单细胞转录组学整合靶向cis-Mendelian随机化识别骨关节炎候选治疗靶点

《Genes》：Single-Cell Transcriptomics Integrated with Targeted cis-Mendelian Randomization Identifies Candidate Targets for Osteoarthritis

背景：改善疾病的骨关节炎药物（DMOADs）仍是紧迫且未被满足的临床需求。本研究整合单细胞转录组学（scRNA-seq）、靶向cis-蛋白质数量性状位点（cis-pQTL）Mendelian随机化（cis-MR）与贝叶斯共定位，以识别并分层骨关节炎（OA）的遗传学支持治疗靶点。方法：分析人OA关节scRNA-seq图谱（GSE152805）以描绘细胞亚群及标志基因。采用来自deCODE Genetics与UK Biobank PPP的先导cis-pQTL工具，对FinnGen OA GWAS汇总统计量进行两样本cis-MR。使用Steiger滤波确认因果方向性，并通过贝叶斯共定位（coloc）评估共享因果变异。将候选靶点映射至药理学数据库以获得假设生成性重定位机会。结果：scRNA-seq描绘出 distinct 细胞群体，显示炎症性滑膜成纤维细胞中THY1与IL6显著上调、肥大软骨细胞中MMP13上调、衰老软骨细胞中CDKN2A上调。cis-MR证明，遗传预测的较高循环THY1水平（MR Beta = 0.626，OR = 1.87，p = 1.76 × 10^?7）、CDKN2A（MR Beta = 0.432，OR = 1.54，p = 5.38 × 10^?5）与MMP13（MR Beta = 0.405，OR = 1.50，p = 4.33 × 10^?5）与膝OA风险增加显著相关；IL6呈现提示性关联。通过共定位的证据分层确认THY1位点具有强共享因果变异支持（PP.H4 = 91.2%），而CDKN2A与MMP13呈现不确定或仅为提示性共定位。计算机映射识别出假设生成性干预措施，包括针对THY1的潜在Notch抑制剂及靶向CDKN2A相关衰老细胞存活通路的senolytics。结论：将单细胞分辨率与cis-MR整合为OA发病机制提供了分层人类遗传学证据。THY1成为由MR与强共定位共同佐证的稳健支持候选靶点，而CDKN2A与MMP13目前缺乏充分共定位支持并需进一步精细定位。这些发现为需未来实验确认的假设生成性治疗开发提供了遗传学知情框架。

研究背景：骨关节炎（OA）是全球最常见的慢性关节疾病，影响超过5亿人，是老龄化人群中慢性疼痛、身体受限和社会经济负担的主要驱动因素。病理生理学上，OA不再被视为单纯的软骨“磨损”，而是以关节软骨破坏、软骨下骨硬化、骨赘形成和复发性滑膜炎症（滑膜炎）为特征的全关节疾病。尽管临床影响巨大，当前药物治疗仍严重局限于姑息性镇痛药（如对乙酰氨基酚、非甾体抗炎药[NSAIDs]和度洛西汀）及关节内治疗（包括皮质类固醇和透明质酸）。这些干预可能提供短暂症状缓解，但长期疗效往往不足，且均不能阻止或逆转潜在的结构性关节退化，关节置换手术仍是终末期OA不可避免的最后手段。这一关键缺口凸显了靶向因果病理机制的改善疾病骨关节炎药物（DMOADs）的紧迫临床需求。此类DMOADs的发现历来受制于深刻的关节细胞异质性以及临床前动物来源靶点在人体临床试验中频繁的转化失败。单细胞RNA测序（scRNA-seq）的进展通过以单细胞分辨率解析软骨和滑膜组织的转录景观，革新了对关节生物学的理解，揭示了包括THY1⁺（CD90⁺）炎症性滑膜成纤维细胞和CDKN2A⁺（p16^INK4a）衰老软骨细胞在内的病理改变细胞亚群。然而，观察性单细胞转录组学本质上捕捉的是关联而非因果关系，疾病簇中上调的基因可能代表原发致病驱动因素，也可能仅是组织破坏引起的反应性下游旁观现象。为弥合转录组相关性与临床因果关系之间的关键缺口，人类遗传学方法——特别是Mendelian随机化（MR）——提供了强大的分析框架。MR利用单核苷酸多态性（SNPs）作为工具变量（IVs）推断可改变暴露（如循环蛋白水平）对疾病结局的潜在因果效应。由于等位基因在受孕时随机分配，MR从根本上抵抗传统环境混杂和反向因果。此外，大规模蛋白质数量性状位点（pQTL）数据集——尤其是来自deCODE Genetics和UK Biobank Pharma Proteomics Project（UKB-PPP）的数据——使得候选蛋白的靶向cis-Mendelian随机化（cis-MR）成为可能。严格聚焦于cis变异（位于蛋白编码基因近端）显著降低水平多效性可能性，并直接模拟靶向药理学干预。然而，cis-MR信号仍可能被连锁不平衡（LD）混杂，即蛋白丰度和疾病风险的不同因果变异位于同一单倍型上。整合贝叶斯共定位分析（coloc）通过评估两性状共享相同因果遗传变异的后验概率来严格应对这一挑战，从而过滤假阳性LD假象。将验证的因果靶点与药物-基因相互作用数据库（如DrugBank和DGIdb）耦合，创建了高效的计算机药物重定位管道，以识别现有、安全性已验证的药物用于新临床适应证。本研究建立了整合多组学框架，结合人OA关节软骨和滑膜scRNA-seq数据集与大规模cis-pQTL蛋白质组学数据、来自FinnGen联盟的广泛OA GWAS汇总统计量以及计算机药物映射。关键的是，研究人员实施了严格的证据分层，基于MR与共定位联合支持对靶点分类，旨在优先考虑稳健的、假设生成性OA治疗机会。

研究人员开展的研究与结论：研究人员开展了整合单细胞转录组学、靶向cis-MR与贝叶斯共定位的多组学分析，以识别和分层OA的遗传学支持治疗靶点。研究得出，THY1成为由MR与强共定位（PP.H4 = 91.2%）共同佐证的稳健支持候选靶点；CDKN2A与MMP13虽具显著MR信号但缺乏确定性共定位支持，需谨慎解读并进一步精细定位。将分层遗传学见解与精确药理学机制（如senolytics对CDKN2A相关细胞的凋亡清除）对齐，提供了生物学基础的、假设生成性路线图，以指导未来需严格实验确认的改善疾病骨关节炎策略。其重要意义在于为OA发病机制提供了分层人类遗传学证据框架，并为假设生成性治疗开发提供了遗传学知情的基础。

关键技术方法：研究获取GEO公开人OA关节软骨与滑膜scRNA-seq数据集（GSE152805，来自3例独立临床膝OA供体，初始含9个样本文库的10,640个滑膜细胞和26,192个软骨细胞），经Scanpy与Seurat质控、Harmony整合、Leiden聚类与UMAP可视化，并采用供体水平 pseudobulk 与DESeq2进行差异表达分析。遗传工具取自deCODE Genetics（N = 35,559）与UKB-PPP（N = 54,219）血浆蛋白质组GWAS，按±500 kb cis窗口、p < 5 × 10^?8、LD clumping（r² < 0.001）、F > 10严格筛选先导cis-pQTL。结局取自FinnGen Release 11膝OA端点（M13_KNEE_ARTHROSIS，56,213例病例与311,463例对照）。采用TwoSampleMR包Wald比值估计量行两样本cis-MR，Steiger滤波确认方向性，Bonferroni校正（p < 0.0071）；以coloc包行贝叶斯共定位（PP.H4 > 0.80为强证据）；并将候选靶点查询DGIdb v4.0与DrugBank v5.1.9进行药理相互作用映射。样本队列来源已在文中注明。

研究结果：

3.1 单细胞转录组异质性：通过对人OA关节组织（n = 3供体，12,480个过滤细胞）的无监督聚类，解析出主要亚群：稳态软骨细胞（32.4%）、炎症性滑膜成纤维细胞（19.6%）、肥大软骨细胞（17.5%）、成纤维性软骨细胞（15.3%）与衰老软骨细胞（15.2%）。供体水平pseudobulk差异表达分析显示病理亚簇中特定标志物相对大体完整内部对照组织显著富集：THY1（平均log₂FC = 1.85）与IL6（log₂FC = 1.75）在炎症性滑膜成纤维细胞中稳健上调；MMP13（log₂FC = 2.10）在肥大软骨细胞中占主导；衰老软骨细胞表现衰老标志物CDKN2A（log₂FC = 1.95）与CDKN1A（log₂FC = 1.45）显著表达。在pseudobulk分析中未达严格倍数变化阈值的DLK1被排除于下游遗传优先排序之外。

3.2 cis-pQTL工具与靶向cis-MR：严格±500-kb cis窗口过滤识别出高置信度遗传工具，所有工具均满足F > 10（范围47.5–216.3）。针对FinnGen Release 11膝OA的cis-MR显示，遗传预测较高THY1表达显著升高膝OA风险（MR Beta = 0.626，SE = 0.120，OR = 1.87，95% CI：1.48–2.37，p = 1.76 × 10^?7，Bonferroni p = 1.23 × 10^?6）；CDKN2A升高（MR Beta = 0.432，SE = 0.107，OR = 1.54，95% CI：1.25–1.90，p = 5.38 × 10^?5）与MMP13升高（MR Beta = 0.405，SE = 0.099，OR = 1.50，95% CI：1.24–1.82，p = 4.33 × 10^?5）在多检验校正后与膝OA风险增加因果相关；循环IL6呈名义显著关联（MR Beta = 0.293，SE = 0.115，OR = 1.34，95% CI：1.07–1.68，p = 0.0109），严格归类为提示性遗传证据。Steiger方向性滤波确认遗传工具在蛋白暴露中解释的方差显著大于OA结局（所有Steiger p < 0.001），表明主要因果方向为蛋白至疾病。

3.3 贝叶斯共定位证据分层：在THY1位点（chr11q23），分析共3421个变异，贝叶斯共定位对单一共享因果变异产生强支持（PP.H4 = 91.2%，PP.H3 = 4.8%），低PP.H3概率排除LD混杂，稳健验证THY1因果信号；不同p12先验（10^?6至10^?4）的敏感性分析确认H4后验概率稳定性。相反，MMP13（PP.H4 = 34.2%，PP.H3 = 45.1%）与IL6呈现仅为提示性或缺乏共定位。关键的是，CDKN2A表现复杂区域连锁结构（不确定PP.H4 < 15%），表明其MR信号虽统计学显著，但需进一步精细定位验证，在当前模型下不能明确宣称为共享因果位点。

3.4 假设生成性治疗机会：候选治疗药物基于生物学机制映射，并严格按衍生的遗传证据层级分层以准确反映转化准备度。THY1（Tier 1证据）：由MR与强共定位共同支持，假设性Notch信号调节剂与研究中抗CD90抗体的临床前评估代表主要假设生成方向。CDKN2A（Tier 2证据）：映射至Senolytic方案（Dasatinib + Quercetin；Navitoclax）；这些药物并非直接酶促抑制p16^INK4a蛋白，而是通过抑制抗凋亡存活通路（如BCL-2）靶向CDKN2A相关衰老细胞群体的潜在senolytic药物；dasatinib的FDA批准状态严格适用于肿瘤适应证而非OA。IL6（Tier 3证据）：仅由提示性MR证据支持，FDA批准的IL-6受体拮抗剂（Tocilizumab、Sarilumab）代表高度炎症性OA内型的潜在重定位候选，但被归入较低遗传证据层级。

讨论部分总结：本研究建立了证据分层多组学管道以评估OA细胞类型特异性遗传学支持候选靶点，通过综合单细胞转录组学、靶向cis-MR与共定位，提供了分层人类遗传学证据等级，为假设生成性治疗开发提供稳健理据。最重要的发现是遗传预测循环THY1水平与膝OA风险升高之间强Tier-1潜在因果关联，并由稳健贝叶斯共定位独特强化（PP.H4 = 91.2%）；单细胞分辨率将THY1特异定位至炎症性滑膜成纤维细胞，共定位数据表明调节THY1表达的遗传变异可能影响OA发病机制，将滑膜成纤维细胞活化定位为遗传学支持的候选致病节点。分析亦识别CDKN2A（p16^INK4a）与MMP13的显著MR因果链接，但缺乏确定性共定位支持（Tier 2）；软骨细胞肥大、细胞衰老与低度炎症的复杂作用日益被视为OA发病机制的协同驱动因素，衰老软骨细胞永久性增殖停滞并分泌破坏性衰老相关分泌表型（SASP），直接触发邻近基质降解；药理映射将该通路链接至senolytic方案，这些药物并非p16直接酶阻断剂，而是选择性破坏衰老细胞促存活通路（如BCL-2）触发凋亡；尽管实验模型显示软骨细胞中整体p16^INK4a失活或缺失不直接改变OA发展，但靶向凋亡清除衰老细胞已显示关节保护效应；然而MR结果支持CDKN2A通路的生物学相关性，不确定共定位意味着特定因果变异仍待精细定位，关系需谨慎实验功能验证。代谢功能障碍、低度炎症与OA严重度的临床相互作用值得审慎考虑，近期临床研究（包括STEP 9试验相关分析和代谢队列）强调全身炎症介质深刻影响关节结构损伤与患者报告临床结局，为靶向炎症通路提供概念基础；然而IL6的MR结果仅产生提示性证据（Tier 3），且近期评估IL-6受体抑制剂tocilizumab在手OA的随机试验未能证明优于安慰剂的疼痛缓解，凸显OA炎症靶向的复杂性，提示广泛IL-6抑制对泛化OA人群可能不足，而靶向高度特异性炎症内型（或靶向上游滑膜成纤维细胞调节因子如THY1）的精准医学方法可能是更可行治疗途径。关键方法学考量涉及组织水平基因转录、全身循环蛋白水平（所用MR暴露）与局部关节内蛋白丰度之间的生物学区别；虽然滑膜血供丰富促进血液-关节屏障双向交换，全身pQTL信号捕捉的全身遗传易感性不能无缝等同于局部关节组织转录，因此发现支持全身循环蛋白作为遗传易感性代理，但推断精确关节内蛋白丰度需未来局部滑液pQTL数据集。本研究具有若干关键优势，特别是通过cis限制严格控制水平多效性及暴露与结局队列零重叠；此外scRNA-seq分析中通过DESeq2的pseudobulk聚合有效控制供体水平伪重复。然而必须承认关键局限：第一，scRNA-seq数据集仅含3例生物学OA供体，严重限制广泛临床代表性，这些发现主要作为生物学定位工具而非独立临床证据；第二，CDKN2A与MMP13缺乏确定性共定位阻止其明确分类为遗传学验证靶点；此外依赖单一遗传变异（k = 1）限制利用常规MR多效性算法的能力，cis限制与Steiger检验均不能完全消除复杂反向因果或LD混杂的可能性；最后MR估计反映终身累积遗传暴露，不能直接预测成熟成人中急性短期治疗干预的疗效。

研究结论翻译：将单细胞转录组学与靶向cis-Mendelian随机化及共定位整合，为骨关节炎发病机制提供了分层人类遗传学证据框架。THY1成为由MR与强共定位共同佐证的稳健支持候选靶点。相反，CDKN2A与MMP13的显著MR信号缺乏确定性共定位，需谨慎解读并进一步精细定位验证。将这些分层遗传学见解与精确药理学机制——如senolytics对CDKN2A相关细胞的凋亡清除——对齐，提供了生物学基础的、假设生成性路线图，以指导未来需严格实验确认的改善疾病骨关节炎策略。