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.

A robust and standardized method to isolate and expand mesenchymal stromal cells from human umbilical cord.

BACKGROUND AIMS: Human umbilical cord-derived mesenchymal stromal cells (hUC-MSCs) are increasingly used in research and therapy. To obtain hUC-MSCs, a diversity of isolation and expansion methods are applied. Here, we report on a robust and standardized method for hUC-MSC isolation and expansion. METHODS: Using 90 hUC donors, we compared and optimized critical variables during each phase of the multi-step procedure involving UC collection, processing, MSC isolation, expansion and characterization. Furthermore, we assessed the effect of donor-to-donor variability regarding UC morphology and donor attributes on hUC-MSC characteristics. RESULTS: We demonstrated robustness of our method across 90 UC donors at each step of the procedure. With our method, UCs can be collected up to 6 h after birth, and UC-processing can be initiated up to 48 h after collection without impacting on hUC-MSC characteristics. The removal of blood vessels before explant cultures improved hUC-MSC purity. Expansion in Minimum essential medium α supplemented with human platelet lysate increased reproducibility of the expansion rate and MSC characteristics as compared with Dulbecco's Modified Eagle's Medium supplemented with fetal bovine serum. The isolated hUC-MSCs showed a purity of ∼98.9%, a viability of >97% and a high proliferative capacity. Trilineage differentiation capacity of hUC-MSCs was reduced as compared with bone marrow-derived MSCs. Functional assays indicated that the hUC-MSCs were able to inhibit T-cell proliferation demonstrating their immune-modulatory capacity. CONCLUSIONS: We present a robust and standardized method to isolate and expand hUC-MSCs, minimizing technical variability and thereby lay a foundation to advance reliability and comparability of results obtained from different donors and different studies.

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.

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

RNA sequencing data integration reveals an miRNA interactome of osteoarthritis cartilage.

OBJECTIVE: To uncover the microRNA (miRNA) interactome of the osteoarthritis (OA) pathophysiological process in the cartilage. METHODS: We performed RNA sequencing in 130 samples (n=35 and n=30 pairs for messenger RNA (mRNA) and miRNA, respectively) on macroscopically preserved and lesioned OA cartilage from the same patient and performed differential expression (DE) analysis of miRNA and mRNAs. To build an OA-specific miRNA interactome, a prioritisation scheme was applied based on inverse Pearson's correlations and inverse DE of miRNAs and mRNAs. Subsequently, these were filtered by those present in predicted (TargetScan/microT-CDS) and/or experimentally validated (miRTarBase/TarBase) public databases. Pathway enrichment analysis was applied to elucidate OA-related pathways likely mediated by miRNA regulatory mechanisms. RESULTS: We found 142 miRNAs and 2387 mRNAs to be differentially expressed between lesioned and preserved OA articular cartilage. After applying prioritisation towards likely miRNA-mRNA targets, a regulatory network of 62 miRNAs targeting 238 mRNAs was created. Subsequent pathway enrichment analysis of these mRNAs (or genes) elucidated that genes within the 'nervous system development' are likely mediated by miRNA regulatory mechanisms (familywise error=8.4×10-5). Herein NTF3 encodes neurotrophin-3, which controls survival and differentiation of neurons and which is closely related to the nerve growth factor. CONCLUSIONS: By an integrated approach of miRNA and mRNA sequencing data of OA cartilage, an OA miRNA interactome and related pathways were elucidated. Our functional data demonstrated interacting levels at which miRNA affects expression of genes in the cartilage and exemplified the complexity of functionally validating a network of genes that may be targeted by multiple miRNAs.

Isolation and tracing of matrix-producing notochordal and chondrocyte cells using ACAN-2A-mScarlet reporter human iPSC lines.

The development of human induced pluripotent stem cell (iPSC)-based regenerative therapies is challenged by the lack of specific cell markers to isolate differentiated cell types and improve differentiation protocols. This issue is particularly critical for notochordal-like cells and chondrocytes, which are crucial in treating back pain and osteoarthritis, respectively. Both cell types produce abundant proteoglycan aggrecan (ACAN), crucial for the extracellular matrix. We generated two human iPSC lines containing an ACAN-2A-mScarlet reporter. The reporter cell lines were validated using CRISPR-mediated transactivation and functionally validated during notochord and cartilage differentiation. The ability to isolate differentiated cell populations producing ACAN enables their enrichment even in the absence of specific cell markers and allows for comprehensive studies and protocol refinement. ACAN's prevalence in various tissues (e.g., cardiac and cerebral) underscores the reporter's versatility as a valuable tool for tracking matrix protein production in diverse cell types, benefiting developmental biology, matrix pathophysiology, and regenerative medicine.

A human in vitro 3D neo-cartilage model to explore the response of OA risk genes to hyper-physiological mechanical stress.

Objective: Due to the complexity and heterogeneity of osteoarthritis (OA) pathophysiology, studying the interaction between intrinsic molecular changes in chondrocytes after hyper-physiological mechanical stress (MS) and aberrant signalling of OA risk genes remains a challenge. In this study we set out to set up an in vitro 3D neo cartilage pellet model that enables us to explore the responses of OA risk genes to hyper-physiological MS. Design: Human primary chondrocyte neo-cartilage pellets were exposed for 2 days to 2 â€‹× â€‹10 â€‹min of hyper-physiological dynamic MS attained by a 20% strain and a frequency of 5 â€‹Hz. In order to assess cartilage damage, sulphated glycosaminoglycan (sGAG) content in the neo-cartilage was quantified using Alcian blue staining and a dimethyl methylene blue (DMMB) assay, while cleavage of aggrecan was visualized by immunohistochemical staining of aggrecan neo-epitope NITEGE. In addition, changes in expression levels of catabolic, anabolic and hypertrophic genes, and of three OA risk genes; IL11, MGP and TGFA were determined. Results: Hyper-physiological MS induced cartilage damage, as reflected by decreased sGAG content. mRNA levels of aggrecanase ADAMTS5 were increased, while hypertrophic gene RUNX2 was downregulated. MS increased expression of pro-apoptotic marker NOXA. Furthermore, 20% MS led to increased expression of all three OA risk genes IL11, MGP and TGFA. Conclusions: We established a human in vitro model in which hyper-physiological MS induced cartilage damage and catabolic signalling. Next, we demonstrated its usage to study OA risk genes and their response to the mechanical aspects of OA pathophysiology.

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.

Human Osteochondral Explants: Reliable Biomimetic Models to Investigate Disease Mechanisms and Develop Personalized Treatments for Osteoarthritis.

INTRODUCTION: Likely due to ignored heterogeneity in disease pathophysiology, osteoarthritis (OA) has become the most common disabling joint disease, without effective disease-modifying treatment causing a large social and economic burden. In this study we set out to explore responses of aged human osteochondral explants upon different OA-related perturbing triggers (inflammation, hypertrophy and mechanical stress) for future tailored biomimetic human models. METHODS: Human osteochondral explants were treated with IL-1ß (10 ng/ml) or triiodothyronine (T3; 10 nM) or received 65% strains of mechanical stress (65% MS). Changes in chondrocyte signalling were determined by expression levels of nine genes involved in catabolism, anabolism and hypertrophy. Breakdown of cartilage was measured by sulphated glycosaminoglycans (sGAGs) release, scoring histological changes (Mankin score) and mechanical properties of cartilage. RESULTS: All three perturbations (IL-1ß, T3 and 65% MS) resulted in upregulation of the catabolic genes MMP13 and EPAS1. IL-1ß abolished COL2A1 and ACAN gene expression and increased cartilage degeneration, reflected by increased Mankin scores and sGAGs released. Treatment with T3 resulted in a high and significant upregulation of the hypertrophic markers COL1A1, COL10A1 and ALPL. However, 65% MS increased sGAG release and detrimentally altered mechanical properties of cartilage. CONCLUSION: We present consistent and specific output on three different triggers of OA. Perturbation with the pro-inflammatory IL-1ß mainly induced catabolic chondrocyte signalling and cartilage breakdown, while T3 initiated expression of hypertrophic and mineralization markers. Mechanical stress at a strain of 65% induced catabolic chondrocyte signalling and changed cartilage matrix integrity. The major strength of our ex vivo models was that they considered aged, preserved, human cartilage of a heterogeneous OA patient population. As a result, the explants may reflect a reliable biomimetic model prone to OA onset allowing for development of different treatment modalities.

RNA Sequencing Reveals Interacting Key Determinants of Osteoarthritis Acting in Subchondral Bone and Articular Cartilage: Identification of IL11 and CHADL as Attractive Treatment Targets.

OBJECTIVE: To identify key determinants of the interactive pathophysiologic processes in subchondral bone and cartilage in osteoarthritis (OA). METHODS: We performed RNA sequencing on macroscopically preserved and lesional OA subchondral bone from patients in the Research Arthritis and Articular Cartilage study who underwent joint replacement surgery due to OA (n = 24 sample pairs: 6 hips and 18 knees). Unsupervised hierarchical clustering and differential expression analyses were conducted. Results were combined with data on previously identified differentially expressed genes in cartilage (partly overlapping samples) as well as data on recently identified OA risk genes. RESULTS: We identified 1,569 genes that were significantly differentially expressed between lesional and preserved subchondral bone, including CNTNAP2 (fold change [FC] 2.4, false discovery rate [FDR] 3.36 × 10-5 ) and STMN2 (FC 9.6, FDR 2.36 × 10-3 ). Among these 1,569 genes, 305 were also differentially expressed, and with the same direction of effect, in cartilage, including the recently recognized OA susceptibility genes IL11 and CHADL. Upon differential expression analysis with stratification for joint site, we identified 509 genes that were exclusively differentially expressed in subchondral bone of the knee, including KLF11 and WNT4. These genes that were differentially expressed exclusively in the knee were enriched for involvement in epigenetic processes, characterized by, e.g., HIST1H3J and HIST1H3H. CONCLUSION: IL11 and CHADL were among the most consistently differentially expressed genes OA pathophysiology-related genes in both bone and cartilage. As these genes were recently also identified as robust OA risk genes, they classify as attractive therapeutic targets acting on 2 OA-relevant tissues.

Elucidating Epigenetic Regulation by Identifying Functional cis-Acting Long Noncoding RNAs and Their Targets in Osteoarthritic Articular Cartilage.

OBJECTIVE: To identify robustly differentially expressed long noncoding RNAs (lncRNAs) with osteoarthritis (OA) pathophysiology in cartilage and to explore potential target messenger RNA (mRNA) by establishing coexpression networks, followed by functional validation. METHODS: RNA sequencing was performed on macroscopically lesioned and preserved OA cartilage from patients who underwent joint replacement surgery due to OA (n = 98). Differential expression analysis was performed on lncRNAs that were annotated in GENCODE and Ensembl databases. To identify potential interactions, correlations were calculated between the identified differentially expressed lncRNAs and the previously reported differentially expressed protein-coding genes in the same samples. Modulation of chondrocyte lncRNA expression was achieved using locked nucleic acid GapmeRs. RESULTS: By applying our in-house pipeline, we identified 5,053 lncRNAs that were robustly expressed, of which 191 were significantly differentially expressed (according to false discovery rate) between lesioned and preserved OA cartilage. Upon integrating mRNA sequencing data, we showed that intergenic and antisense differentially expressed lncRNAs demonstrate high, positive correlations with their respective flanking sense genes. To functionally validate this observation, we selected P3H2-AS1, which was down-regulated in primary chondrocytes, resulting in the down-regulation of P3H2 gene expression levels. As such, we can confirm that P3H2-AS1 regulates its sense gene P3H2. CONCLUSION: By applying an improved detection strategy, robustly differentially expressed lncRNAs in OA cartilage were detected. Integration of these lncRNAs with differential mRNA expression levels in the same samples provided insight into their regulatory networks. Our data indicates that intergenic and antisense lncRNAs play an important role in regulating the pathophysiology of OA.

Suramin inhibits chikungunya virus replication through multiple mechanisms.

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus that causes severe and often persistent arthritis. In recent years, millions of people have been infected with this virus for which registered antivirals are still lacking. Using our recently established in vitro assay, we discovered that the approved anti-parasitic drug suramin inhibits CHIKV RNA synthesis (IC50 of ∼5µM). The compound inhibited replication of various CHIKV isolates in cell culture with an EC50 of ∼80µM (CC50>5mM) and was also active against Sindbis virus and Semliki Forest virus. In vitro studies hinted that suramin interferes with (re)initiation of RNA synthesis, whereas time-of-addition studies suggested it to also interfere with a post-attachment early step in infection, possibly entry. CHIKV (nsP4) mutants resistant against favipiravir or ribavirin, which target the viral RNA polymerase, did not exhibit cross-resistance to suramin, suggesting a different mode of action. The assessment of the activity of a variety of suramin-related compounds in cell culture and the in vitro assay for RNA synthesis provided more insight into the moieties required for antiviral activity. The antiviral effect of suramin-containing liposomes was also analyzed. Its approved status makes it worthwhile to explore the use of suramin to prevent and/or treat CHIKV infections.