中国全科医学 ›› 2025, Vol. 28 ›› Issue (33): 4159-4165.DOI: 10.12114/j.issn.1007-9572.2024.0395
收稿日期:
2024-09-05
修回日期:
2025-01-15
出版日期:
2025-11-20
发布日期:
2025-09-17
通讯作者:
徐积兄
作者贡献:
张书岑提出主要研究目标,负责研究的构思与设计,研究的实施,进行数据的收集与整理,统计学处理,图、表的绘制与展示,撰写论文;徐积兄负责文章的质量控制与审查,对文章整体负责,监督管理。
基金资助:
Received:
2024-09-05
Revised:
2025-01-15
Published:
2025-11-20
Online:
2025-09-17
Contact:
XU Jixiong
摘要: 背景 根据现有的观察性研究结果,多发性硬化症(MS)与自身免疫性甲状腺疾病(AITD)的发病风险之间存在潜在关联,但结果仍存在争议。 目的 旨在采用孟德尔随机化(MR)方法评估MS与AITD之间的因果关系。 方法 本研究使用全基因组关联研究(GWAS)汇总MS(病例组=47 429,对照组=68 374)、自身免疫性甲状腺功能亢进症(病例组=1 991,对照组=305 175)、Graves病(病例组=3 176,对照组=409 005)和自身免疫性甲状腺炎(病例组=539,对照组=349 717)的数据,使用双向MR分析来评估可能的潜在因果关系,MR分析主要采用随机效应逆方差加权法,辅以MR-Egger回归法、加权中位数法、加权模型以及校正最大似然法-多效性分析、稳健校正轮廓评分和去偏反方差加权法等新的MR分析方法。此外,还使用了敏感性检验来确保结果的稳健性,如Cochran Q检验、MR-Presso、MR-Egger回归和留一法检验。 结果 MR结果显示MS与自身免疫性甲状腺炎存在因果关联(OR=1.144,95%CI=1.029~1.272,P=0.013,PFDR=0.039),遗传预测的MS可能会使自身免疫性甲状腺炎的风险增加14%,经错误发现率校正后的结果一致,且结果无异质性(Q-pval=0.356)和水平多效性(P=0.359),留一分析没有显示任何异常单核苷酸多态性值,且Outlier-corrected分析结果显示未检测到离群值。反向MR分析显示AITD与MS之间不存在因果关联。 结论 MS患者发生自身免疫性甲状腺炎的风险增高,但与自身免疫性甲状腺功能亢进症和Graves病的发病风险无关联。
中图分类号:
图1 MR基本假设示意图注:MS=多发性硬化症,AITD=自身免疫性甲状腺疾病,AIT=自身免疫性甲状腺炎,AIH=自身免疫性甲状腺功能亢进症,GD=Graves病,SNPs=单核苷酸多态性。
Figure 1 Schematic diagram of Mendelian randomization's three basic hypotheses
特征 | 数据库 | 种族 | 病例组(例) | 对照组(例) | 样本量(例) | 数据库网址 |
---|---|---|---|---|---|---|
多发性硬化症 | 国际多发性硬化症遗传学联盟 | 欧洲 | 47 429 | 68 374 | 115 803 | http://gwas.mrcieu.ac.uk/ |
自身免疫性甲状腺功能亢进症 | 芬兰数据库 | 欧洲 | 1 991 | 305 175 | 307 166 | www.finngen.fi/en |
Graves病 | 芬兰数据库 | 欧洲 | 3 176 | 409 005 | 412 181 | www.finngen.fi/en |
自身免疫性甲状腺炎 | 芬兰数据库 | 欧洲 | 539 | 349 717 | 350 256 | www.finngen.fi/en |
表1 暴露因素与结局变量的数据来源
Table 1 Data sources for exposure factors and outcome variables
特征 | 数据库 | 种族 | 病例组(例) | 对照组(例) | 样本量(例) | 数据库网址 |
---|---|---|---|---|---|---|
多发性硬化症 | 国际多发性硬化症遗传学联盟 | 欧洲 | 47 429 | 68 374 | 115 803 | http://gwas.mrcieu.ac.uk/ |
自身免疫性甲状腺功能亢进症 | 芬兰数据库 | 欧洲 | 1 991 | 305 175 | 307 166 | www.finngen.fi/en |
Graves病 | 芬兰数据库 | 欧洲 | 3 176 | 409 005 | 412 181 | www.finngen.fi/en |
自身免疫性甲状腺炎 | 芬兰数据库 | 欧洲 | 539 | 349 717 | 350 256 | www.finngen.fi/en |
暴露 | 结局 | Nsnp(个) | 方法 | β | SE | OR(95%CI) | P值 | PFDR值 | Q-pval值 | 水平多效性检验 | |
---|---|---|---|---|---|---|---|---|---|---|---|
PMR-Egger值 | PMR-Presso值 | ||||||||||
MS | AIH | 60 | 逆方差加权法 | 0.072 | 0.043 | 1.074(0.988~1.168) | 0.094 | 0.140 | 5.536×10-8 | 0.089 | — |
MR-Egger回归法 | -0.024 | 0.070 | 0.976(0.852~1.119) | 0.073 | 0.790 | ||||||
加权中位数法 | 0.076 | 0.054 | 1.079(0.972~1.199) | 0.155 | 0.232 | ||||||
加权模型 | 0.130 | 0.111 | 1.138(0.917~1.414) | 0.246 | 0.369 | ||||||
稳健校正轮廓评分 | 0.112 | 0.036 | 1.119(1.043~1.200) | 0.002 | 0.005 | ||||||
去偏反方差加权法 | 0.073 | 0.040 | 1.076(0.995~1.163) | 0.067 | 0.100 | ||||||
校正最大似然法-多效性分析 | 0.151 | 0.049 | 1.163(1.056~1.282) | 0.002 | 0.007 | ||||||
MS | GD | 57 | 逆方差加权法 | 0.038 | 0.037 | 1.039(0.966~1.116) | 0.303 | 0.303 | 1.739×10-10 | 0.262 | — |
MR-Egger回归法 | -0.016 | 0.060 | 0.984(0.875~1.107) | 0.790 | 0.790 | ||||||
加权中位数法 | 0.068 | 0.044 | 1.070(0.981~1.168) | 0.126 | 0.232 | ||||||
加权模型 | 0.108 | 0.071 | 1.114(0.968~1.281) | 0.137 | 0.369 | ||||||
稳健校正轮廓评分 | 0.080 | 0.029 | 1.083(1.023~1.146) | 0.006 | 0.009 | ||||||
去偏反方差加权法 | 0.039 | 0.034 | 1.039(0.973~1.110) | 0.254 | 0.254 | ||||||
校正最大似然法-多效性分析 | 0.104 | 0.046 | 1.110(1.015~1.214) | 0.023 | 0.023 | ||||||
MS | AIT | 62 | 逆方差加权法 | 0.134 | 0.054 | 1.144(1.029~1.272) | 0.013 | 0.039 | 0.356 | 0.359 | 0.349 |
MR-Egger回归法 | 0.068 | 0.090 | 1.070(0.897~1.277) | 0.454 | 0.790 | ||||||
加权中位数法 | 0.086 | 0.088 | 1.090(0.916~1.296) | 0.332 | 0.332 | ||||||
加权模型 | 0.031 | 0.083 | 1.031(0.876~1.213) | 0.714 | 0.714 | ||||||
稳健校正轮廓评分 | 0.131 | 0.055 | 1.139(1.023~1.275) | 0.018 | 0.018 | ||||||
去偏反方差加权法 | 0.137 | 0.054 | 1.146(1.031~1.275) | 0.012 | 0.036 | ||||||
校正最大似然法-多效性分析 | 0.131 | 0.056 | 1.140(1.022~1.273) | 0.019 | 0.023 |
表2 多发性硬化症与自身免疫性甲状腺疾病的正向MR分析
Table 2 Positive MR analysis of multiple sclerosis on autoimmune thyroid diseases
暴露 | 结局 | Nsnp(个) | 方法 | β | SE | OR(95%CI) | P值 | PFDR值 | Q-pval值 | 水平多效性检验 | |
---|---|---|---|---|---|---|---|---|---|---|---|
PMR-Egger值 | PMR-Presso值 | ||||||||||
MS | AIH | 60 | 逆方差加权法 | 0.072 | 0.043 | 1.074(0.988~1.168) | 0.094 | 0.140 | 5.536×10-8 | 0.089 | — |
MR-Egger回归法 | -0.024 | 0.070 | 0.976(0.852~1.119) | 0.073 | 0.790 | ||||||
加权中位数法 | 0.076 | 0.054 | 1.079(0.972~1.199) | 0.155 | 0.232 | ||||||
加权模型 | 0.130 | 0.111 | 1.138(0.917~1.414) | 0.246 | 0.369 | ||||||
稳健校正轮廓评分 | 0.112 | 0.036 | 1.119(1.043~1.200) | 0.002 | 0.005 | ||||||
去偏反方差加权法 | 0.073 | 0.040 | 1.076(0.995~1.163) | 0.067 | 0.100 | ||||||
校正最大似然法-多效性分析 | 0.151 | 0.049 | 1.163(1.056~1.282) | 0.002 | 0.007 | ||||||
MS | GD | 57 | 逆方差加权法 | 0.038 | 0.037 | 1.039(0.966~1.116) | 0.303 | 0.303 | 1.739×10-10 | 0.262 | — |
MR-Egger回归法 | -0.016 | 0.060 | 0.984(0.875~1.107) | 0.790 | 0.790 | ||||||
加权中位数法 | 0.068 | 0.044 | 1.070(0.981~1.168) | 0.126 | 0.232 | ||||||
加权模型 | 0.108 | 0.071 | 1.114(0.968~1.281) | 0.137 | 0.369 | ||||||
稳健校正轮廓评分 | 0.080 | 0.029 | 1.083(1.023~1.146) | 0.006 | 0.009 | ||||||
去偏反方差加权法 | 0.039 | 0.034 | 1.039(0.973~1.110) | 0.254 | 0.254 | ||||||
校正最大似然法-多效性分析 | 0.104 | 0.046 | 1.110(1.015~1.214) | 0.023 | 0.023 | ||||||
MS | AIT | 62 | 逆方差加权法 | 0.134 | 0.054 | 1.144(1.029~1.272) | 0.013 | 0.039 | 0.356 | 0.359 | 0.349 |
MR-Egger回归法 | 0.068 | 0.090 | 1.070(0.897~1.277) | 0.454 | 0.790 | ||||||
加权中位数法 | 0.086 | 0.088 | 1.090(0.916~1.296) | 0.332 | 0.332 | ||||||
加权模型 | 0.031 | 0.083 | 1.031(0.876~1.213) | 0.714 | 0.714 | ||||||
稳健校正轮廓评分 | 0.131 | 0.055 | 1.139(1.023~1.275) | 0.018 | 0.018 | ||||||
去偏反方差加权法 | 0.137 | 0.054 | 1.146(1.031~1.275) | 0.012 | 0.036 | ||||||
校正最大似然法-多效性分析 | 0.131 | 0.056 | 1.140(1.022~1.273) | 0.019 | 0.023 |
暴露 | 结局 | Nsnp(个) | 方法 | β | SE | OR(95%CI) | P值 | PFDR值 | Q-pval值 | 水平多效性检验 | |
---|---|---|---|---|---|---|---|---|---|---|---|
PMR-Egger值 | PMR-Presso值 | ||||||||||
AIH | MS | 51 | 逆方差加权法 | -0.006 | 0.017 | 0.994(0.961~1.029) | 0.731 | 0.741 | 0.088 | 0.095 | — |
MR-Egger回归法 | 0.056 | 0.040 | 1.058(0.977~1.145) | 0.171 | 0.171 | ||||||
加权中位数法 | 0.010 | 0.026 | 1.010(0.960~1.063) | 0.687 | 0.811 | ||||||
加权模型 | 0.069 | 0.036 | 1.071(0.999~1.148) | 0.059 | 0.176 | ||||||
稳健校正轮廓评分 | -0.007 | 0.018 | 0.993(0.959~1.028) | 0.697 | 0.846 | ||||||
去偏反方差加权法 | -0.006 | 0.018 | 0.994(0.959~1.030) | 0.733 | 0.741 | ||||||
校正最大似然法-多效性分析 | -0.010 | 0.025 | 0.990(0.943~1.038) | 0.669 | 0.694 | ||||||
GD | MS | 57 | 逆方差加权法 | -0.036 | 0.020 | 0.964(0.927~1.003) | 0.068 | 0.205 | 0.330 | 0.103 | — |
MR-Egger回归法 | -0.123 | 0.056 | 0.884(0.792~0.987) | 0.032 | 0.096 | ||||||
加权中位数法 | -0.039 | 0.029 | 0.962(0.908~1.018) | 0.181 | 0.543 | ||||||
加权模型 | 0.021 | 0.056 | 1.021(0.916~1.139) | 0.705 | 0.849 | ||||||
稳健校正轮廓评分 | -0.041 | 0.023 | 0.960(0.917~1.004) | 0.077 | 0.230 | ||||||
去偏反方差加权法 | -0.038 | 0.021 | 0.963(0.924~1.003) | 0.068 | 0.204 | ||||||
校正最大似然法-多效性分析 | -0.040 | 0.022 | 0.960(0.919~1.003) | 0.070 | 0.211 | ||||||
AIT | MS | 19 | 逆方差加权法 | 0.006 | 0.018 | 1.006(0.971~1.043) | 0.741 | 0.741 | 0.628 | 0.073 | 0.615 |
MR-Egger回归法 | 0.118 | 0.061 | 1.125(0.998~1.268) | 0.072 | 0.107 | ||||||
加权中位数法 | 0.006 | 0.027 | 1.006(0.955~1.061) | 0.811 | 0.811 | ||||||
加权模型 | 0.007 | 0.034 | 1.007(0.942~1.076) | 0.849 | 0.849 | ||||||
稳健校正轮廓评分 | 0.004 | 0.020 | 1.004(0.965~1.045) | 0.846 | 0.846 | ||||||
去偏反方差加权法 | 0.006 | 0.019 | 1.006(0.969~1.045) | 0.741 | 0.741 | ||||||
校正最大似然法-多效性分析 | 0.008 | 0.019 | 1.008(0.971~1.046) | 0.694 | 0.694 |
表3 多发性硬化症与自身免疫性甲状腺疾病的反向MR分析
Table 3 Reverse MR analysis of multiple sclerosis on autoimmune thyroid diseases
暴露 | 结局 | Nsnp(个) | 方法 | β | SE | OR(95%CI) | P值 | PFDR值 | Q-pval值 | 水平多效性检验 | |
---|---|---|---|---|---|---|---|---|---|---|---|
PMR-Egger值 | PMR-Presso值 | ||||||||||
AIH | MS | 51 | 逆方差加权法 | -0.006 | 0.017 | 0.994(0.961~1.029) | 0.731 | 0.741 | 0.088 | 0.095 | — |
MR-Egger回归法 | 0.056 | 0.040 | 1.058(0.977~1.145) | 0.171 | 0.171 | ||||||
加权中位数法 | 0.010 | 0.026 | 1.010(0.960~1.063) | 0.687 | 0.811 | ||||||
加权模型 | 0.069 | 0.036 | 1.071(0.999~1.148) | 0.059 | 0.176 | ||||||
稳健校正轮廓评分 | -0.007 | 0.018 | 0.993(0.959~1.028) | 0.697 | 0.846 | ||||||
去偏反方差加权法 | -0.006 | 0.018 | 0.994(0.959~1.030) | 0.733 | 0.741 | ||||||
校正最大似然法-多效性分析 | -0.010 | 0.025 | 0.990(0.943~1.038) | 0.669 | 0.694 | ||||||
GD | MS | 57 | 逆方差加权法 | -0.036 | 0.020 | 0.964(0.927~1.003) | 0.068 | 0.205 | 0.330 | 0.103 | — |
MR-Egger回归法 | -0.123 | 0.056 | 0.884(0.792~0.987) | 0.032 | 0.096 | ||||||
加权中位数法 | -0.039 | 0.029 | 0.962(0.908~1.018) | 0.181 | 0.543 | ||||||
加权模型 | 0.021 | 0.056 | 1.021(0.916~1.139) | 0.705 | 0.849 | ||||||
稳健校正轮廓评分 | -0.041 | 0.023 | 0.960(0.917~1.004) | 0.077 | 0.230 | ||||||
去偏反方差加权法 | -0.038 | 0.021 | 0.963(0.924~1.003) | 0.068 | 0.204 | ||||||
校正最大似然法-多效性分析 | -0.040 | 0.022 | 0.960(0.919~1.003) | 0.070 | 0.211 | ||||||
AIT | MS | 19 | 逆方差加权法 | 0.006 | 0.018 | 1.006(0.971~1.043) | 0.741 | 0.741 | 0.628 | 0.073 | 0.615 |
MR-Egger回归法 | 0.118 | 0.061 | 1.125(0.998~1.268) | 0.072 | 0.107 | ||||||
加权中位数法 | 0.006 | 0.027 | 1.006(0.955~1.061) | 0.811 | 0.811 | ||||||
加权模型 | 0.007 | 0.034 | 1.007(0.942~1.076) | 0.849 | 0.849 | ||||||
稳健校正轮廓评分 | 0.004 | 0.020 | 1.004(0.965~1.045) | 0.846 | 0.846 | ||||||
去偏反方差加权法 | 0.006 | 0.019 | 1.006(0.969~1.045) | 0.741 | 0.741 | ||||||
校正最大似然法-多效性分析 | 0.008 | 0.019 | 1.008(0.971~1.046) | 0.694 | 0.694 |
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] | |
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|
[31] |
|
[32] |
|
[33] |
|
[34] |
|
[35] |
|
[36] |
|
[37] |
|
[38] |
|
[39] |
|
[40] |
|
[41] |
|
[42] |
|
[43] |
CONSORTIUM I M S G,Wellcome Consortium 2,
|
[44] |
|
[45] |
|
[46] |
|
[1] | 马文源, 祁烁, 商建伟, 陈晓珩, 李哲, 李会龙, 户蕊, 李璐, 司新颖, 丁治国. 甲状腺功能正常的桥本甲状腺炎证候疗效评价量表的初步研制:基于专家咨询和临床调查[J]. 中国全科医学, 2025, 28(28): 3590-3600. |
[2] | 王爽, 吴树法, 令垚, 谭茜蔚, 曹汝岱, 曾慧婷, 孔丹莉, 丁元林, 于海兵. 基于代谢组学探究非脂质代谢物在肥胖与糖尿病视网膜病变间的中介作用:孟德尔随机化研究[J]. 中国全科医学, 2025, 28(21): 2625-2634. |
[3] | 谭文彬, 李佳, 刘明玉, 路永欣, 程雅欣. 神经系统疾病及相关治疗药物对骨质疏松症影响的研究进展[J]. 中国全科医学, 2025, 28(17): 2092-2100. |
[4] | 曾佳玲, 蒙艳, 邓婷婷, 李锦华, 赵萍. 2型糖尿病患者非酒精性脂肪肝与内脏脂肪面积及甲状腺结节的相关性研究[J]. 中国全科医学, 2025, 28(09): 1072-1083. |
[5] | 易芬, 王勇, 徐爱晖. 外周血绝对嗜酸性粒细胞计数水平对肺癌预后的评估价值研究[J]. 中国全科医学, 2024, 27(32): 4001-4008. |
[6] | 张超, 龙燕, 付溪娜. 双胎妊娠孕妇妊娠期甲状腺功能指标参考区间的探讨研究[J]. 中国全科医学, 2022, 25(27): 3404-3408. |
[7] | 加孜热亚·再依拿提,邢淑清,热那姑丽·艾克拜尔,张凯迪,张玉媛,邢蓓,王新玲,郭艳英. 围绝经期和绝经后女性甲状腺疾病现况及特征研究[J]. 中国全科医学, 2021, 24(3): 310-315. |
[8] | 邓臣前,陈树春. 2020年美国内分泌协会《临床实践指南:内分泌疾病患者的血脂管理》解读[J]. 中国全科医学, 2021, 24(21): 2646-2654. |
[9] | 李婧,邓垂文,徐腾达. 体检人群甲状腺自身抗体与抗谷氨酸脱羧酶抗体的相关性研究[J]. 中国全科医学, 2021, 24(18): 2285-2290. |
[10] | 施根灵,方慧,戴茜茜,高明喜,苏工. 正常甲状腺病态综合征对慢性心力衰竭患者远期预后的影响研究[J]. 中国全科医学, 2021, 24(12): 1476-1481. |
[11] | 张亚琴,陈明卫,王佑民,邓大同,陈婷婷,杨芳. 合肥地区人群甲状腺自身抗体特点及其影响因素研究[J]. 中国全科医学, 2019, 22(24): 2991-2995. |
[12] | 刘欣,王淑萍,吕娜,綦才辉,王振竞. 妊娠期糖尿病与妊娠期甲状腺疾病关系的研究进展[J]. 中国全科医学, 2018, 21(36): 4429-4433. |
[13] | 薛霞1,徐旭娟2*,顾志峰2,沈毅3,董庆2,周静1,吴爱娟2. 中文版甲状腺疾病患者生活质量问卷简明版的信效度检验[J]. 中国全科医学, 2017, 20(17): 2111-2117. |
[14] | 王永森*,赵琳. 内陆至沿海迁移人群甲状腺疾病的调查分析[J]. 中国全科医学, 2017, 20(12): 1486-1490. |
[15] | 赵红,王苏平,宋品,李莉,刘赞华,李淑敏. 双侧甲状腺切除加双侧颈部淋巴结清扫术后引发吉兰-巴雷综合征一例报告并文献复习[J]. 中国全科医学, 2015, 18(7): 850-852. |
阅读次数 | ||||||
全文 |
|
|||||
摘要 |
|
|||||