Skip to main content

Omeprazole and risk of osteoarthritis: insights from a mendelian randomization study in the UK Biobank

Abstract

Background

A former cohort study has raised concern regarding the unanticipated hazard of omeprazole in expediting osteoarthritis (OA) advancement. The precise nature of their causal evidence, however, remains undetermined. The present research endeavors to investigate the underlying causal link between omeprazole and OA through the application of mendelian randomization (MR) analysis.

Methods

The study incorporated the ukb-a-106 and ukb-b-14,486 datasets. The investigation of causal effects employed methodologies such as MR-Egger, Weighted median, Inverse variance weighted (IVW) with multiplicative random effects, and IVW (fixed effects). The IVW approach was predominantly considered for result interpretation. Sensitivity analysis was conducted, encompassing assessments for heterogeneity, horizontal pleiotropy, and the Leave-one-out techniques.

Results

The outcomes of the MR analysis indicated a causal relationship between omeprazole and OA, with omeprazole identified as a contributing risk factor for OA development (IVW model: OR = 1.2473, P < 0.01 in ukb-a-106; OR = 1.1288, P < 0.05 in ukb-b-14,486). The sensitivity analysis underscored the robustness and dependability of the above-mentioned analytical findings.

Conclusion

This study, employing MR, reveals that omeprazole, as an exposure factor, elevates the risk of OA. Considering the drug’s efficacy and associated adverse events, clinical practitioners should exercise caution regarding prolonged omeprazole use, particularly in populations with heightened OA risks. Further robust and high-quality research is warranted to validate our findings and guide clinical practice.

Graphical Abstract

Introduction

Osteoarthritis (OA), the most prevalent form of arthritis [1], affects around 7% of the global population [2]. With the aging population and increasing obesity rates, the incidence of OA is projected to rise [3], imposing a significant socioeconomic burden on nations worldwide [4]. To date, there has been no remedy capable of effectively halting the progression of OA. Clinical intervention primarily focuses on symptom management through approaches such as oral administration of non-steroidal anti-inflammatory drugs (NSAIDs) [5]. However, the use of NSAIDs is linked to a range of adverse effects [6, 7], such as gastrointestinal complications [7]. It is recommended to prescribe proton pump inhibitors (PPIs) along with NSAIDs as a cost-effective treatment for OA patients with moderate or high gastrointestinal risk [8,9,10,11]. As a result, PPIs are now one of the most frequently prescribed medications for OA patients.

While omeprazole, one of the most extensively prescribed PPIs, is generally considered safe and well-tolerated, concerns have arisen regarding its potential adverse effect of hypomagnesemia [12,13,14]. Magnesium, an essential ion, plays a crucial role in supporting and maintaining health [15]. Several observational studies in humans have demonstrated a negative correlation between either dietary magnesium intake or serum magnesium levels and the prevalence of radiographic knee OA [16,17,18,19]. Furthermore, a prospective cohort study found that decreased dietary magnesium intake was linked to poorer pain and function in individuals with radiographic knee OA [20]. Notably, a recent cohort study based on the general population has brought attention to the unexpected potential of omeprazole in accelerating the progression of OA [21]. Given confounding variables and the potential for reverse causation, the observational study can only establish an association between exposure and an event without being able to establish causation. At the pinnacle of the evidence hierarchy for studying therapeutic impacts on diseases stands the randomized controlled trial (RCT) [22]. Despite their effectiveness, RCTs can pose significant financial and temporal burdens, and in certain scenarios, the allocation of exposure may be deemed unethical or impractical [23,24,25].

Since the publication of a seminal article by Smith and Ebrahim at the beginning of the 21st century, mendelian randomization (MR) has garnered growing interest as a distinctive epidemiological method for inferring potential causal relationships using observational data [26]. It employs single nucleotide polymorphisms (SNPs) as instrumental variables for the exposure variable, enabling the estimation of a causal linkage between exposure and outcomes [27]. As the alleles of specific SNPs are randomly assigned at birth, genetic variations remain unaffected by potential confounders [28]. Furthermore, genetic variance is established prior to the disease’s onset, mitigating the possibility of reverse causation [29]. Evidence derived from genetically proxied exposure and disease outcome associations can be utilized to infer a causal link between exposure and outcome. Over the last decade, a plethora of MR studies have been undertaken to pinpoint causal factors linked to OA. For instance, MR investigations have offered backing for causal relationships between OA and elevated adiposity, coffee intake, bone mineral density, and sleep disruption, alongside diminished levels of serum calcium and low-density lipoprotein cholesterol [30]. These findings emphasize the potential advantages of decreasing weight and enhancing sleep quality to lower the risk of OA, while also highlighting the necessity for gaining a deeper understanding of the correlation between coffee intake and serum calcium concentrations regarding OA risk.

To address the identified research gap, this study aims to elucidate the causal relationship between omeprazole and OA through a two-sample MR analysis.

Materials and methods

Data sources

The ukb-a-106 and ukb-b-14,486 datasets were acquired through querying the Integrative Epidemiology Unit Open Genome-Wide Association Study (GWAS) database (https://gwas.mrcieu.ac.uk/) using the search term “osteoarthritis”. The ukb-a-106 dataset incorporated SNPs data from 337,159 samples with OA (10,894,596 SNPs), while the ukb-b-14,486 dataset encompassed 9,851,867 SNPs sourced from a pool of 462,933 OA samples. As for GWAS data concerning omeprazole (ukb-a-129 dataset), it encompassed 337,159 samples and 10,894,596 SNPs. In this study, OA was considered the outcome variable, while omeprazole was examined as the exposure variable. Figure 1 provides a general overview of the present study design.

Fig. 1
figure 1

The overall design of the present Mendelian randomization analysis. Abbreviations: GWAS, Genome-Wide Association Study; IEU, integrative epidemiology unit; MR, mendelian randomization; SNP, single nucleotide polymorphism

Pre-processing of data

The MR studies are rooted in three essential principles: (1) a robust and statistically significant correlation between instrumental variables (IVs) and the focal exposure; (2) the IVs’ independence from potential confounders; (3) the exclusive effect of IVs on outcomes through the designated exposure, without alternative pathways. The exposure factors were assessed and IVs were screened using the “extract_instruments” feature within the TwoSampleMR package (P < 5 × 10− 6) [31]. Then, ‘clump’ was set to TRUE to remove the IVs for linkage disequilibrium analysis (R2 = 0.001 and kb = 10,000). An R2 value of 0 indicated complete linkage equilibrium between the two SNPs, meaning the assignment of these two SNPs was completely random. IVs displaying a significant correlation with the exposure factors were selected. Outcome variables were extracted using the ‘extract_outcome_data’ function from the R package. IVs from the exposure factors were merged with screening metrics, with the proxy setting set to TRUE. Furthermore, weak IVs were omitted based on the F-statistic for each SNP.

$$F=\frac{{R}^{2}}{1-{R}^{2}}\times \frac{N-K-1}{K}$$

Within this formula, R2 signified the cumulative explained variance attributed to the selected IVs concerning the exposure factors. N represented the sample size of the GWAS, while K indicated the number of screened SNPs for the exposure. An F-statistic of ≥ 10 indicated the absence of a weak instrumental bias, thereby highlighting the robust predictive capability of the IVs for the outcome.

MR analysis

The ‘harmonize_data’ function within the TwoSampleMR package was utilized for the normalization of effect estimates. The key MR techniques included MR-Egger [32], Weighted median [33], and Inverse Variance Weighted (IVW), utilizing both multiplicative random effects and fixed effects [34]. Emphasis in the primary analysis was placed on the IVW method, which requires SNPs to fully comply with the three principles of MR research in order to obtain correct causal estimates. MR-Egger adds an additional intercept term, the main purpose of which is to determine the presence or absence of horizontal pleiotropy. The Weighted Median method assesses the presence or absence of causality using the majority of SNPs (majority of genetic variants). If there is no additional heterogeneity in the causal estimates for a given variable in IVW, then the results of the random effects and fixed effects will be the same, and therefore there is no loss of precision. In cases where the IVW method yields a noteworthy outcome (P < 0.05), it can be considered positive even if other methods’ results are not significant, given consistent directionality of β-values across methods and the absence of pleiotropy or heterogeneity. Subsequently, the findings were visually depicted using scatter plots, forest plots, and funnel plots. In the scatter plot, the IVW method remained the primary focus, with a very small intercept indicating that confounding factors had minimal influence and did not affect the reliability of the results. A positive slope of the line indicated a risk factor, while a negative slope indicated a protective factor. The forest plot was designed to assess the diagnostic efficacy of each SNP locus in predicting exposure factors for outcome diagnostics, with solid dots on the left indicating decreasing values and solid dots on the right indicating increasing values, focusing on the position of the IVW. Funnel plots were constructed to assess randomness, and if the IVs were symmetrically distributed along both sides of the IVW line, it indicated that MR conformed to Mendel’s second law of random grouping.

Sensitivity analysis

To assess the robustness of the MR analysis outcomes, a sensitivity analysis was conducted. nitially, a heterogeneity test was conducted. A Cochran’s Q test resulting in a P-value above 0.05 indicated the lack of heterogeneity. Subsequently, a test for horizontal pleiotropy was conducted, and a P-value surpassing 0.05 suggested the absence of horizontal pleiotropic effects, indicating no confounding factors in the study. Finally, a Leave-One-Out (LOO) analysis was performed by iteratively removing each SNP. Consistency in the influence of the remaining SNPs on the outcome variable validated the reliability of the MR analysis results.

Results

Omeprazole is causally associated with OA

After careful screening, a total of 27 SNPs were identified as IVs, with 14 in the ukb-a-106 dataset and 13 in the ukb-b-14,486 dataset (Table 1). Subsequently, the impact of omeprazole on OA in the UK population was evaluated through MR analysis, with omeprazole as the exposure factor and OA as the outcome. MR outcomes from the four distinct methods consistently revealed a positive causal relationship between omeprazole and OA, indicating its role as a risk factor (IVW model: OR = 1.2473, P = 0.0009 in ukb-a-106; OR = 1.1288, P = 0.0392 in ukb-b-14,486) (Table 2; Fig. 2). To clarify the connection between exposure factors and outcomes, scatter plots of SNPs were generated. These scatter plots depicted a positive linear trend for omeprazole, indicating that heightened omeprazole expression corresponded to an increased likelihood of OA development (Fig. 3).

Table 1 Detailed genome-wide association study (GWAS) data on exposure factors and outcomes
Table 2 Statistical data of mendelian randomization analysis
Fig. 2
figure 2

A forest plot illustrating the causal effect of omeprazole on the risk of osteoarthritis, generated using inverse variance weighting. Abbreviations: nsnp, number of single nucleotide polymorphisms; OR, odds ratio; CI, confidence interval, pval, p-value

Fig. 3
figure 3

Scatter plots illustrating the causal relationships between omeprazole and osteoarthritis. Abbreviations: SNP, single nucleotide polymorphism

To assess the predictive efficacy of each SNP locus regarding exposure factors and outcomes, forest plots were constructed. In these plots, solid dots on the left represented lower risk, while solid dots on the right represented higher risk. The forest plot findings consistently positioned solid dots on the right side, indicating that elevated exposure factors increased the risk of disease onset according to the IVW approach (Fig. 4). Simultaneously, assessments of instrumental variable randomness were conducted and visualized through funnel plots. These plots showed a symmetrical distribution of IVs on both sides of the IVW line, confirming that the MR analysis adhered to the principles of MR grouping (Fig. 5).

Fig. 4
figure 4

Single SNP effect combination forest plot of the causal relationships between omeprazole and osteoarthritis

Fig. 5
figure 5

Funnel plots illustrating the causal relationships between omeprazole and osteoarthritis. Abbreviations: βiv, β-value of the instrumental variable

Sensitivity analysis

To ensure the robustness of the derived conclusions, a comprehensive sensitivity analysis was undertaken. First, all Q_pval values resulting from the heterogeneity tests exceeded 0.05, signifying the absence of heterogeneity among the UK samples and emphasizing the predominant employment of the IVW (fixed effects) method in the MR analysis (Table 3). Subsequently, MR-Egger and MR-PRESSO regression tests were implemented to monitor potential horizontal pleiotropy effects stemming from genetic IVs. The horizontal pleiotropy tests consistently revealed the absence of such effects involving omeprazole and OA (P > 0.05) (Tables 4 and 5). Furthermore, a LOO analysis was executed to assess the reliability of the outcomes. LOO, which involves the iterative exclusion of one SNP while utilizing the remainder for MR, gauges the extent to which a given SNP significantly influences the results. The primary objective of LOO is to ascertain whether the line connecting the black dots remains smooth, devoid of pronounced bias points. The outcomes of the LOO analysis demonstrated no significant bias points in the figures, thus affirming the reliability of the results (Fig. 6).

Table 3 Measures of heterogeneity using Cochran’s Q test
Table 4 Examination of horizontal pleiotropy effects using MR-Egger regression tests
Table 5 Examination of horizontal pleiotropy effects using MR-PRESSO regression tests
Fig. 6
figure 6

Forest plots depicting the leave-one-out analysis of the causal relationships between omeprazole and osteoarthritis

Discussion

With the aging population and dietary changes, the prevalence and incidence of OA are escalating annually [35]. Elderly patients suffering from OA may require the simultaneous or separate long-term use of medications like NSAIDs, corticosteroids, anticoagulants, and others. Clinicians frequently prescribe omeprazole as part of the treatment regimen for these patients, with the intent of providing gastric protection [36]. Recent evidence, however, suggests a potential unexpected risk of omeprazole related to accelerated OA progression [21].

The current investigation, utilizing publicly available GWAS data resources and employing a two-sample MR analysis, unravels, for the first time, a potential causal link between omeprazole and OA.

Possible explanations

Hypomagnesemia

Magnesium ions, as a fundamental component within the human body, demonstrate a robust correlation with OA, impacting disease development and progression [16]. In a multicenter randomized controlled clinical trial involving 392 participants, researchers found that increased magnesium intake could potentially confer protective benefits against knee osteoarthritis (KOA) [37]. Conversely, a low intake of magnesium ions correlates with the gravity of OA, heightening the susceptibility to OA development and exacerbating knee joint discomfort [20, 37]. Magnesium insufficiency or diminished levels can induce an inflammatory state, thereby promoting the expression of cellular inflammatory factors such as tumor necrosis factor-α (TNF-α) and interleukin-6 (IL-6) [38, 39]. In animal experimentation, intra-articular injection of magnesium chloride to enhance magnesium ion levels acts as a mediator for hypoxia-inducible factor-1α, leading to reduced oxidative stress levels and effectively slowing the progression of KOA in rat models [40].

Omeprazole effectively halts gastric acid secretion by binding irreversibly to the hydrogen-potassium ATPase pump. Consequently, an increased pH within segments of the small intestine, induced by omeprazole, could potentially decrease the solubility of Mg2 + and hinder its absorption [14, 41]. Another presumed mechanism for the diminished absorption of magnesium by intestinal epithelial cells involves the PPIs-induced inhibition of transient receptor potential melastatin-6 (TRPM6) and TRPM7 channels [42]. A meta-analysis spanning nine observational studies, encompassing a total cohort of 109,798 patients, has unearthed a significant correlation between PPI use and the risk of hypomagnesemia [43]. Although the risk of hypomagnesemia is more pronounced in long-term PPIs users (typically exceeding one year), instances have also been documented within the first year of therapy initiation [44]. Given this potential risk, it is advisable to regularly assess serum magnesium levels in individuals on prolonged omeprazole regimens.

Vitamin B12 malabsorption

Significant causal relationships between serum vitamin B12 levels and musculoskeletal conditions, such as OA, have been observed [45]. Vitamin B12, by reducing blood homocysteine levels, has demonstrated potential protective effects against OA, particularly in weight-bearing joints and among females [46]. There are plausible mechanisms by which inhibiting gastric acid secretion may impair vitamin B12 absorption, ultimately resulting in deficiency. Parietal cells synthesize gastric acid and produce intrinsic factor, which is necessary for vitamin B12 absorption. Additionally, gastric acid is required to separate vitamin B12 from ingested dietary protein in order to facilitate its absorption. Consequently, epidemiological evidence supports a causal relationship between the use of PPIs and decreased levels of vitamin B12 [47].

Multiple studies have demonstrated the impact of omeprazole on vitamin B12 levels. For example, a trial involving healthy volunteers who underwent a Schilling test before and after two weeks of daily omeprazole therapy revealed a decrease in vitamin B12 absorption. The absorption decreased from 3.2 to 0.9% for those on a 20 mg daily dose, and from 3.4 to 0.4% for those on a 40 mg daily dose [48]. In a case-control study of individuals aged 65 years or older, the prolonged use (> 1 year) of PPIs was connected to an elevated risk of vitamin B12 deficiency [49]. Similarly, a more extensive case-control study involving community-dwelling adults highlighted an increased risk of vitamin B12 deficiency in individuals with prolonged PPI usage [50]. The utilization of omeprazole is additionally linked to decreased average serum levels of Vitamin B12 [51]. Given these findings, it is advisable to periodically assess vitamin B12 levels in individuals on long-term PPI regimens.

Intestinal flora disturbance

Gastric acid plays a crucial role in the immune system by effectively eliminating ingested pathogens within the digestive tract. In individuals with normal gastric acid physiology, the gastrointestinal microenvironment becomes unfavorable for the growth of various bacterial and fungal. However, omeprazole use can alkalize the microenvironment, potentially disrupting the gut microbiome and leading to gut microbiota disorders [47, 52].

The intricate interplay between OA and gut microbiota disturbance has emerged as a focal point in recent years [53]. According to the prevailing viewpoint, lipopolysaccharides (LPS) can traverse from the gut into the systemic circulation. They initiate mild inflammation by activating Toll-like receptors on macrophages within the joint cavity, potentially instigating OA [54]. LPS are distinctive constituents of the cell walls of gram-negative bacteria, primarily sourced from gram-negative bacteria in the distal ileum and colon. A prior study revealed a relatively elevated abundance of gram-negative rods in the feces of OA patients [55]. The causal relationship between gut microbiota and the risk of OA has been further validated by a two-sample MR analysis, which identified potentially linking gram-negative bacteria that are causally related to OA [56].

Another perspective suggests that LPS contribute to obesity by inducing metabolic syndrome and insulin resistance [57,58,59,60,61]. For weight-bearing joints like the knees and hips, excess weight inevitably imparts greater stress upon these joints. Macrophages in adipose tissue synthesize pro-inflammatory cytokines such as IL-6 and TNF-α, which might foster a low-grade systemic pro-inflammatory state, thereby triggering OA.

Clinical implications

PPIs are commonly prescribed for patients with OA, with 44.0% of patients using them [21]. Among them, omeprazole is the most commonly prescribed. Building upon previous cohort studies, this study employed MR analysis to confirm a causal relationship between omeprazole as an exposure factor and the incidence of OA. The findings have significant implications for current clinical practice. Prior research has demonstrated gender and age differences in the prescription patterns of PPIs. Females and the elderly over 60 years old have a higher frequency of taking PPIs [62]. Gender differences in PPI prescription patterns may arise because females are more likely to suffer from symptoms like gastric reflux and dyspepsia, prompting clinicians to prescribe omeprazole more frequently for female patients [63, 64]. Age differences might stem from the greater need for gastric protection through acid-suppressing drugs for patients aged 60 and above who are undergoing NSAIDs or dual antiplatelet therapy [36].

Age is a robust predictor of OA, with incidences of hand, hip, and KOA increasing notably after the age of 50 [65,66,67]. Female gender is linked to a heightened prevalence and severity of OA. Extensive meta-analyses have shown that females have a higher susceptibility to prevalent and emerging knee and hand OA, as well as incident hip OA, compared to males [68]. Furthermore, females, particularly post-menopause, exhibit pronounced severity of KOA [65, 68]. In light of these factors, we underscore the imperative of granting clinicians due consideration in their prescription practices and engaging in comprehensive evaluations of potential hazards and merits for their patients.

Strengths and limitations

Strengths

Most of the previous studies that report on the link between the use of PPIs and adverse events are observational [47]. Nevertheless, the analysis of such observational data entails a potential for erroneous associations due to factors like confounding and inadequate study design. The validity of connections between omeprazole uses and adverse events is of particular concern, given that omeprazole is frequently prescribed to individuals suffer from multiple diseases simultaneously, and these conditions themselves could be the actual cause of the adverse outcome. This is commonly known as confounding.

In this study, the application of MR effectively mitigated the interference of confounding factors such as social environment and lifestyle, as well as the impact of reverse causation on the outcomes. Sensitivity analysis further enhanced the reliability and stability of the results. The data utilized in the study were derived from a European ancestry population, which showed minimal heterogeneity. Moreover, the GWAS data originated from two independent samples, resulting in a relatively larger sample size, thus maximizing statistical power.

Limitations

First, it’s worth noting that the two-sample MR employed in this study only addresses linear associations between exposure and outcome, precluding analyses of non-linear relationships. Secondly, the study only included individuals of European descent, which limits its external validity. As a result, the effects of omeprazole on OA within other ethnic groups remain to be confirmed. Lastly, due to the absence of gender or age-stratified GWAS data, this study could not validate if the association between omeprazole and OA holds homogeneity across different genders or age groups, impeding the further elucidation of its underlying biological mechanisms.

Conclusion

The study conducted utilizing MR illustrates that omeprazole, as an exposure factor, increases the risk of OA. Clinicians should be vigilant about the prolonged use of omeprazole, especially in populations with already heightened OA risks. However, the effects and underlying mechanisms of omeprazole on OA remain poorly understood. Future research in this field should prioritize long-term, prospective studies to establish the causal relationship between omeprazole use and OA outcomes more definitively. These studies should encompass diverse population groups and investigate variations in omeprazole dosages and treatment durations. Additionally, there is a significant need for mechanistic studies to elucidate how omeprazole influences OA at the molecular and cellular levels, potentially leading to new preventive strategies or alternative treatments. Furthermore, incorporating real-world data to evaluate the impact of omeprazole across various clinical scenarios will provide a more comprehensive understanding of its effects on OA. This approach will be particularly valuable in identifying patient groups more susceptible to OA associated with omeprazole use.

Data availability

Corresponding authors can be contacted to request access to the datasets used or analyzed in this study.

Abbreviations

GWAS:

Genome-Wide Association Study

IEU:

Integrative epidemiology unit

MR:

Mendelian randomization

SNP:

Single nucleotide polymorphism

References

  1. Hunter D, Bierma-Zeinstra S. Osteoarthritis. Lancet (London England). 2019;393:1745–59.

    Article  CAS  PubMed  Google Scholar 

  2. Hunter D, March L, Chew M. Osteoarthritis in 2020 and beyond: a Lancet Commission. Lancet (London England). 2020;396:1711–2.

    Article  PubMed  Google Scholar 

  3. O’Neill T, McCabe P, McBeth J. Update on the epidemiology, risk factors and disease outcomes of osteoarthritis. Best Pract Res Clin Rheumatol. 2018;32:312–26.

    Article  PubMed  Google Scholar 

  4. Peat G, Thomas MJ. Osteoarthritis year in review 2020: epidemiology & therapy. Osteoarthritis Cartilage. 2021;29:180–9.

    Article  CAS  PubMed  Google Scholar 

  5. Zhang Y, Niu J. Editorial: shifting gears in Osteoarthritis Research toward Symptomatic Osteoarthritis. Arthritis Rheumatol. 2016;68:1797–800.

    Article  PubMed  PubMed Central  Google Scholar 

  6. McGettigan P, Henry D. Cardiovascular risk and inhibition of cyclooxygenase: a systematic review of the observational studies of selective and nonselective inhibitors of cyclooxygenase 2. JAMA. 2006;296:1633–44.

    Article  CAS  PubMed  Google Scholar 

  7. Gabriel SE, Jaakkimainen L, Bombardier C. Risk for serious gastrointestinal complications related to use of nonsteroidal anti-inflammatory drugs. A meta-analysis. Ann Intern Med. 1991;115:787–96.

    Article  CAS  PubMed  Google Scholar 

  8. McAlindon TE, Bannuru RR, Sullivan MC, Arden NK, Berenbaum F, Bierma-Zeinstra SM, Hawker GA, Henrotin Y, Hunter DJ, Kawaguchi H, et al. OARSI guidelines for the non-surgical management of knee osteoarthritis. Osteoarthritis Cartilage. 2014;22:363–88.

    Article  CAS  PubMed  Google Scholar 

  9. Hochberg MC, Altman RD, April KT, Benkhalti M, Guyatt G, McGowan J, Towheed T, Welch V, Wells G, Tugwell P. American College of Rheumatology 2012 recommendations for the use of nonpharmacologic and pharmacologic therapies in osteoarthritis of the hand, hip, and knee. Arthritis Care Res (Hoboken). 2012;64:465–74.

    Article  CAS  PubMed  Google Scholar 

  10. Latimer N, Lord J, Grant RL, O’Mahony R, Dickson J, Conaghan PG. Cost effectiveness of COX 2 selective inhibitors and traditional NSAIDs alone or in combination with a proton pump inhibitor for people with osteoarthritis. BMJ. 2009;339:b2538.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Katz JN, Smith SR, Collins JE, Solomon DH, Jordan JM, Hunter DJ, Suter LG, Yelin E, Paltiel AD, Losina E. Cost-effectiveness of nonsteroidal anti-inflammatory drugs and opioids in the treatment of knee osteoarthritis in older patients with multiple comorbidities. Osteoarthritis Cartilage. 2016;24:409–18.

    Article  CAS  PubMed  Google Scholar 

  12. Seah S, Tan YK, Teh K, Loh WJ, Tan PT, Goh LC, Malakar RD, Aw TC, Lau CS, Dhalliwal T, et al. Proton-pump inhibitor use amongst patients with severe hypomagnesemia. Front Pharmacol. 2023;14:1092476.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Park CH, Kim EH, Roh YH, Kim HY, Lee SK. The association between the use of Proton pump inhibitors and the risk of hypomagnesemia: a systematic review and meta-analysis. PLoS ONE. 2014;9:e112558.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Gommers LMM, Hoenderop JGJ, de Baaij JHF. Mechanisms of proton pump inhibitor-induced hypomagnesemia. Acta Physiol (Oxf). 2022;235:e13846.

    Article  CAS  PubMed  Google Scholar 

  15. Barbagallo M, Veronese N, Dominguez LJ. Magnesium in aging, Health and diseases. Nutrients 2021, 13.

  16. Qin B, Shi X, Samai PS, Renner JB, Jordan JM, He K. Association of dietary magnesium intake with radiographic knee osteoarthritis: results from a population-based study. Arthritis Care Res (Hoboken). 2012;64:1306–11.

    Article  CAS  PubMed  Google Scholar 

  17. Zeng C, Li H, Wei J, Yang T, Deng ZH, Yang Y, Zhang Y, Yang TB, Lei GH. Association between Dietary Magnesium Intake and Radiographic knee osteoarthritis. PLoS ONE. 2015;10:e0127666.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Zeng C, Wei J, Li H, Yang T, Zhang FJ, Pan D, Xiao YB, Yang TB, Lei GH. Relationship between serum magnesium concentration and radiographic knee osteoarthritis. J Rheumatol. 2015;42:1231–6.

    Article  CAS  PubMed  Google Scholar 

  19. Hunter DJ, Hart D, Snieder H, Bettica P, Swaminathan R, Spector TD. Evidence of altered bone turnover, vitamin D and calcium regulation with knee osteoarthritis in female twins. Rheumatology (Oxford). 2003;42:1311–6.

    Article  CAS  PubMed  Google Scholar 

  20. Shmagel A, Onizuka N, Langsetmo L, Vo T, Foley R, Ensrud K, Valen P. Low magnesium intake is associated with increased knee pain in subjects with radiographic knee osteoarthritis: data from the Osteoarthritis Initiative. Osteoarthritis Cartilage. 2018;26:651–8.

    Article  CAS  PubMed  Google Scholar 

  21. Zeng C, Neogi T, Chan AT, Wei J, Misra D, Lu N, Choi HK, Lei G, Zhang Y. Proton pump inhibitor therapy and risk of knee replacement surgery: a general population-based cohort study. Osteoarthritis Cartilage. 2022;30:559–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Bothwell LE, Greene JA, Podolsky SH, Jones DS. Assessing the Gold Standard–lessons from the history of RCTs. N Engl J Med. 2016;374:2175–81.

    Article  CAS  PubMed  Google Scholar 

  23. Larsson SC, Butterworth AS, Burgess S. Mendelian randomization for cardiovascular diseases: principles and applications. Eur Heart J. 2023;44:4913–24.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Sekula P, Del Greco MF, Pattaro C, Köttgen A. Mendelian randomization as an Approach to assess causality using Observational Data. J Am Soc Nephrol. 2016;27:3253–65.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23:R89–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Smith GD, Ebrahim S. Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32:1–22.

    Article  PubMed  Google Scholar 

  27. Sekula P, Pattaro MFDG, Köttgen C. Mendelian randomization as an Approach to assess causality using Observational Data. J Am Soc Nephrol. 2016;27:3253–65.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Burgess S, Butterworth A, Malarstig A, Thompson SG. Use of mendelian randomisation to assess potential benefit of clinical intervention. BMJ. 2012;345:e7325.

    Article  PubMed  Google Scholar 

  29. Davies NM, Holmes MV, Davey Smith G. Reading mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. 2018;362:k601.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Alhassan E, Nguyen K, Hochberg MC, Mitchell BD. Causal factors for Osteoarthritis: a scoping review of mendelian randomization studies. Arthritis Care Res (Hoboken). 2024;76:366–75.

    Article  PubMed  Google Scholar 

  31. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, Laurin C, Burgess S, Bowden J, Langdon R et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife 2018, 7.

  32. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–25.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some Invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40:304–14.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG. Using published data in mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30:543–52.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Safiri S, Kolahi AA, Smith E, Hill C, Bettampadi D, Mansournia MA, Hoy D, Ashrafi-Asgarabad A, Sepidarkish M, Almasi-Hashiani A, et al. Global, regional and national burden of osteoarthritis 1990–2017: a systematic analysis of the global burden of Disease Study 2017. Ann Rheum Dis. 2020;79:819–28.

    Article  PubMed  Google Scholar 

  36. Benmassaoud A, McDonald EG, Lee TC. Potential harms of proton pump inhibitor therapy: rare adverse effects of commonly used drugs. CMAJ. 2016;188:657–62.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Zhang Y, Chen T, Luo P, Li S, Zhu J, Xue S, Cao P, Zhu Z, Li J, Wang X et al. Associations of Dietary Macroelements with knee joint structures, symptoms, Quality of Life, and Comorbid conditions in people with symptomatic knee osteoarthritis. Nutrients 2022, 14.

  38. Weglicki WB, Phillips TM, Freedman AM, Cassidy MM, Dickens BF. Magnesium-deficiency elevates circulating levels of inflammatory cytokines and endothelin. Mol Cell Biochem. 1992;110:169–73.

    Article  CAS  PubMed  Google Scholar 

  39. Song Y, Li TY, van Dam RM, Manson JE, Hu FB. Magnesium intake and plasma concentrations of markers of systemic inflammation and endothelial dysfunction in women. Am J Clin Nutr. 2007;85:1068–74.

    Article  CAS  PubMed  Google Scholar 

  40. Yao H, Xu JK, Zheng NY, Wang JL, Mok SW, Lee YW, Shi L, Wang JY, Yue J, Yung SH, et al. Intra-articular injection of magnesium chloride attenuates osteoarthritis progression in rats. Osteoarthritis Cartilage. 2019;27:1811–21.

    Article  CAS  PubMed  Google Scholar 

  41. Hess MW, Hoenderop JG, Bindels RJ, Drenth JP. Systematic review: hypomagnesaemia induced by Proton pump inhibition. Aliment Pharmacol Ther. 2012;36:405–13.

    Article  CAS  PubMed  Google Scholar 

  42. Perazella MA. Proton pump inhibitors and hypomagnesemia: a rare but serious complication. Kidney Int. 2013;83:553–6.

    Article  CAS  PubMed  Google Scholar 

  43. Cheungpasitporn W, Thongprayoon C, Kittanamongkolchai W, Srivali N, Edmonds PJ, Ungprasert P, O’Corragain OA, Korpaisarn S, Erickson SB. Proton pump inhibitors linked to hypomagnesemia: a systematic review and meta-analysis of observational studies. Ren Fail. 2015;37:1237–41.

    Article  CAS  PubMed  Google Scholar 

  44. Hoorn EJ, van der Hoek J, de Man RA, Kuipers EJ, Bolwerk C, Zietse R. A case series of proton pump inhibitor-induced hypomagnesemia. Am J Kidney Dis. 2010;56:112–6.

    Article  PubMed  Google Scholar 

  45. Fu L, Wang Y, Hu YQ. Causal effects of B vitamins and homocysteine on obesity and musculoskeletal diseases: a mendelian randomization study. Front Nutr. 2022;9:1048122.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Hong H, Chen L, Zhong Y, Yang Z, Li W, Song C, Leng H. Associations of Homocysteine, Folate, and vitamin B12 with osteoarthritis: a mendelian randomization study. Nutrients 2023, 15.

  47. Elias E, Targownik LE. The Clinician’s guide to Proton pump inhibitor related adverse events. Drugs. 2019;79:715–31.

    Article  PubMed  Google Scholar 

  48. Marcuard SP, Albernaz L, Khazanie PG. Omeprazole therapy causes malabsorption of cyanocobalamin (vitamin B12). Ann Intern Med. 1994;120:211–5.

    Article  CAS  PubMed  Google Scholar 

  49. Valuck RJ, Ruscin JM. A case-control study on adverse effects: H2 blocker or proton pump inhibitor use and risk of vitamin B12 deficiency in older adults. J Clin Epidemiol. 2004;57:422–8.

    Article  PubMed  Google Scholar 

  50. Lam JR, Schneider JL, Zhao W, Corley DA. Proton pump inhibitor and histamine 2 receptor antagonist use and vitamin B12 deficiency. JAMA. 2013;310:2435–42.

    Article  CAS  PubMed  Google Scholar 

  51. Hoxha. Q-BH, Sadiku R, Bajraktari. S, MS IH, Thaçi. K, Thaçi. S, Bahtiri E. Proton Pump Inhibitors Intake and Iron and vitamin B12 status: a prospective comparative study with a follow up of 12 months. Open Access Maced J Med Sci. 2018;6:442–6.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Imhann F, Bonder MJ, Vich Vila A, Fu J, Mujagic Z, Vork L, Tigchelaar EF, Jankipersadsing SA, Cenit MC, Harmsen HJ, et al. Proton pump inhibitors affect the gut microbiome. Gut. 2016;65:740–8.

    Article  CAS  PubMed  Google Scholar 

  53. Ramires LC, Santos GS, Ramires RP, da Fonseca LF, Jeyaraman M, Muthu S, Lana AV, Azzini G, Smith CS, Lana JF. The Association between Gut Microbiota and Osteoarthritis: does the disease begin in the gut? Int J Mol Sci 2022, 23.

  54. Huang Z, Kraus VB. Does lipopolysaccharide-mediated inflammation have a role in OA? Nat Rev Rheumatol. 2016;12:123–9.

    Article  CAS  PubMed  Google Scholar 

  55. Liu S, Li G, Zhu Y, Xu C, Yang Q, Xiong A, Weng J, Yu F, Zeng H. Analysis of gut microbiome composition, function, and phenotype in patients with osteoarthritis. Front Microbiol. 2022;13:980591.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Yu XH, Yang YQ, Cao RR, Bo L, Lei SF. The causal role of gut microbiota in development of osteoarthritis. Osteoarthritis Cartilage. 2021;29:1741–50.

    Article  PubMed  Google Scholar 

  57. Ali S, Malloci M, Safiedeen Z, Soleti R, Vergori L, Vidal-Gómez X, Besnard C, Dubois S, Le Lay S, Boursier J, et al. LPS-enriched small extracellular vesicles from metabolic syndrome patients trigger endothelial dysfunction by activation of TLR4. Metabolism. 2021;118:154727.

    Article  CAS  PubMed  Google Scholar 

  58. Croci S, D’Apolito LI, Gasperi V, Catani MV, Savini I. Dietary strategies for management of metabolic syndrome: role of gut microbiota metabolites. Nutrients 2021, 13.

  59. Cani PD, Neyrinck AM, Fava F, Knauf C, Burcelin RG, Tuohy KM, Gibson GR, Delzenne NM. Selective increases of bifidobacteria in gut microflora improve high-fat-diet-induced diabetes in mice through a mechanism associated with endotoxaemia. Diabetologia. 2007;50:2374–83.

    Article  CAS  PubMed  Google Scholar 

  60. Kootte RS, Vrieze A, Holleman F, Dallinga-Thie GM, Zoetendal EG, de Vos WM, Groen AK, Hoekstra JB, Stroes ES, Nieuwdorp M. The therapeutic potential of manipulating gut microbiota in obesity and type 2 diabetes mellitus. Diabetes Obes Metab. 2012;14:112–20.

    Article  CAS  PubMed  Google Scholar 

  61. Saad MJ, Santos A, Prada PO. Linking gut microbiota and inflammation to obesity and insulin resistance. Physiol (Bethesda). 2016;31:283–93.

    CAS  Google Scholar 

  62. Abrahami D, McDonald EG, Schnitzer M, Azoulay L. Trends in acid suppressant drug prescriptions in primary care in the UK: a population-based cross-sectional study. BMJ Open. 2020;10:e041529.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Kim YS, Kim N, Kim GH. Sex and gender differences in Gastroesophageal Reflux Disease. J Neurogastroenterol Motil. 2016;22:575–88.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Piessevaux H, De Winter B, Louis E, Muls V, De Looze D, Pelckmans P, Deltenre M, Urbain D, Tack J. Dyspeptic symptoms in the general population: a factor and cluster analysis of symptom groupings. Neurogastroenterol Motil. 2009;21:378–88.

    Article  CAS  PubMed  Google Scholar 

  65. Johnson VL, Hunter DJ. The epidemiology of osteoarthritis. Best Pract Res Clin Rheumatol. 2014;28:5–15.

    Article  PubMed  Google Scholar 

  66. Allen KD, Golightly YM. State of the evidence. Curr Opin Rheumatol. 2015;27:276–83.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Neogi T, Zhang Y. Epidemiology of osteoarthritis. Rheum Dis Clin North Am. 2013;39:1–19.

    Article  PubMed  Google Scholar 

  68. Srikanth VK, Fryer JL, Zhai G, Winzenberg TM, Hosmer D, Jones G. A meta-analysis of sex differences prevalence, incidence and severity of osteoarthritis. Osteoarthritis Cartilage. 2005;13:769–81.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to express our gratitude to Home for Researchers (https://www.home-for-researchers.com) for providing language assistance and to all organizations that provided funding for this work. Figdraw (https://www.figdraw.com/static/index.html#/) was utilized for creating the graphic abstract.

Funding

This study received funding from various grants: the National Natural Science Foundation of China (Nos. 82172432, 82102568, 82001319); the National & Local Joint Engineering Research Centre of Orthopaedic Biomaterials (No. XMHT20190204007); the Guangdong Basic and Applied Basic Research Foundation (Nos. 2021A1515012586, 2021A1515220054, 2022B1515120046, 2022A1515220038, 2022A1515220111, 2021A1515220037, 2022A1515220165); the Shenzhen High-Level Hospital Construction Fund; the Shenzhen Key Laboratory of Orthopaedic Diseases and Biomaterials Research (No. ZDSYS20220606100602005); the Shenzhen Key Medical Discipline Construction Fund (No. SZXK023); the Sanming Project of Medicine in Shenzhen (No. SZSM202211038); the Research and Development Projects of Shenzhen (Nos. JCYJ20210324110214040, JCYJ20220531094214032, JCYJ20190809152409606, JCYJ20210324105806016, JCYJ20220531094406015); the Shenzhen Sustainable Development Project (No. KCXFZ20201221173411031); and the Scientific Research Foundation of Peking University Shenzhen Hospital (No. KYQD2021099).

Author information

Authors and Affiliations

Authors

Contributions

SYC: Conceptualisation, data analysis and manuscript writing. YHW: Data curation, investigation and manuscript writing. GQL, HLW and YHY: Data curation, reviewing and editing. PL, JJL, and QCW: Methodology, software and supervision. AX, FY and HZ: Project administration and funding acquisition.

Corresponding authors

Correspondence to Fei Yu, Ao Xiong or Hui Zeng.

Ethics declarations

Ethics approval and consent to participate

Informed consent and ethical approval are not required for Mendelian randomization studies.

Consent for publication

Not applicable.

Competing Interests

The authors declare no conflicts of interest related to commercial or financial relationships in the conduct of this research.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cao, S., Wei, Y., Yue, Y. et al. Omeprazole and risk of osteoarthritis: insights from a mendelian randomization study in the UK Biobank. J Transl Med 22, 504 (2024). https://doi.org/10.1186/s12967-024-05255-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-024-05255-y

Keywords