β2-糖蛋白I(APOH)中的编码单核苷酸多态性与自身抗体水平升高但静脉血栓栓塞风险降低相关

A Coding Single Nucleotide Polymorphism in β2-Glycoprotein I (APOH) is Associated with Increased Autoantibody Levels but Reduced Venous Thromboembolism Risk

作者信息Christophe Lalaurie, Lili Liu, Atlas Khan, Chen Wang, Stephen S Rich, R Graham Barr, Elana J Bernstein, Krzysztof Kiryluk, Thomas McDonnell, Yiming Luo
PMID42706482
发布时间2026-09-07
DOI10.1002/art.70333

实验完整度

高

多祖先GWAS(5,969人)与外部VTE GWAS(81,190病例)两独立证据层,含MR、共定位、精细定位、功能基因组及分子动力学模拟

主要模型

MESA多民族队列(5,969名参与者) VTE GWAS汇总数据(81,190例病例和11,419,671例对照) 分子动力学模拟的野生型和W335S β2GPI蛋白模型

重点核对

MESA队列中总抗β2GPI抗体水平经ELISA定量(TheraTest Labs Inc试剂盒),筛选阈值>10.5单位 GWAS分析使用REGENIE v4.1,校正年龄、性别、研究地点、吸烟状况、教育水平和前10个祖先主成分 MR分析使用rs1801690作为工具变量(F-statistic=42),Wald比值法估计VTE风险 贝叶斯共定位分析使用coloc v5.2.3,PP4=0.97 分子动力学模拟从圆形和线性(PDB 1C1Z)模型起始,分析PCA聚类、RMSD、RG和SASA

摘要

目的:抗β2-糖蛋白I(抗β2GPI)抗体是抗磷脂综合征的核心。APOH位点与抗β2GPI相关,但其因果变异和血栓风险意义仍不清楚。方法:我们对来自动脉粥样硬化多种族研究(MESA)的5,969名参与者的定量总抗β2GPI抗体水平进行了多祖先全基因组关联研究(GWAS)。使用两样本孟德尔随机化(MR)和贝叶斯共定位评估遗传决定的抗β2GPI水平与静脉血栓栓塞(VTE)风险之间的关系。通过精细定位和整合功能基因组分析优先考虑候选因果变异。使用分子动力学模拟研究结构效应。结果:我们在APOH位点识别出一个全基因组显著关联(领先SNP rs1801690-G,β = 0.21,p = 1.08 × 10−10)。矛盾的是,APOH位点与较高抗β2GPI水平相关的遗传变异与VTE风险降低相关(β = −0.25,p = 3.95 × 10−6;PP4 = 0.97),可能反映水平多效性而非真正的保护性抗体效应。精细定位和功能基因组学将错义变异rs1801690(W335S)优先为最可能的因果变异。分子动力学模拟支持双重效应机制,即W335S损害Domain V中的磷脂结合,同时增加Domains I-II中的表位暴露。结论:与较高抗β2GPI水平相关的APOH遗传信号矛盾地降低了VTE风险,可能由W335S通过表位可及性和磷脂结合破坏介导。这些发现提供了人类遗传学证据,表明DV的磷脂结合对于抗β2GPI介导的血栓形成至关重要。APOH基因型的临床效用值得进一步研究。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
APOH位点的因果变异及其与抗β2GPI抗体水平和静脉血栓栓塞(VTE)风险的关系是什么?
核心机制
错义变异rs1801690(W335S)通过损害DV的磷脂结合减少血栓风险,同时通过增加DI-II的表位暴露促进自身抗体产生。
主要证据
多祖先GWAS(5,969人)发现APOH位点显著关联;MR和共定位(PP4=0.97)显示遗传决定的抗β2GPI水平与VTE风险负相关;精细定位和功能基因组分析优先rs1801690(W335S);分子动力学模拟显示W335S降低磷脂结合环的SASA并增加抗体反应性基序的暴露。
研究意义
提供了人类遗传学证据支持DV的磷脂结合对于抗β2GPI介导的血栓形成至关重要;APOH基因型的临床效用需要在不同临床背景下的抗β2GPI阳性个体中进一步研究。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

多祖先GWAS发现APOH位点关联

识别与总抗β2GPI抗体水平相关的遗传位点

在MESA队列(5,969人)中对总抗β2GPI抗体水平进行多祖先GWAS meta分析,使用REGENIE v4.1和METAL。

2

MR和共定位评估抗β2GPI水平与VTE风险关系

评估遗传决定的抗β2GPI抗体水平是否因果影响VTE风险

使用rs1801690作为工具变量,Wald比值法进行单变量MR;使用coloc进行贝叶斯共定位。

3

精细定位和功能基因组分析优先因果变异

识别APOH位点内的候选因果变异

使用SuSiE和MESuSiE进行精细定位;评估非编码变异的调控潜力(功能表观基因组学、QTL、深度学习);使用9种预测工具评估错义变异。

4

分子动力学模拟研究W335S结构效应

研究W335S对β2GPI结构和功能的影响

对野生型和W335S β2GPI进行分子动力学模拟,从圆形和线性(PDB 1C1Z)模型起始,分析PCA聚类、RMSD、RG和SASA。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
抗β2GPI抗体ELISA试剂盒TheraTest Labs Inc--
Affymetrix Human SNP array 6.0Affymetrix--
REGENIE v4.1----
METAL----
SuSiE----
MESuSiE----
coloc v5.2.3----
TwoSampleMR v0.6.6----
纳米尺度分子动力学2----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
GWAS分析
样本量:5,969人;抗体测量方法:ELISA(TheraTest Labs Inc);协变量:年龄、性别、研究地点、吸烟状况、教育水平、前10个祖先主成分
阅读提示:Methods部分'Genome-wide association analysis'小节和'Phenotype and genotype measurement'小节
MR分析
工具变量:rs1801690;VTE GWAS样本量:81,190病例和11,419,671对照;MR方法:Wald比值法
阅读提示:Methods部分'Mendelian randomization and colocalization'小节
共定位分析
共定位方法:coloc;后验概率阈值:PP4=0.97
阅读提示:Methods部分'Mendelian randomization and colocalization'小节及Figure 3
精细定位
精细定位方法:SuSiE和MESuSiE;可信集大小:18个变异
阅读提示:Methods部分'Statistical fine-mapping'小节和Table 1
分子动力学模拟
起始结构:圆形模型和PDB 1C1Z;模拟参数:文中未明确说明(在Supplementary Methods中)
阅读提示:Methods部分'Molecular dynamics'小节和Supplementary Methods