Analysis of causal relationship between gut microbiota and hyperthyroidism based on Mendelian randomization
-
摘要:
目的 探索肠道菌群(GM)与甲状腺功能亢进症(以下简称甲亢)的相关性。 方法 以MiBioGen数据库的肠道菌群为暴露因素,以甲状腺功能亢进症为结局变量。从IEU OpenGWAS数据库下载甲状腺功能亢进症的UKB数据,以逆方差加权法(IVW)为主,辅助以MR-Egger回归和加权众数法(WM)等进行双样本孟德尔随机化分析,采用β值、OR值等指标评估GM与甲亢之间的因果关系,并进行异质性检验、多效性检验和敏感性分析,然后验证MR分析结果的稳定性和可靠性。 结果 经孟德尔随机化(MR)分析,发现在119个属水平的肠道菌群中,与甲亢有关且呈负相关的肠道细菌分别为戴阿利斯特杆菌属(OR=0.998,95% CI: 0.996~1.000,P=0.037,β=-0.002)、肠杆菌科的志贺菌属(OR=0.998,95% CI: 0.995~1.000,P=0.026,β=-0.002)和瘤胃球菌科UCG014属(OR=0.998,95% CI: 0.996~1.000,P=0.032,β=-0.002),其β值均<0,OR值均<1且P<0.05,说明以上3个菌属可能是患者发生甲亢风险的保护性因素。异质性检验(P>0.05)、水平基因多效性检验(P>0.05)和敏感性分析结果证明结果可靠。 结论 戴阿利斯特杆菌属、肠杆菌科的志贺菌属、瘤胃球菌科UCG014属可能减少甲亢发病风险,初步发现了肠道菌群和甲亢间的潜在因果关系。 Abstract:Objective To investigate the relationship between human gut microbiota and hyperthyroidism. Methods The gut microbiota of the MiBioGen database was used as the exposure factor, and hyperthyroidism was used as the outcome variable. The GWAS data of hyperthyroidism were downloaded from the IEU Open GWAS database. Two-sample Mendelian randomization analysis was performed using regression models such as inverse variance weighting method, MR-Egger regression and weighted median method. The causal relationship between GM and intestinal flora was evaluated by β-value and odds ratio of effect index and so on. The stability and reliability of the results were verified by leave-one-out method, heterogeneity test and horizontal gene pleiotropic test. Results Mendelian randomization (MR) analysis revealed that, among 119 gut microbiota at the genus level, three bacterial taxa showed a negative association with hyperthyroidism: Dialister (OR=0.998, 95% CI: 0.996-1.000, P=0.037, β=-0.002), Escherichia. Shigella (OR=0.998, 95% CI: 0.996-1.000, P=0.032, β=-0.002), and Ruminococcaceae UCG014 (OR=0.998, 95% CI: 0.996-1.000, P=0.032, β=-0.002). All three genera exhibited β < 0, OR < 1, and P < 0.05, suggesting they may serve as protective factors against hyperthyroidism risk. The results of heterogeneity test (P>0.05), horizontal gene pleiotropic test (P>0.05) and sensitivity analysis proved that the formal results were reliable. Conclusion Dialister, Escherichia. Shigella and Ruminococcaceae UCG014 can reduce the risk of hyperthyroidism, the potential relationship between gut microbiota and hyperthyroidism is found preliminarily. -
表 1 筛选后代表性工具变量SNPs
Table 1. Representative instrumental variable SNPs after screening
工具变量(SNPs) 其他效应位点 效应位点 暴露因素 R2 F rs34583783 T G Actinomyces 0.002 22.23 rs2715439 T C Actinomyces 0.002 20.51 rs7680684 T C Adlercreutzia 0.002 24.36 rs2147798 G C Adlercreutzia 0.002 23.07 rs11729256 C T Akkermansia 0.002 24.97 rs4936098 G A Akkermansia 0.001 22.81 rs67281112 C G Alistipes 0.001 21.18 rs67705352 G T Alistipes 0.001 23.05 rs11769002 A G Alistipes 0.001 23.36 表 2 孟德尔随机化分析后获得3个菌属的分析结果
Table 2. Mendelian randomization analysis yielded results for three bacterial genera
肠道菌属 方法 工具变量数量 β P值 OR(95% CI) Dialister Inverse variance weighted 12 -0.002 0.037 0.998(0.996~1.000) MR Egger 12 -0.009 0.057 0.991(0.982~0.999) Weighted median 12 -0.002 0.185 0.998(0.995~1.001) Simple mode 12 -0.002 0.400 0.997(0.993~1.002) Weighted mode 12 0.000 0.944 1.000(0.995~1.005) Escherichia. Shigella Inverse variance weighted 15 -0.002 0.026 0.998(0.995~1.00) MR Egger 15 -0.003 0.311 0.997(0.990~1.002) Weighted median 15 -0.002 0.158 0.998(0.995~1.000) Simple mode 15 -0.002 0.392 0.998(0.993~1.002) Weighted mode 15 -0.002 0.366 0.998(0.993~1.002) RuminococcaceaeUCG014 Inverse variance weighted 18 -0.002 0.032 0.998(0.996~1.000) MR Egger 18 0.003 0.425 1.000(0.996~1.008) Weighted median 18 -0.002 0.185 0.998(0.995~1.000) Simple mode 18 -0.002 0.355 0.998(0.993~1.002) Weighted mode 18 -0.002 0.425 0.998(0.994~1.002) 表 3 甲亢相关肠道细菌的异质性和多效性分析结果
Table 3. Heterogeneity and pleiotropy analysis of GM related with hyperthyroidism
肠道菌属 方法 异质性检验 多效性检验 Q P值 RSS P值 Egger_截距 P值 Dialister Inverse variance weighted 10.546 0.482 MR Egger 7.680 0.660 0.000 0.121 MR-PRESSO 12.580 0.515 Escherichia. Shigella Inverse variance weighted 18.673 0.178 MR Egger 18.524 0.137 0.000 0.751 MR-PRESSO 21.445 0.197 RuminococcaceaeUCG014 Inverse variance weighted 25.318 0.088 MR Egger 21.615 0.156 0.000 0.117 MR-PRESSO 28.091 0.120 -
[1] 中华医学会内分泌学分会, 中国医师协会内分泌代谢科医师分会, 中华医学会核医学分会, 等. 中国甲状腺功能亢进症和其他原因所致甲状腺毒症诊治指南[J]. 中华内分泌代谢杂志, 2022, 38(8): 700-748.Chinese Society of Endocrinology, Endocrinology and Metabolism Branch of Chinese Medical Doctor Association, Chinese Society of Nuclear Medicine, et al. Guidelines for diagnosis and management of hyperthyroidism and other causes of thyrotoxicosis in China[J]. Chinese Journal of Endocrinology and Metabolism, 2022, 38(8): 700-748. [2] 李雪彧, 刘庆阳. 甲状腺功能亢进症的治疗现状与展望[J]. 中外医学研究, 2022, 20(6): 170-172.LI X Y, LIU Q Y. Treatment status and prospect of hyperthyroidism[J]. Chinese and foreign medical research, 2022, 20(6): 170-172. [3] 周文月, 梁家浩, 李泽豪, 等. 肠道菌群与抽动障碍的因果关联: 孟德尔随机化研究[J]. 华西医学, 2024, 39(8): 1195-1203.ZHOU W Y, LIANG J H, LI Z H, et al. Causal association between intestinal flora and tic disorder: Mendelian randomization study[J]. Huaxi Medicine, 2024, 39(8): 1195-1203. [4] 杨璇, 李岩志, 马伟, 等. 基于两样本孟德尔随机化的肺功能与新型冠状病毒肺炎病死风险的因果关系[J]. 山东大学学报(医学版), 2021, 59(7): 104-111.YANG X, LI Y Z, MA W, et al. Causal relationship between lung function based on two-sample Mendelian randomization and mortality risk of COVID-19[J]. Journal of Shandong University(Medical Edition), 2021, 59(7): 104-111. [5] 郑乾, 余彩遥, 李浩云, 等. 系统性红斑狼疮与重症肌无力: 一项双向孟德尔随机化研究[J]. 中华全科医学, 2024, 22(12): 2020-2024. doi: 10.16766/j.cnki.issn.1674-4152.003788ZHENG Q, YU C Y, LI H Y, et al. Systemic lupus erythematosus and myasthenia gravis : a two-way Mendelian randomized study[J]. Chinese Journal of General Practice, 2024, 22(12): 2020-2024. doi: 10.16766/j.cnki.issn.1674-4152.003788 [6] SHIN H. XGBoost regression of the most significant photoplethysmogram features for assessing vascularaging[J]. IEEE J Biomed Health Inform, 2022, 26(7): 3354-3361. [7] 李金菊, 杨浩, 阮诺冰, 等. 肠道菌群与糖尿病神经病变的因果关系: 双样本孟德尔随机化分析[J]. 中国微生态学杂志, 2024, 36(7): 761-768.LI J J, YANG H, RUAN N B, et al. Causal relationship between gut microbiota and diabetic neuropathy: two-sample Mendelian randomization analysis[J]. Chinese Journal of Microecology, 2024, 36(7): 761-768. [8] 赵文生, 李孝林, 彭昌华, 等. 肠道菌群与骨质疏松性骨折[J]. 中国组织工程研究, 2025, 29(6): 1296-1304.ZHAO W S, LI X L, PENG C H, et al. Intestinal flora and osteoporotic fractures[J]. Chinese Tissue Engineering Research, 2025, 29(6): 1296-1304. [9] KURILSHIKOV A, MEDINA-GOMEZ C, BACIGALUPE R, et al. Large-scale association analyses identify host factors influencing human gut microbiome composition[J]. Nat Genet, 2021, 53(2): 156-165. [10] YANG Y, XIAN W, WU D D, et al. The role of obesity, type 2 diabetes, and metabolic factors in gout: a Mendelian randomization study[J]. Front Endocrinol(Lausanne), 2022, 13(8): 917056. DOI: 10.3389/fendo.2022.917056. [11] 陈继鑫, 周沁心, 郭天赐, 等. 白细胞介素-1受体拮抗剂与骨关节炎及亚型的孟德尔随机化研究[J]. 医学研究杂志, 2024, 53(4): 46-51.CHEN J X, ZHOU Q X, GUO T C, et al. Mendelian randomization study of interleukin-1 receptor antagonists and osteoarthritis and its subtypes[J]. Medical Research Journal, 2024, 53(4): 46-51. [12] 黎征鹏, 章晓云, 武瑞骐, 等. 甲状腺功能障碍与骨质疏松症的因果关系: 双向两样本孟德尔随机化研究[J]. 华西医学, 2024, 39(10): 1617-1624.LI Z P, ZHANG X Y, WU R Q, et al. Causality between thyroid dysfunction and osteoporosis: a bidirectional two-sample Mendelian randomization study[J]. Huaxi Medicine, 2024, 39(10): 1617-1624. [13] XIANG K, WANG P, XU Z W, et al. Causal effects of gut microbiome on systemic lupus erythematosus: a two-sample Mendelian randomization study[J]. Front Immunol, 2021, 12(9): 667097. DOI: 10.3389/fimmu.2021.667097. [14] XUE Y, ZHANG L Z, CHEN Y J, et al. Gut microbiota and atopic dermatitis: a two-sample Mendelian randomization study[J]. Front Med, 2023, 10(6): 1174331. DOI: 10.3389/fmed.2023.1174331. [15] 唐冬梅, 郭艳杰, 关怀, 等. 妊娠期合并甲状腺功能减退患者肠道菌群的研究[J]. 中国微生态学杂志, 2022, 34(2): 179-182.TANG D M, GUO Y J, GUAN H, et al. Study on intestinal flora in patients with hypothyroidism during pregnancy[J]. Chinese Journal of Microecology, 2022, 34(2): 179-182. [16] 冯雨. 肠道菌群介导体液免疫参与Graves病的机制研究[D]. 郑州: 郑州大学, 2022.FENG Y. The mechanism of intestinal flora-mediated humoral immunity involved in Graves'disease[D]. Zhengzhou: Zhengzhou University, 2022. [17] 李岩, 张宇霞, 单海燕. 基于两样本孟德尔随机化分析空腹血糖与冠心病的因果关系[J]. 中华全科医学, 2025, 23(1): 50-54. doi: 10.16766/j.cnki.issn.1674-4152.003833LI Y, ZHANG Y X, SHAN H Y. The causal relationship between fasting blood glucose and coronary heart disease was analyzed based on two-sample Mendelian randomization[J]. Chinese Journal of General Practice, 2025, 23(1): 50-54. doi: 10.16766/j.cnki.issn.1674-4152.003833 [18] 陈桂蓉, 宋锦彬, 潘丹旎, 等. 高校女性日常能耗与肠道菌群的相关性研究[J]. 中国微生态学杂志, 2024, 36(3): 255-263.CHEN G R, SONG J B, PAN D Y, et al. Study on the correlation between daily energy consumption and intestinal flora in college women[J]. Chinese Journal of Microecology, 2024, 36(3): 255-263. [19] 廖驰真, 曾竣. 志贺菌流行耐药现状及溯源研究进展[J]. 中国卫生检验杂志, 2023, 33(16): 2046-2048.LIAO C Z, ZENG J. Research progress on epidemic resistance status and traceability of Shigella[J]. Chinese Journal of Health Inspection, 2023, 33(16): 2046-2048. [20] 谭艳, 贾思锋, 左新河, 等. 从肠道微生物论健脾益气法治疗桥本甲状腺炎[J]. 中国微生态学杂志, 2024, 36(5): 617-621.TAN Y, JIA S F, ZUO X H, et al. Treatment of Hashimoto's thyroiditis with Jianpi Yiqi method from the perspective of intestinal microorganisms[J]. Chinese Journal of Microecology, 2024, 36(5): 617-621. -