Causal Association of Coffee Consumption and Total, Knee, Hip and Self-Reported Osteoarthritis: A Mendelian Randomization Study

Causal factors for osteoarthritis risk revealed by mendelian randomization analysis

Osteoarthritis (OA), a prevalent chronic disease among the elderly, presents a complex pathogenesis and currently lacks effective treatment. Traditional observational studies are time-consuming, labor-intensive, susceptible to confounding factors, and cannot establish causal relationships. Mendelian randomization (MR) analysis, leveraging genetic variation to assess causal associations between exposures and outcomes, offers a cost-effective and efficient alternative. Over the past decade, large-scale genome-wide association studies have identified numerous genetic variants linked to OA risk factors, facilitating MR study design. In this review, we systematically identified 52 MR studies meeting specific criteria and evaluated their quality, exploring the impact of lifestyle, nutrition, comorbidities, circulating metabolites, plasma proteins, and other health factors on OA risk. We discuss the results and potential mechanisms of MR findings, addressing conflicting evidence based on existing literature and our prior research. With the ongoing expansion of genome-wide association data, we anticipate MR's role in future OA studies to broaden, particularly in drug development research using targeted MR approaches. We thus aim for this paper to offer valuable insights for researchers and clinicians in related fields.

A causal relationship between hypothyroidism and rheumatoid arthritis, but not hyperthyroidism: evidence from the mendelian randomization study

OBJECTIVE

To investigate the genetic level causal association among hyperthyroidism, hypothyroidism, and rheumatoid arthritis (RA).

METHODS

We utilized the genome-wide association studies (GWAS) summary data for exposure (hyperthyroidism and hypothyroidism) and outcome (RA) from the IEU OpenGWAS database. We used two different sets of data (test cohort and validation cohort) for causal assessment of exposure and outcome. To establish a causal relationship between these conditions, we conducted a two-sample Mendelian randomization (MR) analysis. Subsequently, we evaluated the MR analysis results for heterogeneity, horizontal pleiotropy, and outliers, aiming to assess the validity and reliability of the findings. Moreover, we conducted additional analyses to examine the robustness of the MR results, including a "Leave one out" analysis and the MR robust adjusted profile score (MR-RAPS) method, ensuring the robustness and adherence to normal distribution assumptions.

RESULTS

The findings from the test cohort indicated that hyperthyroidism did not exhibit a genetic causal association with RA (P = 0.702, odds ratio [OR] 95% confidence interval [CI] = 1.021 [0.918-1.135]). Conversely, hypothyroidism displayed a positive genetic causal relationship with RA (P < 0.001, OR 95% CI = 1.239 [1.140-1.347]). The analysis results of the validation cohort are consistent with those of the test cohort. Notably, our MR analysis results demonstrated no evidence of heterogeneity, horizontal pleiotropy, or outliers. Furthermore, our MR analysis results remained unaffected by any single nucleotide polymorphism (SNP) and exhibited a normal distribution.

CONCLUSION

The results of this study showed that hypothyroidism was positively correlated with RA, while hyperthyroidism was not causally correlated with RA. Hypothyroidism may as a risk factor of RA should be paid attention to in clinical work. Future studies are needed to further confirm this finding.

The role of dietary preferences in osteoarthritis: a Mendelian randomization study using genome-wide association analysis data from the UK Biobank

Background

To understand the impact of individual preferences for specific dietary items on OA, and to help inform the development of effective and targeted OA prevention and management strategies, we performed a Mendelian randomization analysis between dietary preferences and osteoarthritis.

Methods

This study utilized genetic data from the UK Biobank to investigate the association between OA and 21 different common dietary items. Instrumental variables representing European populations were carefully selected based on their genetic significance and linkage disequilibrium. In cases where a dietary item had few relevant genetic markers, a more lenient selection threshold was applied. To prevent bias, the analysis excluded single nucleotide polymorphisms (SNPs) associated with factors such as body mass index (BMI) and cholesterol. Using inverse-variance weighting (IVW) and Mendelian randomization, significant associations were detected between certain dietary items and OA.

Results

Using Mendelian randomization to examine the relationship between 21 different dietary items and OA, significant associations were found for coffee, peas, watercress, and cheese, where the first two had a promoting effect and the last two an inhibiting effect on OA. Due to heterogeneity in the test results for cheese, a random IVW representation was used. The results of sensitivity analysis showed no significant heterogeneity or horizontal pleiotropy in the selected SNPS, demonstrating the reliability of Mendelian randomization analysis.

Conclusion

This study identified coffee, peas, watercress, and cheese as food items that may have significant dietary effects on osteoarthritis. This information may be useful to consider in the development of OA management strategies.

Osteoarthritis and hypertension: observational and Mendelian randomization analyses

BACKGROUND

The association between osteoarthritis (OA) and hypertension is a subject of ongoing debate in observational research, and the underlying causal relationship between them remains elusive.

METHODS

This study retrospectively included 24,871 participants in the National Health and Nutrition Examination Survey (NHANES) from 2013 to 2020. Weighted logistic regression was performed to investigate the connection between OA and hypertension. Additionally, Mendelian randomization (MR) analysis was conducted to explore the potential causal relationship between OA and hypertension.

RESULTS

In the NHANES data, after adjusting for multiple confounding factors, there was no significant relationship between OA and hypertension (OR 1.30, 95% CI, 0.97-1.73, P = 0.089). However, among males, OA appeared to be associated with a higher risk of hypertension (OR 2.25, 95% CI, 1.17-4.32, P = 0.019). Furthermore, MR results indicate no relationship between multiple OA phenotypes and hypertension: knee OA (IVW, OR 1.024, 95% CI: 0.931-1.126, P = 0.626), hip OA (IVW, OR 0.990, 95% CI: 0.941-1.042, P = 0.704), knee or hip OA (IVW, OR 1.005, 95% CI: 0.915-1.105, P = 0.911), and OA from UK Biobank (IVW, OR 0.796, 95% CI: 0.233-2.714, P = 0.715). Importantly, these findings remained consistent across different genders and in reverse MR.

CONCLUSIONS

Our study found that OA patients had a higher risk of hypertension only among males in the observational study. However, MR analysis did not uncover any causal relationship between OA and hypertension.

A protective role of genetically predicted sex hormone-binding globulin on stroke

Introduction

The role of sex hormone-binding globulin (SHBG) on stroke has been investigated in several observational studies. To provide the causal estimates of SHBG on stroke and its subtypes, bi-directional and multivariable Mendelian randomization (MR) analyses are performed.

Methods

The genetic instruments of SHBG were obtained from the UK Biobank. Outcome datasets for stroke and its subtypes were taken from the MEGASTROKE Consortium. The main analysis used in this study is the inverse variance weighting, complemented by other sensitivity approaches to verify the conformity of findings.

Results

We found that the risk of stroke grew by 13% (odd ratio [OR] = 0.87, 95% confidence interval [CI] = 0.79-0.95, P = 0.0041) and the risk of ischemic stroke grew by 15% (OR = 0.85, 95%CI = 0.77-0.95, P = 0.0038) caused by genetically predicted SHBG. The causal association remains robust in the reverse MR and multivariable MR analyses for stroke (reverse MR: all P > 0.01 for the IVW method; MVMR: OR = 0.72, 95%CI = 0.59-0.87, P = 0.0011) and ischemic stroke (reverse MR: all P > 0.01 for IVW; MVMR: OR = 0.70, 95%CI = 0.56-0.86, P = 0.0007).

Conclusion

Our MR study provides novel evidence that SHBG has an inverse association with stroke and ischemic stroke, exerting protective effects on stroke.

Causal Factors for Osteoarthritis: A Scoping Review of Mendelian Randomization Studies

OBJECTIVE

Mendelian randomization (MR) has increasingly been utilized as a tool for establishing causal relations between modifiable exposures and osteoarthritis (OA). The goal of this review was to summarize available MR studies of OA that evaluate the causal role of modifiable risk factors on OA.

METHODS

This review was performed following the Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) Extension for Scoping Reviews model. We performed a literature search for relevant studies published before December 2021 across multiple databases using the search terms "osteoarthritis" and ("Mendelian randomization" or "polygenic risk score"). We reported the MR estimates of causal associations between exposures and OA and then assessed methodologic quality of abstracted studies according to their efforts to validate the three key MR assumptions.

RESULTS

Our search identified 45 studies reporting on 141 exposure-association analyses. All studies performed a formal instrumental variable analysis to estimate the causal effect of exposure on OA. Causal associations (P < 0.05) were reported in 60 of these analyses representing 36 unique publications, and MR-Egger sensitivity analyses were performed in 45 of these analyses. MR studies provided support for causal associations of OA with increased levels of adiposity, coffee consumption, bone mineral density, and sleep disturbance, and decreased levels of serum calcium and low-density lipoprotein cholesterol.

CONCLUSION

These results highlight the potential benefits of weight reduction and improvement of sleep quality to reduce the risk of OA and call for a better understanding of the relations of coffee consumption and serum calcium to OA risk.

Genetic investigation into the broad health implications of caffeine: evidence from phenome-wide, proteome-wide and metabolome-wide Mendelian randomization

BACKGROUND

Caffeine is one of the most utilized drugs in the world, yet its clinical effects are not fully understood. Circulating caffeine levels are influenced by the interplay between consumption behaviour and metabolism. This study aimed to investigate the effects of circulating caffeine levels by considering genetically predicted variation in caffeine metabolism.

METHODS

Leveraging genetic variants related to caffeine metabolism that affect its circulating levels, we investigated the clinical effects of plasma caffeine in a phenome-wide association study (PheWAS). We validated novel findings using a two-sample Mendelian randomization framework and explored the potential mechanisms underlying these effects in proteome-wide and metabolome-wide Mendelian randomization.

RESULTS

Higher levels of genetically predicted circulating caffeine among caffeine consumers were associated with a lower risk of obesity (odds ratio (OR) per standard deviation increase in caffeine = 0.97, 95% confidence interval (CI) CI: 0.95-0.98, p = 2.47 × 10-4), osteoarthrosis (OR = 0.97, 95% CI: 0.96-0.98, P=1.10 × 10-8) and osteoarthritis (OR: 0.97, 95% CI: 0.96 to 0.98, P = 1.09 × 10-6). Approximately one third of the protective effect of plasma caffeine on osteoarthritis risk was estimated to be mediated through lower bodyweight. Proteomic and metabolomic perturbations indicated lower chronic inflammation, improved lipid profiles, and altered protein and glycogen metabolism as potential biological mechanisms underlying these effects.

CONCLUSIONS

We report novel evidence suggesting that long-term increases in circulating caffeine may reduce bodyweight and the risk of osteoarthrosis and osteoarthritis. We confirm prior genetic evidence of a protective effect of plasma caffeine on risk of overweight and obesity. Further clinical study is warranted to understand the translational relevance of these findings before clinical practice or lifestyle interventions related to caffeine consumption are introduced.

Coffee and the risk of osteoarthritis: a two-sample, two-step multivariable Mendelian randomization study

Purpose: To investigate the potential causal relationship between coffee consumption and osteoarthritis (OA), and to disentangle whether body mass index (BMI) and Bone mineral density (BMD) mediate this relationship. Methods: We performed two-sample and two-step Mendelian randomization (MR) analyses utilizing publicly available genome-wide association studies (GWAS) summary statistics to estimate the association between coffee intake and OA risk (including knee OA, hip OA, knee or hip OA, and total OA), as well as the possible mediating effects of BMI and BMD. In addition, data of different coffee types (decaffeinated coffee, instant coffee, ground coffee-including espresso, filter, etc., and other coffee types) were used to explore the effect of coffee type on the risk of OA. Results: In two-sample MR, coffee intake increased the risk of OA in various sites, with the most significant impact observed in knee osteoarthritis (KOA) (odds ratio [OR] 2.03, 95% confidence interval [CI] 1.57-2.61, p < 0.001). The effect on self-reported OA was minimal (OR 1.03, 95% CI 1.01-1.05, p = 0.006). Further analysis of different types of coffee revealed that only decaffeinated coffee was causally associated with both KOA (OR 4.40, 95% CI 1.71-11.33, p = 0.002) and self-reported OA (OR 1.13, 95% CI 1.02-1.26, p = 0.022). In two-step MR, BMI explained over half of the coffee intake-all OA risk association, while BMD accounted for less than 5% of the mediation effect. Conclusion: Our study suggests that coffee intake increase the risk of OA, with BMI playing a significant mediating role. Decaffeinated coffee appears to have the greatest impact on OA risk compared to other types of coffee. Therefore, managing BMI and selecting appropriate types of coffee should be included in the health management of individuals who frequently consume coffee.

Causal relationship between diet and knee osteoarthritis: A Mendelian randomization analysis

BACKGROUND

Knee osteoarthritis (KOA) is a common disabling joint disease that affects millions of people worldwide. Diet may play a role in the etiology and progression of KOA, but evidence for a causal relationship is limited. We aimed to investigate the causal impact of dietary intake on KOA risk using Mendelian randomization (MR).

METHODS

We used summary-level data from genome-wide association studies (GWAS) including dietary intake (n = 335, 394-462, 342), and KOA (n = 403, 124). We selected 6-77 genetic variants as instrumental variables for 18 dietary factors, including processed meat, poultry, beef, oily fish, non-oily fish, pork, lamb, frequency of alcohol intake, alcoholic beverages, tea, coffee, dried fruit, cereals, cheese, bread, cooked vegetables, salad/raw vegetables, and fresh fruit. We performed univariate and multivariate MR analyses to estimate the causal effect of each dietary factor on KOA risk. We also performed some sensitivity analyses to assess the validity of the MR hypothesis.

RESULTS

We found that higher coffee intake was associated with increased KOA risk, whereas higher intake of dried fruits, grains, cheese, and oily fish was associated with reduced KOA risk. After multivariate adjustment, we found that coffee and oily fish intake may affect KOA through obesity, body mass index (BMI), diabetes, hypertension, and prolonged standing. Sensitivity analyses did not reveal any evidence of pleiotropy.

CONCLUSIONS

Our study provides new causal evidence that dietary intake may influence KOA risk. Specifically, we suggest that increased intake of dried fruits, grains, cheese, and oily fish and decreased coffee intake may be beneficial in preventing and mitigating KOA. further studies are needed to elucidate the underlying mechanisms and to confirm our findings in different populations.

The causal relationship between osteoarthritis and bladder cancer: A Mendelian randomization study

OBJECTIVE

The causal association between osteoarthritis (OA) and bladder cancer remains unclear. This Mendelian randomization (MR) study was carried out to assess the potential causal effects of any OA, knee OA and hip OA, and bladder cancer.

METHOD

Genome-wide association study (GWAS) summary data for OA and bladder cancer were obtained in GWAS CATALOG, UK Biobank, and FinnGen Consortium. Inverse-variance weighted (IVW) approach was primarily conducted to evaluate the causal relationships between OA and bladder cancer, and MR-Egger intercept and Cochran's Q test were mainly used to estimate heterogeneity and pleiotropy. MR-PRESSO was used to test the presence of horizontal outliers. Leave-one-out analysis was utilized to ensure the reliability of the results.

RESULTS

A higher genetic predisposition to any OA has a causal association with bladder cancer risk, while neither knee OA nor hip OA is causally linked to bladder cancer. MR-Egger intercept analysis exhibited that any OA and knee OA had no pleiotropic effect on the risk of bladder cancer, and Cochran's Q test showed that any OA, knee OA and hip OA had no heterogeneity on bladder cancer risk. Neither MR PRESSO analysis nor leave-one-out analysis revealed any outlier SNPs.

CONCLUSIONS

This MR study exhibited a positive cause-and-effect relationship between any type of OA and bladder cancer risk, but not between site-specific OA, knee OA and hip OA, and bladder cancer. Attention should be paid to the screening and prevention of bladder cancer in OA patients at any site.

The Causal Association between Alcohol, Smoking, Coffee Consumption, and the Risk of Arthritis: A Meta-Analysis of Mendelian Randomization Studies

Objective: To evaluate the genetic causality between alcohol intake, smoking, coffee consumption, and arthritis. Methods: Mendelian randomization (MR) studies with alcohol, smoking, and coffee consumption behaviors as exposures, and osteoarthritis (OA) and rheumatoid arthritis (RA) as outcomes were retrieved from up to July 2023. Two researchers with relevant professional backgrounds independently assessed the quality and extracted data from the included studies. Meanwhile, we applied MR analyses of four lifestyle exposures and five arthritis outcomes (two for OA and three for RA) with gene-wide association study (GWAS) data that were different from the included studies, and the results were also included in the meta-analysis. Statistical analyses were performed using Stata 16.0 and R software version 4.3.1. Results: A total of 84 studies were assessed. Of these, 11 were selected for meta-analysis. As a whole, the included studies were considered to be at a low risk of bias and were of high quality. Results of the meta-analysis showed no significant genetic causality between alcohol intake and arthritis (odds ratio (OR): 1.02 (0.94-1.11)). Smoking and arthritis had a positive genetic causal association (OR: 1.44 (1.27-1.64)) with both OA (1.44 (1.22-1.71)) and RA (1.37 (1.26-1.50)). Coffee consumption and arthritis also had a positive genetic causal association (OR: 1.02 (1.01-1.03)). Results from the subgroup analysis showed a positive genetic causality between coffee consumption and both OA (OR: 1.02 (1.00-1.03)) and RA (OR: 1.56 (1.19-2.05)). Conclusion: There is positive genetic causality between smoking and coffee consumption and arthritis (OA and RA), while there is insufficient evidence for genetic causality between alcohol intake and arthritis.

Causal association between cardiovascular diseases and erectile dysfunction, a Mendelian randomization study

Background

Cardiovascular diseases (CVD), including coronary heart disease (CHD), heart failure, ischemic heart disease (IHD), and atrial fibrillation, are prevalent in the aged. However, the influence of CVD on ED is less investigated. This study was performed to clarify the causal association between CVD and ED.

Materials and methods

Genome-wide association studies (GWAS) datasets targeting CHD, heart failure, IHD, and atrial fibrillation were downloaded to retrieve single nucleotide polymorphisms (SNPs). Further, single-variable Mendelian randomization and multivariable Mendelian randomization (MVMR) were adopted to explore the causal association between CVD and ED.

Results

Genetically predicted CHD and heart failure were found to increase the risks of ED (OR = 1.09, P < 0.05 and OR = 1.36, P < 0.05, respectively). However, no causal association was disclosed among IHD, atrial fibrillation and ED (all P > 0.05). These findings remained consistent in sensitivity analyses. After controlling for body mass index, alcohol, low density lipoprotein, smoking and total cholesterol levels, the results of MVMR support the causal role of CHD on ED (P < 0.05). Similarly, the direct causal effect estimates of heart failure on ED were significant in MVMR analyses (P < 0.05).

Conclusion

Using genetic data, this study revealed that genetically predicted CHD and heart failure may predict better ED compared with atrial fibrillation and IHD. The results should be interpreted with caution and the insignificant causal inference of IHD still needs further verification in future studies.

Genetically predicted insomnia causally increases the risk of erectile dysfunction

Sleep has attracted extensive attention due to its significance in health. However, its association with erectile dysfunction (ED) is insufficiently investigated. To investigate the potential causal links between sleep traits (insomnia, sleep duration, and chronotype) and ED, this study was performed. The single-nucleotide polymorphisms (SNPs) associated with insomnia, sleep duration, and chronotype were retrieved from previous genome-wide association studies (GWAS). A conventional two-sample Mendelian randomization (MR) was used to estimate the causal links between sleep traits and ED. The summary statistics of ED were from individuals of European ancestry (6175 cases vs 217 630 controls). As shown by the random effect inverse-variance-weighting (IVW) estimator, genetically predicted insomnia was causally associated with a 1.15-fold risk of ED (95% confidence interval: 1.07-1.23, P < 0.001). Sleep duration and morningness were not causally associated with ED, as indicated by the IVW (all P > 0.05). These findings were consistent with the results of sensitivity analyses. Based on genetic data, this study provides causal evidence that genetically predicted insomnia increases the risk of ED, whereas sleep duration and chronotype do not.

Causal Association Between Tea Consumption and Bone Health: A Mendelian Randomization Study

Background

Much observational research reported that tea consumption decreases the risk of osteoarthritis (OA), rheumatoid arthritis (RA), and osteoporosis (OP) which are the three major bone disorders. However, the observed correlation is inconclusive. To determine the causal relationship between genetically predicted tea intake and OA, RA, and OP, we performed a two-sample Mendelian randomization (MR) study based on large samples.

Methods

The European population's genome-wide association meta-analysis (GWAS) dataset identified SNPs associated with tea consumption was obtained from Neale Lab's analysis of UK Biobank data that comprised 349,376 participants of European ancestry. We extracted genetic data for knee OA (17,885 controls and 4,462 cases), hip OA (50,898 controls and 12,625 cases), and RA (43,923 controls and 14,361 cases) from the UK Biobank and OP cases (93083 controls and 1,175 cases) from FinnGen Data Freeze 2. A MR study was conducted to examine the effect of selected single nucleotide polymorphisms (SNPs) and OA, RA, and OP risk. Several sensitivity analyses were performed with weighted median and inverse-variance weighted methods for estimating the causal effects.

Results

In this MR study, the genetically predicted per one cup increase of tea consumption was not associated with knee OA (OR 1.11,95% CI: 0.79-1.55) using IVW with random effect. Genetic predisposition to tea consumption was not associated with hip OA (OR: 1.20, 95% CI: 0.84-1.71), RA (OR: 1.24 95% CI: 0.81-1.91), and OP (OR: 1.11, 95% CI: 0.89, 1.39). Following the sensitivity analysis, there was no potential pleiotropy.

Conclusion

According to our study, According to our study, there was no statistical power to confirm a causal relationship between tea consumption and the risk of knee OA, hip OA, RA, and OP.