跳到主要内容

感谢您访问nature.com。您使用的是对CSS支持有限的浏览器版本。为了获得最好的体验,我们建议您使用最新的浏览器(或关闭Internet Explorer的兼容性模式)。同时,为了确保持续的支持,我们将在没有样式和JavaScript的情况下显示站点。

多祖先全基因组关联分析确定类风湿关节炎的新遗传机制

摘要

类风湿关节炎(RA)是一种高度遗传的复杂疾病,病因不明。RA的多祖先遗传研究有望提高遗传信号检测能力、精细定位分辨率和多基因风险评分(PRS)性能。在这里,我们提出了一项大规模的RA全基因组关联研究(GWAS),其中包括来自五个祖先群体的276,020个样本。我们进行了多祖先荟萃分析,确定了124个基因座(P< 5 × 10−8),其中34部为新作品。新位点上的候选基因表明了免疫系统的重要作用(例如,TNIP2而且TNFRSF11A)和关节组织(例如,WISP1)在类风湿关节炎病因学。多祖先精细映射通过生物学见解确定了假定的因果变异(例如,LEF1).此外,基于多血统GWAS的PRS优于基于单血统GWAS的PRS,并且在欧洲和东亚血统人群之间具有相当的性能。我们的研究为类风湿关节炎的病因学提供了一些见解,并提高了类风湿关节炎的遗传可预测性。

这是订阅内容的预览,通过你所在的机构访问

访问选项

租或购买这篇文章

只要这篇文章,只要你需要它

39.95美元

价格可能受当地税收的影响,在结账时计算

图1:GWAS参与者的不同祖先背景。
图2:精细映射分析确定了候选因果变量。
图3:的拼接和总表达IL6R共同导致类风湿性关节炎风险。
图4:的拼接PADI4会导致类风湿性关节炎风险
图5:S-LDSC分析显示,EUR-GWAS和EAS-GWAS的因果变异分布相似。
图6:IMPACT得分高的精细映射变体。
图7:功能注释和多祖先GWAS提高了PRS性能。

数据可用性

我们在GWAS目录中保存了6组汇总统计数据(所有RA和血清阳性RA的多血统、eu -GWAS和EAS-GWAS),登录号为GCST90132222、GCST90132223、GCST90132224、GCST90132225、GCST90132226和GCST90132227。我们沉积了PRS模型(multi-ancestry PRS with CD4+T细胞T-bet IMPACT注释)到多基因评分(PGS)目录;出版物ID是PGP000357,分数ID是PGS002745。概要统计数字和PRS模型也可在https://data.cyverse.org/dav-anon/iplant/home/kazuyoshiishigaki/ra_gwas/ra_gwas-10-28-2021.tar.英国生物库分析通过申请号进行。47821.

代码的可用性

有关守则可于本署网页(https://github.com/immunogenomics/RA_GWAS),并存档于Zenodo (https://doi.org/10.5281/zenodo.6999289).

参考文献

  1. Ajeganova, S. & Huizinga, T. W. J.血清阴性和血清阳性RA:相似但不同?风湿醇。11, 8-9(2015)。

    文章PubMed谷歌学者

  2. 冈田,Y.等。类风湿关节炎的遗传学有助于生物学和药物的发现。自然506, 376-381(2014)。

    文章中科院PubMed谷歌学者

  3. 石垣,K.等。在日本人群中进行的大规模全基因组关联研究确定了不同疾病的新易感性位点。Nat,麝猫。52, 669-679(2020)。

    文章中科院PubMed公共医学中心谷歌学者

  4. 麦格雷戈,A. J.等人。利用双胞胎数据描述类风湿关节炎的定量遗传贡献。关节炎感冒。43, 30-37(2000)。

    文章中科院PubMed谷歌学者

  5. 石垣,K.等。细胞特异性通路上的多基因负担是类风湿关节炎风险的基础。Nat,麝猫。49, 1120-1125(2017)。

    文章中科院PubMed谷歌学者

  6. Westra周宏儒。et al。精细定位和功能研究强调了类风湿性关节炎和1型糖尿病的潜在因果变异。Nat,麝猫。50, 1366-1374(2018)。

    文章中科院PubMed公共医学中心谷歌学者

  7. Trynka, G.等。染色质标记为精细映射复杂性状变异识别关键细胞类型。Nat,麝猫。45, 124-130(2013)。

    文章中科院PubMed谷歌学者

  8. Asgari, S.等。积极选择FBN1误义变异降低秘鲁个体的身高。自然582, 234-239(2020)。

    文章中科院PubMed公共医学中心谷歌学者

  9. SIGMA 2型糖尿病联盟。中低频变体的关联HNF1A患有2型糖尿病的拉丁裔人群《美国医学会杂志》311, 2305-2314(2014)。

    文章谷歌学者

  10. 毛奇,我等人。普通的格陵兰人TBC1D4变异导致肌肉胰岛素抵抗和2型糖尿病。自然512, 190-193(2014)。

    文章中科院PubMed谷歌学者

  11. 小山,S.等。人群特异性和跨血统全基因组分析确定不同和共享的冠状动脉疾病遗传风险位点。Nat,麝猫。52, 1169-1177(2020)。

    文章中科院PubMed谷歌学者

  12. 陈,M.-H。et al。来自全球5个人群的746,667人的跨种族和祖先特异性血细胞遗传学。细胞182, 1198 - 1213。e14灯头(2020)。

    文章中科院PubMed公共医学中心谷歌学者

  13. 黄,H.等。将炎症性肠病基因座精确定位到单变异分辨率。自然547, 173-178(2017)。

    文章中科院PubMed公共医学中心谷歌学者

  14. Laufer, v.a.等。非裔美国人类风湿关节炎易感性的遗传影响。嗡嗡声。摩尔,麝猫。28, 858-874(2019)。

    文章中科院PubMed谷歌学者

  15. 马丁,a.r.等人。目前多基因风险评分的临床应用可能会加剧健康差异。Nat,麝猫。51, 584-591(2019)。

    文章中科院PubMed公共医学中心谷歌学者

  16. 阮,Y.等。改进祖先多样性群体的多基因预测。Nat,麝猫。54, 573-580(2022)。

    文章中科院PubMed公共医学中心谷歌学者

  17. Márquez-Luna, C.等。多种族多基因风险评分改善了不同人群的风险预测。麝猫。论文。41, 811-823(2017)。

    文章PubMed公共医学中心谷歌学者

  18. 愣,R.-X。et al。类风湿关节炎新易感位点的鉴定。安。感冒。说。79, 1565-1571(2020)。

    文章中科院PubMed谷歌学者

  19. Kochi, Y.等人。一种功能变体FCRL3编码Fc受体样3,与类风湿关节炎和几种自身免疫相关。Nat,麝猫。37, 478-485(2005)。

    文章中科院PubMed公共医学中心谷歌学者

  20. 铃木,A.等。的功能单倍型PADI4,编码瓜氨酸酶肽基精氨酸脱亚胺酶4,与类风湿关节炎相关。Nat,麝猫。34, 395-402(2003)。

    文章中科院PubMed谷歌学者

  21. 冈田,Y.等。荟萃分析在日本人群中确定了9个与类风湿关节炎相关的新位点。Nat,麝猫。44, 511-516(2012)。

    文章中科院PubMed谷歌学者

  22. Diogo, D.等。TYK2蛋白质编码变体可以预防类风湿性关节炎和自身免疫,没有证据表明对非自身免疫复合物性状有主要的多效性作用。《公共科学图书馆•综合》10, e0122271(2015)。

    文章PubMed公共医学中心谷歌学者

  23. 特雷勒等人。类风湿关节炎放射损伤的遗传相关性:7项2775例全基因组相关性研究的荟萃分析《公共科学图书馆•综合》14, e0223246(2019)。

    文章中科院PubMed公共医学中心谷歌学者

  24. Márquez, A.等。四种自身免疫性疾病的免疫芯片数据的元分析揭示了新的单一疾病和跨表型关联。基因组医学。10, 97(2018)。

    文章PubMed公共医学中心谷歌学者

  25. 魏,W.-H。,Viatte, S., Merriman, T. R., Barton, A. & Worthington, J. Genotypic variability based association identifies novel non-additive lociDHCR7而且IRF4在血清阴性类风湿性关节炎中。科学。代表。7, 5261(2017)。

    文章PubMed公共医学中心谷歌学者

  26. Márquez, A.等。一项联合的大规模元分析确定COG6作为一种新的共同风险轨迹用于风湿性关节炎和系统性红斑狼疮。安。感冒。说。76, 286-294(2017)。

    文章PubMed谷歌学者

  27. Bossini-Castillo, L.等人。无瓜氨酸肽抗体的类风湿关节炎全基因组关联研究。安。感冒。说。74, e15(2015)。

    文章中科院PubMed谷歌学者

  28. 爱,S.等。高密度基因定位发现了新的类风湿关节炎易感位点。Nat,麝猫。44, 1336-1340(2012)。

    文章中科院PubMed公共医学中心谷歌学者

  29. 阿科斯塔-埃雷拉,等人。全基因组元分析揭示了共享的新位点在全身性血清阳性风湿病中。安。感冒。说。78, 311-319(2019)。

    文章中科院PubMed谷歌学者

  30. Frisell, T.等人。类风湿关节炎的家族风险和遗传性:类风湿因子/抗瓜氨酸蛋白抗体状态的作用、受影响亲属的数量和类型、性别和年龄。关节炎感冒。65, 2773-2782(2013)。

    文章中科院PubMed谷歌学者

  31. Padyukov, L.等。一项全基因组关联研究表明,acpa阳性和acpa阴性类风湿关节炎存在明显的相关性。安。感冒。说。70, 259-265(2011)。

    文章PubMed谷歌学者

  32. 斐努凯恩,H. K.等。利用全基因组关联汇总统计数据,通过功能注释划分遗传力。Nat,麝猫。47, 1228-1235(2015)。

    文章中科院PubMed公共医学中心谷歌学者

  33. 威康信托病例控制联盟。3种常见疾病14个位点关联信号的贝叶斯细化。Nat,麝猫。44, 1294-1301(2012)。

    文章谷歌学者

  34. Kichaev, G. & Pasaniuc, B.在跨种族精细测绘研究中利用功能注释数据。点。j .的嗡嗡声。麝猫。97, 260-271(2015)。

    文章中科院PubMed公共医学中心谷歌学者

  35. 尤尔什,J. C.等。人类造血在单细胞和单变异分辨率的问题。Nat,麝猫。51, 683-693(2019)。

    文章中科院PubMed公共医学中心谷歌学者

  36. 傅,W.等。一个多重冗余的基因开关“锁定”了调节性T细胞的转录特征。Immunol Nat。13, 972-980(2012)。

    文章中科院PubMed公共医学中心谷歌学者

  37. 陈,L.等。人类免疫细胞表观遗传和转录变异的遗传驱动因素。细胞167, 1398 - 1414。e24(2016)。

    文章中科院PubMed公共医学中心谷歌学者

  38. Giambartolomei, C.等。使用汇总统计对遗传关联研究之间共定位的贝叶斯检验。公共科学图书馆麝猫。10, e1004383(2014)。

    文章PubMed公共医学中心谷歌学者

  39. Ferreira, R. C.等。功能IL6R358Ala等位基因破坏经典IL-6受体信号,影响多种炎症疾病的风险。公共科学图书馆麝猫。9, e1003444(2013)。

    文章中科院PubMed公共医学中心谷歌学者

  40. 冈田,Y.等。mirna -靶基因网络对人类复杂性状遗传的重要影响科学。代表。6, 22223(2016)。

    文章中科院PubMed公共医学中心谷歌学者

  41. Schellekens, g.a., de Jong, b.a.w, van den Hoogen, f.h.j., van de Putte, l.b.a.和van Venrooij, w.j.瓜氨酸是类风湿关节炎特异性自身抗体识别的抗原决定因子的重要组成部分。j .中国。投资。101, 273-281(1998)。

    文章中科院PubMed公共医学中心谷歌学者

  42. 铃木,A.等。肽基精氨酸脱亚胺酶4型敲除小鼠实验性自身免疫性关节炎严重程度降低。BMC Musculoskelet。Disord。17, 205(2016)。

    文章PubMed公共医学中心谷歌学者

  43. Seri, Y.等人。肽精氨酸脱亚胺酶4型在葡萄糖-6-磷酸异构酶诱导的关节炎模型中,缺乏葡萄糖-6-磷酸异构酶降低了关节炎的严重程度。科学。代表。5, 13041(2015)。

    文章中科院PubMed公共医学中心谷歌学者

  44. 有田,K.等。Ca的结构基础2 +-诱导人类PAD4的激活。Nat。结构。摩尔。杂志。11, 777-783(2004)。

    文章中科院PubMed谷歌学者

  45. 南达,S. K.等。ABIN2功能通过tpl2独立机制抑制dss诱导的结肠炎。j . Immunol。201, 3373-3382(2018)。

    文章中科院PubMed谷歌学者

  46. Matmati, M.等人。髓系细胞A20 (TNFAIP3)缺失引发类似类风湿性关节炎的糜糜性多发性关节炎。Nat,麝猫。43, 908-912(2011)。

    文章中科院PubMed谷歌学者

  47. GTEx财团。人类组织中基因表达的遗传效应。自然550, 204-213(2017)。

    文章公共医学中心谷歌学者

  48. 朱,Z.等。整合GWAS和eQTL研究的总结数据预测复杂的性状基因靶点。Nat,麝猫。48, 481-487(2016)。

    文章中科院PubMed谷歌学者

  49. 詹姆斯,E. A.等人。类风湿性关节炎中瓜氨酸特异性Th1细胞增加,其频率受疾病持续时间和治疗的影响。关节炎Rheumatol。66, 1712-1722(2014)。

    文章中科院PubMed公共医学中心谷歌学者

  50. Takayanagi, H.等人。RANKL通过c- fos依赖性的诱导维持骨稳态干扰素-β自然416, 744-749(2002)。

    文章中科院PubMed谷歌学者

  51. Takeuchi, T.等。抗rankl抗体denosumab对用常规合成疾病修饰抗风湿药物治疗的类风湿关节炎患者关节结构损伤的影响(理想的研究):随机、双盲、安慰剂对照阶段安。感冒。说。78, 899-907(2019)。

    文章中科院PubMed谷歌学者

  52. Nakatsuka, K., Nishizawa, Y. & Ralston, S. H.骨早期发病Paget病的表型特征TNFRSF11A基因。挖骨者。Res。18, 1381-1385(2003)。

    文章中科院PubMed谷歌学者

  53. 格雷尼,m.m.等人。人类破骨细胞不足骨石化伴低γ球蛋白血症,由于TNFRSF11A排名)突变。点。j .的嗡嗡声。麝猫。83, 64-76(2008)。

    文章中科院PubMed公共医学中心谷歌学者

  54. French, D. M.等人。WISP-1是一种在骨骼发育和骨折修复过程中表达的成骨细胞调节因子。点。j .分册。165, 855-867(2004)。

    文章中科院PubMed公共医学中心谷歌学者

  55. Maeda, A.等。wnt1诱导的分泌蛋白-1 (WISP1),一种新的骨转换和Wnt信号的调节因子。生物。化学。290, 14004-14018(2015)。

    文章中科院PubMed公共医学中心谷歌学者

  56. 张,F.等。通过整合单细胞转录组学和大规模细胞术定义类风湿关节炎关节滑膜组织中的炎症细胞状态。Immunol Nat。20., 928-942(2019)。

    文章PubMed公共医学中心谷歌学者

  57. 拉莫斯,m.i.p.等。缺乏fms样酪氨酸激酶3配体(Flt3L)信号可以预防胶原诱导的关节炎。安。感冒。说。74, 211-219(2015)。

    文章中科院PubMed谷歌学者

  58. Saevarsdottir, S.等人。FLT3停止突变增加FLT3配体水平和自身免疫性甲状腺疾病的风险。自然584, 619-623(2020)。

    文章中科院PubMed谷歌学者

  59. 茂木,T.等。稀有编码变异的鉴定TYK2对日本人群类风湿性关节炎的保护作用以及它们对细胞因子信号的影响。安。感冒。说。78, 1062-1069(2019)。

    文章中科院PubMed谷歌学者

  60. Brown, B. C.,亚洲遗传流行病学网络2型糖尿病联盟,Ye, C. J., Price, A. L. & Zaitlen, N.根据汇总统计的跨种族遗传相关性估计。点。j .的嗡嗡声。麝猫99, 76-88(2016)。

  61. Amariuta, T.等。通过对预测的细胞类型特异性调控元件中的变异进行优先排序,提高多基因风险评分的跨祖先可移植性。Nat,麝猫。52, 1346-1354(2020)。

    文章中科院PubMed公共医学中心谷歌学者

  62. Amariuta, T.等。影响:从结合转录因子的表观基因组推断的细胞状态特异性调控元件的基因组注释。点。j .的嗡嗡声。麝猫。104, 879-895(2019)。

    文章中科院PubMed公共医学中心谷歌学者

  63. 表观基因组学联盟路线图。111个参考人类表观基因组的综合分析。自然518, 317-330(2015)。

    文章公共医学中心谷歌学者

  64. 斐努凯恩,H. K.等。特异性表达基因的遗传力富集可识别疾病相关组织和细胞类型。Nat,麝猫。50, 621-629(2018)。

    文章中科院PubMed公共医学中心谷歌学者

  65. Chun, S.等人。在三种主要免疫细胞类型中eqtl和自身免疫疾病相关位点共享遗传效应的统计证据有限。Nat,麝猫。49, 600-605(2017)。

    文章中科院PubMed公共医学中心谷歌学者

  66. van der Wijst, m.g.p.等。单细胞RNA测序鉴定细胞类型特异性顺式eqtl和共表达qtl。Nat,麝猫。50, 493-497(2018)。

    文章PubMed公共医学中心谷歌学者

  67. Kumasaka, N, Knights, A. J. & Gaffney, D. J.开放染色质区域之间假定因果相互作用的高分辨率遗传图谱。Nat,麝猫。51, 128-137(2019)。

    文章中科院PubMed谷歌学者

  68. 古铁雷斯-阿塞勒斯等人。在HLA和其他自身免疫位点中,等位基因特异性表达在T细胞激活过程中动态变化。Nat,麝猫。52, 247-253(2020)。

    文章中科院PubMed公共医学中心谷歌学者

  69. Baglaenko, Y., Macfarlane, D., Marson, A., Nigrovic, P. A. & Raychaudhuri, S.基因组编辑定义风湿疾病风险位点和变异的功能。风湿醇。17, 462-474(2021)。

    文章中科院PubMed谷歌学者

  70. 阿内特,f.c.等。美国风湿病协会1987年修订了类风湿关节炎的分类标准。关节炎感冒。31, 315-324(1988)。

    文章中科院PubMed谷歌学者

  71. Aletaha, D.等。2010年类风湿关节炎分类标准:美国风湿病学会/欧洲抗风湿病联盟合作倡议。关节炎感冒。62, 2569-2581(2010)。

    文章PubMed谷歌学者

  72. 秋山等人。描述日本人群中罕见和低频的身高相关变异。Commun Nat。10, 4393(2019)。

    文章中科院PubMed公共医学中心谷歌学者

  73. 张伟等。通过对35人的全基因组测序,可以了解韩国人的遗传结构。BMC Bioinf。15, s6(2014)。

    文章谷歌学者

  74. Lan, T.等。90个汉族人基因组深度全基因组测序。GigaScience6, gix067(2017)。

    文章谷歌学者

  75. Wong L.-P。et al。100个东南亚马来人的深度全基因组测序。点。j .的嗡嗡声。麝猫。92, 52-66(2013)。

    文章中科院PubMed公共医学中心谷歌学者

  76. 遗传流行病学研究中错误发现概率的贝叶斯度量。点。j .的嗡嗡声。麝猫。81, 208-227(2007)。

    文章中科院PubMed公共医学中心谷歌学者

  77. 刘志刚,王志刚,王志刚。髓系细胞中PAD表达与活性的比较分析。j . Immunol。198, 75.18(2017)。

    谷歌学者

  78. 皮切利等人。使用Smart-seq2从单细胞获得全长RNA-seq。Protoc Nat。9, 171-181(2014)。

    文章中科院PubMed谷歌学者

下载参考

确认

我们感谢马来西亚卫生部局长支持南亚(SAS)人口中描述的工作:马来西亚类风湿性关节炎流行病学调查(MyEIRA)研究。MyEIRA研究由马来西亚卫生部(NMRR-08-820-1975)和瑞典国家研究委员会(DNR-348-2009-6468)资助。GENRA研究和CARDERA遗传队列基因分型是由Versus Arthritis(授予文献19739到I.C.S.)资助的。护士健康研究(NHS队列)由美国国立卫生研究院(NIH)资助(R01 AR049880, UM1 CA186107, R01 CA49449, U01 CA176726和R01 CA67262)。瑞典的EIRA研究得到了瑞典研究委员会(对l.k., L.P.和L.A.)的支持。S.S.部分得到Mochida医学和药学研究纪念基金会、Kanae医学科学促进基金会、Astellas代谢紊乱研究基金会、JCR促进基础风湿病赠款和Manabe过敏性和风湿病奖学金赠款的支持。I.C.S.由国家卫生与保健研究所(NIHR)高级研究奖学金(授予参考NIHR300826)资助。本文仅代表作者个人观点,并不代表国家卫生研究院或卫生与社会保障部。K.A.S.由谢尔曼家族基因组医学主席、加拿大卫生研究基金会拨款(FDN 148457)和安大略省研究基金(RE-09-090)和加拿大创新基金会(33374)资助。s - c.b.由教育部基础科学研究基金资助(NRF- 2021r1a6a1a03038899)。 R.P.K. and J.C.E. are funded by NIH (UL1 TR003096). C.M.L. is partly funded by the NIHR Maudsley Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King’s College London. T. Arayssi was partially supported by the National Priorities Research Program (grant 4-344-3-105 from the Qatar National Research Fund, a member of Qatar Foundation). M. Kerick and J.M. are funded by Rheumatology Cooperative Research Thematic Network program RD16/0012/0013 from the Instituto de Salud Carlos III (Spanish Ministry of Science and Innovation). Y.O. is funded by JSPS KAKENHI (19H01021 and 20K21834), AMED (JP21km0405211, JP21ek0109413, JP21ek0410075, JP21gm4010006 and JP21km0405217), JST Moonshot R&D (JPMJMS2021 and JPMJMS2024), Takeda Science Foundation, and the Bioinformatics Initiative of Osaka University Graduate School of Medicine. Y. Kochi is funded by grants from Nanken-Kyoten, TMDU and Medical Research Center Initiative for High Depth Omics. S.R. is supported by UH2AR067677, U01HG009379, R01AR063759 and U01HG012009.

作者信息

作者及隶属关系

作者

财团

贡献

石垣、s.s.、c.t.、Y.O.和S.R.构思并设计了这项研究。石垣圭理根据s.s.、c.t.、y.l.、Y.O.和s.r.k .石垣圭理进行元分析,所有的GWAS下游分析都是在s.s.、c.t.、T. Amariuta、y.l.、Y.O.和s.r.k. Yamaguchi和Y. Kochi的帮助下进行的PADI4长读测序和PADI4sQTL分析。M. Koido, k.t., Y. Kamatani和C.T.对特定人群参考面板的构建做出了贡献。k .石垣岛,砂岩,成分,k, V.A.L,智能卡,S.V,一般,有机,G.X, J.Z,中情局,工作者,R.J.C。·,m . Kerick F.M。m . Traylor C.M.L, H.X, R.S。t . Arayssi赵硕,,L.K Y.O.和S.R. GWAS进行分析。成分,C.L.T V.A.L, S.V。m .高桥X.W, L.L。T.L,一般,A.B。G.O,有机,伺服电动机,K.P.L, R.J.C, E.W.K。k .松尾F.M,东南部,H.X。k . Ikari P.K.G,共和国l.p. Y.O.和导致的基因实验。成分,k。,C.L.T V.A.L,智能卡,S.V。好的,点,M.H,我方。m . Hammoudeh S.A.E, B.K.M, H.H。台北,I.W.U, X.W, L.L。T.L,一般,A.B。G.O, S.M.M.V, A.J.M, S.H, H.T,外星人,其子as, Y.M,伺服电动机,山本贤G.X, J.Z,中情局,工作者,g.w. I.v.d.H.-B。李鸿源K.P.L。,R.J.C H.-S.L, S.-Y.B。·,N.d.V,洛杉矶,S.R.-D。,E.W.K., S.-C.B., R.P.K., J.C.E., X.M., T.H., P.D., M.S., M. Kerick, J.C.D., The BioBank Japan Project, K. Matsuda, K. Matsuo, T.M., F.M., K.F., Y.T., A.K., C.M.L., S.E., H.X., R.S., T. Arayssi, K. Ikari, M. Harigai, P.K.G., Kazuhiko Yamamoto, S.L.B., L.P., J.M., L.K., Y.O. and S.R. contributed to the collection of samples and management of genotype data and clinical information. J.C.D.’s involvement in this project was primarily as faculty at Vanderbilt University Medical Center prior to joining the NIH.

相应的作者

对应到Yukinori冈田克也Soumya Raychaudhuri

道德声明

相互竞争的利益

作者声明没有利益竞争。

同行评审

同行评审信息

自然遗传学感谢匿名审稿人对本工作的同行评议所作的贡献。同行评审报告是可用的。

额外的信息

出版商的注意施普林格自然对出版的地图和机构从属关系中的管辖权主张保持中立。

扩展数据

扩展数据图1 1KG Phase 3 PCA和UMAP图

一个,我们将每个个体的估算基因型投影到PC空间中,该空间使用1KG阶段3的所有个体计算。b,我们进一步使用所有GWAS样本和所有1KG阶段3的个体的前20个PC评分进行UMAP分析。我们绘制了1KG Phase 3中的所有样本,或分别绘制了每个祖先中的样本。所有GWAS样本的UMAP图如图所示。1 c

图2血清阳性和阴性RA的效应量异质性。

一个,绘制血清RA阳性和阴性的比值比及其95%置信区间。在122个显著常染色体位点的先导变异中,我们绘制了118个1KG期3期EUR常见变异(MAF > 0.01)。我们介绍了多血统GWAS和EUR-GWAS的结果。通过Cochran’s Q检验评估多血统GWAS结果中的效应量异质性(Phet).给出了Pearson相关数据。b,血清阳性和阴性RA之间的效应大小大小的差异(delta)。正δ值表明血清阳性RA比血清阴性RA的效应值更大。我们用以下公式定义delta的标准误差(s.e.):\ \√{S.E._{血清阳性}^2 + S.E._{血清阴性}^2}\).我们提供2 × s.e.的delta作为误差条。在血清阳性RA GWAS中,我们有27,448例和240,149例对照;在血清阴性RA GWAS中,我们有4515例病例和240149例对照。

扩展数据图3在18个血细胞群的开放染色质区域富集高PIP变体。

高浓度富集皮普用gchromVAR软件分析了18个血细胞群开放染色质区域内的变异。水平虚线表示Bonferroni校正P值阈值(0.05/18 = 0.0027)。P通过在gchromVAR中进行富集试验来估计值。

扩展数据图4条件分析结果IL2RA而且TNFAIP3位点。

一个b,在每个队列中进行条件分析,并使用逆方差加权固定效应模型对结果进行meta分析(一个IL2RA轨迹;bTNFAIP3轨迹)。我们使用多祖先GWAS结果。调查结果TYK2轨迹在扩展数据图中提供。6.LD的变体与先导变体(r2> 0.6(欧元和EAS祖先)在每一轮条件分析中用不同的颜色突出显示。

扩展数据图5单倍型水平的关联检验IL6R轨迹。

单倍型水平的关联结果在IL6R轨迹。我们对单倍型的定义如表所示,并使用分阶段估算的基因型数据估计了四种单倍型的剂量。我们在每个EUR队列中进行了多元逻辑回归检验(ngydF4y2Ba= 97,173),结果采用逆方差加权固定效应模型进行meta分析。效应量估计和2倍s.e.显示在右侧面板中。

扩展数据图6在TYK2轨迹。

一个b,在每个队列中进行条件分析,并使用逆方差加权固定效应模型对结果进行meta分析(一个EUR-GWAS;bEAS-GWAS)。LD的变体与先导变体(r2在每一轮条件分析中,> 0.6(欧元或EAS祖先)用不同的颜色突出显示。

扩展数据图7使用组蛋白标记注释的S-LDSC分析。

一个给出了396个组蛋白标记注释解释的遗传比例估计值及其95%置信区间。置信区间和P通过LDSC软件中实现的block-jackknife(单侧检验)估计非负tau值。当遗传力富集显著时(P< 0.05/396 = 1.3 ×10−4),该注释是由GWAS的类型着色的。bPMA-I中最佳组蛋白标记注释(H3K4me1)刺激了初级CD4+T细胞)和最好的IMPACT注释(CD4+T细胞T-bet)在S-LDSC分析中联合建模。以总变异遗传力归一化后的Tau估计值(每个注释中的每个变异遗传力)为中心,以其95%置信区间为误差条。采用EUR-GWAS结果。cP值由LDSC软件实现的block-jackknife计算,并表示每个注释的非负tau(每变体遗传力)的显著性(单侧检验)。每个组蛋白标记注释都根据其细胞类型类别着色。水平虚线表示bonferroni校正P-value threshold (0.05/396 = 1.3 ×10−4).上面的面板显示了在不控制最佳IMPACT注释(CD4+T细胞T-bet),底部面板显示了控制该效应的结果。

扩展数据图8以前的RA GWAS与此GWAS的PRS性能比较。

负债规模R2未纳入之前RA GWAS的15个队列。我们使用了冈田报道的多血统GWAS结果.和这个GWAS开发PRS模型。我们使用LOCO方法对基于该GWAS的PRS模型进行评估。通过双面配对Wilcoxon文本(ngydF4y2Ba= 15)。在每个箱线图中,水平线反映中位数,每个箱线图的顶部和底部反映四分位间距(IQR),胡须反映每组内的最大值和最小值,距离铰链不超过1.5 × IQR。

扩展数据图9不同祖先群体的PRS性能。

PRS绩效(责任规模R2),以显示每组PRS模型的组合及应用PRS的群组(ngydF4y2Ba左起分别为25、8、4)。采用双面配对Wilcoxon检验分析成绩最佳组与次优组的差异。在每个箱线图中,水平线反映中位数,每个箱线图的顶部和底部反映四分位间距(IQR),胡须反映每组内的最大值和最小值,距离铰链不超过1.5 × IQR。我们在适用的地方使用了LOCO方法。

扩展数据图10病例与对照组之间的PRS分布差异。

病例和对照组之间的PRS分布差异。带有CD4的多血统PRS+T细胞T-bet IMPACT注释。我们使用了LOCO方法。在每个队列中,使用对照样本的平均值和s.d.对PRS进行缩放,并在一个祖先组的队列中合并个体水平数据。

补充信息

补充信息

补充说明,补充图1-11和补充数据。

报告总结

同行评审文件

补充表

补充表1-14。

权利和权限

根据与作者或其他权利持有人签订的出版协议,自然或其许可方(例如,社会或其他合作伙伴)对本文拥有排他性权利;作者对这篇文章接受的手稿版本的自我存档仅受此类出版协议的条款和适用法律的约束。

转载及权限

关于本文

通过CrossMark验证货币和真实性

引用本文

石垣,K,樱,S,寺尾,C。et al。多祖先全基因组关联分析确定类风湿关节炎的新遗传机制。Nat麝猫54, 1640-1651(2022)。https://doi.org/10.1038/s41588-022-01213-w

下载引用

  • 收到了

  • 接受

  • 发表

  • 发行日期

  • DOIhttps://doi.org/10.1038/s41588-022-01213-w

搜索

快速链接

自然简报

报名参加自然简报时事通讯-什么重要的科学,免费到您的收件箱每天。

获取当天最重要的科学故事,免费在您的收件箱。 注册《自然简报》
Baidu
map