Stratification of knee osteoarthritis: two major patient subgroups identified by genome-wide expression analysis of articular cartilage

Deletion of Bmal1 in aggrecan-expressing cells leads to mouse temporomandibular joint osteoarthritis

INTRODUCTION

Articular cartilage is the major affected tissue during the development of osteoarthritis (OA) in temporomandibular joint (TMJ). The core circadian rhythm molecule Bmal1 regulates chondrocyte proliferation, differentiation and apoptosis; however, its roles in condylar cartilage function and in TMJ OA have not been fully elucidated.

MATERIALS AND METHODS

TMJ OA mouse model was induced by unilateral anterior crossbite (UAC) and Bmal1 protein expression in condylar cartilage were examined by western blot analysis. To determine the role of Bmal1 in TMJ OA, we generated cartilage-specific Bmal1 conditional knockout (cKO) mice (Bmal1Agc1CreER mice) and hematoxylin and eosin staining, toluidine blue and Safranin O/fast green, immunohistochemistry, TUNEL assay, real-time PCR analysis and Western blot assay were followed.

RESULTS

Bmal1 expression was reduced in condylar cartilage in a TMJ OA mouse model induced by UAC. The Bmal1 cKO mice displayed decreased cartilage matrix synthesis, reduced chondrocyte proliferation, increased chondrocyte hypertrophy and apoptosis as well as the upregulation of YAP expression in TMJ condylar cartilage.

CONCLUSIONS

We demonstrated that Bmal1 was essential for TMJ tissue homeostasis and loss-of-function of Bmal1 in chondrocytes leads to the development of TMJ OA.

DSCR1-1 attenuates osteoarthritis-associated chondrocyte injury by regulating the CREB1/ALDH2/Wnt/β-catenin axis: An in vitro and in vivo study

The progression of osteoarthritis (OA) includes the initial inflammation, subsequent degradation of the extracellular matrix (ECM), and chondrocyte apoptosis. Down syndrome candidate region 1 (DSCR1) is a stress-responsive gene and expresses in varied types of cells, including chondrocytes. Bioinformatics analysis of GSE103416 and GSE104739 datasets showed higher DSCR1 expression in the inflamed cartilage tissues and chondrocytes of OA. DSCR1 had two major isoforms, isoform 1 (DSCR1-1) and isoform 4 (DSCR1-4). We found that DSCR1-1 had a faster (in vitro) and higher expression (in vivo) response to OA compared to DSCR1-4. IL-1β-induced apoptosis, inflammation, and ECM degradation in chondrocytes were attenuated by DSCR1-1 overexpression. DSCR1-1 triggered the phosphorylation of cAMP response element-binding 1 (CREB1) at 133 serine sites by decreasing calcineurin activity. Moreover, activated CREB1 moved into the cell nucleus and combined in the promoter regions of aldehyde dehydrogenase 2 (ALDH2), thus enhancing its gene transcription. ALDH2 could recover Wnt/β-catenin signaling transduction by enhancing phosphorylation of β-catenin at 33/37 serine sites and inhibiting the migration of β-catenin protein from the cellular matrix to the nucleus. In vivo, adenoviruses (1 × 108 PFU) overexpressing DSCR1-1 were injected into the articular cavity of C57BL/6 mice with medial meniscus surgery-induced OA, and it showed that DSCR1-1 overexpression ameliorated cartilage injury. Collectively, our study demonstrates that DSCR1-1 may be a potential therapeutic target of OA.

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.

PTH treatment before cyclic joint loading improves cartilage health and attenuates load-induced osteoarthritis development in mice

Osteoarthritis (OA) treatment is limited by the lack of effective nonsurgical interventions to slow disease progression. Here, we examined the contributions of the subchondral bone properties to OA development. We used parathyroid hormone (PTH) to modulate bone mass before OA initiation and alendronate (ALN) to inhibit bone remodeling during OA progression. We examined the spatiotemporal progression of joint damage by combining histopathological and transcriptomic analyses across joint tissues. The additive effect of PTH pretreatment before OA initiation and ALN treatment during OA progression most effectively attenuated load-induced OA pathology. Individually, PTH directly improved cartilage health and slowed the development of cartilage damage, whereas ALN primarily attenuated subchondral bone changes associated with OA progression. Joint damage reflected early transcriptomic changes. With both treatments, the structural changes were associated with early modulation of immunoregulation and immunoresponse pathways that may contribute to disease mechanisms. Overall, our results demonstrate the potential of subchondral bone-modifying therapies to slow the progression of OA.

Data-driven identification of predictive risk biomarkers for subgroups of osteoarthritis using interpretable machine learning

Osteoarthritis (OA) is increasing in prevalence and has a severe impact on patients' lives. However, our understanding of biomarkers driving OA risk remains limited. We developed a model predicting the five-year risk of OA diagnosis, integrating retrospective clinical, lifestyle and biomarker data from the UK Biobank (19,120 patients with OA, ROC-AUC: 0.72, 95%CI (0.71-0.73)). Higher age, BMI and prescription of non-steroidal anti-inflammatory drugs contributed most to increased OA risk prediction ahead of diagnosis. We identified 14 subgroups of OA risk profiles. These subgroups were validated in an independent set of patients evaluating the 11-year OA risk, with 88% of patients being uniquely assigned to one of the 14 subgroups. Individual OA risk profiles were characterised by personalised biomarkers. Omics integration demonstrated the predictive importance of key OA genes and pathways (e.g., GDF5 and TGF-β signalling) and OA-specific biomarkers (e.g., CRTAC1 and COL9A1). In summary, this work identifies opportunities for personalised OA prevention and insights into its underlying pathogenesis.

Three decades of advancements in osteoarthritis research: insights from transcriptomic, proteomic, and metabolomic studies

OBJECTIVE

Osteoarthritis (OA) is a complex disease involving contributions from both local joint tissues and systemic sources. Patient characteristics, encompassing sociodemographic and clinical variables, are intricately linked with OA rendering its understanding challenging. Technological advancements have allowed for a comprehensive analysis of transcripts, proteomes and metabolomes in OA tissues/fluids through omic analyses. The objective of this review is to highlight the advancements achieved by omic studies in enhancing our understanding of OA pathogenesis over the last three decades.

DESIGN

We conducted an extensive literature search focusing on transcriptomics, proteomics and metabolomics within the context of OA. Specifically, we explore how these technologies have identified individual transcripts, proteins, and metabolites, as well as distinctive endotype signatures from various body tissues or fluids of OA patients, including insights at the single-cell level, to advance our understanding of this highly complex disease.

RESULTS

Omic studies reveal the description of numerous individual molecules and molecular patterns within OA-associated tissues and fluids. This includes the identification of specific cell (sub)types and associated pathways that contribute to disease mechanisms. However, there remains a necessity to further advance these technologies to delineate the spatial organization of cellular subtypes and molecular patterns within OA-afflicted tissues.

CONCLUSIONS

Leveraging a multi-omics approach that integrates datasets from diverse molecular detection technologies, combined with patients' clinical and sociodemographic features, and molecular and regulatory networks, holds promise for identifying unique patient endophenotypes. This holistic approach can illuminate the heterogeneity among OA patients and, in turn, facilitate the development of tailored therapeutic interventions.

Identification of therapeutic targets in osteoarthritis by combining heterogeneous transcriptional datasets, drug-induced expression profiles, and known drug-target interactions

BACKGROUND

Osteoarthritis (OA) is a multifactorial, hypertrophic, and degenerative condition involving the whole joint and affecting a high percentage of middle-aged people. It is due to a combination of factors, although the pivotal mechanisms underlying the disease are still obscure. Moreover, current treatments are still poorly effective, and patients experience a painful and degenerative disease course.

METHODS

We used an integrative approach that led us to extract a consensus signature from a meta-analysis of three different OA cohorts. We performed a network-based drug prioritization to detect the most relevant drugs targeting these genes and validated in vitro the most promising candidates. We also proposed a risk score based on a minimal set of genes to predict the OA clinical stage from RNA-Seq data.

RESULTS

We derived a consensus signature of 44 genes that we validated on an independent dataset. Using network analysis, we identified Resveratrol, Tenoxicam, Benzbromarone, Pirinixic Acid, and Mesalazine as putative drugs of interest for therapeutics in OA for anti-inflammatory properties. We also derived a list of seven gene-targets validated with functional RT-qPCR assays, confirming the in silico predictions. Finally, we identified a predictive subset of genes composed of DNER, TNFSF11, THBS3, LOXL3, TSPAN2, DYSF, ASPN and HTRA1 to compute the patient's risk score. We validated this risk score on an independent dataset with a high AUC (0.875) and compared it with the same approach computed using the entire consensus signature (AUC 0.922).

CONCLUSIONS

The consensus signature highlights crucial mechanisms for disease progression. Moreover, these genes were associated with several candidate drugs that could represent potential innovative therapeutics. Furthermore, the patient's risk scores can be used in clinical settings.

Elevated secretion of pro-collagen I-alpha and vascular endothelial growth factor as biomarkers of acetabular labrum degeneration and calcification in hip osteoarthritis: An explant study

Background

Hip osteoarthritis (OA) involves structural degeneration of different joint compartments, including femoral head cartilage, periarticular ligaments and the acetabular labrum. However, the molecular mechanisms underlying labrum degeneration in hip OA remain poorly understood.

Aim

To assess secretion of putative biomarkers for OA from explanted human labrum tissues under basal and inflammatory conditions and to determine whether these could differentiate between OA and calcification status compared to fracture controls.

Methods

Intact labrum specimens were collected from patients undergoing joint arthroplasty for primary hip OA (n ​= ​15, mean age 70) or non-OA femoral neck fracture (n ​= ​5, mean age 64). Tissues were dissected in equal-sized samples and explanted for one week. To mimic activation of inflammatory signaling by endogenous damage-associated molecular patterns (DAMP) tissue were stimulated with a toll-like receptor 4 (TLR4) agonist (1 ​μg/mL LPS). The involvement of transforming growth factor-beta (TGF-beta) signaling was evaluated by treatment with a TGF-beta type 1 receptor inhibitor (10 ​μM SB-505124). Secretion of aggrecan (ACAN), pro-collagen-I alpha (Pro-Col-Iα), cartilage oligomeric matrix protein (COMP), interleukin-6 (IL-6) and vascular endothelial growth factor (VEGF) was assessed by enzyme-linked immunosorbent assay (ELISA). Labrum calcification was evaluated by 3D whole mount fluorescent microscopy of ethyl cinnamate-based optically cleared tissues stained with Alcian blue/Alizarin red.

Results

Whole mount microscopy revealed non-OA fracture controls were non-calcified, whereas six OA labra (40%) were partially calcified or ossified. Basal secretion of Pro-Col-Iα and VEGF was increased four-fold in OA versus non-OA labra. Pro-Col-Iα levels were correlated with those of VEGF (r ​= ​0.65) and COMP (r ​= ​0.54). Stimulation of DAMP signaling through TLR4 affected secretion of IL-6, VEGF, COMP and Pro-Col-Iα, with distinct responses between non-OA and OA tissues. Inhibition of TGF-beta signaling specifically reduced elevated secretion of Pro-Col- Iα and VEGF in calcified OA labrum.

Conclusions

Secretion of the putative OA biomarkers Pro-Col-Iα and VEGF is elevated in degenerated human acetabular labrum and may serve as indicators of OA and calcification status. Secretion of both factors was partially regulated by TGF-beta signaling in calcified OA labrum tissues.The Translational potential of this article:Our findings suggest that a biomarker panel consisting of Pro-Col-Iα/VEGF/COMP may be valuable for assessing subradiographic labrum degeneration and calcification in hip OA. Targeting TGF-beta signaling may offer a means to reduce vascular invasion and fibrosis in acetabular labrum tissue.

Human Umbilical Cord Mesenchymal Stem Cells Alleviate Rat Knee Osteoarthritis via Activating Wnt/ β-catenin Signaling Pathway

BACKGROUND

Osteoarthritis (OA) is a chronic disease characterized by joint cartilage degeneration, destruction, and osteogenic hyperplasia. Human umbilical cord mesenchymal stem cells (hUCMSCs) have attracted increasing research interest due to their high clonogenic, proliferative, and migratory potential, as well as their improved secretion of relevant chondrogenic factors. This study evaluated the therapeutic potential and underlying mechanism of hUC-MSCs in alleviating pathological symptoms of OA.

METHODS

For the in vivo study, OA rats were established by the Hulth method to observe the therapeutic effect of intra-articular injection of hUC-MSCs. X-ray tests, gross observations, and histological and immunohistochemical assessments were conducted in rats. Levels of interleukin-1 beta (IL-1β), IL-6, matrix metalloproteinase-13 (MMP-13), and tissue inhibitor matrix metalloproteinase-1 in rats' synovial fluid were measured using enzyme-linked immunosorbent assay kits. For the in vitro study, hUC-MSCs and chondrocytes were cultured to explore the effect and underlying mechanisms of hUC-MSCs on OA. Apoptosis, proliferation, and glycosaminoglycan (GAG) were measured in the chondrocytes. The relative expression of aggrecan, COL-2, and SOX-9 mRNA was quantified by real-time polymerase chain reaction. Expressions of Wnt/β-catenin signaling molecules were measured by Western blot.

RESULTS

We found that intra-articular injection of hUC-MSCs reduced the combined score, increased the expression of collagen II, and decreased the expression of MMP-13, IL-1β, and IL-6 in rat knee joints. Additionally, hUC-MSCs increased the content of GAGs, inhibited chondrocyte apoptosis, and promoted chondrocyte proliferation. The expression of aggrecan, COL-2, and SOX-9 mRNA in chondrocytes was promoted by hUC-MSCs via activation of the Wnt/β-catenin signaling pathway.

CONCLUSION

Overall, this study demonstrated that hUC-MSCs induce the secretion of some cytokines via the paracrine function to activate the Wnt/β-catenin signaling pathway to reduce the pathological condition of OA and maintain the proper expression of cytokines and extracellular matrix proteins.

Single-cell RNA sequencing in osteoarthritis

Osteoarthritis is a progressive and heterogeneous joint disease with complex pathogenesis. The various phenotypes associated with each patient suggest that better subgrouping of tissues associated with genotypes in different phases of osteoarthritis may provide new insights into the onset and progression of the disease. Recently, single-cell RNA sequencing was used to describe osteoarthritis pathogenesis on a high-resolution view surpassing traditional technologies. Herein, this review summarizes the microstructural changes in articular cartilage, meniscus, synovium and subchondral bone that are mainly due to crosstalk amongst chondrocytes, osteoblasts, fibroblasts and endothelial cells during osteoarthritis progression. Next, we focus on the promising targets discovered by single-cell RNA sequencing and its potential applications in target drugs and tissue engineering. Additionally, the limited amount of research on the evaluation of bone-related biomaterials is reviewed. Based on the pre-clinical findings, we elaborate on the potential clinical values of single-cell RNA sequencing for the therapeutic strategies of osteoarthritis. Finally, a perspective on the future development of patient-centred medicine for osteoarthritis therapy combining other single-cell multi-omics technologies is discussed. This review will provide new insights into osteoarthritis pathogenesis on a cellular level and the field of applications of single-cell RNA sequencing in personalized therapeutics for osteoarthritis in the future.

Tick tock, the cartilage clock

Osteoarthritis (OA) is the most common age-related joint disease, affecting articular cartilage and other joint structures, causing severe pain and disability. Due to a limited understanding of the underlying disease pathogenesis, there are currently no disease-modifying drugs for OA. Circadian rhythms are generated by cell-intrinsic timekeeping mechanisms which are known to dampen during ageing, increasing disease risks. In this review, we focus on one emerging area of chondrocyte biology, the circadian clocks. We first provide a historical perspective of circadian clock discoveries and the molecular underpinnings. We will then focus on the expression and functions of circadian clocks in articular cartilage, including their rhythmic target genes and pathways, links to ageing, tissue degeneration, and OA, as well as tissue niche-specific entrainment pathways. Further research into cartilage clocks and ageing may have broader implications in the understanding of OA pathogenesis, the standardization of biomarker detection, and the development of novel therapeutic routes for the prevention and management of OA and other musculoskeletal diseases.

PSD95 as a New Potential Therapeutic Target of Osteoarthritis: A Study of the Identification of Hub Genes through Self-Contrast Model

Osteoarthritis (OA) is a worldwide joint disease. However, the precise mechanism causing OA remains unclear. Our primary aim was to identify vital biomarkers associated with the mechano-inflammatory aspect of OA, providing potential diagnostic and therapeutic targets for OA. Thirty OA patients who underwent total knee arthroplasty were recruited, and cartilage samples were obtained from both the lateral tibial plateau (LTP) and medial tibial plateau (MTP). GO and KEGG enrichment analyses were performed, and the protein-protein interaction (PPI) assessment was conducted for hub genes. The effect of PSD95 inhibition on cartilage degeneration was also conducted and analyzed. A total of 1247 upregulated and 244 downregulated DEGs were identified. Significant differences were observed between MTP and LTP in mechanical stress-related genes and activated sensory neurons based on a self-contrast model of human knee OA. Cluster analysis identified DLG4 as the hub gene. Cyclic loading stress increased PSD95 (encoded by DLG4) expression in LTP cartilage, and PSD95 inhibitors could alleviate OA progression. This study suggests that inhibiting PSD95 could be a potential therapeutic strategy for preventing articular cartilage degradation.

The circadian clock and extracellular matrix homeostasis in aging and age-related diseases

The extracellular matrix (ECM) is the noncellular scaffolding component present within all tissues and organs. It provides crucial biochemical and biomechanical cues to instruct cellular behavior and has been shown to be under circadian clock regulation, a highly conserved cell-intrinsic timekeeping mechanism that has evolved with the 24-hour rhythmic environment. Aging is a major risk factor for many diseases, including cancer, fibrosis, and neurodegenerative disorders. Both aging and our modern 24/7 society disrupt circadian rhythms, which could contribute to altered ECM homeostasis. Understanding the daily dynamics of ECM and how this mechanism changes with age will have a profound impact on tissue health, disease prevention, and improving treatments. Maintaining rhythmic oscillations has been proposed as a hallmark of health. On the other hand, many hallmarks of aging turn out to be key regulators of circadian timekeeping mechanisms. In this review, we summarize new work linking the ECM with circadian clocks and tissue aging. We discuss how the changes in the biomechanical and biochemical properties of ECM during aging may contribute to circadian clock dysregulation. We also consider how the dampening of clocks with age could compromise the daily dynamic regulation of ECM homeostasis in matrix-rich tissues. This review aims to encourage new concepts and testable hypotheses about the two-way interactions between circadian clocks and ECM in the context of aging.

RNA-binding proteins that are highly expressed and enriched in healthy cartilage but suppressed in osteoarthritis

Objectives: RNA-binding proteins (RBPs) have diverse and essential biological functions, but their role in cartilage health and disease is largely unknown. The objectives of this study were (i) map the global landscape of RBPs expressed and enriched in healthy cartilage and dysregulated in osteoarthritis (OA); (ii) prioritize RBPs for their potential role in cartilage and in OA pathogenesis and as therapeutic targets. Methods: Our published bulk RNA-sequencing (RNA-seq) data of healthy and OA human cartilage, and a census of 1,542 RBPs were utilized to identify RBPs that are expressed in healthy cartilage and differentially expressed (DE) in OA. Next, our comparison of healthy cartilage RNA-seq data to 37 transcriptomes in the Genotype-Tissue Expression (GTEx) database was used to determine RBPs that are enriched in cartilage. Finally, expression of RBPs was analyzed in our single cell RNA-sequencing (scRNA-seq) data from healthy and OA human cartilage. Results: Expression of RBPs was higher than nonRBPs in healthy cartilage. In OA cartilage, 188 RBPs were differentially expressed, with a greater proportion downregulated. Ribosome biogenesis was enriched in the upregulated RBPs, while splicing and transport were enriched in the downregulated. To further prioritize RBPs, we selected the top 10% expressed RBPs in healthy cartilage and those that were cartilage-enriched according to GTEx. Intersecting these criteria, we identified Tetrachlorodibenzodioxin (TCDD) Inducible Poly (ADP-Ribose) Polymerase (TIPARP) as a candidate RBP. TIPARP was downregulated in OA. scRNA-seq data revealed TIPARP was most significantly downregulated in the "pathogenic cluster". Conclusion: Our global analyses reveal expression patterns of RBPs in healthy and OA cartilage. We also identified TIPARP and other RBPs as novel mediators in OA pathogenesis and as potential therapeutic targets.

Combining biomechanical stimulation and chronobiology: a novel approach for augmented chondrogenesis?

The unique structure and composition of articular cartilage is critical for its physiological function. However, this architecture may get disrupted by degeneration or trauma. Due to the low intrinsic regeneration properties of the tissue, the healing response is generally poor. Low-grade inflammation in patients with osteoarthritis advances cartilage degradation, resulting in pain, immobility, and reduced quality of life. Generating neocartilage using advanced tissue engineering approaches may address these limitations. The biocompatible microenvironment that is suitable for cartilage regeneration may not only rely on cells and scaffolds, but also on the spatial and temporal features of biomechanics. Cell-autonomous biological clocks that generate circadian rhythms in chondrocytes are generally accepted to be indispensable for normal cartilage homeostasis. While the molecular details of the circadian clockwork are increasingly well understood at the cellular level, the mechanisms that enable clock entrainment by biomechanical signals, which are highly relevant in cartilage, are still largely unknown. This narrative review outlines the role of the biomechanical microenvironment to advance cartilage tissue engineering via entraining the molecular circadian clockwork, and highlights how application of this concept may enhance the development and successful translation of biomechanically relevant tissue engineering interventions.

MSdb: An integrated expression atlas of human musculoskeletal system

The global prevalence and burden of musculoskeletal (MSK) disorders are immense. Advancements in next-generation sequencing (NGS) have generated vast amounts of data, accelerating the research of pathological mechanisms and the development of therapeutic approaches for MSK disorders. However, scattered datasets across various repositories complicate uniform analysis and comparison. Here, we introduce MSdb, a database for visualization and integrated analysis of next-generation sequencing data from human musculoskeletal system, along with manually curated patient phenotype data. MSdb provides various types of analysis, including sample-level browsing of metadata information, gene/miRNA expression, and single-cell RNA-seq dataset. In addition, MSdb also allows integrated analysis for cross-samples and cross-omics analysis, including customized differentially expressed gene/microRNA analysis, microRNA-gene network, scRNA-seq cross-sample/disease integration, and gene regulatory network analysis. Overall, systematic categorizing, standardized processing, and freely accessible knowledge features MSdb a valuable resource for MSK research community.

Osteoarthritic chondrocytes undergo a glycolysis-related metabolic switch upon exposure to IL-1b or TNF

BACKGROUND

Osteoarthritis is an age-related disease that currently faces a lack of symptomatic treatment. Inflammation, which is mainly sustained by pro-inflammatory cytokines such as IL-1b, TNF, and IL-6, plays an important role in osteoarthritis progression. In this context, pro-inflammatory cytokines are widely used to mimic the inflammatory component of osteoarthritis in vitro. However, the therapeutic failures of clinical trials evaluating anti-cytokines drugs highlight the lack of overall understanding of the effects of these cytokines on chondrocytes.

METHODS

Here, we generated a comprehensive transcriptomic and proteomic dataset of osteoarthritic chondrocytes treated with these cytokines to describe their pro-inflammatory signature and compare it to the transcriptome of non-osteoarthritic chondrocytes. Then, the dysregulations highlighted at the molecular level were functionally confirmed by real-time cellular metabolic assays.

RESULTS

We identified dysregulation of metabolic-related genes in osteoarthritic chondrocytes but not in non-osteoarthritic chondrocytes. A metabolic shift, toward increased glycolysis at the expense of mitochondrial respiration, was specifically confirmed in osteoarthritic chondrocytes treated with IL-1b or TNF.

CONCLUSION

These data show a strong and specific association between inflammation and metabolism in osteoarthritic chondrocytes, which was not found in non-osteoarthritic chondrocytes. This indicates that the link between inflammation and metabolic dysregulation may be exacerbated during chondrocyte damage in osteoarthritis. Video Abstract.

Ozone Therapy versus Hyaluronic Acid Injections for Pain Relief in Patients with Knee Osteoarthritis: Preliminary Findings on Molecular and Clinical Outcomes from a Randomized Controlled Trial

Ozone therapy (OT) is used for the treatment of multiple musculoskeletal disorders. In recent years, there has been a growing interest in its use for the treatment of osteoarthritis (OA). The aim of this double-blind randomized controlled trial was to evaluate the efficacy of OT compared with hyaluronic acid (HA) injections for pain relief in patients with knee OA. Patients with knee OA for at least three months were included and randomly assigned to receive three intra-articular injections of ozone or HA (once a week). Patients were assessed at baseline and at 1, 3, and 6 months after the injections for pain, stiffness, and function using the WOMAC LK 3.1, the NRS, and the KOOS questionnaire. Out of 55 patients assessed for eligibility, 52 participants were admitted to the study and randomly assigned into the 2 groups of treatment. During the study, eight patients dropped out. Thus, a total of 44 patients, reached the endpoint of the study at 6 months. Both Group A and B consisted of 22 patients. At 1-month follow-up after injections, both treatment groups improved statistically significantly from baseline in all outcomes measured. At 3 months, improvements remained similarly consistent for Group A and Group B. At 6-month follow-up, the outcomes were comparable between the 2 groups, showing only a worsening trend in pain. No significant differences were found between the two groups in pain scores. Both therapies have proven to be safe, with the few recorded adverse events being mild and self-limiting. OT has demonstrated similar results to HA injections, proving to be a safe approach with significant effects on pain control in patients affected by knee OA. Due to its anti-inflammatory and analgesic effects, ozone might be considered as a potential treatment for OA.

Transcriptomic changes in porcine articular cartilage one year following disruption of the anterior cruciate ligament

To determine the transcriptomic changes seen in early- to mid-stage posttraumatic osteoarthritis (PTOA) development, 72 Yucatan minipigs underwent transection of the anterior cruciate ligament. Subjects were randomized to no further intervention, ligament reconstruction, or ligament repair, followed by articular cartilage harvesting and RNA-sequencing at three different postoperative timepoints (1, 4, and 52 weeks). Six additional subjects received no ligament transection and provided cartilage tissue to serve as controls. Differential gene expression analysis between post-transection cartilage and healthy cartilage revealed an initial increase in transcriptomic differences at 1 and 4 weeks followed by a stark reduction in transcriptomic differences at 52 weeks. This analysis also showed how different treatments genetically modulate the course of PTOA following ligament disruption. Specific genes (e.g., MMP1, POSTN, IGF1, PTGFR, HK1) were identified as being upregulated in the cartilage of injured subjects across all timepoints regardless of treatment. At the 52-week timepoint, 4 genes (e.g., A4GALT, EFS, NPTXR, ABCA3) that-as far as we know-have yet to be associated with PTOA were identified as being concordantly differentially expressed across all treatment groups when compared to controls. Functional pathway analysis of injured subject cartilage compared to control cartilage revealed overarching patterns of cellular proliferation at 1 week, angiogenesis, ECM interaction, focal adhesion, and cellular migration at 4 weeks, and calcium signaling, immune system activation, GABA signaling, and HIF-1 signaling at 52 weeks.

Re-thinking osteoarthritis pathogenesis: what can we learn (and what do we need to unlearn) from mouse models about the mechanisms involved in disease development

Efforts to develop effective disease-modifying drugs to treat osteoarthritis have so far proved unsuccessful with a number of promising drug candidates from pre-clinical studies failing to show efficacy in clinical trials. It is therefore timely to re-evaluate our current understanding of osteoarthritis pathogenesis and the similarities and differences in disease development between commonly used pre-clinical mouse models and human patients. There is substantial heterogeneity between patients presenting with osteoarthritis and mounting evidence that the pathways involved in osteoarthritis (e.g. Wnt signalling) differ between patient sub-groups. There is also emerging evidence that the pathways involved in osteoarthritis differ between the STR/ort mouse model (the most extensively studied mouse model of spontaneously occurring osteoarthritis) and injury-induced osteoarthritis mouse models. For instance, while canonical Wnt signalling is upregulated in the synovium and cartilage at an early stage of disease in injury-induced osteoarthritis mouse models, this does not appear to be the case in the STR/ort mouse. Such findings may prove insightful for understanding the heterogeneity in mechanisms involved in osteoarthritis pathogenesis in human disease. However, it is important to recognise that there are differences between mice and humans in osteoarthritis pathogenesis. A much more extensive array of pathological changes are evident in osteoarthritic joints in individual mice with osteoarthritis compared to individual patients. There are also specified differences in the pathways involved in disease development. For instance, although increased TGF-β signalling is implicated in osteoarthritis development in both mouse models of osteoarthritis and human disease, in mice, this is mainly mediated through TGF-β3 whereas in humans, it is through TGF-β1. Studies in other tissues have shown TGF-β1 is more potent than TGF-β3 in inducing the switch to SMAD1/5 signalling that occurs in osteoarthritic cartilage and that TGF-β1 and TGF-β3 have opposing effects on fibrosis. It is therefore possible that the relative contribution of TGF-β signalling to joint pathology in osteoarthritis differs between murine models and humans. Understanding the similarities and differences in osteoarthritis pathogenesis between mouse models and humans is critical for understanding the translational potential of findings from pre-clinical studies.

Analysis of RNA Polyadenylation in Healthy and Osteoarthritic Human Articular Cartilage

Polyadenylation (polyA) defines the 3' boundary of a transcript's genetic information. Its position can vary and alternative polyadenylation (APA) transcripts can exist for a gene. This causes variance in 3' regulatory domains and can affect coding sequence if intronic events occur. The distribution of polyA sites on articular chondrocyte transcripts has not been studied so we aimed to define their transcriptome-wide location in age-matched healthy and osteoarthritic knee articular cartilage. Total RNA was isolated from frozen tissue samples and analysed using the QuantSeq-Reverse 3' RNA sequencing approach, where each read runs 3' to 5' from within the polyA tail into the transcript and contains a distinct polyA site. Differential expression of transcripts was significant altered between healthy and osteoarthritic samples with enrichment for functionalities that were strongly associated with joint pathology. Subsequent examination of polyA site data allowed us to define the extent of site usage across all the samples. When comparing healthy and osteoarthritic samples, we found that differential use of polyadenylation sites was modest. However, in the genes affected, there was potential for the APA to have functional relevance. We have characterised the polyadenylation landscape of human knee articular chondrocytes and conclude that osteoarthritis does not elicit a widespread change in their polyadenylation site usage. This finding differentiates knee osteoarthritis from pathologies such as cancer where APA is more commonly observed.

Development of a Clinical Prediction Rule for Treatment Success with Transcranial Direct Current Stimulation for Knee Osteoarthritis Pain: A Secondary Analysis of a Double-Blind Randomized Controlled Trial

The study’s objective was to develop a clinical prediction rule that predicts a clinically significant analgesic effect on chronic knee osteoarthritis pain after transcranial direct current stimulation treatment. This is a secondary analysis from a double-blind randomized controlled trial. Data from 51 individuals with chronic knee osteoarthritis pain and an impaired descending pain inhibitory system were used. The intervention comprised a 15-session protocol of anodal primary motor cortex transcranial direct current stimulation. Treatment success was defined by the Western Ontario and McMaster Universities’ Osteoarthritis Index pain subscale. Accuracy statistics were calculated for each potential predictor and for the final model. The final logistic regression model was statistically significant (p < 0.01) and comprised five physical and psychosocial predictor variables that together yielded a positive likelihood ratio of 14.40 (95% CI: 3.66−56.69) and an 85% (95%CI: 60−96%) post-test probability of success. This is the first clinical prediction rule proposed for transcranial direct current stimulation in patients with chronic pain. The model underscores the importance of both physical and psychosocial factors as predictors of the analgesic response to transcranial direct current stimulation treatment. Validation of the proposed clinical prediction rule should be performed in other datasets.

Characterization of Genetic Risk of End-Stage Knee Osteoarthritis Treated with Total Knee Arthroplasty: A Genome-Wide Association Study

BACKGROUND

End-stage knee osteoarthritis (OA) is a highly debilitating disease for which total knee arthroplasty (TKA) serves as an effective treatment option. Although a genetic component to OA in general has been described, evaluation of the genetic contribution to end-stage OA of the knee is limited. To this end, we present a genome-wide association study involving patients undergoing TKA for primary knee OA to characterize the genetic features of severe disease on a population level.

METHODS

Individuals with the diagnosis of knee OA who underwent primary TKA were identified in the U.K. Biobank using administrative codes. The U.K. Biobank is a data repository containing prospectively collected clinical and genomic data for >500,000 patients. A genome-wide association analysis was performed using the REGENIE software package. Logistic regression was also used to compare the total genetic risk between subgroups stratified by age and body mass index (BMI).

RESULTS

A total of 16,032 patients with end-stage knee OA who underwent primary TKA were identified. Seven genetic loci were found to be significantly associated with end-stage knee OA. The odds ratio (OR) for developing end-stage knee OA attributable to genetics was 1.12 (95% confidence interval [CI], 1.10 to 1.14), which was lower than the OR associated with BMI (OR = 1.81; 95% CI, 1.78 to 1.83) and age (OR = 2.38; 95% CI, 2.32 to 2.45). The magnitude of the OR for developing end-stage knee OA attributable to genetics was greater in patients <60 years old than in patients ≥60 years old (p = 0.002).

CONCLUSIONS

This population-level genome-wide association study of end-stage knee OA treated with primary TKA was notable for identifying multiple significant genetic variants. These loci involve genes responsible for cartilage development, cartilage homeostasis, cell signaling, and metabolism. Age and BMI appear to have a greater impact on the risk of developing end-stage disease compared with genetic factors. The genetic contribution to the development of severe disease is greater in younger patients.

LEVEL OF EVIDENCE

Prognostic Level III . See Instructions for Authors for a complete description of levels of evidence.

Oscillations of the circadian clock protein, BMAL-1, align to daily cycles of mechanical stimuli: a novel means to integrate biological time within predictive in vitro model systems

PURPOSE

In vivo, the circadian clock drives 24-h rhythms in human physiology. Isolated cells in vitro retain a functional clockwork but lack necessary timing cues resulting in the rapid loss of tissue-level circadian rhythms. This study tests the hypothesis that repeated daily mechanical stimulation acts as a timing cue for the circadian clockwork. The delineation and integration of circadian timing cues into predictive in vitro model systems, including organ-on-a-chip (OOAC) devices, represent a novel concept that introduces a key component of in vivo physiology into predictive in vitro model systems.

METHODS

Quiescent bovine chondrocytes were entrained for 3 days by daily 12-h bouts of cyclic biaxial tensile strain (10%, 0.33 Hz, Flexcell) before sampling during free-running conditions. The core clock protein, BMAL-1, was quantified from normalised Western Blot signal intensity and the temporal oscillations characterised by Cosinor linear fit with 24-h period.

RESULTS

Following entrainment, the cell-autonomous oscillations of the molecular clock protein, BMAL-1, exhibited circadian (24 h) periodicity (p < 0.001) which aligned to the diurnal mechanical stimuli. A 6-h phase shift in the mechanical entrainment protocol resulted in an equivalent shift of the circadian clockwork. Thus, repeated daily mechanical stimuli synchronised circadian rhythmicity of chondrocytes in vitro.

CONCLUSION

This work demonstrates that daily mechanical stimulation can act as a timing cue that is sufficient to entrain the peripheral circadian clock in vitro. This discovery may be exploited to induce and sustain circadian physiology within into predictive in vitro model systems, including OOAC systems. Integration of the circadian clock within these systems will enhance their potential to accurately recapitulate human diurnal physiology and hence augment their predictive value as drug testing platforms and as realistic models of human (patho)physiology.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s44164-022-00032-x.

Complement C3-α and C3-β Levels in Synovial Fluid But Not in Blood Correlate With the Severity of Osteoarthritis Research Society International Histopathological Grades in Primary Knee Osteoarthritis

BACKGROUND

Primary osteoarthritis (OA) is the most common cause of knee arthritis worldwide. The knee synovial fluid complement C3-β chain levels have been shown to correlate with clinical symptoms of knee OA. It is not known whether the complement C3 in the synovial fluid is derived from the circulation or is produced locally in the knee.

METHODS

Fifty primary OA patients undergoing a total knee arthroplasty procedure were evaluated for biochemical analyses of C3-α and C3-β chains in the synovial fluid and blood plasma. These levels were corelated with the severity of corresponding knee OA based on the Osteoarthritis Research Society International (OARSI) grade.

RESULTS

Both synovial C3-α and C3-β levels correlated significantly with the severity of OA. Neither plasma C3-α levels nor C3-β levels significantly correlated with OARSI grading. Neither synovial C3-α levels nor C3-β correlated significantly with plasma C3-α or C3-β levels, respectively. Synovial C3-α chain and C3-β chain levels were significantly higher in the grade >6 group. In plasma, neither C3-α chain levels nor C3-β chain levels were significantly different between the groups. Neither synovial C3-α nor C3-β levels significantly correlated with plasma erythrocyte sedimentation rate or C-reactive protein levels.

CONCLUSION

In knee primary OA, C3 seems to be produced and released locally into the synovial fluid instead of being derived from blood in the circulation. Synovial C3 levels, but not blood plasma C3, correlate with the histopathological severity of primary OA in the knee. Synovial C3 may be an important factor in the pathogenesis of primary OA clinical symptoms and a potential target for therapeutic intervention.

The circadian clock has roles in mesenchymal stem cell fate decision

The circadian clock refers to the intrinsic biological rhythms of physiological functions and behaviours. It synergises with the solar cycle and has profound effects on normal metabolism and organismal fitness. Recent studies have suggested that the circadian clock exerts great influence on the differentiation of stem cells. Here, we focus on the close relationship between the circadian clock and mesenchymal stem cell fate decisions in the skeletal system. The underlying mechanisms include hormone signals and the activation and repression of different transcription factors under circadian regulation. Additionally, the clock interacts with epigenetic modifiers and non-coding RNAs and is even involved in chromatin remodelling. Although the specificity and safety of circadian therapy need to be further studied, the circadian regulation of stem cells can be regarded as a promising candidate for health improvement and disease prevention.

Osteoarthritis Pathophysiology: Therapeutic Target Discovery may Require a Multifaceted Approach

Molecular understanding of osteoarthritis (OA) has greatly increased through careful analysis of tissue samples, preclinical models, and large-scale agnostic "-omic" studies. There is broad acceptance that systemic and biomechanical signals affect multiple tissues of the joint, each of which could potentially be targeted to improve patient outcomes. In this review six experts in different aspects of OA pathogenesis provide their independent view on what they believe to be good tractable approaches to OA target discovery. We conclude that molecular discovery has been high but future transformative studies require a multidisciplinary holistic approach to develop therapeutic strategies with high clinical efficacy.

[Osteoarthritis: what's new?]

Osteoarthritis (OA) is a very common disease. As a consequence of the ageing society, osteoarthritis prevalence will further increase. Age itself, trauma, unequal load distribution and overweight are risk factors. Cellular senescence and overweight have been in the focus of scientific interest for the last few years. Both risk factors are able to facilitate joint inflammation, independent of a mechanical approach. Senescent chondrocytes as well as adipocytes can produce increased amounts of inflammatory cytokines. Cornerstones of the therapy are patient education including information on the character/course of the disease and intentional weight loss. Although NSAIDs can be recommended as analgesics, their contraindications limit the widespread use. Alternatively, acetaminophen or low-potency opioids such as tramadol might be considered. Topical NSAIDs and intraarticular glucocorticoid injections can be helpful in pain reduction particularly in knee osteoarthritis. There is still no general recommendation for nutritional supplements including chondroitin or glycosaminoglycan, but they might be considered as an accompanying therapy. With the current non-approval of the nerve growth factor (NGF)-antibody tanezumab, a new therapeutical option for OA suffered a setback. Unfortunately, the results of the phase 2 study on the Wnt inhibitor lorecivivint are barely encouraging. However, the results of the according phase 3 study are eagerly awaited.

OA foundations - experimental models of osteoarthritis

Osteoarthritis (OA) is increasingly recognised as a disease of diverse phenotypes with variable clinical presentation, progression, and response to therapeutic intervention. This same diversity is readily apparent in the many animal models of OA. However, model selection, study design, and interpretation of resultant findings, are not routinely done in the context of the target human (or veterinary) patient OA sub-population or phenotype. This review discusses the selection and use of animal models of OA in discovery and therapeutic-development research. Beyond evaluation of the different animal models on offer, this review suggests focussing the approach to OA-animal model selection on study objective(s), alignment of available models with OA-patient sub-types, and the resources available to achieve valid and translatable results. How this approach impacts model selection is discussed and an experimental design checklist for selecting the optimal model(s) is proposed. This approach should act as a guide to new researchers and a reminder to those already in the field, as to issues that need to be considered before embarking on in vivo pre-clinical research. The ultimate purpose of using an OA animal model is to provide the best possible evidence if, how, when and where a molecule, pathway, cell or process is important in clinical disease. By definition this requires both model and study outcomes to align with and be predictive of outcomes in patients. Keeping this at the forefront of research using pre-clinical OA models, will go a long way to improving the quality of evidence and its translational value.

Osteoarthritis year in review: genetics, genomics, epigenetics

OBJECTIVE

In this review, we have highlighted the advances over the past year in genetics, genomics and epigenetics in the field of osteoarthritis (OA).

METHODS

A literature search of PubMed was performed using the criteria: "osteoarthritis" and one of the following terms "genetic(s), genomic(s), epigenetic(s), polymorphism, noncoding ribonucleic acid (RNA), microRNA, long noncoding RNA, lncRNA, circular RNA, RNA sequencing (RNA-seq), single cell sequencing, transcriptomics, or deoxyribonucleic acid (DNA) methylation between April 01, 2020 and April 30, 2021.

RESULTS

In total we identified 765 unique publications, which eventually reduced to 380 of relevance to the field as judged by two assessors. Many of these studies included multiple search terms. We summarised advances relating to genetics, functional genetics, genomics and epigenetics, focusing on our personal key papers during the year.

CONCLUSIONS

This year few studies have identified new genetic variants contributing to OA susceptibility, but a focus has been on refining risk loci or their functional validation. The use of new technologies together with investigating the cross-talk between multiple tissue types, greater sample sizes and/or better patient classification (OA subtypes) will continue to increase our knowledge of disease mechanisms and progress towards understanding and treating OA.

A molecular map of long non-coding RNA expression, isoform switching and alternative splicing in osteoarthritis

Osteoarthritis is a prevalent joint disease and a major cause of disability worldwide with no curative therapy. Development of disease-modifying therapies requires a better understanding of the molecular mechanisms underpinning disease. A hallmark of osteoarthritis is cartilage degradation. To define molecular events characterizing osteoarthritis at the whole transcriptome level, we performed deep RNA sequencing in paired samples of low- and high-osteoarthritis grade knee cartilage derived from 124 patients undergoing total joint replacement. We detected differential expression between low- and high-osteoarthritis grade articular cartilage for 365 genes and identified a 38-gene signature in osteoarthritis cartilage by replicating our findings in an independent dataset. We also found differential expression for 25 novel long non-coding RNA genes (lncRNAs) and identified potential lncRNA interactions with RNA-binding proteins in osteoarthritis. We assessed alterations in the relative usage of individual gene transcripts and identified differential transcript usage for 82 genes, including ABI3BP, coding for an extracellular matrix protein, AKT1S1, a negative regulator of the mTOR pathway and TPRM4, coding for a transient receptor potential channel. We further assessed genome-wide differential splicing, for the first time in osteoarthritis, and detected differential splicing for 209 genes, which were enriched for extracellular matrix, proteoglycans and integrin surface interactions terms. In the largest study of its kind in osteoarthritis, we find that isoform and splicing changes, in addition to extensive differences in both coding and non-coding sequence expression, are associated with disease and demonstrate a novel layer of genomic complexity to osteoarthritis pathogenesis.

The Homeostasis of Cartilage Matrix Remodeling and the Regulation of Volume-Sensitive Ion Channel

Degenerative joint diseases of the hips and knees are common and are accompanied by severe pain and movement disorders. At the microscopic level, the main characteristics of osteoarthritis are the continuous destruction and degeneration of cartilage, increased cartilage extracellular matrix catabolism, decreased anabolism, increased synovial fluid, and decreased osmotic pressure. Cell volume stability is mainly regulated by ion channels, many of which are expressed in chondrocytes. These ion channels are closely related to pain regulation, volume regulation, the inflammatory response, cell proliferation, apoptosis, and cell differentiation. In this review, we focus on the important role of volume control-related ion channels in cartilage matrix remodeling and summarize current views. In addition, the potential mechanism of the volume-sensitive anion channel LRRC8A in the early occurrence of osteoarthritis is discussed.

Assessment of Clinical, Tissue, and Cell-Level Metrics Identify Four Biologically Distinct Knee Osteoarthritis Patient Phenotypes

OBJECTIVE

Clinical heterogeneity of primary osteoarthritis (OA) is a major challenge in understanding pathogenesis and development of targeted therapeutic strategies. This study aims to (1) identify OA patient subgroups phenotypes and (2) determine predictors of OA severity and cartilage-derived stem/progenitor concentration using clinical-, tissue-, and cell- level metrics.

DESIGN

Cartilage, synovium (SYN) and infrapatellar fatpad (IPFP) were collected from 90 total knee arthroplasty patients. Clinical metrics (patient demographics, radiograph-based joint space width (JSW), Kellgren and Lawrence score (KL)), tissue metrics (cartilage histopathology grade, glycosaminoglycans (GAGs)) and cell-based metrics (cartilage-, SYN-, and IPFP-derived cell concentration ([Cell], cells/mg), connective tissue progenitor (CTP) prevalence (PCTP, CTPs/million cells plated), CTP concentration, [CTP], CTPs/mg)) were assessed using k-mean clustering and linear regression model.

RESULTS

Four patient subgroups were identified. Clusters 1 and 2 comprised of younger, high body mass index (BMI) patients with healthier cartilage, where Cluster 1 had high CTP in cartilage, SYN, and IPFP, and Cluster 2 had low [CTP] in cartilage, SYN, and IPFP. Clusters 3 and 4 comprised of older, low BMI patients with diseased cartilage where Cluster 3 had low [CTP] in SYN, IPFP but high [CTP] in cartilage, and Cluster 4 had high [CTP] in SYN, IPFP but low [CTP] in cartilage. Age (r = 0.23, P = 0.026), JSW (r = 0.28, P = 0.007), KL (r = 0.26, P = 0.012), GAG/mg cartilage tissue (r = -0.31, P = 0.007), and SYN-derived [Cell] (r = 0.25, P = 0.049) were weak but significant predictors of OA severity. Cartilage-derived [Cell] (r = 0.38, P < 0.001) and PCTP (r = 0.9, P < 0.001) were moderate/strong predictors of cartilage-derived [CTP].

CONCLUSION

Initial findings suggests the presence of OA patient subgroups that could define opportunities for more targeted patient-specific approaches to prevention and treatment.

Development of human cartilage circadian rhythm in a stem cell-chondrogenesis model

The circadian clock in murine articular cartilage is a critical temporal regulatory mechanism for tissue homeostasis and osteoarthritis. However, translation of these findings into humans has been hampered by the difficulty in obtaining circadian time series human cartilage tissues. As such, a suitable model is needed to understand the initiation and regulation of circadian rhythms in human cartilage. Methods: We used a chondrogenic differentiation protocol on human embryonic stem cells (hESCs) as a proxy for early human chondrocyte development. Chondrogenesis was validated using histology and expression of pluripotency and differentiation markers. The molecular circadian clock was tracked in real time by lentiviral transduction of human clock gene luciferase reporters. Differentiation-coupled gene expression was assessed by RNAseq and differential expression analysis. Results: hESCs lacked functional circadian rhythms in clock gene expression. During chondrogenic differentiation, there was an expected reduction of pluripotency markers (e.g., NANOG and OCT4) and a significant increase of chondrogenic genes (SOX9, COL2A1 and ACAN). Histology of the 3D cartilage pellets at day 21 showed a matrix architecture resembling human cartilage, with readily detectable core clock proteins (BMAL1, CLOCK and PER2). Importantly, the circadian clocks in differentiating hESCs were activated between day 11 (end of the 2D stage) and day 21 (10 days after 3D differentiation) in the chondrogenic differentiation protocol. RNA sequencing revealed striking differentiation coupled changes in the expression levels of most clock genes and a range of clock regulators. Conclusions: The circadian clock is gradually activated through a differentiation-coupled mechanism in a human chondrogenesis model. These findings provide a human 3D chondrogenic model to investigate the role of the circadian clock during normal homeostasis and in diseases such as osteoarthritis.

Low-Input RNA-Sequencing in Patients with Cartilage Lesions, Osteoarthritis, and Healthy Cartilage

OBJECTIVE

To analyze and compare cartilage samples from 3 groups of patients utilizing low-input RNA-sequencing.

DESIGN

Cartilage biopsies were collected from patients in 3 groups (n = 48): Cartilage lesion (CL) patients had at least ICRS grade 2, osteoarthritis (OA) samples were taken from patients undergoing knee replacement, and healthy cartilage (HC) was taken from ACL-reconstruction patients without CLs. RNA was isolated using an optimized protocol. RNA samples were assessed for quality and sequenced with a low-input SmartSeq2 protocol.

RESULTS

RNA isolation yielded 48 samples with sufficient quality for sequencing. After quality control, 13 samples in the OA group, 9 in the HC group, and 9 in the CL group were included in the analysis. There was a high degree of co-clustering between the HC and CL groups with only 6 genes significantly up- or downregulated. OA and the combined HC/CL group clustered significantly separate from each other, yielding 659 significantly upregulated and 1,369 downregulated genes. GO-term analysis revealed that genes matched to cartilage and connective tissue development terms.

CONCLUSION

The gene expression profiles from the 3 groups suggest that there are no major differences in gene expression between cartilage from knees with a cartilage injury and knees without an apparent cartilage injury. OA cartilage, as expected, showed markedly different gene expression from the other 2 groups. The gene expression profiles resulting from this low-input RNA-sequencing study offer opportunities to discover new pathways not previously recognized that may be explored in future studies.

The non-coding RNA interactome in joint health and disease

Non-coding RNAs have distinct regulatory roles in the pathogenesis of joint diseases including osteoarthritis (OA) and rheumatoid arthritis (RA). As the amount of high-throughput profiling studies and mechanistic investigations of microRNAs, long non-coding RNAs and circular RNAs in joint tissues and biofluids has increased, data have emerged that suggest complex interactions among non-coding RNAs that are often overlooked as critical regulators of gene expression. Identifying these non-coding RNAs and their interactions is useful for understanding both joint health and disease. Non-coding RNAs regulate signalling pathways and biological processes that are important for normal joint development but, when dysregulated, can contribute to disease. The specific expression profiles of non-coding RNAs in various disease states support their roles as promising candidate biomarkers, mediators of pathogenic mechanisms and potential therapeutic targets. This Review synthesizes literature published in the past 2 years on the role of non-coding RNAs in OA and RA with a focus on inflammation, cell death, cell proliferation and extracellular matrix dysregulation. Research to date makes it apparent that 'non-coding' does not mean 'non-essential' and that non-coding RNAs are important parts of a complex interactome that underlies OA and RA.

Comprehensive Analysis of Key Genes, Signaling Pathways and miRNAs in Human Knee Osteoarthritis: Based on Bioinformatics

Background: Osteoarthritis (OA) is one of the main causes of disability in the elderly population, accompanied by a series of underlying pathologic changes, such as cartilage degradation, synovitis, subchondral bone sclerosis, and meniscus injury. The present study aimed to identify key genes, signaling pathways, and miRNAs in knee OA associated with the entire joint components, and to explain the potential mechanisms using computational analysis.

Methods: The differentially expressed genes (DEGs) in cartilage, synovium, subchondral bone, and meniscus were identified using the Gene Expression Omnibus 2R (GEO2R) analysis based on dataset from GSE43923, GSE12021, GSE98918, and GSE51588, respectively and visualized in Volcano Plot. Venn diagram analyses were performed to identify the overlapping DEGs (overlapping DEGs) that expressed in at least two types of tissues mentioned above. Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, protein-protein interaction (PPI) analysis, and module analysis were conducted. Furthermore, qRT-PCR was performed to validate above results using our clinical specimens.

Results: As a result, a total of 236 overlapping DEGs were identified, of which 160 were upregulated and 76 were downregulated. Through enrichment analysis and constructing the PPI network and miRNA-mRNA network, knee OA-related key genes, such as HEY1, AHR, VEGFA, MYC, and CXCL12 were identified. Clinical validation by qRT-PCR experiments further supported above computational results. In addition, knee OA-related key miRNAs such as miR-101, miR-181a, miR-29, miR-9, and miR-221, and pathways such as Wnt signaling, HIF-1 signaling, PI3K-Akt signaling, and axon guidance pathways were also identified. Among above identified knee OA-related key genes, pathways and miRNAs, genes such as AHR, HEY1, MYC, GAP43, and PTN, pathways like axon guidance, and miRNAs such as miR-17, miR-21, miR-155, miR-185, and miR-1 are lack of research and worthy for future investigation.

Conclusion: The present informatic study for the first time provides insight to the potential therapeutic targets of knee OA by comprehensively analyzing the overlapping genes differentially expressed in multiple joint components and their relevant signaling pathways and interactive miRNAs.

The circadian clock: a central mediator of cartilage maintenance and osteoarthritis development?

The circadian clock is a specialized cell signalling pathway present in all cells. Loss of clock function leads to tissue degeneration and premature ageing in animal models demonstrating the fundamental importance of clocks for cell, tissue and organism health. There is now considerable evidence that the chondrocyte circadian clock is altered in OA. The purpose of this review is to summarize current knowledge regarding the nature of the change in the chondrocyte clock in OA and the implications of this change for disease development. Expression of the core clock component, BMAL1, has consistently been shown to be lower in OA chondrocytes. This may contribute to changes in chondrocyte differentiation and extracellular matrix turnover in disease. Circadian clocks are highly responsive to environmental factors. Mechanical loading, diet, inflammation and oxidative insult can all influence clock function. These factors may contribute to causing the change in the chondrocyte clock in OA.

Elucidating mechano-pathology of osteoarthritis: transcriptome-wide differences in mechanically stressed aged human cartilage explants

Background

Failing of intrinsic chondrocyte repair after mechanical stress is known as one of the most important initiators of osteoarthritis. Nonetheless, insight into these early mechano-pathophysiological processes in age-related human articular cartilage is still lacking. Such insights are needed to advance clinical development. To highlight important molecular processes of osteoarthritis mechano-pathology, the transcriptome-wide changes following injurious mechanical stress on human aged osteochondral explants were characterized.

Methods

Following mechanical stress at a strain of 65% (65%MS) on human osteochondral explants (n_65%MS = 14 versus n_control = 14), RNA sequencing was performed. Differential expression analysis between control and 65%MS was performed to determine mechanical stress-specific changes. Enrichment for pathways and protein-protein interactions was analyzed with Enrichr and STRING.

Results

We identified 156 genes significantly differentially expressed between control and 65%MS human osteochondral explants. Of note, IGFBP5 (FC = 6.01; FDR = 7.81 × 10⁻³) and MMP13 (FC = 5.19; FDR = 4.84 × 10⁻²) were the highest upregulated genes, while IGFBP6 (FC = 0.19; FDR = 3.07 × 10⁻⁴) was the most downregulated gene. Protein-protein interactions were significantly higher than expected by chance (P = 1.44 × 10⁻¹⁵ with connections between 116 out of 156 genes). Pathway analysis showed, among others, enrichment for cellular senescence, insulin-like growth factor (IGF) I and II binding, and focal adhesion.

Conclusions

Our results faithfully represent transcriptomic wide consequences of mechanical stress in human aged articular cartilage with MMP13, IGF binding proteins, and cellular senescence as the most notable results. Acquired knowledge on the as such identified initial, osteoarthritis-related, detrimental responses of chondrocytes may eventually contribute to the development of effective disease-modifying osteoarthritis treatments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13075-021-02595-8.

Linking chondrocyte and synovial transcriptional profile to clinical phenotype in osteoarthritis

OBJECTIVES

To determine how gene expression profiles in osteoarthritis joint tissues relate to patient phenotypes and whether molecular subtypes can be reproducibly captured by a molecular classification algorithm.

METHODS

We analysed RNA sequencing data from cartilage and synovium in 113 osteoarthritis patients, applying unsupervised clustering and Multi-Omics Factor Analysis to characterise transcriptional profiles. We tested the association of the molecularly defined patient subgroups with clinical characteristics from electronic health records.

RESULTS

We detected two patient subgroups in low-grade cartilage (showing no/minimal degeneration, cartilage normal/softening only), with differences associated with inflammation, extracellular matrix-related and cell adhesion pathways. The high-inflammation subgroup was associated with female sex (OR 4.12, p=0.0024) and prescription of proton pump inhibitors (OR 4.21, p=0.0040). We identified two independent patient subgroupings in osteoarthritis synovium: one related to inflammation and the other to extracellular matrix and cell adhesion processes. A seven-gene classifier including MMP13, APOD, MMP2, MMP1, CYTL1, IL6 and C15orf48 recapitulated the main axis of molecular heterogeneity in low-grade knee osteoarthritis cartilage (correlation ρ=-0.88, p<10^-10) and was reproducible in an independent patient cohort (ρ=-0.85, p<10^-10).

CONCLUSIONS

These data support the reproducible stratification of osteoarthritis patients by molecular subtype and the exploration of new avenues for tailored treatments.

Identification and characterization of two consistent osteoarthritis subtypes by transcriptome and clinical data integration

OBJECTIVE

To identify OA subtypes based on cartilage transcriptomic data in cartilage tissue and characterize their underlying pathophysiological processes and/or clinically relevant characteristics.

METHODS

This study includes n = 66 primary OA patients (41 knees and 25 hips), who underwent a joint replacement surgery, from which macroscopically unaffected (preserved, n = 56) and lesioned (n = 45) OA articular cartilage were collected [Research Arthritis and Articular Cartilage (RAAK) study]. Unsupervised hierarchical clustering analysis on preserved cartilage transcriptome followed by clinical data integration was performed. Protein-protein interaction (PPI) followed by pathway enrichment analysis were done for genes significant differentially expressed between subgroups with interactions in the PPI network.

RESULTS

Analysis of preserved samples (n = 56) resulted in two OA subtypes with n = 41 (cluster A) and n = 15 (cluster B) patients. The transcriptomic profile of cluster B cartilage, relative to cluster A (DE-AB genes) showed among others a pronounced upregulation of multiple genes involved in chemokine pathways. Nevertheless, upon investigating the OA pathophysiology in cluster B patients as reflected by differentially expressed genes between preserved and lesioned OA cartilage (DE-OA-B genes), the chemokine genes were significantly downregulated with OA pathophysiology. Upon integrating radiographic OA data, we showed that the OA phenotype among cluster B patients, relative to cluster A, may be characterized by higher joint space narrowing (JSN) scores and low osteophyte (OP) scores.

CONCLUSION

Based on whole-transcriptome profiling, we identified two robust OA subtypes characterized by unique OA, pathophysiological processes in cartilage as well as a clinical phenotype. We advocate that further characterization, confirmation and clinical data integration is a prerequisite to allow for development of treatments towards personalized care with concurrently more effective treatment response.

Flightless I is a catabolic factor of chondrocytes that promotes hypertrophy and cartilage degeneration in osteoarthritis

Synovial macrophages that are activated by cartilage fragments initiate synovitis, a condition that promotes hypertrophic changes in chondrocytes leading to cartilage degeneration in OA. In this study, we analyzed the molecular response of chondrocytes under condition of this type of stimulation to identify a molecular therapeutic target. Stimulated macrophages promoted hypertrophic changes in chondrocytes resulting in production of matrix-degrading enzymes of cartilage. Among the top-upregulated genes, FliI was found to be released from activated chondrocytes and exerted autocrine/paracrine effects on chondrocytes leading to an increase in expression of catabolic and hypertrophic factors. Silencing FliI in stimulated cells significantly reduced expression of catabolic and hypertrophic factors in cocultured chondrocytes. Our further results demonstrated that the FliI-TLR4-ERK1/2 axis is involved in the hypertrophic signaling of chondrocytes and catabolism of cartilage. Our findings provide a new insight into the pathogenesis of OA and identify a potentially new molecular target for diagnostics and therapeutics.

•

Activated macrophages promote the secretion of FliI from chondrocytes

•

FliI acts as a DAMP-triggering molecule in cartilage

•

FliI promotes chondrocyte hypertrophy and cartilage catabolism

•

FliI represents attractive target for therapeutic intervention

Two reactive behaviors of chondrocytes in an IL-1β-induced inflammatory environment revealed by the single-cell RNA sequencing

Objective: To investigate the heterogeneous responses of in vitro expanded chondrocytes, which were cultured in an interleukin (IL)-1β -induced inflammatory environment.

Method: Human articular chondrocytes were expanded, in vitro, for 13 days and treated with IL-1β for 0, 24, and 48 h. Cells were collected and subjected to single-cell RNA sequencing. Multiple bioinformatics tools were used to determine the signatures that define chondrocyte physiology.

Results: Two major cell clusters with distinct expression patterns were identified at the initial phase and were with heterogeneous variation that coincides with inflammation progress. They transformed into two terminal cell clusters one of which exhibited OA-phenotype and proinflammatory characteristics through two paths, “response-to-inflammation” and “atypical response-to-inflammation”, respectively. The involved cell clusters exhibited intrinsic relationship with cell types within native cartilage from OA patients. Genes controlling cell transformation to OA-phenotype were relating to the tumor necrosis factor (TNF) signaling pathway via NFKB, up-regulated KRAS signaling and the IL2/STAT5 signaling pathway and pathways relating to apoptosis and reactive oxygen species.

Conclusion: The in vitro expanded chondrocytes under IL-1β-induced inflammatory progression behave heterogeneously. One of the initial cell clusters could transform into a proinflammatory subpopulation through a termed response-to-inflammation path, which may serve as the core target to alleviate OA progression.

Identification of key regulators responsible for dysregulated networks in osteoarthritis by large-scale expression analysis

BACKGROUND

Osteoarthritis (OA) is a worldwide musculoskeletal disorder. However, disease-modifying therapies for OA are not available. Here, we aimed to characterize the molecular signatures of OA and to identify novel therapeutic targets and strategies to improve the treatment of OA.

METHODS

We collected genome-wide transcriptome data performed on 132 OA and 74 normal human cartilage or synovium tissues from 7 independent datasets. Differential gene expression analysis and functional enrichment were performed to identify genes and pathways that were dysregulated in OA. The computational drug repurposing method was used to uncover drugs that could be repurposed to treat OA.

RESULTS

We identified several pathways associated with the development of OA, such as extracellular matrix organization, inflammation, bone development, and ossification. By protein-protein interaction (PPI) network analysis, we prioritized several hub genes, such as JUN, CDKN1A, VEGFA, and FOXO3. Moreover, we repurposed several FDA-approved drugs, such as cardiac glycosides, that could be used in the treatment of OA.

CONCLUSIONS

We proposed that the hub genes we identified would play a role in cartilage homeostasis and could be important diagnostic and therapeutic targets. Drugs such as cardiac glycosides provided new possibilities for the treatment of OA.

[Possibilities and limits of conservative treatment for osteoarthritis : Sport advice, training therapy, orthotics and cartilage therapeutics]

Hintergrund

Arthrose – die Degeneration von Gelenken – ist ein weit verbreitetes Problem durch alle Bevölkerungsschichten, das im zunehmenden Alter vermehrt auftritt und die häufigste Ursache für mobilitätseinschränkende Schmerzen am Bewegungsapparat ist. Etwa 70–80 % der über 70-Jährigen zeigen Zeichen einer Gelenksdegeneration. Insgesamt sind bis zu 25 % der Gesamtbevölkerung davon betroffen, aufgrund der generellen Alterung der Bevölkerung mit steigender Tendenz. Die Inzidenz der Arthrose steigt aber schon ab dem 40 Lebensjahr, wobei besonders posttraumatische und sekundäre Arthroseformen zum Tragen kommen.

Anspruch

Der Wunsch nach hoher Mobilität und Sport zieht sich als Phänomen ebenfalls durch alle Altersgruppe. Dies ist mit hohen Gelenkbelastungen verbunden und stellt damit eine große Herausforderung an vor allem früh degenerativ veränderte Gelenksstrukturen dar. In diesem Zusammenhang ist der orthopädisch tätige Arzt gefordert, die Belastbarkeit von geschädigten Gelenken abzuschätzen und so früh wie möglich präventive Schritte sowie gegebenenfalls konservative Therapien einzuleiten, um die Progression der Arthrose zu verhindern und damit den eventuell notwendigen Gelenkersatz möglichst weit nach hinten zu schieben.

Regulation of TGFβ Signalling by TRPV4 in Chondrocytes

The growth factor TGFβ and the mechanosensitive calcium-permeable cation channel TRPV4 are both important for the development and maintenance of many tissues. Although TRPV4 and TGFβ both affect core cellular functions, how their signals are integrated is unknown. Here we show that pharmacological activation of TRPV4 significantly increased the canonical response to TGFβ stimulation in chondrocytes. Critically, this increase was only observed when TRPV4 was activated after, but not before TGFβ stimulation. The increase was prevented by pharmacological TRPV4 inhibition or knockdown and is calcium/CamKII dependent. RNA-seq analysis after TRPV4 activation showed enrichment for the TGFβ signalling pathway and identified JUN and SP1 as key transcription factors involved in this response. TRPV4 modulation of TGFβ signalling represents an important pathway linking mechanical signalling to tissue development and homeostasis.