Graves眼病(Graves’ ophthalmopathy,GO)又称甲状腺相关性眼病(thyroid-associated ophthalmopathy,TAO),是一种复杂的器官特异性自身免疫性疾病,是全身性甲状腺相关性自身免疫性疾病的眼部表现,即Graves病(Graves’ disease,GD)最常见的甲状腺外表现,可导致眼球突出、眼部软组织红肿、眼睑退缩、复视、暴露性角膜病变和压迫性视神经病变,严重者可引起视力丧失[1-2]。据2021年欧洲甲状腺学会报道,女性GO的发病率为每年(2.67~3.30)/10万,男性为每年(0.54~0.90)/10万[3]。GO的发病机制尚未完全阐明,普遍认为是是一种由免疫、遗传、环境等多因素参与的疾病[4]。近年来,许多遗传学研究[5]和分子研究[6]表明补体系统在GO进展中起至关重要的作用,包括遗传变异、替代途径的过度激活、炎症、氧化应激等。目前,针对补体核心成分3(component 3,C3)和成分5(component 5,C5)及补体调节因子[补体因子D(complement factor D,FD)和补体因子I(complement factor I,FI)等]的抑制剂已超过14种,这些补体抑制剂在年龄相关性黄斑变性中已经完成或正在进行的临床试验近40项[7-10]。
孟德尔随机化(Mendelian randomization,MR)是一种基于全基因组关联研究(genome-wide association studies,GWAS)数据,利用单核苷酸多态性(single nucleotide polymorphism,SNP)作为变量工具(instrumental variable,IV),用于揭示因果关系的新型流行病学方法,效力等同于随机对照试验,其优势在于规避了先前观察性研究的缺点,可最大限度地减少反向因果关系和混淆偏倚的机会[11]。本研究利用MR方法,探究补体C5/C5a与GO的因果关联,以期寻找GO的早期生物标志物及治疗靶点,为其诊断与治疗提供新的理论依据。
资料与方法
研究设计
用于MR分析的IVs需满足3个主要假设:1)从工具变量Z到暴露因素X再到结局变量Y的路径存在强关联,混杂变量C分别与暴露因素X和结局变量Y相关。2)混杂变量C与工具变量Z不存在路径关系。3)从工具变量Z到结果变量Y除了通过暴露外没有其他路径[12-13]。选择血浆蛋白补体C5/C5a作为暴露因素,GO作为结局因素,使用两样本双向MR方法进行因果关联分析,通过共定位分析检测补体C5与GO的遗传关联(图1)。

数据来源
补体C5/C5a蛋白质数量性状位点(protein quantitative trait loci,pQTL)的GWAS汇总数据分别来源于英国生物样本库药物蛋白质组学项目(UK Biobank Pharma Proteomics Project,UKB-PPP)和deCODE数据库。UKB-PPP机构使用基于抗体的 Olink蛋白质组学测定法测试了1 610万个序列变异与54 306名英国生物样本库(UK Biobank)参与者血浆中2 922种蛋白质水平的遗传关联[14],补体C5测定所纳入的有效样本量为33 995;deCODE机构使用基于Soma scan的蛋白质组学测定法测试了2 700万个序列变异与35 559名冰岛人血浆中4 719种蛋白质水平的遗传关联[15],C5a测定所纳入的有效样本量为35 376。GO的汇总数据来自FinnGenR10数据库,共有691例欧洲血统病例和411 490例欧洲血统对照参与了本次GWAS(表1)。
表型 | 机构 | 病例人数 | 对照人数 | PMID | 调整 |
---|---|---|---|---|---|
C5 | UKB-PPP | 33 995 | ─ | 37 794 186 | 年龄、性别、批次、UKB中心、UKB基因阵列、血液采样和测量之间的时间以及前20个遗传主成分 |
C5a | deCODE | 35 376 | ─ | 34 857 953 | 年龄、性别等前20个遗传主成分 |
GO | FinnGenR10 | 691 | 411 490 | ─ | 年龄、性别等前20个遗传主成分 |
本研究采用的所有GWAS汇总数据都是公开的,因此不需要额外的伦理审批。
工具变量
首先,选择与补体C5/C5a显著性相关的cis-pQTL作为工具变量(P<5×10-8)[16-17]。cis-pQTL指的是距离编码对应蛋白质基因的上下游1 Mb区域内的SNP。使用1000Genomes欧洲人群去除存在连锁不平衡的SNP(clump_kb=10 000,r2=0.001)。其次,使用PhenoScanner数据库去除混杂SNP,使用MR-PRESSO测试以识别和去除具有多效性的SNP。最后,为了避免弱工具变量,所有工具变量的F值>10,计算公式:F=β2/SE2,式中β为SNP对暴露的效应量,SE为β的标准误[18-20]。最终纳入UKB-PPP来源的1个SNP、deCODE来源的2个SNP用于MR分析(表2)。
SNP | 蛋白质 | Chr: Pos | EA | OA | EAF | β | SE | P | F |
---|---|---|---|---|---|---|---|---|---|
rs4837810 | UKB-PPP-C5 | 9:121118921 | C | T | 0.577 | -0.143 | 0.008 | 6.55×10-78 | 319.515 |
rs17220750 | deCODE-C5a | 9:121025721 | A | G | 0.107 | 0.186 | 0.013 | 3.78×10-10 | 204.710 |
rs928406 | deCODE-C5a | 9:120365967 | G | C | 0.394 | 0.051 | 0.008 | 4.11×10-48 | 40.641 |
MR分析
对单个SNP使用Wald比率法(Wald ratio)作为主要分析方法。对2个及以上的SNP使用逆方差加权法(inverse variance weighted,IVW)作为主要分析方法[21-23],如果检测到异质性,则选择IVW的随机效应作为结果;否则,IVW的固定效应是首选。对3个及以上的SNP,采用MR-Egger法、加权中位数法(weighted median estimator,WME)作为IVW的补充方法,当MR-Egger、WME的β值与IVW的β值方向一致时,提示IVW结果稳定可靠。此外,本研究采用了最新的cML-MA法,cML-MA是一种基于约束最大似然和模型平均的MR新方法,用于控制相关和不相关的多效性效应,能够提供更可靠和稳健的因果效应估计,比MR-Egger更为强大[24-25]。
敏感性分析
反向MR分析
将GO作为暴露因素,补体C5/C5a作为结局因素,进行反向MR分析,旨在排除补体C5/C5a与GO之间反向因果关联的影响,增加结果的可靠性。
共定位分析
采用colocR软件包进行共定位分析,测试补体C5与GO之间的关联是否由连锁不平衡驱动[29]。默认先验概率为:P1=1E-5,P2=1E-5,P12=1E-5。P1、P2和P12是预定义的概率,即测试区域中的SNP与补体C5表达、GO风险或两者有实质性联系。从共定位分析得出假设的后验概率(posterior probability of hypothesis,PPH),对应如下[30]:1)PPH0,SNP与任何一个表型都无关;2)PPH1,SNP与补体C5表达相关,但与GO风险无关;3)PPH2,SNP与GO风险相关,但与补体C5表达无关;4)PPH3,SNP与GO风险和补体C5表达均相关,但2种表型由不同的SNP驱动;5)PPH4,SNP与GO风险和补体C5表达均相关,且2种表型由共同的SNP驱动。共定位的显著性阈值设定为PPH4>0.80[31],与GO共定位被视为潜在的药物靶蛋白。
统计学处理
所有数据分析和可视化均在R 4.3.2软件中进行,本研究使用了R软件包“TwoSampleMR”“MRPRESSO”“coloc”“ggplot2”“locuscomparer”模块。检验水准α=0.05,P<0.05为差异有统计学意义。
结 果
MR分析结果
Wald比率法模型显示补体C5与GO存在显著正向因果关联(OR=4.109,95% CI 1.990~8.486,P=1.34×10-4),IVW模型显示C5a与G同样存在正向因果关联(OR=2.901, 95% CI 1.225~6.869,P=0.015);cML-MA法也观察到了一致的因果效应(P<0.05,表3)。
暴露 | 结局 | SNP | 方法 | OR | 95% CI | P |
---|---|---|---|---|---|---|
C5 | GO | 1 | Wald ratio | 4.109 | 1.990~8.486 | 1.34×10-4 |
1 | cML-MA | 4.108 | 2.078~8.120 | 4.81×10-5 | ||
C5a | GO | 2 | IVW | 2.901 | 1.225~6.869 | 0.015 |
2 | cML-MA | 2.755 | 1.132~6.705 | 0.025 |
敏感性分析结果
Steiger方向性检验显示:补体C5与GO(P=6.99×10-64)、C5a与GO(P=1.13×10-47)的因果关系均为正向。补体C5作为暴露时,仅有1个SNP纳入分析,未进行Cochran Q检验;C5a作为暴露时,Cochran Q检验未发现SNP间存在异质性(P=0.440)。
反向MR分析
以GO作为暴露因素,补体C5/C5a作为结局因素,最终纳入3个SNP用于反向MR分析(表4)。IVW分析结果提示GO与补体C5(P=0.063)和C5a(P=0.340)均不具有因果关系(表5)。
SNP | Chr: Pos | EA | OA | EAF | β | SE | P | F |
---|---|---|---|---|---|---|---|---|
rs12198492 | 6:33900632 | C | G | 0.059 | 0.581 | 0.094 | 5.53×10-10 | 38.203 |
rs143190692 | 7:9088452 | A | G | 0.012 | 0.984 | 0.180 | 4.39×10-8 | 29.884 |
rs1794530 | 6:32704023 | G | T | 0.101 | 0.784 | 0.068 | 1.70×10-30 | 132.927 |
暴露 | 结局 | SNP | 方法 | OR | 95% CI | P |
---|---|---|---|---|---|---|
GO | C5 | 3 | IVW | 0.979 | 0.957~1.001 | 0.063 |
3 | MR Egger | 0.999 | 0.858~1.162 | 0.988 | ||
3 | Weighted median | 0.978 | 0.955~1.002 | 0.072 | ||
GO | C5a | 3 | IVW | 1.014 | 0.985~1.044 | 0.340 |
3 | MR Egger | 0.954 | 0.799~1.140 | 0.697 | ||
3 | Weighted median | 1.011 | 0.979~1.044 | 0.502 |
共定位分析
共定位分析有强有力的证据(PPH4=0.81)表明补体C5和GO在给定的基因区域内(1 000 kb)共享同一个因果变异(rs7036980),加强了MR分析中确定的致病作用,提示补体C5在GO的发病机制中起重要作用,是GO的潜在药物靶点(图2)。

讨 论
本研究利用两样本MR方法,基于GWAS数据探究补体C5/C5a与GO的因果关联,发现补体C5/C5a与GO存在正向因果关联,共定位分析结果显示补体C5与GO在给定的遗传窗口内共享同一个SNP(rs7036980),预示补体C5有望成为GO早期生物标志物及治疗靶点。
迄今为止,在GO的治疗中出现了许多针对补体成分的创新治疗尝试[32],但仍然缺乏关于补体系统参与GO发病机制的全面研究。补体系统由Jules Bordet和Paul Ehrlich在19世纪后期首次描述为辅助系统,是一个由50多种血清蛋白组成的复杂系统,在先天免疫和适应性免疫中发挥重要作用[33-34],而免疫失调和炎症反应是GO发病机制的潜在介质,故补体系统可能参与GO的发病过程[35]。补体级联反应通过经典(抗体依赖性)、替代性(抗体非依赖性)和凝集素途径激活,3个级联反应汇聚产生补体 C3 转化酶,该酶切割补体C3进而产生补体C5转化酶,该转化酶将补体C5切割为C5a和C5b,这是补体级联反应的关键末端效应子成分,可导致细胞死亡[36]。C5是补体级联反应的核心成分,抑制补体级联反应中的这一关键步骤可防止关键末端片段(C5a和C5b)的形成[37]。C5a可能在炎症小体启动和激活中发挥重要作用,导致细胞焦亡和细胞死亡;C5b可能参与膜攻击复合物(C5b-9)的形成,导致细胞裂解和细胞死亡[38-39]。抑制这些C5介导的复合物活性可能具有减缓GO进展的潜力。Jaffe等[40]的一项II期临床试验表明特异性靶向补体C5的抑制剂Avacincaptad pegol(一种聚乙二醇化的RNA适配体)减缓了年龄相关性黄斑变性患者12个月内地理萎缩的进展,与C3抑制治疗的效果相比,补体C5抑制治疗效果起效更早,这可能是由于补体级联反应在C5水平的下游得到进一步抑制,且C5抑制理论上保留了C3活性,因此可能提供了额外的安全优势。此外,传统的药物开发的主要问题是由于II期和III期缺乏疗效而导致的失败[41],而MR研究的结果为蛋白质药物靶点与疾病之间的因果关系提供了令人信服的证据,增加了候选药物在III期临床试验中成功的可能性[42]。在一些情况下,MR研究在临床试验测试之前可预测候选药物的治疗效果。
本研究存在一定的局限性:首先,本研究使用的GWAS数据均来自于欧洲人群,结论推广到其他人群中需谨慎;其次,本研究中使用的pQTL是使用Olink和Somascan 2个不同的平台获得的,可能会引入异质性;最后,本研究仅从遗传学角度推断补体C5/C5a与GO存在因果关系,具体的生物学机制有待进一步研究。
综上所述,本研究利用MR方法证明了补体C5/C5a与GO存在正向因果关系,补体C5/C5a的高水平显著增加GO的发生风险;共定位分析结果表明补体C5是GO的潜在靶蛋白,通过靶向补体C5的抑制剂可以有效降低GO的发生风险。未来仍需要大样本、多中心的随机对照研究进一步证实补体C5抑制剂治疗GO的有效性和安全性。
Endokrine Orbitopathie [Graves’ ophthalmopathy]
[J]. Ophthalmologe, 2016, 113(4):349-364, 465-466. https://doi.org/10.1007/s00347- 016-0239-3.Graves’ ophthalmopathy
[J]. N Engl J Med, 2010, 362(8): 726-738. https://doi.org/10.1056/NEJMra0905750.The 2021 European Group on Graves’ orbitopathy (EUGOGO) clinical practice guidelines for the medical management of Graves’ orbitopathy
[J]. Eur J Endocrinol, 2021, 185(4): G43-G67. https://doi.org/10. 1530/EJE-21-0479.Current and promising therapies based on the pathogenesis of Graves’ ophthalmopathy
[J]. Front Pharmacol, 2023, 14: 1217253. https://doi.org/10.3389/fphar. 2023.1217253.Complement inhibitors for advanced dry age-related macular degeneration (geographic atrophy): some light at the end of the tunnel?
[J]. J Clin Med, 2023, 12(15): 5131. https://doi.org/10.3390/jcm12155131.Complement inhibitors for age-related macular degeneration
[J]. Cochrane Database Syst Rev, 2014, 2014(1):Complement activation levels are related to disease stage in AMD
[J]. Invest Ophthalmol Vis Sci, 2020, 61(3): 18. https://doi. org/10.1167/iovs.61.3.18.Early local activation of complement in aqueous humour of patients with age-related macular degeneration
[J]. Eye, 2019, 33(12): 1859-1864. https://doi. org/10.1038/s41433-019-0501-4.Local complement activation in aqueous humor in patients with age-related macular degeneration
[J]. Eye, 2017, 31(5): 810-813. https://doi.org/10.1038/eye.2016.328.Risk factors for incident age-related macular degeneration: pooled findings from 3 continents
[J]. Ophthalmology, 2004, 111(7): 1280-1287. https://doi.org/10.1016/j.ophtha.2003.11.010.Detecting and correcting for bias in Mendelian randomization analyses using Gene-by-Environment interactions
[J]. Int J Epidemiol, 2019, 48(3): 702-712. https://doi.org/10.1093/ije/dyy204.MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data
[J]. Int J Epidemiol, 2017, 46(6): 1734-1739. https://doi.org/10.1093/ije/dyx034.Sensitivity analyses for robust causal inference from Mendelian randomization analyses with multiple genetic variants
[J]. Epidemiology, 2017, 28(1): 30-42. https://doi.org/10.1097/EDE.0000000000000559.Plasma proteomic associations with genetics and health in the UK Biobank
[J]. Nature, 2023, 622(7982): 329-338. https://doi.org/10.1038/s41586- 023-06592-6.Large-scale integration of the plasma proteome with genetics and disease
[J]. Nat Genet, 2021, 53(12): 1712-1721. https://doi.org/10. 1038/s41588-021-00978-w.Circulating cytokines and venous thromboembolism: a bidirectional two-sample Mendelian randomization study
[J]. Thromb Haemost, 2024, 124(5): 471-481. https://doi.org/10.1055/s-0043-1777351.XGBoost regression of the most significant photoplethysmogram features for assessing vascular aging
[J]. IEEE J Biomed Health Inform, 2022, 26(7): 3354-3361. https://doi. org/10.1109/JBHI.2022.3151091.Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians
[J]. BMJ, 2018, 362:孟德尔随机化探究白细胞介素与便秘的因果效应
[J]. 西安交通大学学报(医学版), 2023, 44(5): 737-745. https://doi.org/10.7652/jdyxb202305012.胃食管反流病增加慢性阻塞性肺疾病的患病风险: 一项两样本双向孟德尔随机化研究
[J]. 重庆医科大学学报, 2023, 48(12): 1439-1445. https://doi.org/10. 13406/j.cnki.cyxb.003396.A large-scale plasma proteome Mendelian randomization study identifies novel causal plasma proteins related to primary biliary cholangitis
[J]. Front Immunol, 2023, 14: 1052616. https://doi.org/10.3389/fimmu. 2023. 1052616.Novel causal plasma proteins for hypothyroidism: a large-scale plasma proteome Mendelian randomization analysis
[J]. J Clin Endocrinol Metab, 2023, 108(2): 433-442. https://doi.org/10.1210/clinem/dgac575.Plasma proteins and onset of type 2 diabetes and diabetic complications: Proteome-wide Mendelian randomization and colocalization analyses
[J]. Cell Rep Med, 2023, 4(9): 101174. https://doi.org/10.1016/j.xcrm.2023.101174.Does co-localization analysis reinforce the results of Mendelian randomization?
[J/OL]. Brain, 2024, 147(1): e7-e8[Association between gut microbiota and preeclampsia-eclampsia: a two-sample Mendelian randomization study
[J]. BMC Med, 2022, 20(1): 443. https://doi. org/10.1186/s12916-022-02657-x.Calcium homeostasis and psychiatric disorders: a Mendelian randomization study
[J]. Nutrients, 2023, 15(18): 4051. https://doi.org/10.3390/nu15184051.Hypothyroidism has a protective causal association with hepatocellular carcinoma: a two-sample Mendelian randomization study
[J]. Front Endocrinol, 2022, 13: 987401. https://doi.org/10.3389/fendo.2022.987401.Phosphodiesterase and psychiatric disorders: a two-sample Mendelian randomization study
[J]. J Transl Med, 2023, 21(1): 560. https://doi.org/10. 1186/s12967-023-04368-0.Bayesian test for colocalisation between pairs of genetic association studies using summary statistics
[J/OL]. PLoS Genet, 2014, 10(5): e1004383[A fast and efficient colocalization algorithm for identifying shared genetic risk factors across multiple traits
[J]. Nat Commun, 2021, 12(1): 764. https://doi.org/10.1038/s41467-020-20885-8.Therapeutic targets for inflammatory bowel disease: proteome-wide Mendelian randomization and colocalization analyses
[J]. EBioMedicine, 2023, 89: 104494. https://doi.org/10.1016/j.ebiom.2023.104494.Clinical promise of next-generation complement therapeutics
[J]. Nat Rev Drug Discov, 2019, 18(9): 707-729. https://doi.org/10.1038/s41573-019-0031-6.New insights into the immune functions of complement
[J]. Nat Rev Immunol, 2019, 19(8): 503-516. https://doi.org/10.1038/s41577- 019-0168-x.New insights into the immune functions of podocytes: the role of complement
[J]. Mol Cell Pediatr, 2023, 10(1): 3. https://doi.org/10.1186/s40348- 023-00157-3.Rituximab for thyroid-associated ophthalmopathy
[J]. Cochrane Database Syst Rev, 2022, 6(6):Complement: tapping into new sites and effector systems
[J]. Nat Rev Immunol, 2014, 14(12): 811-820. https://doi.org/10.1038/nri3761.The NLRP3 inflammasome: molecular activation and regulation to therapeutics
[J]. Nat Rev Immunol, 2019, 19(8): 477-489. https://doi.org/10.1038/s41577- 019-0165-0.Local factor H production by human choroidal endothelial cells mitigates complement deposition: implications for macular degeneration
[J]. J Pathol, 2022, 257(1): 29-38. https://doi.org/10.1002/path. 5867.Retinal pigment epithelial cells mitigate the effects of complement attack by endocytosis of C5b-9
[J]. J Immunol, 2015, 195(7): 3382-3389. https://doi.org/10.4049/jimmunol.1500937.C5 inhibitor avacincaptad pegol for geographic atrophy due to age-related macular degeneration: a randomized pivotal phase 2/3 trial
[J]. Ophthalmology, 2021, 128(4): 576-586. https://doi.org/10.1016/j. ophtha.2020.08.027.Integrating genomics with biomarkers and therapeutic targets to invigorate cardiovascular drug development
[J]. Nat Rev Cardiol, 2021, 18(6): 435-453. https://doi.org/10.1038/s41569-020-00493-1.The support of human genetic evidence for approved drug indications
[J]. Nat Genet, 2015, 47(8): 856-860. https://doi.org/10.1038/ng.3314.作者声称无任何利益冲突。
朱敏, 吴冰萱, 陈赵昌辞, 陈海燕, 熊炜. 补体C5/C5a在Graves眼病诊断与治疗中的价值: 一项孟德尔随机化研究[J]. 中南大学学报(医学版), 2024, 49(10): 1633-1641. DOI:10.11817/j.issn.1672-7347.2024.240062
ZHU Min, WU Bingxuan, CHEN-ZHAO Changci, CHEN Haiyan, XIONG Wei. Complement C5/C5a in the diagnosis and treatment of Graves’ ophthalmopathy: A Mendelian randomized study[J]. Journal of Central South University. Medical Science, 2024, 49(10): 1633-1641. DOI:10.11817/j.issn.1672-7347. 2024.240062
http://xbyxb.csu.edu.cn/xbwk/fileup/PDF/2024101633.pdf
http://dx.chinadoi.cn/10.11817/j.issn.1672-7347.2024.240062