Genetic insights into serum cathepsins as diagnostic and therapeutic targets in knee and hip osteoarthritis.

Osteoarthritis (OA) is a chronic disease due to the deterioration of cartilage structure and function, involving the progressive degradation of the cartilage extracellular matrix. Cathepsins, lysosomal cysteine proteases, play pivotal roles in various biological and pathological processes, particularly in protein degradation. Excess cathepsins levels are reported to contribute to the development of OA. However, the causal relationship between the cathepsin family and knee and hip OA remains uncertain. Therefore, this study utilized bidirectional Mendelian Randomization (MR) analyses to explore this causal association. Our results indicated that elevated serum levels of cathepsin O increase the overall risk of knee OA, while increased serum levels of cathepsin H enhance the risk of hip OA. Conversely, the reverse MR analyses did not reveal a reverse causal relationship between them. In summary, OA in different anatomical locations may genetically result from pathological elevations in different serum cathepsin isoforms, which could be utilized as diagnostic and therapeutic targets in clinical practice.

Alterations in DNA methylation machinery in a rat model of osteoarthritis of the hip.

BACKGROUND: This study aimed to validate alterations in the gene expression of DNA methylation-related enzymes and global methylation in the peripheral blood mononuclear cell (PBMC) and synovial tissues of animal hip osteoarthritis (OA) models. METHODS: Animals were assigned to the control (no treatment), sham (25 µL of sterile saline), and OA (25 µL of sterile saline and 2 mg of monoiodoacetate) groups. Microcomputed tomography scan, histopathological assessment and pain threshold measurement were performed after induction. The mRNA expression of the DNA methylation machinery genes and global DNA methylation in the PBMC and hip synovial tissue were evaluated. RESULTS: The OA group presented with hip joint OA histopathologically and radiologically and decreased pain threshold. The mRNA expression of DNA methyltransferase (Dnmt 3a), ten-eleven translocation (Tet) 1 and Tet 3 in the synovial tissue of the OA group was significantly upregulated. Global DNA methylation in the synovial tissue of the OA group was significantly higher than that of the control and sham groups. CONCLUSIONS: The intra-articular administration of monoiodoacetate induced hip joint OA and decreased pain threshold. The DNA methylation machinery in the synovial tissues of hip OA was altered.

A systematic review and cross-database analysis of single nucleotide polymorphisms underlying hip morphology and osteoarthritis reveals shared mechanisms.

OBJECTIVE: Understanding the mechanisms of hip disease, such as osteoarthritis (OA), is crucial to advance their treatment. Such hip diseases often involve specific morphological changes. Genetic variations, called single nucleotide polymorphisms (SNPs), influence various hip morphological parameters. This study investigated the biological relevance of SNPs correlated to hip morphology in genome-wide association studies (GWAS). The SNP-associated genes were compared to genes associated with OA in other joints, aiming to see if the same genes play a role in both hip development and the risk of OA in other lower limb joints. METHODOLOGY: A systematic literature review was conducted to identify SNPs correlated with hip morphology, based on the Population, Intervention, Comparison, Outcome, and Study (PICOS) framework. Afterwards, Gene Ontology (GO) analysis was performed, using EnrichR, on the SNP-associated genes and compared with non-hip OA-associated genes, across different databases. RESULTS: Reviewing 49 GWAS identified 436 SNPs associated with hip joint morphology, encompassing variance in bone size, structure and shape. Among the SNP-associated genes, SOX9 plays a pivotal role in size, GDF5 impacts bone structure, and BMP7 affects shape. Overall, skeletal system development, regulation of cell differentiation, and chondrocyte differentiation emerged as crucial processes influencing hip morphology. Eighteen percent of GWAS-identified genes related to hip morphology were also associated with non-hip OA. CONCLUSION: Our findings indicate the existence of multiple shared genetic mechanisms across hip morphology and OA, highlighting the necessity for more extensive research in this area, as in contrast to the hip, the genetic background on knee or foot morphology remains largely understudied.

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

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.

Ultra-processed food consumption, genetic susceptibility, and the risk of hip/knee osteoarthritis.

BACKGROUND: The associations between ultra-processed food (UPF) consumption, genetic susceptibility, and the risk of osteoarthritis (OA) remain unknown. This study was to examine the effect of UPF consumption, genetic susceptibility, and their interactions on hip/knee OA. METHODS: Cohort analyses included 163,987 participants from the UK Biobank. Participants' UPF consumption was derived from their 24-h dietary recall using a questionnaire. Genetic risk scores (GRSs) of 70 and 83 single nucleotide polymorphisms (SNPs) for hip and knee OA were constructed. FINDINGS: After 1,461,447 person-years of follow-up, 11,540 patients developed OA. After adjustments, compared to participants in the low quartile of UPF consumption, those in the high quartile had a 10 % (hazard ratio [HR], 1.10; 95% confidence interval [CI], 1.03-1.18) increased risk of knee OA. No significant association was found between UPF consumption and hip OA. Replacing 20% of UPF diet weight with an equivalent proportion of unprocessed or minimally processed food caused a 6% (HR, 0.94; 95% CI, 0.89-0.98) decreased risk of knee OA, respectively. A significant interaction was found between UPF consumption, genetic predisposition, and the risk of knee OA (P = 0.01). Participants with lower OA-GRS scores experienced higher knee OA risks due to UPF consumption. INTERPRETATION: UPF consumption was associated with a higher risk of knee OA but not hip OA, particularly in those with lower genetic susceptibility. These results highlight the importance of reducing UPF consumption to prevent knee OA.

Circulating cytokines levels and osteoarthritis: A Mendelian randomization study.

BACKGROUND: Previous traditional observational studies have suggested the contribution of several cytokines and growth factors to the development of osteoarthritis (OA). This study aimed to determine the association of circulating cytokine and growth factor levels with OA. METHODS: We used two-sample Mendelian randomization (MR) to explore the causality between circulating cytokine and growth factor levels and OA [including knee or hip OA (K/HOA), knee OA (KOA), and hip OA (HOA)]. Summary level data for circulating cytokine and growth factor levels were sourced from a genome-wide association study (GWAS) involving 8,293 participants of Finnish ancestry. Single-nucleotide polymorphisms related to K/HOA (39,427 cases and 378,169 controls), KOA (24,955 cases and 378,169 controls), and HOA (15,704 cases and 378,169 controls) were obtained from a previous GWAS. The inverse variance weighted (IVW) method was primarily used for our MR analysis. For exposures to only one relevant SNP as IV, we used the Wald ratio as the major method to assess causal effects. We also conducted a series of sensitivity analyses to improve the robustness of the results. RESULTS: Circulating vascular endothelial growth factor levels were suggestively associated with an increased risk of K/HOA (odds ratio (OR) = 1.034; 95 % confidence interval (CI) = 1.013-1.055; P = 0.001), KOA (OR = 1.034; 95 % CI = 1.014-1.065; P = 0.002), and HOA (OR = 1.039; 95 % CI = 1.003-1.067; P = 0.034). Circulating interleukin (IL)-12p70 levels was suggestively associated with K/HOA (OR = 1.047; 95 % CI = 1.018-1.077; P = 0.001), KOA (OR = 1.058; 95 % CI = 1.022-1.095; P = 0.001), and HOA (OR = 1.044; 95 % CI = 1.000-1.091; P = 0.048). Circulating IL-18 levels were suggestively associated with HOA (OR = 1.068; 95 % CI = 1.014-1.125; P = 0.012). However, limited evidence exists to support causal genetic relationships between other circulating cytokines, growth factor levels and K/HOA, KOA, and HOA. CONCLUSIONS: Our MR analysis provides suggestive evidence of causal relationships between circulating cytokines and growth factors levels and OA, providing new insights into the etiology of OA.

Eligibility for knee arthroplasty is associated with increased risk of acquired hallux valgus - a Mendelian randomized study.

OBJECTIVE: Clinically, it has been found that patients undergoing knee replacement have a high incidence of concomitant hallux valgus. In this study, we analyzed whether patients with osteoarthritis who underwent surgery and those patient who did not have surgery had an increased risk of hallux valgus by Mendelian randomization and performed reverse causal analysis. DESIGN: Genomewide association study (GWAS) data for osteoarthritis, categorized by knee arthritis with joint replacement, knee arthritis without joint replacement, hip arthritis with joint replacement, and hip arthritis without joint replacement.And acquired hallux valgus were downloaded for Mendelian randomized studies. MR analysis was performed using inverse variance-weighted (IVW), weighted median, and MR-Egger methods. MR-egger regression, MR pleiotropic residuals and outliers (MR-presso), and Cochran's Q statistical methods were used to evaluate heterogeneity and pleiotropy. RESULTS: The IVW results indicate that, compared to healthy individuals, patients who meet the criteria for knee osteoarthritis joint replacement surgery have a significantly higher risk of acquired hallux valgus. There were no significant causal relationships found for the remaining results. No significant heterogeneity or multiplicity was observed in all the Mr analyses. CONCLUSION: Our study supports the increased risk of acquired hallux valgus in patients eligible for knee replacement. There is necessary for clinicians to be concerned about the hallux valgus status of patients undergoing knee arthroplasty.

Assessing the causal associations of different types of statins use and knee/hip osteoarthritis: A Mendelian randomization study.

OBJECTIVE: This study comprehensively evaluated the causal relationship between different types of statins use and knee/hip osteoarthritis (OA) using a two-sample and multivariate Mendelian randomization (MR) method. METHODS: MR analysis was conducted using publicly available summary statistics data from genome-wide association studies (GWAS) to assess the causal associations between total statins use (including specific types) and knee/hip OA. The primary analysis utilized the inverse variance-weighted (IVW) method, with sensitivity analysis conducted to assess robustness. Multivariable MR (MVMR) analysis adjusted for low-density lipoprotein cholesterol (LDL-C), intermediate-density lipoprotein cholesterol (IDL-C), high-density lipoprotein cholesterol (HDL-C), and body mass index (BMI). RESULTS: The MR analysis revealed a significant inverse association between genetically predicted total statins use and the risk of knee OA (OR = 0.950, 95%CI: 0.920-0.982, p = 0.002) as well as hip OA (OR = 0.932, 95%CI: 0.899-0.966, p <0.001). Furthermore, this study highlighted a reduced risk of knee/hip OA with the use of atorvastatin and simvastatin. Rosuvastatin use was associated with a decreased risk of hip OA but showed no association with knee OA. MVMR results indicated no correlation between exposure factors and outcomes after adjusting for LDL-C or IDL-C. HDL-C may not significantly contribute to statin-induced osteoarthritis, while BMI may play an important role. CONCLUSION: This study provides compelling evidence of the close relationship between statin use and a reduced risk of knee/hip OA, particularly with atorvastatin and simvastatin. LDL-C and IDL-C may mediate these effects. These findings have important implications for the clinical prevention and treatment of knee/hip OA.

Proteome-Wide Mendelian Randomization and Colocalization Analysis Identify Therapeutic Targets for Knee and Hip Osteoarthritis.

Osteoarthritis (OA) is a common degenerative disease. Although some biomarkers and drug targets of OA have been discovered and employed, limitations and challenges still exist in the targeted therapy of OA. Mendelian randomization (MR) analysis has been regarded as a reliable analytic method to identify effective therapeutic targets. Thus, we aimed to identify novel therapeutic targets for OA and investigate their potential side effects based on MR analysis. In this study, two-sample MR, colocalization analysis, summary-data-based Mendelian randomization (SMR) and Mendelian randomization phenome-wide association study (MR-PheWAS) were conducted. We firstly analyzed data from 4907 plasma proteins to identify potential therapeutic targets associated with OA. In addition, blood expression quantitative trait loci (eQTLs) data sources were used to perform additional validation. A protein-protein interaction (PPI) network was also constructed to delve into the interactions among identified proteins. Then, MR-PheWASs were utilized to assess the potential side effects of core therapeutic targets. After MR analysis and FDR correction, we identified twelve proteins as potential therapeutic targets for knee OA or hip OA. Colocalization analysis and additional validation supported our findings, and PPI networks revealed the interactions among identified proteins. Finally, we identified MAPK3 (OR = 0.855, 95% CI: 0.791-0.923, p = 6.88 × 10-5) and GZMK (OR = 1.278, 95% CI: 1.131-1.444, p = 8.58 × 10-5) as the core therapeutic targets for knee OA, and ITIH1 (OR = 0.847, 95% CI: 0.784-0.915, p = 2.44 × 10-5) for hip OA. A further MR phenome-wide association study revealed the potential side effects of treatments targeting MAPK3, GZMK, and ITIH1. This comprehensive study indicates twelve plasma proteins with potential roles in knee and hip OA as therapeutic targets. This advancement holds promise for the progression of OA drug development, and paves the way for more efficacious treatments of OA.

Indirect influence of microRNA-146a on the association of IL-6 and TNF-α genetic polymorphisms with the increased risk of hip osteoarthritis.

Primary osteoarthritis (POA) is a complex hereditary disease that involves the interplay between genetics and epigenetics. MicroRNA molecules play important roles in epigenetic mechanisms. MicroRNA-146a (miR-146a) is a negative regulator of the immune response in osteoarthritis (OA). So, variations in the miR-146a gene could affect OA risk. The aim of this study was to investigate the relationships between single nucleotide polymorphisms (SNPs) in the miR-146a, interleukin-6 (IL-6), Toll-like receptor 10 (TLR10), and tumor necrosis factor-alpha (TNFA) genes and the risk for development of advanced-stage primary hip osteoarthritis (PHOA) and primary knee osteoarthritis (PKOA) in the Croatian population. A total of 609 POA patients and 656 healthy donors were genotyped for SNPs in the miR-146a (rs2910164, G>C). Since we used same patients and controls as two studies before us, we already had information about IL-6 (rs1800795, C>G), TLR10 (rs11096957, C>T), and TNFA (rs1800629, C>T) genotypes of our subjects. None of the differences were statistically significant comparing either allelic or genotypic frequencies of miR-146a SNP rs2910164 (G>C) between the PHOA and PKOA patients and controls. However, we found a significant association with risk to PHOA for the combination of genotypes (stratified miR-146a genotype with the IL-6, and stratified miR-146a genotype with the TNFA). In a multifactorial disease such as POA, we have shown the indirect relevance of a second modifying factor (miR-146a), which apparently contributes to the overall risk of PHOA. There was no risk association with the PKOA, indicating that these two localities (hip and knee) might have different risk-modifying factors.

Causality between Ankylosing Spondylitis and osteoarthritis in European ancestry: a bidirectional Mendelian randomization study.

Objective: To explore the bidirectional causal relationship between Ankylosing Spondylitis (AS) and Osteoarthritis (OA) at the genetic level within the European ancestry. Methods: We implemented a series of quality control steps to select instrumental variables (IVs) related to the exposure. We conducted two-sample Mendelian randomization (MR) using the inverse-variance weighted method as the primary approach. We adjusted significance levels using Bonferroni correction, assessed heterogeneity using Cochrane's Q test. Sensitivity analysis was conducted through leave-one-out method. Additionally, external datasets and relaxed IV selection criteria were employed, and multivariate MR analyses were performed for validation purposes. Finally, Bayesian colocalization (COLOC) analysis identified common genes, validating the MR results. Results: The investigation focused on the correlation between OA and AS in knee, hip, and hand joints. MR results revealed that individuals with AS exhibit a decreased risk of knee OA (OR = 0.9882, 95% CI: 0.9804-0.9962) but no significant increase in the risk of hip OA (OR = 0.9901, 95% CI: 0.9786-1.0018). Conversely, AS emerged as a risk factor for hand OA (OR = 1.0026, 95% CI: 1.0015-1.0036). In reverse-direction MR analysis, OA did not significantly influence the occurrence of AS. Importantly, minimal heterogeneity was observed in our MR analysis results (p > 0.05), and the robustness of these findings was confirmed through sensitivity analysis and multivariate MR analysis. COLOC analysis identified four colocalized variants for AS and hand OA (rs74707996, rs75240935, rs181468789, and rs748670681). Conclusion: In European population, individuals with AS have a relatively lower risk of knee OA, whereas AS serves as a risk factor for hand OA. However, no significant causal relationship was found between AS and hip OA. Additionally, it offers novel insights into genetic research on AS and OA.

Genomic heterozygosity is associated with a lower risk of osteoarthritis.

BACKGROUND: Genomic heterozygosity has been shown to confer a health advantage in humans and play a protective role in complex diseases. Given osteoarthritis (OA) is a highly polygenic disease, we set out to determine if an association exists between OA and genomic heterozygosity. RESULTS: End-stage knee and hip OA patients and healthy controls were recruited from the Newfoundland and Labrador (NL) population. The Arthritis Research UK Osteoarthritis Genetics (arcOGEN) consortium database was utilized as a replication cohort. DNA was extracted from blood samples and genotyped. Individual rates of observed heterozygosity (HetRate) and heterozygosity excess (HetExcess) relative to the expected were mathematically derived, and standardized to a z-score. Logistic regression modeling was used to examine the association between OA and HetRate or HetExcess. A total of 559 knee and hip OA patients (mean age 66.5 years, body mass index (BMI) 33.7 kg/m2, and 55% females) and 118 healthy controls (mean age 56.4 years, BMI 29.5 kg/m2, and 59% female) were included in the NL cohort analysis. We found that OA had an inverse relationship with HetRate and HetExcess with odds ratios of 0.64 (95% CI: 0.45-0.91) and 0.65 (95% CI: 0.45-0.93) per standard deviation (SD), respectively. The arcOGEN data included 2,019 end-stage knee and hip OA patients and 2,029 healthy controls, validating our findings with HetRate and HetExcess odds ratios of 0.60 (95% CI: 0.56-0.64) and 0.44 (95% CI: 0.40-0.47) per SD, respectively. CONCLUSIONS: Our results are the first to clearly show evidence, from two separate cohorts, that reduced genomic heterozygosity confers a risk for the future development of OA.

Characteristics of the synovial microenvironment and synovial mesenchymal stem cells with hip osteoarthritis of different bone morphologies.

BACKGROUND: Variations in bone morphology in patients with hip osteoarthritis (HOA) can be broadly categorized into three types: atrophic, normotrophic, and hypertrophic. Despite the investigations examining clinical elements, such as bone morphology, pain, and range of motion, our understanding of the pathogenesis of HOA remains limited. Previous studies have suggested that osteophytes typically originate at the interface of the joint cartilage, periosteum, and synovium, potentially implicating synovial mesenchymal stem cells (SMSCs) in the process. This study aimed to investigate the potential factors that drive the development of bone morphological features in HOA by investigating the characteristics of the synovium, differentiation potential of SMSCs, and composition of synovial fluid in different types of HOA. METHODS: Synovial tissue and fluid were collected from 30 patients who underwent total hip arthroplasty (THA) with the variable bone morphology of HOA patients. RNA sequencing analysis and quantitative reverse transcription-polymerase chain reaction (RT-qPCR) were performed to analyse the genes in the normotrophic and hypertrophic synovial tissue. SMSCs were isolated and cultured from the normotrophic and hypertrophic synovial tissues of each hip joint in accordance with the variable bone morphology of HOA patients. Cell differentiation potential was compared using differentiation and colony-forming unit assays. Cytokine array was performed to analyse the protein expression in the synovial fluid. RESULTS: In the RNA sequencing analysis, 103 differentially expressed genes (DEGs) were identified, predominantly related to the interleukin 17 (IL-17) signalling pathway. Using a protein-protein interaction (PPI) network, 20 hub genes were identified, including MYC, CXCL8, ATF3, NR4A1, ZC3H12A, NR4A2, FOSB, and FOSL1. Among these hub genes, four belonged to the AP-1 family. There were no significant differences in the tri-lineage differentiation potential and colony-forming capacity of SMSCs. However, RT-qPCR revealed elevated SOX9 expression levels in synovial tissues from the hypertrophic group. The cytokine array demonstrated significantly higher levels of CXCL8, MMP9, and VEGF in the synovial fluid of the hypertrophic group than in the normotrophic group, with CXCL8 and MMP9 being significantly expressed in the hypertrophic synovium. CONCLUSION: Upregulation of AP-1 family genes in the synovium and increased concentrations of CXCL8, MMP9, and VEGF were detected in the synovial fluid of the hypertrophic group of HOA patients, potentially stimulating the differentiation of SMSCs towards the cartilage and thereby contributing to severe osteophyte formation.

Family History of Developmental Dysplasia of the Hip is a Risk Factor for the Progression of Hip Osteoarthritis.

BACKGROUND: Developmental dysplasia of the hip (DDH) is considered to have genetic predisposition and presents many intrafamilial occurrences. However, there is no report that evaluates the effect of DDH family history on the progression after the onset of hip osteoarthritis (OA). METHODS: Medical interviews about detailed clinical information including family history were conducted on 298 consecutive patients who had undergone surgery for OA due to DDH. Clinical or radiographic items that are associated with the severity of DDH (total hip arthroplasty [THA], involvement of bilateral DDH, onset age of hip pain, and three radiological indices of DDH: center-edge angle, sharp angle, and acetabular roof obliquity) were collected and evaluated in multivariate analyses for their associations with DDH family history in a qualitative or quantitative manner. Survival time analyses for THA as the endpoint was also performed to evaluate the effects of DDH family history on the progression of OA. RESULTS: The DDH family history showed significant associations with bilateral involvement of DDH (odds ratio = 2.09 [95% confidence interval {CI} 1.05 to 4.16]; P = .037), early onset of hip pain (P = .0065), and radiological severity of DDH (P = .016). The DDH family history showed a significant association with undergoing THA (odds ratio = 2.25 [95% CI 1.09 to 4.66]; P = .029), further supported by the Cox regression analyses (hazards ratio = 1.56 [95% CI 1.15 to 2.11]; P = .0044). CONCLUSION: A DDH family history is a risk factor for the progression of hip OA. Stronger genetic predisposition to DDH leads to faster onset and progression of hip OA.

Genetic Influence on Osteoarthritis Versus Other Rheumatic Diseases.

OBJECTIVE: We aimed to compare the genetic contribution to osteoarthritis (OA) versus other rheumatic/musculoskeletal diseases (RMDs) in the same population and to explore the role for any shared genetics between OA and other RMDs. METHODS: In 59,970 Swedish twins aged 35 years or older, we estimated the heritability (the variance explained by genetic factors) of OA in peripheral joints, back and neck pain, shoulder pain (adhesive capsulitis, impingement syndrome, etc), rheumatoid arthritis, spondyloarthritis (SpA) and psoriatic arthritis, myalgia, and osteoporosis diagnosed in specialist and inpatient care. We also studied how much covariance between OA and each of the RMDs could be explained by genetics by studying phenotypic correlations in bivariate classical twin models. RESULTS: Any-site OA and hip OA (50% and 64%) were among the most heritable RMDs (as compared with 23% for fibromyalgia [lowest] and 63% for SpA [highest]). The highest phenotypic correlations were between OA (any joint site) and shoulder pain in the same individual (r = 0.33, 95% confidence interval 0.31-0.35), of which 70% (95% confidence interval 52-88) could be explained by shared genetics. The phenotypic correlation between OA and back/neck pain was r = 0.25, with 25% to 75% explained by genetics. Phenotypic correlations between OA and each of the other RMDs were lower (r ~ 0.1 to r ~ 0.2), with inconclusive sources of variation. CONCLUSION: OA has relatively large heritability as compared with other RMDs. The coexistence of OA and shoulder pain, as well as back pain, was common and could often be explained by genetic factors. Findings imply similar etiologies of OA and several pain conditions.

Causal relationships between anthropometric traits, bone mineral density, osteoarthritis and spinal stenosis: A Mendelian randomisation investigation.

OBJECTIVE: Spinal stenosis is a common condition among older individuals, with significant morbidity attached. Little is known about its risk factors but degenerative conditions, such as osteoarthritis (OA) have been identified for their mechanistic role. This study aims to explore causal relationships between anthropometric risk factors, OA, and spinal stenosis using Mendelian randomisation (MR) techniques. DESIGN: We applied two-sample MR to investigate the causal relationships between genetic liability for select risk factors and spinal stenosis. Next, we examined the genetic relationship between OA and spinal stenosis with linkage disequilibrium score regression and Causal Analysis Using Summary Effect estimates MR method. Finally, we used multivariable MR (MVMR) to explore whether OA and body mass index (BMI) mediate the causal pathways identified. RESULTS: Our analysis revealed strong evidence for the effect of higher BMI (odds ratio [OR] = 1.54, 95%CI: 1.41-1.69, p-value = 2.7 × 10-21), waist (OR = 1.43, 95%CI: 1.15-1.79, p-value = 1.5 × 10-3) and hip (OR = 1.50, 95%CI: 1.27-1.78, p-value = 3.3 × 10-6) circumference on spinal stenosis. Strong evidence of causality was also observed for higher bone mineral density (BMD): total body (OR = 1.21, 95%CI: 1.12-1.29, p-value = 1.6 × 10-7), femoral neck (OR = 1.35, 95%CI: 1.09-1.37, p-value = 7.5×10-7), and lumbar spine (OR = 1.38, 95%CI: 1.25-1.52, p-value = 4.4 × 10-11). We detected high genetic correlations between spinal stenosis and OA (rg range: 0.47-0.66), with Causal Analysis Using Summary Effect estimates results supporting a causal effect of OA on spinal stenosis (ORallOA = 1.6, 95%CI: 1.41-1.79). Direct effects of BMI, BMD on spinal stenosis remained after adjusting for OA in the MVMR. CONCLUSIONS: Genetic susceptibility to anthropometric risk factors, particularly higher BMI and BMD can increase the risk of spinal stenosis, independent of OA status. These results may inform preventative strategies and treatments.

Osteoporosis and osteoarthritis: a bi-directional Mendelian randomization study.

OBJECTIVE: To investigate the causal relationship between low bone mineral density (BMD) and osteoarthritis (OA) using Mendelian randomization (MR) design. METHODS: Two-sample bi-directional MR analyses were performed using summary-level information on OA traits from UK Biobank and arcOGEN. Sensitivity analyses including MR-Egger, simple median, weighted median, MR pleiotropy residual sum, and outlier approaches were utilized in conjunction with inverse variance weighting (IVW). Gene ontology (GO) enrichment analyses and expression quantitative trait locus (eQTL) colocalization analyses were used to investigate the potential mechanism and shared genes between osteoporosis (OP) and OA. RESULTS: The IVW method revealed that genetically predicted low femoral neck BMD was significantly linked with hip (ß = 0.105, 95% CI: 0.023-0.188) and knee OA (ß = 0.117, 95% CI: 0.049-0.184), but not with other site-specific OA. Genetically predicted low lumber spine BMD was significantly associated with OA at any sites (ß = 0.048, 95% CI: 0.011-0.085), knee OA (ß = 0.101, 95% CI: 0.045-0.156), and hip OA (ß = 0.150, 95% CI: 0.077-0.224). Only hip OA was significantly linked with genetically predicted reduced total bone BMD (ß = 0.092, 95% CI: 0.010-0.174). In the reverse MR analyses, no evidence for a causal effect of OA on BMD was found. GO enrichment analysis and eQTL analysis illustrated that DDN and SMAD-3 were the most prominent co-located genes. CONCLUSIONS: These findings suggested that OP may be causally linked to an increased risk of OA, indicating that measures to raise BMD may be effective in preventing OA. More research is required to determine the underlying processes via which OP causes OA.

Causal effect of physical activity and sedentary behaviors on the risk of osteoarthritis: a univariate and multivariate Mendelian randomization study.

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

Assessment of causal effects of physical activity on the risk of osteoarthritis: a two-sample Mendelian randomization study.

BACKGROUND: Growing evidence supports an association between physical activity (PA) and the risk of osteoarthritis (OA), but this may be influenced by confounding and reverse causality. Therefore, we performed a two-sample Mendelian randomization (MR) analysis to reveal the causal relationship between PA and OA. METHODS: MR was performed to explore the causation of PA and OA with genetic variants as instrumental variables. The genetic variants were derived from the summary statistics of a large genome-wide association study meta-analysis based on the European population (n = 661,399), including self-reported leisure screen time (LST) and moderate-to-vigorous physical activity (MVPA), and Arthritis Research UK Osteoarthritis Genetics Consortium cohorts (417,596, 393,873 and 403,124 for overall, hip and knee OA, respectively). The major MR analysis used in this work was the inverse variance weighted (IVW) approach, and sensitivity, pleiotropy, and heterogeneity studies were performed to evaluate the validity of the findings. RESULTS: IVW estimates indicated that LST had a risk effect on overall OA (odds ratio (OR) = 1.309, 95% confidence interval (CI): 1.198-1.430, P = 2.330 × 10-9), hip OA (OR = 1.132, 95% CI: 1.009-1.269, P = 0.034) and knee OA (OR = 1.435. 95% CI: 1.286-1.602, P = 1.225 × 10-10). In contrast, no causal relationship was found between MVPA and OA (overall OA: OR = 0.895, 95% CI: 0.664-1.205, P = 0.465; hip OA: OR = 1.189, 95% CI: 0.792-1.786, P = 0.404; knee OA: OR = 0.707, 95% CI: 0.490 -1.021, P = 0.064). In addition, we observed significant heterogeneity in instrumental variables, but no horizontal pleiotropy was detected. CONCLUSIONS: Recent findings demonstrated a protective impact of reducing LST on OA, independent of MVPA. This provides valuable insights into the role of physical activity in OA and offers lifestyle recommendations, such as reducing recreational sedentary behaviors and promoting appropriate exercise, for individuals at risk of OA.

Causal relationships between CD25 on immune cells and hip osteoarthritis.

Objectives: Previous research has indicated a potential association between immune factors and osteoarthritis (OA), but the causal relationship between CD25 expression on immune cells and hip OA remains enigmatic. To shed light on this relationship, this study utilized the two-sample Mendelian Randomization (MR) method. Methods: Leveraging genome-wide association studies (GWAS) data from the UK Biobank and arcOGEN, the investigation encompasses a substantial European cohort comprising 15,704 hip OA cases and 378,169 controls. Genetic insights into CD25 stem from a subgroup of 3,757 individuals with European ancestry, encompassing 77 CD25-related traits. Several MR methods were applied, and robustness was assessed through heterogeneity and sensitivity analysis. Results: Among the 77 traits examined, 66 shared the same single nucleotide polymorphisms (SNPs) with hip OA. Of these, 7 CD25-related traits were found to be causally associated with hip OA (adjusted P><0.05), with F-statistics ranging from 33 to 122. These traits are specifically related to CD4+CD25+ T cells, exhibiting odds ratios (OR) and 95% confidence intervals (CI) less than 1. Notably, no causal link was discerned with the CD8+CD25+ T cell subset. Within absolute count (AC) and relative count (RC) trait types, a significant causal relationship was observed solely between CD4+CD25+ T cells and hip OA, without subtype localization. A more intricate examination of CD25 expression levels within the CD4+CD25+ T cell subset revealed a correlation with the CD39+ regulatory T (Treg) subset and hip OA, particularly within the CD39+ activated Treg subset. Furthermore, a notable causal relationship emerged between CD25 expression levels in the CD45RA- not Treg subset and hip OA. However, no significant causal link was established with any subsets of B cells. Conclusion: The genetic prediction suggests that CD25, particularly within the realm of CD4+CD25+ T cells, may exert a protective influence against the development of hip OA. These findings provide a novel therapeutic approach for the prevention and treatment of hip OA.