Objective To systematically assess the potential causal relationship between circulating inflammatory cytokines and osteoporosis (OP) using a bidirectional Mendelian randomization (MR) approach. Methods Genetic association data for 91 inflammatory cytokines and OP were obtained from genome-wide association studies (GWAS). The causal effects were estimated using several two-sample MR methods, including inverse-variance weighted (IVW), MR-Egger regression, weighted median, and mode-based methods. Heterogeneity and pleiotropy were assessed using Cochran's Q test, MR-Egger intercept, and MR-PRESSO. Results The IVW method identified five cytokines significantly associated with OP risk: CXCL1 (OR=1.100, P=0.020), CXCL11 (OR=1.150, P=0.005), IL-18 (OR=1.087, P=0.012), TNF (OR=1.108, P=0.014) as risk factors, and ARTN (OR=0.895, P=0.015) as a protective factor. Sensitivity analyses supported these causal associations. Interestingly, OP was found to potentially reduce the circulating level of FGF21 (β=-0.076, P=0.003). Conclusion This study provides genetic evidence that CXCL1, CXCL11, IL-18, and TNF may causally increase the risk of OP, while ARTN may be a novel protective factor. Moreover, OP itself may aggravate bone loss by suppressing FGF21 secretion. These findings offer new insights into the pathogenesis of OP and potential intervention targets for its prevention and treatment. |