Advertisement
Research Article|Articles in Press

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

  • A.E.M. Jørgensen
    Correspondence
    Address correspondence and reprint requests to: Institute of Sports Medicine Copenhagen, Department of Orthopaedic Surgery M81, Copenhagen University Hospital of Bispebjerg and Frederiksberg, Copenhagen, Denmark.
    Affiliations
    Institute of Sports Medicine Copenhagen, Department of Orthopaedic Surgery M81, Copenhagen University Hospital of Bispebjerg and Frederiksberg, Copenhagen, Denmark

    Center for Healthy Aging, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
    Search for articles by this author
  • P. Schjerling
    Affiliations
    Institute of Sports Medicine Copenhagen, Department of Orthopaedic Surgery M81, Copenhagen University Hospital of Bispebjerg and Frederiksberg, Copenhagen, Denmark

    Center for Healthy Aging, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
    Search for articles by this author
  • B. DellaValle
    Affiliations
    Department of Endocrinology, Copenhagen University Hospital of Bispebjerg and Frederiksberg, Copenhagen, Denmark

    Copenhagen Center for Translational Research, Copenhagen University Hospital of Bispebjerg and Frederiksberg, Copenhagen, Denmark

    GLX Analytix ApS, Copenhagen, Denmark
    Search for articles by this author
  • J. Rungby
    Affiliations
    Department of Endocrinology, Copenhagen University Hospital of Bispebjerg and Frederiksberg, Copenhagen, Denmark

    Copenhagen Center for Translational Research, Copenhagen University Hospital of Bispebjerg and Frederiksberg, Copenhagen, Denmark
    Search for articles by this author
  • M. Kjær
    Affiliations
    Institute of Sports Medicine Copenhagen, Department of Orthopaedic Surgery M81, Copenhagen University Hospital of Bispebjerg and Frederiksberg, Copenhagen, Denmark

    Center for Healthy Aging, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
    Search for articles by this author
Open AccessPublished:January 28, 2023DOI:https://doi.org/10.1016/j.joca.2023.01.317

      Summary

      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/m2. 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.

      Keywords

      Introduction

      Osteoarthritis (OA) is a painful and debilitating disease of the joint
      • Syx D.
      • Tran P.B.
      • Miller R.E.
      • Malfait A.M.
      Peripheral mechanisms contributing to osteoarthritis pain.
      , and its prevalence increases with age
      • Shane Anderson A.
      • Loeser R.F.
      Why is osteoarthritis an age-related disease?.
      . OA is ultimately treated effectively by joint replacement surgery however, functional outcomes vary, and the lifespan of the prostheses is limited
      • Glyn-Jones S.
      • Palmer A.J.R.
      • Agricola R.
      • Price A.J.
      • Vincent T.L.
      • Weinans H.
      • et al.
      Osteoarthritis.
      . Many non-surgical treatments are recommended, including pain-reducing medication, weight loss, and exercise
      • Kolasinski S.L.
      • Neogi T.
      • Hochberg M.C.
      • Oatis C.
      • Guyatt G.
      • Block J.
      • et al.
      2019 American college of rheumatology/arthritis foundation guideline for the management of osteoarthritis of the hand, hip, and knee.
      , and systematic reviews show that exercise is clinically beneficial
      • Fransen M.
      • McConnell S.
      • Harmer A.R.
      • Van der Esch M.
      • Simic M.
      • Bennell K.L.
      Exercise for osteoarthritis of the knee.
      and safe for the cartilage
      • Bricca A.
      Exercise does not “wear down my knee”: systematic reviews and meta-analyses.
      .
      As sampling from human cartilage is highly invasive and damaging, measurements of cartilage biomarkers have been performed in blood, urine, and synovial fluid (SF) in connection to exercise with months of regular exercise generally reducing markers of cartilage breakdown and inflammation in OA
      • Bricca A.
      • Struglics A.
      • Larsson S.
      • Steultjens M.
      • Juhl C.B.
      • Roos E.M.
      Impact of exercise therapy on molecular biomarkers related to cartilage and inflammation in people at risk of, or with established, knee osteoarthritis: a systematic review and meta-analysis of randomized controlled trials.
      . However, systemic sampling likely does not reflect the cartilage in a single joint. Human cartilage samples can be obtained with joint replacement surgery. This single sampling time point enables a prior stimulus but does not allow for longitudinal observation.
      Studies subjecting human chondrocytes or cartilage explants to mechanical load have shown mechanically regulated genes within groups of matrix proteins, cytokines, growth factors, and degrading enzymes
      • Wong M.
      • Carter D.R.
      Articular cartilage functional histomorphology and mechanobiology: a research perspective.
      . However, conflicting results have been found likely from variations in tissue preparation and loading protocol which further complicates the extrapolation to in vivo conditions to the clinic
      • Wong M.
      • Carter D.R.
      Articular cartilage functional histomorphology and mechanobiology: a research perspective.
      . Thus, despite the comprehensive knowledge, the effect of exercise on human cartilage itself is unknown. Therefore, this study aimed to examine the effect of a single bout of resistance exercise on gene expression and glycosaminoglycan (GAG) content and dynamics in human OA knee cartilage. As the symptomatic benefit of exercise is well described
      • Kraus V.B.
      • Sprow K.
      • Powell K.E.
      • Buchner D.
      • Bloodgood B.
      • Piercy K.
      • et al.
      Effects of physical activity in knee and hip osteoarthritis: a systematic umbrella review.
      , we hypothesise that exercise will provide a beneficial, anti-catabolic and/or anti-inflammatory response at the cellular level.

      Methods

      Trial design

      This was a single day, two-arm, randomised clinical trial of resistance exercise compared to no exercise in patients with OA of the knee scheduled for joint replacement surgery. Primary outcomes were changes in gene expression and GAG content. No study was available for precise sample size calculation, however, a human exercise study on OA patients with direct sampling from the joint found that 29 patients with knee OA divided in two groups were sufficient for detecting a difference in intra articular IL-10 concentrations
      • Helmark I.C.
      • Mikkelsen U.R.
      • Børglum J.
      • Rothe A.
      • Petersen M.C.H.
      • Andersen O.
      • et al.
      Exercise increases interleukin-10 levels both intraarticularly and peri-synovially in patients with knee osteoarthritis: a randomized controlled trial.
      . Accounting for possible dropouts at 10% we estimated the total need of 32 participants. The study was performed at Bispebjerg Hospital, Copenhagen, Denmark. Participants gave written informed content, and the study was approved by the Ethical Committee of the Capital Region of Denmark (H-17015563) in accordance with the Helsinki Declaration of 1975/83. The trial was registered on ClinicalTrials.gov (NCT03410745) – initially as a single trial on collagen synthesis containing 4 weeks of isotope labelling with deuterated water combined with 4 weeks of exercise intervention. However, due to recruitment issues, the trial was split into two studies, where the 4-week prospective deuterium labelling study on collagen synthesis
      • Jørgensen A.E.M.
      • Agergaard J.
      • Schjerling P.
      • Heinemeier K.M.
      • van Hall G.
      • Kjaer M.
      The regional turnover of cartilage collagen matrix in late-stage human knee osteoarthritis.
      is not a trial. We have deliberately not changed the collagen synthesis registration on ClinicalTrials.gov partly for transparency on the initial single study design, and partly to highlight that recruiting patients for long-term exercise intervention was impossible in our local setting.

      Participants

      Patients were recruited from the outpatient clinic at Bispebjerg Hospital orthopaedic department. Patients with OA of the knee as classified by the American College of Rheumatology
      • Altman R.
      • Asch E.
      • Bloch D.
      • Bole G.
      • Borenstein D.
      • Brandt K.
      • et al.
      Development of criteria for the classification and reporting of osteoarthritis: classification of osteoarthritis of the knee.
      scheduled for a knee replacement surgery by an orthopaedic surgeon were considered for participation. Patients were eligible if they were 40–90 years old, had primary OA, were able to perform physical exercise, had a body mass index (BMI) within 18.5–40 kg/m2, and were non-smokers.
      Exclusion criteria included other diseases of the joint (rheumatoid arthritis, gout, psoriasis arthritis), inflammatory diseases, diabetes, cancer, previous knee surgery, recent trauma to the knee, participation in regular strenuous exercise, or required an interpreter. Furthermore, use of non-steroid anti-inflammatory drugs (NSAIDs), oral corticoid-steroids or intraarticular injection within 3 months, disease modifying anti-rheumatic drugs (DMARDs), chemotherapy, or other immuno-suppressive medicine deemed the patient ineligible.

      Intervention

      Exercise was supervised by a physician and began with a light warm-up on an ergometer bike for 5–10 min. In a leg-press machine (Techno Gym, Cesena, Italy), participants were seated upright with the hip bended 90°. To ensure proper load on the knee of interest, one-legged exercise was performed with the working leg extending from 100° to near fully stretched against a vertical plate, and the resting leg positioned on a horizontal plate. Due to safety concerns of the overnight fasted participants with knee pain, we performed a 5 repetition maximum (RM) test as an indirect measure of maximum strength to calculate the corresponding 1 RM
      • Brzycki M.
      Strength testing – predicting a one-rep max from reps-to-fatigue.
      . Three sets of eight repetitions at 70% 1 RM and three sets of four repetitions at 80% 1 RM were performed at a controlled pace.

      Randomisation

      Participants were randomly allocated 1:1 to one of the two groups using the minimisation software MinimPy version 0.3 (Python Software Foundation, Beaverton, OR, USA)
      • Saghaei M.
      • Saghaei S.
      Implementation of an open-source customizable minimization program for allocation of patients to parallel groups in clinical trials.
      , stratified by sex, age (<70 and 70≤ years), and BMI (18.5–25, 25.1–30, and 30.1–40 kg/m2). Due to the nature of the intervention, the participants were unblinded, however, to the best of our knowledge the outcome measurements cannot be influenced by own will. Surgeons were unaware of the allocated group during surgery, and messenger ribonucleic acid (mRNA) and GAG data were obtained blinded.

      Tissue collection

      During knee replacement surgery, the leftover tibia plateau was collected and immediately transferred on ice to the laboratory. From the medial tibia condyle, several full-thickness cartilage biopsies (30–150 mg) were sampled centrally and under the meniscus. By cutting with a scalpel, the subchondral bone was removed from biopsies. From osteophytes, slices (10–120 mg) were sampled using a scalpel. The cartilage was snap-frozen within 15 min from surgery in liquid nitrogen and kept on −80°C until further analysis.

      Gene expression analysis

      From the cartilage biopsies and slices, mRNA expression was measured as previously described
      • Jørgensen A.E.M.
      • Agergaard J.
      • Schjerling P.
      • Heinemeier K.M.
      • van Hall G.
      • Kjaer M.
      The regional turnover of cartilage collagen matrix in late-stage human knee osteoarthritis.
      . In brief, total RNA was isolated from cartilage and osteophyte tissue by homogenisation in a guanidine thiocyanate/phenol solution, before being precipitated, and further extracted using the RNeasy MinElute Cleanup Kit (Qiagen, Hilden, Germany) as per manufacture's protocol. The total RNA concentration was determined by the RiboGreen assay (Molecular Probes Inc., Eugene, OR, USA). 50 ng total RNA was reversed transcribed to complementary deoxyribonucleic acid (cDNA), and target mRNAs were quantified using real-time reverse transcriptase polymerase chain reaction (RT-PCR). Primers are given in Supplementary Table S1. Ct values were converted to cDNA molecules using a standard curve of DNA oligos (Ultramer, Integrated DNA Technologies Inc., Leuven, Belgium) with a DNA sequence corresponding to the expected PCR product. No sample specific inhibition was seen comparing efficiency of the individual sample PCR reactions with the standard curve using the three-parameter logistic model (LRE)
      • Tellinghuisen J.
      Estimating real-time qPCR amplification efficiency from single-reaction data.
      . Expressed values were normalised to internal control of large ribosomal protein P0 (RPLP0), which was then validated (no difference between groups and low variation) using glyceraldehyde 3-phosphate dehydrogenase (GAPDH). The RPLP0-normalized values are expressed as ratio to the total geometric mean for all samples of the central region.

      Extraction and quantification of glycosaminoglycans

      The central and sub-meniscal cartilage biopsies (all osteophyte tissue was used for RNA) were sliced, weighed, freeze-dried, and weighed again (dry weight). Next, GAGs were extracted by overnight treatment with hyaluronidase (H3506, Sigma, Darmstadt, Germany) [5 U/ml in 0.05 M sodium acetate and 0.15 M NaCl (pH = 6)] at 37°C. The samples were centrifuged, and the supernatant collected, aliquoted, and stored at −80°C until further analyses.
      GAG content was measured using the dimethyl-methylene blue (DMMB) assay. The hyaluronidase-treated cartilage biopsies were washed with isotonic NaCl, and then treated overnight at 60°C with 0.125 mg/ml papain (P3125, Sigma, Darmstadt, Germany) in PBE buffer (100 mM sodium phosphate buffer, 10 mM Na2-EDTA) with 10 mM l-cysteine at pH = 6.5. Both the resulting papain digest and the hyaluronidase supernatant were diluted in PBE buffer and mixed with a 1,9-dimethylmethylene blue (341088, Sigma, Darmstadt, Germany) solution containing 38 μM DMMB in 40 mM NaCl, 40 mM glycine at pH = 3. The absorbances were read at 595 nm and 540 nm wavelengths (subtracted) and compared to a known standard curve of chondroitin sulphate (CS) C (C4384, Sigma, Darmstadt, Germany).
      The detailed composition of the hyaluronidase extracted GAGs was determined using immuno-assays of antibodies visualised by chemo-luminescence. On positively charged nylon membranes (Hybond N+, RPN203B, Amersham, Broendby, Denmark), 2 μl was dotted in duplicate and incubated in blocking buffer [5% skimmed milk powder in tris-buffered saline (TBS) with 0.05% Tween20 (TBS-T) (Sigma, Darmstadt, Germany)] for 1 h at room temperature with diluted primary antibodies [chondroitin sulphate (CS) (Sigma–Aldrich, Broendby, Denmark), keratan sulphate (KS) (Santa Cruz, clone 4D3/D10, Dallas, TX, US). Next, the membranes were washed in TBS-T, before incubation with horseradish peroxidase-conjugated secondary antibodies [anti-mouse and anti-goat (DAKO, Glostrup, Denmark)] diluted in blocking buffer (TBS-T) for 1 h at room temperature. Membranes were washed with TBS-T and TBS, before incubation with SuperSignal Femto Reagent (Pierce, Thermo Fischer Scientific, Waltham, MA, USA) for visualisation and imaged by CCD camera (Bio-Rad Laboratories, Hercules, CA, USA). Chemiluminescence was quantified with densitometry after normalizing to background with ImageJ software and normalised to initial dry weight of the cartilage sample. All samples are expressed relative to the geometric mean of the central region.

      Statistics

      Statistical analyses were made in GraphPad Prism version 8.4.3 (GraphPad Software, San Diego, CA, USA) on logarithmic values tested for Gaussian distribution using Shapiro–Wilk's test. Effects of exercise intervention and regional location was compared by mixed-effect analysis with patient as random effect and location + intervention + location∗intervention as fixed effects using Tukey's correction for multiple comparisons test. This mixed model uses a compound symmetry covariance matrix and is fit using Restricted Maximum Likelihood (REML). Normality and homoscedasticity were evaluated by residual- and QQ-plots. Statistical significance was defined as P-values < 0.05 for all analyses with 95% confidence interval (CI) of difference shown in brackets. Values are reported as means ± SD for participant characteristics and RNA extraction; geometric means ± back transformed SEM for gene expression and GAG immuno-assay data; mean with 95% CI for GAG measurements data; and median with 95% CI for 5RM data. An overview of the statistical analyses including 95% CI of differences is reported in Table S3.

      Results

      The trial was completed as the 32 participants were recruited between May 2018 and September 2019 (Fig. 1). There are similar characteristics in the two groups (Table I). No harms were reported during the exercise intervention.
      Fig. 1
      Fig. 1CONSORT diagram. The flow of participants through the study.
      Table IParticipants' characteristics
      ParameterControl (n = 15)Exercise (n = 16)
      Age (years ± SD)67.8 ± 8.069.5 ± 7.2
      Body mass index (kg/m2 ± SD)29.7 ± 3.229.2 ± 5.4
      Sex (percentage female)67%50%
      Kellgren–Lawrence grade (score [range])3.6 [3–4]3.4 [3–4]
      5 repetition maximum (kg median with 95% CI)N/A110 [90–120]
      Time from exercise to tissue freezing (h:mm ± SD)N/A3:24 ± 0:37
      N/A: Not applicable. Values are means if not stated otherwise.

      Gene expression and loading intervention

      In the intervention group, the loaded weight reported as the 5RM test showed a median of 110 kg [95% CI: 90–120 kg] (Fig. S1). For both groups, low amounts of RNA were extracted for articular cartilage with means ± SD of 3.1 ± 1.6, 2.6 ± 1.5, 2.9 ± 1.6, and 3.2 ± 2.0 ng/mg tissue, while more RNA was extracted for osteophytes 16.7 ± 6.8 and 17.3 ± 28.8 ng/mg tissue (Fig. S2). Four samples were of insufficient quality (two in each group), and in total, 23 genes were assessed of which 13 were detected sufficiently for analyses (Table S2). We found no difference in GAPDH, thus validating the use of RPLP0 for normalisation [Fig. 2(A)]. On gene expression, the exercise intervention marked with striped fill pattern on bars in Fig. 2 did not show a significant effect on any target. As it was unexpected that no targets were affected by acute loading, the immediate early genes fos proto-oncogene (FOS) and early growth response 1 (EGR1) were subsequently analysed. They both showed a minor but statistical difference after the loading intervention for FOS [95% CI: 0.21 to 0.99; P = 0.048] and for EGR1 [0.24 to 0.95; P = 0.034], supporting that late-stage OA cartilage is able to respond to the acute loading intervention.
      Fig. 2
      Fig. 2mRNA expression of cartilage tissues. The fold change of every region is shown for sufficiently detected targets relative to central region normalised to RPLP0. Regions are shown in white (central), light grey (submenisc), and dark grey (osteophyte) colours and the exercise intervention by striped fill pattern on top. A: Matrix components and breakdown enzymes. B: Growth factors. Values are individual data points and bars represent geometric means with 95% CI (back transformed) on a logarithmic scale (n = 10–16). Regional differences combining control and exercise together for every region are indicated at P < 0.05 by “a” for central compared to submenisc, “b” for central compared to osteophytes, and “c” submenisc compared to osteophytes. See for complete statistical details.

      Regional gene expression

      In Fig. 2, the submeniscal cartilage region (light grey coloured bars) and osteophytes (dark grey coloured bars) are compared to the central condyle region (white coloured bars). Compared to the central region we found significantly higher fold expression of matrix metalloproteinase 3 (MMP-3) at 2.9 fold [1.7 to 5.0; P < 0.0001], MMP-13 at 3.9 fold [1.4 to 10.7; P = 0.0054], insulin-like growth factor 1 (IGF-1Ea) at 2.6 fold [1.2 to 5.6; P = 0.0127], and connective tissue growth factor (CTGF) at 1.8 fold [1.3 to 2.4; P = 0.0002], combined with lower fold expression of lubricin (PRG4) at 0.4 fold [0.2 to 0.7; P = 0.0004] and cartilage oligomeric protein (COMP) at 0.6 fold [0.4 to 0.9; P = 0.004] in the submeniscal region (Fig. 2).
      Compared to the central region of articular cartilage, osteophytes showed higher fold expression of MMP-1 at 4.8 fold [2.6 to 9.2; P < 0.0001], 78.1 fold MMP-13 [25.8 to 236.1; P < 0.0001], 35.1 fold IGF-1Ea [15.0 to 82.3; P < 0.0001], 2.2 fold TGF- β2 [1.2 to 4.1; P = 0.0067] and 2.8 fold TGF-β3 [1.7 to 4.6; P < 0.0001] together with lower fold expression of aggrecan (ACAN) at 3.2 fold [2.2 to 4.6; P < 0.0001], 0.23 fold COMP [0.15 to 0.35; P < 0.0001], and 0.5 fold fibroblast growth factor 2 (FGF-2) [0.4 to 0.7; P < 0.0001] (Fig. 2).

      Glycosaminoglycans

      From the central region, hyaluronidase treatment extracted a mean of 17.5 μg/mg dry weight (d.w.) GAGs compared to 10.9 μg/mg d.w. from submeniscal regional [1.5 to 11.7; P = 0.012] with no effect of exercise intervention [Fig. 3(A)]. Levels of GAGs remaining in the cartilage biopsies were not significantly different with 95.9 μg/mg d.w. centrally and 107.5 μg/mg d.w. submeniscally and with no effect of exercise intervention [Fig. 3(B)]. Total GAGs were not significantly different with 113.4 μg/mg d.w. centrally and 118.4 μg/mg d.w. submeniscally with no effect of exercise intervention [Fig. 3(C)]. In the supernatant after hyaluronidase treatment, CS content showed higher extractability after loading [1.01 to 2.46; P = 0.0486] with higher amount extracted to the supernatant from submeniscal region compared to central region [1.24 to 2.41; P = 0.0023] [Fig. 4(A)]. KS showed significant interaction [0.81 to 0.99; P = 0.045) (Fig. 4(B)).
      Fig. 3
      Fig. 3Dimethyl-methylene blue assay. A: GAGs extracted to supernatant by overnight treatment with hyaluronidase. B: GAGs remaining in cartilage biopsies after the hyaluronidase treatment. Each data point is means of duplicate analyses and error bars represent means with 95% CI.
      Fig. 4
      Fig. 4Immuno-assay. GAG measurements in the supernatant after hyaluronidase treatment. All samples are expressed relative to the geometric mean of the central region. A: Keratan sulphate B: Chondroitin sulphate. Each data point is means of duplicate analyses and error bars represent geometric means with 95% CI (back transformed) shown on logarithmic scale.

      Discussion

      Acute loading has minor influence on chondrocyte gene expression

      This randomised controlled trial examined the effect of acute exercise in vivo on gene expression from human OA chondrocytes from three regions of the tibia plateau. As we found no cartilage genes responsive to the intervention, we examined two immediate early genes, which both did respond to the loading intervention. This is somewhat in contrast to upregulation of expression in multiple genes that has been reported from a human in vitro study
      • Scholtes S.
      • Krämer E.
      • Weisser M.
      • Roth W.
      • Luginbühl R.
      • Grossner T.
      • et al.
      Global chondrocyte gene expression after a single anabolic loading period: time evolution and re-inducibility of mechano-responses.
      . On the other hand, this is the first study to examine the effect of exercise on gene expression directly in joint cartilage in humans in vivo, and thus, we are unable to compare our results directly. Our data suggest that in humans with OA, chondrocyte gene expression of cartilage relevant genes is in fact not regulated in a detectable manner by this single bout of exercise setup within the given timeframe.
      In humans subjected to acute exercise, biomarker studies on serum
      • Roberts H.M.
      • Law R.J.
      • Thom J.M.
      The time course and mechanisms of change in biomarkers of joint metabolism in response to acute exercise and chronic training in physiologic and pathological conditions.
      and SF
      • Mazor M.
      • Best T.M.
      • Cesaro A.
      • Lespessailles E.
      • Toumi H.
      Osteoarthritis biomarker responses and cartilage adaptation to exercise: a review of animal and human models.
      exist, and although COMP has been widely examined for breakdown of cartilage, COMP is also expressed in synovial cells
      • Tseng S.
      • Reddi A.H.
      • Di Cesare P.E.
      Cartilage oligomeric matrix protein (COMP): a biomarker of arthritis.
      , complicating the measurements on SF. We did not detect any effect of loading on chondrocyte gene expression. Decreased levels of SF COMP has been reported after acute exercise
      • Helmark I.C.
      • Petersen M.C.H.
      • Christensen H.E.
      • Kjaer M.
      • Langberg H.
      Moderate loading of the human osteoarthritic knee joint leads to lowering of intraarticular cartilage oligomeric matrix protein.
      , although baseline measurements were made 3 months earlier, and an effect of time cannot be ruled out. Further, with acute exercise, we did not find any effect upon aggrecan expression. Similarly, in OA chondrocytes, loading did not change aggrecan expression, whereas an increase was observed in healthy chondrocytes
      • Millward-Sadler S.J.
      • Wright M.O.
      • Davies L.W.
      • Nuki G.
      • Salter D.M.
      Mechanotransduction via integrins and interleukin-4 results in altered aggrecan and matrix metalloproteinase 3 gene expression in normal, but not osteoarthritic, human articular chondrocytes.
      . Using delayed gadolinium-enhanced magnetic resonance imaging of cartilage (dGEMRIC) as a surrogate marker for GAG content, long-term exercise provided no change in tibiofemoral cartilage with mild OA
      • Multanen J.
      • Nieminen M.T.
      • Häkkinen A.
      • Kujala U.M.
      • Jämsä T.
      • Kautiainen H.
      • et al.
      Effects of high-impact training on bone and articular cartilage: 12-month randomized controlled quantitative MRI study.
      • Munukka M.
      • Waller B.
      • Rantalainen T.
      • Häkkinen A.
      • Nieminen M.T.
      • Lammentausta E.
      • et al.
      Efficacy of progressive aquatic resistance training for tibiofemoral cartilage in postmenopausal women with mild knee osteoarthritis: a randomised controlled trial.
      • Multanen J.
      • Rantalainen T.
      • Kautiainen H.
      • Ahola R.
      • Jämsä T.
      • Nieminen M.T.
      • et al.
      Effect of progressive high-impact exercise on femoral neck structural strength in postmenopausal women with mild knee osteoarthritis: a 12-month RCT.
      . Thus, in OA, an exercise-induced increase in GAG content seems unattainable.
      In the current study, no change in lubricin gene expression of chondrocytes with exercise was found. Synoviocytes and superficial zone chondrocytes express lubricin, which reduces the sheer stress to the surface and protects against chondrocyte apoptosis
      • Waller K.A.
      • Zhang L.X.
      • Elsaid K.A.
      • Fleming B.C.
      • Warman M.L.
      • Jay G.D.
      Role of lubricin and boundary lubrication in the prevention of chondrocyte apoptosis.
      . In a rat OA model, lubricin injections inhibited toll like receptors
      • Iqbal S.M.
      • Leonard C.
      • C Regmi S.
      • De Rantere D.
      • Tailor P.
      • Ren G.
      • et al.
      Lubricin/proteoglycan 4 binds to and regulates the activity of toll-like receptors in vitro.
      , thus providing cartilage protection. In healthy rodents, low and moderate intensity running increased lubricin gene expression in chondrocytes after 4–8 weeks
      • Ni G.X.
      • Lei L.
      • Zhou Y.Z.
      • Su J.
      • Schumacher B.
      • Lindley K.
      • et al.
      Intensity-dependent effect of treadmill running on lubricin metabolism of rat articular cartilage.
      ,
      • Ogawa H.
      • Kozhemyakina E.
      • Hung H.H.
      • Grodzinsky A.J.
      • Lassar A.B.
      Mechanical motion promotes expression of Prg4 in articular cartilage via multiple CREB-dependent, fluid flow shear stress-induced signaling pathways.
      . A single study has measured the loading effect in vivo in rats with a damaged joint, mimicking early post-traumatic OA. Here, an ACL lesion reduced lubricin expression in chondrocytes, which was then further reduced by forced running
      • Elsaid K.A.
      • Zhang L.
      • Waller K.
      • Tofte J.
      • Teeple E.
      • Fleming B.C.
      • et al.
      The impact of forced joint exercise on lubricin biosynthesis from articular cartilage following ACL transection and intra-articular lubricin's effect in exercised joints following ACL transection.
      . Thus, in OA, an exercise-induced lubricin secretion from chondrocytes are uncertain at best. Further, a study applying in vivo stretching of synovial fibroblasts, found increased lubricin gene expression in healthy but not OA donors
      • Jamal J.
      • Roebuck M.M.
      • Lee S.Y.
      • Frostick S.P.
      • Abbas A.A.
      • Merican A.M.
      • et al.
      Modulation of the mechanical responses of synovial fibroblasts by osteoarthritis-associated inflammatory stressors.
      .
      In our present study, similar to others
      • Dunn S.L.
      • Soul J.
      • Anand S.
      • Schwartz J.M.
      • Boot-Handford R.P.
      • Hardingham T.E.
      Gene expression changes in damaged osteoarthritic cartilage identify a signature of non-chondrogenic and mechanical responses.
      , we did not detect any mRNA from interleukins (ILs). This raises the question of whether cells from other regions (i.e., synovial lining or subchondral bone) could in fact be the population responsible for IL production. Another possibility is differing tissue preparation, as one study found a marked increase in IL expression in vitro compared to in vivo
      • Tsuchida A.I.
      • Beekhuizen M.
      • t Hart M.C.
      • Radstake T.R.D.J.
      • Dhert W.J.A.
      • Saris D.B.F.
      • et al.
      Cytokine profiles in the joint depend on pathology, but are different between synovial fluid, cartilage tissue and cultured chondrocytes.
      . In patients with moderate OA (Kellgren–Lawrence 2–3), increased SF and perisynovial IL-10 was found in the recovery after acute exercise, however, IL-10 levels were not measured before exercise
      • Helmark I.C.
      • Mikkelsen U.R.
      • Børglum J.
      • Rothe A.
      • Petersen M.C.H.
      • Andersen O.
      • et al.
      Exercise increases interleukin-10 levels both intraarticularly and peri-synovially in patients with knee osteoarthritis: a randomized controlled trial.
      . In our setup with direct visualisation of the joint capsule during surgery, too few SF samples were collectable to perform feasible SF analysis and further conclusions on IL-10 with exercise remains speculative.
      In healthy human chondrocytes that are exposed to load, MMP-1 is either unchanged
      • Millward-Sadler S.J.
      • Wright M.O.
      • Davies L.W.
      • Nuki G.
      • Salter D.M.
      Mechanotransduction via integrins and interleukin-4 results in altered aggrecan and matrix metalloproteinase 3 gene expression in normal, but not osteoarthritic, human articular chondrocytes.
      , or decreases together with MMP-3 and MMP-13 due to increasing IL-4
      • Salter D.M.
      • Millward-Sadler S.J.
      • Nuki G.
      • Wright M.O.
      Differential responses of chondrocytes from normal and osteoarthritic human articular cartilage to mechanical stimulation.
      • Leong D.J.
      • Li Y.H.
      • Gu X.I.
      • Sun L.
      • Zhou Z.
      • Nasser P.
      • et al.
      Physiological loading of joints prevents cartilage degradation through CITED2.
      • He Z.
      • Leong D.J.
      • Xu L.
      • Hardin J.A.
      • Majeska R.J.
      • Schaffler M.B.
      • et al.
      CITED2 mediates the cross-talk between mechanical loading and IL-4 to promote chondroprotection.
      . However, in OA there was no change in MMP-1 or MMP-3
      • Millward-Sadler S.J.
      • Wright M.O.
      • Davies L.W.
      • Nuki G.
      • Salter D.M.
      Mechanotransduction via integrins and interleukin-4 results in altered aggrecan and matrix metalloproteinase 3 gene expression in normal, but not osteoarthritic, human articular chondrocytes.
      ,
      • Salter D.M.
      • Millward-Sadler S.J.
      • Nuki G.
      • Wright M.O.
      Differential responses of chondrocytes from normal and osteoarthritic human articular cartilage to mechanical stimulation.
      , similar to our data. Further, in the present study we barely detected MMP-1, consistent with one study
      • Bau B.
      • Gebhard P.M.
      • Haag J.
      • Knorr T.
      • Bartnik E.
      • Aigner T.
      Relative messenger RNA expression profiling of collagenases and aggrecanases in human articular chondrocytes in vivo and in vitro.
      but in contrast to others
      • Snelling S.
      • Rout R.
      • Davidson R.
      • Clark I.
      • Carr A.
      • Hulley P.A.
      • et al.
      A gene expression study of normal and damaged cartilage in anteromedial gonarthrosis, a phenotype of osteoarthritis.
      ,
      • Jessberger S.
      • Högger P.
      • Genest F.
      • Salter D.M.
      • Seefried L.
      Cellular pharmacodynamic effects of Pycnogenol® in patients with severe osteoarthritis: a randomized controlled pilot study.
      . Thus, our human exercise data are in line with the in vitro studies, and MMPs are not influenced within our setup. Finally, ADAMTS-4 was undetected, and ADAMTS-5 was found in too low amount contrary to others
      • Bau B.
      • Gebhard P.M.
      • Haag J.
      • Knorr T.
      • Bartnik E.
      • Aigner T.
      Relative messenger RNA expression profiling of collagenases and aggrecanases in human articular chondrocytes in vivo and in vitro.
      ,
      • Jessberger S.
      • Högger P.
      • Genest F.
      • Salter D.M.
      • Seefried L.
      Cellular pharmacodynamic effects of Pycnogenol® in patients with severe osteoarthritis: a randomized controlled pilot study.
      . From our data, additional conclusions on ADAMTS are speculative. Currently, exercise did not lead to increased expression of breakdown enzymes nor cytokines, and consequently, in patients with late-stage OA, the performed exercise intervention appears safe for cartilage.
      One limitation of the study is that only 23% of the patients screened were eligible for participation, thus limiting the generalisability of the study. However, the comprehensive eligibility criteria optimised the possibility of identifying an effect of exercise, and the trial participants are believed to represent the most likely to respond to exercise. Hence, an effect in post-traumatic OA cannot be excluded. Another limitation is the intervention itself. The effect of exercise may emerge later – although the genes of interest are reported to change after the 3-h timeframe
      • Scholtes S.
      • Krämer E.
      • Weisser M.
      • Roth W.
      • Luginbühl R.
      • Grossner T.
      • et al.
      Global chondrocyte gene expression after a single anabolic loading period: time evolution and re-inducibility of mechano-responses.
      . Our intervention could be inadequate as a stimulus as the participants did not perform an all-out resistance exercise bout. However, the performed exercise uses moderate to hard load, and thus, the cartilage should be loaded sufficiently, which is supported by the loading effect in FOS and EGR1. Further, despite our participants being overnight fasted and untrained, all completed the protocol without any harm, as has been reported previously
      • Skoffer B.
      • Dalgas U.
      • Maribo T.
      • Søballe K.
      • Mechlenburg I.
      No exacerbation of knee joint pain and effusion following preoperative progressive resistance training in patients scheduled for total knee arthroplasty: secondary analyses from a randomized controlled trial.
      , thus making the intervention clinically possible and applicable. Finally, a single bout of resistance exercise might not be enough, as months of exercise tended to show a beneficial response on cartilage biomarkers
      • Bricca A.
      • Struglics A.
      • Larsson S.
      • Steultjens M.
      • Juhl C.B.
      • Roos E.M.
      Impact of exercise therapy on molecular biomarkers related to cartilage and inflammation in people at risk of, or with established, knee osteoarthritis: a systematic review and meta-analysis of randomized controlled trials.
      . Potentially, the slight and insignificant difference in sex distribution between the two groups would have influenced the results as women in general are found to have more severe OA. In counter to this argument, the analysis of KL-score between males and females in the present study were similar. Further, in the present study we did determine the outcome result before and after exercise, but due to the inability to sample more than one cartilage sample within a single individual, we cannot specifically know what the true baseline values in cartilage were. Given this drawback in human studies we still find it important to measure cartilage characteristic as done in the present study.

      Regional differences in chondrocyte gene expression

      As the submeniscal area on the medial plateau experience higher load than central
      • Wong M.
      • Carter D.R.
      Articular cartilage functional histomorphology and mechanobiology: a research perspective.
      , we sampled from both regions. Our data showed increased MMP-3, MMP-13, IGF-1Ea, and CTGF, and decreased lubricin and COMP in the submeniscal area. Work from our lab confirms similar regional differences for MMP-3 and lubricin
      • Jørgensen A.E.M.
      • Agergaard J.
      • Schjerling P.
      • Heinemeier K.M.
      • van Hall G.
      • Kjaer M.
      The regional turnover of cartilage collagen matrix in late-stage human knee osteoarthritis.
      . Thus, it appears that breakdown, fibrosis, and reduced anabolism dominate submeniscally. In the current study, cartilage fibrillation was universally found centrally, whereas submeniscal cartilage conditions varied from appearing healthy to a complete loss. This appearance is also reported by others, who showed increased MMP-3, MMP-13, and IGF-1Ea in undamaged cartilage (corresponding to submeniscal) from the medial tibia plateau
      • Snelling S.
      • Rout R.
      • Davidson R.
      • Clark I.
      • Carr A.
      • Hulley P.A.
      • et al.
      A gene expression study of normal and damaged cartilage in anteromedial gonarthrosis, a phenotype of osteoarthritis.
      , and decreased COMP
      • Dunn S.L.
      • Soul J.
      • Anand S.
      • Schwartz J.M.
      • Boot-Handford R.P.
      • Hardingham T.E.
      Gene expression changes in damaged osteoarthritic cartilage identify a signature of non-chondrogenic and mechanical responses.
      ,
      • Snelling S.
      • Rout R.
      • Davidson R.
      • Clark I.
      • Carr A.
      • Hulley P.A.
      • et al.
      A gene expression study of normal and damaged cartilage in anteromedial gonarthrosis, a phenotype of osteoarthritis.
      . From femur condyles, results varies: some studies found differences in MMP-3
      • Dunn S.L.
      • Soul J.
      • Anand S.
      • Schwartz J.M.
      • Boot-Handford R.P.
      • Hardingham T.E.
      Gene expression changes in damaged osteoarthritic cartilage identify a signature of non-chondrogenic and mechanical responses.
      , and MMP-13
      • Sato T.
      • Konomi K.
      • Yamasaki S.
      • Aratani S.
      • Tsuchimochi K.
      • Yokouchi M.
      • et al.
      Comparative analysis of gene expression profiles in intact and damaged regions of human osteoarthritic cartilage.
      , while others did not
      • Yagi R.
      • McBurney D.
      • Laverty D.
      • Weiner S.
      • Horton W.E.
      Intrajoint comparisons of gene expression patterns in human osteoarthritis suggest a change in chondrocyte phenotype.
      ,
      • Brew C.J.
      • Clegg P.D.
      • Boot-Handford R.P.
      • Andrew J.G.
      • Hardingham T.
      Gene expression in human chondrocytes in late osteoarthritis is changed in both fibrillated and intact cartilage without evidence of generalised chondrocyte hypertrophy.
      . We found no difference in expression of aggrecan or COL2A1, similar to many
      • Yagi R.
      • McBurney D.
      • Laverty D.
      • Weiner S.
      • Horton W.E.
      Intrajoint comparisons of gene expression patterns in human osteoarthritis suggest a change in chondrocyte phenotype.
      • Brew C.J.
      • Clegg P.D.
      • Boot-Handford R.P.
      • Andrew J.G.
      • Hardingham T.
      Gene expression in human chondrocytes in late osteoarthritis is changed in both fibrillated and intact cartilage without evidence of generalised chondrocyte hypertrophy.
      • Fukui N.
      • Ikeda Y.
      • Ohnuki T.
      • Tanaka N.
      • Hikita A.
      • Mitomi H.
      • et al.
      Regional differences in chondrocyte metabolism in osteoarthritis: a detailed analysis by laser capture microdissection.
      , although decreased aggrecan has been found centrally
      • Dunn S.L.
      • Soul J.
      • Anand S.
      • Schwartz J.M.
      • Boot-Handford R.P.
      • Hardingham T.E.
      Gene expression changes in damaged osteoarthritic cartilage identify a signature of non-chondrogenic and mechanical responses.
      .
      Osteophytes are bony outgrowths covered by fibro-cartilage
      • van der Kraan P.M.
      • van den Berg W.B.
      Osteophytes: relevance and biology.
      , and compared to articular cartilage, we found decreased aggrecan, COMP, and FGF-2 and increased MMP-1, MMP-13, TGF-β2, TGF-β3, and IGF-1Ea compatible with osteophyte growth and ossification capability. Very limited work has been done in humans on chondrocyte gene expression from osteophytes, but similar increases in TGF-β and unchanged FGF-2 have been found
      • Uchino M.
      • Izumi T.
      • Tominaga T.
      • Wakita R.
      • Minehara H.
      • Sekiguchi M.
      • et al.
      Growth factor expression in the osteophytes of the human femoral head in osteoarthritis.
      , while IGF-1 was found expressed in osteophytes, the level was not compared to articular cartilage
      • Middleton J.
      • Arnott N.
      • Walsh S.
      • Beresford J.
      Osteoblasts and osteoclasts in adult human osteophyte tissue express the mRNAs for insulin-like growth factors I and II and the type 1 IGF receptor.
      . Thus, osteophyte gene expression differs from articular cartilage, and the pattern conforms with osteophyte growth, remodelling, and ossification capability supporting the underlying bony growth. Hence, care should be made during cartilage tissue sampling of an entire joint surface for RNA analysis.

      Glycosaminoglycans are slightly affected by acute loading and regionally different

      The total content of GAGs in cartilage and SF is unlikely to change sufficiently for robust detection, however, we examined the potential influence of acute exercise upon GAG dynamics. We used hyaluronidase in rather low amounts to assess the extractable amounts of GAGs as a marker for loosened matrix. We examined the GAG composition of the hyaluronidase supernatant further and found minor effects of exercise with increased extractability of CS and an interaction of KS. Thus, acute loading could potentially affect some GAG dynamics, and we cannot exclude that more pronounced dynamic changes would have been observed if we had sampled sometime later than the 3.5 h.
      Regionally, more GAGs were extracted from the central region, which showed macroscopic fibrillation universally, and thus, likely having more loosely bound matrix. Others too, have shown less GAGs left in the damaged medial femur compared to the lateral
      • Dunn S.L.
      • Soul J.
      • Anand S.
      • Schwartz J.M.
      • Boot-Handford R.P.
      • Hardingham T.E.
      Gene expression changes in damaged osteoarthritic cartilage identify a signature of non-chondrogenic and mechanical responses.
      . Aggrecan is linked to HA and contains a smaller KS rich section followed by two large CS rich sections
      • Jørgensen A.E.M.
      • Kjær M.
      • Heinemeier K.M.
      The effect of aging and mechanical loading on the metabolism of articular cartilage.
      , and as aggrecan gene expression is regionally similar, breakdown could be the dominant factor. Finally, more CS were extracted submeniscally with equal KS extractability. Combined with the results on regional gene expression, a careful suggestion could be that matrix breakdown by MMPs, ossification due to IGF-1 and CTGF signalling, and reduced cartilage matrix formation of COMP and lubricin dominates under the meniscus in areas where cartilage still exists. Enzymatic degradation of aggrecan would release CS rich sections sequentially while keeping most KS attached
      • Vynios D.H.
      Metabolism of cartilage proteoglycans in health and disease.
      • Dudhia J.
      • Davidson C.M.
      • Wells T.M.
      • Vynios D.H.
      • Hardingham T.E.
      • Bayliss M.T.
      Age-related changes in the content of the C-terminal region of aggrecan in human articular cartilage.
      • Struglics A.
      • Hansson M.
      MMP proteolysis of the human extracellular matrix protein aggrecan is mainly a process of normal turnover.
      , and in support, others have found less CS left in areas with more damage
      • Yagi R.
      • McBurney D.
      • Laverty D.
      • Weiner S.
      • Horton W.E.
      Intrajoint comparisons of gene expression patterns in human osteoarthritis suggest a change in chondrocyte phenotype.
      . As submeniscal cartilage varied from healthy looking to complete subchondral induration, this process seems the likely future outcome for all submeniscal cartilage had surgery not been performed. Our data suggest regional differences and a possible acute effect of in vivo loading despite the absolute amounts of GAGs involved being rather small, and thus, further studies are necessary to uncover exercise-induced GAG dynamics.
      Taken together, we found minor effects on GAG extractability. Thus, the beneficial response found in healthy cartilage appears not to be reflected in late-stage human OA within the given timeframe. This suggests clinical improvements of exercise are likely a result of the known extra-articular musculo-skeletal adaptations. Nevertheless, from the lack of a major effect, exercise appears safe in late-stage OA, albeit not beneficial per se. Finally, regional differences in chondrocyte gene expression indicate an unfavourable environment under the meniscus, possibly leading to subchondral exposure frequently observed is this region. In conclusion, acute exercise did not influence expression of cartilage relevant proteins and had minor effects on GAG dynamics, indicative of a relatively stable pathological cartilage in OA. The major regional cartilage differences demonstrated in this study indicates a chronic mechanical influence upon joint cartilage in humans.

      Contributions

      AJ, PS & MK concepted and designed the study. AJ acquired data and drafted the article. PS & AJ performed the analysis and interpretation of RNA data and DMMB assay. BDV performed the analysis of immune assay. All authors revised the manuscript critically for important intellectual content and approved the final version to be submitted. AJ and MK take responsibility for the integrity of the work as a whole, from inception to finished article.

      Funding

      This study has received financial support from The Danish Rheumatism Association (AJ & MK), Bispebjerg Hospital Research Grant (AJ), The Nordea Foundation (Healthy Aging Grant) (MK), The Novo Nordisk Foundation (MK), and The Lundbeck Foundation (AJ & MK).

      Declaration of competing interest

      BDV and JR are owners of GLX Analytix, a company that studies glycocalyx shedding in relation to disease. GLX Analytix has had no influence on the design, data handling and decision to publish. The remaining authors report no conflicts of interest.

      Acknowledgements

      The authors kindly acknowledge technical assistance from Anja Jokipii-Utzon and Ann-Christina Reimann, support from the participating surgeons and other staff at orthopaedic department M50 and the outpatient clinic.

      Appendix A. Supplementary data

      The following are the Supplementary data to this article.
      Fig. S1
      Fig. S1Loaded weight. The amount of weight lifted during the one-legged leg-press reported as the 5 repeated maximum test. Values are individual data points and bars represent median with 95% CI.
      Fig. S2
      Fig. S2RNA extraction. The total amount of extracted RNA from each location of cartilage tissue. Values are individual data points and bars represent means with SD.
      Fig. S3
      Fig. S3Raw mRNA expression of cartilage tissues. The fold change of every region is shown for sufficiently detected targets normalised to RPLP0. Regions are shown in white (central), light grey (submenisc), and dark grey (osteophyte) colours and the exercise intervention by striped fill pattern on top. A: Matrix components and breakdown enzymes. B: Growth factors. Values are individual data points and bars represent geometric means with 95% CI (back transformed) on a logarithmic scale (n = 10–16).

      References

        • Syx D.
        • Tran P.B.
        • Miller R.E.
        • Malfait A.M.
        Peripheral mechanisms contributing to osteoarthritis pain.
        Curr Rheumatol Rep. 2018; 20https://doi.org/10.1007/s11926-018-0716-6
        • Shane Anderson A.
        • Loeser R.F.
        Why is osteoarthritis an age-related disease?.
        Best Pract Res Clin Rheumatol. 2010; 24: 15-26https://doi.org/10.1016/j.berh.2009.08.006
        • Glyn-Jones S.
        • Palmer A.J.R.
        • Agricola R.
        • Price A.J.
        • Vincent T.L.
        • Weinans H.
        • et al.
        Osteoarthritis.
        Lancet. 2015; 386: 376-387https://doi.org/10.1016/S0140-6736(14)60802-3
        • Kolasinski S.L.
        • Neogi T.
        • Hochberg M.C.
        • Oatis C.
        • Guyatt G.
        • Block J.
        • et al.
        2019 American college of rheumatology/arthritis foundation guideline for the management of osteoarthritis of the hand, hip, and knee.
        Arthritis Care Res. 2020; 72: 149-162https://doi.org/10.1002/acr.24131
        • Fransen M.
        • McConnell S.
        • Harmer A.R.
        • Van der Esch M.
        • Simic M.
        • Bennell K.L.
        Exercise for osteoarthritis of the knee.
        in: Fransen M. Cochrane Database of Systematic Reviews Published online January 9, 2015. 2015https://doi.org/10.1002/14651858.CD004376.pub3
        • Bricca A.
        Exercise does not “wear down my knee”: systematic reviews and meta-analyses.
        Br J Sports Med. 2018; 52: 1591-1592https://doi.org/10.1136/bjsports-2018-099705
        • Bricca A.
        • Struglics A.
        • Larsson S.
        • Steultjens M.
        • Juhl C.B.
        • Roos E.M.
        Impact of exercise therapy on molecular biomarkers related to cartilage and inflammation in people at risk of, or with established, knee osteoarthritis: a systematic review and meta-analysis of randomized controlled trials.
        Arthritis Care Res. 2018; (Published online)https://doi.org/10.1002/acr.23786
        • Wong M.
        • Carter D.R.
        Articular cartilage functional histomorphology and mechanobiology: a research perspective.
        Bone. 2003; 33: 1-13https://doi.org/10.1016/S8756-3282(03)00083-8
        • Kraus V.B.
        • Sprow K.
        • Powell K.E.
        • Buchner D.
        • Bloodgood B.
        • Piercy K.
        • et al.
        Effects of physical activity in knee and hip osteoarthritis: a systematic umbrella review.
        Med Sci Sports Exerc. 2019; 51: 1324-1339https://doi.org/10.1249/MSS.0000000000001944
        • Helmark I.C.
        • Mikkelsen U.R.
        • Børglum J.
        • Rothe A.
        • Petersen M.C.H.
        • Andersen O.
        • et al.
        Exercise increases interleukin-10 levels both intraarticularly and peri-synovially in patients with knee osteoarthritis: a randomized controlled trial.
        Arthritis Res Ther. 2010; 12: R126https://doi.org/10.1186/ar3064
        • Jørgensen A.E.M.
        • Agergaard J.
        • Schjerling P.
        • Heinemeier K.M.
        • van Hall G.
        • Kjaer M.
        The regional turnover of cartilage collagen matrix in late-stage human knee osteoarthritis.
        Osteoarthr Cartil. 2022; 30: 886-895https://doi.org/10.1016/j.joca.2022.03.007
        • Altman R.
        • Asch E.
        • Bloch D.
        • Bole G.
        • Borenstein D.
        • Brandt K.
        • et al.
        Development of criteria for the classification and reporting of osteoarthritis: classification of osteoarthritis of the knee.
        Arthr Rheum. 1986; 29: 1039-1049https://doi.org/10.1002/art.1780290816
        • Brzycki M.
        Strength testing – predicting a one-rep max from reps-to-fatigue.
        J Phys Educ Recreat Dance. 1993; 64: 88-90https://doi.org/10.1080/07303084.1993.10606684
        • Saghaei M.
        • Saghaei S.
        Implementation of an open-source customizable minimization program for allocation of patients to parallel groups in clinical trials.
        J Biomed Sci Eng. 2011; 4: 734-739https://doi.org/10.4236/jbise.2011.411090
        • Tellinghuisen J.
        Estimating real-time qPCR amplification efficiency from single-reaction data.
        Life. 2021; 11https://doi.org/10.3390/LIFE11070693
        • Scholtes S.
        • Krämer E.
        • Weisser M.
        • Roth W.
        • Luginbühl R.
        • Grossner T.
        • et al.
        Global chondrocyte gene expression after a single anabolic loading period: time evolution and re-inducibility of mechano-responses.
        J Cell Physiol. 2018; 233: 699-711https://doi.org/10.1002/jcp.25933
        • Roberts H.M.
        • Law R.J.
        • Thom J.M.
        The time course and mechanisms of change in biomarkers of joint metabolism in response to acute exercise and chronic training in physiologic and pathological conditions.
        Eur J Appl Physiol. 2019; 119: 2401-2420https://doi.org/10.1007/s00421-019-04232-4
        • Mazor M.
        • Best T.M.
        • Cesaro A.
        • Lespessailles E.
        • Toumi H.
        Osteoarthritis biomarker responses and cartilage adaptation to exercise: a review of animal and human models.
        Scand J Med Sci Sports. 2019; 2913435https://doi.org/10.1111/sms.13435
        • Tseng S.
        • Reddi A.H.
        • Di Cesare P.E.
        Cartilage oligomeric matrix protein (COMP): a biomarker of arthritis.
        Biomark Insights. 2009; 4: 33-44https://doi.org/10.4137/bmi.s645
        • Helmark I.C.
        • Petersen M.C.H.
        • Christensen H.E.
        • Kjaer M.
        • Langberg H.
        Moderate loading of the human osteoarthritic knee joint leads to lowering of intraarticular cartilage oligomeric matrix protein.
        Rheumatol Int. 2012; 32: 1009-1014https://doi.org/10.1007/s00296-010-1716-7
        • Millward-Sadler S.J.
        • Wright M.O.
        • Davies L.W.
        • Nuki G.
        • Salter D.M.
        Mechanotransduction via integrins and interleukin-4 results in altered aggrecan and matrix metalloproteinase 3 gene expression in normal, but not osteoarthritic, human articular chondrocytes.
        Arthr Rheum. 2000; 43: 2091-2099https://doi.org/10.1002/1529-0131(200009)43:9<2091::AID-ANR21>3.0.CO;2-C
        • Multanen J.
        • Nieminen M.T.
        • Häkkinen A.
        • Kujala U.M.
        • Jämsä T.
        • Kautiainen H.
        • et al.
        Effects of high-impact training on bone and articular cartilage: 12-month randomized controlled quantitative MRI study.
        J Bone Miner Res. 2014; 29: 192-201https://doi.org/10.1002/jbmr.2015
        • Munukka M.
        • Waller B.
        • Rantalainen T.
        • Häkkinen A.
        • Nieminen M.T.
        • Lammentausta E.
        • et al.
        Efficacy of progressive aquatic resistance training for tibiofemoral cartilage in postmenopausal women with mild knee osteoarthritis: a randomised controlled trial.
        Osteoarthr Cartil. 2016; 24: 1708-1717https://doi.org/10.1016/j.joca.2016.05.007
        • Multanen J.
        • Rantalainen T.
        • Kautiainen H.
        • Ahola R.
        • Jämsä T.
        • Nieminen M.T.
        • et al.
        Effect of progressive high-impact exercise on femoral neck structural strength in postmenopausal women with mild knee osteoarthritis: a 12-month RCT.
        Osteoporos Int. 2017; 28: 1323-1333https://doi.org/10.1007/s00198-016-3875-1
        • Waller K.A.
        • Zhang L.X.
        • Elsaid K.A.
        • Fleming B.C.
        • Warman M.L.
        • Jay G.D.
        Role of lubricin and boundary lubrication in the prevention of chondrocyte apoptosis.
        Proc Natl Acad Sci USA. 2013; 110: 5852-5857https://doi.org/10.1073/pnas.1219289110
        • Iqbal S.M.
        • Leonard C.
        • C Regmi S.
        • De Rantere D.
        • Tailor P.
        • Ren G.
        • et al.
        Lubricin/proteoglycan 4 binds to and regulates the activity of toll-like receptors in vitro.
        Sci Rep. 2016; 618910https://doi.org/10.1038/srep18910
        • Ni G.X.
        • Lei L.
        • Zhou Y.Z.
        • Su J.
        • Schumacher B.
        • Lindley K.
        • et al.
        Intensity-dependent effect of treadmill running on lubricin metabolism of rat articular cartilage.
        Arthr Res Ther. 2012; 14: R256https://doi.org/10.1186/ar4101
        • Ogawa H.
        • Kozhemyakina E.
        • Hung H.H.
        • Grodzinsky A.J.
        • Lassar A.B.
        Mechanical motion promotes expression of Prg4 in articular cartilage via multiple CREB-dependent, fluid flow shear stress-induced signaling pathways.
        Genes Develop. 2014; 28: 127-139https://doi.org/10.1101/gad.231969.113
        • Elsaid K.A.
        • Zhang L.
        • Waller K.
        • Tofte J.
        • Teeple E.
        • Fleming B.C.
        • et al.
        The impact of forced joint exercise on lubricin biosynthesis from articular cartilage following ACL transection and intra-articular lubricin's effect in exercised joints following ACL transection.
        Osteoarthr Cartil. 2012; 20: 940-948https://doi.org/10.1016/j.joca.2012.04.021
        • Jamal J.
        • Roebuck M.M.
        • Lee S.Y.
        • Frostick S.P.
        • Abbas A.A.
        • Merican A.M.
        • et al.
        Modulation of the mechanical responses of synovial fibroblasts by osteoarthritis-associated inflammatory stressors.
        Int J Biochem Cell Biol. 2020; 126105800https://doi.org/10.1016/j.biocel.2020.105800
        • Dunn S.L.
        • Soul J.
        • Anand S.
        • Schwartz J.M.
        • Boot-Handford R.P.
        • Hardingham T.E.
        Gene expression changes in damaged osteoarthritic cartilage identify a signature of non-chondrogenic and mechanical responses.
        Osteoarthr Cartil. 2016; 24: 1431-1440https://doi.org/10.1016/j.joca.2016.03.007
        • Tsuchida A.I.
        • Beekhuizen M.
        • t Hart M.C.
        • Radstake T.R.D.J.
        • Dhert W.J.A.
        • Saris D.B.F.
        • et al.
        Cytokine profiles in the joint depend on pathology, but are different between synovial fluid, cartilage tissue and cultured chondrocytes.
        Arthr Res Ther. 2014; 16: 441https://doi.org/10.1186/s13075-014-0441-0
        • Salter D.M.
        • Millward-Sadler S.J.
        • Nuki G.
        • Wright M.O.
        Differential responses of chondrocytes from normal and osteoarthritic human articular cartilage to mechanical stimulation.
        Biorheology. 2002; 39: 97-108
        • Leong D.J.
        • Li Y.H.
        • Gu X.I.
        • Sun L.
        • Zhou Z.
        • Nasser P.
        • et al.
        Physiological loading of joints prevents cartilage degradation through CITED2.
        FASEB J. 2011; 25: 182-191https://doi.org/10.1096/fj.10-164277
        • He Z.
        • Leong D.J.
        • Xu L.
        • Hardin J.A.
        • Majeska R.J.
        • Schaffler M.B.
        • et al.
        CITED2 mediates the cross-talk between mechanical loading and IL-4 to promote chondroprotection.
        Ann N Y Acad Sci. 2019; 1442: 128-137https://doi.org/10.1111/nyas.14021
        • Bau B.
        • Gebhard P.M.
        • Haag J.
        • Knorr T.
        • Bartnik E.
        • Aigner T.
        Relative messenger RNA expression profiling of collagenases and aggrecanases in human articular chondrocytes in vivo and in vitro.
        Arthr Rheum. 2002; 46: 2648-2657https://doi.org/10.1002/art.10531
        • Snelling S.
        • Rout R.
        • Davidson R.
        • Clark I.
        • Carr A.
        • Hulley P.A.
        • et al.
        A gene expression study of normal and damaged cartilage in anteromedial gonarthrosis, a phenotype of osteoarthritis.
        Osteoarthr Cartil. 2014; 22: 334-343https://doi.org/10.1016/j.joca.2013.12.009
        • Jessberger S.
        • Högger P.
        • Genest F.
        • Salter D.M.
        • Seefried L.
        Cellular pharmacodynamic effects of Pycnogenol® in patients with severe osteoarthritis: a randomized controlled pilot study.
        BMC Compl Altern Med. 2017; 17: 537https://doi.org/10.1186/s12906-017-2044-1
        • Skoffer B.
        • Dalgas U.
        • Maribo T.
        • Søballe K.
        • Mechlenburg I.
        No exacerbation of knee joint pain and effusion following preoperative progressive resistance training in patients scheduled for total knee arthroplasty: secondary analyses from a randomized controlled trial.
        PM&R. 2018; 10: 687-692https://doi.org/10.1016/j.pmrj.2017.11.002
        • Sato T.
        • Konomi K.
        • Yamasaki S.
        • Aratani S.
        • Tsuchimochi K.
        • Yokouchi M.
        • et al.
        Comparative analysis of gene expression profiles in intact and damaged regions of human osteoarthritic cartilage.
        Arthr Rheum. 2006; 54: 808-817https://doi.org/10.1002/art.21638
        • Yagi R.
        • McBurney D.
        • Laverty D.
        • Weiner S.
        • Horton W.E.
        Intrajoint comparisons of gene expression patterns in human osteoarthritis suggest a change in chondrocyte phenotype.
        J Orthop Res. 2005; 23: 1128-1138https://doi.org/10.1016/j.orthres.2004.12.016
        • Brew C.J.
        • Clegg P.D.
        • Boot-Handford R.P.
        • Andrew J.G.
        • Hardingham T.
        Gene expression in human chondrocytes in late osteoarthritis is changed in both fibrillated and intact cartilage without evidence of generalised chondrocyte hypertrophy.
        Ann Rheum Dis. 2010; 69: 234-240https://doi.org/10.1136/ard.2008.097139
        • Fukui N.
        • Ikeda Y.
        • Ohnuki T.
        • Tanaka N.
        • Hikita A.
        • Mitomi H.
        • et al.
        Regional differences in chondrocyte metabolism in osteoarthritis: a detailed analysis by laser capture microdissection.
        Arthr Rheum. 2008; 58: 154-163https://doi.org/10.1002/art.23175
        • van der Kraan P.M.
        • van den Berg W.B.
        Osteophytes: relevance and biology.
        Osteoarthr Cartil. 2007; 15: 237-244https://doi.org/10.1016/j.joca.2006.11.006
        • Uchino M.
        • Izumi T.
        • Tominaga T.
        • Wakita R.
        • Minehara H.
        • Sekiguchi M.
        • et al.
        Growth factor expression in the osteophytes of the human femoral head in osteoarthritis.
        Clin Orthop Relat Res. 2000; 377: 119-125https://doi.org/10.1097/00003086-200008000-00017
        • Middleton J.
        • Arnott N.
        • Walsh S.
        • Beresford J.
        Osteoblasts and osteoclasts in adult human osteophyte tissue express the mRNAs for insulin-like growth factors I and II and the type 1 IGF receptor.
        Bone. 1995; 16: 287-293https://doi.org/10.1016/8756-3282(94)00040-9
        • Jørgensen A.E.M.
        • Kjær M.
        • Heinemeier K.M.
        The effect of aging and mechanical loading on the metabolism of articular cartilage.
        J Rheumatol. 2017; 44: 410-417https://doi.org/10.3899/jrheum.160226
        • Vynios D.H.
        Metabolism of cartilage proteoglycans in health and disease.
        BioMed Res Int. 2014; 2014452315https://doi.org/10.1155/2014/452315
        • Dudhia J.
        • Davidson C.M.
        • Wells T.M.
        • Vynios D.H.
        • Hardingham T.E.
        • Bayliss M.T.
        Age-related changes in the content of the C-terminal region of aggrecan in human articular cartilage.
        Biochem J. 1996; 313: 933-940
        • Struglics A.
        • Hansson M.
        MMP proteolysis of the human extracellular matrix protein aggrecan is mainly a process of normal turnover.
        Biochem J. 2012; 446: 213-223https://doi.org/10.1042/BJ20120274