Search

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of the British Society for Rheumatology.)

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of the British Society for Rheumatology.)

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of the British Society for Rheumatology.)

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of the British Society for Rheumatology.)

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-8), AC005165.1 (FC = 0.44, FDR = 2.37 × 10-6) and empty spiracles homeobox 2 opposite strand RNA (EMX2OS, FC = 0.41, FDR = 7.64 × 10-3). 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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the British Society for Rheumatology.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the British Society for Rheumatology.)

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., (© The Author(s) 2020. Published by Oxford University Press on behalf of the British Society for Rheumatology.)

Objectives: Previously, we have shown the involvement of Wnt-activated protein Wnt-1-induced signaling protein 1 (WISP1) in the development of OA in mice. Here, we aimed to characterize the relation between WISP1 expression and human OA and its regulatory epigenetic determinants., Methods: Preserved and lesioned articular cartilage from end-stage OA patients and non-OA-diagnosed individuals was collected. WISP1 expression was determined using immunohistochemistry and damage was classified using Mankin scoring. RNA expression and DNA methylation were assessed in silico from genome-wide datasets (microarray analysis and RNA sequencing, and 450 k-methylationarrays, respectively). Effects of WISP1 were tested in pellet cultures of primary human chondrocytes., Results: WISP1 expression in cartilage of OA patients was increased compared with non-OA-diagnosed controls and, within OA patients, WISP1 was even higher in lesioned compared with preserved regions, with expression strongly correlating with Mankin score. In early symptomatic OA patients with disease progression, higher synovial WISP1 expression was observed as compared with non-progressors. Notably, increased WISP1 expression was inversely correlated with methylation levels of a positional CpG-dinucleotide (cg10191240), with lesioned areas showing strong hypomethylation for this CpG as compared with preserved cartilage. Additionally, we observed that methylation levels were allele-dependent for an intronic single-nucleotide polymorphism nearby cg10191240. Finally, addition of recombinant WISP1 to pellets of primary chondrocytes strongly inhibited deposition of extracellular matrix as reflected by decreased pellet circumference, proteoglycan content and decreased expression of matrix components., Conclusion: Increased WISP1 expression is found in lesioned human articular cartilage, and appears epigenetically regulated via DNA methylation. In vitro assays suggest that increased WISP1 is detrimental for cartilage integrity., (© The Author(s) 2019. Published by Oxford University Press on behalf of the British Society for Rheumatology. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Objectives: To investigate metric properties of four hand mobility tests in hand OA patients, using the OMERACT filter., Methods: Trained assessors examined the Hand Mobility in Scleroderma test (HAMIS), fingertip-to-palm distance (FPD), modified Kapandji index (MKI) and number of hand joints with limited mobility in participants from two cohorts [Genetics ARthrosis and Progression (n = 207) and Hand OSTeoArthritis in Secondary care (n = 174)]. Validity was appraised by assessment of correlations with other outcome measures, and ability to measure thumb vs finger mobility specifically, using cumulative probability plots. The proportion of participants changing in hand mobility based on the smallest detectable difference was calculated for responsiveness. Intraclass correlation coefficients (ICCs) for intra- and interobserver reliability, and feasibility (time to perform tests) were studied in a random sample (n = 20)., Results: Participants displayed large variation in mobility scores. Strongest correlations were observed with structural damage (rs = 0.43-0.52) and bony swelling (rs = 0.46-0.58); correlation patterns were similar among tests. HAMIS, FPD and MKI could all measure finger mobility specifically, but only HAMIS measured thumb mobility particularly. Interobserver reliability was best for HAMIS, ICC 0.90 (95% CI: 0.76, 0.96); intraobserver reliability was excellent for all (ICCs 0.94-0.97). In 2 years, little change was observed; HAMIS was the most sensitive-to-change (smallest detectable difference 3.7% of maximum score). The mean performance time ranged from 0.7 (s.d. 0.5, for FPD) to 5.7 (s.d. 1.3, for HAMIS) min., Conclusion: HAMIS, FPD, MKI and number of joints with limited mobility are all valid, reliable and feasible measures for assessing hand mobility in hand OA, although HAMIS had slightly more favourable properties. Studies assessing sensitivity-to-change in a clinical trial setting are warranted., (© The Author(s) 2017. Published by Oxford University Press on behalf of the British Society for Rheumatology. All rights reserved. For Permissions, please email: journals.permissions@oup.com)