Search

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.

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.

ObjectivesThe etiology of osteoarthritis revolves around the interplay between genetic predisposition and perturbing environmental cues, such as mechanical stress. The pericellular matrix, with its hallmark proteins collagen type VI and fibronectin, surrounds chondrocytes and is critical in transducing the biomechanical cues. The objective is to study the functional effects of an OA disease-risk mutation inCOL6A3in interaction with hyper-physiological mechanical cues in a tailored human induced pluripotent stem cells (hiPSCs) derived cartilage organoid model.MethodTo identify pathogenic OA mutations exome sequencing in symptomatic OA patients was performed. To study functional effects, CRISPR-Cas9 genome engineering was used to introduce the mutation in our established human induced pluripotent stem cell-derived in-vitro neo-cartilage organoid model in interaction with hyper-physiological mechanical loading conditions.ResultsA high-impact mutation inCOL6A3was identified that resulted in significantly lower binding between the PCM proteins COLVI and fibronectin (FN) and provoked an osteoarthritic chondrocyte state. Moreover, aberrant function of the PCM, secondary to theCOL6A3mutation, abolished the initial stress responses marked particularly by upregulation ofPTGS2encoding cyclooxygenase-2 (COX-2), after hyper-physiological mechanical loading conditions.ConclusionThese findings demonstrate that ablating the characteristic transient COX-2 response after injurious mechanical cues may have a direct negative impact on chondrocyte health.What is already knownThe etiology of osteoarthritis revolves around the interplay between genetic predisposition and perturbing environmental cues, such as mechanical stress.The pericellular matrix, with its hallmark proteins collagen type VI and fibronectin, surrounds the chondrocytes and is critical in transducing biomechanical cues from the extracellular matrix to chondrocytes henceforth it determines the chondrocyte mechanical environment.The mechanical environment of the chondrocytes is a critical factor that influences chondrocyte health as it determines the balance between synthesis and degradation of the articular cartilage extracellular matrix.What this study addsA sustainable human induced pluripotent stem cell-derived in-vitro neo-cartilage organoid model that is tailored to study detailed biologic effects of mechanical cues to chondrocytes.An OA disease-risk mutation inCOL6A3reduces the binding between collagen type VI to fibronectin and provoked an osteoarthritic chondrocyte state.Upon hyper-physiological mechanical loading, aberrant function of the pericellular matrix, secondary to theCOL6A3mutation, ablates the initial transient inflammatory response, characterized particularly byPTGS2encoding cyclooxygenase-2 (COX-2).How this study might affect research practice or policyInhibiting COX-2, as an important transient inflammatory response after hyper-physiological mechanical cues, could worsen the loss of structural integrity of the cartilage in osteoarthritis patients. Henceforth, prescription of COX-2 inhibitors as pain treatment for OA patients should be reconsidered.

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 genome-wide association studies (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 two 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.7×10−4) and ILF3-DT (meta-φ = 0.6, FDR = 1.75×10−5). Moreover, 12 drugs were interacting with seven genes displaying AEI, of which seven drugs have been already approved. Conclusions By prioritizing proxy transcript SNPs that mark AEI in cartilage and/or subchondral bone at loci harbouring 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.

ObjectivesIn 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.MethodsWe 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.ResultsWe 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.ConclusionTargeting KDM7A/B histone demethylases could be an attractive strategy to protect joints against osteoarthritis.

ObjectiveTo 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.MethodsCo-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.ResultsUpon 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.ConclusionsMimicking 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.

Objectives OA is a complex genetic disease with different risk factors contributing to its development. One of the genes, TNFRSF11B, previously identified with gain-of-function mutation in a family with early-onset OA with chondrocalcinosis, is among the highest upregulated genes in lesioned OA cartilage (RAAK-study). Here, we determined the role of TNFRSF11B overexpression in development of OA. Methods Human primary articular chondrocytes (9 donors RAAK study) were transduced using lentiviral particles with or without TNFRSF11B. Cells were cultured for 1 week in a 3 D in-vitro chondrogenic model. TNFRSF11B overexpression was confirmed by RT-qPCR, immunohistochemistry and ELISA. Effects of TNFRSF11B overexpression on cartilage matrix deposition, matrix mineralization, and genes highly correlated to TNFRSF11B in RNA-sequencing dataset (r >0.75) were determined by RT-qPCR. Additionally, glycosaminoglycans and collagen deposition were visualized with Alcian blue staining and immunohistochemistry (COL1 and COL2). Results Overexpression of TNFRSF11B resulted in strong upregulation of MMP13, COL2A1 and COL1A1. Likewise, mineralization and osteoblast characteristic markers RUNX2, ASPN and OGN showed a consistent increase. Among 30 genes highly correlated to TNFRSF11B, expression of only eight changed significantly, with BMP6 showing the highest increase (9-fold) while expression of RANK and RANKL remained unchanged indicating previously unknown downstream pathways of TNFRSF11B in cartilage. Conclusion Results of our 3D in vitro chondrogenesis model indicate that upregulation of TNFRSF11B in lesioned OA cartilage may act as a direct driving factor for chondrocyte to osteoblast transition observed in OA pathophysiology. This transition does not appear to act via the OPG/RANK/RANKL triad common in bone remodeling.