Abstract The prevalence of thyroid dysfunction diseases (TDFDs) and osteoporosis (OP) is high. Previous studies have indicated a potential association between TDFDs and OP, yet the causal direction remains unclear. This study aimed to investigate the potential causal relationship between TDFDs and the risk of developing OP and related fractures. We obtained pooled data from genome-wide association studies (GWASs) conducted on TDFDs and OP in European populations and identified single-nucleotide polymorphisms (SNPs) with genome-wide significance levels associated with exposure to TDFDs as instrumental variables. Inverse variance weighted (IVW) was employed as the primary method for Mendelian randomization (MR) analysis, supplemented by MR‒Egger, weighted median, simple mode and weighted mode methods. Sensitivity analyses were conducted to evaluate the robustness of the findings. The IVW method demonstrated an increased risk of OP in patients with TDFDs, including hyperthyroidism and hypothyroidism (TDFDs: OR = 1.11; 95% CI 1.09, 1.13; hypothyroidism: OR = 1.14; 95% CI 1.10, 1.17; hyperthyroidism: OR = 1.09; 95% CI 1.06, 1.12). These findings were supported by supplementary analysis, which revealed a positive correlation between TDFDs and the risk of OP. Multiple sensitivity analyses confirmed the absence of horizontal pleiotropy in the study, thus indicating the robustness of our results. The causal relationship between TDFDs and increased risk of OP implies the need for early bone mineral density (BMD) screening and proactive prevention and treatment strategies for individuals with TDFDs.

Investigating the causal relationship between thyroid dysfunction diseases and osteoporosis: a two-sample Mendelian randomization analysis