A gene expression study of normal and damaged cartilage in anteromedial gonarthrosis, a phenotype of osteoarthritis

A versatile platform based on matrix metalloproteinase-sensitive peptides for novel diagnostic and therapeutic strategies in arthritis

Matrix metalloproteinases (MMPs), coupled with other proteinases and glycanases, can degrade proteoglycans, collagens, and other extracellular matrix (ECM) components in inflammatory and non-inflammatory arthritis, making them important pathogenic molecules and ideal disease indicators and pharmaceutical intervention triggers. For MMP responsiveness, MMP-sensitive peptides (MSPs) are among the most easily synthesized and cost-effective substrates, with free terminal amine and/or carboxyl groups extensively employed in multiple designs. We hereby provide a comprehensive review over the mechanisms and advances in MSP applications for the management of arthritis. These applications include early and precise diagnosis of MMP activity via fluorescence probe technologies; acting as nanodrug carriers to enable on-demand drug release triggered by pathological microenvironments; and facilitating cartilage engineering through MMP-mediated degradation, which promotes cell migration, matrix synthesis, and tissue integration. Specifically, the ultra-sensitive MSP diagnostic probes could significantly advance the early diagnosis and detection of osteoarthritis (OA), while MSP-based drug carriers for rheumatoid arthritis (RA) can intelligently release anti-inflammatory drugs effectively during flare-ups, or even before symptoms manifest. The continuous progress in MSP development may acceleratedly lead to novel management regimens for arthropathy in the future.

Cell culture expansion media choice affects secretory, protective and immuno-modulatory features of adipose mesenchymal stromal cell-derived secretomes for orthopaedic applications

Introduction

Mesenchymal stromal cells (MSCs) gained attention for their anti-inflammatory and trophic properties, with musculoskeletal diseases and osteoarthritis (OA) being among the most studied conditions. Alongside cells, their released factors and extracellular vesicles (EVs), overall termed "secretome", are actively sifted being envisioned as the main therapeutic actors. In addition to standard supplementation given by foetal bovine serum (FBS) or human platelet lysate (hPL), new good manufacturing practice (GMP)-compliant serum/xeno (S/X)-free media formulations have been proposed, although their influence on MSCs phenotype and potential is scarcely described. The aim of this study is therefore to evaluate, in the OA context, the differences in secretome composition and potential after adipose-MSCs (ASCs) cultivation in both standard (FBS and hPL) and two next generation (S/X) GMP-ready supplements.

Methods

Immunophenotype and secretory ability at soluble protein and EV-related levels, including embedded miRNAs, were analysed in the secretomes by means of flow cytometry, nanoparticle tracking analysis, high throughput ELISA and qRT-PCR arrays. Secretomes effect was tested in in vitro models of chondrocytes, lymphocytes and monocytes to mimic the OA microenvironment.

Results

Within a conserved molecular signature, a divergent fingerprint emerged for ASCs' secretomes collected after expansion in standard FBS/hPL or next-generation S/X formulations. Regarding soluble factors, a less protective feature for those in the secretome collected after ASCs were cultured in S/X media emerged. Moreover, the overall message for EV-miRNAs was characterized by a preponderance of protective signals in FBS and hPL conditions in a context of general safeguard given by ASCs released molecules. This dichotomy was reflected on secretomes' potential in vitro, with expansion in hPL resulting in the most effective secretome for chondrocytes and in FBS for immune cells.

Conclusions

These data open the question about the implications from using new media for MSCs expansion for clinical application. Although the undeniable advantages for GMP compliant processes, this study results suggest that new media formulations would deserve a deep characterization to drive the choice of the most effective one tailored to each specific application.

Dipeptidyl Peptidase 4 (DPP4) Exacerbates Osteoarthritis Progression in an Enzyme-Independent Manner

Chondrocyte senescence is a key driver of osteoarthritis (OA). Mitochondrial dysfunction and oxidative stress can induce chondrocyte senescence. However, the specific mechanisms by which senescence contributes to OA progression are not fully understood. Here, it is attested that Dipeptidyl peptidase 4 (DPP4) is significantly upregulated in osteoarthritic chondrocytes in both humans and mice. DPP4 promotes oxidative stress and cellular senescence in chondrocytes through excessive mitochondrial fission in an enzyme-independent manner. Intra-articular injection of adeno-associated virus 2 to upregulate DPP4 in chondrocytes promotes post-traumatic and aging-induced OA in mice in an enzyme-independent manner. Mechanistically, DPP4 competitively binds to Myosin heavy chain 9 (MYH9), interfering with its E3 ubiquitin ligase Carboxyl terminus of Hsc70-interacting protein (CHIP), and thereby upregulates MYH9 expression. Finally, a small molecule, 4,5-Dicaffeoylquinic acid is identified, which disrupts the interaction between DPP4 and MYH9, thereby ameliorating post-traumatic and aging-induced OA in mice caused by DPP4 upregulation. The study indicates that the non-enzymatic activity of DPP4 is a promising target for OA treatment.

Primary cartilage transcriptional signatures reflect cell-type-specific molecular pathways underpinning osteoarthritis

Translational efforts in osteoarthritis are hampered by a gap in our understanding of disease processes at the molecular level. Here, we present evidence of pronounced transcriptional changes in high- and low-disease-grade cartilage tissue, pointing to embryonic processes involved in disease progression. We identify shared transcriptional programs between osteoarthritis cartilage and cell populations in the human embryonic and fetal limb, pointing to increases in pre-hypertrophic chondrocytes' transcriptional programs in low-grade cartilage and increases in osteoblastic signatures in high-grade disease tissue. We find that osteoarthritis genetic risk signals are enriched in six gene co-expression modules and show that these transcriptional signatures reflect cell-type-specific expression along the endochondral ossification developmental trajectory. Using this network approach in combination with causal inference analysis, we present evidence of a causal effect on osteoarthritis risk for variants associated with the expression of ten genes that have not been previously reported as effector genes in genome-wide association studies in osteoarthritis. Our findings point to key molecular pathways as drivers of cartilage degeneration and identify high-value drug targets and repurposing opportunities.

Dissecting SOX9 dynamics reveals its differential regulation in osteoarthritis

The transcription factor SOX9 is integral to tissue homeostasis and is implicated in skeletal malformation, campomelic dysplasia, and osteoarthritis (OA). Despite extensive research, the complete regulatory landscape of SOX9 transcriptional activity, interconnected with signaling pathways (TGFβ, WNT, BMP, IHH, NFκB, and HIF), remains challenging to decipher. This study focuses on elucidating SOX9 signaling in OA pathology using Fluorescence Recovery After Photobleaching (FRAP) to assess SOX9 activity directly in live human primary chondrocytes (hPCs). Single cell FRAP data revealed two distinct subpopulations with differential SOX9 dynamics, showing varied distribution between healthy and OA hPCs. Moreover, inherently elevated SOX9-DNA binding was observed in healthy hPCs compared to preserved and OA counterparts. Anabolic factors (BMP7 and GREM1) and catabolic inhibitors (DKK1 and FRZb) were found to modulate SOX9 transcriptional activity in OA-hPCs. These findings provide valuable insights into the intricate regulation of SOX9 signaling in OA, suggesting potential therapeutic avenues for modulating SOX9 activity in diseased states.

Synoviocytes assist in modulating the effect of Ross River virus infection in micromass-cultured primary human chondrocytes

Introduction. Ross River virus (RRV) is a mosquito-borne virus prevalent in Australia and the islands of the South Pacific, where it causes an arthritogenic illness with a hallmark feature of severe joint pain. The joint space is a unique microenvironment that contains cartilage and synovial fluid. Chondrocytes and synoviocytes are crucial components of the joint space and are known targets of RRV infection.Hypothesis/Gap statement. Understanding the relationship between synoviocytes and chondrocytes during RRV infection will provide further insights into RRV-induced joint pathology.Methodology. To better understand the unique dynamics of these cells during RRV infection, we used primary chondrocytes cultured in physiologically relevant micromasses. We then directly infected micromass chondrocytes or infected primary fibroblast-like synoviocytes (FLS), co-cultured with micromass chondrocytes. Micromass cultures and supernatants were collected and analysed for viral load with a PCR array of target genes known to play a role in arthritis.Results. We show that RRV through direct or secondary infection in micromass chondrocytes modulates the expression of cellular factors that likely contribute to joint inflammation and disease pathology, as well as symptoms such as pain. More importantly, while we show that RRV can infect micromass-cultured chondrocytes via FLS infection, FLS themselves affect the regulation of cellular genes known to contribute to arthritis.Conclusion. Single-cell culture systems lack the complexity of in vivo systems, and understanding the interaction between cell populations is crucial for deciphering disease pathology, including for the development of effective therapeutic strategies.

Inflammatory biomarkers and state of the tibiofemoral joint in the osteoarthritic knee: a narrative review

Background

The healing process is initiated by injurious stimuli in response to cellular damage. Upon recruiting proinflammatory biomarkers to the tissue site of injury, the release of additional biomarkers occurs, including the likes of cytokines, matrix molecules, macrophages, neutrophils, and others. This influx of immune system mediators can occur for chronic periods, and though its intention is for healing the original injurious stimuli, it is also suspected of causing long term cartilage impairment following internal structure damage. The objective of this narrative review is to identify which inflammatory factors have the leading roles in the progression of osteoarthritis (OA) following knee injuries and how they fluctuate throughout the healing process, both acutely and chronically.

Methods

This narrative review was performed following a computerized search of the electronic database on PubMed in May 2023. Abstracts related to the inflammatory biomarkers of the post-traumatic knee were included for review.

Key Content and Findings

The chronic low-level inflammation that leads to OA leads to the destruction of the cartilage extracellular matrix, which new and developing orthopedic research is still attempting to find resolve for. Some of this damage is attributed to the biomechanical alterations that occurs following injury, though with most procedures capable of joint biomechanical restoration, focus has rather been shifted toward the environment of inflammatory biomarkers.

Conclusions

Future studies will be aiming to improve the diagnostics of OA, focusing on a consistent correlation of inflammatory biomarkers with imaging. Additionally, biochemical treatments will need to focus on validating reproducible modulation of signaling molecules, in attempts to lessen the chronic elevations of destructive biomarkers.

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.

Cyaonoside A-loaded composite hydrogel microspheres to treat osteoarthritis by relieving chondrocyte inflammation

Cyaonoside A (CyA), derived from the natural Chinese medicine, Cyathula officinalis Kuan, which was for a long time used to treat knee injuries and relieve joint pain in traditional Chinese medicine, showed an unclear mechanism for protecting cartilage. In addition, CyA was poorly hydrosoluble and incapable of being injected directly into the joint cavity, which limited its clinical application. This study reveals that CyA resisted IL-1β-mediated chondrogenic inflammation and apoptosis. Next, transcriptome sequencing is used to explore the potential mechanisms underlying CyA regulation of MSC chondrogenic differentiation. Based on these findings, CyA-loaded composite hydrogel microspheres (HLC) were developed and they possessed satisfactory loading efficiency, a suitable degradation rate and good biocompatibility. HLC increased chondrogenic anabolic gene (Acan, COL2A, and SOX9) expression, while downregulating the expression of the catabolic marker MMP13 in vitro. In the osteoarthritis mouse model, HLC demonstrated promising therapeutic capabilities by protecting the integrity of articular cartilage. In conclusion, this study provides insights into the regulatory mechanisms of CyA for chondrocytes and proposes a composite hydrogel microsphere-based advanced therapeutic strategy for osteoarthritis.

DDX5 inhibits hyaline cartilage fibrosis and degradation in osteoarthritis via alternative splicing and G-quadruplex unwinding

Hyaline cartilage fibrosis is typically considered an end-stage pathology of osteoarthritis (OA), which results in changes to the extracellular matrix. However, the mechanism behind this is largely unclear. Here, we found that the RNA helicase DDX5 was dramatically downregulated during the progression of OA. DDX5 deficiency increased fibrosis phenotype by upregulating COL1 expression and downregulating COL2 expression. In addition, loss of DDX5 aggravated cartilage degradation by inducing the production of cartilage-degrading enzymes. Chondrocyte-specific deletion of Ddx5 led to more severe cartilage lesions in the mouse OA model. Mechanistically, weakened DDX5 resulted in abundance of the Fn1-AS-WT and Plod2-AS-WT transcripts, which promoted expression of fibrosis-related genes (Col1, Acta2) and extracellular matrix degradation genes (Mmp13, Nos2 and so on), respectively. Additionally, loss of DDX5 prevented the unfolding Col2 promoter G-quadruplex, thereby reducing COL2 production. Together, our data suggest that strategies aimed at the upregulation of DDX5 hold significant potential for the treatment of cartilage fibrosis and degradation in OA.

Articular cartilage corefucosylation regulates tissue resilience in osteoarthritis

This study aimed to investigate the glycan structural changes that occur before histological degeneration in osteoarthritis (OA) and to determine the mechanism by which these glycan conformational changes affect cartilage degeneration. An OA model was established in rabbits using mannosidase injection, which reduced high-mannose type N-glycans and led to cartilage degeneration. Further analysis of glycome in human OA cartilage identified specific corefucosylated N-glycan expression patterns. Inhibition of N-glycan corefucosylation in mice resulted in unrecoverable cartilage degeneration, while cartilage-specific blocking of corefucosylation led to accelerated development of aging-associated and instability-induced OA models. We conclude that α1,6 fucosyltransferase is required postnatally to prevent preosteoarthritic deterioration of articular cartilage. These findings provide a novel definition of early OA and identify glyco-phenotypes of OA cartilage, which may distinguish individuals at higher risk of progression.

Characterization of adult human skeletal cells in different tissues reveals a CD90+CD34+ periosteal stem/progenitor population

The periosteum plays a crucial role in bone healing and is an important source of skeletal stem and progenitor cells. Recent studies in mice indicate that diverse populations of skeletal progenitors contribute to growth, homeostasis and healing. Information about the in vivo identity and diversity of skeletal stem and progenitor cells in different compartments of the adult human skeleton is limited. In this study, we compared non-hematopoietic populations in matched tissues from the femoral head and neck of 21 human participants using spectral flow cytometry of freshly isolated cells. High-dimensional clustering analysis indicated significant differences in marker distribution between periosteum, articular cartilage, endosteum and bone marrow populations, and identified populations that were highly enriched or unique to specific tissues. Periosteum-enriched markers included CD90 and CD34. Articular cartilage, which has very poor regenerative potential, showed enrichment of multiple markers, including the PDPN+CD73+CD164+CD146- population previously reported to represent human skeletal stem cells. We further characterized periosteal populations by combining CD90 with other strongly expressed markers. CD90+CD34+ cells sorted directly from periosteum showed significant colony-forming unit fibroblasts (CFU-F) enrichment, rapid expansion, and consistent multi-lineage differentiation of clonal populations in vitro. In situ, CD90+CD34+ cells include a perivascular population in the outer layer of the periosteum and non-perivascular cells closer to the bone surface. CD90+ cells are also highly enriched for CFU-F in bone marrow and endosteum, but not articular cartilage. In conclusion, our study indicates considerable diversity in the non-hematopoietic cell populations in different tissue compartments within the adult human skeleton, and suggests that periosteal progenitor cells reside within the CD90+CD34+ population.

Is the Synovium the First Responder to Posttraumatic Knee Joint Stress? The Molecular Pathogenesis of Traumatic Cartilage Degeneration

OBJECTIVE

The aim of this study was to evaluate if a similar catabolic and inflammatory gene pattern exists between the synovium, hyaline cartilage, and blood of patients with the knee joint tissues and if one precedes the other.

DESIGN

A total of fifty-eight patients (34 females and 24 males) with a mean age of 44.7 years (range, 18-75) underwent elective knee arthroscopy due to previously diagnosed pathology. Full blood samples were collected preoperatively from synovium and cartilage samples intraoperatively. Real time PCR with spectrophotometric analysis was performed. Following genes taking part in ECM (extracellular matrix) remodeling were selected for analysis: MMP-1, MMP-2, MMP-8, MMP-9, MMP-13, MMP-14, ADAMTS-4 (Agg1) and ADAMTS-5 (Agg2) proteases, TIMP-1, and TIMP-2 - their inhibitors - and IL-1 and TNF-α cytokines.

RESULTS

Analysis revealed a strong and significant correlation between gene expression in synovial and systemic blood cells (p <0.05 for all studied genes) with ADAMTS-4, ADAMTS-5, IL-1, TNF-α and TIMP-2 expression most positively correlated with an R>0.8 for each. An analysis between chondrocytes and systemic blood gene expression shown no significant correlation for all genes. Bivariate correlation of International Cartilage Repair Society grading and genes expression revealed significant associations with synovial MMP-1, MMP-2, MMP-8, MMP-9, IL-1, TNF-α and TIMP-2.

CONCLUSION

We suggest that the synovial tissue is the first responder for knee joint stress factors in correlation with the response of blood cells. The chondrocyte's genetic response must be further investigated to elucidate the genetic program of synovial joints, as an organ, during OA development and progression.

Microscopic and transcriptomic changes in porcine synovium one year following disruption of the anterior cruciate ligament

OBJECTIVE

There is no disease-modifying treatment for posttraumatic osteoarthritis (PTOA). This may be partly due to an incomplete understanding of synovitis, which has been causally linked to PTOA progression. The microscopic and transcriptomic changes in synovium seen in early- to mid-stage PTOA were evaluated to better characterize this knowledge gap.

METHODS

Seventy-two Yucatan minipigs underwent transection of the anterior cruciate ligament (ACL). Subjects were randomized to no further intervention, ligament reconstruction, or ligament repair, followed by microscopic synovium evaluation and RNA-sequencing at 1, 4, and 52 weeks. Six additional subjects received no ligament transection and served as 1- and 4-week controls and 12 contralateral knees served as 52-week controls.

RESULTS

Synovial lining thickness, stromal cellularity, and overall microscopic synovitis reached their highest levels in the first few weeks following injury. Inflammatory infiltration continued to increase over the course of a year. Leaving the ACL transected, reconstructing the ligament, or repairing the ligament did not modulate synovitis development at 1, 4, or 52 weeks. Differential gene expression analysis of PTOA-affected synovium compared to control synovium revealed increased cell proliferation, angiogenesis, collagen breakdown, and diminished lipid metabolism at 1 and 4 weeks, and increased axonogenesis and focal adhesion with reduced immune activation at 52 weeks.

CONCLUSIONS

Synovitis was present one year after ACL injury and was not alleviated by surgical intervention. Gene expression in early synovitis was characterized by cell proliferation, angiogenesis, proteolysis, and reduced lipolysis, which was followed by nerve growth and cellular adhesion with less immune activation at 52 weeks.

The recombinant Link module of human TSG-6 suppresses cartilage damage in models of osteoarthritis: A potential disease-modifying OA drug

OBJECTIVE

To investigate the role of endogenous TSG-6 in human osteoarthritis (OA) and assess the disease-modifying potential of a TSG-6-based biological treatment in cell, explant and animal models of OA.

DESIGN

Knee articular cartilages from OA patients were analyzed for TSG-6 protein and mRNA expression using immunohistochemistry and RNAscope, respectively. The inhibitory activities of TSG-6 and its isolated Link module (Link_TSG6) on cytokine-induced degradation of OA cartilage explants were compared. Human mesenchymal stem/stromal cell-derived chondrocyte pellet cultures were used to determine the effects of Link_TSG6 and full-length TSG-6 on IL-1α-, IL-1β-, or TNF-stimulated ADAMTS4, ADAMTS5, and MMP13 mRNA expression. Link_TSG6 was administered i.a. to the rat ACLTpMMx model; cartilage damage and tactile allodynia were assessed.

RESULTS

TSG-6 is predominantly associated with chondrocytes in regions of cartilage damage where high TSG-6 expression aligns with low MMP13, the major collagenase implicated in OA progression. Link_TSG6 is more potent than full-length TSG-6 at inhibiting cytokine-mediated matrix breakdown in human OA cartilage explants;>50% of donor cartilages, from 59 tested, were responsive to Link_TSG6 treatment. Link_TSG6 also displayed more potent effects in 3D pellet cultures, suppressing ADAMTS4, ADAMTS5, and MMP13 gene expression, which was consistent with reduced aggrecanase and collagenase activities in explant cultures. Link_TSG6 treatment reduced touch-evoked pain behavior and dose-dependently inhibited cartilage damage in a rodent model of surgically-induced OA.

CONCLUSIONS

Link_TSG6 has enhanced chondroprotective activity compared to the full-length TSG-6 protein and shows potential as a disease modifying OA drug via its inhibition of aggrecanase and collagenase activity.

REST, RCOR1 and RCOR2 expression is reduced in osteoarthritic chondrocytes and contributes to increasing MMP13 and ADAMTS5 expression through upregulating HES1

Expression of key transcriptional regulators is altered in chondrocytes in osteoarthritis (OA). This contributes to an increase in production of cartilage-catabolizing enzymes such as MMP13 and ADAMTS5. RCOR1 and RCOR2, binding partners for the transcriptional repressor REST, have previously been found to be downregulated in OA chondrocytes although their function in chondrocytes is unclear. HES1 is a known REST/RCOR1 target gene and HES1 has been shown to promote MMP13 and ADAMTS5 expression in murine OA chondrocytes. The purpose of this study was to determine whether reduced REST/RCOR levels leads to increased HES1 expression in human OA chondrocytes and whether HES1 also promotes ADAMTS5 and MMP13 expression in these cells. Chondrocytes were isolated from osteoarthritic and adjacent macroscopically normal cartilage obtained from patients undergoing total knee arthroplasty. RNA and protein levels of REST, RCOR1 and RCOR2 were lower, but levels of HES1 higher, in chondrocytes isolated from osteoarthritic compared to macroscopically normal cartilage. Over-expression of either REST, RCOR1 or RCOR2 resulted in reduced HES1 levels in OA chondrocytes whereas knockdown of REST, RCOR1 or RCOR2 led to increased HES1 expression in chondrocytes from macroscopically normal cartilage. In OA chondrocytes, ADAMTS5 and MMP13 expression were reduced following HES1 knockdown, but further enhanced following HES1 over-expression. Levels of phosphorylated CaMKII were higher in chondrocytes from OA cartilage consistent with previous findings that HES1 only promotes gene transcription in the presence of active CaMKII. These findings identify the REST/RCOR/HES1 pathway as a contributing factor leading to increased ADAMTS5 and MMP13 expression in OA chondrocytes.

Differential Effects of Hypoxia versus Hyperoxia or Physoxia on Phenotype and Energy Metabolism in Human Chondrocytes from Osteoarthritic Compared to Macroscopically Normal Cartilage

Chondrocyte phenotype and energy metabolism are altered in osteoarthritis (OA). However, most studies characterising the change in human chondrocyte behaviour in OA have been conducted in supraphysiological oxygen concentrations. The purpose of this study was to compare phenotype and energy metabolism in chondrocytes from macroscopically normal (MN) and OA cartilage maintained in 18.9% (standard tissue culture), 6% (equivalent to superficial zone of cartilage in vivo) or 1% oxygen (equivalent to deep zone of cartilage in vivo). MMP13 production was higher in chondrocytes from OA compared to MN cartilage in hyperoxia and physoxia but not hypoxia. Hypoxia promoted SOX9, COL2A1 and ACAN protein expression in chondrocytes from MN but not OA cartilage. OA chondrocytes used higher levels of glycolysis regardless of oxygen availability. These results show that differences in phenotype and energy metabolism between chondrocytes from OA and MN cartilage differ depending on oxygen availability. OA chondrocytes show elevated synthesis of cartilage-catabolising enzymes and chondrocytes from MN cartilage show reduced cartilage anabolism in oxygenated conditions. This is relevant as a recent study has shown that oxygen levels are elevated in OA cartilage in vivo. Our findings may indicate that this elevated cartilage oxygenation may promote cartilage loss in OA.

Acute loading has minor influence on human articular cartilage gene expression and glycosaminoglycan composition in late-stage knee osteoarthritis: a randomised controlled trial

OBJECTIVE

Osteoarthritis (OA) remains clinically challenging. Regular physical exercise improves symptoms though it is unclear whether exercise influences cartilage at the molecular level. Thus, we aimed to determine the effect of acute loading on gene expression and glycosaminoglycan (GAG) content in human OA cartilage.

DESIGN

Patients with primary knee OA participated in this single-blind randomised controlled trial initiated 3.5 h prior to scheduled joint replacement surgery with or without loading by performing one bout of resistance exercise (one-legged leg press). Cartilage from the medial tibia condyle was sampled centrally, under the meniscus, and from peripheral osteophytes. Samples were analysed for gene expression by real-time reverse transcriptase polymerase chain reaction, and hyaluronidase-extracted matrix was analysed for GAG composition by immuno- and dimethyl-methylene blue assays.

RESULTS

Of 32 patients randomised, 31 completed the intervention: mean age 69 ± 7.5 years (SD), 58% female, BMI 29.4 ± 4.4 kg/m². Exercise increased chondroitin sulphate extractability [95% CI: 1.01 to 2.46; P = 0.0486] but cartilage relevant gene expression was unchanged. Regionally, the submeniscal area showed higher MMP-3, MMP-13, IGF-1Ea, and CTGF, together with lower lubricin and COMP expression compared to the central condylar region. Further, osteophyte expression of MMP-1, MMP-13, IGF-1Ea, and TGF-β3 was higher than articular cartilage and lower for aggrecan, COMP, and FGF-2. Hyaluronidase-extracted matrix from central condylar cartilage contained more GAGs but less chondroitin sulphate compared to submeniscal cartilage.

CONCLUSION

Acute exercise had minor influence on cartilage GAG dynamics, indicating that osteoarthritic cartilage is not significantly affected by acute exercise. However, the regional differences suggest a chronic mechanical influence on human cartilage.

CLINICALTRIALS

GOV REGISTRATION NUMBER

NCT03410745.

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 regional turnover of cartilage collagen matrix in late-stage human knee osteoarthritis

OBJECTIVE

Cartilage collagen has very limited repair potential, though some turnover and incorporation has not been fully excluded. We aim to determine the regional turnover of human osteoarthritis cartilage.

DESIGN

Patients scheduled for knee joint replacement surgery due to osteoarthritis were recruited in this prospective study of four weeks duration. Deuterium oxide (D₂O) was administered orally by weekly boluses at 70% D₂O, initially 150 ml followed by three boluses of 50 ml. Cartilage from the medial tibia plateau was sampled centrally, under the meniscus, and from osteophytes and treated enzymatically with hyaluronidase and trypsin. Samples were analysed for deuterium incorporation in alanine using mass spectrometry and for gene expression by real-time reverse transcriptase polymerase chain reaction.

RESULTS

Twenty participants completed the study: mean (SD) age 64 ± 9.1 years, 45% female, BMI 29.5 ± 4.8 kg/m². Enzymatically treated cartilage from central and submeniscal regions showed similar enrichments at 0.063% APE, while osteophytes showed significantly greater enrichment at 0.072% APE (95% confidence interval of difference) [0.004-0.015]). Fractional synthesis rates were similar for central 0.027%/day and submeniscal cartilage 0.022%/day but 10-fold higher in osteophytes 0.22%/day [0.098-0.363]. When compared to central cartilage, submeniscal cartilage had increased gene expression of MMP-3 and decreased lubricin expression. Untreated cartilage had higher turnover (enrichments at 0.073% APE) than enzymatically treated cartilage (0.063% APE).

CONCLUSIONS

In OA, despite regional differences in gene expression, the turnover of the articular cartilage matrix across the entire joint surface is very limited, but higher turnover was observed in osteophyte cartilage.

Protein tyrosine phosphatases in skeletal development and diseases

Skeletal development and homeostasis in mammals are modulated by finely coordinated processes of migration, proliferation, differentiation, and death of skeletogenic cells originating from the mesoderm and neural crest. Numerous molecular mechanisms are involved in these regulatory processes, one of which is protein posttranslational modifications, particularly protein tyrosine phosphorylation (PYP). PYP occurs mainly through the action of protein tyrosine kinases (PTKs), modifying protein enzymatic activity, changing its cellular localization, and aiding in the assembly or disassembly of protein signaling complexes. Under physiological conditions, PYP is balanced by the coordinated action of PTKs and protein tyrosine phosphatases (PTPs). Dysregulation of PYP can cause genetic, metabolic, developmental, and oncogenic skeletal diseases. Although PYP is a reversible biochemical process, in contrast to PTKs, little is known about how this equilibrium is modulated by PTPs in the skeletal system. Whole-genome sequencing has revealed a large and diverse superfamily of PTP genes (over 100 members) in humans, which can be further divided into cysteine (Cys)-, aspartic acid (Asp)-, and histidine (His)-based PTPs. Here, we review current knowledge about the functions and regulatory mechanisms of 28 PTPs involved in skeletal development and diseases; 27 of them belong to class I and II Cys-based PTPs, and the other is an Asp-based PTP. Recent progress in analyzing animal models that harbor various mutations in these PTPs and future research directions are also discussed. Our literature review indicates that PTPs are as crucial as PTKs in supporting skeletal development and homeostasis.

Transcriptome Analysis Reveals Novel Genes Associated with Cartilage Degeneration in Posttraumatic Osteoarthritis Progression

OBJECTIVE

The current therapeutic strategy for posttraumatic osteoarthritis (PTOA) focuses on early intervention to attenuate disease progression, preserve joint function, and defer joint replacement timing. Sequential transcriptomic changes of articular cartilage in a rat model were investigated to explore the molecular mechanism in early PTOA progression.

DESIGN

Anterior cruciate ligament transection and medial meniscectomy (ACLT + MMx)-induced PTOA model was applied on male Wistar rats. Articular cartilages were harvested at time 0 (naïve), 2 week, and 4 weeks after surgery. Affymetrix Rat genome 230 2.0 array was utilized to analyze the gene expression changes of articular cartilages.

RESULTS

We identified 849 differentially expressed genes (DEGs) at 2 weeks and 223 DEGs at 4 weeks post-ACLT + MMx surgery compared with time 0 (naïve group). Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to gain further insights from these DEGs. 22 novel genes and 1 novel KEGG pathway (axon guidance) in cartilage degeneration of osteoarthritis were identified. Axon guidance molecules-Gnai1, Sema4d, Plxnb1, and Srgap2 commonly dysregulated in PTOA progression. Gnai1 gene showed a concordant change in protein expression by immunohistochemistry staining.

CONCLUSIONS

Our study identified 22 novel dysregulated genes and axon guidance pathway associated with articular cartilage degeneration in PTOA progression. These findings provide the potential candidates of biomarkers and therapeutic targets for further investigation.

General Caspase Inhibition in Primary Chondrogenic Cultures Impacts Their Transcription Profile Including Osteoarthritis-Related Factors

OBJECTIVE

The knowledge about functions of caspases, usually associated with cell death and inflammation, keeps expanding also regarding cartilage. Active caspases are present in the growth plate, and caspase inhibition in limb-derived chondroblasts altered the expression of osteogenesis-related genes. Caspase inhibitors were reported to reduce the severity of cartilage lesions in osteoarthritis (OA), and caspase-3 might represent a promising biomarker for OA prognosis. The objective of this investigation was to decipher the transcriptomic regulation of caspase inhibition in chondrogenic cells.

DESIGN

Limb-derived chondroblasts were cultured in the presence of 2 different inhibitors: Z-VAD-FMK (FMK) and Q-VD-OPH (OPH). A whole transcriptome RNA sequencing was performed as the key analysis.

RESULTS

The analysis revealed a statistically significant increase in the expression of 252 genes in the FMK samples and 163 genes in the OPH samples compared with controls. Conversely, there was a significant decrease in the expression of 290 genes in the FMK group and 188 in the OPH group. Among the top up- and downregulated genes (more than 10 times changed), almost half of them were associated with OA. Both inhibitors displayed the highest upregulation of the inflammatory chemokine Ccl5, the most downregulated gene was the one for mannose receptors Mrc1.

CONCLUSIONS

The obtained datasets pointed to a significant impact of caspase inhibition on the expression of several chondro-/osteogenesis-related markers in an in vitro model of endochondral ossification. Notably, the list of these genes included some encoding for factors associated with cartilage/bone pathologies such as OA.

A Synchronized Circadian Clock Enhances Early Chondrogenesis

OBJECTIVE

Circadian rhythms in cartilage homeostasis are hypothesized to temporally segregate and synchronize the activities of chondrocytes to different times of the day, and thus may provide an efficient mechanism by which articular cartilage can recover following physical activity. While the circadian clock is clearly involved in chondrocyte homeostasis in health and disease, it is unclear as to what roles it may play during early chondrogenesis.

DESIGN

The purpose of this study was to determine whether the rhythmic expression of the core circadian clock was detectable at the earliest stages of chondrocyte differentiation, and if so, whether a synchronized expression pattern of chondrogenic transcription factors and developing cartilage matrix constituents was present during cartilage formation.

RESULTS

Following serum shock, embryonic limb bud-derived chondrifying micromass cultures exhibited synchronized temporal expression patterns of core clock genes involved in the molecular circadian clock. We also observed that chondrogenic marker genes followed a circadian oscillatory pattern. Clock synchronization significantly enhanced cartilage matrix production and elevated SOX9, ACAN, and COL2A1 gene expression. The observed chondrogenesis-promoting effect of the serum shock was likely attributable to its synchronizing effect on the molecular clockwork, as co-application of small molecule modulators (longdaysin and KL001) abolished the stimulating effects on extracellular matrix production and chondrogenic marker gene expression.

CONCLUSIONS

Results from this study suggest that a functional molecular clockwork plays a positive role in tissue homeostasis and histogenesis during early chondrogenesis.

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.

SHP2 inhibition attenuates osteoarthritis by maintaining homeostasis of cartilage metabolism via the DOK1/UPP1/uridine cascade

OBJECTIVE

Protein tyrosine kinases (PTKs) regulate osteoarthritis (OA) progression by activating a series of signal transduction pathways. However, the roles of protein tyrosine phosphatases (PTPs) in OA remain obscure.

METHODS

The expression of 107 PTP genes in human OA cartilage was analyzed based on a single-cell sequencing dataset. The enzyme activity of the PTP SHP2 was detected in primary chondrocytes after interleukin (IL)-1β treatment and in human OA cartilage. Destabilized medial meniscus (DMM) model and IL-1β-stimulated primary mouse chondrocytes were treated with an SHP2 inhibitor and celecoxib (a clinical drug for the treatment of OA). The function of SHP2 in OA pathogenesis was further verified in Aggrecan-Cre^ERT ; SHP2 ^flox/flox mice. The downstream protein expression profile and dephosphorylated substrate of SHP2 were examined by tandem mass tag (TMT) labeling-based global proteomic and stable isotope labeling using amino acids in cell culture (SILAC)-labeled tyrosine phosphoproteomic analysis, respectively.

RESULTS

SHP2 enzyme activity significantly increased in human OA samples with serious articular cartilage injury and in IL-1β-stimulated chondrocytes. Pharmacological inhibition or genetic deletion of SHP2 ameliorated OA progression. SHP2 inhibitors dramatically reduced the expression of cartilage degradation-related genes and simultaneously promoted the expression of cartilage synthesis-related genes. Mechanistically, SHP2 inhibition suppressed the dephosphorylation of DOK1 and subsequently reduced the expression of uridine phosphorylase 1 and increased uridine level, thereby contributing to the homeostasis of cartilage metabolism.

CONCLUSIONS

SHP2 is a novel accelerator of the imbalance in the cartilage homeostasis. Specific inhibition of SHP2 may ameliorate OA by maintaining the anabolic and catabolic balance.

Distinct whole-blood transcriptome profile of children with metabolic healthy overweight/obesity compared to metabolic unhealthy overweight/obesity

BACKGROUND

Youth populations with overweight/obesity (OW/OB) exhibit heterogeneity in cardiometabolic health phenotypes. The underlying mechanisms for those differences are still unclear. This study aimed to analyze the whole-blood transcriptome profile (RNA-seq) of children with metabolic healthy overweight/obesity (MHO) and metabolic unhealthy overweight/obesity (MUO) phenotypes.

METHODS

Twenty-seven children with OW/OB (10.1 ± 1.3 years, 59% boys) from the ActiveBrains project were included. MHO was defined as having none of the following criteria for metabolic syndrome: elevated fasting glucose, high serum triglycerides, low high-density lipoprotein-cholesterol, and high systolic or diastolic blood pressure, while MUO was defined as presenting one or more of these criteria. Inflammatory markers were additionally determined. Total blood RNA was analyzed by 5'-end RNA-sequencing.

RESULTS

Whole-blood transcriptome analysis revealed a distinct pattern of gene expression in children with MHO compared to MUO children. Thirty-two genes differentially expressed were linked to metabolism, mitochondrial, and immune functions.

CONCLUSIONS

The identified gene expression patterns related to metabolism, mitochondrial, and immune functions contribute to a better understanding of why a subset of the population remains metabolically healthy despite having overweight/obesity.

IMPACT

A distinct pattern of whole-blood transcriptome profile (RNA-seq) was identified in children with metabolic healthy overweight/obesity (MHO) compared to metabolic unhealthy overweight/obesity (MUO) phenotype. The most relevant genes in understanding the molecular basis underlying the MHO/MUO phenotypes in children could be: RREB1, FAM83E, SLC44A1, NRG1, TMC5, CYP3A5, TRIM11, and ADAMTSL2. The identified whole-blood transcriptome profile related to metabolism, mitochondrial, and immune functions contribute to a better understanding of why a subset of the population remains metabolically healthy despite having overweight/obesity.

OATargets: a knowledge base of genes associated with osteoarthritis joint damage in animals

OBJECTIVES

To collate the genes experimentally modulated in animal models of osteoarthritis (OA) and compare these data with OA transcriptomics data to identify potential therapeutic targets.

METHODS

PubMed searches were conducted to identify publications describing gene modulations in animal models. Analysed gene expression data were retrieved from the SkeletalVis database of analysed skeletal microarray and RNA-Seq expression data. A network diffusion approach was used to predict new genes associated with OA joint damage.

RESULTS

A total of 459 genes were identified as having been modulated in animal models of OA, with ageing and post-traumatic (surgical) models the most prominent. Ninety-eight of the 143 genes (69%) genetically modulated more than once had a consistent effect on OA joint damage severity. Several discrepancies between different studies were identified, providing lessons on interpretation of these data. We used the data collected along with OA gene expression data to expand existing annotations and prioritise the most promising therapeutic targets, which we validated using the latest reported associations. We constructed an online database OATargets to allow researchers to explore the collated data and integrate it with existing OA and skeletal gene expression data.

CONCLUSIONS

We present a comprehensive survey and online resource for understanding gene regulation of animal model OA pathogenesis.

Endoglin Level Is Critical for Cartilage Tissue Formation In Vitro by Passaged Human Chondrocytes

Transforming growth factor beta (TGFβ) signaling is required for in vitro chondrogenesis. In animal models of osteoarthritis (OA), TGFβ receptor alterations are detected in chondrocytes in severe OA cartilage. It is not known whether such changes are dependent on the grade of human OA and if they affect chondrogenesis. Thus, the purpose of this study was to determine if human OA chondrocytes obtained from low-grade or high-grade disease could form cartilage tissue and to assess the role of the co-receptors, endoglin (ENG) and TGFβ receptor 3 (TGFBRIII), in the regulation of this tissue generation in vitro. We hypothesized that the grade of OA disease would not affect the ability of cells to form cartilage tissue and that the TGFβ co-receptor, ENG, would be critical to regulating tissue formation. Chondrocytes isolated from low-grade OA or high-grade OA human articular cartilage (AC) were analyzed directly (P0) or passaged in monolayer to P2. Expression of the primary TGFβ receptor ALK5, and the co-receptors ENG and TGFβRIII, was assessed by image flow cytometry. To assess the ability to form cartilaginous tissue, cells were placed in three-dimensional culture at high density and cultured in chondrogenic media containing TGFβ3. ENG knockdown was used to determine its role in regulating tissue formation. Overall, grade-specific differences in expression of ALK5, ENG, and TGFβRIII in primary or passaged chondrocytes were not detected; however, ENG expression increased significantly after passaging. Despite the presence of ALK5, P0 cells did not form cartilaginous tissue. In contrast, P2 cells derived from low-grade and high-grade OA AC formed hyaline-like cartilaginous tissues of similar quality. Knockdown of ENG in P2 cells inhibited cartilaginous tissue formation compared to controls indicating that the level of ENG protein expression is critical for in vitro chondrogenesis by passaged articular chondrocytes. This study demonstrates that it is not the grade of OA, but the levels of ENG in the presence of ALK5 that influences the ability of human passaged articular chondrocytes to form cartilaginous tissue in vitro in 3D culture. This has implications for cartilage repair therapies. Impact statement These findings are important clinically, given the limited availability of osteoarthritis (OA) cartilage tissue. Being able to use cells from all grades of OA will increase our ability to obtain sufficient cells for cartilage repair. In addition, it is possible that endoglin (ENG) levels, in the presence of ALK5 expression, may be suitable to use as biomarkers to identify cells able to produce cartilage.

Rheometer enabled study of cartilage frequency-dependent properties

Despite the well-established dependence of cartilage mechanical properties on the frequency of the applied load, most research in the field is carried out in either load-free or constant load conditions because of the complexity of the equipment required for the determination of time-dependent properties. These simpler analyses provide a limited representation of cartilage properties thus greatly reducing the impact of the information gathered hindering the understanding of the mechanisms involved in this tissue replacement, development and pathology. More complex techniques could represent better investigative methods, but their uptake in cartilage research is limited by the highly specialised training required and cost of the equipment. There is, therefore, a clear need for alternative experimental approaches to cartilage testing to be deployed in research and clinical settings using more user-friendly and financial accessible devices. Frequency dependent material properties can be determined through rheometry that is an easy to use requiring a relatively inexpensive device; we present how a commercial rheometer can be adapted to determine the viscoelastic properties of articular cartilage. Frequency-sweep tests were run at various applied normal loads on immature, mature and trypsinased (as model of osteoarthritis) cartilage samples to determine the dynamic shear moduli (G*, G′ G″) of the tissues. Moduli increased with increasing frequency and applied load; mature cartilage had generally the highest moduli and GAG depleted samples the lowest. Hydraulic permeability (K_H) was estimated from the rheological data and decreased with applied load; GAG depleted cartilage exhibited higher hydraulic permeability than either immature or mature tissues. The rheometer-based methodology developed was validated by the close comparison of the rheometer-obtained cartilage characteristics (G*, G′, G″, K_H) with results obtained with more complex testing techniques available in literature. Rheometry is relatively simpler and does not require highly capital intensive machinery and staff training is more accessible; thus the use of a rheometer would represent a cost-effective approach for the determination of frequency-dependent properties of cartilage for more comprehensive and impactful results for both healthcare professional and R&D.

Widespread regulation of gene expression by glucocorticoids in chondrocytes from patients with osteoarthritis as determined by RNA-Seq

Background

Intra-articular glucocorticoid (GC) injections are widely used as a symptomatic treatment for osteoarthritis (OA). However, there are also concerns about their potentially harmful effects, and their detailed effects on chondrocyte phenotype remain poorly understood.

Methods

We studied the effects of dexamethasone on gene expression in OA chondrocytes with RNA-Seq. Chondrocytes were isolated from the cartilage from OA patients undergoing knee replacement surgery and cultured with or without dexamethasone for 24 h. Total RNA was isolated and sequenced, and functional analysis was performed against the Gene Ontology (GO) database. Results for selected genes were confirmed with RT-PCR. We also investigated genes linked to OA in recent genome-wide expression analysis (GWEA) studies.

Results

Dexamethasone increased the expression of 480 and reduced that of 755 genes with a fold change (FC) 2.0 or greater. Several genes associated with inflammation and cartilage anabolism/catabolism as well as lipid and carbohydrate metabolism were among the most strongly affected genes. In the GO analysis, genes involved in the extracellular matrix organization, cell proliferation and adhesion, inflammation, and collagen synthesis were enriched among the significantly affected genes. In network analysis, NGF, PI3KR1, and VCAM1 were identified as central genes among those most strongly affected by dexamethasone.

Conclusions

This is the first study investigating the genome-wide effects of GCs on the gene expression in OA chondrocytes. In addition to clear anti-inflammatory and anticatabolic effects, GCs affect lipid and glucose metabolism in chondrocytes, an observation that might be particularly important in the metabolic phenotype of OA.

A Microarray Study of Articular Cartilage in Relation to Obesity and Severity of Knee Osteoarthritis

OBJECTIVE

To query the transcript-level changes in the medial and lateral tibial plateau cartilage in tandem with obesity in patients with end-stage osteoarthritis (OA).

DESIGN

Cartilage was obtained from 23 patients (20 obese [body mass index > 30 kg/m²], 3 overweight [body mass index < 30 kg/m²]) at the time of total knee replacement. Cartilage integrity was assessed using Outerbridge scale, while radiographic changes were scored on preoperative X-rays using Kellgren-Lawrence (K-L) classification. RNA was probed for differentially expressed transcripts between medial and lateral compartments using Affymetrix Gene 2.0 ST Array and validated via real-time polymerase chain reaction. Gene ontology and pathway analyses were also queried.

RESULTS

Scoring of cartilage integrity by the Outerbridge scale indicated that the medial and lateral compartments were similar, while scoring by the K-L classification indicated that the medial compartment was more severely damaged than the lateral compartment. We observed a distinct transcript profile with >50% of transcripts unique between medial and lateral compartments. MMP13 and COL2A1 were more highly expressed in medial versus lateral compartment. Polymerase chain reaction confirmed expression of 4 differentially expressed transcripts. Numerous transcripts, biological processes, and pathways were significantly different between overweight and obese patients with a differential response of obesity on medial and lateral compartments.

CONCLUSIONS

Our findings support molecular differences between medial and lateral compartments reflective of the greater severity of OA in the medial compartment. The K-L system better reflected the molecular results than did the Outerbridge. Moreover, the molecular effect of obesity was different between the medial and lateral compartments of the same knee plausibly reflecting the molecular effects of differential biomechanical loading.

Gene Expression Profiling Studies Using Microarray in Osteoarthritis: Genes in Common and Different Conditions

Osteoarthritis (OA), which is characterized mainly by cartilage degradation, is the most prevalent joint disorder worldwide. Although OA is identified as a major cause of joint pain, disability, and socioeconomic burden, the etiology of OA is still not clearly known. Recently, gene microarray analysis has become an efficient method for the research of complex diseases and has been employed to determine what genes and pathways are involved in the pathological process of OA. In this review, OA study results over the last decade are summarized for gene expression profiling of various tissues, such as cartilage, subchondral bone, and synovium in human OA and mouse OA models. Many differentially expressed genes, which mainly involve matrix metabolism, bone turnover, and inflammation pathways, were identified in diseased compared with "normal" tissues. Nevertheless, rare common genes were reported from studies using different tissue sources, microarray chips, and research designs. Thus, future novel and carefully designed microarray studies are required to elucidate underlying genetic mechanisms in the pathogenesis of OA as well as new directions for potential OA-targeted pharmaceutical therapies.

Secreted Frizzled-related proteins (sFRPs) in osteo-articular diseases: much more than simple antagonists of Wnt signaling?

Osteo-articular diseases are characterized by a dysregulation of joint and/or bone homeostasis. These include diseases affecting the joints originally, such as osteoarthritis and rheumatoid arthritis, or the bone, such as osteoporosis. Inflammation and the involvement of Wingless-related integration site (Wnt) signaling pathways are key pathophysiological features of these diseases resulting in tissue degradation by matrix-degrading enzymes, namely matrix metalloproteinases (MMPs) and a disintegrin and metalloproteinases with thrombospondin motifs (ADAMTs), secreted by the joint resident cells and/or by infiltrating immune cells. Activation of Wnt signaling pathways is modulated by different families of proteins, including Dickkopfs and the secreted Frizzled-related proteins (sFRPs). The sFRP family is composed of five secreted glycoproteins in mammals that regulate Wnt signaling in the extracellular compartment. Indeed, sFRPs are able to bind both to the soluble Wnt ligands and to their cell membrane receptors, the Frizzled proteins. Their expression profile is altered in osteo-articular diseases, suggesting that they could account for the abnormal activation of Wnt pathways. In the present article, we review how sFRPs are more than simple antagonists of the Wnt signaling pathways and discuss their pathophysiological relevance in the context of osteo-articular diseases. We detail their Wnt-dependent and their Wnt-independent roles, with a particular emphasis on their ability to modulate the inflammatory response and extracellular matrix (ECM) remodeling. We also discuss their potential therapeutic use with a focus on bone remodeling, osteo-articular cancers, and tissue engineering.

Neural EGFL like 1 as a potential pro-chondrogenic, anti-inflammatory dual-functional disease-modifying osteoarthritis drug

Arthritis, an inflammatory condition that causes pain and cartilage destruction in joints, affects over 54.4 million people in the US alone. Here, for the first time, we demonstrated the emerging role of neural EGFL like 1 (NELL-1) in arthritis pathogenesis by showing that Nell-1-haploinsufficient (Nell-1^+/6R) mice had accelerated and aggravated osteoarthritis (OA) progression with elevated inflammatory markers in both spontaneous primary OA and chemical-induced secondary OA models. In the chemical-induced OA model, intra-articular injection of interleukin (IL)1β induced more severe inflammation and cartilage degradation in the knee joints of Nell-1^+/6R mice than in wildtype animals. Mechanistically, in addition to its pro-chondrogenic potency, NELL-1 also effectively suppressed the expression of inflammatory cytokines and their downstream cartilage catabolic enzymes by upregulating runt-related transcription factor (RUNX)1 in mouse and human articular cartilage chondrocytes. Notably, NELL-1 significantly reduced IL1β-stimulated inflammation and damage to articular cartilage in vivo. In particular, NELL-1 administration markedly reduced the symptoms of antalgic gait observed in IL1β-challenged Nell-1^+/6R mice. Therefore, NELL-1 is a promising pro-chondrogenic, anti-inflammatory dual-functional disease-modifying osteoarthritis drug (DMOAD) candidate for preventing and suppressing arthritis-related cartilage damage.

lnc-SAMD14-4 can regulate expression of the COL1A1 and COL1A2 in human chondrocytes

Osteoarthritis (OA) is the most common motor system disease in aging people, characterized by matrix degradation, chondrocyte death, and osteophyte formation. OA etiology is unclear, but long noncoding RNAs (lncRNAs) that participate in numerous pathological and physiological processes may be key regulators in the onset and development of OA. Because profiling of lncRNAs and their biological function in OA is not understood, we measured lncRNA and mRNA expression profiles using high-throughput microarray to study human knee OA. We identified 2,042 lncRNAs and 2,011 mRNAs that were significantly differentially expressed in OA compared to non-OA tissue (>2.0- or < - 2.0-fold change; p < 0.5), including 1,137 lncRNAs that were upregulated and 905 lncRNAs that were downregulated. Also, 1,386 mRNA were upregulated and 625 mRNAs were downregulated. QPCR was used to validate chip results. Gene Ontology analysis and the Kyoto Encyclopedia of Genes and Genomes was used to study the biological function enrichment of differentially expressed mRNA. Additionally, coding-non-coding gene co-expression (CNC) network construction was performed to explore the relevance of dysregulated lncRNAs and mRNAs. Finally, the gain/loss of function experiments of lnc-SAMD14-4 was implemented in IL-1β-treated human chondrocytes. In general, this study provides a preliminary database for further exploring lncRNA-related mechnisms in OA.

Whole-transcriptome sequencing of knee joint cartilage from osteoarthritis patients

Objectives

The aim of this study was to provide a comprehensive understanding of alterations in messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs) in cartilage affected by osteoarthritis (OA).

Methods

The expression profiles of mRNAs, lncRNAs, and circRNAs in OA cartilage were assessed using whole-transcriptome sequencing. Bioinformatics analyses included prediction and reannotation of novel lncRNAs and circRNAs, their classification, and their placement into subgroups. Gene ontology and pathway analysis were performed to identify differentially expressed genes (DEGs), differentially expressed lncRNAs (DELs), and differentially expressed circRNAs (DECs). We focused on the overlap of DEGs and targets of DELs previously identified in seven high-throughput studies. The top ten DELs were verified by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) in articular chondrocytes, both in vitro and in vivo.

Results

In total, 739 mRNAs, 1152 lncRNAs, and 42 circRNAs were found to be differentially expressed in OA cartilage tissue. Among these, we identified 18 overlapping DEGs and targets of DELs, and the top ten DELs were screened by expression profile analysis as candidate OA-related genes. WISP2, ATF3, and CHI3L1 were significantly increased in both normal versus OA tissues and normal versus interleukin (IL)-1β-induced OA-like cell models, while ADAM12, PRELP, and ASPN were shown to be significantly decreased. Among the identified DELs, we observed higher expression of ENST00000453554 and MSTRG.99593.3, and lower expression of MSTRG.44186.2 and NONHSAT186094.1 in normal versus OA cells and tissues.

Conclusion

This study revealed expression patterns of coding and noncoding RNAs in OA cartilage, which added sets of genes and noncoding RNAs to the list of candidate diagnostic biomarkers and therapeutic agents for OA patients.

Cite this article: H. Li, H. H. Yang, Z. G. Sun, H. B. Tang, J. K. Min. Whole-transcriptome sequencing of knee joint cartilage from osteoarthritis patients. Bone Joint Res 2019;8:290–303. DOI: 10.1302/2046-3758.87.BJR-2018-0297.R1.

Hydrogen and Water Bonding between Glycosaminoglycans and Phospholipids in the Synovial Fluid: Molecular Dynamics Study

Synovial fluid is a lubricant of the synovial joint that shows remarkable tribological properties. These properties originate in the synergy between its components, with two of its major components, glycosaminoglycans (GAGs) and phospholipids (PLs), playing a major role in boundary and mixed lubrication regimes. All-atom molecular dynamic simulations were performed to investigate the way these components bond. Hyaluronic acid (HA) and chondroitin sulphate (CS) bonding with three types of lipids was tested. The results show that both glycosaminoglycans bind lipids at a similar rate, except for 1,2-d-ipalmitoyl-sn-glycero-3-phosphoethanolamine lipids, which bind to chondroitin at a much higher rate than to hyaluronan. The results suggest that different synovial fluid lipids may play a different role when binding to both hyaluronan and chondroitin sulphate. The presented results may help in understanding a process of lubrication of articular cartilage at a nanoscale level.

Distinct degenerative phenotype of articular cartilage from knees with meniscus tear compared to knees with osteoarthritis

OBJECTIVE

To compare the transcriptome of articular cartilage from knees with meniscus tears to knees with end-stage osteoarthritis (OA).

DESIGN

Articular cartilage was collected from the non-weight bearing medial intercondylar notch of knees undergoing arthroscopic partial meniscectomy (APM; N = 10, 49.7 ± 10.8 years, 50% females) for isolated medial meniscus tears and knees undergoing total knee arthroplasty (TKA; N = 10, 66.0 ± 7.6 years, 70% females) due to end-stage OA. Ribonucleic acid (RNA) preparation was subjected to SurePrint G3 human 8 × 60K RNA microarrays to probe differentially expressed transcripts followed by computational exploration of underlying biological processes. Real-time polymerase chain reaction amplification was performed on selected transcripts to validate microarray data.

RESULTS

We observed that 81 transcripts were significantly differentially expressed (45 elevated, 36 repressed) between APM and TKA samples (≥ 2 fold) at a false discovery rate of ≤ 0.05. Among these, CFD, CSN1S1, TSPAN11, CSF1R and CD14 were elevated in the TKA group, while CHI3L2, HILPDA, COL3A1, COL27A1 and FGF2 were highly expressed in APM group. A few long intergenic non-coding RNAs (lincRNAs), small nuclear RNAs (snoRNAs) and antisense RNAs were also differentially expressed between the two groups. Transcripts up-regulated in TKA cartilage were enriched for protein localization and activation, chemical stimulus, immune response, and toll-like receptor signaling pathway. Transcripts up-regulated in APM cartilage were enriched for mesenchymal cell apoptosis, epithelial morphogenesis, canonical glycolysis, extracellular matrix organization, cartilage development, and glucose catabolic process.

CONCLUSIONS

This study suggests that APM and TKA cartilage express distinct sets of OA transcripts. The gene profile in cartilage from TKA knees represents an end-stage OA whereas in APM knees it is clearly earlier in the degenerative process.

Altered expression of the core circadian clock component PERIOD2 contributes to osteoarthritis-like changes in chondrocyte activity

In osteoarthritis, chondrocytes undergo a phenotype shift characterised by reduced expression of SOX9 (sry-box 9) and increased production of cartilage-degrading enzymes, e.g. MMP13 (matrix metalloproteinase 13) and ADAMTS5 (a disintegrin and metalloproteinase with thrombospondin motifs 5). The chondrocyte clock is also altered. Specifically, the peak level of PER2 is elevated, but peak level of BMAL1 reduced in osteoarthritic chondrocytes. The purpose of this study was to determine whether increased PER2 expression causes disease-associated changes in chondrocyte activity and to identify whether known risk factors for osteoarthritis induce changes in PER2 and BMAL1 expression. Primary human chondrocytes isolated from macroscopically normal cartilage were serum-starved overnight then re-fed with serum-replete media with/without interleukin 1β (IL-1β) (10 ng/mL), hydrogen peroxide (100 µM) or basic calcium phosphate (BCP) crystals (50 µg/mL). Peak level of BMAL1 was lower, whereas PER2 levels remained elevated for longer, in chondrocytes treated with IL-1β, hydrogen peroxide or BCP crystals compared to untreated cells. Levels of SOX9 were lower, whereas levels of ADAMTS5 and MMP13 were higher, in chondrocytes exposed to any of the three treatments compared to untreated cells. Knockdown of PER2 using siRNA partially abrogated the effects of each treatment on chondrocyte phenotype marker expression. Similarly, in chondrocytes isolated from osteoarthritic cartilage PER2 knockdown was associated with increased SOX9, reduced ADAMTS5 and reduced RNA and protein levels of MMP13 indicating partial mitigation of the osteoarthritic phenotype. Conversely, further ablation of BMAL1 expression in osteoarthritic chondrocytes resulted in a further reduction in SOX9 and increase in MMP13 expression. Overexpression of PER2 in the H5 chondrocyte cell line led to increased ADAMTS5 and MMP13 and decreased SOX9 expression. Localised inflammation, oxidative stress and BCP crystal deposition in osteoarthritic joints may contribute to disease pathology by inducing changes in the chondrocyte circadian clock.

Identification of TGFβ-related genes regulated in murine osteoarthritis and chondrocyte hypertrophy by comparison of multiple microarray datasets

OBJECTIVE

Osteoarthritis (OA) is a joint disease characterized by progressive degeneration of articular cartilage. Some features of OA, including chondrocyte hypertrophy and focal calcification of articular cartilage, resemble the endochondral ossification processes. Alterations in transforming growth factor β (TGFβ) signaling have been associated with OA as well as with chondrocyte hypertrophy. Our aim was to identify novel candidate genes implicated in chondrocyte hypertrophy during OA pathogenesis by determining which TGFβ-related genes are regulated during murine OA and endochondral ossification.

METHODS

A list of 580 TGFβ-related genes, including TGFβ signaling pathway components and TGFβ-target genes, was generated. Regulation of these TGFβ-related genes was assessed in a microarray of murine OA cartilage: 1, 2 and 6 weeks after destabilization of the medial meniscus (DMM). Subsequently, genes regulated in the DMM model were studied in two independent murine microarray datasets on endochondral ossification: the growth plate and transient embryonic cartilage (joint development).

RESULTS

A total of 106 TGFβ-related genes were differentially expressed in articular cartilage of DMM-operated mice compared to sham-control. From these genes, 43 were similarly regulated during chondrocyte hypertrophy in the growth plate or embryonic joint development. Among these 43 genes, 18 genes have already been associated with OA. The remaining 25 genes were considered as novel candidate genes involved in OA pathogenesis and endochondral ossification. In supplementary data of published human OA microarrays we found indications that 15 of the 25 novel genes are indeed regulated in articular cartilage of human OA patients.

CONCLUSION

By focusing on TGFβ-related genes during OA and chondrocyte hypertrophy in mice, we identified 18 known and 25 new candidate genes potentially implicated in phenotypical changes in chondrocytes leading to OA. We propose that 15 of these candidates warrant further investigation as therapeutic target for OA as they are also regulated in articular cartilage of OA patients.

Identification of transcription factors responsible for dysregulated networks in human osteoarthritis cartilage by global gene expression analysis

OBJECTIVE

Osteoarthritis (OA) is the most prevalent joint disease. As disease-modifying therapies are not available, novel therapeutic targets need to be discovered and prioritized for their importance in mediating the abnormal phenotype of cells in OA-affected joints. Here, we generated a genome-wide molecular profile of OA to elucidate regulatory mechanisms of OA pathogenesis and to identify possible therapeutic targets using integrative analysis of mRNA-sequencing data obtained from human knee cartilage.

DESIGN

RNA-sequencing (RNA-seq) was performed on 18 normal and 20 OA human knee cartilage tissues. RNA-seq datasets were analysed to identify genes, pathways and regulatory networks that were dysregulated in OA.

RESULTS

RNA-seq data analysis revealed 1332 differentially expressed (DE) genes between OA and non-OA samples, including known and novel transcription factors (TFs). Pathway analysis identified 15 significantly perturbed pathways in OA with ECM-related, PI3K-Akt, HIF-1, FoxO and circadian rhythm pathways being the most significantly dysregulated. We selected DE TFs that are enriched for regulating DE genes in OA and prioritized these TFs by creating a cartilage-specific interaction subnetwork. This analysis revealed eight TFs, including JUN, Early growth response (EGR)1, JUND, FOSL2, MYC, KLF4, RELA, and FOS that both target large numbers of dysregulated genes in OA and are themselves suppressed in OA.

CONCLUSIONS

We identified a novel subnetwork of dysregulated TFs that represent new mediators of abnormal gene expression and promising therapeutic targets in OA.

Prenatal caffeine exprosure increases adult female offspring rat's susceptibility to osteoarthritis via low-functional programming of cartilage IGF-1 with histone acetylation

Our previous in vivo studies showed that prenatal caffeine exposure (PCE) could restrain the development of chondrogenesis, which may delay fetal articular cartilage development and increase susceptibility to osteoarthritis in adults. So, the goal of the current study is to clarify theincreasing susceptibility to adult osteoarthritis in caffeine-exposed female offspring and its'mechanism. Pregnant rats were treated with 120 mg/kg·d caffeine or equal volumes of saline from gestational day (GD) 9 to 20. knee joints were collected from GD20 female fetuses and 18-week old female offspring which was treated with strenuous running for 6 weeks (55 min/day at 20 m/min) load to induce osteoarthritis. Knee joints from GD20 fetuses and adult offspring were collected for histochemistry and immunohistochemistry. Next, chondrocytes were isolated from 1-day-old newborn rats and in vitro studies were conducted where the cells in primary culture were exposed to 1, 10, and 100 μM caffeine and 250, 500, and 1,250 nM corticosterone. Insulin-like growth factor 1 (IGF-1) signal pathway genes' expression levels in fetal chondrocytes were studied, and IGF-1 histone acetylation was detected in vitro. Immunohistochemical results showed low expression levels of IGF-1 signaling genes (IGF-1, IRS-1, AKT, and COL2A1) both in fetal and adult cartilage with PCE. For adult offspring, histological results and Mankin score revealed increased cartilage destruction and accelerated osteoarthritis progression in PCE group with strenuous running exercise. Analysis in vitro revealed that caffeine and corticosterone impeded the expression of IGF-1 signaling pathway aggrecan and COL2A1 genes, but only corticosterone decreased H3K9 and H3K27 acetylation in the IGF-1 promoter region. In concluson, PCE low functional programmed cartilage IGF-1 by histone acetylation modification via overexposure to corticosterone and delayed articular cartilage development from fetus to adults. Then, the delayed cartilage development increased susceptibility to osteoarthritis in offsprings.

A whole-genome transcriptome analysis of articular chondrocytes in secondary osteoarthritis of the hip

Objective

To date, exhaustive gene expression analyses of chondrocytes in hip osteoarthritis (OA) have yielded specific gene expression patterns. No study has reported on the exhaustive transcriptome of secondary hip OA based on acetabular dysplasia in a Japanese population, while previous reports have focused on primary or idiopathic hip OA in Caucasian populations. This study aims to search for specific gene expression patterns of secondary hip OA chondrocytes by transcriptome analysis.

Design

Human articular cartilage was obtained from femoral heads following hemiarthroplasty for femoral neck fracture (N = 8; non-OA) and total hip arthroplasty for secondary hip OA (N = 12). Total RNA was extracted from the articular cartilage and submitted for microarray analysis. The obtained data were used to perform gene expression analysis, GO enrichment analysis and pathway analysis and were compared with data from primary hip OA in Caucasian populations in the literature.

Results

We identified 888 upregulated (fold change: FC ≥ 2) and 732 downregulated (FC ≤ 0.5) genes in hip OA versus non-OA chondrocytes, respectively. Only 10% of upregulated genes were common between the secondary and primary OA. The newly found genes prominently overexpressed in the secondary hip OA chondrocytes were DPT, IGFBP7, and KLF2. Pathway analysis revealed extracellular matrix (ECM)-receptor interaction as an OA-related pathway, which was similar to previous reports in primary hip OA.

Conclusions

This is the first study to report the genome-wide transcriptome of secondary hip OA chondrocytes and demonstrates new potential OA-related genes. Gene expression patterns were different between secondary and primary hip OA, although the results of pathway and functional analysis were similar.

Quality Assessment of Surgical Disc Samples Discriminates Human Annulus Fibrosus and Nucleus Pulposus on Tissue and Molecular Level

A discrimination of the highly specialised annulus fibrosus (AF) and nucleus pulposus (NP) cells in the mature human intervertebral disc (IVD) is thus far still not possible in a reliable way. The aim of this study was to identify molecular markers that distinguish AF and NP cells in human disc tissue using microarray analysis as a screening tool. AF and NP samples were obtained from 28 cervical discs. First, all samples underwent quality sorting using two novel scoring systems for small-sized disc tissue samples including macroscopic, haptic and histological evaluation. Subsequently, samples with clear disc characteristics of either AF or NP that were free from impurities of foreign tissue (IVD score) and with low signs of disc degeneration on cellular level (DD score) were selected for GeneChip analysis (HGU1332P). The 11 AF and 9 NP samples showed distinctly different genome-wide transcriptomes. The majority of differentially expressed genes (DEGs) could be specifically assigned to the AF, whereas no DEG was exclusively expressed in the NP. Nevertheless, we identified 11 novel marker genes that clearly distinguished AF and NP, as confirmed by quantitative gene expression analysis. The novel established scoring systems and molecular markers showed the identity of AF and NP in disc starting material and are thus of great importance in the quality assurance of cell-based therapeutics in regenerative treatment of disc degeneration.

Use of integrative epigenetic and mRNA expression analyses to identify significantly changed genes and functional pathways in osteoarthritic cartilage

Aim

Osteoarthritis (OA) is caused by complex interactions between genetic and environmental factors. Epigenetic mechanisms control the expression of genes and are likely to regulate the OA transcriptome. We performed integrative genomic analyses to define methylation-gene expression relationships in osteoarthritic cartilage.

Patients and Methods

Genome-wide DNA methylation profiling of articular cartilage from five patients with OA of the knee and five healthy controls was conducted using the Illumina Infinium HumanMethylation450 BeadChip (Illumina, San Diego, California). Other independent genome-wide mRNA expression profiles of articular cartilage from three patients with OA and three healthy controls were obtained from the Gene Expression Omnibus (GEO) database. Integrative pathway enrichment analysis of DNA methylation and mRNA expression profiles was performed using integrated analysis of cross-platform microarray and pathway software. Gene ontology (GO) analysis was conducted using the Database for Annotation, Visualization and Integrated Discovery (DAVID).

Results

We identified 1265 differentially methylated genes, of which 145 are associated with significant changes in gene expression, such as DLX5, NCOR2 and AXIN2 (all p-values of both DNA methylation and mRNA expression < 0.05). Pathway enrichment analysis identified 26 OA-associated pathways, such as mitogen-activated protein kinase (MAPK) signalling pathway (p = 6.25 × 10-4), phosphatidylinositol (PI) signalling system (p = 4.38 × 10-3), hypoxia-inducible factor 1 (HIF-1) signalling pathway (p = 8.63 × 10-3 pantothenate and coenzyme A (CoA) biosynthesis (p = 0.017), ErbB signalling pathway (p = 0.024), inositol phosphate (IP) metabolism (p = 0.025), and calcium signalling pathway (p = 0.032).

Conclusion

We identified a group of genes and biological pathwayswhich were significantly different in both DNA methylation and mRNA expression profiles between patients with OA and controls. These results may provide new clues for clarifying the mechanisms involved in the development of OA.

Cite this article: A. He, Y. Ning, Y. Wen, Y. Cai, K. Xu, Y. Cai, J. Han, L. Liu, Y. Du, X. Liang, P. Li, Q. Fan, J. Hao, X. Wang, X. Guo, T. Ma, F. Zhang. Use of integrative epigenetic and mRNA expression analyses to identify significantly changed genes and functional pathways in osteoarthritic cartilage. Bone Joint Res 2018;7:343–350. DOI: 10.1302/2046-3758.75.BJR-2017-0284.R1.

Hyaluronan-Chondroitin Sulfate Anomalous Crosslinking Due to Temperature Changes

Glycosaminoglycans are a wide class of biopolymers showing great lubricating properties due to their structure and high affinity to water. Two of them, hyaluronic acid and chondroitin sulfate, play an important role in articular cartilage lubrication. In this work, we present results of the all-atom molecular dynamics simulations of both molecules placed in water-based solution. To mimic changes of the physiological conditions, especially temperature, of the synovial fluid in joints under successive load (e.g., walking, jogging, jumping), simulations have been performed at different physiological temperatures in the range of 300 to 320 Kelvin (normal intra-articular temperature is 305 K). The stability of the biopolymeric network at equilibrium (isothermal and isobaric) conditions has been studied. To understand the process of physical crosslinking, the dynamics of intra- and intermolecular hydrogen bonds forming and breaking have been studied. The results show that following addition of chondroitin sulfate, hyaluronan creates more intermolecular hydrogen bonds than when in homogeneous solution. The presence of chondroitin in a hyaluronan network is beneficial as it may increase its stability. Presented data show hyaluronic acid and chondroitin sulfate as viscosity modifiers related to their crosslinking properties in different physicochemical conditions.

The Anomalies of Hyaluronan Structures in Presence of Surface Active Phospholipids-Molecular Mass Dependence

Interactions between hyaluronan (A-) and phospholipids play a key role in many systems in the human body. One example is the articular cartilage system, where the synergistic effect of such interactions supports nanoscale lubrication. A molecular dynamics simulation has been performed to understand the process of formation of hydrogen bonds inside the hyaluronan network, both in the presence and absence of phospholipids. Additionally, the effect of the molecular mass of (A-) was analyzed. The main finding of this work is a robust demonstration of the optimal parameters (H-bond energy, molecular mass) influencing the facilitated lubrication mechanism of the articular cartilage system. Simulation results show that the presence of phospholipids has the greatest influence on hyaluronan at low molecular mass. We also show the specific sites of H-bonding between chains. Simulation results can help to understand how hyaluronan and phospholipids interact at several levels of articular cartilage system functioning.

Transcriptome comparison of meniscus from patients with and without osteoarthritis

OBJECTIVE

To assess the impact of osteoarthritis (OA) on the meniscus by comparing transcripts and biological processes in the meniscus between patients with and without OA.

DESIGN

RNA microarrays were used to identify transcripts differentially expressed (DE) in meniscus obtained from 12 OA and 12 non-OA patients. The non-OA specimens were obtained at the time of arthroscopic partial meniscectomy. Real-time PCR was performed on selected transcripts. Biological processes and gene-networking was examined computationally. Transcriptome signatures were mapped with 37 OA-related transcripts to evaluate how meniscus gene expression relates to that of OA cartilage.

RESULTS

We identified 168 transcripts significantly DE between OA (75 elevated, 93 repressed) and non-OA samples (≥1.5-fold). Among these, CSN1S1, COL10A1, WIF1, and SPARCL1 were the most prominent transcripts elevated in OA meniscus, POSTN and VEGFA were most highly repressed in OA meniscus. Transcripts elevated in OA meniscus represented response to external stimuli, cell migration and cell localization while those repressed in OA meniscus represented histone deacetylase activity (related to epigenetics) and skeletal development. Numerous long non-coding RNAs (lncRNAs) were DE between the two groups. When segregated by OA-related transcripts, two distinct clustering patterns appeared: OA meniscus appeared to be more inflammatory while non-OA meniscus exhibited a "repair" phenotype.

CONCLUSIONS

Numerous transcripts with potential relevance to the pathogenesis of OA are DE in OA and non-OA meniscus. These data suggest an involvement of epigenetically regulated histone deacetylation in meniscus tears as well as expression of lncRNAs. Patient clustering based on transcripts related to OA in articular cartilage confirmed distinct phenotypes between injured (non-OA) and OA meniscus.