Background aimsHuman 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... Show moreBackground aimsHuman 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.MethodsUsing 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.ResultsWe 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.ConclusionsWe 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. Show less
Tuerlings, M.; Janssen, G.M.C.; Boone, I.; Hoolwerff, M. van; Ruiz, A.R.; Houtman, E.; ... ; Meulenbelt, I. 2023
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... Show moreObjective: 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 (| r| > 0.7) were identified. Among these genes, we identified GJA1, GDF10, STC2, WDR1, and WNK4. Sub-sequent 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, upregu-lation 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 detri-mental 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.(c) 2022 The Author(s). Published by Elsevier Ltd on behalf of Osteoarthritis Research Society International. This is an open access article under the CC BY license (http://creativecommons.org/ licenses/by/4.0/). Show less
The aim of this thesis was to combine transcriptomics, genetics and human disease modelling to obtain further insight into molecular processes underlying osteoarthritis. More specifically, we aimed... Show moreThe aim of this thesis was to combine transcriptomics, genetics and human disease modelling to obtain further insight into molecular processes underlying osteoarthritis. More specifically, we aimed to elucidate the role of long noncoding RNAs expression changes as aberrant epigenetic mechanism in regulating gene expression in chondrocytes. We identified previously unknown long noncoding RNAs associated with the osteoarthritic process and showed enrichment for cis¬-regulation of these long noncoding RNAs with target messenger RNAs.To provide insight in the etiology of osteoarthritis, causal pathways can be identified by unravelling the substantial genetic component. To this end, we investigated the biological functionality of the high-impact, pathogenic mutation identified in the gene fibronectin1 in an early-onset osteoarthritis family. 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.Finally, the common function of fibronectin1 was investigated in cartilage and what changes occur at the transcript level of fibronectin1 with osteoarthritis. Down-regulation of full-length fibronectin was unbeneficial for in vitro chondrogenesis, we hypothesize that this was caused by decreased availability of the classical integrin binding site of fibronectin. Show less
Ruiz, A.R.; Hoolwerff, M. van; Sprangers, S.; Suchiman, E.; Schoenmaker, T.; Dibbets-Schneider, P.; ... ; Meulenbelt, I. 2022
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... Show moreObjectives: 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. Show less
Hoolwerff, M. van; Tuerlings, M.; Wijnen, I.J.L.; Suchiman, H.E.D.; Cats, D.; Mei, H.L.; ... ; Meulenbelt, I. 2022
Objective: To identify FN1 transcripts associated with OA pathophysiology and investigate the downstream effects of modulating FN1 expression and relative transcript ratio. Methods: FN1... Show moreObjective: 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. Show less
Hoolwerff, M. van; Ruiz, A.R.; Bouma, M.; Suchiman, H.E.D.; Koning, R.I.; Jost, C.R.; ... ; Meulenbelt, I. 2021
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... Show moreOsteoarthritis 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. Show less
Tuerlings, M.; Hoolwerff, M. van; Houtman, E.; Suchiman, E.H.E.D.; Lakenberg, N.; Mei, H.L.; ... ; Meulenbelt, I. 2021
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... Show moreObjective 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 x 10(-5)) and STMN2 (FC 9.6, FDR 2.36 x 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. Show less
Houtman, E.; Hoolwerff, M. van; Lakenberg, N.; Suchiman, E.H.D.; Zwaag, E.V.V. van der van der; Nelissen, R.G.H.H.; ... ; Meulenbelt, I. 2021
Introduction: Likely due to ignored heterogeneity in disease pathophysiology, osteoarthritis (OA) has become the most common disabling joint disease, without effective disease-modifying treatment... Show moreIntroduction: 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 beta (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 beta, T3 and 65% MS) resulted in upregulation of the catabolic genes MMP13 and EPAS1. IL-1 beta 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 beta 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. Show less
Hoolwerff, M. van; Metselaar, P.I.; Tuerlings, M.; Suchiman, H.E.D.; Lakenberg, N.; Ramos, Y.F.M.; ... ; Meulenbelt, I. 2020
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... Show moreObjective 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 selectedP3H2-AS1, which was down-regulated in primary chondrocytes, resulting in the down-regulation ofP3H2gene expression levels. As such, we can confirm thatP3H2-AS1regulates its sense geneP3H2. 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. Show less
Almeida, R.C. de; Ramos, Y.F.M.; Mahfouz, A.; Hollander, W. den; Lakenberg, N.; Houtman, E.; ... ; Meulenbelt, I. 2019
