Thyroid hormone induces ossification and terminal maturation in a preserved OA cartilage biomimetic model.

OBJECTIVE: To characterize aspects of triiodothyronine (T3) induced chondrocyte terminal maturation within the molecular osteoarthritis pathophysiology using the previously established T3 human ex vivo osteochondral explant model. DESIGNS: RNA-sequencing was performed on explant cartilage obtained from OA patients (n = 8), that was cultured ex vivo with or without T3 (10 ng/ml), and main findings were validated using RT-qPCR in an independent sample set (n = 22). Enrichment analysis was used for functional clustering and comparisons with available OA patient RNA-sequencing and GWAS datasets were used to establish relevance for OA pathophysiology by linking to OA patient genomic profiles. RESULTS: Besides the upregulation of known hypertrophic genes EPAS1 and ANKH, T3 treatment resulted in differential expression of 247 genes with main pathways linked to extracellular matrix and ossification. CCDC80, CDON, ANKH and ATOH8 were among the genes found to consistently mark early, ongoing and terminal maturational OA processes in patients. Furthermore, among the 37 OA risk genes that were significantly affected in cartilage by T3 were COL12A1, TNC, SPARC and PAPPA. CONCLUSIONS: RNA-sequencing results show that metabolic activation and recuperation of growth plate morphology are induced by T3 in OA chondrocytes, indicating terminal maturation is accelerated. The molecular mechanisms involved in hypertrophy were linked to all stages of OA pathophysiology and will be used to validate disease models for drug testing.

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.

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.

Variants of FOXO3 and RPA3 genes affecting IGF-1 levels alter the risk of development of primary osteoarthritis.

INTRODUCTION: Pathologically high growth hormone (GH) and insulin-like growth factor-1 (IGF-1) levels in patients with acromegaly are associated with arthropathy. Several studies highlight the potential role of the GH/IGF-1 axis in primary osteoarthritis (OA). We aimed to disentangle the role of IGF-1 levels in primary OA pathogenesis. METHODS: Patients from the Genetics osteoARthritis and Progression (GARP) Study with familial, generalized, symptomatic OA (n = 337, mean age: 59.8 ± 7.4 years, 82% female) were compared to Leiden Longevity Study (LLS) controls (n = 456, mean age: 59.8 ± 6.8 years, 51% female). Subjects were clinically and radiographically assessed, serum IGF-1 levels were measured, and 10 quantitative trait loci (QTL) in the FOXO3, IGFBP3/TNS3, RPA3, SPOCK2 genes, previously related to serum IGF-1 levels, were genotyped. Linear or binary logistic generalized estimating equation models were performed. RESULTS: Serum IGF-1 levels were increased in OA patients, with male patients exhibiting the strongest effect (males OR = 1.10 (1.04-1.17), P=0.002 vs females OR = 1.04 (1.01-1.07), P = 0.02). Independent of the increased IGF-1 levels, male carriers of the minor allele of FOXO3 QTL rs4946936 had a lower risk to develop hip OA (OR = 0.41 (0.18-0.90), P = 0.026). Additionally, independent of IGF-1 levels, female carriers of the minor alleles of RPA3 QTL rs11769597 had a higher risk to develop knee OA (OR = 1.90 (1.20-2.99), P = 0.006). CONCLUSION: Patients with primary OA had significantly higher IGF-1 levels compared to controls. Moreover, SNPs in the FOXO3 and RPA3 genes were associated with an altered risk of OA. Therefore, altered IGF-1 levels affect the development of OA, and are potentially the result of the pathophysiological OA process.

Neo-cartilage engineered from primary chondrocytes is epigenetically similar to autologous cartilage, in contrast to using mesenchymal stem cells.

OBJECTIVES: To compare the epigenetic landscape of 3D cell models of human primary articular chondrocytes (hPACs) and human bone-marrow derived mesenchymal stem cells (hBMSCs) and their respective autologous articular cartilage. DESIGN: Using Illumina Infinium HumanMethylation450 BeadChip arrays, the DNA methylation landscape of the different cell sources and autologous cartilage was determined. Pathway enrichment was analyzed using DAVID. RESULTS: Principal Component Analysis (PCA) of methylation data revealed separate clustering of hBMSC samples. Between hBMSCs and autologous cartilage 86,881 cytosine-phosphate-guanine dinucleotides (CpGs) (20.2%), comprising 3,034 differentially methylated regions (DMRs; Δß > 0.1; with the same direction of effect), were significantly differentially methylated. In contrast, between hPACs and autologous cartilage only 5,706 CpGs (1.33%) were differentially methylated. Of interest was the finding of the transcriptionally active, hyper-methylation of a Cartilage Intermediate Layer Protein (CILP) annotated DMR (Δß = 0.16) in PAC-cartilage, corresponding to a profound decrease in CILP expression after in vitro culturing of hPACs as compared to autologous cartilage. CONCLUSIONS: In vitro engineered neo-cartilage tissue from primary chondrocytes, hPACs, exhibits a DNA methylation landscape that is almost identical (99% similarity) to autologous cartilage, in contrast to neo-cartilage engineered from bone marrow-derived mesenchymal stem cells (MSCs). Although hBMSCs are widely used for cartilage engineering purposes the effects of these vast differences on cartilage regeneration and long term consequences of implantation, are not known. The use of hBMSCs or hPACs for future cartilage tissue regeneration purposes should therefore be investigated in more depth in future endeavors to better understand the consequences of the differential methylome on neo-cartilage.

Knee and hip articular cartilage have distinct epigenomic landscapes: implications for future cartilage regeneration approaches.

OBJECTIVES: To elucidate the functional epigenomic landscape of articular cartilage in osteoarthritis (OA) affected knee and hip joints in relation to gene expression. METHODS: Using Illumina Infinium HumanMethylation450 BeadChip arrays, genome-wide DNA methylation was measured in 31 preserved and lesioned cartilage sample pairs (14 knees and 17 hips) from patients who underwent a total joint replacement due to primary OA. Using previously published genome-wide expression data of 33 pairs of cartilage samples, of which 13 pairs were overlapping with the current methylation dataset, we assessed gene expression differences in differentially methylated regions (DMRs). RESULTS: Principal component analysis of the methylation data revealed distinct clustering of knee and hip samples, irrespective of OA pathophysiology. A total of 6272 CpG dinucleotides were differentially methylated between the two joints, comprising a total of 357 DMRs containing 1817 CpGs and 245 unique genes. Enrichment analysis of genes proximal of the DMRs revealed significant enrichment for developmental pathways and homeobox (HOX) genes. Subsequent transcriptomic analysis of DMR genes exposed distinct knee and hip expression patterns. CONCLUSIONS: Our findings reveal consistent DMRs between knee and hip articular cartilage that marked transcriptomic differences among HOX genes, which were not reflecting the temporal sequential HOX expression pattern during development. This implies distinct mechanisms for maintaining cartilage integrity in adulthood, thereby contributing to our understanding of cartilage homeostasis and future tissue regeneration approaches.