Tissue catabolism and donor-specific dexamethasone response in a human osteochondral model of post-traumatic osteoarthritis
Arthritis Research & Therapy volume 24, Article number: 137 (2022)
Post-traumatic osteoarthritis (PTOA) does not currently have clinical prognostic biomarkers or disease-modifying drugs, though promising candidates such as dexamethasone (Dex) exist. Many challenges in studying and treating this disease stem from tissue interactions that complicate understanding of drug effects. We present an ex vivo human osteochondral model of PTOA to investigate disease effects on cartilage and bone homeostasis and discover biomarkers for disease progression and drug efficacy.
Human osteochondral explants were harvested from normal (Collins grade 0–1) ankle talocrural joints of human donors (2 female, 5 male, ages 23–70). After pre-equilibration, osteochondral explants were treated with a single-impact mechanical injury and TNF-α, IL-6, and sIL-6R ± 100 nM Dex for 21 days and media collected every 2–3 days. Chondrocyte viability, tissue DNA content, and glycosaminoglycan (sGAG) percent loss to the media were assayed and compared to untreated controls using a linear mixed effects model. Mass spectrometry analysis was performed for both cartilage tissue and pooled culture medium, and the statistical significance of protein abundance changes was determined with the R package limma and empirical Bayes statistics. Partial least squares regression analyses of sGAG loss and Dex attenuation of sGAG loss against proteomic data were performed.
Injury and cytokine treatment caused an increase in the release of matrix components, proteases, pro-inflammatory factors, and intracellular proteins, while tissue lost intracellular metabolic proteins, which was mitigated with the addition of Dex. Dex maintained chondrocyte viability and reduced sGAG loss caused by injury and cytokine treatment by 2/3 overall, with donor-specific differences in the sGAG attenuation effect. Biomarkers of bone metabolism had mixed effects, and collagen II synthesis was suppressed with both disease and Dex treatment by 2- to 5-fold. Semitryptic peptides associated with increased sGAG loss were identified. Pro-inflammatory humoral proteins and apolipoproteins were associated with lower Dex responses.
Catabolic effects on cartilage tissue caused by injury and cytokine treatment were reduced with the addition of Dex in this osteochondral PTOA model. This study presents potential peptide biomarkers of early PTOA progression and Dex efficacy that can help identify and treat patients at risk of PTOA.
Though millions of patients worldwide suffer from osteoarthritis (OA), no disease-modifying OA drug (DMOAD) has been approved due to many challenges in the drug development pipeline, and promising candidates often fail at the level of clinical trials [1,2,3]. Post-traumatic OA (PTOA) is an important disease target since the time of the disease onset (the injury) is known. Many in vitro models of PTOA using human cartilage and isolated chondrocytes offer insights into disease progression, but do not capture the complexity of a full joint that has vasculature, bone, and synovial cross-talk . Models of OA progression using osteochondral explants exist, but have used either diseased knee arthroplasty discards or stem cell-derived engineered organoids that do not completely recapitulate an in vivo healthy control state [5,6,7,8]. Therefore, it is desirable to use a model starting with healthy primary human cartilage and bone to better understand disease progression and drug effects on healthy as well as diseased joint tissues.
Outcome measures of drug efficacy often focus on late-stage pain or structural changes, while early tissue breakdown events can progress years before macroscopic changes are observed. Molecular biomarkers of disease progression can offer insight into early disease processes that could help in drug development as therapeutic endpoints, as well as potential prognostic measurements of patients for the earliest stages of PTOA. In a previous model of PTOA using human knee cartilage explants, mass spectrometry analysis of culture media identified the time-dependent release of proteases and catabolic signaling processes within days after mechanical injury and inflammatory cytokine exposure, hypothesized to be potential biomarkers of disease progression . However, this in vitro model did not incorporate tissue crosstalk that could affect signaling processes and only focused on changes in the media proteome and not the cartilage tissue. Another set of potential biomarker targets is protein fragments produced by enzymatic activity not present in healthy cartilage. A neoepitope of cartilage oligomeric protein (COMP) produced by proteolytic activity was discovered in the synovial fluid of patients and later validated in a human knee cartilage explant model of PTOA, and there is potential to identify fragments of other proteins that are generated under disease stress [10, 11].
In the search for DMOADs for patients at risk of PTOA development, corticosteroids have been identified as promising therapeutics. One of the most well-studied is dexamethasone (Dex), the most potent of the corticosteroid family, currently used for its analgesic effects . In human and animal cartilage tissue models of PTOA, Dex exerted anti-catabolic effects on matrix breakdown and protease production [9, 13,14,15]. However, there are concerns over the safety of Dex use, as in some models using isolated chondrocytes, higher doses of Dex have been shown to cause chondrocyte death [16, 17]. In models using cartilage explants, negative effects on cartilage viability are not commonly observed; however, Dex can exacerbate changes in immune signaling late into disease treatment [9, 14]. Long-term systemic corticosteroid treatment has also been shown to increase incidence of osteoporosis, leading to concerns over effects on bone homeostasis . It is of interest to identify which patients might be the most or least responsive to Dex to target treatment, as there is variability in the response to both disease and Dex at the cellular level [19, 20].
In the present study, we use a human ankle osteochondral model of PTOA, starting with initially normal donor joints, to study early disease progression and to better understand catabolic and anabolic responses in the context of both cartilage and bone. Ankle OA is estimated to have an ~8% prevalence, rarer than hip or knee OA but non-negligible, and up to 90% of the incidence of ankle OA is post-traumatic [21, 22]. This discovery-mode analysis of both the cartilage and media proteomes allows identification of potential biomarkers of PTOA progression and proteolytic events and analysis of donor-specific responses to Dex. We hypothesize that, in this model of PTOA that incorporates cartilage-bone crosstalk, (1) the cartilage and bone tissues will undergo catabolic degradation, (2) anabolic and homeostatic cellular processes will differ from control, and (3) donors will respond differently to disease and drug treatment due to donor-specific differences in cartilage and bone biology.
Materials and methods
Explant harvest and treatment
Human osteochondral explants (3.5 mm diameter, full-thickness cartilage, and ~4 mm bone) were harvested from ankle talocrural joints of seven human donors (62F, 66M, 66M, 44F, 23M, 39M, 70M, Collins grade 0–1) obtained postmortem within 48 h of death through the Gift of Hope Organ and Tissue Donor Network (Itasca, IL). Explants were pre-equilibrated for two days in high glucose phenol red-free Dulbecco’s Modified Eagle Medium (DMEM) (Thermo Fisher) before switching to low-glucose phenol red-free DMEM, supplemented as described . After pre-equilibration, osteochondral explants were treated for 21 days ± a single-impact mechanical injury (60% final strain at 300%/s strain rate; both followed by immediate release at the same rate [11, 23]) and inflammatory cytokines (25 ng/mL tumor necrosis factor alpha (TNF-α) + 50 ng/mL interleukin-6 (IL-6) + 250 ng/mL soluble IL-6 receptor (sIL-6R); treatment IC), as well as with 100 nM Dex alone (D) or treatment IC plus 100 nM Dex (ICD) (Fig. 1A). The use of these cytokines to simulate early inflammation in PTOA has been well-established in previous studies [13, 24]. All donors provided both left and right ankles, and explants from one ankle were used for proteomic analysis while the other ankle was used for sGAG and DNA biochemical analysis. Viability analysis was performed on samples randomized from both ankles. Culture medium was collected and stored at −20°C until analysis.
Biochemical and viability analysis
Cartilage tissue was removed from the underlying bone and digested in 1 mg/mL proteinase K (Sigma). Percent loss of sGAG content from digested cartilage tissue into the medium was determined using the dimethylmethylene blue (DMMB) assay , and tissue DNA content was quantified with the Quanti-iT PicoGreen dsDNA kit (Thermo) according to manufacturer instructions. Significance for biochemical measurements was determined fitting the data to a linear mixed effects model with donor as a random effect followed by a least squares means test, using the R package lmerTest. To determine chondrocyte viability within the cartilage, 100–200 μm vertical slices were cut from intact cartilage and stained with fluorescein diacetate and propidium iodide (Sigma) as previously described .
Mass spectrometry preparation and identification
Culture medium (50 μL) was prepared for mass spectrometry (MS) analysis as described [11, 15]. Cartilage tissue samples were removed from the underlying bone and prepared for MS analysis as described . Discovery MS was performed using a quadrupole Orbitrap benchtop mass spectrometer (Q-Exactive HFX, Thermo Scientific) with prior separation of peptides using a liquid chromatography system (EASY-nLC 1000, Thermo Scientific) on an analytical column (PepMap RSLC C18, 75μm × 25cm, Thermo Scientific) coupled on-line using a nano-electrospray ion source with a column temperature at +45°C (EASY-Spray, Thermo Scientific) using a flow rate of 300nL/min. Protein identification was performed in Proteome Discoverer 2.5 (Thermo Scientific) using two search engines in parallel: a tryptic search against the UniProt human (UP000005640 from 2021-01) sequence database combined with an MSPep spectral search against the NIST_human_Orbitrap_HCD_20160923 library (mass tolerance: 10 and 20 ppm in MS1, MS2 respectively). Other Sequest search settings were modifications: carbamidomethylation (fixed: C), oxidation (variable: M, P) missed cleavages (max 2), mass tolerance (MS1-10ppm, MS2-0.02Da). Label-free protein abundance quantification was obtained by averaging peak area intensities from the top three unique peptides for each protein. To determine individual peptide abundances, we performed a semi-tryptic database search to enable identification of non-tryptic cleavages within the dataset. This was performed using the same combined searches as above but in series. The protein false discovery rate (FDR) was 0.01 for both searches.
For both the peptide and protein MS data sets, proteins were filtered out if they were exogenous or not identified and quantified in at least 70% of samples, and missing values were imputed using the k-nearest neighbor method [9, 27]. Cartilage tissue data was normalized to the DNA content per wet weight of cartilage tissue to adjust for different cell densities between donors (Supplemental Figure 1A). The media data had significant batch effects due to two collections of donors (Supplemental Figure 1B), so data used for principle component and regression analyses were batch corrected using the R limma package function “removeBatchEffect” . Protein and peptide abundance data were log2-transformed and scaled, and principle component analysis (PCA) was performed using the “prcomp” function . Pairwise comparisons between treatments were performed on the individual peptide and protein abundances. Statistical analysis on proteomic and peptide data was performed using limma and MATLAB (MathWorks). Protein and peptide abundances were regressed using partial least squares regression against total sGAG loss or the percent of sGAG loss that was attenuated by Dex treatment, which yielded a dot product of the first two loading vectors for each protein or peptide. The proteins and peptides were ranked by their dot product and analyzed using Gene Set Enrichment Analysis (GSEA) c5 gene sets with the Human UniProt IDs chip and 1,000 repeats for enrichment score distributions [29, 30]. Enrichment analysis for biological processes was performed using the PANTHER database and STRING analysis as previously described [15, 31].
Cartilage viability and sGAG loss during disease and drug treatment
Osteochondral plugs from seven pairs of Collins grade 0-1 human donor ankles were successfully cultured for 3 weeks. Without any treatment, the cartilage tissue lost around 20% sGAG content and maintained chondrocyte viability, consistent with other studies using isolated cartilage (Fig. 1B, Supplemental Figure 2). Low-dose treatment of only Dex had no effect on cell viability or sGAG loss, while treatment with injury and cytokines (IC) caused widespread cell death in the cartilage and increased sGAG loss significantly, up to 54%. The addition of Dex to IC treatment (ICD) ameliorated sGAG loss by 2/3, though not back to control levels, and reduced the amount of chondrocyte death. There were donor-specific differences in how much sGAG loss each individual donor experienced, as well as the degree to which Dex attenuated sGAG loss, if at all (Supplemental Figure 3).
Mass spectrometry identification of proteins in media and cartilage tissue and disease and Dex effects on media and cartilage proteomes
MS analysis identified 18,913 peptides with 2041 identified and quantified proteins in the media, and 8732 peptides corresponding to 1389 proteins in the cartilage tissue. The raw data are available via ProteomeXchange with identifier PXD032213 . After filtering as described in the methods, the media contained 8718 peptides corresponding to 1451 proteins, and the cartilage tissue 7082 peptides and 1141 proteins. Filtered imputed data for individual peptides and proteins can be found in Supplemental File 1. The biological processes with the greatest representation in proteins identified in the media were cellular process, metabolic process, and biological regulation (Supplemental Figure 4A), with the same processes represented in the tissue proteome as well (Supplemental Figure 4B). PCA revealed that injury and cytokine treatment had the strongest contribution to sample variability in the media, while the donor effect was more significant for tissue samples (Fig. 2).
IC treatment caused an increase in the release of extracellular matrix (ECM) components (e.g., COMP, collagens III, V, VI, lumican), proteases (e.g., MMPs −1, −2, −3, −9, −10, −13, −14), pro-inflammatory factors (e.g., lipopolysaccharide-binding protein, serum amyloid A1, complement subcomponents 1R and 1S), and intracellular proteins (e.g., many ribosomal proteins, fructose-bisphosphate aldolase A, heat shock proteins) to the media (Table 1, Supplemental File 2). The major effect on the tissue proteome was a loss of intracellular metabolic proteins, and an increase in MMP2, MMP3, SOD2, and collagen II. Adding Dex to IC-treated osteochondral plugs reduced the media release of ECM proteins including aggrecan, collagen VII, lumican, and fibrillin-2, as well as proteases including MMPs −1, −3, −9, and −10. Many pro-inflammatory factors including interleukin-11, CXCL1, and complement C2, had increased levels with ICD treatment compared to IC, while in the tissue there was no effect on protein abundances. Dex alone reduced the release of some collagens (I, II, IX, XI, XII) and MMPs (−1, −13, −14) into the media, with mixed effects on protease inhibitors. Dex had little effect on the tissue proteome normalized to DNA content (Table 2), with a notable increase in tissue levels of matrix metalloproteinase 3 (MMP3) and superoxide dismutase 2 (SOD2), and a decrease in tissue levels of procollagen-lysine, 2-oxoglutarate 5-dioxygenase 1 and 2 (PLOD1 and 2), and collagen XI. Proteins and peptides found in only one or two treatment conditions, classified as those with quantified abundances in at least six of the seven donors for media samples or all four donors for tissue samples are listed in Supplemental File 3. Statistical results for both DNA normalized and non-DNA normalized cartilage tissue data are included in Supplemental Files 4 and 5, respectively.
Disease and Dex effects on cartilage tissue and bone homeostasis
To investigate proteins identified in both the tissue and media proteomes and their behavior in each compartment under disease stress, we compared the proteins with a significant effect of IC treatment in the tissue to their corresponding changes in the media with IC treatment (Fig. 3A). Some proteins had a decrease in both media and tissue, including PLOD1, collagen IX, and collagen XI. The majority had a decrease in the tissue and a corresponding increase in the media, which had many intracellular metabolic and homeostatic proteins such as ribosomal proteins, alcohol dehydrogenase, and protein disulfide-isomerase A6. Proteins with increased levels in both tissue and the media included MMPs −1, −2, and −3, serpin family E member 1 and member 2 (SERPINE1 and -2), and SOD2.
Media levels of osteomodulin (OMD), osteopontin (SPP1), osteonectin (SPARC), sclerostin (SOST), and alkaline phosphatase (ALPL) were decreased with IC treatment up to 16-fold (Fig. 3B). With ICD treatment, media levels of SPARC and ALPL were increased compared to IC treatment, but SPP1 was further decreased. Dex alone only caused the increase of OMD abundance in the media. Osteoprotegerin (TNFRSF11B) did not experience a change with Dex or IC treatment alone and was slightly decreased by Dex addition to IC treatment compared to IC only.
Analysis of matrix breakdown at the peptide level
Statistical results for individual peptides can be found in Supplemental Files 6 and 7 for media and non-DNA normalized cartilage tissue, respectively. Aggrecan (ACAN, Fig. 4A) had peptides identified in the G1, G2, and G3 globular domains that were mostly increased with IC treatment compared to control in the media. The tryptic peptides identified in the media for collagen II (COL2A1, Fig. 4B) were decreased with IC treatment in the N- and C- terminal regions 2- to 5-fold but increased in the center of the protein. In the tissue, all tryptic peptides identified in the G3 domain of aggrecan had decreased release compared to control under injury and disease stress (Fig. 4C), and the effect on collagen II was only significant in the C-terminal region (Fig. 4D). Semitryptic peptides found only in the IC-treated condition were found from proteins including fibronectin-1 (FN1), COMP, biglycan (BGN), and collagen VI in the tissue, and collagen II, MMP3, MMP13, fibrillin-1 (FBN1), and FN1 in the media (Supplemental File 3).
Dex had mixed effects on semitryptic peptides, those generated by endogenous proteases and not through trypsin digestion during MS preparation. Many proteins with several identified semitryptic peptides had opposite effects, such as aggrecan (Fig. 4E), where semitryptic peptides found in the G1 domain had increased release into the media with Dex treatment, but semitryptic peptides from the G2 and G3 domains were decreased compared to control. Dex broadly caused a decrease in media levels of peptides from collagen II (Fig. 4F), including in the N- and C-terminal regions, with only three having increased abundances compared to control. Collagen I (COL1A1) peptides in both the N- and C-terminal propeptide regions were decreased in the media with both Dex and IC treatment (Supplemental Figure 5).
To investigate which semitryptic peptides may be biomarkers of disease progression, the abundances of each semitryptic peptide were averaged across three biological replicates within each treatment condition and donor, and then regressed against the total percent sGAG loss, an analog for disease severity, for that treatment condition and donor. From the top 50 peptides with the greatest association with increased sGAG loss, the most represented proteins were MMP1, ACAN, lumican (LUM), COMP, and MMP3, with other ECM components fibromodulin (FMOD), and FBN1 notably among the top 50 as well (Fig. 5).
Donor-specific effects on sGAG loss and media proteome
To investigate what protein biomarkers in the media correlated with Dex response, we calculated the percent of the increased sGAG loss caused by IC treatment compared to control that was attenuated by Dex for each donor, and regressed the batch-corrected media proteomic data against this percent reduction. When regressing all samples, apolipoproteins, proteins involved in the complement response, and pro-inflammatory factors were all associated with a lower response to Dex, i.e. less of an attenuation of increased sGAG loss with IC treatment (Fig. 6A). To look for potentially predictive markers of a lack of a Dex response, we regressed only the control samples against the Dex percent reduction (Fig. 6B). As with the regression against all media samples, higher levels of proteins associated with humoral immune responses was associated with a lack of a Dex response, as well as proteasomes 20S and 26S.
In this ex vivo model of human ankle PTOA, we sought to answer how initial stages of cartilage catabolism and anabolism occur in the presence of both cartilage and underlying bone, and what donor-specific responses to both disease and drug treatment would occur. This system used healthy primary osteochondral plugs, rather than total joint replacement discards, as is often the default for in vitro osteochondral studies . IC treatment increased the release of ECM components, proteases, pro-inflammatory factors, and intracellular proteins to the media, with mixed effect on protease inhibitors. This agrees with previous models using isolated human knee cartilage [9, 15]. The addition of Dex to injury and cytokine-treated osteochondral explants reduced the release of many ECM proteins and proteases, consistent with previous conclusions on the anti-catabolic actions of Dex in models of PTOA [9, 12, 14, 15]. Many intracellular proteins identified in both the tissue and media proteomes had decreased levels in tissue and corresponding increases in the media under injury and cytokine stress, and under IC stress, there was a significant amount of chondrocyte death compared to control. As previously hypothesized, the presence of these proteins in the media could be markers of early necrotic cell death [9, 15].
This model allowed for the analysis of specific mechanisms of matrix breakdown by observing the behavior of individual regions of proteins such as aggrecan. Under disease stress, the G2 and G3 globular domains saw significant levels of cleavage and release into the media. The aggrecan found in the tissue was depleted for peptides from the G3 region, which supports previous findings that aggrecan is degraded primarily from the G3 domain in the N-terminal direction towards the other interglobular regions [11, 33, 34]. Decreased abundances of N- and C- terminal domain peptides in the tissue and media were also observed: the N- and C- terminal domains of collagen II are cleaved during maturation and incorporation into the matrix, and lower levels of these peptides in the media and tissue suggest decreases in the synthesis of new collagen II and suppression of anabolic repair . Dex also suppressed the synthesis of new collagen II as indicated by lower levels of the N- and C-terminal peptides in the media. Additionally, Dex affected the activity of proteases and protease inhibitors in healthy tissue, as based on increases in the release of specific semitryptic peptides from aggrecan and decreasing the release of others.
Bone also undergoes changes due to disease and Dex treatment. OMD, SPP1, SPARC, SOST, and ALPL are common biomarkers of bone health and disease responses. All were decreased in the media with IC treatment compared to controls. Some of the selected biomarkers could be protective under osteoarthritic stress: OMD stimulates osteogenesis, and SOST was protective of cartilage viability in a mouse model of OA/PTOA [36,37,38]. The synthesis of collagen I, an indicator of osteoblast metabolism, was decreased under disease stress based on levels of its propeptides in the media . The results presented here contradict previous reports on the presence of these biomarkers under OA stress. SPP1, ALPL, and SPARC have all been shown to be increased in models of OA or associated with OA severity in patient serum [40,41,42]. This could be due to the lack of synovial and other joint tissues or the shorter time span of this model compared to long-term PTOA progression [40, 43]. Corticosteroids are also expected to affect bone biology, as at certain doses and durations are considered detrimental to bone health [18, 44]. However, Dex rescued decreased levels of ALPL and SPARC when added to IC-treated osteochondral explants, but further decreased the levels of SPP1 and collagen I synthesis. Dex alone only caused a change in secreted OMD, increasing the amount of this osteogenic biomarker. The relatively short duration (3 weeks) and low dose of Dex used in our study (100 nM) may mitigate against strong catabolic effects on the bone.
As semitryptic peptides are generated by endogenous protease activity and not by trypsin treatment during MS preparation, these peptides can serve as potential biomarkers for the identification of PTOA progression . Proteases including MMPs are synthesized and released at the earliest stages of PTOA progression, so semitryptic peptides are promising candidates for prognostic PTOA biomarkers . Regression of semitryptic peptides identified in this study against sGAG loss identified semitryptic peptides from aggrecan, COMP, FBN1, FMOD, LUM, MMP1, and MMP3 as those most predictive of disease severity. These are promising candidates for PTOA biomarkers that should be investigated further in clinical samples to determine their potential prognostic use.
Biomarkers are also desirable to determine which patients might respond best to a therapeutic such as Dex [45, 46]. Proteomic data were regressed against the amount of sGAG lost under disease stress that was attenuated by Dex treatment for each donor in order to determine which proteins best predicted anti-catabolic activity. Humoral immune proteins and the pro-inflammatory apolipoproteins A, B, C, and H were all inversely associated with a Dex attenuation of sGAG loss, while the anti-inflammatory apolipoprotein E was associated with a Dex attenuation effect [47,48,49]. However, apolipoproteins A and B have also been shown to have anti-inflammatory activity, inhibiting T cell activation and the innate immune response . When regressing only the untreated control condition, most analogous to healthy patients or very early stages of the disease, a lack of a Dex attenuation response was also associated with proteasomes 20S and 26S. Proteasomes play a role in T cell activation and immune responses, and circulating proteasomes have been associated with immunological activity and cellular damage [51, 52]. These results suggest that a highly elevated initial inflammatory state could prevent the anti-catabolic activity of Dex, which may help to distinguish patients that may best respond to Dex. Variability in donor responses to disease as well as the baseline inflammatory state of a joint are likely tied to sex, diet, BMI, lifestyle, genetic disposition, cartilage thickness, and many other factors that cannot be explored with such a small number of donors . Further investigation of clinical samples from many Dex-treated patients is necessary to follow up on these results.
PTOA of the talocrural joint is a relevant model system to study due to the high incidence of ankle PTOA, but these results must be validated in samples from knee and hip cartilage to confirm their translatability to those joints that may have different treatment and disease responses [22, 54]. Cartilage thickness of each osteochondral plug varied within and between donors, though normalizing to cartilage weight and DNA content addressed differences in both the amount of tissue and number of cells within the tissue. Osteochondral plugs were randomized between treatment conditions to further address location-specific variation. Osteochondral plugs had differences in bone density, vascularity, and cellularity, and while the bone could be cut to approximately the same length, it did not always cut perfectly. Our present results confirm several of those previous findings of disease and Dex effects , but notably used a single concentration of and continuous exposure to Dex, and these results must be explored in further experiments that compare different timings and doses of Dex treatment. The current approach to intra-articular injection of corticosteroids involves a much higher initial dose that is maintained for a much shorter time due to rapid clearance from the joint. However, novel methods of drug delivery are currently under investigation to replicate low-dose Dex delivery over a longer durations targeted to cartilage [55,56,57]. The addition of synovial tissue to this system would broaden applicability to knee and other joints, but would introduce uncontrolled and inconsistent levels of inflammatory cytokines compared to controlled addition of exogenous cytokines . Our proteomic analysis used MS identification based on trypsin digestion of media and tissue proteins, which removes the ability to identify some endogenously produced peptides. We utilized a high threshold of confidence for identification and strict filtering for inclusion to compensate for any single-peptide identifications or quantifications. The tissue proteome only includes proteins that could be solubilized and extracted from the tissue. Proteins identified in the media by MS could have been released from both cartilage and bone, so it is difficult to explicitly say what tissue changes in media levels of proteins originate from.
This osteochondral human ex vivo model of PTOA allowed for the analysis of disease progression and Dex effects in the context of both cartilage and bone starting with healthy tissue. The combination of a single mechanical injury, TNF-α, and IL-6 treatment caused catabolic changes and a suppression of cartilage anabolism, in line with other studies using isolated cartilage. This model of early PTOA caused some changes to bone metabolism as well, which did not extend to full bone breakdown but suppressed normal homeostatic processes. The addition of Dex to this PTOA model confirmed the anti-catabolic effects seen in previous studies using isolated cartilage, and peptidomic analyses revealed a suppression of collagen II synthesis and changes to endogenous protease activity, even with healthy cartilage. Elevated levels of humoral proteins and apolipoproteins were associated with a lack of a Dex response, a promising set of predictive biomarkers for patient Dex response that could allow for targeted, personalized drug treatment that can be combined with pro-anabolic drugs to offset potential off-target effects . With Dex as a promising DMOAD candidate for early PTOA intervention, prognostic biomarkers are critical to determine patients most at risk of developing PTOA as well as those most likely to respond to Dex. Semitryptic peptide biomarkers associated with PTOA progression from several ECM proteins and proteases presented in this study should be validated against longitudinal clinical samples and used in tandem with other biomarker candidates, as with patient heterogeneity there is likely not only one single sufficient biomarker.
Cartilage oligomeric matrix protein
100 nM dexamethasone treatment
Dulbecco’s modified Eagle’s medium
Disease-modifying osteoarthritis drug
False discovery rate
Injury and cytokine treatment
Injury, cytokine, and 100 nM Dex treatment
Procollagen-lysine, 2-oxoglutarate 5-dioxygenase
Serpin family E member
Soluble IL-6 receptor
Superoxide dismutase 2
Tumor necrosis factor alpha
Oo WM, Little C, Duong V, Hunter DJ. The development of disease-modifying therapies for osteoarthritis (DMOADs): The evidence to date. Drug Des Devel Ther. 2021;15:2921–45.
Ghouri A, Conaghan PG. Update on novel pharmacological therapies for osteoarthritis. Ther Adv Musculoskelet Dis. 2019;11:1759720X1986449.
Vincent TL. Of mice and men: converging on a common molecular understanding of osteoarthritis. Lancet Rheumatol. 2020;2:e633–45.
Jiang A, Xu P, Sun S, Zhao Z, Tan Q, Li W, et al. Cellular alterations and crosstalk in the osteochondral joint in osteoarthritis and promising therapeutic strategies. Connect Tissue Res. 2021;62:709–19.
Geurts J, Jurić D, Müller M, Schären S, Netzer C. Novel ex vivo human osteochondral explant model of knee and spine osteoarthritis enables assessment of inflammatory and drug treatment responses. Int J Mol Sci. 2018;19:1314.
Huntley JS, Bush PG, McBirnie JM, Simpson AH, Hall AC. Chondrocyte death associated with human femoral osteochondral harvest as performed for mosaicplasty. J Bone Jt Surg - Ser A. 2005;87:351–60.
Lin H, Lozito TP, Alexander PG, Gottardi R, Tuan RS. Stem cell-based microphysiological osteochondral system to model tissue response to interleukin-1β. Mol Pharm. 2014;11:2203–12.
Lozito TP, Alexander PG, Lin H, Gottardi R, Cheng AWM, Tuan RS. Three-dimensional osteochondral microtissue to model pathogenesis of osteoarthritis. Stem Cell Res Ther. 2013;4:S6.
Black RM, Wang Y, Struglics A, Lorenzo P, Chubinskaya S, Grodzinsky AJ, et al. Proteomic clustering reveals the kinetics of disease biomarkers in bovine and human models of post-traumatic osteoarthritis. Osteoarthr Cartil Open. 2021;3:100191.
Åhrman E, Lorenzo P, Holmgren K, Grodzinsky AJ, Dahlberg LE, Saxne T, et al. Novel cartilage oligomeric matrix protein (COMP) neoepitopes identified in synovial fluids from patients with joint diseases using affinity chromatography and mass spectrometry. J Biol Chem. 2014;289:20908–16.
Wang Y, Li Y, Khabut A, Chubinskaya S, Grodzinsky AJ, Önnerfjord P. Quantitative proteomics analysis of cartilage response to mechanical injury and cytokine treatment. Matrix Biol. 2017;63:11–22.
Black R, Grodzinsky AJ. Dexamethasone: Chondroprotective corticosteroid or catabolic killer? Eur Cells Mater. 2019;38:246–63.
Lu YCS, Evans CH, Grodzinsky AJ. Effects of short-term glucocorticoid treatment on changes in cartilage matrix degradation and chondrocyte gene expression induced by mechanical injury and inflammatory cytokines. Arthritis Res Ther. 2011;13:R142.
Li Y, Wang Y, Chubinskaya S, Schoeberl B, Florine E, Kopesky P, et al. Effects of insulin-like growth factor-1 and dexamethasone on cytokine-challenged cartilage: Relevance to post-traumatic osteoarthritis. Osteoarthr Cartil. 2015;23:266–74.
Black RM, Wang Y, Struglics A, Lorenzo P, Tillgren V, Rydén M, et al. Proteomic analysis reveals dexamethasone rescues matrix breakdown but not anabolic dysregulation in a cartilage injury model. Osteoarthr Cartil Open. 2020;2:100099.
Zaman F, Chrysis D, Huntjens K, Chagin A, Takigawa M, Fadeel B, et al. Dexamethasone differentially regulates Bcl-2 family proteins in human proliferative chondrocytes: Role of pro-apoptotic Bid. Toxicol Lett. 2014;224:196–200.
Shen C, Cai GQ, Peng JP, Chen XD. Autophagy protects chondrocytes from glucocorticoids-induced apoptosis via ROS/Akt/FOXO3 signaling. Osteoarthr Cartil. 2015;23:2279–87.
Weinstein RS. Glucocorticoid-induced csteoporosis and csteonecrosis. Endocrinol Metab Clin North Am. 2012;41:595–611.
Chan MWY, Gomez-Aristizábal A, Mahomed N, Gandhi R, Viswanathan S. A tool for evaluating novel osteoarthritis therapies using multivariate analyses of human cartilage-synovium explant co-culture. Osteoarthr Cartil. 2022;30:147–59.
Siddappa R, Licht R, van Blitterswijk C, de Boer J. Donor variation and loss of multipotency during in vitro expansion of human mesenchymal stem cells for bone tissue engineering. J Orthop Res. 2007;25:1029–41.
Murray C, Marshall M, Rathod T, Bowen CJ, Menz HB, Roddy E. Population prevalence and distribution of ankle pain and symptomatic radiographic ankle osteoarthritis in community dwelling older adults: a systematic review and cross-sectional study. PLoS One. 2018;13:e0193662.
Delco ML, Kennedy JG, Bonassar LJ, Fortier LA. Post-traumatic osteoarthritis of the ankle: A distinct clinical entity requiring new research approaches. J Orthop Res. 2017;35:440–53.
Li Y, Frank E, Wang Y, Chubinskaya S, Huang H-H, Grodzinsky AJ. Moderate dynamic compression inhibits pro-catabolic response of cartilage to mechancial injury, TNF-alpha and IL-6, but accentuates degradation above a strain threshold. Osteoarthr Cartil. 2013;21:1–17.
Sui Y, Lee JH, DiMicco MA, Vanderploeg EJ, Blake SM, Hung HH, et al. Mechanical injury potentiates proteoglycan catabolism induced by interleukin-6 with soluble interleukin-6 receptor and tumor necrosis factor α in immature bovine and adult human articular cartilage. Arthritis Rheum. 2009;60:2985–96.
Farndale RW, Buttle DJ, Barrett AJ. Improved quantitation and discrimination of sulphated glycosaminoglycans by use of dimethylmethylene blue. Biochim Biophys Acta. 1986;883:173–7.
Hsueh MF, Khabut A, Kjellström S, Önnerfjord P, Kraus VB. Elucidating the molecular composition of cartilage by proteomics. J Proteome Res. 2016;15:374–88.
Lazar C, Gatto L, Ferro M, Bruley C, Burger T. Accounting for the multiple natures of missing values in label-free quantitative proteomics data sets to compare imputation strategies. J Proteome Res. 2016;15:1116–25.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.
Daly MJ, Patterson N, Mesirov JP, Golub TR, Tamayo P, Spiegelman B. PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34:267–73.
Mi H, Ebert D, Muruganujan A, Mills C, Albou LP, Mushayamaha T, et al. PANTHER version 16: A revised family classification, tree-based classification tool, enhancer regions and extensive API. Nucleic Acids Res. 2021;49:D394–403.
Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, et al. The PRIDE database and related tools and resources in 2019: Improving support for quantification data. Nucleic Acids Res. 2019;47:D442–50.
Patwari P, Kurz B, Sandy JD, Grodzinsky AJ. Mannosamine inhibits aggrecanase-mediated changes in the physical properties and biochemical composition of articular cartilage. Arch Biochem Biophys. 2000;374:79–85.
Swärd P, Wang Y, Hansson M, Lohmander LS, Grodzinsky AJ, Struglics A. Coculture of bovine cartilage with synovium and fibrous joint capsule increases aggrecanase and matrix metalloproteinase activity. Arthritis Res Ther. 2017;19:1–12.
Kraus VB, Karsdal MA. Osteoarthritis: Current molecular biomarkers and the way forward. Calcif Tissue Int. 2021;109:329–38.
Zappia J, Van der Cruyssen R, Sanchez C, Lausberg C, Lambert C, Florin A, et al. Osteomodulin impacts positively the bone remodeling process in osteoarthritis. Bone Rep. 2021;14:100904.
Lin W, Zhu X, Gao L, Mao M, Gao D, Huang Z. Osteomodulin positively regulates osteogenesis through interaction with BMP2. Cell Death Dis. 2021;12:147.
Chang JC, Christiansen BA, Murugesh DK, Sebastian A, Hum NR, Collette NM, et al. SOST/Sclerostin improves posttraumatic osteoarthritis and inhibits MMP2/3 expression after injury. J Bone Miner Res. 2018;33:1105–13.
Eastell R, Hannon RA. Biochemical markers of bone turnover. In: Lobo RA, editor. Treatment of the Postmenopausal Woman (Third Edition): Basic and Clinical Aspects: Academic Press; 2007. p. 337–49.
Cheng C, Gao S, Lei G. Association of osteopontin with osteoarthritis. Rheumatol Int. 2014;34:1627–31.
Park HM, Lee JH, Lee YJ. Positive association of serum alkaline phosphatase level with severe knee osteoarthritis: a nationwide population-based study. Diagnostics. 2020;10:1016.
Nakamura S, Kamihagi K, Satakeda H, Katayama M, Pan H, Okamoto H, et al. Enhancement of SPARC (osteonectin) synthesis in arthritic cartilage: Increased levels in synovial fluids from patients with rheumatoid arthritis and regulation by growth factors and cytokines in chondrocyte cultures. Arthritis Rheum. 1996;39:539–51.
Nanba Y, Nishida K, Yoshikawa T, Sato T, Inoue H, Kuboki Y. Expression of osteonectin in articular cartilage of osteoarthritic knees. Acta Med Okayama. 1997;51:239–43.
Liu Y, Chen Y, Zhao H, Zhong L, Wu L, Cui L. Effects of different doses of dexamethasone on bone qualities in rats. J Biomed Eng. 2011;28:737–43.
Huebner KD, Shrive NG, Frank CB. Dexamethasone inhibits inflammation and cartilage damage in a new model of post-traumatic osteoarthritis. J Orthop Res. 2014;32:566–72.
Saklatvala J. Glucocorticoids: Do we know how they work? Arthritis Res. 2002;4:146–50.
De Seny D, Cobraiville G, Charlier E, Neuville S, Lutteri L, Le GC, et al. Apolipoprotein-A1 as a damage-associated molecular patterns protein in osteoarthritis: Ex vivo and in vitro pro-inflammatory properties. PLoS One. 2015;10:1–17.
Lee JY, Kang MJ, Choi JY, Park JS, Park JK, Lee EY, et al. Apolipoprotein B binds to enolase-1 and aggravates inflammation in rheumatoid arthritis. Ann Rheum Dis. 2018;77:1480–9.
Zhang H, Wu LM, Wu J. Cross-talk between apolipoprotein E and cytokines. Mediators Inflamm. 2011;(949072):1–10.
Cho NH, Seong SY. Apolipoproteins inhibit the innate immunity activated by necrotic cells or bacterial endotoxin. Immunology. 2009;128:479–86.
Majetschak M, Perez M, Sorell LT, Lam J, Maldonado ME, Hoffman RW. Circulating 20S proteasome levels in patients with mixed connective tissue disease and systemic lupus erythematosus. Clin Vaccine Immunol. 2008;15:1489–93.
Verbrugge EE, Scheper RJ, Lems WF, de Gruijl TD, Jansen G. Proteasome inhibitors as experimental therapeutics of autoimmune diseases. Arthritis Res Ther. 2015;17:1–10.
Webb RJ, Mazidi M, Lip GYH, Kengne AP, Banach M, Davies IG. The role of adiposity, diet and inflammation on the discordance between LDL-C and apolipoprotein B. Nutr Metab Cardiovasc Dis. 2022;32:605–15.
Hsueh MF, Önnerfjord P, Bolognesi MP, Easley ME, Kraus VB. Analysis of “old” proteins unmasks dynamic gradient of cartilage turnover in human limbs. Sci Adv. 2019;5:1–10.
Krishnan Y, Rees HA, Rossitto CP, Kim SE, Hung HHK, Frank EH, et al. Green fluorescent proteins engineered for cartilage-targeted drug delivery: Insights for transport into highly charged avascular tissues. Biomaterials. 2018;183:218–33.
Bajpayee AG, De la Vega RE, Scheu M, Varady NH, Yannatos IA, Brown LA, et al. Sustained intra-cartilage delivery of low dose dexamethasone using a cationic carrier for treatment of post traumatic osteoarthritis. Eur Cell Mater. 2017;34:341–64.
Bajpayee AG, Quadir MA, Hammond PT, Grodzinsky AJ. Charge based intra-cartilage delivery of single dose dexamethasone using Avidin nano-carriers suppresses cytokine-induced catabolism long term. Osteoarthr Cartil. 2016;49:48–56.
Dwivedi G, Flaman L, Frank E, Geishecker E, Rosen V, Chubinskaya S, et al. Human cartilage-bone-synovium microphysiological system to study PTOA pathogenesis and treatment on earth and in space. Osteoarthr Cartil. 2019;27:S167.
The authors acknowledge Gift of Hope Organ & Tissue Donor Network (Itasca, IL) and donors' families.
Supported by grants NIH-NCATS UG3/UH3 TR002186 (AJG), the 2022 Siebel Scholars award (RMB), the Rush Klaus Kuettner Endowed Chair (SC), the Swedish Rheumatism Association (PÖ), the Alfred Österlund Foundation (PÖ), the Greta & Johan Kock and Crafoord Foundations (PÖ), the Olle Engkvist Foundation (PÖ). Instrument funding (MS) from the IAB Lundberg Foundation (PÖ).
Ethics approval and consent to participate
All procedures were approved by Rush University Medical IRB and the MIT COUHES committee.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Donor differences in DNA content and batch effect in media proteome. Figure S2. Fluorescent imaging assessment of cartilage viability. Figure S3. GAG loss within individual donors. Figure S4. Proteomic identification of media and extracted tissue proteins. Figure S5. Changes in media abundances of collagen I tryptic peptides with Dex and disease treatment.
About this article
Cite this article
Black, R.M., Flaman, L.L., Lindblom, K. et al. Tissue catabolism and donor-specific dexamethasone response in a human osteochondral model of post-traumatic osteoarthritis. Arthritis Res Ther 24, 137 (2022). https://doi.org/10.1186/s13075-022-02828-4
- Post-traumatic osteoarthritis
- Mass spectrometry
- Cartilage matrix