IF50.5｜196万人大规模骨关节炎基因组学揭示700个效应基因

迄今最大规模骨关节炎GWAS，发现513个新关联，整合多组学锁定700个效应基因。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

骨关节炎是全球增长最快的致残疾病之一，预计2050年患者将达10亿。目前尚无改变病情的药物，亟需深入理解其发病机制。遗传因素在骨关节炎中起重要作用，既往GWAS已发现约150个风险变异，但样本量有限，且多数变异位于非编码区，难以直接指向因果基因。骨关节炎具有明显的关节异质性，不同关节部位的遗传结构可能不同，但既往研究多聚焦于总体表型，缺乏关节特异性分析。此外，骨关节炎的遗传风险如何通过特定细胞类型和生物学通路发挥作用，仍不清楚。

科学问题

卡在哪里： 此前最大规模骨关节炎GWAS仅纳入约74万样本，发现的关联位点有限，且对非欧洲人群覆盖不足。更重要的是，从风险变异到因果基因的转化存在巨大鸿沟：绝大多数变异位于非编码区，缺乏功能注释，难以确定其靶基因和生物学效应。

本文要回答： 通过大规模多祖先GWAS meta分析，系统鉴定骨关节炎风险变异，并整合多组学数据，定位效应基因，揭示关键生物学通路，为药物重定位提供遗传学证据。

技术路线

作者主要用了：GWAS meta分析、fGWAS、单细胞多组学、精细定位、基因负荷分析、pQTL共定位。

作者首先对11个骨关节炎表型进行跨87个队列、近196万人的GWAS meta分析，鉴定出962个独立关联信号。随后利用人类骨骼发育图谱的单细胞ATAC和RNA-seq数据进行fGWAS富集分析，发现信号富集于软骨细胞、成骨细胞等特定细胞状态。进一步通过精细定位、3D染色质互作、转录因子结合分析等24条证据线，对286个风险位点内的基因进行打分，鉴定出700个效应基因。最后对效应基因进行罕见变异负荷分析、通路富集分析和药物靶点注释，揭示关键生物学过程和潜在治疗机会。

图 1　Fig. 1。图1展示了962个独立关联变异的OR值随风险等位基因频率的分布，每个点的大小代表表型方差解释量（VEP）。不同颜色对应不同骨关节炎表型。可见大多数变异为常见变异（MAF≥5%），效应量较小（OR<1.2），VEP普遍较低；少数低频变异（MAF 1-5%）效应量略大。该图直观呈现了骨关节炎的遗传结构：多基因、微效、常见变异为主。

核心亮点

亮点 1｜962个独立关联，513个新发现

在P≤1.3×10-8阈值下，共鉴定出962个独立关联信号，其中513个为首次报道。这些信号分布在286个基因组位点，176个为新位点。95%的关联变异MAF≥5%，效应量较小（OR 1.016-1.186），另有49个低频变异（MAF 1-5%）信号。此外，还发现3个女性特异和1个男性特异的关联。遗传风险评分（GRS）预测能力有限，最佳AUC为髋骨关节炎的58.6%，加入BMI后提升至66%。

图 2　Fig. 2。图2展示了骨关节炎GWAS信号在骨骼发育图谱30种细胞状态中的fGWAS富集结果。点的大小代表效应量（log OR），颜色代表显著性（FDR<0.1）。在软骨细胞谱系中，成熟、肥大、关节软骨细胞和DLK1+软骨细胞显著富集；成骨谱系中骨细胞、成骨细胞和软骨膜细胞也富集，尤其在髋和手指骨关节炎。该图表明胚胎骨骼发育中的特定细胞状态与成年后骨关节炎风险相关。

亮点 2｜信号富集于骨骼发育细胞状态

fGWAS分析显示，骨关节炎GWAS信号在软骨细胞谱系中显著富集，包括成熟、肥大、关节软骨细胞和DLK1+软骨细胞，以及高细胞周期活性的软骨细胞。在成骨谱系中，成熟骨细胞、成骨细胞和软骨膜细胞也显示富集，尤其在髋和手指骨关节炎中。此外，肌腱细胞富集于髋骨关节炎和全髋关节置换信号，提示肌腱发育在疾病中的作用。

亮点 3｜700个效应基因，ALDH1A2得分最高

整合24条正交证据线，对286个位点内的8785个基因进行打分，鉴定出700个效应基因（得分≥3），覆盖88%的位点。最高分基因ALDH1A2获得11条证据支持。小鼠和人类肌肉骨骼及疼痛表型、软骨细胞Hi-C和差异甲基化是信息贡献最高的证据类型。70%的位点含多个效应基因，提示多基因共同作用。

亮点 4｜罕见变异负荷效应量更大

对700个效应基因进行LOF和LOF+MIS负荷分析，发现ADAMTSL3、VIT、COL27A1、IL11、PMVK等基因的罕见变异与骨关节炎显著相关，其中ADAMTSL3和VIT的髋骨关节炎关联达到全基因组显著。效应量普遍大于常见变异，如ADAMTSL3的LOF负荷OR为1.44，而常见变异OR仅1.05。ADAMTS6、SPRY2和COLGALT2的LOF+MIS变异具有保护作用。

亮点 5｜八条生物学通路汇聚效应基因

通路富集分析揭示八条相互关联的生物学过程：视黄酸信号、TGFβ信号、BMP信号、WNT信号、FGF信号、ECM组织、昼夜节律和胶质细胞相关通路。其中ECM和WNT信号贡献的遗传力最高。例如，ALDH1A2和CYP26B1调控视黄酸合成与降解；TGFB1和SMAD6在软骨下骨显示等位基因表达不平衡；SFRP4表达与髋骨关节炎风险等位基因相关。昼夜节律通路首次在基因组水平与骨关节炎关联。

亮点 6｜10%效应基因是已批准药物靶点

在700个效应基因中，69个基因的蛋白产物是473种已批准药物的靶点，其中5个基因与疼痛表型相关。例如，CYP26B1抑制剂taralazole正在进行拇指基底骨关节炎临床试验；FGF18配体sprifermin在膝骨关节炎中显示改善软骨厚度；SOST抑制剂romosozumab用于骨质疏松；NR3C1激动剂糖皮质激素用于骨关节炎疼痛。这些发现为药物重定位提供了直接遗传学依据。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： GWAS meta分析、fGWAS、单细胞多组学、精细定位、基因负荷分析、pQTL共定位
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事骨关节炎遗传学、功能基因组学或药物重定位研究的科研人员。

带走一句

大规模GWAS结合多组学功能注释，可将风险变异转化为效应基因和通路，为复杂疾病机制解析和药物重定位提供可复用的分析框架。