Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 08 January 2024
Sec. Statistical Genetics and Methodology
This article is part of the Research Topic Mendelian Randomization: the Study of Causal Relationships between Tumors and Underlying Pathogenic Factors View all 14 articles

Genetic causality and site-specific relationship between sarcopenia and osteoarthritis: a bidirectional Mendelian randomization study

Xue-Min Jia,Xue-Min Jia1,2Ting-Ting DengTing-Ting Deng3Hang SuHang Su4Hao-Jun ShiHao-Jun Shi5Hao QinHao Qin2Gong-Chang Yu
Gong-Chang Yu1*Ying Yin
Ying Yin6*Fan-Jie Liu,
Fan-Jie Liu1,2*Bin Shi,
Bin Shi1,2*
  • 1Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, Shandong, China
  • 2Neck-Shoulder and Lumbocrural Pain Hospital of Shandong First Medical University, Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, Shandong, China
  • 3College of Traditional Chinese Medicine, Shandong University of Traditional Chinese Medicine, Jinan, China
  • 4College of Rehabilitation, Shandong University of Traditional Chinese Medicine, Jinan, Shandong, China
  • 5School of TCM, Macau University of Science and Technology, Macau, China
  • 6Affiliated Hospital of Shandong University of Traditional Chinese Medicine, Jinan, Shandong, China

Background: Previous studies demonstrated a controversial relationship between sarcopenia (SP) and osteoarthritis (OA) and their genetic causality is unclear. Thus, we conducted a Mendelian randomization (MR) analysis to evaluate the possible causal association between sarcopenia-related traits (appendicular lean mass (ALM), grip strength, usual walking pace) and OA.

Method: We used pooled genetic data from the UK Biobank for ALM(n = 450,243), left-hand grip strength (n = 461,026), right-hand grip strength (n = 461,089) and usual walking pace (n = 459,915). Moreover, summary statistics for OA were obtained from the latest study conducted by the Genetics of Osteoarthritis Consortium, including all OA (n = 826,690), hand OA (n = 303,7782), hip OA (n = 353,388) and knee OA (n = 396,054). The primary method for estimating causal effects was the inverse-variance weighted (IVW) method, with the utilizing of false discovery rate adjusted p values (PFDR). Additional MR methods such as MR-Egger regression, MR pleiotropy residual sum and outlier (MR-PRESSO), weighted median were employed as supplementary analyses.

Results: We discovered ALM (odds ratio (OR) = 1.103, 95% confidence interval (CI) = 1.052–1.156, PFDR = 2.87E-04), hand grip strength (left, IVW OR = 0.823, 95% CI = 0.712 to 0.952, PFDR = 0.020; right, OR = 0.826, 95% CI = 0.718 to 0.950, PFDR = 0.020), and usual walking pace (OR = 0.339, 95% CI = 0.204 to 0.564, PFDR = 2.38E-04) were causally associated with OA risk. In the reverse MR analysis, we identified a causal effect of OA on ALM (β = −0.258, 95% CI = −0.369 to 0.146, PFDR = 0.6.07E-06), grip strength (left, β = −0.064, 95% CI = −0.104 to 0.024, PFDR = 0.002; right, β = −0.055, 95% CI = −0.095 to 0.014, PFDR = 0.008), and usual walking pace (β = −0.104, 95% CI = −0.147 to 0.061, PFDR = 1.61E-05).

Conclusion: This present study suggests an obvious causality of SP on OA, with condition exhibiting site-specific effects, while evidence was also provided for the causal effect of OA on SP.

Introduction

Sarcopenia (SP) and osteoarthritis (OA) as degenerative musculoskeletal diseases (DMD) emerged as major challenges for the aging population (Yin et al., 2023). SP is a muscle disease (muscle failure) (Cruz-Jentoft et al., 2019), characterized by an accelerated loss of muscle mass and function (Cruz-Jentoft and Sayer, 2019). Currently, about 50 million people worldwide suffer from SP, and its prevalence increases with age (Hida et al., 2014). Some studies have shown that the global prevalence of SP in people over 60 years old ranges from 10.00% to 27.00%, and in people over 80 years old is as high as 50.00% (Therakomen et al., 2020; Petermann-Rocha et al., 2022). OA is a degenerative disease with clinical manifestations of chronic pain, joint stiffness, and swelling (Cho et al., 2021), which afflicts more than 500 million people worldwide and has become the leading cause of chronic pain and disability in older adults (Wen and Xiao, 2022). As SP and OA are often diagnosed as comorbidities clinically, epidemiologic studies are increasingly examining the relationship between these two prevalent diseases.

Several studies have demonstrated a significant interaction between SP and OA. SP or its related traits are likely to be associated with outcomes in OA (predominantly OA of the knee). James S Andrews et al. showed that ALM and grip strength may be related to the development of KOA in older men (Andrews et al., 2021). On the other hand, a systematic review and meta-analysis showed that the prevalence of sarcopenia was more than twice as high in patients with osteoarthritis of the knee compared with controls (Pegreffi et al., 2023). However, there are also some studies do not agree with the aforementioned notion, and they consider the interconnection of the two to be highly controversial (Jones et al., 2021; Mezian et al., 2021; Tzartza et al., 2023). In addition, for ethical and practical purposes, a causal association between the two diseases could not be proved by a randomized controlled trial (RCT).

Mendelian randomization (MR) is a data analysis technique used to evaluate etiologic inferences in epidemiologic studies. The technique utilizes genetic variation as an instrumental variable (IV) to estimate reliable causal associations between exposures and outcomes (Burgess et al., 2019; Richmond and Davey Smith, 2022). Based on a two-sample bidirectional MR framework, we examined the potential causality of all OA, hand OA, hip OA, and knee OA with SP-related traits (appendicular lean mass (ALM), hand grip strength (left), hand grip strength (right), and usual walking pace).

Materials and methods

Study design overview

Figure 1 illustrates the design of our bidirectional MR study. We first estimated the causal effect of SP-related traits on OA and then assessed the causal effect of OA on SP-related traits. Genetic variants were considered as IVs only if they met the following three strict core assumptions. First, genetic variants were highly correlated with exposure. Second, genetic variants are not associated with confounding factors. Finally, genetic variation cannot act directly on the outcome, its effect on the outcome can only be reflected by exposure.

FIGURE 1
www.frontiersin.org

FIGURE 1. Generalized diagram of the study process. Abbreviations: OA: osteoarthritis; MR: Mendelian randomization;PFDR: The inverse-variance weighted (IVW) method with false discovery rate (FDR) adjusted p values.

Data sources for sarcopenia-related traits

All summary-level genetic data for three SP-related traits were obtained from the UK Biobank (UKB). UKB is a large-scale repository of biomedical genes and information resources, containing about half a million people. The repository includes samples of volunteers’ genetic information, lifestyle choices, and pedigree data (Sudlow et al., 2015). Identification of ALM, hand grip strength, and usual walking pace as consensus diagnostic criteria for SP was based on the report of the European Working Group on Sarcopenia in Older People (EWGSOP) (Cruz-Jentoft et al., 2010). In terms of the ALM, pooled data were analyzed for 450,243 UKB cohort participants, and ALM-related values were quantified and adjusted (Pei et al., 2020). Grip strength data were obtained from UKB’s summary of hand grip strength (left and right) for 461,026 and 461,089 European ethnicity (Mitchell et al., 2019), calibrated to hand size, and adjusted for factors such as age and sex. Finally, for the “usual walking pace”, data was similarly summarized from the GWAS summary data collected by the UKB for 459,915 European populations (Mitchell et al., 2019).

Data sources for all and site-specific OA

Summary statistics for all OA and its specific sites were obtained from the Genetics of Osteoarthritis (GO) Consortium’s GWAS meta-analysis involving 826,690 individuals (177,517 patients with OA and 649,173 controls) from nine different populations (Boer et al., 2021). OA was defined by the GO based on self-reported status, hospital diagnosis, the 10th edition of the International Classification of Diseases (ICD-10) code, or TREAT-OA Consortium-defined radiology definition. The study identified 100 independently associated risk variants in 11 OA phenotypes and is the most recent and comprehensive GWAS analysis known for OA. We selected four phenotypes for analysis based on research needs: all OA (n = 826,690), hand OA (n = 303,782), hip OA (n = 353,388) and knee OA (n = 396,054).

Table 1 shows additional details such as phenotypes of all study participants.

TABLE 1
www.frontiersin.org

TABLE 1. Data sources used in this study.

Selection of IVs

We first ensured that the SNPs of the three SP-related traits were genome-wide significant (p < 5E-08). Furthermore, to exclude the SNPs of strong linkage disequilibrium, we carried out an aggregation with r2 < 0.001 and kb = 10,000 processes. Then, we also assessed the strength of IVs based on calculated R2 [R2 = 2 * MAF*(1-MAF) *beta (Cruz-Jentoft et al., 2019)] and F-statistics [F = R (Cruz-Jentoft et al., 2019)/(1-R2) * (N-K-1)/K] for each SNP (Burgess et al., 2011). In order to reduce bias, the IVs with F-statistics less than 10 were excluded (Hemani et al., 2018a), and the SNPs and their related data such as F-statistics and R2 for subsequent causal analysis are listed in Supplementary Table S1. In the reverse MR analysis, the SNP screening process was consistent with the aforementioned procedure.

MR analysis

This study used R software (Version 4.3.0) and the two-sample MR package (Version 0.5.7) for data analysis (Yavorska and Burgess, 2017). We used the random-effects IVW method as the primary analytical method for MR estimation. This method was considered the most reliable in the absence of indications of directional pleiotropy in the selected IVs (Holmes et al., 2017). In addition, several sensitivity analyses were carried out, including weighted median (WM), MR-Egger, MR pleiotropy residual sum, and outlier (MR-PRESSO) test. Within this framework, the WM selected the median estimate to calculate the causal effects (Burgess et al., 2017). The MR-Egger regression method effectively tests the null causality hypothesis and gives consistent estimates of causality even if no genetic variation is valid. Additionally, the MR-Egger regression method is robust to horizontal pleiotropy (Bowden et al., 2015). The MR-PRESSO method detects pleiotropy, removes potentially pleiotropic IVs (outliers), and provides outlier-adjusted estimates (Verbanck et al., 2018). In order to correct the problem of multiple testing relatively gently, the p-value was adjusted using the false discovery rate (FDR) which is also called q value in the main IVW MR analysis, and the significance of causal inference was set to less than 0.05 (Chen et al., 2021). In addition, the statistical power was calculated with an online sample size and power calculator for MR (https://sb452.shinyapps.io/power/), and the results are shown in Supplementary Table S4.

Heterogeneity was assessed using Cochran’s Q test, and p < 0.05 was considered statistically significant (Bowden et al., 2016). Moreover, we performed pleiotropic tests using MR-Egger intercept test and MR-PRESSO global testing to ensure that IVs do not influence the risk of outcome through other confounding factors or other biological pathways unrelated to exposure (Hemani et al., 2018b). Furthermore, we performed leave-one-out analyses to ensure the reliability of associations with individual SNPs.

Results

Overview

Based on the inverse-variance weighted (IVW) method, we observed significant evidence of a bi-directionally causal relationship between SP and OA. After FDR correction, most of the meaningful results were still retained. Additionally, the results were corroborated using other MR analysis methods.

The causal effects of SP-related traits on OA

After accounting for the independence of the genetic variation, none of the IVs for ALM, hand grip strength (left), hand grip strength (right), and usual walking pace that we obtained were without linkage disequilibrium (kb > 10,000 and r2 < 0.001). Furthermore, all IVs reached genome-wide significance (p < 5E-08). Additionally, IVs with an F-statistic of less than 10 were considered weak instruments and were omitted from the MR analysis thus reducing the bias in the estimation of the IVs. Ultimately, we selected 583 SNPs, 127 SNPs, 144 SNPs, and 47 SNPs as IVs for ALM, hand grip strength (left), hand grip strength (right), and usual walking pace, respectively. Details of the IVs for SP-related traits are displayed in Supplementary Tables S1, S2.

As displayed in Table 2, a portion of Cochran’s Q test was used to detect heterogeneity (p < 0.05), hence we used the IVW method in the random effects model. In the IVW model, adjusted for the FDR, genetically elevated SP-related traits levels were causally associated with all OA, ALM (Odds ratio (OR) = 1.103, 95% confidence interval (CI) = 1.052–1.156, PFDR = 2.87E-04), hand grip strength (left, OR = 0.823, 95% CI = 0.712 to 0.952, PFDR = 0.020; right, OR = 0.826, 95% CI = 0.718 to 0.950, PFDR = 0.020), and usual walking pace (OR = 0.339, 95% CI = 0.204 to 0.564, PFDR = 2.38E-04). The results of the WM and MR-PRESSO were consistent with the IVW (Table 3). No evidence of directional pleiotropy in the MR-Egger intercept test was observed for all IVs (p for intercept >0.05) (Table 2). Leave-one-out test indicated that SNPs without large effect sizes skewed the estimates (Supplementary Figure S1).

TABLE 2
www.frontiersin.org

TABLE 2. Sensitivity and IVW analysis of the causal effects of sarcopenia-related traits on osteoarthritis.

TABLE 3
www.frontiersin.org

TABLE 3. Weighted median, MR Egger and MR PRESSO analysis of the causal effects of sarcopenia-related traits on osteoarthritis.

In the site-specific OA analysis, considering the results of the heterogeneity test, the IVW in the random effects model was used. In the IVW model, we observed that genetically determined levels of ALM were causally associated with knee OA (OR = 1.246, 95% CI = 1.151 to 1.350, PFDR = 9.71E-07), usual walking pace was negatively correlated with hip OA (OR = 0.474, 95% CI = 0.272 to 0.827, PFDR = 0.020) and knee OA (OR = 0.299, 95% CI = 0.136 to 0.659, PFDR = 0.011) (Table 2). Moreover, WM and MR-PRESSO reached similar causal conclusions, and the MR-Egger effect estimate was in the same direction as IVW, which was considered supportive (Table 3). Further, no directional pleiotropy was detected for our selected IVs on MR-Egger analysis (Table 2). Leave-one-out test results were consistent with the above (Supplementary Figures S2–S4).

The causal effect of OA on SP-related traits

After similar screening criteria, 25 SNPs, 8 SNPs, 28 SNPs, and 22 SNPs were obtained as IVs for all OA, hand OA, hip OA, and knee OA, respectively (Supplementary Tables S1, S3). As Table 4 illustrates, heterogeneity was observed (p < 0.05) between the IVs of all the selected OA and the traits associated with SP. Therefore, we used the random effects model of IVW. The results of our analysis showed that the onset and progression of OA could lead to the worsening of sarcopenia-related traits [ALM (IVW β = −0.258, 95% CI = −0.369 to −0.146, PFDR = 0.6.07E-06), grip strength (left, β = −0.064, 95% CI = −0.104 to −0.024, PFDR = 0.002; right, β = −0.055, 95% CI = −0.095 to −0.014, PFDR = 0.008), and usual walking pace (β = −0.104, 95% CI = −0.147 to −0.061, PFDR = 1.61E-05)]. The WM and MR-PRESSO suggested similar findings, and the beta values of the MR-Egger are in the same direction (Table 5). MR-Egger intercept test did not demonstrate directional pleiotropy for our selected IVs (Table 4). The results of the leave-one-out- test displayed that SNPs without large effect sizes resulted in biases in the estimates (Supplementary Figure S5).

TABLE 4
www.frontiersin.org

TABLE 4. Sensitivity and IVW analysis of the causal effects of osteoarthritis on sarcopenia-related traits.

TABLE 5
www.frontiersin.org

TABLE 5. Weighted median, MR Egger and MR PRESSO analysis of the causal effects of osteoarthritis on sarcopenia-related traits.

In the site-specific analysis, the results of using the IVW method in the random effects model pointed to a negative causal relationship between hand OA and grip strength (left, β = −0.102, 95% CI = −0.149 to −0.055, PFDR = 7.23E-05; right, β = -0.104, 95% CI = −0.152 to −0.057, PFDR = 7.23E-05). Hip OA (β = −0.025, 95% CI = −0.038 to −0.013, PFDR = 2.62E-04) and knee OA (β = −0.064, 95% CI = −0.087 to −0.04, PFDR = 1.37E-06), remained significantly causally and negatively correlated with usual walking pace (Table 4). In addition, we identified that only knee OA was causally related to ALM in the available site-specific cohorts (β = −0.179, 95% CI = −0.281 to −0.076, PFDR = 0.001), the results of all other MR analysis methods were similar to those of the IVW (Table 5). Heterogeneity for the associations between the selected IVs of the knee OA and ALM was not observed, as in the study above no directional pleiotropy was detected (Table 4). The results of the leave-one-out test were also consistent with the aforementioned results (Supplementary Figures S6–S8).

Discussion

We investigated the potential causality between SP-related traits and OA. Following the FDR correction, the results demonstrated a possible causal relationship between SP-related traits and all OA, the usual walking pace and hip OA, the usual walking pace and knee OA, the ALM and knee OA, respectively. In addition, we found evidence of their previous reverse causation with each other. The results also indicated a causal relationship between hand OA and hand grip strength.

Firstly, SP increases the risk of OA. In recent years, studies have demonstrated a correlation of varying intensity between SP and OA (Jin et al., 2017; Vlietstra et al., 2019; Dalle and Koppo, 2020; Godziuk et al., 2021). Decreased muscle strength is the main feature of SP and previous experiments have demonstrated that decreased muscle strength or muscle weakness is a risk factor for the development and progression of OA (Tanaka et al., 2019; Xu et al., 2020). Andrews et al. discovered that in men only, the likelihood of knee OA was reduced for each standard deviation (SD) reduction in ALM (OR per SD reduction: 0.68; 95% CI: 0.47–0.97) (Andrews et al., 2021). Similarly, another study noted significant ALM or total lean mass and increased fat mass (FM) were associated with radiographic knee OA (Azuma et al., 2017). Optimizing medial femoral size is important in clinical management to reduce the progression of OA and subsequent knee arthroplasty, which in part reflects the important impact of SP on the development of OA (Wang et al., 2012).

In addition, OA has been noted as a risk factor for increasing the incidence of SP (Kemmler et al., 2015; Dharmakulsakti et al., 2022). A previous foundational experimental study provided some mechanistic insights, indicating that knee osteoarthritis induced by anterior cruciate ligament transection promotes remodeling and atrophy in the neuromuscular junctions of the quadriceps and tibialis anterior muscles. These changes were associated with signs of inflammation and alterations in muscle gene and protein expression (Cunha et al., 2019). Besides, some studies have indicated a strong relationship between OA and low skeletal muscle mass (Berenbaum and van den Berg, 2015; Jeon et al., 2019). A significant association between OA of the knee and walking pace was observed (p < 0.001, OR:0.073) (Nakamura and Ogata, 2016). A recent study conducted a more comprehensive analysis, revealing that the OA group had statistically significantly worse SP parameters than the control group, with lower appendicular skeletal muscle mass (p = 0.041), impaired performance on the 40-m fast walk test (p = 0.020), and reduced right (p < 0.01) and left (p < 0.01) hand grip strength. The findings suggest an early onset of sarcopenia in these individuals (de Almeida et al., 2020).

According to previous studies, bone, muscle, and fat tissue are connected and interact with each other through molecules. And SP seems to have a bidirectional relationship with the maintenance or destruction of joint structures (Spanoudaki et al., 2023). SP and OA elevate the risk of mutual development, which is one of the most important reasons for their frequent coexistence (Peng and Zeng, 2022). A new concept of “sarcopenic knee OA” has also been proposed (Iijima and Aoyama, 2021). The study by Jiyong Yang et al. identified a common network of genetic interactions between KOA and SP, including 14 common differentially expressed genes, 4 hub genes, and 10 potential chemical compounds, among other important findings, which have updated the research results of the mechanism between OA and SP (Yang et al., 2023).

These previous studies support our view to some extent. Our study’s methodology thus has some advantages. Firstly, the MR method can effectively avoid the drawbacks of traditional observational research methods such as residual confounding uncertainty and reverse causality. Secondly, the IVs for SP-related traits were obtained from the existing large GWAS, and the IVs for OA were obtained from the most recent GWAS. This allowed for a more precise assessment of effect sizes than would be possible with individual-level data or findings from studies with limited sample sizes. Lastly, we performed an analysis of the specific relationships between individual SP-related traits and the different sites of OA, which led to a more comprehensive understanding of the potential link between them.

Nevertheless, our study also has some limitations. Firstly, due to the lack of data on the large GWAS concerning SP, we could only use the related traits to analyze the relationship with OA. Secondly, selecting SNPs from the different large-sample GWAS studies may increase the risk of sample overlap between exposure and outcome variables, which may bias the results. Furthermore, due to the complexity of biological systems, bidirectional MR assumes that causality happens in one direction, and feedback loops may exist between the exposure and the outcome, which might affect the accuracy of the results. Finally, since a majority of the participants included in the study were of European ancestry and were not representative of other racial groups, further magnetic resonance studies are needed to verify causality.

Conclusion

This present study suggests an obvious causality of SP on OA, with condition exhibiting site-specific effects, while evidence was also provided for the causal effect of OA on SP. It presents some evidence of reciprocal interaction between SP and OA, which may facilitate the development of novel treatment strategies for both diseases. However, the causal relationship between the two conditions still necessitates further investigation and substantiation through a multitude of studies.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author contributions

X-MJ: Writing–original draft, Writing–review and editing. T-TD: Methodology, Supervision, Validation, Writing–review and editing. HS: Resources, Visualization, Writing–review and editing. H-JS: Validation, Writing–review and editing. HQ: Validation, Writing–review and editing. G-CY: Conceptualization, Validation, Writing–review and editing. YY: Writing–review and editing. F-JL: Methodology, Supervision, Writing–review and editing. BS: Funding acquisition, Resources, Supervision, Validation, Writing–review and editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The key R&D Program of Shandong Province (2022CXGC020510) and the Academic Promotion Program of Shandong First Medical University (2019QL003) funded this study.

Acknowledgments

We would like to thank Editage (www.editage.cn) for English language editing.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1340245/full#supplementary-material

References

Andrews, J. S., Gold, L. S., Nevitt, M., Heagerty, P. J., and Cawthon, P. M. (2021). Appendicular lean mass, grip strength, and the development of knee osteoarthritis and knee pain among older adults. ACR open Rheumatol. 3 (8), 566–572. doi:10.1002/acr2.11302

PubMed Abstract | CrossRef Full Text | Google Scholar

Azuma, K., Sera, Y., Shinjo, T., Takayama, M., Shiomi, E., Momoshima, S., et al. (2017). Maintenance of the paraspinal muscles may protect against radiographic knee osteoarthritis. Open access rheumatology Res. Rev. 9, 151–158. doi:10.2147/oarrr.S130688

PubMed Abstract | CrossRef Full Text | Google Scholar

Berenbaum, F., and van den Berg, W. B. (2015). Inflammation in osteoarthritis: changing views. Osteoarthritis Cartilage. Nov 23 (11), 1823–1824. doi:10.1016/j.joca.2015.09.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Boer, C. G., Hatzikotoulas, K., Southam, L., Stefánsdóttir, L., Zhang, Y., Coutinho de Almeida, R., et al. (2021). Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations. Cell 184 (18), 4784–4818. doi:10.1016/j.cell.2021.07.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 44 (2), 512–525. doi:10.1093/ije/dyv080

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Del Greco, M. F., Minelli, C., Davey Smith, G., Sheehan, N. A., and Thompson, J. R. (2016). Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: the role of the I2 statistic. Int. J. Epidemiol. 45 (6), 1961–1974. doi:10.1093/ije/dyw220

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgess, S., Bowden, J., Fall, T., Ingelsson, E., and Thompson, S. G. (2017). Sensitivity analyses for robust causal inference from mendelian randomization analyses with multiple genetic variants. Epidemiol. Camb. Mass) 28 (1), 30–42. doi:10.1097/ede.0000000000000559

CrossRef Full Text | Google Scholar

Burgess, S., Davey Smith, G., Davies, N. M., Dudbridge, F., Gill, D., Glymour, M. M., et al. (2019). Guidelines for performing Mendelian randomization investigations. Wellcome open Res. 4, 186. doi:10.12688/wellcomeopenres.15555.2

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgess, S., and Thompson, S. G.CRP CHD Genetics Collaboration (2011). Avoiding bias from weak instruments in Mendelian randomization studies. Int. J. Epidemiol. 40 (3), 755–764. doi:10.1093/ije/dyr036

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Robinson, D. G., and Storey, J. D. (2021). The functional false discovery rate with applications to genomics. Biostat. Oxf. Engl. 22 (1), 68–81. doi:10.1093/biostatistics/kxz010

CrossRef Full Text | Google Scholar

Cho, Y., Jeong, S., Kim, H., Kang, D., Lee, J., Kang, S. B., et al. (2021). Disease-modifying therapeutic strategies in osteoarthritis: current status and future directions. Exp. Mol. Med. 53 (11), 1689–1696. doi:10.1038/s12276-021-00710-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Cruz-Jentoft, A. J., Baeyens, J. P., Bauer, J. M., Boirie, Y., Cederholm, T., Landi, F., et al. (2010). Sarcopenia: European consensus on definition and diagnosis: report of the European working group on sarcopenia in older people. Age ageing 39 (4), 412–423. doi:10.1093/ageing/afq034

PubMed Abstract | CrossRef Full Text | Google Scholar

Cruz-Jentoft, A. J., Bahat, G., Bauer, J., Boirie, Y., Bruyère, O., Cederholm, T., et al. (2019). Sarcopenia: revised European consensus on definition and diagnosis. Age ageing 48 (1), 16–31. doi:10.1093/ageing/afy169

PubMed Abstract | CrossRef Full Text | Google Scholar

Cruz-Jentoft, A. J., and Sayer, A. A. (2019). Sarcopenia. Lancet (London, Engl. 393 (10191), 2636–2646. doi:10.1016/s0140-6736(19)31138-9

CrossRef Full Text | Google Scholar

Cunha, J. E., Barbosa, G. M., Castro, P., Luiz, B. L. F., Silva, A. C. A., Russo, T. L., et al. (2019). Knee osteoarthritis induces atrophy and neuromuscular junction remodeling in the quadriceps and tibialis anterior muscles of rats. Sci. Rep. 9 (1), 6366. doi:10.1038/s41598-019-42546-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Dalle, S., and Koppo, K. (2020). Is inflammatory signaling involved in disease-related muscle wasting? Evidence from osteoarthritis, chronic obstructive pulmonary disease and type II diabetes. Exp. Gerontol. 137, 110964. doi:10.1016/j.exger.2020.110964

PubMed Abstract | CrossRef Full Text | Google Scholar

de Almeida, A. C., Aily, J., Pedroso, M., Bergamaschi, J. R., de Carvalho Silva, G., and Mattiello, S. (2020). Middle-aged adults with knee osteoarthritis present worse sarcopenia parameters, symptoms and functional performance. Osteoarthr. Cartil. 28, S399–S400. doi:10.1016/j.joca.2020.02.623

CrossRef Full Text | Google Scholar

Dharmakulsakti, P., Roopsawang, I., and Aree-Ue, S. (2022). Sarcopenia among older adults with knee osteoarthritis: a cross-sectional study of prevalence and its associated factors. Pac. Rim Int. J. Nurs. Res. 26 (1), 121–134.

Google Scholar

Godziuk, K., Prado, C. M., Woodhouse, L. J., and Forhan, M. (2021). Associations between self-reported weight history and sarcopenic obesity in adults with knee osteoarthritis. Obes. (Silver Spring, Md) 29 (2), 302–307. doi:10.1002/oby.23074

CrossRef Full Text | Google Scholar

Hemani, G., Bowden, J., and Davey Smith, G. (2018b). Evaluating the potential role of pleiotropy in Mendelian randomization studies. Hum. Mol. Genet. 27 (R2), R195–R208. doi:10.1093/hmg/ddy163

PubMed Abstract | CrossRef Full Text | Google Scholar

Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., Baird, D., et al. (2018a). The MR-Base platform supports systematic causal inference across the human phenome. eLife 7, e34408. doi:10.7554/eLife.34408

PubMed Abstract | CrossRef Full Text | Google Scholar

Hida, T., Harada, A., Imagama, S., and Ishiguro, N. (2014). Managing sarcopenia and its related-fractures to improve quality of life in geriatric populations. Aging Dis. 5 (4), 226–237. doi:10.14336/ad.2014.0500226

PubMed Abstract | CrossRef Full Text | Google Scholar

Holmes, M. V., Ala-Korpela, M., and Smith, G. D. (2017). Mendelian randomization in cardiometabolic disease: challenges in evaluating causality. Nat. Rev. Cardiol. 14 (10), 577–590. doi:10.1038/nrcardio.2017.78

PubMed Abstract | CrossRef Full Text | Google Scholar

Iijima, H., and Aoyama, T. (2021). Increased recurrent falls experience in older adults with coexisting of sarcopenia and knee osteoarthritis: a cross-sectional study. BMC Geriatr. 21 (1), 698. doi:10.1186/s12877-021-02654-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeon, H., Lee, S. U., Lim, J. Y., Chung, S. G., Lee, S. J., and Lee, S. Y. (2019). Low skeletal muscle mass and radiographic osteoarthritis in knee, hip, and lumbar spine: a cross-sectional study. Aging Clin. Exp. Res. 31 (11), 1557–1562. doi:10.1007/s40520-018-1108-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, W. S., Choi, E. J., Lee, S. Y., Bae, E. J., Lee, T. H., and Park, J. (2017). Relationships among obesity, sarcopenia, and osteoarthritis in the elderly. J. Obes. metabolic syndrome 26 (1), 36–44. doi:10.7570/jomes.2017.26.1.36

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, T. L., Esa, M. S., Li, K. H. C., Krishnan, S. R. G., Elgallab, G. M., Pearce, M. S., et al. (2021). Osteoporosis, fracture, osteoarthritis and sarcopenia: a systematic review of circulating microRNA association. Bone 152, 116068. doi:10.1016/j.bone.2021.116068

PubMed Abstract | CrossRef Full Text | Google Scholar

Kemmler, W., Teschler, M., Goisser, S., Bebenek, M., von Stengel, S., Bollheimer, L. C., et al. (2015). Prevalence of sarcopenia in Germany and the corresponding effect of osteoarthritis in females 70 years and older living in the community: results of the FORMoSA study. Clin. interventions aging 10, 1565–1573. doi:10.2147/cia.S89585

PubMed Abstract | CrossRef Full Text | Google Scholar

Mezian, K., Angerová, Y., Kara, M., and Özçakar, L. (2021). Obesity paradox in sarcopenia and knee osteoarthritis: comment on the article by Andrews et al. ACR open Rheumatol. 3 (11), 812–813. doi:10.1002/acr2.11330

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, R., Elsworth, B., Mitchell, R., Raistrick, C. A., Paternoster, L., Hemani, G., et al. (2019). MRC IEU UK Biobank GWAS pipeline version 2. Bristol, UK: University of Bristol.

Google Scholar

Nakamura, K., and Ogata, T. (2016). Locomotive syndrome: definition and management. Clin. Rev. bone mineral metabolism 14 (2), 56–67. doi:10.1007/s12018-016-9208-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Pegreffi, F., Balestra, A., De Lucia, O., Smith, L., Barbagallo, M., and Veronese, N. (2023). Prevalence of sarcopenia in knee osteoarthritis: a systematic review and meta-analysis. J. Clin. Med. 12 (4), 1532. doi:10.3390/jcm12041532

PubMed Abstract | CrossRef Full Text | Google Scholar

Pei, Y. F., Liu, Y. Z., Yang, X. L., Zhang, H., Feng, G. J., Wei, X. T., et al. (2020). The genetic architecture of appendicular lean mass characterized by association analysis in the UK Biobank study. Commun. Biol. 3 (1), 608. doi:10.1038/s42003-020-01334-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, H., and Zeng, Y. (2022). Research progress on the correlation between sarcopenia and osteoarthritis. Zhongguo xiu fu chong jian wai ke za zhi = Zhongguo xiufu chongjian waike zazhi = Chin. J. reparative Reconstr. Surg. 36 (12), 1549–1557. doi:10.7507/1002-1892.202209015

CrossRef Full Text | Google Scholar

Petermann-Rocha, F., Balntzi, V., Gray, S. R., Lara, J., Ho, F. K., Pell, J. P., et al. (2022). Global prevalence of sarcopenia and severe sarcopenia: a systematic review and meta-analysis. J. cachexia, sarcopenia muscle 13 (1), 86–99. doi:10.1002/jcsm.12783

PubMed Abstract | CrossRef Full Text | Google Scholar

Richmond, R. C., and Davey Smith, G. (2022). Mendelian randomization: concepts and scope. Cold Spring Harb. Perspect. Med. 12 (1), a040501. doi:10.1101/cshperspect.a040501

PubMed Abstract | CrossRef Full Text | Google Scholar

Spanoudaki, M., Giaginis, C., Mentzelou, M., Bisbinas, A., Solovos, E., Papadopoulos, K., et al. (2023). Sarcopenia and sarcopenic obesity and osteoarthritis: a discussion among muscles, fat, bones, and aging. Life (Basel, Switz. 13 (6), 1242. doi:10.3390/life13061242

CrossRef Full Text | Google Scholar

Sudlow, C., Gallacher, J., Allen, N., Beral, V., Burton, P., Danesh, J., et al. (2015). UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 12 (3), e1001779. doi:10.1371/journal.pmed.1001779

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka, R., Hirohama, K., and Ozawa, J. (2019). Can muscle weakness and disability influence the relationship between pain catastrophizing and pain worsening in patients with knee osteoarthritis? A cross-sectional study. Braz. J. Phys. Ther. 23 (3), 266–272. doi:10.1016/j.bjpt.2018.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Therakomen, V., Petchlorlian, A., and Lakananurak, N. (2020). Prevalence and risk factors of primary sarcopenia in community-dwelling outpatient elderly: a cross-sectional study. Sci. Rep. 10 (1), 19551. doi:10.1038/s41598-020-75250-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Tzartza, C. L., Karapalis, N., Voulgaridou, G., Zidrou, C., Beletsiotis, A., Chatziprodromidou, I. P., et al. (2023). Comparative analysis on the effect of sarcopenia in patients with knee osteoarthritis before and after total knee arthroplasty. Dis. (Basel, Switz. 11 (1), 36. doi:10.3390/diseases11010036

CrossRef Full Text | Google Scholar

Verbanck, M., Chen, C. Y., Neale, B., and Do, R. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50 (5), 693–698. doi:10.1038/s41588-018-0099-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Vlietstra, L., Stebbings, S., Meredith-Jones, K., Abbott, J. H., Treharne, G. J., and Waters, D. L. (2019). Sarcopenia in osteoarthritis and rheumatoid arthritis: the association with self-reported fatigue, physical function and obesity. PloS one 14 (6), e0217462. doi:10.1371/journal.pone.0217462

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Wluka, A. E., Berry, P. A., Siew, T., Teichtahl, A. J., Urquhart, D. M., et al. (2012). Increase in vastus medialis cross-sectional area is associated with reduced pain, cartilage loss, and joint replacement risk in knee osteoarthritis. Arthritis rheumatism 64 (12), 3917–3925. doi:10.1002/art.34681

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, C., and Xiao, G. (2022). Advances in osteoarthritis research in 2021 and beyond. J. Orthop. Transl. 32, A1–A2. doi:10.1016/j.jot.2022.02.011

CrossRef Full Text | Google Scholar

Xu, J., She, G., Gui, T., Hou, H., Li, J., Chen, Y., et al. (2020). Knee muscle atrophy is a risk factor for development of knee osteoarthritis in a rat model. J. Orthop. Transl. 22, 67–72. doi:10.1016/j.jot.2019.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Jiang, T., Xu, G., Wang, S., and Liu, W. (2023). Exploring molecular mechanisms underlying the pathophysiological association between knee osteoarthritis and sarcopenia. Osteoporos. sarcopenia 9 (3), 99–111. doi:10.1016/j.afos.2023.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Yavorska, O. O., and Burgess, S. (2017). MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data. Int. J. Epidemiol. 46 (6), 1734–1739. doi:10.1093/ije/dyx034

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, P., Jiang, Y., Fang, X., Wang, D., Li, Y., Chen, M., et al. (2023). Cell-based therapies for degenerative musculoskeletal diseases. Adv. Sci. (Weinheim, Baden-Wurttemberg, Ger. 10, e2207050. doi:10.1002/advs.202207050

CrossRef Full Text | Google Scholar

Keywords: sarcopenia, osteoarthritis, Mendelian randomization, degenerative musculoskeletal diseases, causal relationship

Citation: Jia X-M, Deng T-T, Su H, Shi H-J, Qin H, Yu G-C, Yin Y, Liu F-J and Shi B (2024) Genetic causality and site-specific relationship between sarcopenia and osteoarthritis: a bidirectional Mendelian randomization study. Front. Genet. 14:1340245. doi: 10.3389/fgene.2023.1340245

Received: 17 November 2023; Accepted: 18 December 2023;
Published: 08 January 2024.

Edited by:

Ziheng Wang, University of Macau, China

Reviewed by:

Zhengxiao Wei, Public Health Clinical Center of Chengdu, China
Lin Chen, Fujian Medical University, China

Copyright © 2024 Jia, Deng, Su, Shi, Qin, Yu, Yin, Liu and Shi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Gong-Chang Yu, eXVnb25nY2hhbmdAc2RmbXUuZWR1LmNu; Ying Yin, NTYzMjk4MDk4QHFxLmNvbQ==; Fan-Jie Liu, NTAxODA0NzlAcXEuY29t; Bin Shi, YnNoaUBzZGZtdS5lZHUuY24=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.