Skip to main content
Download PDF

Causal relationship between sarcopenia with osteoarthritis and the mediating role of obesity: a univariate, multivariate, two-step Mendelian randomization study

BMC Geriatrics volume 24, Article number: 469 (2024) Cite this article

Abstract

Background

Recent genetic evidence supports a causal role for sarcopenia in osteoarthritis, which may be mediated by the occurrence of obesity or changes in circulating inflammatory protein levels. Here, we leveraged publicly available genome-wide association study data to investigate the intrinsic causal relationship between sarcopenia, obesity, circulating inflammatory protein levels, and osteoarthritis.

Methods

In this study, we used Mendelian randomization analyses to explore the causal relationship between sarcopenia phenotypes (Appendicular lean mass [ALM], Low hand-grip strength [LHG], and usual walking pace [UWP]) and osteoarthritis (Knee osteoarthritis [KOA], and Hip osteoarthritis [HOA]). Univariable Mendelian randomization (UVMR) analyses were performed using the inverse variance weighted (IVW) method, MR-Egger, weighted median method, simple mode, and weighted mode, with the IVW method being the primary analytical technique. Subsequently, the independent causal effects of sarcopenia phenotype on osteoarthritis were investigated using multivariate Mendelian randomization (MVMR) analysis. To further explore the mechanisms involved, obesity and circulating inflammatory proteins were introduced as the mediator variables, and a two-step Mendelian randomization analysis was used to explore the mediating effects of obesity and circulating inflammatory proteins between ALM and KOA as well as the mediating proportions.

Results

UVMR analysis showed a causal relationship between ALM, LHG, UWP and KOA [(OR = 1.151, 95% CI: 1.087–1.218, P = 1.19 × 10–6, PFDR = 7.14 × 10–6) (OR = 1.215, 95% CI: 1.004–1.470; P = 0.046, PFDR = 0.055) (OR = 0.503, 95% CI: 0.292–0.867; P = 0.013, PFDR = 0.027)], and a causal relationship between ALM, UWP and HOA [(OR = 1.181, 95% CI: 1.103–1.265, P = 2.05 × 10–6, PFDR = 6.15 × 10–6) (OR = 0.438, 95% CI: 0.226–0.849, P = 0.014, PFDR = 0.022)]. In the MVMR analyses adjusting for confounders (body mass index, insomnia, sedentary behavior, and bone density), causal relationships were observed between ALM, LHG, UWP and KOA [(ALM: OR = 1.323, 95%CI: 1.224- 1.431, P = 2.07 × 10–12), (LHG: OR = 1.161, 95%CI: 1.044- 1.292, P = 0.006), (UWP: OR = 0.511, 95%CI: 0.290- 0.899, P = 0.020)], and between ALM and HOA (ALM: OR = 1.245, 95%CI: 1.149- 1.348, P = 7.65 × 10–8). In a two-step MR analysis, obesity was identified to play a potential mediating role in ALM and KOA (proportion mediated: 5.9%).

Conclusions

The results of this study suggest that decreased appendicular lean mass, grip strength, and walking speed increase the risk of KOA and decreased appendicular lean mass increases the risk of HOA in patients with sarcopenia in a European population. Obesity plays a mediator role in the occurrence of KOA due to appendicular lean body mass reduction.

Peer Review reports

Introduction

Osteoarthritis is a chronic degenerative joint disease, with knee osteoarthritis (KOA) and hip osteoarthritis (HOA) being prevalent types. Clinically, the main manifestations include limited movements, arthralgia, and joint deformity, etc., which severely impact the patients' quality of life and may even lead to disability as the disease worsens [1, 2]. The pathological mechanisms of these two diseases are mainly related to the imbalance of cartilage matrix repair, but the causes of the two diseases remain unknown, typically involving the interaction of multiple factors, including genetic, mechanical and chemical factors [3, 4]. Mechanically, as important weight-bearing joints, the decline in lower limb muscle strength may lead to decreased joint stability and changes in lower limb force lines, resulting in increased joint stress, increasing the risk of osteoarthritis [5, 6]. Chemically, the reduction of skeletal muscle may cause hormonal regulatory disorders, thereby promoting the occurrence and deterioration of osteoarthritis [7]. Sarcopenia is an age-related chronic degenerative disease characterized by a decrease in overall muscle mass and strength [8]. The decline of the skeletal muscle system in sarcopenia patients may trigger the above mechanical and chemical factors, thereby increasing the risk of KOA and HOA.

Skeletal muscle has metabolic functions, and its reduction can lead to metabolic disorder in the body, affecting energy metabolism and inflammation regulation [9, 10], potentially causing obesity and exacerbating inflammatory responses. Obesity and chronic inflammation are considered as important promotive factors for osteoarthritis [11, 12]. Elevated levels of inflammatory biomarkers, including cytokines and soluble receptors, are closely associated with the onset of osteoarthritis [13]. Additionally, obesity is closely related with the occurrence and deterioration of KOA [14]. Although observational studies show a strong association between sarcopenia and KOA [15], it is unclear whether circulating inflammatory proteins and obesity mediate between sarcopenia and osteoarthritis. Additionally, observational studies have certain limitations in causal inference, making it difficult to establish causality for this association [16]. Therefore, new models are needed to explore the potential risk factors for osteoarthritis, in order to develop more precise and scientific intervention strategies.

Mendelian randomization (MR) is a statistical strategy that employs genetic variation as instrumental variables. Since genetic variation is typically undisturbed by environmental and behavioral factors, this method effectively reduces bias in estimating causal relationships [17]. MR analysis includes both Univariate Mendelian randomization (UVMR) and Multivariable Mendelian randomization (MVMR). MVMR, as an extension of UVMR, can explore the independent causal relationships between exposure factors and outcome factors, and it can also address the issue of multiple exposure factors and the mutual influence of similar exposure factors [18, 19]. This study collected large-scale data from genome-wide association study (GWAS) for UVMR and MVMR analyses, to access the causality of sarcopenia and related phenotypes with KOA and HOA, providing theoretical basis for their prevention and occurrence mechanism studies. Simultaneously, a Two-step Mendelian randomization (Two-step MR) model was used to investigate the mediating effect of obesity and circulating inflammatory proteins between sarcopenia and osteoarthritis, and to calculate the proportion of obesity and circulating inflammatory proteins in this possible mechanism.

Materials and methods

Study design

Due to the lack of published GWAS data for sarcopenia, according to the diagnostic criteria for sarcopenia recommended by the European Working Group on Sarcopenia in Older People (EWGSOP), appendicular lean mass (ALM), low hand-grip strength (LHG), and usual walking pace (UWP) were used instead of sarcopenia for MR analysis in this study [20].

In this study, we used UVMR and MVMR analyses to infer causal relationships between myasthenia gravis phenotypes and KOA and HOA. To be a valuable tool for causal inference in MR studies, genetic variation must satisfy three fundamental criteria: Assumption 1, genetic variation as an instrumental variable must be genuinely associated with exposure (sarcopenia phenotypes); assumption 2, exposure-outcome confounders have no effect on genetic variation; and assumption 3, genetic variation affects outcome through exposure (sarcopenia phenotypes) only, independent of other pathways. Previous studies have shown that KOA and HOA can be affected by obesity [21], insomnia [22], sedentary behavior [23], and bone mineral density [24]. To exclude the potential influence of these factors, this study used MVMR analysis to determine the causal associations of sarcopenia phenotypes with KOA and HOA independent of body mass index (BMI), insomnia, TV watching time, and femoral-neck bone mineral density (FN-BMD).

To further explore the mechanism between sarcopenia and osteoarthritis, the mediator variables obesity (BMI) and circulating inflammatory proteins were introduced to investigate their mediating effects between sarcopenia and osteoarthritis. The main characteristic of patients with sarcopenia is the decrease in muscle mass, therefore ALM was chosen as the exposure factor for the mediator analysis. Then, considering that KOA is the most common type of osteoarthritis, the study analyzed KOA as the outcome factor. In the mediation analysis, UVMR analysis was first used to explore the causal relationship between circulating inflammatory proteins and KOA, screening out inflammatory proteins that have a causal relationship with KOA. Subsequently, the mediation effects of obesity and the selected inflammatory proteins between ALM and KOA were investigated. Traditionally Two-Step Method was utilized for analyzing the mediation effects. we used the “product of coefficients approach” to evaluate the indirect impact of ALM on KOA through each potential mediator [25]. The indirect effect and proportion were obtained by the delta method [26]. In the first step, the causal relationship between the exposure factor (ALM) and mediator variables (BMI, circulating inflammatory proteins) was assessed. Gene variants strongly associated with the exposure factor and not directly related to the mediator and outcome variables were selected as instrumental variables. The Inverse Variance Weighted (IVW) method was employed as the primary analytical method. In the second step, the causal relationship between mediator variables (BMI, circulating inflammatory proteins) and the outcome variable (KOA) was evaluated. Gene variants strongly associated with the mediator variables and not directly related to the outcome variable were selected as instrumental variables. Finally, a comprehensive analysis was conducted by combining the results of the first and second steps to estimate the mediation effects of the mediator variables (BMI, circulating inflammatory proteins) between the exposure factor (ALM) and the outcome variable (KOA). A flowchart of the study design is shown in Fig. 1.

Fig. 1
figure 1

Mendelian randomization analysis flowchart

Full size image

Data source

All data used in this study are available in the public domain without additional ethical approval. The ALM data came from a GWAS study published by Pei YF et al. ALM was quantified in 450,243 participants of the UK Biobanking cohort using the Tanita BC 418ma body fat analyzer and its measurement accuracy was verified using DEXA method [27]; the LHG data was obtained from a GWAS meta-analysis study published by Garan G et al., which included 256,523 European descendants aged 60 years or older. Grip strength was measured using the Jamar J00105 hydraulic hand dynamometer, with maximum grip strength recorded in kilograms. Low grip strength was defined as < 30 kg for males and < 20 kg for females [28]; UWP data were sourced from the publicly available UK BioBank database, containing 459,915 European participants. Participants were assessed based on their responses to the following question: "How do you classify your regular walking speed? " (Slow Pace was considered below 3 mph, moderate pace corresponded to 3–4 mph, and speeds exceeding 4 mph were classified as fast) [29].

BMI data came from a study published by Sakaue S et al. which included a total of 359,983 samples [30]; Circulating inflammatory protein data were sourced from a study published by Zhao JH et al. involving 14,824 participants measured for 91 plasma proteins using the linkTarget platform for genome-wide protein quantitative trait locus (pQTL) analysis. Catlog numbers for the 91 circulating inflammatory proteins were GCST90274758 ~ GCST90274848 [31]; Television viewing time data were obtained from the UK BioBank database, including 437,887 European participants. Participants were assessed based on their responses to the following question: "How many hours do you spend watching television in a day? " [29]; the insomnia data were also resourced from the UK Biobank database, comprising 462,341 samples. Participants were assessed based on their responses to the following question: " Do you have difficulty falling asleep or waking up at night?". Response options included "Never/Rarely", "Sometimes", "Usually", or "Prefer not to answer". Participants who answered " Prefer not to answer " were defined as missing [29]; FN-BMD data were derived from a GWAS meta-analysis conducted by the Genetic Factors for Osteoporosis (GEFOS), involving 32,735 samples [32].

The KOA and HOA data were sourced from the largest-scale genome-wide meta-analysis to date. Inclusion and exclusion criteria for KOA and HOA were based on ICD codes. KOA data included 62,497 cases and 333,557 controls, while HOA data included 36,445 cases and 316,943 controls [33] (Table 1).

Table 1 Data information used in the MR analyses for this study
Full size table

Criteria for IV selection

The SNPs included in this analysis were selected based on the following criteria. To screen highly correlated SNPs, all selected instrumental variables reached genome-wide significance (significance threshold: P < 5 × 10–8). To avoid the impact of linkage disequilibrium, SNPs that did not meet the associated criteria (r2 = 0.001; 10,000 kb distance) were removed using R software. The F value, which represents a criterion for the strength of the relationship between instrumental variables and exposure, for all selected SNPs was greater than 10 and SNPs with F < 10 were excluded. In this study, the F statistic was calculated by the following formula: F = β exposure2/SE exposure2 [34]. The remaining SNPs were detected by the Phenoscanner database (http://www.phenoscanner.medschl.cam.ac.uk/), and SNPs associated with KOA, HOA confounders (e.g., BMI, percent body fat) were manually removed based on the results of the analysis.

Causal effect estimation and sensitivity analyses

In this study, the random-effects inverse-variance weighting (IVW) method, MR‒Egger method, weighted median method, simple mode, and weighted mode were used for the UVMR analysis [35,36,37].

Pleiotropy was tested by MR‒Egger regression. MR‒Egger regression analysis was used to detect whether the selected SNP had horizontal pleiotropy, and its regression intercept was used to evaluate the magnitude of pleiotropy. The null hypothesis of this hypothesis test is that the intercept term is zero, that is, there is no pleiotropy, and the rejection of the null hypothesis indicates that there is pleiotropy between the instrumental variables and the outcomes [36]. Heterogeneity between each SNP estimate was assessed using Cochran’s Q test, and if the results of the Q test were statistically significant and proved that the analysis results were significantly heterogeneous, we focused on the results of the random-effects IVW method [38]. The leave-one-out test was used to gradually eliminate each SNP, and the effect size change of the remaining SNPs was calculated [39]. If the results are greatly changed after excluding a single SNP, this SNP has a great impact on the results and is a sensitive SNP, so it needs to be eliminated, and MR analysis should be performed again. In addition, MR-PRESSO analysis was used to detect and remove outliers, and MR analysis was repeated after removing outliers [40].

To explore the independent causal associations of the three sarcopenia phenotypes with KOA and HOA, MVMR analysis was further performed using the IVW, MR‒Egger and weighted median methods. As sarcopenia is often accompanied by a loss of muscle mass and a decrease in physical function, in the MVMR analysis of this study, the sarcopenia phenotype combined with confounding factors was analyzed separately for patients with KOA and those with HOA. IVW was used as the primary method for the MVMR analysis. In addition, LASSO regression analysis was used to ensure the robustness of the results.

Because the outcome was a dichotomous variable, the MR analysis results needed to be converted into odds ratios (ORs) to evaluate the causal relationships of the sarcopenia phenotypes with KOA and HOA, and P < 0.05 was used as the threshold for determining statistical significance. In addition, to improve statistical rigor, this study applied False Discovery Rate (FDR) correction to adjust for multiple comparisons, thereby reducing the risk of Type I errors. Adjusted P-value of < 0.05 indicates a significant causal relationship. If the unadjusted p-value is < 0.05 but the FDR-adjusted p-value is > 0.05, it suggests a suggestive association [41].

Results

Results of MR analysis of sarcopenia phenotype with KOA and HOA

Results of UVMR analysis of sarcopenia phenotype with KOA and HOA

According to the screening criteria, in the UVMR analysis of the sarcopenia phenotype with KOA, a total of 504 SNPs were included in the analysis of ALM with KOA for analysis, with an F-statistic range of 29.75–708.02; 4 SNPs were included in the analysis of LHG with KOA for analysis, with an F-statistic range of 29.55–39.69; 14 SNPs were included in the analysis of UWP with KOA for analysis The F-statistic range was 29.79–42.61. In the UVMR analyses of the sarcopenia phenotype with HOA, the analyses of ALM with HOA incorporated 502 SNPs for analysis, with an F-statistic range of 29.75–609.31; and the analyses of LHG with HOA incorporated 10 SNPs for analysis, with an F-statistic range of 29.55–50.80; The analysis of UWP versus HOA included 16 SNPs for analysis, with an F-statistic range of 29.79–42.61. In this study, the F-statistic for each instrumental variable was > 10, indicating that the instrumental variables used in this study had a low bias. (Supplementary Material 1).

In the UVMR IVW analysis of ALM and KOA revealed a causal relationship between those variables (OR = 1.151, 95% CI: 1.087–1.218, P = 1.19 × 10–6, PFDR = 7.14 × 10–6), and the weighted median method was consistent with the IVW results. UVMR analysis of LHG and KOA revealed a causal relationship between the two variables (OR = 1.215, 95% CI: 1.004–1.470; P = 0.046, PFDR = 0.055). UVMR analysis of UWP and KOA revealed a causal relationship between those variables (OR = 0.503, 95% CI: 0.292–0.867; P = 0.013, PFDR = 0.027).

The UVMR IVW analysis revealed a causal relationship between ALM and HOA (OR = 1.181, 95% CI: 1.103–1.265, P = 2.05 × 10–6, PFDR = 6.15 × 10–6). The results of the weighted median, simple mode, and weighted mode analyses were consistent with the IVW results. The UVMR IVW analysis revealed no causal relationship between LHG and HOA (OR = 1.206, 95% CI: 0.909–1.600, P = 0.194, PFDR = 0.194) and a causal relationship between UWP and HOA (OR = 0.438, 95% CI: 0.226–0.849, P = 0.014, PFDR = 0.022) (Table 2).

Table 2 Results of univariate mendelian randomization analyses on sarcopenia phenotypes and osteoarthritis in knee and hip joints
Full size table

None of the MR‒Egger pleiotropy tests in the UVMR analysis detected potential horizontal pleiotropy, indicating that instrumental variables did not significantly affect outcomes through pathways other than exposure. According to the Cochran Q test, the analysis results of ALM with KOA, HOA and UWP with HOA showed significant heterogeneity, so the random-effects IVW was used as the main analysis method in this study. (Table 3) Leave-one-out cross-validation analysis revealed that no single-nucleotide polymorphisms affected the overall causal estimation. (Supplementary Material 2) Pleiotropy was found in the MR-PRESSO analysis of ALM and the two outcome factors, but the results still showed statistical significance after excluding abnormal SNPs.

Table 3 Results of MR‒Egger pleiotropy test and heterogeneity test of univariate Mendelian randomization analysis
Full size table

Results of MVMR analysis of sarcopenia phenotype with KOA and HOA

After adjusting for confounding factors, MVMR analysis of KOA revealed a causal relationship between the three sarcopenia phenotypes and KOA (ALM: OR = 1.323, 95%CI: 1.224- 1.431, P = 2.07 × 10–12) (LHG: OR = 1.161, 95%CI: 1.044- 1.292, P = 0.006) (UWP: OR = 0.511, 95%CI: 0.290- 0.899, P = 0.020). Since the results of the UVMR analysis of the association between LHG and HOA showed that there was no causal relationship between them, only ALM and UWP combined with confounding factors were included in the MVMR analysis in this study. In the MVMR IVW analysis, there was a direct causal effect of ALM on HOA (ALM: OR = 1.245, 95%CI: 1.149- 1.348, P = 7.65 × 10–8), but there was no causal effect of UWP on HOA (UWP: OR = 0.825, 95%CI: 0.419- 1.626, P = 0.579). There was no causal relationship between UWP and HOA after adjustment for confounders or internal adjustment between sarcopenia phenotypes. In addition, the MR-LASSO results were consistent with the above results, indicating the robustness of the results. (Table 4).

Table 4 Results of multivariable mendelian randomization analyses on sarcopenia phenotypes and osteoarthritis in knee and hip joints
Full size table

Mediator analysis

Results of UVMR analysis of circulating inflammatory proteins with KOA

In the UVMR analysis of circulating inflammatory proteins with KOA, the IVW analysis results showed that C-C motif chemokine 23 levels (CCL23) (OR= 0.959, 95%CI: 0.923- 0.995, P= 0.028, PFDR= 0.399), Fibroblast growth factor 19 levels (FGF19) (OR= 1.079, 95%CI: 1.010- 1.153, P= 0.025, PFDR= 0.399), Latency-associated peptide transforming growth factor β 1 levels (LAP TGF-β1) (OR= 0.839, 95%CI: 0.761- 0.924, P= 3.89×10-4, PFDR= 0.023), and Leukemia inhibitory factor receptor levels(LIFR) (OR=0.933, 95% CI: 0.883-0.985, P=3.89×10-4, PFDR=0.023) were causally associated with KOA. (See Supplementary Material 3 for details of the results of the analysis of 91 circulating inflammatory proteins with KOA).

Two-step MR analysis

According to the two-step mediation analysis method, the potential mediation effects of BMI, CCL23, FGF19, LAP-TGF-β1, and LIFR were assessed in the previously identified promotion of KOA risk by ALM. Testing the mediation effect of BMI in the relationship between ALM and KOA yielded a p-value of 0.049, with the proportion of mediated effect (IE_div_TE) at 0.059 and a confidence interval of the proportion of mediated effect ranging from 7.083×10-5 to 0.119, indicating the existence of a mediation effect. Testing the mediating effect of CCL23 in the relationship between ALM and KOA resulted in a p-value of 0.240, IE_div_TE of -0.009, and a confidence interval for the proportion of mediated effect ranging from -0.024 to 0.007, indicating no mediating effect. Testing the mediating effect of FGF19 in the relationship between ALM and KOA yielded a p-value of 0.257, IE_div_TE of -0.015, and a confidence interval for the proportion of mediated effect ranging from -0.043 to 0.013, indicating no mediating effect. Testing the mediating effect of LAP-TGF-β1 in the relationship between ALM and KOA yielded a p-value of 0.606, IE_div_TE of 0.013, and a confidence interval for the proportion of mediated effect ranging from -0.041 to 0.068, indicating no mediating effect. Testing the mediating effect of LIFR in the relationship between ALM and KOA resulted in a p-value of 0.128, IE_div_TE of -0.022, and a confidence interval for the proportion of mediated effect ranging from -0.052 to 0.008, indicating no mediating effect. (See Table 5 for analysis results of exposure with outcome, exposure with mediator, and mediator with outcome. See Table 6 for analysis results of total effect, direct effect, and mediated effect in the two-step MR analysis).

Table 5 The causal relationship between exposure and outcome, exposure and mediation, mediation and outcome
Full size table
Table 6 The total effect, direct effect and intermediary effect mediated by BMI
Full size table

Discussion

Amidst the global aging population, sarcopenia, as a degenerative disease of old age, is gradually increasing in incidence. The relationship between sarcopenia and arthritis has drawn increasing attention. Previous MR studies on sarcopenia and osteoarthritis mainly employed UVMR analysis [42,43,44], and they found a positive causal relationship between ALM and KOA, as well as HOA, and a negative causal relationship between UWP and KOA, as well as HOA, consistent with the findings of this study. Furthermore, the causal relationship between LHG and KOA was also found in the UVMR analysis of the present study. Considering that UVMR analysis cannot fully eliminate the effect of horizontal pleiotropy, MVMR was able to efficiently remove the effect of horizontal pleiotropy to access the independent causal relationship of exposure factors on outcome factors. Based on UVMR, this study further used MVMR for analysis. By adjusting for BMI, insomnia, television viewing time and FN-BMD, it investigated the independent causal relationship of sarcopenia phenotype on KOA and HOA. The results of the adjusted MVMR analysis indicated that there was no causal relationship between UWP and HOA, which might be due to the confounding effect modification in the UVMR analysis and the internal effect modification among the sarcopenia phenotypes masking the relationship between UWP and HOA. The findings of this study suggest that ALM, LHG, and UWP have causal relationships with KOA, and ALM has a causal relationship with HOA. In addition, the results of two-step MR analysis indicate that BMI plays a mediating role in the causal relationship between ALM and KOA.

ALM refers to the mass of body appendage muscles and bones, and in sarcopenia, a decrease in ALM mainly manifests as a decrease in appendage muscle mass. The results of this study suggest that the reduction of ALM increases the risk of KOA and HOA. Previous observational studies and animal experiments support our conclusion [45,46,47,48,49,50]. In the aspect of muscle function, the study shows that a decrease in grip strength and walking speed will increase the risk of KOA. Andrews J S et al. found in a European population that for every 1 standard deviation decrease in grip strength, the risk of KOA increases by 1.2 times [15]. Similarly, decreased grip strength has been found to be strongly associated with the development of osteoarthritis in East Asian populations [51]. Studies on the grip strength of elderly people have also shown that the grip strength of elderly people is positively correlated with the strength of other muscles in the whole body [52]. Therefore, a decrease in grip strength can indirectly reflect the deterioration of lower-limb muscles. As one of the most important muscles for lower-limb movement, the quadriceps femoris is crucial to the influence of KOA [53], and a decrease in strength will lead to increased pressure on the knee joint and then accelerate the progression of cartilage tissue degradation and arthritis [54]. In addition, walking speed can reflect the balance ability and nervous system function of the elderly, and is widely used to evaluate the overall physical function of the elderly. The results of this study suggest that a decrease in walking speed leads to an increased risk of KOA, but does not affect the occurrence and development of HOA. Previous observational studies have similarly found a negative correlation between walking speed and KOA, while walking speed was not associated with the development of HOA [55]. In the study of gait characteristics of sarcopenia patients, it has been found that the time of the support phase in the gait cycle of patients is significantly increased, and the delay of the support phase will lead to the joint bearing longer gravity and ground reaction force, thus increasing the risk of joint soft tissue injury [56].

The results of this study suggest that obesity plays a mediator role between the decline in ALM and the onset of KOA. The reduction of skeletal muscle in elderly individuals may trigger obesity through multiple mechanisms, including chronic inflammation, hormonal metabolic disorders, and decreased basal metabolic rate [57]. Additionally, obesity has long been recognized as one of the risk factors for KOA, as confirmed in the study conducted by Wang et al., which also found a positive causal relationship between obesity and KOA [58], consistent with the findings of this study. Obesity increases load on knee joint cartilage tissue and changes in body composition may also lead to metabolic disorders, thereby inducing the occurrence arthritis [59]. For sarcopenic patients, it is important to reduce the incidence of KOA by lowering BMI. It should be emphasized that during the process of reducing BMI, adipose tissue should be reduced as much as possible in sarcopenic patients while preserving muscle mass. The results of this study provide some new insights into the pathogenesis of KOA involving circulating inflammatory proteins. Although circulating inflammatory proteins may not mediate the occurrence of osteoarthritis resulting from skeletal muscle loss, this study found that high levels of CCL23, LAP-TGF-β1, and LIFR might contribute to reduce the risk of KOA, while high levels of FGF19 might increase the risk of KOA. These findings provide new research directions and insights for future studies on the interaction between circulating inflammatory proteins and osteoarthritis.

This study provides potential evidence that the sarcopenia phenotype contributes to the onset or development of KOA and HOA. Compared with observational studies, MR analysis has the advantage of avoiding reverse causation and confounding factors. However, this study has several limitations. First, the GWAS data used in this study came from a European database and were limited to individuals of European ancestry; it is unknown whether these results apply to other ethnic groups. Second, due to the current lack of GWAS data on the phenotype of sarcopenia itself, this study chose ALM, LHG and UWP to replace sarcopenia in exploring its causal relationships with KOA and HOA.

In conclusion, this study suggests that in a European population, the loss of appendicular lean body mass in patients with sarcopenia has a pathogenic effect on KOA and HOA, and a decrease in grip strength and walking speed in patients with sarcopenia may increase the risk of KOA. Meanwhile, obesity plays a mediator role in the occurrence of KOA due to appendicular lean body mass reduction.

Availability of data and materials

Data used in the present study are all publicly available. The Usual walking pace, insomnia and TV watching time data can be found at http://www.nealelab.is/uk-biobank/; Appendicular lean mass, Low hand grip and femoral-neck bone mineral density data can be found at https://gwas.mrcieu.ac.uk/. Knee osteoarthritis and Hip osteoarthritis data can be found at https://msk.hugeamp.org/downloads.html. All the data generated by the MR analysis is in the included Supplementary Material 1.

References

  1. Sharma L. Osteoarthritis of the Knee. N Engl J Med. 2021;384(1):51–9.

    Article  CAS  PubMed  Google Scholar 

  2. Hando BR, Gill NW 3rd. Osteoarthritis of the hip. New Engl J Med. 2008;358(5):534 author reply 534–535.

    Article  CAS  PubMed  Google Scholar 

  3. He Y, Li Z, Alexander PG, Ocasio-Nieves BD, Yocum L, Lin H, Tuan RS. Pathogenesis of Osteoarthritis: Risk Factors, Regulatory Pathways in Chondrocytes, and Experimental Models. Biology. 2020;9(8):194.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Katz JN, Arant KR, Loeser RF. Diagnosis and Treatment of Hip and Knee Osteoarthritis: A Review. JAMA. 2021;325(6):568–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Andriacchi TP. Dynamics of knee malalignment. Orthop Clin North Am. 1994;25(3):395–403.

    Article  CAS  PubMed  Google Scholar 

  6. Tanamas S, Hanna FS, Cicuttini FM, Wluka AE, Berry P, Urquhart DM. Does knee malalignment increase the risk of development and progression of knee osteoarthritis? A systematic review.  Arthritis Rheumatism. 2009;61(4):459–67.

    Article  PubMed  Google Scholar 

  7. Herrmann M, Engelke K, Ebert R, Müller-Deubert S, Rudert M, Ziouti F, Jundt F, Felsenberg D, Jakob F. Interactions between Muscle and Bone-Where Physics Meets Biology. Biomolecules. 2020;10(3):432.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Cruz-Jentoft AJ, Bahat G, Bauer J, Boirie Y, Bruyère O, Cederholm T, Cooper C, Landi F, Rolland Y, Sayer AA, et al. Sarcopenia: revised European consensus on definition and diagnosis. Age Ageing. 2019;48(1):16–31.

    Article  PubMed  Google Scholar 

  9. Umegaki H. Sarcopenia and frailty in older patients with diabetes mellitus. Geriatr Gerontol Int. 2016;16(3):293–9.

    Article  PubMed  Google Scholar 

  10. Petersen KF, Dufour S, Savage DB, Bilz S, Solomon G, Yonemitsu S, Cline GW, Befroy D, Zemany L, Kahn BB, et al. The role of skeletal muscle insulin resistance in the pathogenesis of the metabolic syndrome. Proc Natl Acad Sci USA. 2007;104(31):12587–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Michael JW, Schlüter-Brust KU, Eysel P. The epidemiology, etiology, diagnosis, and treatment of osteoarthritis of the knee. Deutsches Arzteblatt international. 2010;107(9):152–62.

    PubMed  PubMed Central  Google Scholar 

  12. Kim ME, Lee JS. Molecular Foundations of Inflammatory Diseases: Insights into Inflammation and Inflammasomes. Curr Issues Mol Biol. 2024;46(1):469–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Molnar V, Matišić V, Kodvanj I, Bjelica R, Jeleč Ž, Hudetz D, Rod E, Čukelj F, Vrdoljak T, Vidović D, et al. Cytokines and Chemokines Involved in Osteoarthritis Pathogenesis. Int J Mol Sci. 2021;22(17):9208.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Martin KR, Kuh D, Harris TB, Guralnik JM, Coggon D, Wills AK. Body mass index, occupational activity, and leisure-time physical activity: an exploration of risk factors and modifiers for knee osteoarthritis in the 1946 British birth cohort. BMC Musculoskelet Disord. 2013;14:219.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Andrews JS, Gold LS, Nevitt M, Heagerty PJ, Cawthon PM. Appendicular Lean Mass, Grip Strength, and the Development of Knee Osteoarthritis and Knee Pain Among Older Adults. ACR open rheumatology. 2021;3(8):566–72.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Skrivankova VW, Richmond RC, Woolf BAR, Davies NM, Swanson SA, VanderWeele TJ, Timpson NJ, Higgins JPT, Dimou N, Langenberg C, et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration. BMJ (Clinical research ed). 2021;375:n2233.

    PubMed  Google Scholar 

  17. Boef AG, Dekkers OM, le Cessie S. Mendelian randomization studies: a review of the approaches used and the quality of reporting. Int J Epidemiol. 2015;44(2):496–511.

    Article  PubMed  Google Scholar 

  18. Sanderson E, Smith GD, Windmeijer F, Bowden J. Corrigendum to: An examination of multivariable Mendelian randomization in the single-sample and two-sample summary data settings. Int J Epidemiol. 2020;49(3):1057.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Burgess S, Thompson SG. Multivariable Mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. 2015;181(4):251–60.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Sha T, Wang Y, Zhang Y, Lane NE, Li C, Wei J, Zeng C, Lei G. Genetic Variants, Serum 25-Hydroxyvitamin D Levels, and Sarcopenia: A Mendelian Randomization Analysis. JAMA Netw Open. 2023;6(8):e2331558.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Zhang L, Zhang W, Wu X, Cui H, Yan P, Yang C, Zhao X, Xiao J, Xiao C, Tang M, et al. A sex- and site-specific relationship between body mass index and osteoarthritis: evidence from observational and genetic analyses. Osteoarthritis Cartilage. 2023;31(6):819–28.

    Article  CAS  PubMed  Google Scholar 

  22. Ni J, Zhou W, Cen H, Chen G, Huang J, Yin K, Sui C. Evidence for causal effects of sleep disturbances on risk for osteoarthritis: a univariable and multivariable Mendelian randomization study. Osteoarthritis Cartilage. 2022;30(3):443–50.

    Article  CAS  PubMed  Google Scholar 

  23. Li X, Wang S, Liu W, Wu H, Zhu Y. Causal effect of physical activity and sedentary behaviors on the risk of osteoarthritis: a univariate and multivariate Mendelian randomization study. Sci Rep. 2023;13(1):19410.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Qu Y, Chen S, Han M, Gu Z, Zhang Y, Fan T, Zeng M, Ruan G, Cao P, Yang Q, et al. Osteoporosis and osteoarthritis: a bi-directional Mendelian randomization study. Arthritis Res Ther. 2023;25(1):242.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. VanderWeele TJ. Mediation Analysis: A Practitioner’s Guide. Annu Rev Public Health. 2016;37:17–32.

    Article  PubMed  Google Scholar 

  26. Carter AR, Gill D, Davies NM, Taylor AE, Tillmann T, Vaucher J, Wootton RE, Munafò MR, Hemani G, Malik R, et al. Understanding the consequences of education inequality on cardiovascular disease: mendelian randomisation study. BMJ (Clinical research ed). 2019;365:1855.

    PubMed  Google Scholar 

  27. Pei YF, Liu YZ, Yang XL, Zhang H, Feng GJ, Wei XT, Zhang L. The genetic architecture of appendicular lean mass characterized by association analysis in the UK Biobank study. Commun Biol. 2020;3(1):608.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Jones G, Trajanoska K, Santanasto AJ, Stringa N, Kuo CL, Atkins JL, Lewis JR, Duong T, Hong S, Biggs ML, et al. Genome-wide meta-analysis of muscle weakness identifies 15 susceptibility loci in older men and women. Nat Commun. 2021;12(1):654.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, Downey P, Elliott P, Green J, Landray M, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12(3):e1001779.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Sakaue S, Kanai M, Tanigawa Y, Karjalainen J, Kurki M, Koshiba S, Narita A, Konuma T, Yamamoto K, Akiyama M, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. 2021;53(10):1415–24.

    Article  CAS  PubMed  Google Scholar 

  31. Zhao JH, Stacey D, Eriksson N, Macdonald-Dunlop E, Hedman ÅK, Kalnapenkis A, Enroth S, Cozzetto D, Digby-Bell J, Marten J, et al. Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets. Nat Immunol. 2023;24(9):1540–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Zheng HF, Forgetta V, Hsu YH, Estrada K, Rosello-Diez A, Leo PJ, Dahia CL, Park-Min KH, Tobias JH, Kooperberg C, et al. Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture. Nature. 2015;526(7571):112–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Boer CG, Hatzikotoulas K, Southam L, Stefánsdóttir L, Zhang Y. Coutinho de Almeida R, Wu TT, Zheng J, Hartley A, Teder-Laving M et al: Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations. Cell. 2021;184(18):4784-4818.e4717.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Bowden J, Del Greco MF, Minelli C, Zhao Q, Lawlor DA, Sheehan NA, Thompson J, Davey Smith G. Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumption. Int J Epidemiol. 2019;48(3):728–42.

    Article  PubMed  Google Scholar 

  35. Slob EAW, Burgess S. A comparison of robust Mendelian randomization methods using summary data. Genet Epidemiol. 2020;44(4):313–29.

    Article  PubMed  PubMed Central  Google Scholar 

  36. 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(2):512–25.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 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(4):304–14.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Greco MF, Minelli C, Sheehan NA, Thompson JR. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat Med. 2015;34(21):2926–40.

    Article  Google Scholar 

  39. Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity Analyses for Robust Causal Inference from Mendelian Randomization Analyses with Multiple Genetic Variants. Epidemiology. 2017;28(1):30–42.

    Article  PubMed  Google Scholar 

  40. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125(1–2):279–84.

    Article  CAS  PubMed  Google Scholar 

  42. Zhang L, Zhang C, Zhang J, Liu A, Wang P, Xu J. A Bidirectional Mendelian Randomization Study of Sarcopenia-Related Traits and Knee Osteoarthritis. Clin Interv Aging. 2023;18:1577–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Yang J, Liu P, Wang S, Jiang T, Zhang Y, Liu W. Causal relationship between sarcopenia and osteoarthritis: a bi-directional two-sample mendelian randomized study. Eur J Med Res. 2023;28(1):327.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Jia XM, Deng TT, Su H, Shi HJ, Qin H, Yu GC, Yin Y, Liu FJ, Shi B. Genetic causality and site-specific relationship between sarcopenia and osteoarthritis: a bidirectional Mendelian randomization study. Front Genet. 2023;14:1340245.

    Article  PubMed  Google Scholar 

  45. Suh DH, Han KD, Hong JY, Park JH, Bae JH, Moon YW, Kim JG. Body composition is more closely related to the development of knee osteoarthritis in women than men: a cross-sectional study using the Fifth Korea National Health and Nutrition Examination Survey (KNHANES V-1, 2). Osteoarthritis Cartilage. 2016;24(4):605–11.

    Article  CAS  PubMed  Google Scholar 

  46. Long H, Xie D, Zeng C, Wei J, Wang Y, Yang T, Xu B, Qian Y, Li J, Wu Z, et al. Association between body composition and osteoarthritis: A systematic review and meta-analysis. Int J Rheum Dis. 2019;22(12):2108–18.

    Article  PubMed  Google Scholar 

  47. Chen Y, Sun Y, Pan X, Ho K, Li G. Joint distraction attenuates osteoarthritis by reducing secondary inflammation, cartilage degeneration and subchondral bone aberrant change. Osteoarthritis Cartilage. 2015;23(10):1728–35.

    Article  CAS  PubMed  Google Scholar 

  48. Rehan Youssef A, Longino D, Seerattan R, Leonard T, Herzog W. Muscle weakness causes joint degeneration in rabbits. Osteoarthritis Cartilage. 2009;17(9):1228–35.

    Article  CAS  PubMed  Google Scholar 

  49. Egloff C, Sawatsky A, Leonard T, Hart DA, Valderrabano V, Herzog W. Effect of muscle weakness and joint inflammation on the onset and progression of osteoarthritis in the rabbit knee. Osteoarthritis Cartilage. 2014;22(11):1886–93.

    Article  CAS  PubMed  Google Scholar 

  50. Xu J, She G, Gui T, Hou H, Li J, Chen Y, Zha Z. Knee muscle atrophy is a risk factor for development of knee osteoarthritis in a rat model. J Orthop Translat. 2020;22:67–72.

    Article  CAS  PubMed  Google Scholar 

  51. Chang SY, Han BD, Han KD, Park HJ, Kang S. Relation between Handgrip Strength and Quality of Life in Patients with Arthritis in Korea: The Korea National Health and Nutrition Examination Survey, 2015–2018. Medicina (Kaunas, Lithuania). 2022;58(2):172.

    Article  PubMed  Google Scholar 

  52. Porto JM, Nakaishi APM, Cangussu-Oliveira LM, Freire Júnior RC, Spilla SB, Abreu DCC. Relationship between grip strength and global muscle strength in community-dwelling older people. Arch Gerontol Geriatr. 2019;82:273–8.

    Article  PubMed  Google Scholar 

  53. Takagi S, Omori G, Koga H, Endo K, Koga Y, Nawata A, Endo N. Quadriceps muscle weakness is related to increased risk of radiographic knee OA but not its progression in both women and men: the Matsudai Knee Osteoarthritis Survey. Knee surgery, sports traumatology, arthroscopy : official journal of the ESSKA. 2018;26(9):2607–14.

    Article  PubMed  Google Scholar 

  54. Bennell KL, Hunt MA, Wrigley TV, Lim BW, Hinman RS: Muscle and exercise in the prevention and management of knee osteoarthritis: an internal medicine specialist's guide. Med Clin North Am 2009;93(1):161–177, xii.

  55. Purser JL, Golightly YM, Feng Q, Helmick CG, Renner JB, Jordan JM. Association of slower walking speed with incident knee osteoarthritis-related outcomes. Arthritis Care Res. 2012;64(7):1028–35.

    Article  Google Scholar 

  56. Fan Y, Zhang B, Huang G, Zhang G, Ding Z, Li Z, Sinclair J, Fan Y. Sarcopenia: Body Composition and Gait Analysis. Frontiers in aging neuroscience. 2022;14: 909551.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Müller MJ, Lagerpusch M, Enderle J, Schautz B, Heller M, Bosy-Westphal A. Beyond the body mass index: tracking body composition in the pathogenesis of obesity and the metabolic syndrome. Obesity reviews : an official journal of the International Association for the Study of Obesity. 2012;13(Suppl 2):6–13.

    Article  PubMed  Google Scholar 

  58. Wang C, Zhu Y, Liu Z, Long H, Ruan Z, Zhao S. Causal associations of obesity related anthropometric indicators and body compositions with knee and hip arthritis: A large-scale genetic correlation study. Front Endocrinol. 2022;13:1011896.

    Article  Google Scholar 

  59. Kulkarni K, Karssiens T, Kumar V, Pandit H. Obesity and osteoarthritis. Maturitas. 2016;89:22–8.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to sincerely thank the original GWASs and the related consortiums for sharing and managing the summary statistics.

Funding

Philosophy and Social Science Research Projects, project number: 2023BTY016.

Author information

Authors and Affiliations

  1. College of Physical Education, Henan Normal University, Xinxiang, 453007, China

    Zicheng Jin, Rui Wang, Lishuang Wan & Yuzhou Li

  2. College of Music and Dance, Henan Normal University, Xinxiang, 453007, China

    Linzi Jin

Authors
  1. Zicheng Jin
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Rui Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Linzi Jin
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Lishuang Wan
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yuzhou Li
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Conceptualization and methodology: Zicheng Jin and Yuzhou Li; data collection and analysis: Zicheng Jin and Rui Wang; validation: Linzi Jin and Lishuang Wan; writing—original draft preparation: Zicheng Jin, Rui Wang, Linzi Jin, Lishuang Wan; writing—review and editing: Zicheng Jin and Yuzhou Li; supervision: Yuzhou Li. All the authors have read and approved the final version of the manuscript.

Corresponding author

Correspondence to Yuzhou Li.

Ethics declarations

Ethics approval and consent to participate

All of the studies that were referenced in the genome-wide association studies we cited had obtained the necessary approvals from the appropriate review board and informed consent from participants. And our research obtained ethics approval from the Henan Normal University Ethics Committee and informed consent from all subjects or their legal guardians.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

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

Jin, Z., Wang, R., Jin, L. et al. Causal relationship between sarcopenia with osteoarthritis and the mediating role of obesity: a univariate, multivariate, two-step Mendelian randomization study. BMC Geriatr 24, 469 (2024). https://doi.org/10.1186/s12877-024-05098-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12877-024-05098-8

Keywords

BMC Geriatrics

ISSN: 1471-2318

Contact us