Skip to main content
Download PDF

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

Cell Communication and Signaling volume 21, Article number: 137 (2023) Cite this article

Abstract

Background

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

Methods

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

Results

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

Conclusion

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

Video Abstract

Introduction

Osteoarthritis (OA), the most common joint disease in humans, affects over 500 million people worldwide, including a significant proportion of the elderly population, resulting in pain and disability [1]. The incidence of OA and its socioeconomic costs, representing between 1% and 2.5% of the gross domestic product in developed countries [2], are expected to rise with the increase in global life expectancy. However, there is currently no cure for OA, and management of this debilitating condition only consists of diet and lifestyle changes, pain management, and prosthetic joint replacement at late stages. Thus, the development of disease-modifying osteoarthritis drugs (DMOADs) is a major challenge and requires a better understanding of OA pathophysiology.

Although long considered to be only a cartilage disease, it is now widely recognized that OA affects all joint tissues through inflammatory and degenerative processes [3]. Chronic low-grade inflammation appears to play a key role in OA development, and joint inflammation is largely sustained by the synovium. Indeed, synovitis correlates with the severity of OA [4], and synovium changes can arise before the visible onset of cartilage degradation, thus highlighting the significance of synovial inflammation in OA. There is nonetheless a close relationship between the synovium and cartilage during OA. Chondrocytes, through the production of damage-associated molecular patterns, undergo and support local inflammation of the synovial tissue by producing and releasing pro-inflammatory cytokines. In this chronic pro-inflammatory environment, chondrocytes acquire a pro-catabolic and pro-inflammatory phenotype, leading to the production of a plethora of proteins, such as proteases and pro-inflammatory cytokines [5, 6].

Moreover, IL-1β (IL-1b), TNF, and, to a lesser extent IL-6, the three main pro-inflammatory cytokines found in abundance in OA synovial fluid [7], are commonly used to mimic pro-catabolic and pro-inflammatory chondrocyte phenotypes in vitro. Thus, several mechanisms are altered by these cytokines. They decrease the capacity of chondrocytes to produce extracellular matrix, while broadly enhancing their catabolic activity [8], ultimately leading to extracellular matrix (ECM) degradation. Moreover, pro-inflammatory stimulation of chondrocytes results in IL-1b, TNF, and IL-6 overexpression [5], ultimately promoting deleterious mechanisms involved in OA, such as cellular senescence [9] and oxidative stress [10].

Given the crucial role played by these cytokines in perpetuating the low-grade inflammation characteristic of OA, it appeared relevant to evaluate the effect of anti-TNF (e.g., etanercept, adalimumab), anti-IL-1 (e.g., anakinra, lutikizumab), and anti-IL-6 (tocilizumab) drugs as DMOADs. According to the regulatory requirements, a DMOAD must improve symptoms, such as pain and/or disability, and stop or at least slow the loss of joint space width [11]. However, despite promising results in OA animal models [12,13,14] and efficacy in other types of inflammatory rheumatisms, all of these anti-cytokine drugs failed to relieve pain in patients with OA [15]. Specifically, the recently developed dual anti-IL1α/β molecule lutikizumab was tested in OA patients with synovitis and demonstrated only limited improvement in pain scores and synovial inflammation [16]. Similarly, etanercept [17] and tocilizumab [15] failed to relieve pain in patients with hand OA. Thus, given the therapeutic failure of individual targeting of IL-1b, TNF, or IL-6 signaling, a thorough investigation of shared dysregulated pathways may allow the identification of common clinically relevant therapeutic targets for the development of higher-impact therapies for OA.

Several metabolic dysregulations have been described in an OA context. Among them, mitochondrial defects and alterations in glycolytic enzymes are described in OA chondrocytes [18,19,20]. Interestingly, promising new therapeutic compounds currently under clinical evaluation such as sprifermin (recombinant human FGF18) have recently shown a protective effect against OA-related metabolic alterations in vitro, evidenced by an improvement of mitochondrial function and inhibition of ROS production [21]. It is also increasingly recognized that there is a close link between inflammation and metabolic adaptation [22]. In inflammatory arthritis, deleterious metabolism changes receive increasing attention, particularly in synovial cells [23]. Moreover, impairment of adenosine metabolism by methotrexate, a drug used in rheumatoid arthritis treatment [22] has been shown to significantly improve pain and stiffness in OA patients, while drugs targeting pro-inflammatory cytokines fail to improve it [27], thus confirming that metabolism should be considered a key player in the clinical response to treatment. Nevertheless, the underlying mechanisms of these metabolic dysregulations in the pathophysiology of OA, although present, are still poorly understood and the link between chronic low-grade inflammation and these dysregulations in cartilage remains to be explored.

Here, multi-omics analyses were performed in articular chondrocytes from OA joints (OACs) to decipher common transcriptional and proteomic changes in response to IL-1b, TNF, and IL-6 pro-inflammatory stimuli. To unravel whether the common inflammatory signature identified is specific to OACs, the transcriptome of chondrocytes from non-osteoarthritic joints (NCs) following pro-inflammatory stimuli was used as an internal control and compared to the transcriptome of OACs.

In this study, we highlighted that the metabolic shift usually observed in chondrocytes upon exposure to pro-inflammatory stimuli was a specific feature of OA chondrocytes and was not found in non-OA chondrocytes. This insensitivity of non-OA chondrocytes to metabolic alteration may open new avenues for understanding the dysregulated mechanisms that occur during OA.

Results

Transcriptomic analysis revealed overlapping but non-redundant NC and OAC pro-inflammatory responses

Following RNA-sequencing (Fig. 1), OACs principal component analysis (PCA) showed that the samples treated with either IL-1b or TNF significantly clustered away from the untreated samples (Fig. 2A). By contrast, samples treated with IL-6 overlapped with untreated samples. Interestingly the sample segregation observed on NCs PCA was not as distinct as on OACs PCA, exhibiting only significant segregation of IL-1b-treated samples from untreated ones (Fig. 2B). Altogether these data suggest a more pronounced pro-inflammatory response by OACs compared to NCs.

Fig. 1
figure 1

Overview of the study design. OA chondrocytes (OACs) and non-OA chondrocytes (NCs) were treated for 24 h with or without pro-inflammatory cytokines. A multi-omic (RNA 3’ SRP and protein MS) approach was used in OACs to identify a pro-inflammatory signature. Then, NC transcriptomics data were compared to the OACs pro-inflammatory signature to define disease-specific expression features. Scale bar histology: 500 µm, Scale bar cells culture: 100 µm

Full size image
Fig. 2
figure 2

OAC and NC RNA-seq data. (a, b) Principal component analysis (PCA) was performed on the expression data of OAC (a) and NC samples (b). The first two eigenvalues were plotted with data ellipses for each treatment. (c) Summary of the number of DEGs in each condition and DEGs common between IL-1b and TNF treatments highlighted by Venn diagrams. (d, e, f, g) Volcano plots representing DEGs (p-adjust < 0.10; log2(Fold-Change) >|0.58|) in OACs (d, e) and NCs (f, g) in response to 1 ng/mL IL-1b (d, f) or 25 ng/mL TNF (e, g). (h) Histograms representing normalized MMP1, SOD2, SERPINE2, and Slc39a8 (ZIP8) gene counts in OACs and NCs; **P < 0.01 and ***P < 0.001, versus untreated condition using multiple t-tests with Benjamini & Hochberg correction; N = 5 per group. Values are expressed as means ± SEM

Full size image

As expected from the PCAs, very few differentially expressed genes (DEGs) were identified between IL-6-treated and untreated samples in OACs (56 DEGs), and no DEGs were identified in NCs (data not shown). We identified 3365 DEGs in response to IL-1b and 1607 DEGs in response to TNF in OACs, while only 537 and 499 DEGs were found in response to IL-1b and TNF, respectively, in NCs (Fig. 2C). Indeed, volcano plots indicated that the overall responses of OACs (Fig. 2D-E) to either IL-1b or TNF were greater compared to NCs (Fig. 2F-G). Moreover, OACs exhibited a higher number of DEG in response to IL-1b (3365 DEGs) than to TNF treatment (1607 DEGs). Nevertheless, 78% of the DEGs identified in response to TNF (1239/1607) were common to those identified in response to IL-1b, highlighting strong redundancy of the TNF response with that of IL-1b (Fig. 2C).

Despite the overall differences observed between NC and OAC responses, 408 DEGs in response to IL-1b and 411 DEGs in response to TNF were found to be common in both cell types (Figure S1A, B). Among these common DEGs overexpressed in both cell types in response to IL-1b and TNF, there were several key players of OA, such as MMP1 [24]; SOD2 [25], SERPINE2 [26, 27], and Slc39a8 [28] (Fig. 2H).

Nevertheless, among the 1249 DEGs identified as common to both cytokines in OACs, the DEGs of interest, 1079 (≈86%) were specific to OACs and warranted further examination.

Multi-omics analysis identified a canonical pro-inflammatory signature in OACs

To decipher the specific OAC signature in response to pro-inflammatory stimuli and to complete the description of OACs upon exposure to pro-inflammatory stimuli, in addition to the 3’ SRP analysis, we performed tandem mass spectrometry analysis of the proteins extracted from OAC lysates (5 OA patients) treated the same way as the samples described previously (Fig. 1). This analysis validated and quantified 2688 proteins and, for the entire dataset, we were able to assign 794,053 MS/MS spectra, or an average of 26,468 per sample. The PCA analysis of proteomic samples revealed one significant outlier patient despite the absence of previously obvious differences with other samples (Fig. 3A, in black). Following the removal of this patient from our analysis, the variability between the different samples was drastically decreased, as illustrated by the percentage of variance associated with the first component (PC1), which changed from 67% (N = 5) to 26% (N = 4) (Fig. 3B). Similar to the transcriptomic results, a low number (52) of differentially expressed proteins (DEPs) were identified in samples treated with IL-6 (data not shown). In light of the small number of DEGs and DEPs identified in response to IL-6, our subsequent analyses focused only on samples treated with either IL-1b or TNF. Interestingly, whereas a significant difference in the number of DEGs in response to IL-1b (3365) or TNF (1607) was observed at the transcriptomic level, at the proteomic level a similar number of DEPs was found in response to IL-1b (154 DEP) or TNF (140 DEPs) (Fig. 3C, D). Among these identified DEPs, 89 (Fig. 3C) and 82 (Fig. 3D) DEPs matched with identified DEGs in response to IL-1b or TNF, respectively. The KEGG pathways analysis based on these common DEGs/DEPs revealed key over-represented pathways such as “Antigen processing and presentation” and “NF-kB signaling pathway” that are shared by both the IL-1b and TNF responses (Fig. 3E, F).

Fig. 3
figure 3

Definition of the canonical pro-inflammatory signature in OACs. (a, b) PCAs performed on the protein expression data from (a) 5 OA patients or (b) the 4 selected patients. The first two eigenvalues were plotted with data ellipses for each patient. (c, d) Volcano plots representing DEPs (p-value < 0.05; |log2(Fold-Change)|> 0) and Venn diagram highlighting DEPs and DEGs matched in OAC in response to (c) 1 ng/mL IL-1b, or (d) 25 ng/mL TNF compared to the untreated condition. (e, f) KEGG enrichment analysis based on DEGs/DEPs in response to (e) IL-1b or (f) TNF. (g) Venn diagram highlighting 56 DEGs/DEPs common between IL-1b and TNF responses. (h) KEGG enrichment analysis using the list of 56 matched DEPs/DEGs

Full size image

Thereafter, to define the pro-inflammatory signature in OACs, we mainly focused on the similarities between the IL-1b and TNF responses at both the transcriptomic and proteomic levels. With these criteria, we identified 56 matched DEPs/DEGs found in both the IL-1b and TNF responses, which determined a canonical pro-inflammatory signature in OACs (Fig. 3G, Table 1). To begin the exploration of this pro-inflammatory signature, we performed a KEGG enrichment analysis. Among the KEGG-enriched pathways identified using this list of 56 DEPs/DEGs, some were expected, such as the “NF-κB signaling pathway” [10] and the “HIF-1 signaling pathway” [29], whereas others, such as “Antigen presentation” [30] and “Ferroptosis” [31], were more recently identified in the context of OA inflammation (Fig. 3H).

Table 1 List of the 56 common DEG/DEPS defining the canonical OAC pro-inflammatory signature. Summarized log2(Fold-Change) and associated p-value for the 56 targets that define the OAC pro-inflammatory signature are shown. The first 21 genes, in bold, are also differentially expressed in NCs in response to either IL-1b or TNF
Full size table

OACs exhibit disease-specific metabolic features following pro-inflammatory exposure

Following the description of the canonical pro-inflammatory signature characterizing OACs, we wished to determine whether the same pro-inflammatory signature could also be observed in non-OA chondrocytes. To address this question, the 56 matched DEPs/DEGs defining the canonical pro-inflammatory signature of OACs were searched in the 214 common DEGs identified in NCs in response to IL-1b or TNF. These 214 DEGs comprised 21 of the 56 DEGs/DEPs identified in OACs (Table 1, bold genes) (Fig. 4A). To further characterize the differences in the inflammatory responses of OACs versus NCs, protein–protein networks were constructed, using either the 56 DEGs/DEPs of OACs (Fig. 4B) or the 21 DEGs retrieved in NCs (Fig. 4C). Interestingly, the 56 DEGs/DEPs characterizing the pro-inflammatory signature of OACs revealed three interacting groups of proteins (Fig. 4B): one group was characterized by common pro-inflammatory targets (Blue circle: SOD2, ICAM-1, EGFR, etc.), another was associated with proteins involved in antigen presentation (Yellow circle: HLA-B, B2M, PSME2, etc.), and a third group was composed of metabolism-related proteins (Green box: SLC2A1, HK2, SLC16A3, and GFPT2). A similar analysis using the 21 DEGs characterizing the NCs response to IL-1b and TNF (Fig. 4C) demonstrated that, while the common pro-inflammatory targets (Blue circle) and antigen presentation (Yellow circle) groups were shared between both the NC and OAC pro-inflammatory responses, the group related to metabolism was specific to the OAC response. This metabolism group is composed of four proteins: Glut-1 (Solute Carrier Family 2 Member 1 – Slc2a1), a glucose transporter; MCT-4 (Solute Carrier Family 16 Member 3 – Slc16a3), a lactate/H+ exporter; HK2 (hexokinase 2), the enzyme catalyzing the first step in glycolysis; and GFAT2 (Glutamine-Fructose-6-Phosphate Transaminase 2 – Gfpt2), the rate-limiting enzyme of the hexosamine biosynthesis pathway (HBP). Among these, Gfpt2 expression was also increased in NCs, but the expression of Glut-1, MCT-4, and HK2 was specifically modulated in OACs (Table 1). Using RT-qPCR, we then validated the overexpression of Glut-1, MCT-4, and HK2 in OACs by both cytokine treatments as well as the overexpression of GFAT2 in response to either IL-1b or TNF in OACs and NCs (Fig. 4D).

Fig. 4
figure 4

Identification of a specific feature of OACs. (a) Venn diagram showing the 21 common DEGs among the 56 defining the pro-inflammatory signature of OACs and TNF and IL-1b responses in NCs. (b) Protein–protein network analysis using the pro-inflammatory signature identified in OACs (c) and using the 21 DEGs also found in response to IL-1b and TNF treatments in NCs. (b, c) Three interacting groups were revealed by this analysis: The blue circle group was characterized by common pro-inflammatory targets, the yellow circle group was associated with proteins involved in antigen presentation, and the green box group was composed of metabolism-related proteins. (d) RT-qPCR analysis of Glut-1, MCT-4, HK2, and GFAT2 expression in OACs and NCs after 24 h of IL-1b or TNF treatment; *P < 0.05, ** P < 0.01 using the Mann–Whitney test; N = 5 per group. Values are expressed as means ± SEM

Full size image

Since Glut-1, MCT-4, and HK2 are involved in the key steps controlling glycolytic flux [32], we used the glycolysis/gluconeogenesis-related KEGG pathway annotation (hsa00010) to construct heatmaps and visualize dysregulations occurring in OACs versus NCs in response to pro-inflammatory cytokines (Fig. 5A). This analysis revealed an overall increase in the expression of glycolysis enzymes specifically in OACs treated with either IL-1b or TNF (Fig. 5A). As the final product of glycolysis is pyruvate, which is a metabolite that can initiate the tricarboxylic acid cycle (TCA, or citrate/Krebs cycle), we examined TCA-related genes using unsupervised heatmap analysis based on the TCA-related KEGG annotation (hsa00020). The TCA heatmap showed an overall decrease in TCA-related genes in OACs but not in NCs (Fig. 5B). To strengthen our results, glycolytic and TCA gene dysregulation were also explored in two public datasets: GSE162510, which explored the transcriptomic changes upon IL-1b treatment, and E-MTAB-6266, which compared gene expression between the cartilage of OA knees and the cartilage of healthy knees from multi-organ donors as a control. Interestingly, both datasets confirmed the expression profile highlighted in our study, with overexpression of the genes associated with glycolysis (Fig. 5C), and downregulation of the genes related to TCA (Fig. 5D) in OA conditions compared to the controls (untreated chondrocytes or chondrocytes from non-OA joints).

Fig. 5
figure 5

Evidence of an OAC-specific metabolic switch. (a, b) Heatmaps displaying the log2(Fold-Change) in OACs or NCs for genes annotated in the KEGG database for (a) glycolysis/neo-glucogenesis and (b) the Krebs cycle (TCA). Heatmaps from the public transcriptome dataset (E-MTAB-6266 and GSE162510) highlighting glycolysis (c) and TCA (d) gene modulation. e Lactate measurement in OAC (N = 11) and NC (N = 11) supernatants; ns.: not significant; *P < 0.05 using the Mann–Whitney test; N = 11 per group. Values are expressed as means ± SEM

Full size image

To better understand the impact of this glycolytic increase on the overall cellular metabolism, we also explored other metabolic pathways directly related to glycolysis, such as the HBP and the pentose phosphate pathway. No regulation as clear and specific as the one observed for glycolysis was revealed for these pathways (Figure S2A, B). Moreover, in keeping with the increase observed in the expression of glycolysis enzymes and the decrease in the expression of TCA enzymes, our results revealed overexpression of LDHA and, as previously described, overexpression of MCT-4 (Slc16a3) at both the transcriptomic and proteomic levels in response to either IL-1b or TNF in OACs (Fig. 5A, Table 1). Based on these observations, we examined whether lactate was abnormally secreted in OACs treated with IL-1b or TNF, and we found a significant increase in the extracellular lactate concentrations in treated OACs compared to NCs (Fig. 5E). These data suggest a metabolic shift of OACs in response to pro-inflammatory cytokines.

Real-time cell metabolic analysis confirmed induction of OAC-specific metabolism by pro-inflammatory stimuli

To strengthen our multi-omic analysis, we then examined whether the metabolic shift has a functional impact on cellular metabolism. This was addressed by performing Seahorse® real-time metabolic assays. We evaluated both the extracellular acidification rate (ECAR) during a glycolytic stress test and the oxygen consumption rate (OCR) during a mitochondrial stress test in both OACs and NCs upon IL-1b and TNF treatments. As expected, according to previous results presented in Fig. 5, OACs and NCs had distinct metabolic profiles. Indeed, the ECAR measured during the glycolytic stress assay revealed an altered profile of OACs treated with IL-1b or TNF compared with untreated cells, while this was not observed in NCs (Fig. 6A, B). Unexpectedly, OAC and NC samples did not exhibit differences in glycolysis levels, irrespective of the applied treatments (Fig. 6C). However, pro-inflammatory stimuli significantly decreased the glycolytic reserve of OACs but not NCs (Fig. 6D). Interestingly, non-glycolytic acidification was increased only in treated OAC samples but not in treated NC samples (Fig. 6E), suggesting that other pathways are involved in the acidification of the OAC environment.

Fig. 6
figure 6

Investigation of OAC-specific functional metabolic changes. a-e Seahorse® assay evaluating ECAR measurements in the glycolysis stress test in OACs (a) and NCs (b). c Quantification of glycolysis, (d) glycolytic reserve, and (e) non-glycolytic acidification from the glycolysis stress tests. f-j) Seahorse® assay evaluating OCR measurements in the mitochondrial stress assay in OACs (f) and NCs (g). (h) Quantification of basal respiration, (i) ATP production, and (j) non-mitochondrial oxygen consumption (NMOC) from the mitochondrial stress assay. k Nitrite measurements in OAC and NC supernatants; *P < 0.05, **P < 0.01 or ***P < 0.001 using the Mann–Whitney test; N = 5 (Seahorse® assays) or N = 8 (nitrite measurements) per group. Values are expressed as means ± SEM

Full size image

In parallel, the mitochondrial stress assay revealed altered oxidative phosphorylation (OXPHOS) in response to IL-1b or TNF in OACs (Fig. 6F) compared to the expected NC profiles (Fig. 6G). Consequently, the OAC basal respiration level was significantly decreased by pro-inflammatory stimuli (Fig. 6H), and it was associated with a significant reduction in ATP production by OXPHOS (F ig. 6I). Consistent with our RNA-seq data, pro-inflammatory stimuli did not alter OXPHOS or ATP production in NCs (Fig. 6H-J). Interestingly, non-mitochondrial O2 consumption (NMOC) increased specifically in OACs upon IL-1b and TNF treatments (Fig. 6J). An increased NMOC suggests that OACs treated with cytokines consumed O2 in other ways than OXPHOS, such as NO and/or ROS production. Of note, the nitrite (NO2−) concentration was specifically increased in OAC supernatants by pro-inflammatory stimuli (Fig. 6K). Since nitrite production resulted from NO oxidation, part of the increased NMOC observed in OACs but not in NCs could, therefore, be explained by NO production.

Materials and methods

Experimental design

Primary human chondrocytes, isolated from osteoarthritic knees of advanced OA patients (OACs) and vertebral transverse costal facet joints of non-OA patients (NCs). For histological validation, cartilage samples were paraffin-embedded and stained with Safranin-O. Isolated OA chondrocytes (OACs) and non-OA chondrocytes (NCs) were treated for 24 h with or without 1 ng/mL of IL-1b or 25 ng/mL of TNF or 50 ng/mL of IL-6/IL-6 soluble receptor (IL-6). Sequencing by 3’ SRP was performed on OAC and NC samples. Moreover, mass-spectrometry was performed on OACs. The pro-inflammatory signature was defined with matched transcripts and proteins dysregulated by IL-1b or TNF treatments. NC transcriptomics data were then compared to the OAC pro-inflammatory signature (Fig. 1).

Cell culture

Human chondrocytes

Primary human chondrocytes were isolated from osteoarthritic joints of advanced OA patients undergoing total knee surgery (OACs, 71–81 years old, patient number indicated in the figure legends) and from vertebral transverse costal facet joints of non-OA patients undergoing scoliosis surgery (NCs, patient number indicated in the figure legends, 15–21 years old). According to the Declaration of Helsinki, all human articular cartilage samples were harvested from patients after they had provided informed consent. This study was carried out following the recommendation of the Committee for Person Protection of Pays de La Loire and approved by the French Ministry of Higher Education and Research (registration number: DC-2017–2987). Briefly, cartilage was cut into small slices and placed in Hanks' Balanced Salt Solution (HBSS, L0606-500, Biowest) supplemented with 10% penicillin–streptomycin solution (PS, 1000 U/mL, 15,140,122, Thermo Fisher Scientific) at 37 °C in a humidified atmosphere with 5% CO2 for up to 5 h. Then, small slices were digested for 30 min in 0.2% collagenase (125 U/mg, C6885, Sigma-Aldrich), and chondrocytes were released by overnight digestion (~ 14 h) in 0.03% collagenase in high-glucose Dulbecco's Modified Eagle’s Medium (DMEM, 31,966–021, Gibco) containing 10% fetal calf serum (FCS, Dominique Dutscher) and 1% PS (complete medium) and passed through a 70 µm filter to remove undigested cartilage fragments. Isolated chondrocytes were seeded at a density of 6.0 × 104 cells/cm2 for OACs and 1.5 × 104 cells/cm2 for NCs and used at passage 1 after reaching 90% confluence.

Chondrocyte treatments

Primary human chondrocytes (OACs and NCs) were seeded at 6 × 104 cells/cm2 in complete medium and cultured at 37 °C in a humidified atmosphere with 5% CO2. After 24 h, cytokines were added as follows: 1 ng/mL IL-1b (Merck, IL038) or 25 ng/mL TNF (Miltenyi Biotec 130–094-023) or 50 ng/mL IL-6/IL-6 soluble receptor (R&D, 8954-SR-025) for 24 h. At the end of the treatments, the supernatants were collected and the cells were lysed before being stored immediately at -80 °C until further analysis.

For the transcriptomic analysis, total RNA was extracted from lysed chondrocytes using NucleoSpin RNA XS (740,902.50, Macherey–Nagel) according to the supplier’s instructions. The RNA yield and RNA quality were assessed for the 3’ SRP analysis using a NanoDrop™ spectrophotometer and an Agilent Bioanalyzer system (Agilent Technologies), respectively. The RNA Integrity Number was greater than 9 for each RNA sample.

For the mass spectrometry analysis, cells were trypsinized, washed with PBS, and then pelleted before being snap-frozen and stored at -80 °C until further analysis.

Chondrocytes analysis

qRT-PCR

RNAs were reverse transcribed to cDNA using a Verso™ cDNA Synthesis Kit (Thermo Fisher Scientific, AB-1453/B). The diluted cDNA was then amplified with specific primers (Table 2) using the SYBR® Select Master Mix (Thermo Fisher Scientific, 4,472,897). qPCR was then performed using a CFX96™ Touch Deep Well Real-Time PCR Detection System (Bio-Rad).

Table 2 Specific primers used for qPCR analysis
Full size table

Supernatants analysis

Nitrite concentrations were determined using the Griess reaction. Aliquots of 100 µL of cell supernatants were mixed with 100 µL of Griess solution (v/v mixture of a 0.1% aqueous solution of naphthylenediamine dihydrochloride and a 1% solution of sulfanilamide in 5% H3PO4). The absorbance was then measured at 550 nm, and the nitrite concentration was determined using a standard range of nitrites diluted in complete medium. Lactate concentrations were determined using a Lactate Assay Kit (Sigma-Aldrich, MAK064) with 1 µL of the cell supernatants following the manufacturer’s instructions.

Measurement of cellular metabolism by Seahorse® assay

Primary human chondrocytes (OACs or NCs) were seeded in Seahorse® XFp plates at 15 000 cells per well (1.3 × 105 cells/cm2) and treated with IL-1b (1 ng/mL) or TNF (25 ng/mL). After 24 h, the complete medium was replaced with DMEM (Sigma-Aldrich, 5030) with 1% P/S and glutamine (4 mM) for the glycolytic stress assay or glutamine (4 mM), glucose (25 mM), and pyruvate (1 mM) for the mitochondrial stress test. The plate was then preincubated at 37 °C for 45 min in the absence of CO2. The glycolytic stress assay consisted of the sequential addition of glucose (25 mM), oligomycin (1 µM), and 2-deoxy-D-glucose (2-DG, 50 mM) according to the manufacturer's instructions. The mitochondrial stress test consisted of the sequential addition of oligomycin (1 µM), carbonylcyanide-3-chlorophenylhydrazone (CCCP, 1 µM), and rotenone/antimycin A (1 µM each). Measurements of the extracellular acidification rate (ECAR) and the oxygen consumption rate (OCR) were performed every 5 min before and after the addition of the various drugs as indicated in the figures. The data analysis was performed using Wave software (Agilent).

3’ Seq-RNA Profiling (3’ SRP)

Preparation of the sequencing library

The 3'seq-RNA profiling was performed as described by Charpentier et al., 2021 and (48). Samples were randomly plated and diluted to a concentration of 2.5 ng/μL to construct the library. The library was prepared from 10 ng of total RNA in 4 µL. The mRNA poly-A tails were tagged with universal adapters, well-specific barcodes, and unique molecular identifiers (UMIs) during template-switching reverse transcriptase (Thermo Scientific™, EP0751, Maxima H Minus Reverse Transcriptase). Barcoded cDNAs from 96 samples were then pooled and purified using a Zymo kit. cDNAs were treated with exonuclease I (NEB, M0293S) before being amplified (12 cycles) and fragmented using a transposon-fragmentation approach that enriched for 3′ ends of cDNA. The fragments were size controlled on a 2200 TapeStation system (Agilent Technologies) to form a library of 350–800 bp in length. Finally, the library was validated with fragments of approximately 484 pb length and with 38.5 nM of properly tagged cDNA.

Sequencing and primary analysis

Sequencing and primary analysis were performed according to Charpentier et al., 2021. Briefly, the library was sequenced on a NovaSeq 6000 using NovaSeq 6000 SP Reagent Kit 100 cycles (#20,027,464, Illumina) with 17*-8–105* cycle reads (* addition of one cycle according to Illumina’s recommendation). Bioinformatics steps were performed using a Snakemake (49) pipeline. Sample sequences were demultiplexed using a python script, and raw paired-end fastq files were transformed into a single-end fastq file for each sample and aligned on the Ensembl transcriptome hg38 using bwa.

Public datasets analysis

Two public datasets, E-MTAB-6266 and GSE162510, were uploaded to explore the regulation of metabolism-associated genes. E-MTAB-6266 compared gene expression between the cartilage of 58 OA knees and 10 non-OA knees. GSE162510 explored the transcriptomic changes occurring in chondrocytes from 10 OA patients following 24 h of IL-1b treatment (0.2 ng/mL). A differential expression analysis table of the E-MTAB-6266 public dataset was uploaded from Soul et al. [33], while GSE162510 was uploaded with the R package GEO query (2.62.2). Data visualization was then performed as previously described.

Mass-spectrometry

Proteome sample preparation and tandem mass spectrometry

Cell pellets were dissolved in 25 µL of LDS buffer (26.5 mM Tris HCl, 35.25 mM Tris base, 0.5% LDS, 2.5% glycerol, 0.13 mM EDTA, supplemented with 5% beta-mercaptoethanol) per mg of pellet. Samples were heated for 5 min at 99 °C in a thermomixer (Eppendorf). For each sample, a volume of 20 µL (i.e., 10 µg of protein) was subjected to a 5 min denaturing electrophoresis on a NuPAGE 4%–12% gradient gel with MES SDS as the running buffer (50 mM MES ([2-(N-morpholino) ethane sulfonic acid), 50 mM Tris Base, 0.1% SDS, 1 mM EDTA, pH 7.3). Following electrophoresis, the gel was briefly washed with Milli-Q® water, stained with SimplyBlue™ SafeStain (Thermo Fisher Scientific) for 5 min to visualize proteins, and then washed extensively with Milli-Q® water. Each proteome was extracted as a single polyacrylamide band with a volume of approximately 100 µL. Each sample was processed as previously described (50) and then proteolyzed with Trypsin Gold (Promega) in 50 mM NH4HCO3 in the presence of ProteaseMAX™ detergent (Promega). A volume of 10 µL of the resulting peptide mixture (45 µL) was injected into a nanoscale C18 PepMap™ 100 capillary column (3 µm, 100 Å, 75 µm internal diameter × 50 cm length, LC Packings) mounted with a desalting pre-column and resolved with a 65-min gradient of acetonitrile (4.8%–24.8% in 60 min followed by 24.8%–40% in 5 min), 0.1% formic acid, at a flow rate of 0.25 µL/min. The peptides resolved by reverse-phase chromatography were analyzed by tandem mass spectrometry with an Exploris™ 480 mass spectrometer (Thermo Fisher Scientific) connected directly to the column exit. The instrument was operated in data-dependent acquisition mode, with a full scan of peptide ions acquired at a resolution of 120,000 from m/z 350 to 1500 and with a dynamic exclusion of 10 s. Each MS scan was followed by high-energy collisional dissociation and MS/MS scans at a resolution of 15,000 on the 20 most abundant precursor ions identified within the full scan, selecting only ions with a charge of 2+ or 3+.

Label-free shotgun proteomic interpretation

MS/MS spectra were assigned to peptide sequences by the MASCOT Daemon 2.3.2 search engine (Matrix Science) in follow-up mode using first the contaminant database described by Pereira et al. (51) including the 23 most abundant proteins from fetal calf serum, and then the SwissProt Homo sapiens database (20,396 polypeptide sequences). Standard search parameters included: trypsin as a proteolytic enzyme with two possible miss-cleavages at maximum, tolerances of 5 ppm and 0.02 Da for the MS and MS/MS signals, respectively, oxidation of methionine, deamidation of glutamine and asparagine, and acetylation of N-termini, as variable modifications, carbamidomethylation of cysteine as fixed modification, and a peptide p-value below 0.01. A protein was considered validated when at least two different peptides were detected, resulting in a protein identification false discovery rate below 1% as verified with a reverse decoy database search. To account for possible genome heterogeneity between patients, only peptides common between all patients (at least identified once in a patient sample) were used for the comparative proteomics. Spectral counts corresponding to the number of MS/MS spectra assigned per protein were used as a proxy for the abundance of the proteins in each condition, considering only unambiguous peptides (listed only once in the database).

Omic data processing

Differential expression analyses

Aligned reads were parsed, and UMIs were counted for each gene in each sample to create an expression matrix containing the absolute abundance of mRNAs in all samples. Reads aligned on multiple genes or containing more than three mismatches with the reference were discarded. The expression matrix was normalized and differentially expressed genes (DEG) were searched using the R package DeSeq2 (52) (thresholds used: adjusted p-value < 0.01 and log2(Fold-Change) > 0.58). Supplementary R packages compiled with R version 4.1.1 were used to generate data visualization: ggplot2 (3.3.5), VennDiagram (1.7.3), enrichR (3.0), and ComplexHeatmap (2.8.0).

Protein–protein network

Protein–protein network analysis, using the 56 DEGs/DEPs identified in OACs and the 21 DEGs of NCs, were determined using the Search Tool for the Retrieval of Interacting Genes (STRING) database version 11.5. All available interactions were mapped with an interaction score cutoff of 0.4, and unconnected nodes were removed.

Statistical analysis

The statistical analyses used in the transcriptomic analysis were performed using the R software package DESeq2 (52) with a Benjamini & Hochberg correction. The PCA and biplots were created using the prcomp function from the stats package within R statistics (version 4.1.1) and the factoextra R package (1.0.7). Differences in continuous variables were assessed using the non-parametric Mann–Whitney test (n < 30), and p-values less than 0.05 were considered to be statistically significant. The statistical analyses were performed using GraphPad Prism 8.0v software.

Discussion

Using a multi-omics approach, this study presents, a broad description of IL-1b and TNF responses of chondrocytes from OA joints compared to chondrocytes from non-OA joints. Our multi-omics approach based on a 3’-seq RNA analysis and mass spectrometry analysis of OAC lysates identified 56 targets of interest that define a canonical pro-inflammatory signature of OACs induced in response to either IL-1b or TNF. The protein–protein network analysis of the 56 targets identified in the canonical pro-inflammatory signature of OACs defined three main interactome groups. By comparing the OAC canonical pro-inflammatory signature to NCs, we identified that the group corresponding to metabolism was absent in the pro-inflammatory response of NCs.

This metabolism group comprises four proteins: Glut-1, HK2, MCT-4, and GFAT2. Thus, a special focus on energetic metabolism pathways such as glycolysis and TCA at the transcriptomic level revealed that the expression of glycolytic enzymes and TCA enzymes are inversely regulated in OACs. Indeed, our study, showed, that the expression of most genes involved in glycolysis is increased whereas the expression of those related to TCA is decreased in OACs upon exposure to IL-1b or TNF. This increased expression of glycolytic genes such as Glut-1, HK2, or LDHA has been previously described in the pro-inflammatory response of OA and non-OA chondrocytes [34,35,36,37]. However, in our study, no significant regulation of metabolism-related genes upon exposure to IL-1b or TNF was observed in non-OA chondrocytes (NCs). This finding was confirmed by a public data set obtained from chondrocytes isolated from knees with OA and healthy knees from age-matched multi-organ donors (E-MTAB-6266) [33]. Interestingly, promising studies have also highlighted the potential benefit of glycolysis inhibition on the inflammatory response of chondrocytes [38, 39].

We further examined the bioenergetic profiles of chondrocytes using real-time cell metabolic assays. Interestingly, we confirmed the functional-specific impact on OACs of the metabolic gene modulations identified in our transcriptomic analysis. Indeed, Seahorse® glycolytic stress assay showed that pro-inflammatory conditions forced OACs to use glycolysis at their full capacity, which was not observed in NCs. Furthermore, a significant increase in non-glycolytic acidification was specifically observed in OACs upon treatment with IL-1b or TNF. An acidic microenvironment is a hallmark of inflamed tissue as well as cartilage from OA patients [40] and has been associated with decreased glycosaminoglycan production and increased MMP activity in cartilage [41]. Microenvironment acidification could be due to TCA activity, glycogenolysis, or lactate export (linked to glycolysis) [42]. In our study, the increased expression of LDHA and MCT-4 expression suggests enhanced lactate secretion in OACs following exposure to pro-inflammatory stimuli. We confirmed that lactate concentrations increased in the supernatants of OACs compared to NCs under pro-inflammatory conditions. An increase in synovial fluid lactate levels has already been associated with inflammatory joint diseases, such as rheumatoid arthritis [43], and in OA animal models [44, 45]. Altogether, these data confirmed that disease-specific metabolic dysregulation characterizes OACs.

In parallel, our results showing decreased expression in TCA-related genes, basal respiration, and ATP production suggest an impaired mitochondrial function in OACs upon IL-1b and TNF treatments. Consistent with our data, several metabolic dysregulations have been shown previously in OA chondrocytes, such as OXPHOS impairments [18, 19, 24, 46]. Here, our study provides the additional information that pro-inflammatory stimuli exacerbate OXPHOS dysfunction in OA chondrocytes, whereas NCs appeared to be less sensitive to pro-inflammatory stimuli and thus did not exhibit metabolic alterations associated with glycolysis or OXPHOS. Contrary to our results, the previous study by Eitner et al. showed that IL-1b and TNF treatments induced much more pronounced metabolic alterations in non-OA chondrocytes compared to OA chondrocytes [24] after 48 h of treatment. However, these differences with our results could be explained by differences in treatment conditions (duration and dose), as well as in the anatomical source of chondrocytes (hip).

Interestingly, in OACs but not in NCs, both IL-1b and TNF induced an increase in NMOC associated with an increase in nitrite production, reflecting O2 consumption for NO production. As NO is a known chondrocyte catabolic product found in OA [47, 48], this increased nitrite production strengthens the association of these metabolic dysfunctions with pathogenic OA processes. Our observations, therefore, suggest that the inflammatory and/or catabolic environment chondrocytes experience in OA joints for several years makes them more responsive to metabolic changes induced by pro-inflammatory stimuli. Whether sustained treatment of NCs with low-grade pro-inflammatory stimuli could sensitize them to this metabolic shift merits further analysis.

Although our study focused on glycolysis and TCA due to its specificity to OACs, other pathways, such as the hexosamine biosynthetic pathway (HBP), are regulated by pro-inflammatory stimuli. This HBP produces the monosaccharide UDP-GlcNAc involved in GAG synthesis and the post-translational modification called O-GlcNacylation [49], integrates multiple metabolic pathways such as glucose, glucosamine, amino acid, fatty acid, and nucleotide metabolisms. Thus, it could also be explored in depth. Indeed, our study highlighted that the expression of several players of the HBP are up-regulated in both OAC and NC treated with IL-1b- or TNF. Unlike what is observed for glycolysis and OXPHOS, inflammation can induce a modification of the HBP in chondrocytes regardless of their status, OA or not. This suggests either that this metabolic pathway is not specifically altered in OA or that it may be altered early in the development of OA. Few studies have explored the involvement of this HBP in OA [49, 50], mainly focusing on its interplay with glucosamine, a known but controversial component in the management of OA.

Interestingly our study also highlighted a lack of relevant transcriptomic and proteomic changes by IL-6 exposure in human chondrocytes (osteoarthritic and non-osteoarthritic) under the experimental conditions tested. This result may seem surprising as studies have described a protective effect of IL-6 inhibition in preclinical models of OA [12], suggesting that IL-6 is a key signaling cytokine in joints. However, anti-IL-6 clinical trials have been ineffective in OA, highlighting the species variability in IL-6 signaling. This is consistent with our data, where a transcriptomic response to the same IL-6 treatment was observed in murine chondrocytes (530 DEGs, data not shown) but not in human chondrocytes.

The design of our study does not completely rule out the possibility that the observed differences in response to pro-inflammatory stimuli are related to the use of chondrocytes from different joints. However, all our observations were supported by additional studies using publicly available datasets. Nevertheless, since accessing a healthy knee joint is challenging, supplementation of this study with other animal models such as dogs or horses would be substantial. These future studies should also complement the pro-inflammatory signature with data on secretome and metabolome dysregulations. Despite these limitations, the main interest of our study lies in the multi-omics approach completed by a functional approach to investigate the action of two different pro-inflammatory treatments.

Our study identified a canonical pro-inflammatory signature and its specific signaling pathways in chondrocytes from OA joints compared to chondrocytes from non-OA joints. Although IL-1b or TNF treatments are widely used to mimic the pro-inflammatory component of OA, the specificity of the metabolic alterations observed only in OACs strengthens the need for a thorough study of the link between inflammation and metabolism. Indeed, our study suggests that the pathological history of OA chondrocytes (inflammation, pre-established metabolic alterations) may promote abnormal catabolic responses to stress and thus sustain joint degeneration. Finally, based on the relative resistance of NCs to metabolic shifts induced by proinflammatory stimuli, further studies may help design relevant strategies to efficiently target upstream processes to inflammation-related metabolic change.

Availability of data and materials

The mass spectrometry and proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE (53) partner repository with the dataset identifiers PXD037145 and 10.6019/PXD037145. All sequencing 3' SRP data are publicly available via the NCBI Gene Expression Omnibus (GEO) using the accession number GSE215039.

References

  1. Leifer VP, Katz JN, Losina E. The burden of OA-health services and economics. Osteoarthr Cartil. 2022;30:10–6.

    Article  CAS  Google Scholar 

  2. Hunter DJ, Bierma-Zeinstra S. Osteoarthritis Lancet. 2019;393:1745–59.

    CAS  PubMed  Google Scholar 

  3. Martel-Pelletier J, et al. Osteoarthritis. Nat Rev Dis Prim. 2016;2:16072.

    Article  PubMed  Google Scholar 

  4. Mathiessen A, Conaghan PG. Synovitis in osteoarthritis: Current understanding with therapeutic implications. Arthritis Res Ther. 2017;19:1–9.

    Article  Google Scholar 

  5. Jenei-Lanzl Z, Meurer A, Zaucke F. Interleukin-1β signaling in osteoarthritis – chondrocytes in focus. Cell Signal. 2019;53:212–23.

    Article  CAS  PubMed  Google Scholar 

  6. Ansari MY, Ahmad N, Haqqi TM. Oxidative stress and inflammation in osteoarthritis pathogenesis: role of polyphenols. Biomed Pharmacother. 2020;129:110452.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Chou CH, et al. Synovial cell cross-talk with cartilage plays a major role in the pathogenesis of osteoarthritis. Sci Rep. 2020;10:10868.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Lv M, et al. Identification of chondrocyte genes and signaling pathways in response to acute joint inflammation. Sci Rep. 2019;9:1–12.

    Article  Google Scholar 

  9. Wang Q, et al. Celecoxib prevents tumor necrosis factor-α (TNF-α)-induced cellular senescence in human chondrocytes. Bioengineered. 2021;12:12812–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Lepetsos P, Papavassiliou KA, Papavassiliou AG. Redox and NF-κB signaling in osteoarthritis. Free Radic Biol Med. 2019;132:90–100.

    Article  CAS  PubMed  Google Scholar 

  11. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Latourte A, et al. Systemic inhibition of IL-6/Stat3 signalling protects against experimental osteoarthritis. Ann Rheum Dis. 2017;76:748–55.

    Article  CAS  PubMed  Google Scholar 

  13. Wang F, et al. IL-1β receptor antagonist (IL-1Ra) combined with autophagy inducer (TAT-Beclin1) is an effective alternative for attenuating extracellular matrix degradation in rat and human osteoarthritis chondrocytes. Arthritis Res Ther. 2019;21:171.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Chisari E, Yaghmour KM, Khan WS. The effects of TNF-alpha inhibition on cartilage: a systematic review of preclinical studies. Osteoarthritis and Cartilage at. 2019. https://doi.org/10.1016/j.joca.2019.09.008.

    Article  Google Scholar 

  15. Richette P, et al. Efficacy of tocilizumab in patients with hand osteoarthritis: double blind, randomised, placebo-controlled, multicentre trial. Ann Rheum Dis. 2021;80:349–55.

    Article  CAS  PubMed  Google Scholar 

  16. Chevalier X, Eymard F. Anti-IL-1 for the treatment of OA: dead or alive? Nat Rev Rheumatol. 2019;15:191–2.

    Article  PubMed  Google Scholar 

  17. Kloppenburg M, et al. Etanercept in patients with inflammatory hand osteoarthritis (EHOA): a multicentre, randomised, double-blind, placebo-controlled trial. Ann Rheum Dis. 2018;77:1757–64.

    Article  CAS  PubMed  Google Scholar 

  18. Maneiro E, et al. Mitochondrial respiratory activity is altered in osteoarthritic human articular chondrocytes. Arthritis Rheum. 2003;48:700–8.

    Article  CAS  PubMed  Google Scholar 

  19. Yang X, et al. Pyruvate Kinase M2 Modulates the Glycolysis of Chondrocyte and Extracellular Matrix in Osteoarthritis. DNA Cell Biol. 2018;37:271–7.

    Article  CAS  PubMed  Google Scholar 

  20. Qu J, et al. PFKFB3 modulates glycolytic metabolism and alleviates endoplasmic reticulum stress in human osteoarthritis cartilage. Clin Exp Pharmacol Physiol. 2016;43:312–8.

    Article  CAS  PubMed  Google Scholar 

  21. Yao X, et al. Fibroblast growth factor 18 exerts anti-osteoarthritic effects through PI3K-AKT signaling and mitochondrial fusion and fission. Pharmacol Res. 2019;139:314–24.

    Article  CAS  PubMed  Google Scholar 

  22. Farah H, et al. Differential Metabotypes in Synovial Fibroblasts and Synovial Fluid in Hip Osteoarthritis Patients Support Inflammatory Responses. Int J Mol Sci. 2022;23:3266.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Garcia-Carbonell R, et al. Critical Role of Glucose Metabolism in Rheumatoid Arthritis Fibroblast-like Synoviocytes. Arthritis Rheumatol. 2016;68:1614–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Eitner A, et al. Osteoarthritis-Induced Metabolic Alterations of Human Hip Chondrocytes. Biomedicines. 2022;10:1349.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Bolduc JA, Collins JA, Loeser RF. Reactive oxygen species, aging and articular cartilage homeostasis. Free Radic Biol Med. 2019;132:73–82.

    Article  CAS  PubMed  Google Scholar 

  26. Conde J, et al. Identification of novel adipokines in the joint. Differential expression in healthy and osteoarthritis tissues. PLoS One. 2015;10:e0123601.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Ajekigbe B, et al. Identification of long non-coding RNAs expressed in knee and hip osteoarthritic cartilage. Osteoarthr Cartil. 2019;27:694.

    Article  CAS  Google Scholar 

  28. Kim J-HH, et al. Regulation of the catabolic cascade in osteoarthritis by the zinc-ZIP8-MTF1 axis. Cell. 2014;156:730–43.

    Article  CAS  PubMed  Google Scholar 

  29. Zhang FJ, Luo W, Lei GH. Role of HIF-1α And HIF-2α In osteoarthritis. Jt Bone Spine. 2015;82:144–7.

    Article  CAS  Google Scholar 

  30. Sengprasert P, et al. Upregulation of antigen presentation function and inflammation in chondrocytes by induction of proteoglycan aggrecan peptides (P16–31 and P263–280). Clin Exp Rheumatol. 2022;40:596–607.

    Article  PubMed  Google Scholar 

  31. Lv Z, et al. Single cell RNA-seq analysis identifies ferroptotic chondrocyte cluster and reveals TRPV1 as an anti-ferroptotic target in osteoarthritis. EBioMedicine. 2022;84:104258.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Tanner LB, et al. Four key steps control glycolytic flux in mammalian cells. Cell Syst. 2018;7:49-62.e8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Soul J, et al. Stratification of knee osteoarthritis: Two major patient subgroups identified by genome-wide expression analysis of articular cartilage. Ann Rheum Dis. 2018;77:423–30.

    Article  CAS  PubMed  Google Scholar 

  34. Shikhman AR, Brinson DC, Lotz MK. Distinct pathways regulate facilitated glucose transport in human articular chondrocytes during anabolic and catabolic responses. Am J Physiol - Endocrinol Metab. 2004;286:980–5.

    Article  Google Scholar 

  35. Bustamante MF, et al. Hexokinase 2 as a novel selective Metabolic target for rheumatoid arthritis. Ann Rheum Dis. 2018;77:1636.

    Article  CAS  PubMed  Google Scholar 

  36. Rosa SC, et al. Impaired glucose transporter-1 degradation and increased glucose transport and oxidative stress in response to high glucose in chondrocytes from osteoarthritic versus normal human cartilage. Arthritis Res Ther. 2009;11:R80.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Arra M, et al. LDHA-mediated ROS generation in chondrocytes is a potential therapeutic target for osteoarthritis. Nat Commun. 2020;11:1–16.

    Article  Google Scholar 

  38. Ohashi Y, et al. Metabolic reprogramming in chondrocytes to promote mitochondrial respiration reduces downstream features of osteoarthritis. Sci Rep. 2021;11:1–16.

    Article  Google Scholar 

  39. Yu SM, Kim SJ. Endoplasmic reticulum stress (ER-stress) by 2-deoxy-d-glucose (2DG) reduces cyclooxygenase-2 (COX-2) expression and N-glycosylation and induces a loss of COX-2 activity via a Src kinase-dependent pathway in rabbit articular chondrocytes. Exp Mol Med. 2010;42:777–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Konttinen YT, et al. Acidic cysteine endoproteinase cathepsin K in the degeneration of the superficial articular hyaline cartilage in osteoarthritis. Arthritis Rheum. 2002;46:953–60.

    Article  CAS  PubMed  Google Scholar 

  41. Razaq S, Wilkins RJ, Urban JPG. The effect of extracellular pH on matrix turnover by cells of the bovine nucleus pulposus. Eur Spine J. 2003;12:341.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Pike Winer LS, Wu M. Rapid analysis of glycolytic and oxidative substrate flux of cancer cells in a microplate. PLoS One. 2014;9:e109916.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Haas R, et al. lactate regulates metabolic and pro-inflammatory circuits in control of t cell migration and effector functions. PLoS Biol. 2015;13:e1002202.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Mickiewicz B, et al. Metabolic analysis of knee synovial fluid as a potential diagnostic approach for osteoarthritis. J Orthop Res. 2015;33:1631–8.

    Article  PubMed  Google Scholar 

  45. Damyanovich AZ, Styaples JR, Marshall KW. 1H NMR investigation of changes in the metabolic profile of synovial fluid in bilateral canine osteoarthritis with unilateral joint denervation. Osteoarthr Cartil. 1999;7:165–72.

    Article  CAS  Google Scholar 

  46. Wu X, et al. Dysregulated energy metabolism impairs chondrocyte function in osteoarthritis. Osteoarthr Cartil. 2022. https://doi.org/10.1016/J.JOCA.2022.11.004.

    Article  Google Scholar 

  47. Yang Z, et al. Large-scale generation of functional mRNA-encapsulating exosomes via cellular nanoporation. Nat Biomed Eng. 2019. https://doi.org/10.1038/s41551-019-0485-1.

    Article  PubMed  PubMed Central  Google Scholar 

  48. López-Armada MJ, et al. Mitochondrial activity is modulated by TNFα and IL-1β in normal human chondrocyte cells. Osteoarthr Cartil. 2006;14:1011–22.

    Article  Google Scholar 

  49. Riegger J, Baumert J, Zaucke F, Brenner RE. The hexosamine biosynthetic pathway as a therapeutic target after cartilage trauma: modification of chondrocyte survival and metabolism by glucosamine derivatives and PUGNAc in an Ex Vivo Model. Int J Mol Sci. 2021;22:7247.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Gouze JN, et al. Adenovirus-mediated gene transfer of glutamine: fructose-6-phosphate amidotransferase antagonizes the effects of interleukin-1beta on rat chondrocytes. Osteoarthr Cartil. 2004;12:217–24.

    Article  Google Scholar 

Download references

Acknowledgements

We are most grateful to the Genomics Core Facility GenoA, member of Biogenouest and France Genomique and to the Bioinformatics Core Facility BiRD, member of Biogenouest and the Institut Français de Bioinformatique (IFB) (ANR-11-INBS-0013) for the use of their resources and their technical support.

We thank Jean-Charles Gaillard (ProGénomix Platform, CEA-Marcoule) for help with the mass spectrometry analysis.

Funding

This work was supported by the Agence Nationale de la Recherche (ANR) through the project KLOTHOA (ANR-18-CE14-0024–01; CV), by the region “Pays de la Loire” (RFI Bioregate, GenOA, CV, AD), and by research grants from the Arthritis Foundation and the “Société Française de Rhumatologie” (CV).

Author information

Authors and Affiliations

  1. Nantes Université, Oniris, CHU Nantes, INSERM, Regenerative Medicine and Skeleton, RMeS, UMR 1229, 44000, Nantes, France

    Anais Defois, Nina Bon, Alexandre Charpentier, Melina Georget, Nicolas Gaigeard, Frederic Blanchard, Antoine Hamel, Denis Waast, Yves Maugars, Marie-Astrid Boutet, Jerome Guicheux & Claire Vinatier

  2. Département Médicaments Et Technologies Pour La Santé (DMTS), Université Paris-Saclay, CEA, INRAE, SPI, Bagnols-Sur-Cèze, 30200, France

    Jean Armengaud

  3. Nantes Université, INSERM, CNRS, CRCI2NA, F-44000, Nantes, France

    Ophelie Renoult & Claire Pecqueur

  4. Centre for Experimental Medicine & Rheumatology, William Harvey Research Institute and Barts and The London School of Medicine and Dentistry, Queen Mary University of London, London, UK

    Marie-Astrid Boutet

Authors
  1. Anais Defois
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Nina Bon
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Alexandre Charpentier
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Melina Georget
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Nicolas Gaigeard
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Frederic Blanchard
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Antoine Hamel
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Denis Waast
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Jean Armengaud
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Ophelie Renoult
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Claire Pecqueur
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Yves Maugars
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Marie-Astrid Boutet
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Jerome Guicheux
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Claire Vinatier
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Conceptualization: AD, YM, CV, and JG; Methodology: AD, CP, CV, FB, JA, JG, MAB, and NB; Investigation: AC, AD, JA, MG, NB, NG, and OR; Samples providers: AH and DW; Supervision: AD, CV, JG, and YM; Writing—original draft: AD and CV; All authors reviewed the manuscript. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jerome Guicheux or Claire Vinatier.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Comparison of OACs and NCs pro-inflammatory responses. (a,b) Venn diagram highlighting OAC and NC matched DEGs in response to (a) 1 ng/mL of IL-1b or (b) 25 ng/mL of TNF

Additional file 2: Fig. S2.

Glycolysis-linked metabolic pathways. (a,b) Heatmaps displaying the log2 (Fold-Change) in OACs and NCs for genes associated with (a) hexosamine biosynthesis pathway and (b) pentose phosphate pathway. ns.: not significant.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Defois, A., Bon, N., Charpentier, A. et al. Osteoarthritic chondrocytes undergo a glycolysis-related metabolic switch upon exposure to IL-1b or TNF. Cell Commun Signal 21, 137 (2023). https://doi.org/10.1186/s12964-023-01150-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12964-023-01150-z

Keywords

Cell Communication and Signaling

ISSN: 1478-811X