Genome-wide analyses using UK Biobank data provide insights into the genetic architecture of osteoarthritis

Integrin signalling in joint development, homeostasis and osteoarthritis

Integrins are key regulators of cell-matrix interactions during joint development and joint tissue homeostasis, as well as in the development of osteoarthritis (OA). The signalling cascades initiated by the interactions of integrins with a complex network of extracellular matrix (ECM) components and intracellular adaptor proteins orchestrate cellular responses necessary for maintaining joint tissue integrity. Dysregulated integrin signalling, triggered by matrix degradation products such as matrikines, disrupts this delicate balance, tipping the scales towards an environment conducive to OA pathogenesis. The interplay between integrin signalling and growth factor pathways further underscores the multifaceted nature of OA. Moreover, emerging insights into the role of endocytic trafficking in regulating integrin signalling add a new layer of complexity to the understanding of OA development. To harness the therapeutic potential of targeting integrins for mitigation of OA, comprehensive understanding of their molecular mechanisms across joint tissues is imperative. Ultimately, deciphering the complexities of integrin signalling will advance the ability to treat OA and alleviate its global burden.

The relationship between obstructive sleep apnea and osteoarthritis: evidence from an observational and Mendelian randomization study

Objectives

Obstructive sleep apnea (OSA) and osteoarthritis (OA) are common comorbidities that significantly impact individuals' quality of life. However, the relationship between OSA and OA remains unclear. This study aims to explore the connection between OSA and OA and evaluate causality using Mendelian randomization (MR).

Methods

A total of 12,454 participants from the National Health and Nutrition Examination Survey (2009-2012) were included. OSA participants were identified based on self-reported interviews. The association between OA and OSA was assessed through multivariable logistic regression analysis. A two-sample MR was employed to investigate the relationship between OSA and OA, specifically hip OA and knee OA, utilizing the inverse variance-weighted (IVW) approach.

Results

Based on the observational study, individuals with OSA exhibited a higher risk of OA (OR = 1.67, 95% CI = 1.40-1.98). IVW demonstrated that the risk of OA (OR = 1.13, 95% CI: 1.05-1.21, p = 0.001), hip OA (OR = 1.11, 95% CI: 1.04-1.18, p = 0.002), and knee OA (OR = 1.08, 95% CI: 1.02-1.14, p = 0.005) was significantly associated with OSA. Reverse MR analyses indicated no effect of OA on OSA. Additionally, body mass index (BMI) was found to mediate 36.9% (95% CI, 4.64-73.2%, p = 0.026) of the OSA effects on OA risk.

Conclusion

The cross-sectional observational analysis unveiled noteworthy associations between OSA and OA. Meanwhile, findings from the MR study provide support for a causal role.

Unveiling inflammatory and prehypertrophic cell populations as key contributors to knee cartilage degeneration in osteoarthritis using multi-omics data integration

OBJECTIVES

Single-cell and spatial transcriptomics analysis of human knee articular cartilage tissue to present a comprehensive transcriptome landscape and osteoarthritis (OA)-critical cell populations.

METHODS

Single-cell RNA sequencing and spatially resolved transcriptomic technology have been applied to characterise the cellular heterogeneity of human knee articular cartilage which were collected from 8 OA donors, and 3 non-OA control donors, and a total of 19 samples. The novel chondrocyte population and marker genes of interest were validated by immunohistochemistry staining, quantitative real-time PCR, etc. The OA-critical cell populations were validated through integrative analyses of publicly available bulk RNA sequencing data and large-scale genome-wide association studies.

RESULTS

We identified 33 cell population-specific marker genes that define 11 chondrocyte populations, including 9 known populations and 2 new populations, that is, pre-inflammatory chondrocyte population (preInfC) and inflammatory chondrocyte population (InfC). The novel findings that make this an important addition to the literature include: (1) the novel InfC activates the mediator MIF-CD74; (2) the prehypertrophic chondrocyte (preHTC) and hypertrophic chondrocyte (HTC) are potentially OA-critical cell populations; (3) most OA-associated differentially expressed genes reside in the articular surface and superficial zone; (4) the prefibrocartilage chondrocyte (preFC) population is a major contributor to the stratification of patients with OA, resulting in both an inflammatory-related subtype and a non-inflammatory-related subtype.

CONCLUSIONS

Our results highlight InfC, preHTC, preFC and HTC as potential cell populations to target for therapy. Also, we conclude that profiling of those cell populations in patients might be used to stratify patient populations for defining cohorts for clinical trials and precision medicine.

Causal relationship between depression and metabolic dysfunction-associated steatotic liver disease: a bidirectional Mendelian randomized study

Background

With the global rise in obesity, metabolic dysfunction-associated steatotic liver disease (MASLD) has emerged as the most common chronic liver disease. Concurrently, depression is a highly prevalent mental disorder. As the incidence of MASLD and depression continues to increase, a growing body of research indicates a potential association between the two conditions. However, the direction of causality between depression and MASLD remains uncertain. To address this gap, our study utilizes a two-sample Mendelian randomization (MR) approach to explore the bidirectional causal relationship between depression and MASLD.

Methods

We extracted single nucleotide polymorphisms (SNPs) associated with depression and MASLD from pooled data of genome-wide association studies (GWAS). A comprehensive assessment of possible causality was also performed. Possible mediating effects of liver enzymes on MASLD were also assessed.

Results

A total of three GWAS pooled data on depression as well as GWAS data related to MASLD and GWAS data on four liver enzymes were used in this study. Our findings indicated a strong causal relationship between depression and MASLD (OR, 1.557; 95% CI, 1.097-2.211; P = 0.016). And we found a mediating effect of gamma-glutamyl transferase (GGT), alanine aminotransferase (ALT) and aspartate aminotransferase (AST). ALT 10% (95% CI: 7% - 13%, P< 0.0002). AST, 4.14% (95% CI: 2.34% - 5.94%, P < 0.05). GGT 0.19% (95% CI: 0.15% - 0.22%, P< 0.000000002). However, we did not find a mediating effect of alkaline phosphatase (ALP). Our inverse MR analysis did not reveal any causal relationship between MASLD and depression.

Conclusions

The MR analysis revealed a positive causal relationship between depression and MASLD, while no reverse causal relationship was identified. Liver enzymes may mediate the role between depression and MASLD.

Assessing the causal relationship of birth weight and childhood obesity on osteoarthritis: a Mendelian randomization study

Obesity is associated with osteoarthritis (OA), but few studies have used fetal origin to explore the association. Our study aims to disentangle the causality between birth weight, childhood obesity, and adult OA using Mendelian randomization (MR). We identified single nucleotide polymorphisms (SNPs) related to birth weight (n = 298,142) and childhood obesity (n = 24,160) from two genome-wide association studies contributed by the Early Growth Genetics Consortium. Summary statistics of OA and its phenotypes (knee, hip, spine, hand, thumb, and finger OA) from the Genetics of Osteoarthritis Consortium (n = 826,690) were used to estimate the effects of SNPs on OA. Multivariable MR (MVMR) was conducted to investigate the independent effects of exposures. It turned out that genetically predicted standard deviation increase in birth weight was not associated with OA. In contrast, there was a marginally positive effect of childhood obesity on total [odds ratio (OR) = 1.07, 95% confidence interval (CI) = 1.00, 1.15 using IVW], knee (OR = 1.13, 95% CI = 1.05, 1.22 using weighted median), hip (OR = 1.13, 95% CI = 1.04, 1.24 using IVW), and spine OA (OR = 1.12, 95% CI = 1.03, 1.22 using IVW), but not hand, thumb, or finger OA. MVMR indicated a potential adulthood body mass index-dependent causal pathway between childhood obesity and OA. In conclusion, no association of birth weight with OA was suggested. Childhood obesity, however, showed a causality with OA in weight-bearing joints, which seems to be a general association of obesity with OA.

Dual-energy X-ray absorptiometry derived knee shape may provide a useful imaging biomarker for predicting total knee replacement: Findings from a study of 37,843 people in UK Biobank

Objective

We aimed to create an imaging biomarker for knee shape using knee dual-energy x-ray absorptiometry (DXA) scans and investigate its potential association with subsequent total knee replacement (TKR), independently of radiographic features of knee osteoarthritis and established risk factors.

Methods

Using a 129-point statistical shape model, knee shape (expressed as a B-score) and minimum joint space width (mJSW) of the medial joint compartment (binarized as above or below the first quartile) were derived. Osteophytes were manually graded in a subset of images and an overall score was assigned. Cox proportional hazards models were used to examine the associations of B-score, mJSW and osteophyte score with TKR risk, adjusting for age, sex, height and weight.

Results

The analysis included 37,843 individuals (mean age 63.7 years). In adjusted models, B-score was associated with TKR: each unit increase in B-score, reflecting one standard deviation from the mean healthy shape, corresponded to a hazard ratio (HR) of 2.25 (2.08, 2.43), while a lower mJSW had a HR of 2.28 (1.88, 2.77). Among the 6719 images scored for osteophytes, mJSW was replaced by osteophyte score in the most strongly predictive model for TKR. In ROC analyses, a model combining B-score, osteophyte score, and demographics outperformed a model including demographics alone (AUC ​= ​0.87 vs 0.73).

Conclusions

Using statistical shape modelling, we derived a DXA-based imaging biomarker for knee shape that was associated with kOA progression. When combined with osteophytes and demographic data, this biomarker may help identify individuals at high risk of TKR, facilitating targeted interventions.

Clinical Assessment of Osteoarthritis Pain: Contemporary Scenario, Challenges, and Future Perspectives

The multifaceted nature of osteoarthritis (OA) pain presents a challenge in understanding and managing the condition. The diverse pain experiences, progression rates, individual responses to treatments, and complex disease mechanisms contribute to heterogeneity in the clinical studies outcomes. The lack of a standardized methodology for assessing and classifying OA pain challenges healthcare practitioners. This complicates the establishment of universally applicable protocols or standardized guidelines for treatment. This article explores the heterogeneity observed in clinical studies evaluating OA pain treatments, highlighting the necessity for refined methodologies, personalized patient categorization, and consistent outcome measures. It discusses the role of the multidimensional nature of OA pain, underlying pain mechanisms, and other contributing factors to the heterogeneity in outcome measures. Addressing these variations is crucial to establishing a more consistent framework for evidence-based treatments and advancing care of the patient with OA pain.

Association between omega-3 polyunsaturated fatty acids and osteoarthritis: results from the NHANES 2003-2016 and Mendelian randomization study

BACKGROUND

Omega-3 polyunsaturated fatty acids (omega-3 PUFAs) exhibit potential as therapeutics for a variety of diseases. This observational and Mendelian randomization (MR) study aims to explore the relationship between omega-3 PUFAs and osteoarthritis (OA).

METHODS

Excluding individuals under 20 years old and those with missing data on relevant variables in the National Health and Nutrition Examination Survey (NHANES) spanning from 2003 to 2016, a total of 22 834 participants were included in this cross-sectional study. Weighted multivariable-adjusted logistic regression was used to estimate the association between omega-3 PUFAs and OA in adults. Moreover, restricted cubic splines were utilized to examine the dose-response relationship between omega-3 PUFAs and OA. To further investigate the potential causal relationship between omega-3 PUFAs and OA risk, a two-sample MR study was conducted. Furthermore, the robustness of the findings was assessed using various methods.

RESULTS

Omega-3 PUFAs intake were inversely associated with OA in adults aged 40 ∼ 59 after multivariable adjustment [Formula: see text], with a nonlinear relationship observed between omega-3 PUFAs intake and OA [Formula: see text]. The IVW results showed there was no evidence to suggest a causal relationship between omega-3 PUFAs and OA risk [Formula: see text].

CONCLUSIONS

Omega-3 PUFAs were inversely associated with OA in adults aged 40 ∼ 59. However, MR studies did not confirm a causal relationship between the two.

Global burden and socioeconomic impact of knee osteoarthritis: a comprehensive analysis

Objective

To report the trend changes of the prevalence, incidence, and disability-adjusted life years (DALYs) of knee osteoarthritis (KOA) according to age, sex, sociodemographic index (SDI), and income.

Methods

This analysis used estimates from the broader Global Burden of Disease (GBD) study 2019, which covered 201 countries from 1990 to 2019. National prevalence, incidence, and DALYs were shown by using ggplot2 and maps packages. Five-year intervals were used for age groupings. The Commonwealth and the World Bank income classifications were used for income grouping.

Results

Globally, there were ~364.58 million prevalent cases (females: 225.16 million), 29.51 million incident cases (females: 17.53 million), and 11.53 million DALYs (females: 7.09 million) due to KOA in 2019. The Western Pacific Region had a high endemicity of ~7,319.87 cases per 100,000 people (7.64%). Japan's prevalence rate (12,610.12 cases per 100,000 population) was 10 times that of Somalia (1,178.23) in 2019. In 200 countries (except the United Arab Emirates), the prevalence, incidence, and DALY rates of KOA in 2019 were higher among females than among males and increased with age up to the oldest age group. The prevalence was highest in the high-middle SDI countries, and the incidence and DALYs were highest in the middle SDI countries.

Conclusion

There was a large burden of KOA worldwide, with some notable intercountry variation. Some countries had 10 times the burden of other countries. Increasing population awareness regarding the prevalence, incidence, and DALYs of KOA with a focus on the population in the Western Pacific Region is needed, particularly for older females. informs health policy development, and contributes to improving the efficiency, equity, and effectiveness of healthcare systems.

Identification and Verification of Endoplasmic Reticulum Stress-Related Genes as Novel Signatures for Osteoarthritis Diagnosis and Therapy: A Bioinformatics Analysis-Oriented Pilot Study

BACKGROUND AND PURPOSE

Endoplasmic reticulum stress (ERS) has been reported to be closely associated with the development of osteoarthritis (OA), but the underlying mechanisms are not fully delineated. The present study was designed to investigate the involvement of ERS-related genes in regulating OA progression.

METHODS

The expression profiles of OA patients and normal people were downloaded from the gene expression omnibus (GEO) database. The differentially expressed genes (DEGs) in datasets GSE55457 and GSE55235 were screened and identified by R software with the construction of the protein-protein interaction (PPI) networks. Through the STRING and Venn diagram analysis, hub ERS-related genes were obtained. Gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses were performed. Biomarkers with high diagnostic values of osteoarthritis (OA) were studied. The hematoxylin and eosin (H&E) staining and micro-CT were applied to evaluate the establishment of the OA model. The expression levels of biomarkers were validated with the use of reverse transcription‑quantitative polymerase chain reaction (RT-qPCR) and western blot. Finally, we evaluated the correlations of hub ERS-related genes with the immune infiltration cells via the CIBERSORT algorithm.

RESULTS

A total of 60 downregulated and 52 upregulated DEGs were identified, and the following GO and KEGG pathway analyses verified that those DEGs were mainly enriched in biological process (BP), cellular component (CC), molecular function (MF), and inflammation-associated signal pathways. Interestingly, among all the DEGs, six ER stress-associated genes, including activating transcription factor 3 (ATF3), DEAD-Box Helicase 3 X-Linked (DDX3X), AP-1 transcription factor subunit (JUN), eukaryotic initiation factor 4 (EIF4A1), KDEL endoplasmic reticulum protein retention receptor 3 (KDELR3), and vascular endothelial growth factor A (VEGFA), were found to be closely associated with OA progression, and the following RT-qPCR and Western Blot analysis confirmed that DDX3X, JUN, and VEGFA were upregulated, whereas KDELR3, EIF4A1, and ATF3 were downregulated in OA rats tissues compared to the normal tissues, which were in accordance with our bioinformatics findings. Furthermore, our receiver operating characteristic (ROC) curve analysis verified that the above six ER stress-associated genes could be used as ideal biomarkers for OA diagnosis and those genes also potentially regulated immune responses by influencing the biological functions of mast cells and macrophages.

CONCLUSION

Collectively, the present study firstly identified six ER stress-associated genes (ATF3, DDX3X, JUN, EIF4A1, KDELR3, and VEGFA) that may play critical role in regulating the progression of OA.

Mitonuclear epistasis involving TP63 and haplogroup Uk: Risk of rapid progression of knee OA in patients from the OAI

OBJECTIVE

To investigate genetic interactions between mitochondrial deoxyribonucleic acid (mtDNA) haplogroups and nuclear single nucleotide polymorphisms (nSNPs) to analyze their impact on the development of the rapid progression of knee osteoarthritis (OA).

DESIGN

A total of 1095 subjects from the Osteoarthritis Initiative, with a follow-up time of at least 48-months, were included. Appropriate statistical approaches were performed, including generalized estimating equations adjusting for age, gender, body mass index, contralateral knee OA, Western Ontario and McMaster Universities Osteoarthritis Index pain, previous injury in target knee and the presence of the mtDNA variant m.16519C. Additional genomic data consisted in the genotyping of Caucasian mtDNA haplogroups and eight nSNPs previously associated with the risk of knee OA in robust genome-wide association studies.

RESULTS

The simultaneous presence of the G allele of rs12107036 at TP63 and the haplogroup Uk significantly increases the risk of a rapid progression of knee OA (odds ratio = 1.670; 95% confidence interval [CI]: 1.031-2.706; adjusted p-value = 0.027). The assessment of the population attributable fraction showed that the highest proportion of rapid progressors was under the simultaneous presence of the G allele of rs12107036 and the haplogroup Uk (23.4%) (95%CI: 7.89-38.9; p-value < 0.05). The area under the curve of the cross-validation model (0.730) was very similar to the obtained for the predictive model (0.735). A nomogram was constructed to help clinicians to perform clinical trials or epidemiologic studies.

CONCLUSIONS

This study demonstrates the existence of a mitonuclear epistasis in OA, providing new mechanisms by which nuclear and mitochondrial variation influence the susceptibility to develop different OA phenotypes.

Thyroid hormone induces ossification and terminal maturation in a preserved OA cartilage biomimetic model

OBJECTIVE

To characterize aspects of triiodothyronine (T3) induced chondrocyte terminal maturation within the molecular osteoarthritis pathophysiology using the previously established T3 human ex vivo osteochondral explant model.

DESIGNS

RNA-sequencing was performed on explant cartilage obtained from OA patients (n = 8), that was cultured ex vivo with or without T3 (10 ng/ml), and main findings were validated using RT-qPCR in an independent sample set (n = 22). Enrichment analysis was used for functional clustering and comparisons with available OA patient RNA-sequencing and GWAS datasets were used to establish relevance for OA pathophysiology by linking to OA patient genomic profiles.

RESULTS

Besides the upregulation of known hypertrophic genes EPAS1 and ANKH, T3 treatment resulted in differential expression of 247 genes with main pathways linked to extracellular matrix and ossification. CCDC80, CDON, ANKH and ATOH8 were among the genes found to consistently mark early, ongoing and terminal maturational OA processes in patients. Furthermore, among the 37 OA risk genes that were significantly affected in cartilage by T3 were COL12A1, TNC, SPARC and PAPPA.

CONCLUSIONS

RNA-sequencing results show that metabolic activation and recuperation of growth plate morphology are induced by T3 in OA chondrocytes, indicating terminal maturation is accelerated. The molecular mechanisms involved in hypertrophy were linked to all stages of OA pathophysiology and will be used to validate disease models for drug testing.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Causal Factors for Osteoarthritis: A Scoping Review of Mendelian Randomization Studies

OBJECTIVE

Mendelian randomization (MR) has increasingly been utilized as a tool for establishing causal relations between modifiable exposures and osteoarthritis (OA). The goal of this review was to summarize available MR studies of OA that evaluate the causal role of modifiable risk factors on OA.

METHODS

This review was performed following the Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) Extension for Scoping Reviews model. We performed a literature search for relevant studies published before December 2021 across multiple databases using the search terms "osteoarthritis" and ("Mendelian randomization" or "polygenic risk score"). We reported the MR estimates of causal associations between exposures and OA and then assessed methodologic quality of abstracted studies according to their efforts to validate the three key MR assumptions.

RESULTS

Our search identified 45 studies reporting on 141 exposure-association analyses. All studies performed a formal instrumental variable analysis to estimate the causal effect of exposure on OA. Causal associations (P < 0.05) were reported in 60 of these analyses representing 36 unique publications, and MR-Egger sensitivity analyses were performed in 45 of these analyses. MR studies provided support for causal associations of OA with increased levels of adiposity, coffee consumption, bone mineral density, and sleep disturbance, and decreased levels of serum calcium and low-density lipoprotein cholesterol.

CONCLUSION

These results highlight the potential benefits of weight reduction and improvement of sleep quality to reduce the risk of OA and call for a better understanding of the relations of coffee consumption and serum calcium to OA risk.

Polymorphism rs143384 GDF5 reduces the risk of knee osteoarthritis development in obese individuals and increases the disease risk in non-obese population

BACKGROUND

We investigated the effect of obesity on the association of genome-wide associative studies (GWAS)-significant genes with the risk of knee osteoarthritis (KOA).

METHODS

All study participants (n = 1,100) were divided into 2 groups in terms of body mass index (BMI): BMI ≥ 30 (255 KOA patients and 167 controls) and BMI < 30 (245 KOA and 433 controls). The eight GWAS-significant KOA single nucleotide polymorphisms (SNP) of six candidate genes, such as LYPLAL1 (rs2820436, rs2820443), SBNO1 (rs1060105, rs56116847), WWP2 (rs34195470), NFAT5 (rs6499244), TGFA (rs3771501), GDF5 (rs143384), were genotyped. Logistic regression analysis (gPLINK online program) was used for SNPs associations study with the risk of developing KOA into 2 groups (BMI ≥ 30 and BMI < 30) separately. The functional effects of KOA risk loci were evaluated using in silico bioinformatic analysis.

RESULTS

Multidirectional relationships of the rs143384 GDF5 with KOA in BMI-different groups were found: This SNP was KOA protective locus among individuals with BMI ≥ 30 (OR 0.41 [95%CI 0.20-0.94] recessive model) and was disorder risk locus among individuals with BMI < 30 (OR 1.32 [95%CI 1.05-1.65] allele model, OR 1.44 [95%CI 1.10-1.86] additive model, OR 1.67 [95%CI 1.10-2.52] dominant model). Polymorphism rs143384 GDF5 manifested its regulatory effects in relation to nine genes (GDF5, CPNE1, EDEM2, ERGIC3, GDF5OS, PROCR, RBM39, RPL36P4, UQCC1) in adipose tissue, which were involved in the regulation of pathways of apoptosis of striated muscle cells.

CONCLUSIONS

In summary, the effect of obesity on the association of the rs143384 GDF5 with KOA was shown: the "protective" value of this polymorphism in the BMI ≥ 30 group and the "risk" meaning in BMI < 30 cohort.

Epigenetics in Knee Osteoarthritis: A 2020-2023 Update Systematic Review

Osteoarthritis is a leading cause of disability in the world. The scientific literature highlights the critical importance of epigenetic regulatory effects, intertwined with biomechanical and biochemical peculiar conditions within each musculoskeletal district. While the contribution of genetic and epigenetic factors to knee OA is well-recognized, their precise role in disease management remains an area of active research. Such a field is particularly heterogeneous, calling for regular analysis and summarizing of the data that constantly emerge in the scientific literature, often sparse and scant of integration. The aim of this study was to systematically identify and synthesize all new evidence that emerged in human and animal model studies published between 2020 and 2023. This was necessary because, to the best of our knowledge, articles published before 2019 (and partly 2020) had already been included in systematic reviews that allowed to identify the ones concerning the knee joint. The review was carried out in accordance with Preferential Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. Only peer-reviewed articles were considered for inclusion. A total of 40 studies were identified, showing promising results in terms either of biomarker identification, new insight in mechanism of action or potential therapeutic targets for knee OA. DNA methylation, histone modification and ncRNA were all mechanisms involved in epigenetic regulation of the knee. Most recent evidence suggests that epigenetics is a most promising field with the long-term goal of improving understanding and management of knee OA, but a variety of research approaches need greater consolidation.

A landscape of gene expression regulation for synovium in arthritis

The synovium is an important component of any synovial joint and is the major target tissue of inflammatory arthritis. However, the multi-omics landscape of synovium required for functional inference is absent from large-scale resources. Here we integrate genomics with transcriptomics and chromatin accessibility features of human synovium in up to 245 arthritic patients, to characterize the landscape of genetic regulation on gene expression and the regulatory mechanisms mediating arthritic diseases predisposition. We identify 4765 independent primary and 616 secondary cis-expression quantitative trait loci (cis-eQTLs) in the synovium and find that the eQTLs with multiple independent signals have stronger effects and heritability than single independent eQTLs. Integration of genome-wide association studies (GWASs) and eQTLs identifies 84 arthritis related genes, revealing 38 novel genes which have not been reported by previous studies using eQTL data from the GTEx project or immune cells. We further develop a method called eQTac to identify variants that could affect gene expression by affecting chromatin accessibility and identify 1517 regions with potential regulatory function of chromatin accessibility. Altogether, our study provides a comprehensive synovium multi-omics resource for arthritic diseases and gains new insights into the regulation of gene expression.

The causal effect of two occupational factors on osteoarthritis and rheumatoid arthritis: a Mendelian randomization study

Background

Osteoarthritis (OA) and rheumatoid arthritis (RA) are two common types of arthritis. We conducted a two-sample Mendelian randomization (MR) study to estimate the causal effects of two common occupational factors-job involves heavy manual or physical work and job involves mainly walking or standing-on OA and RA in individuals of European ancestry.

Methods

Instruments were chosen from genome-wide association studies (GWASs) that identified independent single nucleotide polymorphisms (SNPs) robustly linked to job involves heavy manual or physical work (N = 263,615) as well as job involves mainly walking or standing (N = 263,556). Summary statistics for OA and RA were taken from the Integrative Epidemiology Unit (IEU) GWAS database; both discovery and replication GWAS datasets were considered. The primary analysis utilized the inverse variance weighted (IVW) MR method supplemented by various sensitivity MR analyses.

Results

In the IVW model, we found that genetically predicted job involves heavy manual or physical work was significantly associated with OA in both the discovery [odds ratio (OR) = 1.034, 95% confidence interval (CI): 1.016-1.053, P = 2.257 × 10-4] and replication (OR = 1.857, 95% CI: 1.223-2.822, P = 0.004) analyses. The causal associations were supported in diverse sensitivity analyses. MR analyses suggested no causal effect of genetically predicted job involves heavy manual or physical work on RA. Similarly, our data provided no evidence that genetically predicted job involves mainly walking or standing was related to OA and RA.

Conclusions

Our MR study suggests that job involves heavy manual or physical work is a risk factor for OA. It is of utmost importance to create preventive strategies aimed at reducing its impact on OA at such work sites.

Coffee and the risk of osteoarthritis: a two-sample, two-step multivariable Mendelian randomization study

Purpose: To investigate the potential causal relationship between coffee consumption and osteoarthritis (OA), and to disentangle whether body mass index (BMI) and Bone mineral density (BMD) mediate this relationship. Methods: We performed two-sample and two-step Mendelian randomization (MR) analyses utilizing publicly available genome-wide association studies (GWAS) summary statistics to estimate the association between coffee intake and OA risk (including knee OA, hip OA, knee or hip OA, and total OA), as well as the possible mediating effects of BMI and BMD. In addition, data of different coffee types (decaffeinated coffee, instant coffee, ground coffee-including espresso, filter, etc., and other coffee types) were used to explore the effect of coffee type on the risk of OA. Results: In two-sample MR, coffee intake increased the risk of OA in various sites, with the most significant impact observed in knee osteoarthritis (KOA) (odds ratio [OR] 2.03, 95% confidence interval [CI] 1.57-2.61, p < 0.001). The effect on self-reported OA was minimal (OR 1.03, 95% CI 1.01-1.05, p = 0.006). Further analysis of different types of coffee revealed that only decaffeinated coffee was causally associated with both KOA (OR 4.40, 95% CI 1.71-11.33, p = 0.002) and self-reported OA (OR 1.13, 95% CI 1.02-1.26, p = 0.022). In two-step MR, BMI explained over half of the coffee intake-all OA risk association, while BMD accounted for less than 5% of the mediation effect. Conclusion: Our study suggests that coffee intake increase the risk of OA, with BMI playing a significant mediating role. Decaffeinated coffee appears to have the greatest impact on OA risk compared to other types of coffee. Therefore, managing BMI and selecting appropriate types of coffee should be included in the health management of individuals who frequently consume coffee.

Genomic Determinants of Knee Joint Biomechanics: An Exploration into the Molecular Basis of Locomotor Function, a Narrative Review

In recent years, the nexus between genetics and biomechanics has garnered significant attention, elucidating the role of genomic determinants in shaping the biomechanical attributes of human joints, specifically the knee. This review seeks to provide a comprehensive exploration of the molecular basis underlying knee joint locomotor function. Leveraging advancements in genomic sequencing, we identified specific genetic markers and polymorphisms tied to key biomechanical features of the knee, such as ligament elasticity, meniscal resilience, and cartilage health. Particular attention was devoted to collagen genes like COL1A1 and COL5A1 and their influence on ligamentous strength and injury susceptibility. We further investigated the genetic underpinnings of knee osteoarthritis onset and progression, as well as the potential for personalized rehabilitation strategies tailored to an individual's genetic profile. We reviewed the impact of genetic factors on knee biomechanics and highlighted the importance of personalized orthopedic interventions. The results hold significant implications for injury prevention, treatment optimization, and the future of regenerative medicine, targeting not only knee joint health but joint health in general.

Calycosin ameliorates osteoarthritis by regulating the imbalance between chondrocyte synthesis and catabolism

Osteoarthritis (OA) is a severe chronic inflammatory disease. As the main active component of Astragalus mongholicus Bunge, a classic traditional ethnic herb, calycosin exhibits anti-inflammatory action and its mechanism of exact targets for OA have yet to be determined. In this study, we established an anterior cruciate ligament transection (ACLT) mouse model. Mice were randomized to sham, OA, and calycosin groups. Cartilage synthesis markers type II collagen (Col-2) and SRY-Box Transcription Factor 9 (Sox-9) increased significantly after calycosin gavage. While cartilage matrix degradation index cyclooxygenase-2 (COX-2), phosphor-epidermal growth factor receptor (p-EGFR), and matrix metalloproteinase-9 (MMP9) expression were decreased. With the help of network pharmacology and molecular docking, these results were confirmed in chondrocyte ADTC5 cells. Our results indicated that the calycosin treatment significantly improved cartilage damage, this was probably attributed to reversing the imbalance between chondrocyte synthesis and catabolism.

Investigating the causal relationship between physical activity and incident knee osteoarthritis: a two-sample Mendelian randomization study

There is evidence that physical activity (PA) has a long-term positive impact on disease. Whether PA is a risk factor for knee osteoarthritis (OA) is still controversial. The purpose of this study was to explore whether there is a causal relationship between PA and knee OA. We extracted PA and knee OA data from genome-wide association study (GWAS) databases. We used single-nucleotide polymorphisms (SNPs) as instrumental variables. We performed MR analysis by random-effects inverse-variance weighting (IVW), MR‒Egger, weighted median, simple mode, and weighted mode methods. We evaluated the stability and reliability of the results through sensitivity analysis. There was no significant association between PA and knee OA (p > 0.05). We did not detect any pleiotropy (MR‒Egger intercept test et al.: p > 0.05). The sensitivity analysis confirmed our results (p > 0.05). There is no causal relationship between PA and knee OA.

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

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

The causal effect of cytokine cycling levels on osteoarthritis: a bidirectional Mendelian randomized study

Objective

Osteoarthritis (OA) is the most prevalent joint disease globally, serving as a primary cause of pain and disability. However, the pathological processes underlying OA remain incompletely understood. Several studies have noted an association between cytokines and OA, yet the causal link between them remains ambiguous. This study aims to identify cytokines potentially causally related to OA using Mendelian randomization (MR) analysis, informing early clinical diagnosis and treatment decisions.

Methods

We conducted a genome-wide association study (GWAS) on 12 OA traits involving 177,517 cases and 649,173 controls from 9 international cohorts. For discovery MR analysis, we used 103 cytokines from two European populations as instrumental variables (IVs). Concurrently, another European population OA GWAS database (36,185 cases and 135,185 controls) was used to replicate MR analysis, employing the inverse variance weighted (IVW) method as the primary analytic approach. Additional methods tested included MR Egger, Weighted median, and Weighted mode. We merged the MR findings through meta-analysis. Heterogeneity testing, level pleiotropy testing (MR Egger intercept test and MRPRESSO), and sensitivity analysis via Leave One Out (LOO) were conducted to verify result robustness. Lastly, reverse MR analysis was performed.

Results

The meta-analysis merger revealed a correlation between CX3CL1 cycle levels and increased OA risk (OR=1.070, 95% CI: 1.040-1.110; P<0.010). We also observed associations between MCP4 (OR=0.930, 95% CI: 0.890-0.970; P<0.010) and CCL25 (OR=0.930, 95% CI: 0.890-0.970; P<0.010) with reduced OA risk. The sensitivity analysis results corroborate the robustness of these findings.

Conclusion

Our MR analysis indicates a potential causal relationship between CX3CL1, MCP4, CCL25, and OA risk changes. Further research is warranted to explore the influence of cytokines on OA development.

Evaluating the potential of Vitamin D and curcumin to alleviate inflammation and mitigate the progression of osteoarthritis through their effects on human chondrocytes: A proof-of-concept investigation

Osteoarthritis (OA) is a chronic degenerative joint disorder primarily affecting the elderly, characterized by a prominent inflammatory component. The long-term side effects associated with current therapeutic approaches necessitate the development of safer and more efficacious alternatives. Nutraceuticals, such as Vitamin D and curcumin, present promising therapeutic potentials due to their safety, efficacy, and cost-effectiveness. In this study, we utilized a proinflammatory human chondrocyte model of OA to assess the anti-inflammatory properties of Vitamin D and curcumin, with a particular focus on the Protease-Activated Receptor-2 (PAR-2) mediated inflammatory pathway. Employing a robust siRNA approach, we effectively modulated the expression of PAR-2 to understand its role in the inflammatory process. Our results reveal that both Vitamin D and curcumin attenuate the expression of PAR-2, leading to a reduction in the downstream proinflammatory cytokines, such as Tumor Necrosis Factor-alpha (TNF-α), Interleukin 6 (IL-6), and Interleukin 8 (IL-8), implicated in the OA pathogenesis. Concurrently, these compounds suppressed the expression of Receptor Activator of Nuclear Factor kappa-Β Ligand (RANKL) and its receptor RANK, which are associated with PAR-2 mediated TNF-α stimulation. Additionally, Vitamin D and curcumin downregulated the expression of Interferon gamma (IFN-γ), known to elevate RANKL levels, underscoring their potential therapeutic implications in OA. This study, for the first time, provides evidence of the mitigating effect of Vitamin D and curcumin on PAR-2 mediated inflammation, employing an siRNA approach in OA. Thus, our findings pave the way for future research and the development of novel, safer, and more effective therapeutic strategies for managing OA.

Genetic liability for diet-derived circulating antioxidants, oxidative stress, and risk of osteoarthritis: a Mendelian randomization study

Background

Although well-documented, the causal relationships between diet-derived circulating antioxidants, oxidative stress, and osteoarthritis (OA) are equivocal. The objective of this study is to employ two-sample Mendelian randomization (MR) to investigate possible causal relationships among dietary-derived circulating antioxidants, oxidative stress damage indicators, and OA risk.

Methods

Single-nucleotide polymorphisms for diet-derived circulating antioxidants (ascorbate, β-carotene, lycopene, retinol, and α-and γ-tocopherol), assessed as absolute levels and metabolites, as well as oxidative stress injury biomarkers (GSH, GPX, CAT, SOD, albumin, and total bilirubin), were retrieved from the published data and were used as genetic instrumental variables. Summary statistics for gene-OA associations were obtained from publicly available and two relatively large-scale GWAS meta-analyses to date. The inverse-variance weighting method was utilized as the primary MR analysis. Moreover, multivariable MR was used to determine if mediators (BMI and smoking) causally mediated any connection. Furthermore, for each exposure, MR analyses were conducted per outcome database and then meta-analyzed.

Results

Genetically predicted absolute retinol level was causally associated with hip OA risk [odds ratios (ORs) = 0.40, 95% confidence interval (CI) = 0.24-0.68, FDR-corrected p = 0.042]. Moreover, genetically predicted albumin level was causally associated with total OA risk (OR = 0.80, 95% CI = 0.75-0.86, FDR-corrected p = 2.20E-11), as well as the risk of hip OA (OR = 0.75, 95% CI = 0.68-0.84, FDR-corrected p = 1.38E-06) and knee OA (OR = 0.82, 95% CI = 0.76-0.89, FDR-corrected p = 4.49E-06). In addition, MVMR confirmed that the effect of albumin on hip OA is independent of smoking initiation, alcoholic drinks per week, and moderate-to-vigorous physical activity levels but may be influenced by BMI.

Conclusion

Evidence from our study supports a potentially protective effect of high levels of retinol and albumin on OA risk.

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.

Age-dependent genetic regulation of osteoarthritis: independent effects of immune system genes

OBJECTIVES

Osteoarthritis (OA) is a joint disease with a heritable component. Genetic loci identified via genome-wide association studies (GWAS) account for an estimated 26.3% of the disease trait variance in humans. Currently, there is no method for predicting the onset or progression of OA. We describe the first use of the Collaborative Cross (CC), a powerful genetic resource, to investigate knee OA in mice, with follow-up targeted multi-omics analysis of homologous regions of the human genome.

METHODS

We histologically screened 275 mice for knee OA and conducted quantitative trait locus (QTL) mapping in the complete cohort (> 8 months) and the younger onset sub-cohort (8-12 months). Multi-omic analysis of human genetic datasets was conducted to investigate significant loci.

RESULTS

We observed a range of OA phenotypes. QTL mapping identified a genome-wide significant locus on mouse chromosome 19 containing Glis3, the human equivalent of which has been identified as associated with OA in recent GWAS. Mapping the younger onset sub-cohort identified a genome-wide significant locus on chromosome 17. Multi-omic analysis of the homologous region of the human genome (6p21.32) indicated the presence of pleiotropic effects on the expression of the HLA - DPB2 gene and knee OA development risk, potentially mediated through the effects on DNA methylation.

CONCLUSIONS

The significant associations at the 6p21.32 locus in human datasets highlight the value of the CC model of spontaneous OA that we have developed and lend support for an immune role in the disease. Our results in mice also add to the accumulating evidence of a role for Glis3 in OA.

Cam morphology is strongly and consistently associated with development of radiographic hip osteoarthritis throughout 4 follow-up visits within 10 years

OBJECTIVE

To determine the association between cam morphology and the development of radiographic hip osteoarthritis (RHOA) at four time points within 10-year follow-up.

DESIGN

The nationwide prospective Cohort Hip and Cohort Knee study includes 1002 participants aged 45-65 years with 2-, 5-, 8-, and 10-year follow-ups. The associations of cam morphology (alpha angle >60°) and large cam morphology (alpha angle >78°) in hips free of osteoarthritis at baseline (Kellgren & Lawrence (KL) grade <2) with the development of both incident RHOA (KL grade≥2) and end-stage RHOA (KL grade≥3) were estimated using logistic regression with generalized estimating equation at each follow-up and using Cox regression over 10 years, adjusted for age, sex, and body mass index.

RESULTS

Both cam morphology and large cam morphology were associated with the development of incident RHOA at all follow-ups with adjusted Odd Ratios (aORs) ranging from 2.7 (95% Confidence interval 1.8-4.1) to 2.9 (95% CI 2.0-4.4) for cam morphology and ranging from 2.5 (95% CI 1.5-4.3) to 4.2 (95% CI 2.2-8.3) for large cam morphology. For end-stage RHOA, cam morphology resulted in aORs ranging from 4.9 (95% CI 1.8-13.2) to 8.5 (95% CI 1.1-64.4), and aORs for large cam morphology ranged from 6.7 (95% CI 3.1-14.7) to 12.7 (95% CI 1.9-84.4).

CONCLUSIONS

Cam morphology poses the hip at 2-13 times increased odds for developing RHOA within a 10-year follow-up. The association was particularly strong for large cam morphology and end-stage RHOA, while the strength of association was consistent over time.

The effect of immunomodulatory drugs on aortic stenosis: a Mendelian randomisation analysis

There are currently no approved pharmacological treatment options for aortic stenosis (AS), and there are limited identified drug targets for this chronic condition. It remains unclear whether inflammation plays a role in AS pathogenesis and whether immunomodulation could become a therapeutic target. We evaluated the potentially causal association between inflammation and AS by investigating the genetically proxied effects of tocilizumab (IL6 receptor, IL6R, inhibitor), canakinumab (IL1β inhibitor) and colchicine (β-tubulin inhibitor) through a Mendelian randomisation (MR) approach. Genetic proxies for these drugs were identified as single nucleotide polymorphisms (SNPs) in the gene, enhancer or promoter regions of IL6R, IL1β or β-tubulin gene isoforms, respectively, that were significantly associated with serum C-reactive protein (CRP) in a large European genome-wide association study (GWAS; 575,531 participants). These were paired with summary statistics from a large GWAS of AS in European patients (653,867 participants) to then perform primary inverse-variance weighted random effect and sensitivity MR analyses for each exposure. This analysis showed that genetically proxied tocilizumab was associated with reduced risk of AS (OR 0.56, 95% CI 0.45-0.70 per unit decrease in genetically predicted log-transformed CRP). Genetically proxied canakinumab was not associated with risk of AS (OR 0.80, 95% CI 0.51-1.26), and only one suitable SNP was identified to proxy the effect of colchicine (OR 34.37, 95% CI 1.99-592.89). The finding that genetically proxied tocilizumab was associated with reduced risk of AS is concordant with an inflammatory hypothesis of AS pathogenesis. Inhibition of IL6R may be a promising therapeutic target for AS management.

Crystal arthropathies and osteoarthritis-where is the link?

Osteoarthritis (OA) is one of the leading causes of disability worldwide. As our understanding of OA progressively has moved from a purely mechanical "wear and tear" concept toward a complex multi-tissue condition in which inflammation plays a central role, the possible role of crystal-induced inflammation in OA incidence and progression may be relevant. In addition to gout, which affects 4% of the US population, basic calcium phosphate and calcium pyrophosphate deposition both may induce joint inflammation and may play a role in pain in OA. This narrative review article discusses the possible mechanisms underlying the associations between crystal-induced arthropathies and OA, and the important implications of these for clinical practice and future research.

Causal effects of walking pace on osteoarthritis: a two-sample mendelian randomization study

Background: Osteoarthritis (OA) is one of the most common joint diseases worldwide, imposing a substantial burden on individuals and society. Numerous pieces of evidence suggest that walking pace (WP) can serve as a predictive indicator for the risk of various diseases, and observational studies have also found a potential link between WP and the risk of OA. However, the causal relationship between WP and the risk of OA remains unclear. Methods: We conducted a mendelian randomization (MR) study using data from the European Genome-wide Association Study, which included WP (including 459,915 participants), OA (including 10,083 cases and 40,425 controls), knee OA (including 24,955 cases and 378,169 controls), and hip OA (including 15,704 cases and 378,169 controls). Single nucleotide polymorphisms (SNPs) associated with WP were utilized to infer causal associations with OA and its subtypes. The Inverse Variance Weighted (IVW) technique served as the primary causal analysis method. Three auxiliary MR methods - MR-Egger, weighted median, and maximum likelihood - were used to substantiate the IVW results. Sensitivity analyses were performed to examine heterogeneity and pleiotropy. In addition, multivariate MR (MVMR) analysis was used to assess causality after adjustment for three potential confounders. Results: According to the results of the IVW method, every 1 standard deviation increased in genetic WP corresponds to an 89% reduction in the risk of OA (odds ratio (OR) = 0.11; 95% confidence interval (CI), 0. 06-0.19; p = 1.57 × 10-13), an 83% reduction in the risk of knee OA (OR = 0.17; 95% CI, 0.11-0.28; p = 2.78 × 10-13), and a 76% reduction in the risk of hip OA (OR = 0.24; 95% CI, 0.14-0.43; p = 1.51 × 10-6). These results were confirmed by the three additional MR methods and validated by the sensitivity analysis. Ultimately, the MVMR analysis confirmed that the role of WP in reducing the risk of OA and its subtypes remains consistent regardless of potential confounders. Conclusion: The results of our MR study highlight a significant causal association between WP and the susceptibility to OA, including its knee and hip subtypes. These findings propose that WP could be utilized as a potential prognostic factor for OA risk.

Genetic Risk Factors for End-Stage Hip Osteoarthritis Treated With Total Hip Arthroplasty: A Genome-wide Association Study

BACKGROUND

Although a genetic component to hip osteoarthritis (OA) has been described, focused evaluation of the genetic components of end-stage disease is limited. We present a genomewide association study for patients undergoing total hip arthroplasty (THA) to characterize the genetic risk factors associated with end-stage hip osteoarthritis (ESHO), defined as utilization of the procedure.

METHODS

Patients who underwent primary THA for hip OA were identified in a national patient data repository using administrative codes. Fifteen thousand three hundred and fifty-five patients with ESHO and 374,193 control patients were identified. Whole genome regression of genotypic data for patients who underwent primary THA for hip OA corrected for age, sex, and body mass index (BMI) was performed. Multivariate logistic regression models were used to evaluate the composite genetic risk from the identified genetic variants.

RESULTS

There were 13 significant genes identified. Composite genetic factors resulted in an odds ratio 1.04 for ESHO (P < .001). The effect of genetics was lower than that of age (Odds Ratio (OR): 2.38; P < .001) and BMI (1.81; P < .001).

CONCLUSION

Multiple genetic variants, including 5 novel loci, were associated with end-stage hip OA treated with primary THA. Age and BMI were associated with greater odds of developing end-stage disease when compared to genetic factors.

Lifetime risk and genetic predisposition to post-traumatic OA of the knee in the UK Biobank

OBJECTIVE

Acute knee injury is associated with post-traumatic OA (PTOA). Very little is known about the genome-wide associations of PTOA when compared with idiopathic OA (iOA). Our objective was to describe the development of knee OA after a knee injury and its genetic associations in UK Biobank (UKB).

DESIGN

Clinically significant structural knee injuries in those ≤50 years were identified from electronic health records and self-reported data in 502,409 UKB participants. Time-to-first knee osteoarthritis (OA) code was compared in injured cases and age-/sex-matched non-injured controls using Cox Proportional Hazards models. A time-to-OA genome-wide association study (GWAS) sought evidence for PTOA risk variants 6 months to 20 years following injury. Evidence for associations of two iOA polygenic risk scores (PRS) was sought.

RESULTS

Of 4233 knee injury cases, 1896 (44.8%) were female (mean age at injury 34.1 years [SD 10.4]). Over a median of 30.2 (IQR 19.5-45.4) years, 1096 (25.9%) of injured cases developed knee OA. The overall hazards ratio (HR) for knee OA after injury was 1.81 (1.70,1.93), P = 8.9 × 10-74. Female sex and increasing age at injury were associated with knee OA following injury (HR 1.15 [1.02,1.30];1.07 [1,07,1.07] respectively). OA risk was highest in the first 5 years after injury (HR 3.26 [2.67,3.98]), persisting for 40 years. In 3074 knee injury cases included in the time-to-OA GWAS, no variants reached genome-wide significance. iOA PRS was not associated with time-to-OA (HR 0.43 [0.02,8.41]).

CONCLUSIONS

Increasing age at injury and female sex appear to be associated with future development of PTOA in UKB, the risk of which was greatest in the 5 years after injury. Further international efforts towards a better-powered meta-analysis will definitively elucidate genetic similarities and differences of PTOA and iOA.

Identification of candidate genes and chemicals associated with osteoarthritis by transcriptome-wide association study and chemical-gene interaction analysis

BACKGROUND

Osteoarthritis (OA) is a common degenerative joint disease and causes chronic pain and disability to the elderly. Several risk factors are involved, such as aging, obesity, genetic susceptibility, and environmental factors. We conducted a transcriptome-wide association study (TWAS) and chemical-related gene set enrichment analysis (CGSEA) to investigate the susceptibility genes and environmental factors.

METHODS

TWAS analysis was conducted to identify the susceptibility genes by integrating the summary-level genome-wide association study data of knee OA (KOA) and hip OA (HOA) with the precomputed expression weights from the Genotype-Tissue Expression Project (Version 8). The FUSION software was used for both single-tissue and cross-tissue TWAS, which were combined using an aggregate Cauchy association test. The biological function and pathways of the TWAS genes were explored using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases, and the human cartilage mRNA expression profiles were utilized to validate the TWAS genes. CGSEA analysis was performed to scan the OA-associated chemicals by integrating the TWAS results with the chemical-related gene sets.

RESULTS

There were 44 and 93 unique TWAS genes identified in 7 and 11 chromosomes for KOA and HOA, respectively, fourteen and four of which showed significantly differential expression in the mRNA profiles, such as CRHR1, LTBP1, WWP2, LMX1B, and PTHLH. OA-related pathways were found in the KEGG and GO analysis, such as TGF-beta signaling pathway, MAPK signaling pathway, hyaluronan metabolic process, and chondrocyte differentiation. Forty-five OA-associated chemicals were identified, including quercetin, bisphenol A, and cadmium chloride.

CONCLUSIONS

Several candidate OA-associated genes and chemicals were identified through TWAS and CGSEA analysis, which expanded our understanding of the relationship between genes, chemicals, and their impact on OA.

The identification of distinct protective and susceptibility mechanisms for hip osteoarthritis: findings from a genome-wide association study meta-analysis of minimum joint space width and Mendelian randomisation cluster analyses

BACKGROUND

Hip minimum joint space width (mJSW) provides a proxy for cartilage thickness. This study aimed to conduct a genome-wide association study (GWAS) of mJSW to (i) identify new genetic determinants of mJSW and (ii) identify which mJSW loci convey hip osteoarthritis (HOA) risk and would therefore be of therapeutic interest.

METHODS

GWAS meta-analysis of hip mJSW derived from plain X-rays and DXA was performed, stratified by sex and adjusted for age and ancestry principal components. Mendelian randomisation (MR) and cluster analyses were used to examine causal effect of mJSW on HOA.

FINDINGS

50,745 individuals were included in the meta-analysis. 42 SNPs, which mapped to 39 loci, were identified. Mendelian randomisation (MR) revealed little evidence of a causal effect of mJSW on HOA (ORIVW 0.98 [95% CI 0.82-1.18]). However, MR-Clust analysis suggested the null MR estimates reflected the net effect of two distinct causal mechanisms cancelling each other out, one of which was protective, whereas the other increased HOA susceptibility. For the latter mechanism, all loci were positively associated with height, suggesting mechanisms leading to greater height and mJSW increase the risk of HOA in later life.

INTERPRETATIONS

One group of mJSW loci reduce HOA risk via increased mJSW, suggesting possible utility as targets for chondroprotective therapies. The second group of mJSW loci increased HOA risk, despite increasing mJSW, but were also positively related to height, suggesting they contribute to HOA risk via a growth-related mechanism.

FUNDING

Primarily funded by the Medical Research Council and Wellcome Trust.

The association between frailty and osteoarthritis based on the NHANES and Mendelian randomization study

Introduction

The aim of the study was to examine the association between frailty and osteoarthritis.

Methods

We conducted a cross-sectional study in the National Health and Nutrition Examination Survey, while logistic regression was used to explore the association of the two. Mendelian randomization (MR) study was used to explore the causal relationship between the two.

Results

In the cross-sectional study, logistic regression analysis showed that odds ratio (OR) (95% CI) value was 1.07 (1.05, 1.08). In the MR study, the inverse-variance weighting (IVW) results showed OR (95% CI) value of 1.69 (1.01-2.83).

Conclusions

There is both a correlation and a causal relationship between frailty and osteoarthritis, and frailty may be a potentially better response than age to osteoarthritis.

Genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis

Multimorbidity is a rising public health challenge with important implications for health management and policy. The most common multimorbidity pattern is the combination of cardiometabolic and osteoarticular diseases. Here, we study the genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis. We find genome-wide genetic correlation between the two diseases and robust evidence for association-signal colocalization at 18 genomic regions. We integrate multi-omics and functional information to resolve the colocalizing signals and identify high-confidence effector genes, including FTO and IRX3, which provide proof-of-concept insights into the epidemiologic link between obesity and both diseases. We find enrichment for lipid metabolism and skeletal formation pathways for signals underpinning the knee and hip osteoarthritis comorbidities with type 2 diabetes, respectively. Causal inference analysis identifies complex effects of tissue-specific gene expression on comorbidity outcomes. Our findings provide insights into the biological basis for the type 2 diabetes-osteoarthritis disease co-occurrence.

Down-regulation of Jun induces senescence through destabilizing chromatin in osteoarthritis chondrocytes

OBJECTIVE

Osteoarthritis (OA) is the most common degenerative joint disease leading to disability worldwide. Cellular senescence is considered to be a fundamental pathogenic mechanism in the development of OA and has attracted increasing attention. However, regulatory mechanisms underlying chondrocyte senescence in OA remain unclear.

METHODS

Bioinformatic methods were used to screen key genes. Immunohistochemistry and the quantitative reverse transcription polymerase chain reaction were used to evaluate gene expression. RNA intervention experiments were performed to explore the functions of key genes.

RESULTS

We used 494 aging-associated genes provided by the Aging Atlas to identify the co-expression modules associated with age and OA. Thirty age-associated differentially expressed genes (ASDEGs) were identified. Using cytoHubba in Cytoscape, we identified Jun as the hub-ASDEG for OA chondrocytes. We confirmed the downregulation of Jun in OA rats and senescent chondrocytes by immunohistochemistry and quantitative reverse transcription polymerase chain reaction, respectively. Inhibition of proliferation and accelerated senescence were observed in chondrocytes treated with siRNA against Jun. Mechanistically, we observed micronuclei formation and reduced expression of H3K9me3 and heterochromatin protein 1gamma in siRNA-Jun-treated chondrocytes, indicating that destabilization of chromatin occurred during this treatment.

CONCLUSION

Jun plays a crucial role in OA development and causes senescence by destabilizing chromatin in chondrocytes. These findings provide new insights into OA progression and suggest promising therapeutic targets.

Sox, Fox, and Lmx1b binding sites differentially regulate a Gdf5-Associated regulatory region during elbow development

Introduction: The articulating ends of limb bones have precise morphology and asymmetry that ensures proper joint function. Growth differentiation factor 5 (Gdf5) is a secreted morphogen involved in cartilage and bone development that contributes to the architecture of developing joints. Dysregulation of Gdf5 results in joint dysmorphogenesis often leading to progressive joint degeneration or osteoarthritis (OA). The transcription factors and cis-regulatory modules (CRMs) that regulate Gdf5 expression are not well characterized. We previously identified a Gdf5-associated regulatory region (GARR) that contains predicted binding sites for Lmx1b, Osr2, Fox, and the Sox transcription factors. These transcription factors are recognized factors involved in joint morphogenesis and skeletal development. Methods: We used in situ hybridization to Gdf5, Col2A1, and the transcription factors of interest in developing chicken limbs to determine potential overlap in expression. We further analyzed scRNA-seq data derived from limbs and knees in published mouse and chicken datasets, identifying cells with coexpression of Gdf5 and the transcription factors of interest. We also performed site-directed mutatgenesis of the predicted transcription factor binding sites in a GARR-reporter construct and determined any change in activity using targeted regional electroporation (TREP) in micromass and embryonic chicken wing bioassays. Results: Gdf5 expression overlapped the expression of these transcription factors during joint development both by in situ hybridization (ISH) and scRNA-seq analyses. Within the GARR CRM, mutation of two binding sites common to Fox and Sox transcripstion factors reduced enhancer activity to background levels in micromass cultures and in ovo embryonic chicken wing bioassays, whereas mutation of two Sox-only binding sites caused a significant increase in activity. These results indicate that the Fox/Sox binding sites are required for activity, while the Sox-only sites are involved in repression of activity. Mutation of Lmx1b binding sites in GARR caused an overall reduction in enhancer activity in vitro and a dorsal reduction in ovo. Despite a recognized role for Osr2 in joint development, disruption of the predicted Osr2 site did not alter GARR activity. Conclusion: Taken together, our data indicates that GARR integrates positive, repressive, and asymmetrical inputs to fine-tune the expression of Gdf5 during elbow joint development.

A shared genetic contribution to osteoarthritis and COVID-19 outcomes: a large-scale genome-wide cross-trait analysis

Background

Patients with osteoarthritis (OA) are exposed to an increased risk of adverse outcomes of COVID-19, and they tend to experience disruption in access to healthcare services and exercise facilities. However, a deep understanding of this comorbidity phenomenon and the underlying genetic architecture of the two diseases is still unclear. In this study, we aimed to untangle the relationship between OA and COVID-19 outcomes by conducting a large-scale genome-wide cross-trait analysis.

Methods

Genetic correlation and causal relationships between OA and COVID-19 outcomes (critical COVID-19, COVID-19 hospitalization, and COVID-19 infection) were estimated by linkage disequilibrium score regression and Mendelian Randomization approaches. We further applied Multi-Trait Analysis of GWAS and colocalization analysis to identify putative functional genes associated with both OA and COVID-19 outcomes.

Results

Significant positive genetic correlations between OA susceptibility and both critical COVID-19 (rg=0.266, P=0.0097) and COVID-19 hospitalization (rg=0.361, P=0.0006) were detected. However, there was no evidence to support causal genetic relationships between OA and critical COVID-19 (OR=1.17[1.00-1.36], P=0.049) or OA and COVID-19 hospitalization OR=1.08[0.97-1.20], P=0.143). These results were robustly consistent after the removal of obesity-related single nucleotide polymorphisms (SNPs). Moreover, we identified a strong association signal located near the FYCO1 gene (lead SNPs: rs71325101 for critical COVID-19, Pmeta=1.02×10-34; rs13079478 for COVID-19 hospitalization, Pmeta=1.09×10-25).

Conclusion

Our findings further confirmed the comorbidity of OA and COVID-19 severity, but indicate a non-causal impact of OA on COVID-19 outcomes. The study offers an instructive perspective that OA patients did not generate negative COVID-19 outcomes during the pandemic in a causal way. Further clinical guidance can be formulated to enhance the quality of self-management in vulnerable OA patients.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

TRIAL REGISTRATION

No. C278.

Automatic block-wise genotype-phenotype association detection based on hidden Markov model

BACKGROUND

For detecting genotype-phenotype association from case-control single nucleotide polymorphism (SNP) data, one class of methods relies on testing each genomic variant site individually. However, this approach ignores the tendency for associated variant sites to be spatially clustered instead of uniformly distributed along the genome. Therefore, a more recent class of methods looks for blocks of influential variant sites. Unfortunately, existing such methods either assume prior knowledge of the blocks, or rely on ad hoc moving windows. A principled method is needed to automatically detect genomic variant blocks which are associated with the phenotype.

RESULTS

In this paper, we introduce an automatic block-wise Genome-Wide Association Study (GWAS) method based on Hidden Markov model. Using case-control SNP data as input, our method detects the number of blocks associated with the phenotype and the locations of the blocks. Correspondingly, the minor allele of each variate site will be classified as having negative influence, no influence or positive influence on the phenotype. We evaluated our method using both datasets simulated from our model and datasets from a block model different from ours, and compared the performance with other methods. These included both simple methods based on the Fisher's exact test, applied site-by-site, as well as more complex methods built into the recent Zoom-Focus Algorithm. Across all simulations, our method consistently outperformed the comparisons.

CONCLUSIONS

With its demonstrated better performance, we expect our algorithm for detecting influential variant sites may help find more accurate signals across a wide range of case-control GWAS.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Genetically predicted circulating levels of cytokines and the risk of osteoarthritis: A mendelian randomization study

Background: The association between inflammatory cytokines and osteoarthritis (OA) has been reported in several observational studies, but the causal relationship between these two remains unknown. Hence, we performed this two-sample Mendelian randomization (MR) to confirm the causal relationship between circulating levels of inflammatory factors and osteoarthritis risk.

Method: We used genetic variants associated with cytokine circulation levels from a meta-analysis of genome-wide association studies (GWASs) in 8,293 Finns as instrumental variables and obtained OA data from the United Kingdom Biobank, including a total of 345,169 subjects of European ancestry (66,031 diagnosed OA cases and 279,138 controls). Inverse variance weighting (IVW), MR-Egger, Wald Ratio, weighted median, and MR multiplicity residual sums with outliers (MR-PRESSO) were used.

Result: We found a causal relationship between circulating levels of macrophage inflammatory protein-1beta (MIP-1β) and risk of OA (OR = 0.998, 95% CI = 0.996–0.999p = 9.61 × 10−5); tumour necrosis factor beta (TNF-β) was also causally associated with risk of OA (OR = 0.996,95%CI = 0.994–0.999, p = 0.002); finally we found a suggestive association between C-C motif chemokine ligand 5(CCL5, also called Rantes) and OA risk (OR = 1.013, 95%CI = 1.002–1.024,p = 0.016).

Conclusion: Our findings offer promising leads for the development of new therapeutic targets in the treatment of osteoarthritis. By identifying the role of inflammatory cytokines in this debilitating condition through a genetic epidemiological approach, our study contributes to a better understanding of the underlying disease mechanisms. These insights may ultimately pave the way for more effective treatments that improve patient outcomes.

A sex- and site-specific relationship between body mass index and osteoarthritis: evidence from observational and genetic analyses

OBJECTIVE

We primarily aimed to investigate whether there are phenotypic and genetic links underlying body mass index (BMI) and overall osteoarthritis (OA). We then intended to explore whether the relationships differ across sexes and sites.

METHOD

We first evaluated the phenotypic association between BMI and overall OA using data from the UK Biobank. We then investigated the genetic relationship leveraging summary statistics of the hitherto largest genome-wide association studies performed for BMI and overall OA. Finally, we repeated all analyses in a sex- (female, male) and site- (knee, hip, spine) specific manner.

RESULTS

Observational analysis suggested an increased hazard of diagnosed OA per 5 kg/m² increment in BMI (hazard ratio = 1.38, 95% confidence interval (CI) = 1.37-1.39). A positive overall genetic correlation was observed for BMI and OA (r_g = 0.43, P = 4.72 × 10^-133), corroborated by 11 significant local signals. Cross-trait meta-analysis identified 34 pleiotropic loci shared between BMI and OA, of which seven were novel. Transcriptome-wide association study revealed 29 shared gene-tissue pairs, targeting nervous, digestive, and exo/endocrine systems. Mendelian randomization demonstrated a robust BMI-OA causal relationship (odds ratio = 1.47, 95% CI = 1.42-1.52). A similar pattern of effects was observed in sex- and site-specific analyses, with BMI affecting OA comparably in both sexes and most strongly in the knee.

CONCLUSION

Our work demonstrates an intrinsic relationship underlying BMI and overall OA, reflected by a pronounced phenotypic association, significant biological pleiotropy, and a putative causal link. Stratified analysis further reveals that the effects are distinct across sites and comparable across sexes.

Association of Lipid-Lowering Drugs With Risk of Psoriasis: A Mendelian Randomization Study

Importance

Lipid pathways have been implicated in the pathogenesis of psoriasis, and some lipid-lowering drugs, such as statins, are hypothesized to have disease-modifying properties. However, large population-level studies are scarce, and causal interpretation of results from traditional observational designs is limited by confounding.

Objective

To investigate the causal association between genetically proxied lipid-lowering drugs and psoriasis risk.

Design, Setting, and Participants

This 2-sample mendelian randomization study was performed from August to October 2022 and included population-level genome-wide association studies of psoriasis in the UK Biobank and FinnGen studies and low-density lipoprotein (LDL) by the Global Lipids Genetics Consortium. The inverse variance-weighted method was used with pleiotropy robust methods and colocalization as sensitivity analyses.

Exposures

Genetically proxied inhibition of 3-hydroxy-3-methylglutaryl CoA reductase (HMGCR, targeted by statins), Niemann-Pick C1-like 1 (NPC1L1, targeted by ezetimibe), and proprotein convertase subtilisin/kexin type 9 (PCSK9, targeted by, eg, alirocumab), using LDL as the biomarker.

Main Outcomes and Measures

Risk of psoriasis.

Results

Data from 12 116 patients with psoriasis and approximately 1.3 million individuals with LDL measurement were analyzed. Genetically proxied PCSK9 inhibition was associated with reduced risk of psoriasis (odds ratio, 0.69 per standard deviation reduction in LDL; 95% CI, 0.55-0.88; P = .003), which was replicated in FinnGen (odds ratio, 0.71; 95% CI, 0.57-0.88; P = .002). Sensitivity analyses did not provide statistical evidence of bias from pleiotropy or genetic confounding. No robust association was found for HMGCR or NPC1L1 inhibition.

Conclusions and Relevance

This mendelian randomization study suggests that PCSK9 is implicated in psoriasis pathogenesis, and its inhibition is associated with reduced psoriasis risk. These findings potentially pave the way for future studies that may allow personalized selection of lipid-lowering drugs for those at risk of psoriasis.

Gsmtx4 Alleviated Osteoarthritis through Piezo1/Calcineurin/NFAT1 Signaling Axis under Excessive Mechanical Strain

Excessive mechanical strain is the prominent risk factor for osteoarthritis (OA), causing cartilage destruction and degeneration. However, the underlying molecular mechanism contributing to mechanical signaling transduction remains unclear in OA. Piezo type mechanosensitive ion channel component 1 (Piezo1) is a calcium-permeable mechanosensitive ion channel and provides mechanosensitivity to cells, but its role in OA development has not been determined. Herein, we found up-regulated expression of Piezo1 in OA cartilage, and that its activation contributes to chondrocyte apoptosis. The knockdown of Piezo1 could protect chondrocytes from apoptosis and maintain the catabolic and anabolic balance under mechanical strain. In vivo, Gsmtx4, a Piezo1 inhibitor, markedly ameliorated the progression of OA, inhibited the chondrocyte apoptosis, and accelerated the production of the cartilage matrix. Mechanistically, we observed the elevated activity of calcineurin (CaN) and the nuclear transfection of nuclear factor of activated T cells 1 (NFAT1) under mechanical strain in chondrocytes. Inhibitors of CaN or NFAT1 rescued the pathologic changes induced by mechanical strain in chondrocytes. Overall, our findings revealed that Piezo1 was the essential molecule response to mechanical signals and regulated apoptosis and cartilage matrix metabolism via the CaN/NFAT1 signaling axis in chondrocytes, and that Gsmtx4 could be an attractive therapeutic drug for OA treatment.

Genetically predicted obesity and risk of hip osteoarthritis

OBJECTIVES

To determine the causal association between genetically predicted obesity and the risk of hip osteoarthritis.

METHODS

We performed two-sample Mendelian randomization (MR) analysis to analyze the association between body mass index (BMI) and hip osteoarthritis using pooled-level genome-wide association study (GWAS) data. The inverse variance weighted (IVW), MR‒Egger, and weighted median methods were used to estimate the causal association. In addition, we applied the MR Steiger filtering method, MR robust adjusted profile score (MR.RAPS) methods, and the MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO) global test to examine and address potential horizontal pleiotropy.

RESULTS

We found a causal relationship between genetically predicted BMI and the risk of hip osteoarthritis by the IVW method [OR = 1.45, 95% confidence interval (CI) = 1.04-2.00, P = 0.02]. In the sensitivity analysis, the results of the MR‒Egger and weighted median methods revealed similar estimations but with a wide CI with lower precision. The funnel plot, MR-Egger intercept, and MR-PRESSO all indicated the absence of a directional pleiotropic effect. In addition, no heterogeneity was observed in the present analysis. Therefore, the result of IVW is most suitable and reliable for the present MR analysis.

CONCLUSION

There is a causal relationship between obesity and a higher risk of hip osteoarthritis, suggesting that weight management may be an intervention for the prevention and management of hip osteoarthritis.

LEVEL OF EVIDENCE

Bioinformatics, Basic science.