Cardiometabolic traits mediating the effect of education on osteoarthritis risk: a Mendelian randomization study

Causal relationship between modifiable risk factors and knee osteoarthritis: a Mendelian randomization study

Background

While several risk factors for knee osteoarthritis (KOA) have been recognized, the pathogenesis of KOA and the causal relationship between modifiable risk factors and KOA in genetic epidemiology remain unclear. This study aimed to determine the causal relationship between KOA and its risk factors.

Methods

Data were obtained from published Genome-Wide Association study (GWAS) databases. A two-sample Mendelian randomization (MR) analysis was performed with genetic variants associated with risk factors as instrumental variables and KOA as outcome. First, inverse variance weighting was used as the main MR analysis method, and then a series of sensitivity analyses were conducted to comprehensively evaluate the causal relationship between them.

Results

Univariate forward MR analysis revealed that genetically predicted hypothyroidism, hyperthyroidism/thyrotoxicosis, educational level, income level, metabolic syndrome (MS), essential hypertension, height, hot drink temperature, diet (abstaining from sugar-sweetened or wheat products), and psychological and psychiatric disorders (stress, depression, and anxiety) were causally associated with KOA. Reverse MR exhibits a causal association between KOA and educational attainment. Multivariate MR analysis adjusted for the inclusion of potential mediators, such as body mass index (BMI), smoking, alcohol consumption, and sex, exhibited some variation in causal effects. However, hyperthyroidism/thyrotoxicosis had a significant causal effect on KOA, and there was good evidence that height, hypothyroidism, educational level, psychological and psychiatric disorders (stress, depression, and anxiety), and abstaining from wheat products had an independent causal relationship. The mediating effect of BMI as a mediator was also identified.

Conclusion

This study used MR to validate the causal relationship between KOA and its risk factors, providing new insights for preventing and treating KOA in clinical practice and for developing public health policies.

The role of inflammatory proteins in regulating the impact of lipid specifications on deep venous thrombosis: a two sample and mediated Mendelian randomization study

Objective

To investigate the potential mediating role of inflammatory proteins in the association between lipid species and Deep Venous Thrombosis (DVT).

Methods

A comprehensive analysis was conducted using pooled data from genome-wide association studies (GWAS), incorporating double-sample and reverse Mendelian randomization (MR) techniques, to identify the specific inflammatory proteins that act as intermediaries among 91 screened proteins in relation to deep vein thrombosis (DVT). Furthermore, a two-step MR approach was employed to quantify the proportion of DVT risk attributed to lipid effects mediated by these inflammatory proteins.

Results

The MR Analysis revealed that the two inflammatory proteins, as predicted by genetics, served as mediating factors in the impact of five lipids on DVT. No reverse effect of DVT was observed on 179 lipid species and 91 inflammatory proteins. In the case of TAG(58:7) and its influence on DVT, CCL20 played an intermediary role with an estimated proportion of 12.51% (ranging from 12% to 13%). SIRT2 exhibited a masking effect on DVT for PC(17:0/20:4) and PC(18:0/20:4), while CCL20 masked the impact of DVT on PC(14:0/18:2), PC(15:0/18:2), and PC(18:0/20:5).

Conclusions

In our study, we identified CCL20 as a crucial mediator in the association between TAG(58:7) and DVT, with a mediating proportion of 12.51% (12%-13%). Further investigations are warranted to explore other potential risk factors acting as mediators.

Iron, copper, zinc and magnesium on rheumatoid arthritis: a two-sample Mendelian randomization study

This study aimed to elucidate the causal genetic relationships between iron, copper, zinc, magnesium, and rheumatoid arthritis (RA). A two-sample Mendelian randomization (MR) analysis was conducted using the "TwoSampleMR" and "MendelianRandomization" packages in R. The random-effects inverse variance-weighted (IVW) method was used as the primary approach. We performed sensitivity analyses to test the reliability of the results. The random-effects IVW analysis revealed that there was no genetic causal relationship between iron (P = 0.429, odds ratio [OR] 95% confidence interval [CI] = 0.919 [0.746-1.133]), copper (P = 0.313, OR 95% CI = 0.973 [0.921-1.027]), zinc (P = 0.633, OR 95% CI = 0.978 [0.891-1.073]), or magnesium (P = 0.218, OR 95% CI = 0.792 [0.546-1.148]) and RA. Sensitivity analysis verified the reliability of the results. Therefore, there is no evidence to support a causal relationship between iron, copper, zinc, and magnesium intake at the genetic level and the development of RA.

Clinical phenotypes of comorbidities in end-stage knee osteoarthritis: a cluster analysis

OBJECTIVES

Comorbidities, as components of these heterogeneous features, often coexist with knee osteoarthritis, and are particularly prevalent in end-stage knee osteoarthritis. Here, we attempted to identify the different clinical phenotypes of comorbidities in patients with end-stage knee osteoarthritis by cluster analysis.

METHODS

A total of 421 inpatients diagnosed with end-stage knee osteoarthritis who underwent inpatient surgery were included in this cross-sectional study. 23 demographic, comorbidity, inflammatory immune and evaluation scale variables were collected. Systematic clustering after factor analysis and separate two-step cluster analysis were performed for individual comorbidity variables and all variables, respectively, to objectively identify the different clinical phenotypes of the study patients.

RESULTS

Four clusters were finally identified. Cluster 1 had the largest proportion of obese patients (93.8%) and hypertension was common (71.2%). Almost all patients in cluster 2 were depressed (95.8%) and anxiety disorders (94.7%). Cluster 3 combined patients with isolated end-stage knee osteoarthritis and a few comorbidities. Cluster 4 had the highest proportion of patients with rheumatoid arthritis (58.8%).

CONCLUSIONS

Patients with end-stage knee osteoarthritis may be classified into four different clinical phenotypes: "isolated end-stage knee osteoarthritis"; "obesity + hypertension"; "depression + anxiety"; and "rheumatoid arthritis", which may help guide individualized patient care and treatment strategies.

Phenome-wide causal associations between osteoarthritis and other complex traits through the latent causal variable analysis

BACKGROUND

Individuals with osteoarthritis present with comorbidities, and the potential causal associations remain incompletely elucidated. The present study undertook a large-scale investigation about the causality between osteoarthritis and variable traits, using the summary-level data of genome-wide association studies (GWAS).

METHODS

The present study included the summary-level GWS data of knee osteoarthritis, hip osteoarthritis, hip or knee osteoarthritis, hand osteoarthritis, and other 1355 traits. Genetic correlation analysis was conducted between osteoarthritis and other traits through cross-trait bivariate linkage disequilibrium score regression. Subsequently, latent causal variable analysis was performed to explore the causal association when there was a significant genetic correlation. Genetic correlation and latent causal variable analysis were conducted on the Complex Traits Genomics Virtual Lab platform ( https://vl.genoma.io/ ).

RESULTS

We found 133 unique phenotypes showing causal relationships with osteoarthritis. Our results confirmed several well-established risk factors of osteoarthritis, such as obesity, weight, BMI, and meniscus derangement. Additionally, our findings suggested putative causal links between osteoarthritis and multiple factors. Socioeconomic determinants such as occupational exposure to dust and diesel exhaust, extended work hours exceeding 40 per week, and unemployment status were implicated. Furthermore, our analysis revealed causal associations with cardiovascular and metabolic disorders, including heart failure, deep venous thrombosis, type 2 diabetes mellitus, and elevated cholesterol levels. Soft tissue and musculoskeletal disorders, such as hallux valgus, internal derangement of the knee, and spondylitis, were also identified to be causally related to osteoarthritis. The study also identified the putative causal associations of osteoarthritis with digestive and respiratory diseases, such as Barrett's esophagus, esophagitis, and asthma, as well as psychiatric conditions including panic attacks and manic or hyperactive episodes. Additionally, we observed osteoarthritis causally related to pharmacological treatments, such as the use of antihypertensive medications, anti-asthmatic drugs, and antidepressants.

CONCLUSION

Our study uncovered a wide range of traits causally associated with osteoarthritis. Further studies are needed to validate and illustrate the detailed mechanism of those causal associations.

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.

The association between lipid biomarkers and osteoarthritis based on the National Health and Nutrition Examination Survey and Mendelian randomization study

To explore the association between lipid markers and osteoarthritis (OA). First, the National Health and Nutrition Examination Survey (NHANES) database was used to screen participants with lipid markers, OA and relevant covariates, and logistic regression was used to analyze the association between lipid markers and OA; Then, under the theoretical framework of Mendelian randomization (MR), two-sample MR was performed using GWAS data of lipid markers and OA to explore the causal association between the two, which was analyzed by inverse variance weighting (IVW) method. Heterogeneity test, sensitivity analysis and pleiotropy analysis were also performed. The NHANES database screened a total of 3706 participants, of whom 836 had OA and 2870 did not have OA. When lipid markers were used as continuous variables, multivariate logistic results showed an association between HDL, LDL and OA (HDL, OR (95%):1.01 (1.00, 1.01); LDL, OR (95%):1.00 (0.99, 1.00)). When lipid markers were used as categorical variables, multivariate logistic results showed the fourth quartile result of 0.713 (0.513, 0.992) for LDL relative to the first quartile. In MR study, the results of the IVW method for TG, TL, HDL and LDL showed OR (95% CI) of 1.06 (0.97-1.16), 0.95 (0.85-1.06), 0.94 (0.86-1.02) and 0.89 (0.80-0.998) with P-values of 0.21, 0.37. 013, 0.046. The heterogeneity tests and multiplicity analyses showed P-values greater than 0.05, and sensitivity analyses showed no abnormal single nucleotide polymorphisms. Through NHANES database and MR analyses, LDL was found to be a protective factor for OA, while HDL still needs further study. Our results provide new biomarkers for preventive and therapeutic strategies for OA.

Late age at first birth is a protective factor for oesophageal cancer and gastro-oesophageal reflux: the evidence from the genetic study

Objective

The primary objective of this research endeavor was to examine the underlying genetic causality between the age at first birth (AFB) and four prevalent esophageal diseases, namely oesophageal obstruction (OO), oesophageal varices (OV), gastro-oesophageal reflux (GOR), and oesophageal cancer (OC).

Methods

We conducted a two-sample Mendelian randomization (MR) analysis to examine the causal association between AFB and four prevalent esophageal disorders. We employed eight distinct MR analysis techniques to evaluate causal relationships, encompassing random-effects inverse variance weighted (IVW), MR Egger, weighted median, simple mode, weighted mode, maximum likelihood, penalized weighted median, and fixed-effects IVW. The random-effects IVW method served as the primary approach for our analysis. Furthermore, we executed several sensitivity analyses to assess the robustness of the genetic causal inferences.

Results

The random-effects IVW analysis revealed a significant negative genetic causal association between AFB and both GOR (P < 0.001, Odds Ratio [OR] 95% Confidence Interval [CI] = 0.882 [0.828-0.940]) and OC (P < 0.001, OR 95% CI = 0.998 [0.998-0.999]). Conversely, there was insufficient evidence support to substantiate a genetic causal link between AFB and OO (P = 0.399, OR 95% CI = 0.873 [0.637-1.197]) or OV (P = 0.881, OR 95% CI = 0.978 [0.727-1.314]). The results of sensitivity analyses underscore the robustness and reliability of our MR analysis.

Conclusion

The findings of this investigation substantiate the notion that elevated AFB confers a protective effect against GOR and OC. In addition, no causative association was discerned between AFB and OO or OV at the genetic level.

Age at first birth, age at menopause, and risk of ovarian cyst: a two-sample Mendelian randomization study

Background

Increasing observational studies have indicated that hormonal reproductive factors were associated with ovarian cyst, a common gynecological disease. A two-sample Mendelian randomization (MR) was carried out by investigating the causality of reproductive factors including age at first birth (AFB), age at natural menopause (ANM), and age at menarche (AAM), and the risk of ovarian cyst (OC).

Method

Summary statistics were collected from a large genome-wide association study (GWAS), and we used a two-sample MR study to clarify the causal association between the exposure of AFB (N = 542,901), ANM (N = 69,360), and AAM (N = 29,346) and the outcome of the OC (N case = 20,750, N control = 107,564). We separately selected 51, 35, and 6 single-nucleotide polymorphisms (SNPs) as instrumental variables (IVs) for assaying the influence of AFB, ANM, and AAM on OC, respectively. Then, the causal relationship was tested through multiple approaches including an inverse-variance weighted method, an MR-Egger regression, and a weighted median method. In addition, the MR-PRESSO method was also used to verify the horizontal pleiotropy. Subsequently, we adjust the confounders for MR design.

Results

The MR analysis results showed that AFB was negatively associated with the OC (IVW Beta: -0.09, OR: 0.91, 95% CI: 0.86-0.96, p = 0.00185), and the greater AAM decreased the risk of OC (IVW Beta: -0.10, OR: 0.91, 95% CI: 0.82-0.99, p = 0.0376). However, ANM has a positive correlation with the OC (IVW Beta: 0.05, OR: 1.05, 95% CI: 1.03-1.08, p = 8.38 × 10-6). After adjusting BMI, alcohol intake frequency, and ever smoked, we also obtained a negative relationship between AFB and OC (p < 0.005). Meanwhile, we adjusted weight, alcohol intake frequency, and height, and then found a causal relationship between older AMN and an increased risk of OC (p < 0.005).

Conclusion

A causal effect of reproductive factors on the development of OC, affected by AFB, ANM, and AAM, was found convincingly. After adjusting the confounders, we also successfully found the substantial causal effect of younger AFB, younger AAM, and older ANM on an increased risk of OC.

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 effect of physical activity and sedentary behaviors on the risk of osteoarthritis: a univariate and multivariate Mendelian randomization study

There is still a lot of ambiguity about the link between physical activity (PA), sedentary behaviors (SBs) and osteoarthritis (OA). This study aimed to investigate the causal relationship of PA/SBs on the risk of OA. A univariate and multivariate Mendelian randomization (MR) analysis was conducted to investigate the causal effect of five PA phenotypes and three SB phenotypes on overall OA, knee OA, hip OA, total hip arthroplasty, and total knee arthroplasty (TKA). MR methods used were inverse-variance weighting, MR-Egger regressions, and weighted median. Sensitivity analysis examined horizontal pleiotropy and heterogeneity, and confirmed the reliability of the results. After false discovery rate, light do-it-yourself (DIY) activities decreased the risk for overall OA (OR: 0.32, 95% CI 0.16-0.65), and knee OA (OR: 0.26, 95% CI 0.12-0.51). Resulting in a decreased risk of walking for pleasure on overall OA (OR: 0.87, 95% CI 0.70-1.04) and knee OA (OR: 0.14, 95% CI 0.06-0.32) was also observed. Television viewing, however, significantly increased the risk of OA, knee OA, hip OA, and TKA. MVMR findings revealed independent causal impacts of walking for pleasure and watching television on overall and knee OA, taking into account BMI, smoking, and education. This study suggested that light DIY and walking for pleasure were beneficial for preventing OA, and the risk of OA and TKA increased with prolonged television watching.

Mendelian randomization analyses identified bioavailable testosterone mediates the effect of sex hormone-binding globulin on prostate cancer

OBJECTIVE

A better knowledge of the hormonal etiology of prostate cancer is essential for its prevention and treatment. The goal of this study was to provide causal estimates of the connection between sex hormone-binding globulin and prostate cancer and investigate the possible mediating function of other modifiable risk indicators.

METHODS

We used two-step, two-sample multivariable Mendelian randomization using single-nucleotide polymorphisms as instrumental variables for exposure and mediators. Single-nucleotide polymorphisms associated with sex hormone-binding globulin and bioavailable testosterone were screened via a genome-wide association study enrolling European-descent adult male individuals. Summary-level data for prostate cancer (79,148 cases and 61,106 controls) were extracted from the PRACTICAL consortium. The total effect of sex hormone-binding globulin on prostate cancer risk was decomposed into direct and indirect effects through the mediator, bioavailable testosterone. An inverse-variance-weighted method was the primary Mendelian randomization analysis method. Sensitivity analyses were performed via Mendelian randomization-Egger regression, heterogeneity test, pleiotropy test, and leave-one-out test. The directionality that exposure causes the outcome was verified using Mendelian randomization-Steiger test.

RESULTS

In the univariable Mendelian randomization analysis, genetically predicted higher sex hormone-binding globulin levels had a causal association with lower prostate cancer risk (odds ratio = 0.944, 95% confidence interval = 0.897-0.993, p = 0.027) and an inverse association with bioavailable testosterone level (odds ratio = 0.945, 95% confidence interval = 0.926-0.965, p = 1.62E-07) without controlling for other factors. Moreover, an increase of one standard deviation (59.5 pmol/L) in genetically predicted bioavailable testosterone level was significantly associated with a 22.0% increase in the overall prostate cancer risk (odds ratio = 1.220, 95% confidence interval = 1.064-1.398, p = 0.004) after adjusting for sex hormone-binding globulin level. The effect size ratio of bioavailable testosterone-mediated sex hormone-binding globulin to prostate cancer was further analyzed to clarify the importance of the mediating effect. Notably, the mediator bioavailable testosterone explained 19.28% (95% confidence interval = 10.76%, 73.78%) of the total effect of sex hormone-binding globulin level on prostate cancer risk.

CONCLUSION

The results support the potentially protective causal effect of genetically predicted higher sex hormone-binding globulin levels against prostate cancer with mediation by the modifiable risk factor, bioavailable testosterone. More research is needed to determine how this possible sex hormone-binding globulin-bioavailable testosterone-prostate cancer link works. Targeting sex hormone-binding globulin and bioavailable testosterone traits may be a valuable strategy for preventing prostate cancer.

Improved genetic prediction of the risk of knee osteoarthritis using the risk factor-based polygenic score

BACKGROUND

Polygenic risk score (PRS) analysis is used to predict disease risk. Although PRS has been shown to have great potential in improving clinical care, PRS accuracy assessment has been mainly focused on European ancestry. This study aimed to develop an accurate genetic risk score for knee osteoarthritis (OA) using a multi-population PRS and leveraging a multi-trait PRS in the Japanese population.

METHODS

We calculated PRS using PRS-CS-auto, derived from genome-wide association study (GWAS) summary statistics for knee OA in the Japanese population (same ancestry) and multi-population. We further identified risk factor traits for which PRS could predict knee OA and subsequently developed an integrated PRS based on multi-trait analysis of GWAS (MTAG), including genetically correlated risk traits. PRS performance was evaluated in participants of the Nagahama cohort study who underwent radiographic evaluation of the knees (n = 3,279). PRSs were incorporated into knee OA integrated risk models along with clinical risk factors.

RESULTS

A total of 2,852 genotyped individuals were included in the PRS analysis. The PRS based on Japanese knee OA GWAS was not associated with knee OA (p = 0.228). In contrast, PRS based on multi-population knee OA GWAS showed a significant association with knee OA (p = 6.7 × 10-5, odds ratio (OR) per standard deviation = 1.19), whereas PRS based on MTAG of multi-population knee OA, along with risk factor traits such as body mass index GWAS, displayed an even stronger association with knee OA (p = 5.4 × 10-7, OR = 1.24). Incorporating this PRS into traditional risk factors improved the predictive ability of knee OA (area under the curve, 74.4% to 74.7%; p = 0.029).

CONCLUSIONS

This study showed that multi-trait PRS based on MTAG, combined with traditional risk factors, and using large sample size multi-population GWAS, significantly improved predictive accuracy for knee OA in the Japanese population, even when the sample size of GWAS of the same ancestry was small. To the best of our knowledge, this is the first study to show a statistically significant association between the PRS and knee OA in a non-European population.

TRIAL REGISTRATION

No. C278.

Association Between Metabolically Different Adiposity Subtypes and Osteoarthritis: A Mendelian Randomization Study

OBJECTIVE

In this Mendelian randomization (MR) study, the objective was to investigate the causal effect of metabolically different adiposity subtypes on osteoarthritis.

METHODS

We performed 2-sample MR using summary-level data for osteoarthritis (10,083 cases and 40,425 controls) from a genome-wide association using the UK Biobank, and for site-specific osteoarthritis from the Arthritis Research UK Osteoarthritis Genetics consortium. We used 3 classes of genetic instruments, which all increase body mass index but are associated with different metabolic profiles (unfavorable, neutral, and favorable). Primary analysis was performed using inverse variance weight (IVW), with additional sensitivity analysis from different MR methods. We further applied a nonlinear MR using UK Biobank data to understand the nature of the adiposity-osteoarthritis relationship.

RESULTS

Greater metabolically unfavorable and metabolically neutral adiposity were associated with higher odds of osteoarthritis (IVW odds ratio [OR] 1.56 [95% confidence interval (95% CI) 1.31, 1.85] and OR 1.60 [95% CI 1.15, 2.23], respectively). The estimate for the association between metabolically favorable adiposity and osteoarthritis was similar, although with notable imprecision (OR 1.55 [95% CI 0.70, 3.41]). Using site-specific osteoarthritis, metabolically unfavorable, neutral, and favorable adiposity were all associated with higher odds of knee osteoarthritis (OR 1.44 [95% CI 1.04, 1.98], OR 2.28 [95% CI 1.04, 4.99], and OR 6.80 [95% CI 2.08, 22.19], respectively). We found generally consistent estimates with a wider confidence interval crossing the null from other MR methods. The nonlinear MR analyses suggested a nonlinear relationship between metabolically unfavorable adiposity and osteoarthritis (Pnonlinear  = 0.003).

CONCLUSION

Metabolic abnormalities did not explain the association between greater adiposity and the risk of osteoarthritis, which might suggest that the association is largely due to a mechanical effect on the joints.

The impact of obesity: a narrative review

Obesity is a disease with a major negative impact on human health. However, people with obesity may not perceive their weight to be a significant problem and less than half of patients with obesity are advised by their physicians to lose weight. The purpose of this review is to highlight the importance of managing overweight and obesity by discussing the adverse consequences and impact of obesity. In summary, obesity is strongly related to >50 medical conditions, with many of them having evidence from Mendelian randomisation studies to support causality. The clinical, social and economic burdens of obesity are considerable, with these burdens potentially impacting future generations as well. This review highlights the adverse health and economic consequences of obesity and the importance of an urgent and concerted effort towards the prevention and management of obesity to reduce the burden of obesity.

Body mass index, C-reactive protein, and pancreatic cancer: A Mendelian randomization analysis to investigate causal pathways

Aim

To explore whether C-reactive protein (CRP) mediates the risk of body mass index (BMI) in pancreatic cancer (PC) and calculate the mediate proportion of CRP in this possible mechanism.

Methods

Based on two-sample Mendelian randomization (TSMR), a two-step Mendelian randomization (TM) model was conducted to determine whether CRP was a mediator of the causal relationship between BMI and PC. The multivariable Mendelian randomization (MVMR) study was designed for mediating analysis and to calculate the mediating proportion mediated by CRP.

Results

BMI has a positive causal relationship with PC (n = 393 SNPs, OR = 1.484, 95% CI: 1.021-2.157, p< 0.05). BMI has a positive causal relationship with CRP (n = 179 SNPs, OR = 1.393, 95% CI: 1.320-1.469, p< 0.05). CRP has a positive causal relationship with PC (n = 54 SNPs, OR = 1.348, 95% CI: 1.004-1.809, p< 0.05). After adjusting CRP, BMI has no causal relationship with PC (n = 334 SNPs, OR = 1.341, 95% CI: 0.884-2.037, p< 0.05). After adjusting BMI, there was still a positive causal relationship between CRP and PC (n = 334 SNPs, OR = 1.441, 95% CI: 1.064-1.950, p< 0.05). The mediating effect of CRP was 29%.

Conclusions

In clinical practice, while actively advocating for weight loss among obese patients, we should focus on chronic inflammation levels in obese patients as well. In addition, anti-inflammatory dietary patterns and appropriate physical activity are important in preventing PC.

Polycystic ovary syndrome and 25-hydroxyvitamin D: A bidirectional two-sample Mendelian randomization study

BACKGROUND

Accumulating observational studies have indicated that vitamin D deficiency (serum 25-hydroxyvitamin D (25OHD) < 50 nmol/L) is common in women with polycystic ovary syndrome (PCOS). However, the direction and causal nature remain unclear. In this study, we aimed to investigate the causal association between PCOS and 25OHD.

METHODS

A bidirectional two-sample Mendelian randomization (MR) study was used to evaluate the causal association between PCOS and 25OHD. From the publicly available European-lineage genome-wide association studies (GWAS) summary statistics for PCOS (4,890 cases of PCOS and 20,405 controls) and 25OHD (n = 417,580), we selected 11 and 102 single nucleotide polymorphisms (SNPs) as instrumental variables (IVs), respectively. In univariate MR (uvMR) analysis, inverse-variance weighted (IVW) method was employed in the primary MR analysis and multiple sensitivity analyses were implemented. Additionally, a multivariable MR (mvMR) design was carried to adjust for obesity and insulin resistance (IR) as well.

RESULTS

UvMR demonstrated that genetically determined PCOS was negatively associated with 25OHD level (IVW Beta: -0.02, P = 0.008). However, mvMR found the causal effect disappeared when adjusting the influence of obesity and IR. Both uvMR and mvMR analysis didn't support the causal effect of 25OHD deficiency on risk of PCOS (IVW OR: 0.86, 95% CI: 0.66 ~ 1.12, P = 0.280).

CONCLUSION

Our findings highlighted that the casual effect of PCOS on 25OHD deficiency might be mediated by obesity and IR, and failed to find substantial causal effect of 25OHD deficiency on risk of PCOS. Further observational studies and clinical trials are necessary.

Appraising the causal role of smoking in multiple diseases: A systematic review and meta-analysis of Mendelian randomization studies

BACKGROUND

The causal association between cigarette smoking and several diseases remains equivocal. The purpose of this study was to appraise the causal role of smoking in a wide range of diseases by summarizing the evidence from Mendelian randomization (MR) studies.

METHODS

MR studies on genetic liability to smoking initiation or lifetime smoking (composite of smoking initiation, heaviness, duration, and cessation) in relation to circulatory system, digestive system, nervous system, musculoskeletal system, endocrine, metabolic, and eye diseases, and neoplasms published until February 15, 2022, were identified in PubMed. De novo MR analyses were performed using summary statistics data from genome-wide association studies. Meta-analysis was applied to combine study-specific estimates.

FINDINGS

Meta-analyses of findings of 29 published MR studies and 123 de novo MR analyses of 57 distinct primary outcomes showed that genetic liability to smoking (smoking initiation or lifetime smoking) was associated with increased risk of 13 circulatory system diseases, several digestive system diseases (including diverticular, gallstone, gastroesophageal reflux, and Crohn's disease, acute pancreatitis, and periodontitis), epilepsy, certain musculoskeletal system diseases (including fracture, osteoarthritis, and rheumatoid arthritis), endocrine (polycystic ovary syndrome), metabolic (type 2 diabetes) and eye diseases (including age-related macular degeneration and senile cataract) as well as cancers of the lung, head and neck, esophagus, pancreas, bladder, kidney, cervix, and ovaries, and myeloid leukemia. Smoking liability was associated with decreased risk of Parkinson's disease and prostate cancer.

INTERPRETATION

This study found robust evidence that cigarette smoking causes a wide range of diseases.

FUNDING

This work was supported by research grants from the Swedish Cancer Society (Cancerfonden), the Swedish Heart Lung Foundation (Hjärt-Lungfonden, 20210351), the Swedish Research Council for Health, Working Life and Welfare (Forte, 2018-00123), and the Swedish Research Council (Vetenskapsrådet, 2019-00977). Stephen Burgess is supported by Sir Henry Dale Fellowship jointly funded by the Wellcome Trust and the Royal Society (204623/Z/16/Z) and the National Institute for Health Research Cambridge Biomedical Research Centre (BRC-1215-20014). The views expressed are those of the authors and not necessarily those of the National Institute for Health Research or the Department of Health and Social Care.

Causal mediation analysis with multiple causally non-ordered and ordered mediators based on summarized genetic data

Causal mediation analysis investigates the mechanism linking exposure and outcome. Dealing with the impact of unobserved confounders among exposure, mediator and outcome is an issue of great concern. Moreover, when multiple mediators exist, this causal pathway intertwines with other causal pathways, rendering it difficult to estimate the path-specific effects. In this study, we propose a method (PSE-MR) to identify and estimate path-specific effects of an exposure (e.g. education) on an outcome (e.g. osteoarthritis risk) through multiple causally ordered and non-ordered mediators (e.g. body mass index and pack-years of smoking) using summarized genetic data, when the sequential ignorability assumption is violated. Specifically, PSE-MR requires a specific rank condition in which the number of instrumental variables is larger than the number of mediators. Furthermore, we illustrate the utility of PSE-MR by providing guidance for practitioners and exploring the mediation effects of body mass index and pack-years of smoking in the causal pathways from education to osteoarthritis risk. Additionally, the results of simulation reveal that the causal estimates of path-specific effects are almost unbiased with good coverage and Type I error properties. Also, we summarize the least number of instrumental variables for the specific number of mediators to achieve 80% power.

Causal Associations of Circulating Lipids with Osteoarthritis: A Bidirectional Mendelian Randomization Study

Osteoarthritis (OA) imposes an increasing social burden due to global activity limitations, especially among the aged. Links between circulating lipids and OA have been reported; however, confounding data from observational studies have hindered causal conclusions. We used Mendelian randomization (MR) approach to evaluate the genetic causal effects of circulating apolipoproteins and lipoprotein lipids on OA risk. Genetic instruments at the genome-wide significance level (p < 5 × 10−8) were selected from genome-wide association studies (n = 393,193−441,016 individuals). Summary-level OA data were obtained from the UK Biobank (39,427 cases, 378,169 controls). Bidirectional two-sample Mendelian randomization (MR) analyses used MR-Egger, weighted median, and MR-PRESSO for sensitivity analysis. Genetic predisposition to 1-SD increments of Apolipoprotein B (APOB), and low-density lipoprotein (LDL) was associated with a decreased risk of knee or hip OA (KHOA) (odds ratio (OR) = 0.925, 95% confidence interval (95% CI): 0.881−0.972, p = 0.002; OR = 0.898, 95% CI: 0.843−0.957, p = 0.001) and hip OA (HOA) (OR = 0.894; 95% CI: 0.832−0.961, p = 0.002; OR = 0.870 95% CI: 0.797−0.949, p = 0.002). Genetically predicted APOB showed an association with knee OA (KOA) (OR per SD increase, 0.930, 95% CI: 0.876−0.987, p = 0.016). The OR of KOA was 0.899 (95% CI: 0.835−0.968, p = 0.005) for a 1-SD increase in LDL. Apolipoprotein A1, high-density lipoprotein, and triglycerides showed no association. Inverse MR showed no causal effect of KOA, HOA, or KHOA on these serum lipids. Distinct protective genetic-influence patterns were observed for APOB and LDL on OA, offering new insights into relationships between lipids and OA risk and a better understanding of OA etiology.

No Genetic Causal Association Between Periodontitis and Arthritis: A Bidirectional Two-Sample Mendelian Randomization Analysis

Objectives

Periodontitis (PD) has been linked to arthritis in previous epidemiological observational studies; however, the results are inconclusive. It remains unclear whether the association between PD and arthritis is causal. The purpose of this study was to investigate the causal association of PD with arthritis, including rheumatoid arthritis (RA) and osteoarthritis (OA).

Methods

We performed a two-sample bidirectional Mendelian randomization (MR) analysis using publicly released genome-wide association studies (GWAS) statistics. The inverse-variance weighted (IVW) method was used as the primary analysis. We applied four complementary methods, including weighted median, weighted mode, MR-Egger regression and MR pleiotropy residual sum and outlier (MR-PRESSO) to detect and correct for the effect of horizontal pleiotropy.

Results

Genetically determined PD did not have a causal effect on OA (OR = 1.06, 95% CI: 0.99-1.15, P = 0.09) and RA (OR = 0.99, 95% CI: 0.87-1.13, P = 0.89). Furthermore, we did not find a significant causal effect of arthritis on PD in the reverse MR analysis. The results of MR-Egger regression, Weighted Median, and Weighted Mode methods were consistent with those of the IVW method. Horizontal pleiotropy was unlikely to distort the causal estimates according to the sensitivity analysis.

Conclusions

Our MR analysis reveals non-causal association of PD with arthritis, despite observational studies reporting an association between PD and arthritis.

Genetically Predicted Higher Educational Attainment Decreases the Risk of COVID-19 Susceptibility and Severity: A Mendelian Randomization Study

Background: Prior observational studies indicated that lower educational attainment (EA) is associated with higher COVID-19 risk, while these findings were vulnerable to bias from confounding factors. We aimed to clarify the causal effect of EA on COVID-19 susceptibility, hospitalization, and severity using Mendelian randomization (MR).

Methods: We identified genetic instruments for EA from a large genome-wide association study (GWAS) (n = 1,131,881). Summary statistics for COVID-19 susceptibility (112,612 cases and 2,474,079 controls), hospitalization (24,274 cases and 2,061,529 controls), and severity (8,779 cases and 1,001,875 controls) were obtained from the COVID-19 Host Genetics Initiative. We used the single-variable MR (SVMR) and the multivariable MR (MVMR) controlling intelligence, income, body mass index, vigorous physical activity, sedentary behavior, smoking, and alcohol consumption to estimate the total and direct effects of EA on COVID-19 outcomes. Inverse variance weighted was the primary analysis method. All the statistical analyses were performed using R software.

Results: Results from the SVMR showed that genetically predicted higher EA was correlated with a lower risk of COVID-19 susceptibility [odds ratio (OR) 0.86, 95% CI 0.84–0.89], hospitalization (OR 0.67, 95% CI 0.62–0.73), and severity (OR 0.67, 95% CI 0.58–0.79). EA still maintained its effects in most of the MVMR.

Conclusion: Educational attainment is a predictor for susceptibility, hospitalization, and severity of COVID-19 disease. Population with lower EA should be provided with a higher prioritization to public health resources to decrease the morbidity and mortality of COVID-19.

Effect of blood lipids and lipid-lowering therapies on osteoarthritis risk: A Mendelian randomization study

Background

We aimed to investigate the effects of blood lipids and lipid-lowering agents on osteoarthritis (OA) risk.

Materials and methods

We performed Mendelian randomization (MR) analyses to estimate the causal effect of blood low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and triglyceride (TG) levels on knee and hip OA. Single nucleotide polymorphisms (SNPs) were selected from large genome-wide association studies (GWASs) of individuals of European ancestry as genetic instruments for blood lipid levels. The associations of selected genetic instruments with knee and hip OA were estimated in a recent GWAS of the UK Biobank and arcOGEN datasets. Univariate and multivariate MR analyses were performed to detect and adjust for potential pleiotropy. Furthermore, genetic instruments in HMGCR, NPC1L1, and PCSK9 regions were used to mimic LDL-C-lowering effects of statin, ezetimibe, and evolocumab, respectively.

Results

Genetically determined LDL-C increments led to reduced risks of both knee OA (OR = 0.91 per 1-SD increment, 95% CI: 0.86-0.95, P = 6.3 × 10^-5) and hip OA (OR = 0.92, 95% CI: 0.85-0.99, P = 0.027). Multivariate MR analysis proved that the effect was independent of HDL-C, TG, and body mass index. TG increment was associated with reduced risks of hip OA in the univariate MR analysis; however, this was not supported by the multivariate MR analysis. Genetically proxied LDL-C-lowering effects of statins are related to increased risks of knee OA but not hip OA.

Conclusions

The findings suggested that LDL-C increments have independent protective effects on both knee and hip OA. LDL-C-lowering effects of statins may increase the risk of knee OA.

Causal Associations Between Educational Attainment and 14 Urological and Reproductive Health Outcomes: A Mendelian Randomization Study

Background: The impact of educational attainment (EA) on multiple urological and reproductive health outcomes has been explored in observational studies. Here we used Mendelian randomization (MR) to investigate whether EA has causal effects on 14 urological and reproductive health outcomes.

Methods: We obtained summary statistics for EA and 14 urological and reproductive health outcomes from genome-wide association studies (GWAS). MR analyses were applied to explore the potential causal association between EA and them. Inverse variance weighted was the primary analytical method.

Results: Genetically predicted one standard deviation (SD) increase in EA was causally associated with a higher risk of prostate cancer [odds ratio (OR) 1.14, 95% confidence interval (CI) 1.05–1.25, P = 0.003] and a reduced risk of kidney stone (OR 0.73, 95% CI 0.62–0.87, P < 0.001) and cystitis (OR 0.76, 95% CI 0.67–0.86, P < 0.001) after Bonferroni correction. EA was also suggestively correlated with a lower risk of prostatitis (OR 0.76, 95% CI 0.59–0.98, P = 0.037) and incontinence (OR 0.64, 95% CI 0.47–0.87, P = 0.004). For the bioavailable testosterone levels and infertility, sex-specific associations were observed, with genetically determined increased EA being related to higher levels of testosterone in men (β 0.07, 95% CI 0.04–0.10, P < 0.001), lower levels of testosterone in women (β −0.13, 95% CI−0.16 to−0.11, P < 0.001), and a lower risk of infertility in women (OR 0.74, 95% CI 0.64–0.86, P < 0.001) but was not related to male infertility (OR 0.79, 95% CI 0.52–1.20, P = 0.269) after Bonferroni correction. For bladder cancer, kidney cancer, testicular cancer, benign prostatic hyperplasia, and erectile dysfunction, no causal effects were observed.

Conclusions: EA plays a vital role in urological diseases, especially in non-oncological outcomes and reproductive health. These findings should be verified in further studies when GWAS data are sufficient.

Estimating causal effects of atherogenic lipid-related traits on COVID-19 susceptibility and severity using a two-sample Mendelian randomization approach

BACKGROUND

As the number of COVID-19 deaths continues to rise worldwide, the identification of risk factors for the disease is an urgent issue, and it remains controversial whether atherogenic lipid-related traits including serum apolipoprotein B, low-density lipoprotein (LDL)-cholesterol, and triglyceride levels, are risk factors. The aim of this study was to estimate causal effects of lipid-related traits on COVID-19 risk in the European population using a two-sample Mendelian randomization (MR) approach.

METHODS

We used summary statistics from a genome-wide association study (GWAS) that included 441,016 participants from the UK Biobank as the exposure dataset of lipid-related traits and from COVID-19 Host Genetics Initiative GWAS meta-analyses of European ancestry as the outcome dataset for COVID-19 susceptibility (32,494 cases and 1,316,207 controls), hospitalization (8316 cases and 1,549,095 controls), and severity (4792 cases and 1,054,664 controls). We performed two-sample MR analyses using the inverse variance weighted (IVW) method. As sensitivity analyses, the MR-Egger regression, weighted median, and weighted mode methods were conducted as were leave-one-out sensitivity analysis, the MR-PRESSO global test, PhenoScanner searches, and IVW multivariable MR analyses. A P value below 0.0055 with Bonferroni correction was considered statistically significant.

RESULTS

This MR study suggested that serum apolipoprotein B or LDL-cholesterol levels were not significantly associated with COVID-19 risk. On the other hand, we inferred that higher serum triglyceride levels were suggestively associated with higher risks of COVID-19 susceptibility (odds ratio [OR] per standard deviation increase in lifelong triglyceride levels, 1.065; 95% confidence interval [CI], 1.001-1.13; P = 0.045) and hospitalization (OR, 1.174; 95% CI, 1.04-1.33; P = 0.012), and were significantly associated with COVID-19 severity (OR, 1.274; 95% CI, 1.08-1.50; P = 0.004). Sensitivity and bidirectional MR analyses suggested that horizontal pleiotropy and reverse causation were unlikely.

CONCLUSIONS

Our MR study indicates a causal effect of higher serum triglyceride levels on a greater risk of COVID-19 severity in the European population using the latest and largest GWAS datasets to date. However, as the underlying mechanisms remain unclear and our study might be still biased due to possible horizontal pleiotropy, further studies are warranted to validate our findings and investigate underlying mechanisms.

Educational Attainment Decreases the Risk of COVID-19 Severity in the European Population: A Two-Sample Mendelian Randomization Study

Observational studies have reported that the severity of COVID-19 depends not only on physical conditions but also on socioeconomic status, including educational level. Because educational attainment (EA), which measures the number of years of schooling, is moderately heritable, we investigated the causal association of EA on the risk of COVID-19 severity using the Mendelian randomization (MR) approach. A two-sample MR analysis was performed using publicly available summary-level data sets of genome-wide association studies (GWASs). A total of 235 single-nucleotide polymorphisms (SNPs) were extracted as instrumental variables for the exposure of EA from the Social Science Genetic Association Consortium GWAS summary data of 766,345 participants of European ancestry. The effect of each SNP on the outcome of COVID-19 severity risk was obtained from the GWAS summary data of 1,059,456 participants of European ancestry gathered from the COVID-19 Host Genetics Initiative. Using inverse variance weighted method, our MR study shows that EA was significantly associated with a lower risk of COVID-19 severity (odds ratio per one standard deviation increase in years of schooling, 0.540; 95% confidence interval, 0.376–0.777, P = 0.0009). A series of sensitivity analyses showed little evidence of bias. In conclusion, we show for the first time using a two-sample MR approach the associations between higher EA and the lower risk of COVID-19 severity in the European population. However, the genetic or epidemiological mechanisms underlying the association between EA and the risk of COVID-19 severity remain unknown, and further studies are warranted to validate the MR findings and investigate underlying mechanisms.