Genes involved in the osteoarthritis process identified through genome wide expression analysis in articular cartilage; the RAAK study

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

OBJECTIVE

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

DESIGN

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

RESULTS

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

CONCLUSIONS

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

PPAR γ activation in chondrocytes alleviates glucocorticoid-induced oxidative stress, mitochondrial impairment, and pyroptosis via autophagic flow enhancement

Osteoarthritis (OA) is a progressive age-related disease characterised by pathological changes in the synovium, articular cartilage, and subchondral bone, significantly reducing the patients' quality of life. This study investigated the role of glucocorticoids, specifically dexamethasone, in OA progression, with a particular focus on their effects on chondrocytes. Although glucocorticoids are commonly used for OA pain relief, our research demonstrated that high concentrations of dexamethasone may accelerate OA progression by enhancing the ability of reactive oxygen species to inhibit chondrocyte autophagy, resulting in cell death and accelerated cartilage degeneration. Despite reports on the acceleration of pathogenesis and cartilage damage in some patients of OA taking corticosteroids, the mechanism behind the same has not been investigated. This necessitates an investigation of the concentration-dependent changes in the cartilage cells upon dexamethasone administration. In addition, the protective effect of PPAR γ on chondrocytes can prevent the decrease in chondrocyte autophagy and delay cartilage degeneration. Therefore, our study suggests that the therapeutic use of glucocorticoids in OA treatment should be more nuanced considering their potential detrimental effects. Future investigations should focus on the mechanisms underlying the glucocorticoid-mediated modulation of cell death processes, which could provide insights into new therapeutic strategies for OA treatment.

Transcriptomic analyses and machine-learning methods reveal dysregulated key genes and potential pathogenesis in human osteoarthritic cartilage

Aims

This study aimed to explore the biological and clinical importance of dysregulated key genes in osteoarthritis (OA) patients at the cartilage level to find potential biomarkers and targets for diagnosing and treating OA.

Methods

Six sets of gene expression profiles were obtained from the Gene Expression Omnibus database. Differential expression analysis, weighted gene coexpression network analysis (WGCNA), and multiple machine-learning algorithms were used to screen crucial genes in osteoarthritic cartilage, and genome enrichment and functional annotation analyses were used to decipher the related categories of gene function. Single-sample gene set enrichment analysis was performed to analyze immune cell infiltration. Correlation analysis was used to explore the relationship among the hub genes and immune cells, as well as markers related to articular cartilage degradation and bone mineralization.

Results

A total of 46 genes were obtained from the intersection of significantly upregulated genes in osteoarthritic cartilage and the key module genes screened by WGCNA. Functional annotation analysis revealed that these genes were closely related to pathological responses associated with OA, such as inflammation and immunity. Four key dysregulated genes (cartilage acidic protein 1 (CRTAC1), iodothyronine deiodinase 2 (DIO2), angiopoietin-related protein 2 (ANGPTL2), and MAGE family member D1 (MAGED1)) were identified after using machine-learning algorithms. These genes had high diagnostic value in both the training cohort and external validation cohort (receiver operating characteristic > 0.8). The upregulated expression of these hub genes in osteoarthritic cartilage signified higher levels of immune infiltration as well as the expression of metalloproteinases and mineralization markers, suggesting harmful biological alterations and indicating that these hub genes play an important role in the pathogenesis of OA. A competing endogenous RNA network was constructed to reveal the underlying post-transcriptional regulatory mechanisms.

Conclusion

The current study explores and validates a dysregulated key gene set in osteoarthritic cartilage that is capable of accurately diagnosing OA and characterizing the biological alterations in osteoarthritic cartilage; this may become a promising indicator in clinical decision-making. This study indicates that dysregulated key genes play an important role in the development and progression of OA, and may be potential therapeutic targets.

Mechanical stress and inflammation have opposite effects on Wnt signaling in human chondrocytes

Dysregulation of Wingless and Int-1 (Wnt) signaling has been strongly associated with development and progression of osteoarthritis (OA). Here, we set out to investigate the independent effects of either mechanical stress (MS) or inflammation on Wnt signaling in human neocartilage pellets, and to relate this Wnt signaling to OA pathophysiology. OA synovium-conditioned media (OAS-CM) was collected after incubating synovium from human end-stage OA joints for 24 h in medium. Cytokine levels in the OAS-CM were determined with a multiplex immunoassay (Luminex). Human neocartilage pellets were exposed to 20% MS, 2% OAS-CM or 1 ng/mL Interleukin-1β (IL-1β). Effects on expression levels of Wnt signaling members were determined by reverse transcription-quantitative polymerase chain reaction. Additionally, the expression of these members in articular cartilage from human OA joints was analyzed in association with joint space narrowing (JSN) and osteophyte scores. Protein levels of IL-1β, IL-6, IL-8, IL-10, tumor necrosis factor α, and granulocyte-macrophage colony-stimulating factor positively correlated with each other. MS increased noncanonical WNT5A and FOS expression. In contrast, these genes were downregulated upon stimulation with OAS-CM or IL-1β. Furthermore, Wnt inhibitors DKK1 and FRZB decreased in response to OAS-CM or IL-1β exposure. Finally, expression of WNT5A in OA articular cartilage was associated with increased JSN scores, but not osteophyte scores. Our results demonstrate that MS and inflammatory stimuli have opposite effects on canonical and noncanonical Wnt signaling in human neocartilage. Considering the extent to which MS and inflammation contribute to OA in individual patients, we hypothesize that targeting specific Wnt pathways offers a more effective, individualized approach.

CD19, ALDH18A1, and CACNA1G as Significant Hub Genes in End-Stage Osteoarthritis

Background

Osteoarthritis is one of the principal causes of chronic joint disease and may progressively engender disability in elderly individuals. The present study aimed to identify differentially expressed genes and associated signaling pathways in end-stage osteoarthritis.

Methods

Differentially expressed messenger RNAs in the early and end stages of osteoarthritis were examined through gene expression omnibus 2R (GEO2R) in the GSE32317 dataset. Subsequently, gene ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG), and protein-protein interaction (PPI) analyses were conducted. Furthermore, microRNAs targeting hub genes were investigated using the miRcode database. This study was conducted jointly at Bam University of Medical Sciences and Rajaie Cardiovascular, Medical and Research Center on October 2022.

Results

Differentially expressed data demonstrated downregulation in 134 genes and upregulation in 189 genes in end-stage knee osteoarthritis. The results of the enrichment and PPI analyses determined 4 end-stage knee osteoarthritis-related hub genes: IL-1B, CD19, CACNA1G, and ALDH18A1. The knee osteoarthritis-related key genes were involved in the Wnt signaling, B cell receptor signaling, calcium signaling, circadian entrainment, arginine and proline metabolism, axon guidance, and cytokine-cytokine receptor pathways. Additionally, the microRNAs targeting the 4 aforementioned genes were predicted.

Conclusion

The present study is the first to provide fresh insights into the potential therapeutic targets of key genes, namely CD19, CACNA1G, and ALDH18A1, differentially expressed in end-stage osteoarthritis and their relevant signaling pathways and interactive microRNAs.

Technical Report: A Comprehensive Comparison between Different Quantification Versions of Nightingale Health's 1H-NMR Metabolomics Platform

1H-NMR metabolomics data is increasingly used to track health and disease. Nightingale Health, a major supplier of 1H-NMR metabolomics, has recently updated the quantification strategy to further align with clinical standards. Such updates, however, might influence backward replicability, particularly affecting studies with repeated measures. Using data from BBMRI-NL consortium (~28,000 samples from 28 cohorts), we compared Nightingale data, originally released in 2014 and 2016, with a re-quantified version released in 2020, of which both versions were based on the same NMR spectra. Apart from two discontinued and twenty-three new analytes, we generally observe a high concordance between quantification versions with 73 out of 222 (33%) analytes showing a mean ρ > 0.9 across all cohorts. Conversely, five analytes consistently showed lower Spearman's correlations (ρ < 0.7) between versions, namely acetoacetate, LDL-L, saturated fatty acids, S-HDL-C, and sphingomyelins. Furthermore, previously trained multi-analyte scores, such as MetaboAge or MetaboHealth, might be particularly sensitive to platform changes. Whereas MetaboHealth replicated well, the MetaboAge score had to be retrained due to use of discontinued analytes. Notably, both scores in the re-quantified data recapitulated mortality associations observed previously. Concluding, we urge caution in utilizing different platform versions to avoid mixing analytes, having different units, or simply being discontinued.

Identification of candidate genes and chemicals associated with osteoarthritis by transcriptome-wide association study and chemical-gene interaction analysis

BACKGROUND

Osteoarthritis (OA) is a common degenerative joint disease and causes chronic pain and disability to the elderly. Several risk factors are involved, such as aging, obesity, genetic susceptibility, and environmental factors. We conducted a transcriptome-wide association study (TWAS) and chemical-related gene set enrichment analysis (CGSEA) to investigate the susceptibility genes and environmental factors.

METHODS

TWAS analysis was conducted to identify the susceptibility genes by integrating the summary-level genome-wide association study data of knee OA (KOA) and hip OA (HOA) with the precomputed expression weights from the Genotype-Tissue Expression Project (Version 8). The FUSION software was used for both single-tissue and cross-tissue TWAS, which were combined using an aggregate Cauchy association test. The biological function and pathways of the TWAS genes were explored using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases, and the human cartilage mRNA expression profiles were utilized to validate the TWAS genes. CGSEA analysis was performed to scan the OA-associated chemicals by integrating the TWAS results with the chemical-related gene sets.

RESULTS

There were 44 and 93 unique TWAS genes identified in 7 and 11 chromosomes for KOA and HOA, respectively, fourteen and four of which showed significantly differential expression in the mRNA profiles, such as CRHR1, LTBP1, WWP2, LMX1B, and PTHLH. OA-related pathways were found in the KEGG and GO analysis, such as TGF-beta signaling pathway, MAPK signaling pathway, hyaluronan metabolic process, and chondrocyte differentiation. Forty-five OA-associated chemicals were identified, including quercetin, bisphenol A, and cadmium chloride.

CONCLUSIONS

Several candidate OA-associated genes and chemicals were identified through TWAS and CGSEA analysis, which expanded our understanding of the relationship between genes, chemicals, and their impact on OA.

Exposure of Tissue-Engineered Cartilage Analogs to Synovial Fluid Hematoma After Ankle Fracture Is Associated With Chondrocyte Death and Altered Cartilage Maintenance Gene Expression

BACKGROUND

The first stage of fracture healing consists of hematoma formation with recruitment of proinflammatory cytokines and matrix metalloproteinases. Unfortunately, when there is an intra-articular fracture, these inflammatory mediators are not retained at the fracture site, but instead, envelop the healthy cartilage of the entire joint via the synovial fluid fracture hematoma (SFFH). These inflammatory cytokines and matrix metalloproteinases are known factors in the progression of osteoarthritis and rheumatoid arthritis. Despite the known inflammatory contents of the SFFH, little research has been done on the effects of the SFFH on healthy cartilage with regard to cell death and alteration in gene expression that could lead to posttraumatic osteoarthritis (PTOA).

METHODS

SFFH was collected from 12 patients with intraarticular ankle fracture at the time of surgery. Separately, C20A4 immortalized human chondrocytes were 3-dimensionally cultured to create scaffold-free cartilage tissue analogs (CTAs) to simulate healthy cartilage. Experimental CTAs (n = 12) were exposed to 100% SFFH for 3 days, washed, and transferred to complete media for 3 days. Control CTAs (n = 12) were simultaneously cultured in complete medium without exposure to SFFH. Subsequently, CTAs were harvested and underwent biochemical, histological, and gene expression analysis.

RESULTS

Exposure of CTAs to ankle SFFH for 3 days significantly decreased chondrocyte viability by 34% (P = .027). Gene expression of both COL2A1 and SOX9 were significantly decreased after exposure to SFFH (P = .012 and P = .0013 respectively), while there was no difference in COL1A1, RUNX2, and MMP13 gene expression. Quantitative analysis of Picrosirius red staining demonstrated increased collagen I deposition with poor ultrastructural organization in SFFH-exposed CTAs.

CONCLUSION

Exposure of an organoid model of healthy cartilage tissue to SFFH after intraarticular ankle fracture resulted in decreased chondrocyte viability, decreased expression of genes regulating normal chondrocyte phenotype, and altered matrix ultrastructure indicating differentiation toward an osteoarthritis phenotype.

CLINICAL RELEVANCE

The majority of ankle fracture open reduction and internal fixation does not occur immediately after fracture. In fact, typically these fractures are treated several days to weeks later in order to let the swelling subside. This means that the healthy innocent bystander cartilage not involved in the fracture is exposed to SFFH during this time. In this study, the SFFH caused decreased chondrocyte viability and specific altered gene expression that might have the potential to induce osteoarthritis. These data suggest that early intervention after intraarticular ankle fracture could possibly mitigate progression toward PTOA.

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

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

Inhibition of KDM7A/B histone demethylases restores H3K79 methylation and protects against osteoarthritis

OBJECTIVES

In osteoarthritis, methylation of lysine 79 on histone H3 (H3K79me), a protective epigenetic mechanism, is reduced. Histone methylation levels are dynamically regulated by histone methyltransferases and demethylases. Here, we aimed to identify which histone demethylases regulate H3K79me in cartilage and investigate whether their targeting protects against osteoarthritis.

METHODS

We determined histone demethylase expression in human non-osteoarthritis and osteoarthritis cartilage using qPCR. The role of histone demethylase families and subfamilies on H3K79me was interrogated by treatment of human C28/I2 chondrocytes with pharmacological inhibitors, followed by western blot and immunofluorescence. We performed C28/I2 micromasses to evaluate effects on glycosaminoglycans by Alcian blue staining. Changes in H3K79me after destabilisation of the medial meniscus (DMM) in mice were determined by immunohistochemistry. Daminozide, a KDM2/7 subfamily inhibitor, was intra-articularly injected in mice upon DMM. Histone demethylases targeted by daminozide were individually silenced in chondrocytes to dissect their role on H3K79me and osteoarthritis.

RESULTS

We documented the expression signature of histone demethylases in human non-osteoarthritis and osteoarthritis articular cartilage. Inhibition of Jumonji-C demethylase family increased H3K79me in human chondrocytes. Blockade of KDM2/7 histone demethylases with daminozide increased H3K79me and glycosaminoglycans. In mouse articular cartilage, H3K79me decayed rapidly upon induction of joint injury. Early and sustained intra-articular treatment with daminozide enhanced H3K79me and exerted protective effects in mice upon DMM. Individual silencing of KDM7A/B demethylases in human chondrocytes demonstrated that KDM7A/B mediate protective effects of daminozide on H3K79me and osteoarthritis.

CONCLUSION

Targeting KDM7A/B histone demethylases could be an attractive strategy to protect joints against osteoarthritis.

Mechanotransduction pathways in articular chondrocytes and the emerging role of estrogen receptor-α

In the synovial joint, mechanical force creates an important signal that influences chondrocyte behavior. The conversion of mechanical signals into biochemical cues relies on different elements in mechanotransduction pathways and culminates in changes in chondrocyte phenotype and extracellular matrix composition/structure. Recently, several mechanosensors, the first responders to mechanical force, have been discovered. However, we still have limited knowledge about the downstream molecules that enact alterations in the gene expression profile during mechanotransduction signaling. Recently, estrogen receptor α (ERα) has been shown to modulate the chondrocyte response to mechanical loading through a ligand-independent mechanism, in line with previous research showing that ERα exerts important mechanotransduction effects on other cell types, such as osteoblasts. In consideration of these recent discoveries, the goal of this review is to position ERα into the mechanotransduction pathways known to date. Specifically, we first summarize our most recent understanding of the mechanotransduction pathways in chondrocytes on the basis of three categories of actors, namely mechanosensors, mechanotransducers, and mechanoimpactors. Then, the specific roles played by ERα in mediating the chondrocyte response to mechanical loading are discussed, and the potential interactions of ERα with other molecules in mechanotransduction pathways are explored. Finally, we propose several future research directions that may advance our understanding of the roles played by ERα in mediating biomechanical cues under physiological and pathological conditions.

Lysosomal control of senescence and inflammation through cholesterol partitioning

Whereas cholesterol is vital for cell growth, proliferation, and remodeling, dysregulation of cholesterol metabolism is associated with multiple age-related pathologies. Here we show that senescent cells accumulate cholesterol in lysosomes to maintain the senescence-associated secretory phenotype (SASP). We find that induction of cellular senescence by diverse triggers enhances cellular cholesterol metabolism. Senescence is associated with the upregulation of the cholesterol exporter ABCA1, which is rerouted to the lysosome, where it moonlights as a cholesterol importer. Lysosomal cholesterol accumulation results in the formation of cholesterol-rich microdomains on the lysosomal limiting membrane enriched with the mammalian target of rapamycin complex 1 (mTORC1) scaffolding complex, thereby sustaining mTORC1 activity to support the SASP. We further show that pharmacological modulation of lysosomal cholesterol partitioning alters senescence-associated inflammation and in vivo senescence during osteoarthritis progression in male mice. Our study reveals a potential unifying theme for the role of cholesterol in the aging process through the regulation of senescence-associated inflammation.

Identification and functional characterization of imbalanced osteoarthritis-associated fibronectin splice variants

OBJECTIVE

To identify FN1 transcripts associated with OA pathophysiology and investigate the downstream effects of modulating FN1 expression and relative transcript ratio.

METHODS

FN1 transcriptomic data was obtained from our previously assessed RNA-seq dataset of lesioned and preserved OA cartilage samples from the Research osteoArthritis Articular Cartilage (RAAK) study. Differential transcript expression analysis was performed on all 27 FN1 transcripts annotated in the Ensembl database. Human primary chondrocytes were transduced with lentiviral particles containing short hairpin RNA (shRNA) targeting full-length FN1 transcripts or non-targeting shRNA. Subsequently, matrix deposition was induced in our 3D in vitro neo-cartilage model. Effects of changes in the FN1 transcript ratio on sulphated glycosaminoglycan (sGAG) deposition were investigated by Alcian blue staining and dimethylmethylene blue assay. Moreover, gene expression levels of 17 cartilage-relevant markers were determined by reverse transcription quantitative polymerase chain reaction.

RESULTS

We identified 16 FN1 transcripts differentially expressed between lesioned and preserved cartilage. FN1-208, encoding migration-stimulating factor, was the most significantly differentially expressed protein coding transcript. Downregulation of full-length FN1 and a concomitant increased FN1-208 ratio resulted in decreased sGAG deposition as well as decreased ACAN and COL2A1 and increased ADAMTS-5, ITGB1 and ITGB5 gene expression levels.

CONCLUSION

We show that full-length FN1 downregulation and concomitant relative FN1-208 upregulation was unbeneficial for deposition of cartilage matrix, likely due to decreased availability of the classical RGD (Arg-Gly-Asp) integrin-binding site of fibronectin.

WWP2 confers risk to osteoarthritis by affecting cartilage matrix deposition via hypoxia associated genes

OBJECTIVE

To explore the co-expression network of the osteoarthritis (OA) risk gene WWP2 in articular cartilage and study cartilage characteristics when mimicking the effect of OA risk allele rs1052429-A on WWP2 expression in a human 3D in vitro model of cartilage.

METHOD

Co-expression behavior of WWP2 with genes expressed in lesioned OA articular cartilage (N = 35 samples) was explored. By applying lentiviral particle mediated WWP2 upregulation in 3D in vitro pellet cultures of human primary chondrocytes (N = 8 donors) the effects of upregulation on cartilage matrix deposition was evaluated. Finally, we transfected primary chondrocytes with miR-140 mimics to evaluate whether miR-140 and WWP2 are involved in similar pathways.

RESULTS

Upon performing Spearman correlations in lesioned OA cartilage, 98 highly correlating genes (|ρ| > 0.7) were identified. Among these genes, we identified GJA1, GDF10, STC2, WDR1, and WNK4. Subsequent upregulation of WWP2 on 3D chondrocyte pellet cultures resulted in a decreased expression of COL2A1 and ACAN and an increase in EPAS1 expression. Additionally, we observed a decreased expression of GDF10, STC2, and GJA1. Proteomics analysis identified 42 proteins being differentially expressed with WWP2 upregulation, which were enriched for ubiquitin conjugating enzyme activity. Finally, upregulation of miR-140 in 2D chondrocytes resulted in significant upregulation of WWP2 and WDR1.

CONCLUSIONS

Mimicking the effect of OA risk allele rs1052429-A on WWP2 expression initiates detrimental processes in the cartilage shown by a response in hypoxia associated genes EPAS1, GDF10, and GJA1 and a decrease in anabolic markers, COL2A1 and ACAN.

Whole-genome RNA sequencing identifies distinct transcriptomic profiles in impingement cartilage between patients with femoroacetabular impingement and hip osteoarthritis

Femoroacetabular impingement (FAI) has a strong clinical association with the development of hip osteoarthritis (OA); however, the pathobiological mechanisms underlying the transition from focal impingement to global joint degeneration remain poorly understood. The purpose of this study is to use whole-genome RNA sequencing to identify and subsequently validate differentially expressed genes (DEGs) in femoral head articular cartilage samples from patients with FAI and hip OA secondary to FAI. Thirty-seven patients were included in the study with whole-genome RNA sequencing performed on 10 gender-matched patients in the FAI and OA cohorts and the remaining specimens were used for validation analyses. We identified a total of 3531 DEGs between the FAI and OA cohorts with multiple targets for genes implicated in canonical OA pathways. Quantitative reverse transcription-polymerase chain reaction validation confirmed increased expression of FGF18 and WNT16 in the FAI samples, while there was increased expression of MMP13 and ADAMTS4 in the OA samples. Expression levels of FGF18 and WNT16 were also higher in FAI samples with mild cartilage damage compared to FAI samples with severe cartilage damage or OA cartilage. Our study further expands the knowledge regarding distinct genetic reprogramming in the cartilage between FAI and hip OA patients. We independently validated the results of the sequencing analysis and found increased expression of anabolic markers in patients with FAI and minimal histologic cartilage damage, suggesting that anabolic signaling may be increased in early FAI with a transition to catabolic and inflammatory gene expression as FAI progresses towards more severe hip OA. Clinical significance:Cam-type FAI has a strong clinical association with hip OA; however, the cellular pathophysiology of disease progression remains poorly understood. Several previous studies have demonstrated increased expression of inflammatory markers in FAI cartilage samples, suggesting the involvement of these inflammatory pathways in the disease progression. Our study further expands the knowledge regarding distinct genetic reprogramming in the cartilage between FAI and hip OA patients. In addition to differences in inflammatory gene expression, we also identified differential expression in multiple pathways involved in hip OA progression.

Angiopoietin-like 3-derivative LNA043 for cartilage regeneration in osteoarthritis: a randomized phase 1 trial

Osteoarthritis (OA) is a common, debilitating, chronic disease with no disease-modifying drug approved to date. We discovered LNA043-a derivative of angiopoietin-like 3 (ANGPTL3)-as a potent chondrogenesis inducer using a phenotypic screen with human mesenchymal stem cells. We show that LNA043 promotes chondrogenesis and cartilage matrix synthesis in vitro and regenerates hyaline articular cartilage in preclinical OA and cartilage injury models in vivo. LNA043 exerts at least part of these effects through binding to the fibronectin receptor, integrin α5β1 on mesenchymal stem cells and chondrocytes. In a first-in-human (phase 1), randomized, double-blinded, placebo-controlled, single ascending dose, single-center trial ( NCT02491281 ; sponsored by Novartis Pharmaceuticals), 28 patients with knee OA were injected intra-articularly with LNA043 or placebo (3:1 ratio) either 2 h, 7 d or 21 d before total knee replacement. LNA043 met its primary safety endpoint and showed short serum pharmacokinetics, cartilage penetration and a lack of immunogenicity (secondary endpoints). Post-hoc transcriptomics profiling of cartilage revealed that a single LNA043 injection reverses the OA transcriptome signature over at least 21 d, inducing the expression of hyaline cartilage matrix components and anabolic signaling pathways, while suppressing mediators of OA progression. LNA043 is a novel disease-modifying OA drug candidate that is currently in a phase 2b trial ( NCT04864392 ) in patients with knee OA.

Characterizing the secretome of licensed hiPSC-derived MSCs

Although mesenchymal stromal cells (MSCs) from primary tissues have been successfully applied in the clinic, their expansion capabilities are limited and results are variable. MSCs derived from human-induced pluripotent stem cells (hiMSCs) are expected to overcome these limitations and serve as a reproducible and sustainable cell source. We have explored characteristics and therapeutic potential of hiMSCs in comparison to hBMSCs. RNA sequencing confirmed high resemblance, with average Pearson correlation of 0.88 and Jaccard similarity index of 0.99, and similar to hBMSCs the hiMSCs released extracellular vesicles with in vitro immunomodulatory properties. Potency assay with TNFα and IFNγ demonstrated an increase in well-known immunomodulatory genes such as IDO1, CXCL8/IL8, and HLA-DRA which was also highlighted by enhanced secretion in the media. Notably, expression of 125 genes increased more than 1000-fold. These genes were predicted to be regulated by NFΚB signaling, known to play a central role in immune response. Altogether, our data qualify hiMSCs as a promising source for cell therapy and/or cell-based therapeutic products. Additionally, the herewith generated database will add to our understanding of the mode of action of regenerative cell-based therapies and could be used to identify relevant potency markers.

Allelic expression imbalance in articular cartilage and subchondral bone refined genome-wide association signals in osteoarthritis

OBJECTIVES

To present an unbiased approach to identify positional transcript single nucleotide polymorphisms (SNPs) of osteoarthritis (OA) risk loci by allelic expression imbalance (AEI) analyses using RNA sequencing of articular cartilage and subchondral bone from OA patients.

METHODS

RNA sequencing from 65 articular cartilage and 24 subchondral bone from OA patients was used for AEI analysis. AEI was determined for all genes present in the 100 regions reported by the GWAS catalog that were also expressed in cartilage or bone. The count fraction of the alternative allele (φ) was calculated for each heterozygous individual with the risk-SNP or with the SNP in linkage disequilibrium (LD) with it (r2 > 0.6). Furthermore, a meta-analysis was performed to generate a meta-φ (null hypothesis median φ = 0.49) and P-value for each SNP.

RESULTS

We identified 30 transcript SNPs (28 in cartilage and 2 in subchondral bone) subject to AEI in 29 genes. Notably, 10 transcript SNPs were located in genes not previously reported in the GWAS catalog, including two long intergenic non-coding RNAs (lincRNAs), MALAT1 (meta-φ = 0.54, FDR = 1.7x10-4) and ILF3-DT (meta-φ = 0.6, FDR = 1.75x10-5). Moreover, 12 drugs were interacting with 7 genes displaying AEI, of which 7 drugs have been already approved.

CONCLUSIONS

By prioritizing proxy transcript SNPs that mark AEI in cartilage and/or subchondral bone at loci harboring GWAS signals, we present an unbiased approach to identify the most likely functional OA risk-SNP and gene. We identified 10 new potential OA risk genes ready for further, translation towards underlying biological mechanisms.

Human Mesenchymal Stem Cell-Derived Miniature Joint System for Disease Modeling and Drug Testing

Diseases of the knee joint such as osteoarthritis (OA) affect all joint elements. An in vitro human cell-derived microphysiological system capable of simulating intraarticular tissue crosstalk is desirable for studying etiologies/pathogenesis of joint diseases and testing potential therapeutics. Herein, a human mesenchymal stem cell-derived miniature joint system (miniJoint) is generated, in which engineered osteochondral complex, synovial-like fibrous tissue, and adipose tissue are integrated into a microfluidics-enabled bioreactor. This novel design facilitates different tissues communicating while still maintaining their respective phenotypes. The miniJoint exhibits physiologically relevant changes when exposed to interleukin-1β mediated inflammation, which are similar to observations in joint diseases in humans. The potential of the miniJoint in predicting in vivo efficacy of drug treatment is confirmed by testing the "therapeutic effect" of the nonsteroidal anti-inflammatory drug, naproxen, as well as four other potential disease-modifying OA drugs. The data demonstrate that the miniJoint recapitulates complex tissue interactions, thus providing a robust organ chip model for the study of joint pathology and the development of novel therapeutic interventions.

Epigenetic Regulation of Chondrocytes and Subchondral Bone in Osteoarthritis

The aim of this review is to provide an updated review of the epigenetic factors involved in the onset and development of osteoarthritis (OA). OA is a prevalent degenerative joint disease characterized by chronic inflammation, ectopic bone formation within the joint, and physical and proteolytic cartilage degradation which result in chronic pain and loss of mobility. At present, no disease-modifying therapeutics exist for the prevention or treatment of the disease. Research has identified several OA risk factors including mechanical stressors, physical activity, obesity, traumatic joint injury, genetic predisposition, and age. Recently, there has been increased interest in identifying epigenetic factors involved in the pathogenesis of OA. In this review, we detail several of these epigenetic modifications with known functions in the onset and progression of the disease. We also review current therapeutics targeting aberrant epigenetic regulation as potential options for preventive or therapeutic treatment.

Mutation in the CCAL1 locus accounts for bidirectional process of human subchondral bone turnover and cartilage mineralization

OBJECTIVES

To study the mechanism by which the readthrough mutation in TNFRSF11B, encoding osteoprotegerin (OPG) with additional 19 amino acids at its C-terminus (OPG-XL), causes the characteristic bidirectional phenotype of subchondral bone turnover accompanied by cartilage mineralization in chondrocalcinosis patients.

METHODS

OPG-XL was studied by human induced pluripotent stem cells expressing OPG-XL and two isogenic CRISPR/Cas9-corrected controls in cartilage and bone organoids. Osteoclastogenesis was studied with monocytes from OPG-XL carriers and matched healthy controls followed by gene expression characterization. Dual energy X-ray absorptiometry scans and MRI analyses were used to characterize the phenotype of carriers and non-carriers of the mutation.

RESULTS

Human OPG-XL carriers relative to sex- and age-matched controls showed, after an initial delay, large active osteoclasts with high number of nuclei. By employing hiPSCs expressing OPG-XL and isogenic CRISPR/Cas9-corrected controls to established cartilage and bone organoids, we demonstrated that expression of OPG-XL resulted in excessive fibrosis in cartilage and high mineralization in bone accompanied by marked downregulation of MGP, encoding matrix Gla protein, and upregulation of DIO2, encoding type 2 deiodinase, gene expression, respectively.

CONCLUSIONS

The readthrough mutation at CCAL1 locus in TNFRSF11B identifies an unknown role for OPG-XL in subchondral bone turnover and cartilage mineralization in humans via DIO2 and MGP functions. Previously, OPG-XL was shown to affect binding between RANKL and heparan sulphate (HS) resulting in loss of immobilized OPG-XL. Therefore, effects may be triggered by deficiency in the immobilization of OPG-XL Since the characteristic bidirectional pathophysiology of articular cartilage calcification accompanied by low subchondral bone mineralization is also a hallmark of OA pathophysiology, our results are likely extrapolated to common arthropathies.

Inhibiting thyroid activation in aged human explants prevents mechanical induced detrimental signalling by mitigating metabolic processes

OBJECTIVES

To investigate whether the deiodinase inhibitor iopanoic acid (IOP) has chondroprotective properties, a mechanical stress induced model of human aged explants was used to test both repeated dosing and slow release of IOP.

METHODS

Human osteochondral explants subjected to injurious mechanical stress (65%MS) were treated with IOP or IOP encapsulated in poly lactic-co-glycolic acid-polyethylene glycol nanoparticles (NP-IOP). Changes to cartilage integrity and signalling were determined by Mankin scoring of histology, sulphated glycosaminoglycan (sGAG) release and expression levels of catabolic, anabolic and hypertrophic markers. Subsequently, on a subgroup of samples, RNA sequencing was performed on 65%MS (n = 14) and 65%MS+IOP (n = 7) treated cartilage to identify IOP's mode of action.

RESULTS

Damage from injurious mechanical stress was confirmed by increased cartilage surface damage in the Mankin score, increased sGAG release, and consistent upregulation of catabolic markers and downregulation of anabolic markers. IOP and, though less effective, NP-IOP treatment, reduced MMP13 and increased COL2A1 expression. In line with this, IOP and NP-IOP reduced cartilage surface damage induced by 65%MS, while only IOP reduced sGAG release from explants subjected to 65%MS. Lastly, differential expression analysis identified 12 genes in IOP's mode of action to be mainly involved in reducing metabolic processes (INSIG1, DHCR7, FADS1 and ACAT2) and proliferation and differentiation (CTGF, BMP5 and FOXM1).

CONCLUSION

Treatment with the deiodinase inhibitor IOP reduced detrimental changes of injurious mechanical stress. In addition, we identified that its mode of action was likely on metabolic processes, cell proliferation and differentiation.

Novel role of estrogen receptor-α on regulating chondrocyte phenotype and response to mechanical loading

OBJECTIVE

In knee cartilage from patients with osteoarthritis (OA), both preserved cartilage and damaged cartilage are observed. In this study, we aim to compare preserved with damaged cartilage to identify the molecule(s) that may be responsible for the mechanical loading-induced differences within cartilage degradation.

METHODS

Preserved and damaged cartilage were harvested from the same OA knee joint. RNA Sequencing was performed to examine the transcriptomic differences between preserved and damaged cartilage cells. Estrogen receptor-α (ERα) was identified, and its function of was tested through gene knockin and knockout. The role of ERα in mediating chondrocyte response to mechanical loading was examined via compression of chondrocyte-laded hydrogel in a strain-controlled manner. Findings from the studies on human samples were verified in animal models.

RESULTS

Level of estrogen receptor α (ERα) was significantly reduced in damaged cartilage compared to preserved cartilage, which were observed in both human and mice samples. Knockdown of ESR1, the gene encoding ERα, resulted in an upregulation of senescence- and OA-relevant markers in chondrocytes. Conversely, knockin of ESR1 partially reversed the osteoarthritic and senescent phenotype of OA chondrocytes. Using a three-dimensional (3D) culture model, we demonstrated that mechanical overload significantly suppressed ERα level in chondrocytes with concomitant upregulation of osteoarthritic phenotype. When ESR1 expression was suppressed, mechanical loading enhanced hypertrophic and osteogenic transition.

CONCLUSION

Our study demonstrates a new estrogen-independent role of ERα in mediating chondrocyte phenotype and its response to mechanical loading, and suggests that enhancing ERα level may represent a new method to treat osteoarthritis.

Identification and development of the novel 7-genes diagnostic signature by integrating multi cohorts based on osteoarthritis

Background

A chronic progressive degenerative joint disease, such as osteoarthritis (OA) is positively related to age. The medical economy is facing a major burden, because of the high disability rate seen in patients with OA. Therefore, to prevent and treat OA, exploring the diagnostic biomarkers of OA will be of great significance.

Methods

Differentially expressed genes (DEGs) were obtained from the Gene Expression Omnibus database using the RobustRankAggreg R package, and a protein–protein interaction network was constructed. The module was obtained from Cytoscape, and the four algorithms of degree, MNC, closeness, and MCC in CytoHubba were used to identify the hub genes. A diagnostic model was constructed using Support Vector Machines (SVM), and the ability of the model to predict was evaluated by other cohorts.

Results

From normal and OA samples, 136 DEGs were identified, out of which 45 were downregulated in the normal group and 91 were upregulated in the OA group. These genes were associated with the extracellular matrix-receptor interactions, the PI3K-Akt signaling pathway, and the protein digestion and absorption pathway, as per a functional enrichment analysis. Finally, we identified the 7 hub genes (COL6A3, COL1A2, COL1A1, MMP2, COL3A1, POST, and FN1). These genes have important roles and are widely involved in the immune response, apoptosis, inflammation, and bone development. These 7 genes were used to construct a diagnostic model by SVM, and it performed well in different cohorts. Additionally, we verified the methylation expression of these hub genes.

Conclusions

The 7-genes signature can be used for the diagnosis of OA and can provide new ideas in the clinical decision-making for patients with OA.

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

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

CRLF1 and CLCF1 in Development, Health and Disease

Cytokines and their receptors have a vital function in regulating various processes such as immune function, inflammation, haematopoiesis, cell growth and differentiation. The interaction between a cytokine and its specific receptor triggers intracellular signalling cascades that lead to altered gene expression in the target cell and consequent changes in its proliferation, differentiation, or activation. In this review, we highlight the role of the soluble type I cytokine receptor CRLF1 (cytokine receptor-like factor-1) and the Interleukin (IL)-6 cytokine CLCF1 (cardiotrophin-like cytokine factor 1) during development in physiological and pathological conditions with particular emphasis on Crisponi/cold-induced sweating syndrome (CS/CISS) and discuss new insights, challenges and possibilities arising from recent studies.

Proteomics profiling of human synovial fluid suggests increased protein interplay in early-osteoarthritis (OA) that is lost in late-stage OA

The underlying molecular mechanisms in osteoarthritis (OA) development are largely unknown. This study explores the proteome and the pairwise interplay of proteins in synovial fluid from patients with late-stage knee OA (arthroplasty), early knee OA (arthroscopy due to degenerative meniscal tear), and from deceased controls without knee OA. Synovial fluid samples were analyzed using state-of-the-art mass spectrometry with data-independent acquisition. The differential expression of the proteins detected was clustered and evaluated with data mining strategies and a multilevel model. Group-specific slopes of associations were estimated between expressions of each pair of identified proteins to assess the co-expression (i.e., interplay) between the proteins in each group. More proteins were increased in early-OA versus controls than late-stage OA versus controls. For most of these proteins, the fold changes between late-stage OA versus controls and early-stage OA versus controls were remarkably similar suggesting potential involvement in the OA process. Further, for the first time, this study illustrated distinct patterns in protein co-expression suggesting that the interplay between the protein machinery is increased in early-OA and lost in late-stage OA. Further efforts should focus on earlier stages of the disease than previously considered.

•

Synovial fluid proteomics study of different stages of osteoarthritis (OA).

•

Higher catabolic activity is found in both early- and late-stage OA.

•

Imbalance of the metabolic homeostasis in late-stage OA.

•

Understanding early-stage OA may lead to finding better effective therapies.

Identification and development of a novel 5-gene diagnostic model based on immune infiltration analysis of osteoarthritis

Background

Osteoarthritis (OA), which is due to the progressive loss and degeneration of articular cartilage, is the leading cause of disability worldwide. Therefore, it is of great significance to explore OA biomarkers for the prevention, diagnosis, and treatment of OA.

Methods and materials

The GSE129147, GSE57218, GSE51588, GSE117999, and GSE98918 datasets with normal and OA samples were downloaded from the Gene Expression Omnibus (GEO) database. The GSE117999 and GSE98918 datasets were integrated, and immune infiltration was evaluated. The differentially expressed genes (DEGs) were analyzed using the limma package in R, and weighted gene co-expression network analysis (WGCNA) was used to explore the co-expression genes and co-expression modules. The co-expression module genes were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. A protein–protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database, and hub genes were identified by the degree, MNC, closeness, and MCC algorithms. The hub genes were used to construct a diagnostic model based on support vector machines.

Results

The Immune Score in the OA samples was significantly higher than in the normal samples, and a total of 2313 DEGs were identified. Through WGCNA, we found that the yellow module was significantly positively correlated with the OA samples and Immune Score and negatively correlated with the normal samples. The 142 DEGs of the yellow module were related to biological processes such as regulation of inflammatory response, positive regulation of inflammatory response, blood vessel morphogenesis, endothelial cell migration, and humoral immune response. The intersections of the genes obtained by the 4 algorithms resulted in 5 final hub genes, and the diagnostic model constructed with these 5 genes showed good performance in the training and validation cohorts.

Conclusions

The 5-gene diagnostic model can be used to diagnose OA and guide clinical decision-making.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03183-9.

Scientific Developments and Clinical Applications Utilizing Chondrons and Chondrocytes with Matrix for Cartilage Repair

Injuries to articular cartilage of the knee are increasingly common. The operative management of these focal chondral lesions continues to be problematic for the treating orthopedic surgeon secondary to the limited regenerative capacity of articular cartilage. The pericellular matrix (PCM) is a specialized, thin layer of the extracellular matrix that immediately surrounds chondrocytes forming a unit together called the chondron. The advancements in our knowledge base with regard to the PCM/chondrons as well as interterritorial matrix has permeated and led to advancements in product development in conjunction with minced cartilage, marrow stimulation, osteochondral allograft, and autologous chondrocyte implantation (ACI). This review intends to summarize recent progress in chondrocytes with matrix research, with an emphasis on the role the PCM/extracellular matrix (ECM) plays for favorable chondrogenic gene expression, as a barrier/filtration unit, and in osteoarthritis. The bulk of the review describes cutting-edge and evolving clinical developments and discuss these developments in light of underlying basic science applications. Clinical applications of chondrocytes with matrix science include Reveille Cartilage Processor, Cartiform, and ACI with Spherox (which was recently recommended for the treatment of grade III or IV articular cartilage defects over 2 cm² by the National Institute of Health and Care Excellence [NICE] in the United Kingdom). The current article presents a comprehensive overview of both the basic science and clinical results of these next-generation cartilage repair techniques by focusing specifically on the scientific evolution in each category as it pertains with underlying chondrocytes with matrix theory.

A dysfunctional TRPV4-GSK3β pathway prevents osteoarthritic chondrocytes from sensing changes in extracellular matrix viscoelasticity

Changes in the composition and viscoelasticity of the extracellular matrix in load-bearing cartilage influence the proliferation and phenotypes of chondrocytes, and are associated with osteoarthritis. However, the underlying molecular mechanism is unknown. Here we show that the viscoelasticity of alginate hydrogels regulates cellular volume in healthy human chondrocytes (with faster stress relaxation allowing cell expansion and slower stress relaxation restricting it) but not in osteoarthritic chondrocytes. Cellular volume regulation in healthy chondrocytes was associated with changes in anabolic gene expression, in the secretion of multiple pro-inflammatory cytokines, and in the modulation of intracellular calcium regulated by the ion-channel protein transient receptor potential cation channel subfamily V member 4 (TRPV4), which controls the phosphorylation of glycogen synthase kinase 3β (GSK3β), an enzyme with pleiotropic effects in osteoarthritis. A dysfunctional TRPV4-GSK3β pathway in osteoarthritic chondrocytes rendered the cells unable to respond to environmental changes in viscoelasticity. Our findings suggest strategies for restoring chondrocyte homeostasis in osteoarthritis.

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

OBJECTIVE

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

DESIGN

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

RESULTS

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

CONCLUSION

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

High-impact FN1 mutation decreases chondrogenic potential and affects cartilage deposition via decreased binding to collagen type II

Osteoarthritis is the most prevalent joint disease worldwide, yet progress in development of effective disease-modifying treatments is slow because of lack of insight into the underlying disease pathways. Therefore, we aimed to identify the causal pathogenic mutation in an early-onset osteoarthritis family, followed by functional studies in human induced pluripotent stem cells (hiPSCs) in an in vitro organoid cartilage model. We demonstrated that the identified causal missense mutation in the gelatin-binding domain of the extracellular matrix protein fibronectin resulted in significant decreased binding capacity to collagen type II. Further analyses of formed hiPSC-derived neo-cartilage tissue highlighted that mutated fibronectin affected chondrogenic capacity and propensity to a procatabolic osteoarthritic state. Together, we demonstrate that binding of fibronectin to collagen type II is crucial for fibronectin downstream gene expression of chondrocytes. We advocate that effective treatment development should focus on restoring or maintaining proper binding between fibronectin and collagen type II.

Application of Induced Pluripotent Stem Cells for Disease Modeling and 3D Model Construction: Focus on Osteoarthritis

Since their discovery in 2006, induced pluripotent stem cells (iPSCs) have shown promising potential, specifically because of their accessibility and plasticity. Hence, the clinical applicability of iPSCs was investigated in various fields of research. However, only a few iPSC studies pertaining to osteoarthritis (OA) have been performed so far, despite the high prevalence rate of degenerative joint disease. In this review, we discuss some of the most recent applications of iPSCs in disease modeling and the construction of 3D models in various fields, specifically focusing on osteoarthritis and OA-related conditions. Notably, we comprehensively reviewed the successful results of iPSC-derived disease models in recapitulating OA phenotypes for both OA and early-onset OA to encompass their broad etiology. Moreover, the latest publications with protocols that have used iPSCs to construct 3D models in recapitulating various conditions, particularly the OA environment, were further discussed. With the overall optimistic results seen in both fields, iPSCs are expected to be more widely used for OA disease modeling and 3D model construction, which could further expand OA drug screening, risk assessment, and therapeutic capabilities.

Long non-coding RNA expression profiling of subchondral bone reveals AC005165.1 modifying FRZB expression during osteoarthritis

Objective

To gain insight in the expression profile of long non-coding RNAs (lncRNAs) in OA subchondral bone.

Methods

RNA sequencing data of macroscopically preserved and lesioned OA subchondral bone of patients that underwent joint replacement surgery due to OA (N = 22 pairs; 5 hips, 17 knees, Research osteoArthrits Articular Tissue (RAAK study) was run through an in-house pipeline to detect expression of lncRNAs. Differential expression analysis between preserved and lesioned bone was performed. Spearman correlations were calculated between differentially expressed lncRNAs and differentially expressed mRNAs identified previously in the same samples. Primary osteogenic cells were transfected with locked nucleic acid (LNA) GapmeRs targeting AC005165.1 lncRNA, to functionally investigate its potential mRNA targets.

Results

In total, 2816 lncRNAs were well-expressed in subchondral bone and we identified 233 lncRNAs exclusively expressed in knee and 307 lncRNAs exclusively in hip. Differential expression analysis, using all samples (N = 22 pairs; 5 hips, 17 knees), resulted in 21 differentially expressed lncRNAs [false discovery rate (FDR) < 0.05, fold change (FC) range 1.19–7.39], including long intergenic non-protein coding RNA (LINC) 1411 (LINC01411, FC = 7.39, FDR = 2.20 × 10⁻⁸), AC005165.1 (FC = 0.44, FDR = 2.37 × 10⁻⁶) and empty spiracles homeobox 2 opposite strand RNA (EMX2OS, FC = 0.41, FDR = 7.64 × 10⁻³). Among the differentially expressed lncRNAs, five were also differentially expressed in articular cartilage, including AC005165.1, showing similar direction of effect. Downregulation of AC005165.1 in primary osteogenic cells resulted in consistent downregulation of highly correlated frizzled related protein (FRZB).

Conclusion

The current study identified a novel lncRNA, AC005165.1, being dysregulated in OA articular cartilage and subchondral bone. Downregulation of AC005165.1 caused a decreased expression of OA risk gene FRZB, an important member of the wnt pathway, suggesting that AC005165.1 could be an attractive potential therapeutic target with effects in articular cartilage and subchondral bone.

An evolutionary perspective of DNA methylation patterns in skeletal tissues using a baboon model of osteoarthritis

Epigenetic factors, such as DNA methylation, play an influential role in the development of the degenerative joint disease osteoarthritis (OA). These molecular mechanisms have been heavily studied in humans, and although OA affects several other animals in addition to humans, few efforts have taken an evolutionary perspective. This study explores the evolution of OA epigenetics by assessing the relationship between DNA methylation variation and knee OA development in baboons (Papio spp.) and by comparing these findings to human OA epigenetic associations. Genome-wide DNA methylation patterns were identified in bone and cartilage of the right distal femora from 56 pedigreed, adult baboons (28 with and 28 without knee OA) using the Illumina Infinium MethylationEPIC BeadChip. Several significantly differentially methylated positions (DMPs) and regions were found between tissue types. Substantial OA-related differential methylation was also identified in cartilage, but not in bone, suggesting that cartilage epigenetics may be more influential in OA than bone epigenetics. Additionally, some genes containing OA-related DMPs overlap with and display methylation patterns similar to those previously identified in human OA, revealing a mixture of evolutionarily conserved and divergent OA-related methylation patterns in primates. Overall, these findings reinforce the current etiological perspectives of OA and enhance our evolutionary understanding of epigenetic mechanisms associated with OA. This study further establishes baboons as a valuable nonhuman primate model of OA, and continued investigations in baboons will help to disentangle the molecular mechanisms contributing to OA and their evolutionary histories.

Circulating MicroRNAs Highly Correlate to Expression of Cartilage Genes Potentially Reflecting OA Susceptibility-Towards Identification of Applicable Early OA Biomarkers

Objective: To identify and validate circulating micro RNAs (miRNAs) that mark gene expression changes in articular cartilage early in osteoarthritis (OA) pathophysiology process. Methods: Within the ongoing RAAK study, human preserved OA cartilage and plasma (N = 22 paired samples) was collected for RNA sequencing (respectively mRNA and miRNA). Spearman correlation was determined for 114 cartilage genes consistently and significantly differentially expressed early in osteoarthritis and 384 plasma miRNAs. Subsequently, the minimal number of circulating miRNAs serving to discriminate between progressors and non-progressors was assessed by regression analysis and area under receiver operating curves (AUC) was calculated with progression data and plasma miRNA sequencing from the GARP study (N = 71). Results: We identified strong correlations (ρ ≥ |0.7|) among expression levels of 34 unique plasma miRNAs and 21 genes, including 4 genes that correlated with multiple miRNAs. The strongest correlation was between let-7d-5p and EGFLAM (ρ = −0.75, P = 6.9 × 10⁻⁵). Regression analysis of the 34 miRNAs resulted in a set of 7 miRNAs that, when applied to the GARP study, demonstrated clinically relevant predictive value with AUC > 0.8 for OA progression over 2 years and near-clinical value for progression over 5 years- (AUC = 0.8). Conclusions: We show that plasma miRNAs levels reflect gene expression levels in cartilage and can be exploited to represent ongoing pathophysiological processes in articular cartilage. We advocate that identified signature of 7 plasma miRNAs can contribute to direct further studies toward early biomarkers predictive for progression of osteoarthritis over 2 and 5 years.

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Chondrocytes from Osteoarthritis Patients Adopt Distinct Phenotypes in Response to Central TH1/TH2/TH17 Cytokines

Chronic low-grade inflammation plays a central role in the pathogenesis of osteoarthritis (OA), and several pro- and anti-inflammatory cytokines have been implicated to mediate and regulate this process. Out of these cytokines, particularly IFNγ, IL-1β, IL-4 and IL-17 are associated with different phenotypes of T helper (T_H) cells and macrophages, both examples of cells known for great phenotypic and functional heterogeneity. Chondrocytes also display various phenotypic changes during the course of arthritis. We set out to study the hypothesis of whether chondrocytes might adopt polarized phenotypes analogous to T_H cells and macrophages. We studied the effects of IFNγ, IL-1β, IL-4 and IL-17 on gene expression in OA chondrocytes with RNA-Seq. Chondrocytes were harvested from the cartilage of OA patients undergoing knee replacement surgery and then cultured with or without the cytokines for 24 h. Total RNA was isolated and sequenced, and GO (Gene Ontology) functional analysis was performed. We also separately investigated genes linked to OA in recent genome wide expression analysis (GWEA) studies. The expression of more than 2800 genes was significantly altered in chondrocytes treated with IL-1β [in the C(IL-1β) phenotype] with a fold change (FC) > 2.5 in either direction. These included a large number of genes associated with inflammation, cartilage degradation and attenuation of metabolic signaling. The profile of genes differentially affected by IFNγ (the C(IFNγ) phenotype) was relatively distinct from that of the C(IL-1β) phenotype and included several genes associated with antigen processing and presentation. The IL-17-induced C(IL-17) phenotype was characterized by the induction of a more limited set of proinflammatory factors compared to C(IL-1β) cells. The C(IL-4) phenotype induced by IL-4 displayed a differential expression of a rather small set of genes compared with control, primarily those associated with TGFβ signaling and the regulation of inflammation. In conclusion, our results show that OA chondrocytes can adopt diverse phenotypes partly analogously to T_H cells and macrophages. This phenotypic plasticity may play a role in the pathogenesis of arthritis and open new therapeutic avenues for the development of disease-modifying treatments for (osteo)arthritis.

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

Background

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

Methods

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

Results

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

Conclusions

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

Supplementary Information

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

Cripto favors chondrocyte hypertrophy via TGF-β SMAD1/5 signaling during development of osteoarthritis

Chondrocytes in mice developing osteoarthritis (OA) exhibit an aberrant response to the secreted cytokine transforming growth factor (TGF)‐β, consisting in a potentiation of intracellular signaling downstream of the transmembrane type I receptor kinase activin receptor‐like kinase (ALK)1 against canonical TGF‐β receptor ALK5‐mediated signaling. Unfortunately, the underlying mechanisms remain elusive. In order to identify novel druggable targets for OA, we aimed to investigate novel molecules regulating the ALK1/ALK5 balance in OA chondrocytes. We performed gene expression analysis of TGF‐β signaling modulators in joints from three different mouse models of OA and found an upregulated expression of the TGF‐β co‐receptor Cripto (Tdgf1), which was validated in murine and human cartilage OA samples at the protein level. In vitro and ex vivo, elevated expression of Cripto favors the hypertrophic differentiation of chondrocytes, eventually contributing to tissue calcification. Furthermore, we found that Cripto participates in a TGF‐β–ALK1–Cripto receptor complex in the plasma membrane, thereby inducing catabolic SMAD1/5 signaling in chondrocytes. In conclusion, we demonstrate that Cripto is expressed in OA and plays a functional role promoting chondrocyte hypertrophy, thereby becoming a novel potential therapeutic target in OA, for which there is no efficient cure or validated biomarker. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Green-lipped (greenshell™) mussel (Perna canaliculus) extract supplementation in treatment of osteoarthritis: a systematic review

OBJECTIVES

Intervention studies using New Zealand green-lipped or greenshell™ mussel (GSM) (Perna canaliculus) extract in osteoarthritis (OA) patients have shown effective pain relief. This systematic review summarises the efficacy of GSM extracts in the treatment of OA.

METHODS

A literature search of the three databases EMBASE, MEDLINE, and Scopus was performed to identify relevant articles published up to March 2020. Inclusion criteria were clinical trials published in English measuring the effect of supplementation of whole or a lipid extract from GSM on pain and mobility outcomes in OA patients.

RESULTS

A total of nine clinical trials were included in systematic review, from which five studies were considered appropriate for inclusion in a forest plot. Pooled results showed that GSM extracts (lipid extract or whole powder) provide moderate and clinically significant treatment effects on a visual analogue scale (VAS) pain score (effect size: - 0.46; 95% CI - 0.82 to - 0.10; p = 0.01). The whole GSM extract improved gastrointestinal symptoms in OA patients taking anti-inflammatory medications. The GSM extract was considered to be generally well tolerated in most of the studies.

CONCLUSION

The overall analysis showed that GSM provided moderate and clinically meaningful treatment effects on OA pain. However, the current evidence is limited by the number and quality of studies, and further larger and high-quality studies are needed to confirm the effectiveness and to identify the optimal GSM format. Nevertheless, it is worth considering using GSM extracts especially for patients seeking alternative pain relief treatments with fewer side effects compared to conventional treatment.

Characterization of dynamic changes in Matrix Gla Protein (MGP) gene expression as function of genetic risk alleles, osteoarthritis relevant stimuli, and the vitamin K inhibitor warfarin

OBJECTIVE

We here aimed to characterize changes of Matrix Gla Protein (MGP) expression in relation to its recently identified OA risk allele rs1800801-T in OA cartilage, subchondral bone and human ex vivo osteochondral explants subjected to OA related stimuli. Given that MGP function depends on vitamin K bioavailability, we studied the effect of frequently prescribed vitamin K antagonist warfarin.

METHODS

Differential (allelic) mRNA expression of MGP was analyzed using RNA-sequencing data of human OA cartilage and subchondral bone. Human osteochondral explants were used to study exposures to interleukin one beta (IL-1β; inflammation), triiodothyronine (T3; Hypertrophy), warfarin, or 65% mechanical stress (65%MS) as function of rs1800801 genotypes.

RESULTS

We confirmed that the MGP risk allele rs1800801-T was associated with lower expression and that MGP was significantly upregulated in lesioned as compared to preserved OA tissues, mainly in risk allele carriers, in both cartilage and subchondral bone. Moreover, MGP expression was downregulated in response to OA like triggers in cartilage and subchondral bone and this effect might be reduced in carriers of the rs1800801-T risk allele. Finally, warfarin treatment in cartilage increased COL10A1 and reduced SOX9 and MMP3 expression and in subchondral bone reduced COL1A1 and POSTN expression.

DISCUSSION & CONCLUSIONS

Our data highlights that the genetic risk allele lowers MGP expression and upon OA relevant triggers may hamper adequate dynamic changes in MGP expression, mainly in cartilage. The determined direct negative effect of warfarin on human explant cultures functionally underscores the previously found association between vitamin K deficiency and OA.

Cartilage from human-induced pluripotent stem cells: comparison with neo-cartilage from chondrocytes and bone marrow mesenchymal stromal cells

Cartilage has little intrinsic capacity for repair, so transplantation of exogenous cartilage cells is considered a realistic option for cartilage regeneration. We explored whether human-induced pluripotent stem cells (hiPSCs) could represent such unlimited cell sources for neo-cartilage comparable to human primary articular chondrocytes (hPACs) or human bone marrow-derived mesenchymal stromal cells (hBMSCs). For this, chondroprogenitor cells (hiCPCs) and hiPSC-derived mesenchymal stromal cells (hiMSCs) were generated from two independent hiPSC lines and characterized by morphology, flow cytometry, and differentiation potential. Chondrogenesis was compared to hBMSCs and hPACs by histology, immunohistochemistry, and RT-qPCR, while similarities were estimated based on Pearson correlations using a panel of 20 relevant genes. Our data show successful differentiations of hiPSC into hiMSCs and hiCPCs. Characteristic hBMSC markers were shared between hBMSCs and hiMSCs, with the exception of CD146 and CD45. However, neo-cartilage generated from hiMSCs showed low resemblances when compared to hBMSCs (53%) and hPACs (39%) characterized by lower collagen type 2 and higher collagen type 1 expression. Contrarily, hiCPC neo-cartilage generated neo-cartilage more similar to hPACs (65%), with stronger expression of matrix deposition markers. Our study shows that taking a stepwise approach to generate neo-cartilage from hiPSCs via chondroprogenitor cells results in strong similarities to neo-cartilage of hPACs within 3 weeks following chondrogenesis, making them a potential candidate for regenerative therapies. Contrarily, neo-cartilage deposited by hiMSCs seems more prone to hypertrophic characteristics compared to hPACs. We therefore compared chondrocytes derived from hiMSCs and hiCPCs with hPACs and hBMSCs to outline similarities and differences between their neo-cartilage and establish their potential suitability for regenerative medicine and disease modelling.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Single-Cell RNA-Seq Reveals Transcriptomic Heterogeneity and Post-Traumatic Osteoarthritis-Associated Early Molecular Changes in Mouse Articular Chondrocytes

Articular cartilage is a connective tissue lining the surfaces of synovial joints. When the cartilage severely wears down, it leads to osteoarthritis (OA), a debilitating disease that affects millions of people globally. The articular cartilage is composed of a dense extracellular matrix (ECM) with a sparse distribution of chondrocytes with varying morphology and potentially different functions. Elucidating the molecular and functional profiles of various chondrocyte subtypes and understanding the interplay between these chondrocyte subtypes and other cell types in the joint will greatly expand our understanding of joint biology and OA pathology. Although recent advances in high-throughput OMICS technologies have enabled molecular-level characterization of tissues and organs at an unprecedented resolution, thorough molecular profiling of articular chondrocytes has not yet been undertaken, which may be in part due to the technical difficulties in isolating chondrocytes from dense cartilage ECM. In this study, we profiled articular cartilage from healthy and injured mouse knee joints at a single-cell resolution and identified nine chondrocyte subtypes with distinct molecular profiles and injury-induced early molecular changes in these chondrocytes. We also compared mouse chondrocyte subpopulations to human chondrocytes and evaluated the extent of molecular similarity between mice and humans. This work expands our view of chondrocyte heterogeneity and rapid molecular changes in chondrocyte populations in response to joint trauma and highlights potential mechanisms that trigger cartilage degeneration.

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

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

•

Activated macrophages promote the secretion of FliI from chondrocytes

•

FliI acts as a DAMP-triggering molecule in cartilage

•

FliI promotes chondrocyte hypertrophy and cartilage catabolism

•

FliI represents attractive target for therapeutic intervention

Characterizing correlations among disease severity measures in osteochondral tissues from osteoarthritic knees

Osteoarthritis (OA) is a complex disease with biologic, biomechanical, and clinical heterogeneity among patients. Relationships among OA tissue metabolism, histopathology, and extracellular matrix (ECM) composition have not been well characterized. It was hypothesized that moderate (r = .4-.69) to strong (r > .7) correlations exist among these different measures of disease severity in osteochondral tissues from OA knees. Joint surfaces were obtained from patients (n = 6) undergoing total knee arthroplasty. Osteochondral explants (n = 136) were created and cultured for 3 days. Culture media were collected for biomarker analyses, and tissue was assessed for viability, histological scoring, and ECM composition. Correlations among media biomarker concentrations, histological scoring, ECM composition, and viability were determined using a Spearman correlation. GRO-α, IL-6, IL-8, and MCP-1 showed strong positive correlations to each other, and moderate positive correlations to NO, PGE2, and MMP-2. Total MMP activity, MMP-9, and MMP-13 had strong positive correlations to each other, and moderate positive correlations to MMP-1. MMP-2 had a moderate to strong positive correlations to histological scores (total and cartilage structure) and collagen content. MMP-2, IL-6, IL-8, and MCP-1 had moderate negative correlations, and MMP-9 had a moderate positive correlation, to viability. GRO-α, IL-6, IL-8, and MCP-1 had moderate positive correlations to collagen content. MMP-9, MMP-13, and total MMP activity had moderate negative correlations to tissue GAG. The data suggest links among proinflammatory and degradative pathways are present in OA osteochondral tissues. Further characterization of these links have the potential to delineate mechanisms of disease and diagnostic and therapeutic targets for knee OA.

Cartilage Protective and Immunomodulatory Features of Osteoarthritis Synovial Fluid-Treated Adipose-Derived Mesenchymal Stem Cells Secreted Factors and Extracellular Vesicles-Embedded miRNAs

Intra-articular administration of adipose-derived mesenchymal stem cells (ASCs), either in vitro expanded or within adipose tissue-based products obtained at point-of-care, has gained popularity as innovative regenerative medicine approach for osteoarthritis (OA) treatment. ASCs can stimulate tissue repair and immunomodulation through paracrine factors, both soluble and extracellular vesicles (EV) embedded, collectively defining the secretome. Interaction with the degenerative/inflamed environment is a crucial factor in understanding the finely tuned molecular message but, to date, the majority of reports have described ASC-secretome features in resting conditions or under chemical stimuli far from the in vivo environment of degenerated OA joints. In this report, the secretory profile of ASCs treated with native synovial fluid from OA patients was evaluated, sifting 200 soluble factors and 754 EV-embedded miRNAs. Fifty-eight factors and 223 EV-miRNAs were identified, and discussed in the frame of cartilage and immune cell homeostasis. Bioinformatics gave a molecular basis for M2 macrophage polarization, T cell proliferation inhibition and T reg expansion enhancement, as well as cartilage protection, further confirmed in an in vitro model of OA chondrocytes. Moreover, a strong influence on immune cell chemotaxis emerged. In conclusion, obtained molecular data support the regenerative and immunomodulatory properties of ASCs when interacting with osteoarthritic joint environment.

Coordinate regulation of the senescent state by selective autophagy

Cellular senescence is a complex stress response implicated in aging. Autophagy can suppress senescence but is counterintuitively necessary for full senescence. Although its anti-senescence role is well described, to what extent autophagy contributes to senescence establishment and the underlying mechanisms is poorly understood. Here, we show that selective autophagy of multiple regulatory components coordinates the homeostatic state of senescence. We combined a proteomic analysis of autophagy components with protein stability profiling, identifying autophagy substrate proteins involved in several senescence-related processes. Selective autophagy of KEAP1 promoted redox homeostasis during senescence. Furthermore, selective autophagy limited translational machinery components to ameliorate senescence-associated proteotoxic stress. Lastly, selective autophagy of TNIP1 enhanced senescence-associated inflammation. These selective autophagy networks appear to operate in vivo senescence during human osteoarthritis. Our data highlight a caretaker role of autophagy in the stress support network of senescence through regulated protein stability and unravel the intertwined relationship between two important age-related processes.