ObjectiveOsteoarthritis (OA) is the most prevalent joint disease globally, serving as a primary cause of pain and disability. However, the pathological processes underlying OA remain incompletely understood. Several studies have noted an association between cytokines and OA, yet the causal link between them remains ambiguous. This study aims to identify cytokines potentially causally related to OA using Mendelian randomization (MR) analysis, informing early clinical diagnosis and treatment decisions.MethodsWe conducted a genome-wide association study (GWAS) on 12 OA traits involving 177,517 cases and 649,173 controls from 9 international cohorts. For discovery MR analysis, we used 103 cytokines from two European populations as instrumental variables (IVs). Concurrently, another European population OA GWAS database (36,185 cases and 135,185 controls) was used to replicate MR analysis, employing the inverse variance weighted (IVW) method as the primary analytic approach. Additional methods tested included MR Egger, Weighted median, and Weighted mode. We merged the MR findings through meta-analysis. Heterogeneity testing, level pleiotropy testing (MR Egger intercept test and MRPRESSO), and sensitivity analysis via Leave One Out (LOO) were conducted to verify result robustness. Lastly, reverse MR analysis was performed.ResultsThe meta-analysis merger revealed a correlation between CX3CL1 cycle levels and increased OA risk (OR=1.070, 95% CI: 1.040-1.110; P<0.010). We also observed associations between MCP4 (OR=0.930, 95% CI: 0.890-0.970; P<0.010) and CCL25 (OR=0.930, 95% CI: 0.890-0.970; P<0.010) with reduced OA risk. The sensitivity analysis results corroborate the robustness of these findings.ConclusionOur MR analysis indicates a potential causal relationship between CX3CL1, MCP4, CCL25, and OA risk changes. Further research is warranted to explore the influence of cytokines on OA development.