- 1Honghui Hospital, Xi’an Jiaotong University, Xi’an, China
- 2Department of Orthopedic Center, The First Hospital of Jilin University, Changchun, China
Background: Despite extensive research, the association between serum α-klotho levels and osteoarthritis (OA) remains unclear, predominantly relying on findings from OA mouse models. This study used data from the National Health and Nutrition Examination Survey (NHANES) to conduct a cross-sectional study examining the relationship between α-klotho and human OA. In addition, we used Mendelian randomization (MR) to genetically infer a causal relationship between serum α-klotho and the three OA subtypes.
Method: A cohort of 12,037 subjects from NHANES (2007-2016) was analyzed. Multivariate logistic regression was utilized to examine the association between α-klotho concentration and OA, alongside subgroup analysis and interaction tests. Additionally, a two-sample bi-directional MR analysis was conducted to evaluate the relationship between serum α-klotho and three OA subtypes, including all OA, hip OA, and knee OA, employing the inverse variance weighting (IVW) method as the primary approach.
Results: Following adjustment for covariates, a nonlinear negative correlation between serum α-klotho and OA was observed (OR=0.77; 95% CI, 0.68-0.88, p < 0.0001). The IVW method revealed that higher serum α-klotho levels were associated with decreased susceptibility to hip OA (OR = 0.92, 95% CI: 0.87–0.98, P = 9.64×10-3). However, MR analysis did not establish a causal relationship between serum α-klotho and OA or knee OA. Inverse MR also indicated that the three subtypes of OA do not causally affect serum α-klotho concentrations.
Conclusions: In cross-sectional studies, α-klotho showed a nonlinear negative correlation with OA. MR analysis of outcomes was not identical to cross-sectional studies.
Introduction
Osteoarthritis (OA) stands as a significant global contributor to pain and disability, affecting an estimated 10% of men and 18% of women over 60 years of age (1), imposing substantial societal and individual economic burdens (1, 2). There are various risk factors for OA such as age, female gender, genetic predisposition, comorbidities, obesity, and joint injuries (1, 3–8). OA was once thought to be a purely mechanical disease of cartilage degeneration, but it is now known to be a complex disease affecting the entire joint: cartilage, subchondral bone, synovium, and infrapatellar fat pad may all play a role in the pathogenesis of the disease, and it may be associated with systemic inflammation (1, 9–12). Biological pathways within joints are mechanosensitive (13). Prolonged mechanical stress and age-related wear and tear typically lead to cartilage degradation (14). But a breakthrough in non-surgical OA treatment efficacy remains elusive, underscoring the growing importance of investigating OA progression mechanisms (15).
The intricate molecular mechanisms underpinning OA development involve proteomics as a pivotal player (16, 17). Encoded by the Klotho gene, α-Klotho was initially characterized as an anti-aging protein (18). Its functions include inhibiting insulin and Wnt signaling pathways, mitigating oxidative stress, and regulating phosphatase activity and calcium uptake (18). Soluble α-klotho has paracrine and endocrine functions in different organs, and reduced levels of α-klotho in animal models of aging correlate with the initial manifestations of aging (18). Recent research has unveiled the role of α-Klotho’s epigenetic regulation in modulating the extracellular matrix associated with aging, thereby influencing chondrocyte physiology (19). Gu et al. noted a significant elevation of Klotho expression in normal mouse cartilage compared to OA models, alongside heightened activity in the Wnt/β-catenin pathway and its downstream targets in the latter (20). Similarly, Klotho gene transfer into mouse joints impeded cartilage degradation in OA models by suppressing the NOS2 and ZIP8-MMP13 catabolic axes, highlighting its crucial role in preserving cartilage tissue homeostasis (21). However, many existing mouse models of OA exhibit limitations; notably, they predominantly utilize post-traumatic young male mice, diverging from the predominantly age-related nature of human OA (19, 22). Furthermore, there remains uncertainty regarding the translatability of these preclinical models to human disease. Additionally, observational studies encounter inherent challenges in controlling for confounding factors and addressing bias from reverse causation, posing obstacles to elucidating the relationship between α-klotho and OA.
Mendelian randomization (MR) represents an innovative approach in genetic epidemiology, leveraging genome-wide association studies (GWAS) datasets and employing genetic variation as instrumental variables (IVs) to discern the causal effects of exposure on disease outcomes (23). By virtue of its random allocation and stability across germline genotypes, MR mitigates biases stemming from confounding and reverse causation (23, 24). In our current investigation, we initially utilized data sourced from the National Health and Nutrition Examination Survey (NHANES) to scrutinize the association between serum α-klotho concentrations and OA, adjusting for various potential confounders. Subsequently, we delved deeper into elucidating the causal link between serum α-klotho concentration and OA through MR analysis. From those, we wished to validate the effect of changes in human serum α-klotho concentration on OA susceptibility through a combination of observational and MR analyses.
Method
All cross-sectional and GWAS data used in this study were appropriately consented by the participants and ethical approval was obtained in accordance with national legal and institutional requirements.
Study design
In this study, we first examined the relationship between α-klotho and OA using data from the NHANES. We then analyzed summary statistics from GWAS using MR analysis and assessed the bidirectional causal associations between α-klotho and three osteoarthritis subtypes.
Cross-sectional study
Description of data sources
The NHANES, conducted by the Centers for Disease Control and Prevention (CDC), serves as a comprehensive survey of the noninstitutionalized U.S. civilian population, ensuring national representation. The NHANES study protocols received approval from the National Center for Health Statistics Institutional Review Board (IRB), as documented at https://www.cdc.gov/nchs/nhanes/irba98.htm. Prior to participation, all individuals provided both verbal and written consent for their involvement in future research endeavors. Detailed information regarding interviews, medical examinations, and laboratory tests conducted during NHANES is accessible at www.cdc.gov/nchs/nhanes/. Our investigation utilized NHANES data spanning from 2007 to 2016, encompassing a cohort of 50,588 subjects. Following the exclusion of individuals with missing data on serum α-Klotho levels, OA, and covariates, our analysis involved a total of 12,037 participants (Figure 1).
OA defined
Participants with OA were identified based on a documented diagnosis by a healthcare professional. Individuals aged 20 years or older were queried with two specific questions: “Have you ever received a diagnosis of arthritis from a doctor or other healthcare professional?” and “If so, what type of arthritis were you diagnosed with?” Those respondents affirming both questions by indicating a diagnosis of osteoarthritis were classified into the OA group, while those providing different responses were categorized into the non-OA group(www.cdc.gov/nchs/nhanes/).
Serum α-klotho defined
Serum α-klotho served as the primary variable under observation in this study. Soluble α-klotho, a multifaceted protein encompassing paracrine and endocrine functions across various organs, was the focus of investigation. Diminished levels of α-klotho in animal aging models correlate with initial manifestations of aging. Klotho concentration was quantified in NHANES frozen samples using a commercially available enzyme-linked immunosorbent assay (ELISA) kit manufactured by IBL International, Japan (25). The samples were received on dry ice and were scanned in the laboratory receiving area where the data were compared to data received on an electronic manifest and then entered into the laboratory information system. Samples were analyzed in duplicate and averaged to calculate the final value. Results are automatically transferred from the instrument to the laboratory’s Oracle management system and are evaluated by the regional supervisor. Samples with more than 10% duplicate results are flagged for repeat analysis. If the QC sample value is not within 2SD of the specified value, the entire analytical run is rejected and the sample is repeated.
Covariates
Covariates encompassed demographic data such as age, gender, and race. Body mass index (BMI) was assessed by healthcare professionals at the Mobile Examination Center (MEC). Diabetes and hypertension were delineated based on documented physician diagnoses or prescribed medications. Cotinine concentrations served as indicators of both active smoking and exposure to secondhand smoke. Individuals classified as alcohol drinkers were those who reported consuming more than 12 alcoholic beverages within the past year.
Statistical
The statistical analyses were carried out employing Empower Stats (version 2.0). To compare the demographic characteristics of subjects in the OA and non-OA groups, t-tests and chi-square tests were utilized. Beta values and their corresponding 95% confidence intervals were computed using multivariate logistic regression analysis to assess the relationship between α-klotho concentration and OA.
Model 1 was unadjusted, Model 2 incorporated adjustments for age, race, and gender, while Model 3 included additional adjustments for age, gender, race, education level, Diabetes, Hypertension, Alcohol drinking, Cotinine, and BMI. Simultaneous smoothed curve fits were performed after adjusting these variables. Furthermore, subgroup analyses and interaction tests were also undertaken to investigate the relationship between the α-klotho concentration and OA. We analyzed subgroups separately for BMI, gender, and age. Subjects were divided into three groups according to BMI <25, 25-30, and >30; two groups according to gender; and two groups according to age 40-59 and 60-79 for subgroup analyses, respectively, in order to further exclude confounding factors. The threshold for statistical significance was defined as p < 0.05.
Mendelian randomization
Figure 2 presents a flowchart elucidating the bidirectional causal relationship between α-klotho and OA. For an MR study to yield valid results, three key assumptions must be met (1): genetic variants must exhibit correlation with the exposure under investigation (2); these variants should not be linked with confounding factors; and (3) they should solely influence the outcome through the exposure. To conduct bidirectional two-sample MR analyses, we constructed distinct IVs comprising independent single nucleotide polymorphisms (SNPs) associated with α-klotho and the three subtypes of OA, including all OA, knee OA and hip OA.
Figure 2. Schematic representation of Mendelian randomization analysis. α-klotho and osteoarthritis as exposure and outcome, respectively, and the instrumental variables must meet three major assumptions. IVs, instrument variables; SNPs, single-nucleotide polymorphisms.
Data source
A summary-level GWAS dataset including data for α-klotho was extracted from a comprehensive meta-analysis of circulating α-klotho (26), comprising 4,675 European-ancestry participants. The GWAS data for OA consisted of 39,515 cases and 445,083 controls, and all subjects were from UK ancestry (27). OA is defined as a non-inflammatory degenerative joint disease occurring primarily in older adults and characterized by degeneration of articular cartilage, marginal osteophytes and synovial changes. It is associated with pain and stiffness, especially after prolonged activity (27). To increase the generalizability of the study, we additionally added two major OA subtypes, knee OA and hip OA. Pooled GWAS data for these two phenotypes were obtained from UK biobank, containing 403,124 and 393,873 UK ancestry individual data, respectively (28).
Instrumental variables extracting
We conducted SNP screening based on genome-wide significance (p < 5 × 10−8) regarding their correlation with the exposure, followed by further refinement through the elimination of linkage disequilibrium (LD) (clumping r2 = 0.001 and kb = 10,000) to identify instrumental variables (IVs) (29). Subsequently, prior to each MR analysis, IVs associated with the outcome were excluded. Additionally, we utilized PhenoScanner (www.phenoscanner.medschl.cam.ac.uk), an extensive platform providing detailed information on genotype-phenotype relationships, to identify and eliminate potential confounders by screening out secondary phenotypes of SNPs.
Statistical
The inverse-variance weighted (IVW) method (30) served as our primary MR analytical approach. To determine the causal association between OA and α-klotho, we conducted random-effects meta-analysis, combining the ratios of SNP-exposure to SNP-outcome. Additionally, to enhance the precision of IVW estimates, we employed MR-Egger and weighted median methods. For sensitivity analysis, we utilized Cochran’s Q test to assess heterogeneity. Furthermore, we employed the MR-Egger intercept test, “leave-one-out analyses,” funnel plot, and MR-PRESSO method to evaluate horizontal and directional pleiotropy.
All statistical analyses were performed using the TwoSampleMR packages in R software (Version 4.3.0), and p-value < 0.05 was seen as significantly association.
Results
Cross-sectional study
In the cross-sectional study, a total of 12,036 participants were included, of which 48.85% were males and 51.15% were females, with a mean age of 57.86 ± 10.81 years. Among all the participants, 1,686 (14%) had OA. Age, gender, race, education, cotinine, BMI, diabetes, hypertension, and serum α-klotho concentrations were statistically significant between OA and non-OA individuals (p < 0.05). OA patients, compared to those without, were more likely to be female, Non-Hispanic White, be nonsmoker, exhibit advanced age, higher BMI, in our study (p < 0.05). Table 1 presents the clinical and biochemical characteristics of the subjects based on whether they were OA patients.
Our study reveals an association between elevated serum α-klotho concentrations and susceptibility to OA. In the unadjusted model, a 1ug/ml increase in α-klotho concentration was associated with a 34% reduction in the risk of OA development (OR=0.66; 95% CI, 0.55-0.80, p < 0.0001). Similarly, in the minimally adjusted model, each 1ug/ml increase in α-klotho concentration was associated with a 24% reduction in susceptibility to OA (OR=0.76; 95% CI, 0.64-0.92, p = 0.0047). To further verify the robustness of the conclusions, we triple α-klotho and perform a sensitivity analysis. Individuals in the highest tertile of α-klotho concentration had a statistically significant 23% lower risk of developing OA compared to individuals in the lowest tertile of α-klotho concentration (OR=0.77; 95% CI, 0.68-0.88, p < 0.0001). In addition, participants in the middle third of the α-klotho concentration group also exhibited a lower risk of OA compared to the lowest third (OR=0.83; 95% CI, 0.73-0.94, p=0.0034) (Table 2). In addition, an attempt was made to fit the nonlinear relationship between α-klotho concentration and OA into a smooth curve, which showed a nonlinear negative correlation between the two variables (Figure 3).
Figure 3. Smooth curve fittings. Association between α-klotho and osteoarthritis, Gender, Age, Race, Education level, Diabetes, Hypertension, Alcohol drinking, Cotinine, and BMI were adjusted; The solid and dotted lines represent the estimated values and their corresponding 95% CIs, respectively.
We analyzed subgroups separately for BMI, gender, and age. As shown in Table 3, there was no significant correlation between the three subgroups of BMI and the interaction test between the two gender subgroups (P for interaction>0.05). This means that there is no evidence that the association between α-klotho and OA is related to BMI and gender. Notably, when we divided the age into two subgroups, 40-59 and 60-79, for the interaction test, the results suggested that the association between α-klotho and OA might have different effects in different age groups (p=0.031).
Mendelian randomization
Ultimately, we identified 5 IVs for genetic prediction of circulating α-klotho. Additionally, 7 index SNPs were selected for genetic prediction of overall OA, with 10 index SNPs designated for knee OA and 27 index SNPs for hip OA, representing the two OA subtypes. Detailed IVs information is provided in the Supplementary Table S1–S4. Subsequently, we conducted a bidirectional two-sample MR study to elucidate the causal relationship between circulating α-klotho and OA.
In the IVW model, α-klotho demonstrated a significant reduction in the risk of hip OA (OR = 0.92, 95% CI: 0.87–0.98, P = 9.64×10-3) (Figure 4). The Cochran’s Q-derived P-value of 0.56 suggested no evidence of heterogeneity, while the MR-Egger intercept-derived P-value of 0.66 indicated no discernible pleiotropy. Conversely, no causal association between α-klotho and overall OA or knee OA was observed. Furthermore, in the reverse MR analysis, no significant causal effect of the three OA subtypes on circulating α-klotho was evident. The results of all sensitivity analyses are presented in the Supplementary Figure S1.
Figure 4. Forest plot for MR analyses. OR, odds ratio; CI, confidence interval; IVW, Inverse-variance weighted; OA, osteoarthritis.
Discussion
This study provides a comprehensive analysis of the relationship between human serum α-klotho protein and OA using data from a large observational study and an extensive genetic dataset. Investigations across diverse clinical populations have indicated a nonlinear negative correlation between serum α-klotho levels and OA, suggesting a potential protective role for α-klotho in cartilage health. Notably, subgroup analyses conducted across different age groups revealed varying associations between α-klotho concentration and OA, may indicating age-dependent effects. Xu et al. performed similar work and came to conclusions similar to ours (31). In our study, we further utilized MR analysis to stratify subgroups of OA on top of conducting observational studies, and investigated the relationship between klotho and susceptibility to the three OA subtypes at the genetic level. MR studies have further highlighted a significant protective effect of increased serum α-klotho concentration specifically in individuals with the hip OA subtype. However, no significant associations were observed between serum α-klotho levels and overall OA or knee OA. Additionally, reverse MR analysis did not suggest a causal effect of the three OA subtypes on serum α-klotho levels. These findings provide valuable insights into the mechanisms underlying OA and aging processes, shedding light on the potential role of serum α-klotho in these contexts.
The Klotho (KL) gene, initially identified as a potential senescence suppressor gene, has sparked significant interest due to its impact on lifespan regulation (32–34). Mouse studies have shown that overexpression of the KL gene extends lifespan, while mutations in this gene shorten it, thereby deepening our understanding of the aging process (35, 36). The human KL gene encodes the multifunctional α-Klotho protein, which plays a crucial role in regulating phosphate, calcium, and vitamin D metabolism (18). Although research on α-klotho and OA is relatively limited, existing findings consistently suggest its relevance. Notably, dating back to 2007, Zhang et al. examined the association between α-klotho gene expression and OA susceptibility in a female Caucasian population, observing an association between a klotho gene variant and hand OA susceptibility (37). In recent years, numerous mouse models of OA have reinforced the notion that α-klotho serves as a protective factor against OA progression (20, 21). For instance, Gu et al. demonstrated that Klotho acts as an antagonist to endogenous Wnt/β-catenin activity, thereby mitigating cartilage damage in OA (20). Furthermore, studies in rat models have shown that treatment with α-klotho and sTGFβR2 counteracts the development of OA phenotype, offering potential avenues for future cartilage repair drug development (38). Additionally, investigations by Iijima examining OA-related changes in mice across lifespan and gender found that age-related alterations in extracellular matrix properties initiate pathogenic mechanotransduction signals, leading to α-klotho promoter methylation and compromised cellular health, with broader implications for aging research (19). Despite the high degree of sequence similarity (98%) between human and mouse α-Klotho proteins (18), caution must be exercised in directly applying conclusions from animal models to humans. However, our study, based on data from a large-scale human observational study, reaffirms previous findings suggesting that higher serum α-klotho concentrations may reduce susceptibility to OA. This underscores a potential avenue for future OA prevention and non-surgical cartilage treatment strategies.
In contrast to the cross-sectional findings, MR analysis only indicated a protective effect of serum α-klotho causality on the hip OA subtype. The absence of heterogeneity and pleiotropy in the sensitivity analyses, including MR-Egger, Cochran’s Q test, and leave-one-out sensitivity, underscores the reliability of these positive findings. Furthermore, to mitigate the possibility of reverse causation, we conducted bidirectional MR analysis (39). However, MR results did not suggest a causal association between serum α-klotho concentrations and the other two OA subtypes. This disparity from observational studies may arise from the fact that participants in this cross-sectional study were exclusively middle-aged and older adults aged over 50 years, whereas MR analyses estimated exposure effects over a lifetime. It remains unclear whether an increase in serum α-klotho can exclusively reduce susceptibility to OA during specific periods of time, a bias observed in previous klotho studies (40). However, it is imperative to note that MR analysis provides a higher level of evidence. Consequently, we cannot infer that elevated serum α-klotho concentrations would significantly protect against all OA subtypes.
Our study boasts several strengths. Firstly, our cross-sectional investigation benefits from a sizable sample size and employs a diverse array of statistically robust methodologies. Secondly, we conducted a comprehensive two-sample MR analysis to assess the causal relationship between serum α-klotho and various subtypes of OA phenotypes. In contrast to observational studies, MR analysis offers the advantage of mitigating potential confounders and reverse causality. Thirdly, while prior research predominantly relied on mouse OA models, our study utilized human data, thus enhancing the reliability of our conclusions. However, our study also exhibits several limitations. Firstly, our MR analysis only considered a linear relationship between exposure and outcome, thereby failing to account for potential nonlinear relationships. Secondly, due to limitations in the available GWAS data, we were unable to conduct stratified MR analysis by sex, age, OA grade, classification and other factors, complicating the reconciliation of discrepancies with cross-sectional study findings. Third, whether subjects had OA in our observational study came from self-report rather than clinical diagnosis, and thus may have been partially biased. Fourth, our study analyses used only cross-sectional data from the U.S. population and GWAS data from individuals of European ancestry, which limits the generalizability of our findings to all population groups. Future analyses should incorporate data from different ethnic backgrounds to improve the applicability of our findings to different populations. Finally, our study focused solely on investigating the association between serum α-klotho and OA susceptibility, without delving into the specific underlying mechanisms. Further basic research and animal experiments are warranted to elucidate these mechanisms in greater detail.
Conclusion
Our cross-sectional study found an inverse nonlinear relationship between serum α-klotho and OA. MR analyses revealed that the conclusions of the cross-sectional study also applied to serum α-klotho and hip OA subtypes in European pedigrees. However, OA and knee OA subtypes were not found to be causally related to klotho. It is possible that klotho may not have the same effect on OA in different age strata, which may also account for the discrepancy between the findings of MR analyses and cross-sectional 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.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
ZhaL: Writing – original draft, Writing – review & editing. ZhoL: Writing – review & editing. QC: Writing – original draft. XN: Writing – review & editing. YC: Writing – original draft. BD: Writing – original draft. TR: Writing – original draft. YX: Writing – review & editing. TM: Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by Xi'an Science and Technology Program(2022JH-YBYJ-0432) and Shaanxi Provincial Natural Science Basic Research Program (S2024-JC-YB-0174).
Acknowledgments
We gratefully acknowledge the NHANES, the GWAS Catalog, and the UK BioBank for generously providing access to their datasets. Our heartfelt thanks go out to the dedicated researchers involved in the contributing studies, whose invaluable work has laid the foundation for our research. We also extend our sincere appreciation to the participants who willingly contributed their time and efforts to these important investigations.
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/fendo.2024.1450472/full#supplementary-material
Abbreviations
OA, Osteoarthritis; NHANES, National Health and Nutrition Examination Survey; MR, Mendelian randomization; IVW, inverse variance weighting; GWAS, genome-wide association studies; IVs, instrumental variables; CDC, Centers for Disease Control and Prevention; IRB, Institutional Review Board; BMI, Body mass index; MEC, Mobile Examination Center; LD, linkage disequilibrium; KL, Klotho; OR, odds ratio; CI, confidence interval.
References
1. Glyn-Jones S, Palmer AJ, Agricola R, Price AJ, Vincent TL, Weinans H, et al. Osteoarthritis. Lancet. (2015) 386:376–87. doi: 10.1016/S0140-6736(14)60802-3
2. Safiri S, Kolahi AA, Smith E, Hill C, Bettampadi D, Mansournia MA, et al. Global, regional and national burden of osteoarthritis 1990-2017: a systematic analysis of the Global Burden of Disease Study 2017. Ann Rheum Dis. (2020) 79:819–28. doi: 10.1136/annrheumdis-2019-216515
3. Lawrence RC, Felson DT, Helmick CG, Arnold LM, Choi H, Deyo RA, et al. Estimates of the prevalence of arthritis and other rheumatic conditions in the United States. Part II Arthritis Rheum. (2008) 58:26–35. doi: 10.1002/art.23176
4. Muthuri SG, McWilliams DF, Doherty M, Zhang W. History of knee injuries and knee osteoarthritis: a meta-analysis of observational studies. Osteoarthritis Cartilage. (2011) 19:1286–93. doi: 10.1016/j.joca.2011.07.015
5. Conde J, Scotece M, Gomez R, Lopez V, Gomez-Reino JJ, Gualillo O. Adipokines and osteoarthritis: novel molecules involved in the pathogenesis and progression of disease. Arthritis. (2011) 2011:203901. doi: 10.1155/2011/203901
6. Zhang Y, Niu J, Felson DT, Choi HK, Nevitt M, Neogi T. Methodologic challenges in studying risk factors for progression of knee osteoarthritis. Arthritis Care Res (Hoboken). (2010) 62:1527–32. doi: 10.1002/acr.v62:11
7. arc OC, arc OC, Zeggini E, Panoutsopoulou K, Southam L, Rayner NW, et al. Identification of new susceptibility loci for osteoarthritis (arcOGEN): a genome-wide association study. Lancet. (2012) 380:815–23. doi: 10.1016/S0140-6736(12)60681-3
8. Fontanella CG, Belluzzi E, Pozzuoli A, Scioni M, Olivotto E, Reale D, et al. Exploring anatomo-morphometric characteristics of infrapatellar, suprapatellar fat pad, and knee ligaments in osteoarthritis compared to post-traumatic lesions. Biomedicines. (2022) 10:1639. doi: 10.3390/biomedicines10061369
9. Fontanella CG, Belluzzi E, Pozzuoli A, Favero M, Ruggieri P, Macchi V, et al. Mechanical behavior of infrapatellar fat pad of patients affected by osteoarthritis. J Biomech. (2022) 131:110931. doi: 10.1016/j.jbiomech.2021.110931
10. Burleigh A, Chanalaris A, Gardiner MD, Driscoll C, Boruc O, Saklatvala J, et al. Joint immobilization prevents murine osteoarthritis and reveals the highly mechanosensitive nature of protease expression in vivo. Arthritis Rheum. (2012) 64:2278–88. doi: 10.1002/art.34420
11. Weinans H, Siebelt M, Agricola R, Botter SM, Piscaer TM, Waarsing JH. Pathophysiology of peri-articular bone changes in osteoarthritis. Bone. (2012) 51:190–6. doi: 10.1016/j.bone.2012.02.002
12. Scanzello CR, Goldring SR. The role of synovitis in osteoarthritis pathogenesis. Bone. (2012) 51:249–57. doi: 10.1016/j.bone.2012.02.012
13. Vincent TL. Targeting mechanotransduction pathways in osteoarthritis: a focus on the pericellular matrix. Curr Opin Pharmacol. (2013) 13:449–54. doi: 10.1016/j.coph.2013.01.010
14. Abramoff B, Caldera FE. Osteoarthritis: pathology, diagnosis, and treatment options. Med Clin North Am. (2020) 104:293–311. doi: 10.1016/j.mcna.2019.10.007
15. Jiang Y. Osteoarthritis year in review 2021: biology. Osteoarthritis Cartilage. (2022) 30:207–15. doi: 10.1016/j.joca.2021.11.009
16. Ruiz-Romero C, Blanco FJ. Proteomics role in the search for improved diagnosis, prognosis and treatment of osteoarthritis. Osteoarthritis Cartilage. (2010) 18:500–9. doi: 10.1016/j.joca.2009.11.012
17. Boris Chan PM, Zhu L, Wen CY, Chiu KY. Subchondral bone proteomics in osteoarthritis: Current status and perspectives. J Orthop Translat. (2015) 3:71–7. doi: 10.1016/j.jot.2015.02.002
18. Xu Y, Sun Z. Molecular basis of Klotho: from gene to function in aging. Endocr Rev. (2015) 36:174–93. doi: 10.1210/er.2013-1079
19. Iijima H, Gilmer G, Wang K, Bean AC, He Y, Lin H, et al. Age-related matrix stiffening epigenetically regulates alpha-Klotho expression and compromises chondrocyte integrity. Nat Commun. (2023) 14:18. doi: 10.1038/s41467-022-35359-2
20. Gu Y, Ren K, Wang L, Yao Q. Loss of Klotho contributes to cartilage damage by derepression of canonical Wnt/beta-catenin signaling in osteoarthritis mice. Aging (Albany NY). (2019) 11:12793–809. doi: 10.18632/aging.102603
21. Chuchana P, Mausset-Bonnefont AL, Mathieu M, Espinoza F, Teigell M, Toupet K, et al. Secreted alpha-Klotho maintains cartilage tissue homeostasis by repressing NOS2 and ZIP8-MMP13 catabolic axis. Aging (Albany NY). (2018) 10:1442–53. doi: 10.18632/aging.101481
22. Prieto-Alhambra D, Judge A, Javaid MK, Cooper C, Diez-Perez A, Arden NK. Incidence and risk factors for clinically diagnosed knee, hip and hand osteoarthritis: influences of age, gender and osteoarthritis affecting other joints. Ann Rheum Dis. (2014) 73:1659–64. doi: 10.1136/annrheumdis-2013-203355
23. Skrivankova VW, Richmond RC, Woolf BAR, Yarmolinsky J, Davies NM, Swanson SA, et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomization: the STROBE-MR statement. JAMA. (2021) 326:1614–21. doi: 10.1001/jama.2021.18236
24. Burgess S, Foley CN, Allara E, Staley JR, Howson JMM. A robust and efficient method for Mendelian randomization with hundreds of genetic variants. Nat Commun. (2020) 11:376. doi: 10.1038/s41467-019-14156-4
25. Yamazaki Y, Imura A, Urakawa I, Shimada T, Murakami J, Aono Y, et al. Establishment of sandwich ELISA for soluble alpha-Klotho measurement: Age-dependent change of soluble alpha-Klotho levels in healthy subjects. Biochem Biophys Res Commun. (2010) 398:513–8. doi: 10.1016/j.bbrc.2010.06.110
26. Gergei I, Zheng J, Andlauer TFM, Brandenburg V, Mirza-Schreiber N, Muller-Myhsok B, et al. GWAS meta-analysis followed by Mendelian randomization revealed potential control mechanisms for circulating alpha-Klotho levels. Hum Mol Genet. (2022) 31:792–802. doi: 10.1093/hmg/ddab263
27. Donertas HM, Fabian DK, Valenzuela MF, Partridge L, Thornton JM. Common genetic associations between age-related diseases. Nat Aging. (2021) 1:400–12. doi: 10.1038/s43587-021-00051-5
28. Tachmazidou I, Hatzikotoulas K, Southam L, Esparza-Gordillo J, Haberland V, Zheng J, et al. Identification of new therapeutic targets for osteoarthritis through genome-wide analyses of UK Biobank data. Nat Genet. (2019) 51:230–6. doi: 10.1038/s41588-018-0327-1
29. Burgess S, Small DS, Thompson SG. A review of instrumental variable estimators for Mendelian randomization. Stat Methods Med Res. (2017) 26:2333–55. doi: 10.1177/0962280215597579
30. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. (2013) 37:658–65. doi: 10.1002/gepi.2013.37.issue-7
31. Xu Q, Wang J, Li H, Gao Y. Association between serum alpha-Klotho levels and osteoarthritis prevalence among middle-aged and older adults: an analysis of the NHANES 2007-2016. Rev Clin Esp (Barc). (2024) 224:366–78. doi: 10.1016/j.rceng.2024.04.012
32. Buchanan S, Combet E, Stenvinkel P, Shiels PG. Klotho, aging, and the failing kidney. Front Endocrinol (Lausanne). (2020) 11:560. doi: 10.3389/fendo.2020.00560
34. Zhu Z, Xia W, Cui Y, Zeng F, Li Y, Yang Z, et al. Klotho gene polymorphisms are associated with healthy aging and longevity: Evidence from a meta-analysis. Mech Ageing Dev. (2019) 178:33–40. doi: 10.1016/j.mad.2018.12.003
35. Kuro-o M, Matsumura Y, Aizawa H, Kawaguchi H, Suga T, Utsugi T, et al. Mutation of the mouse klotho gene leads to a syndrome resembling ageing. Nature. (1997) 390:45–51. doi: 10.1038/36285
36. Kurosu H, Yamamoto M, Clark JD, Pastor JV, Nandi A, Gurnani P, et al. Suppression of aging in mice by the hormone Klotho. Science. (2005) 309:1829–33. doi: 10.1126/science.1112766
37. Zhang F, Zhai G, Kato BS, Hart DJ, Hunter D, Spector TD, et al. Association between KLOTHO gene and hand osteoarthritis in a female Caucasian population. Osteoarthritis Cartilage. (2007) 15:624–9. doi: 10.1016/j.joca.2006.12.002
38. Martinez-Redondo P, Guillen-Guillen I, Davidsohn N, Wang C, Prieto J, Kurita M, et al. alphaKLOTHO and sTGFbetaR2 treatment counteract the osteoarthritic phenotype developed in a rat model. Protein Cell. (2020) 11:219–26. doi: 10.1007/s13238-019-00685-7
39. Nie X, Liu Z, Xie D, Sun Y. Inflammatory arthritis and eye diseases: a Mendelian randomization study. Front Endocrinol (Lausanne). (2023) 14:1251167. doi: 10.3389/fendo.2023.1251167
Keywords: α-klotho, osteoarthritis, NHANES, Mendelian randomization, causal effect
Citation: Li Z, Li Z, Cheng Q, Nie X, Cui Y, Du B, Ren T, Xu Y and Ma T (2024) α-klotho reduces susceptibility to osteoarthritis: evidence from cross-sectional studies and Mendelian randomization. Front. Endocrinol. 15:1450472. doi: 10.3389/fendo.2024.1450472
Received: 17 June 2024; Accepted: 31 October 2024;
Published: 19 November 2024.
Edited by:
Rita Cassia Menegati Dornelles, São Paulo State University, BrazilReviewed by:
Assunta Pozzuoli, University of Padua, ItalyYang Sun, The University of Auckland, New Zealand
Copyright © 2024 Li, Li, Cheng, Nie, Cui, Du, Ren, Xu and Ma. 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: Yibo Xu, bW9yaWFfMzNAMTYzLmNvbQ==; Teng Ma, Z3VrZW10QDE2My5jb20=
†These authors have contributed equally to this work