Hyper-physiologic mechanical cues, as an osteoarthritis disease-relevant environmental perturbation, cause a critical shift in set points of methylation at transcriptionally active CpG sites in neo-cartilage organoids.

BACKGROUND: Osteoarthritis (OA) is a complex, age-related multifactorial degenerative disease of diarthrodial joints marked by impaired mobility, joint stiffness, pain, and a significant decrease in quality of life. Among other risk factors, such as genetics and age, hyper-physiological mechanical cues are known to play a critical role in the onset and progression of the disease (Guilak in Best Pract Res Clin Rheumatol 25:815-823, 2011). It has been shown that post-mitotic cells, such as articular chondrocytes, heavily rely on methylation at CpG sites to adapt to environmental cues and maintain phenotypic plasticity. However, these long-lasting adaptations may eventually have a negative impact on cellular performance. We hypothesize that hyper-physiologic mechanical loading leads to the accumulation of altered epigenetic markers in articular chondrocytes, resulting in a loss of the tightly regulated balance of gene expression that leads to a dysregulated state characteristic of the OA disease state. RESULTS: We showed that hyper-physiological loading evokes consistent changes in CpGs associated with expression changes (ML-tCpGs) in ITGA5, CAV1, and CD44, among other genes, which together act in pathways such as anatomical structure morphogenesis (GO:0009653) and response to wound healing (GO:0042060). Moreover, by comparing the ML-tCpGs and their associated pathways to tCpGs in OA pathophysiology (OA-tCpGs), we observed a modest but particular interconnected overlap with notable genes such as CD44 and ITGA5. These genes could indeed represent lasting detrimental changes to the phenotypic state of chondrocytes due to mechanical perturbations that occurred earlier in life. The latter is further suggested by the association between methylation levels of ML-tCpGs mapped to CD44 and OA severity. CONCLUSION: Our findings confirm that hyper-physiological mechanical cues evoke changes to the methylome-wide landscape of chondrocytes, concomitant with detrimental changes in positional gene expression levels (ML-tCpGs). Since CAV1, ITGA5, and CD44 are subject to such changes and are central and overlapping with OA-tCpGs of primary chondrocytes, we propose that accumulation of hyper-physiological mechanical cues can evoke long-lasting, detrimental changes in set points of gene expression that influence the phenotypic healthy state of chondrocytes. Future studies are necessary to confirm this hypothesis.

Identified senescence endotypes in aged cartilage are reflected in the blood metabolome.

Heterogeneous accumulation of senescent cells expressing the senescence-associated secretory phenotype (SASP) affects tissue homeostasis which leads to diseases, such as osteoarthritis (OA). In this study, we set out to characterize heterogeneity of cellular senescence within aged articular cartilage and explored the presence of corresponding metabolic profiles in blood that could function as representative biomarkers. Hereto, we set out to perform cluster analyses, using a gene-set of 131 senescence genes (N = 57) in a previously established RNA sequencing dataset of aged articular cartilage and a generated metabolic dataset in overlapping blood samples. Using unsupervised hierarchical clustering and pathway analysis, we identified two robust cellular senescent endotypes. Endotype-1 was enriched for cell proliferating pathways, expressing forkhead box protein O4 (FOXO4), RB transcriptional corepressor like 2 (RBL2), and cyclin-dependent kinase inhibitor 1B (CDKN1B); the FOXO mediated cell cycle was identified as possible target for endotype-1 patients. Endotype-2 showed enriched inflammation-associated pathways, expressed by interleukin 6 (IL6), matrix metallopeptidase (MMP)1/3, and vascular endothelial growth factor (VEGF)C and SASP pathways were identified as possible targets for endotype-2 patients. Notably, plasma-based metabolic profiles in overlapping blood samples (N = 21) showed two corresponding metabolic clusters in blood. These non-invasive metabolic profiles could function as biomarkers for patient-tailored targeting of senescence in OA.

Mechanical stress and inflammation have opposite effects on Wnt signaling in human chondrocytes.

Dysregulation of Wingless and Int-1 (Wnt) signaling has been strongly associated with development and progression of osteoarthritis (OA). Here, we set out to investigate the independent effects of either mechanical stress (MS) or inflammation on Wnt signaling in human neocartilage pellets, and to relate this Wnt signaling to OA pathophysiology. OA synovium-conditioned media (OAS-CM) was collected after incubating synovium from human end-stage OA joints for 24 h in medium. Cytokine levels in the OAS-CM were determined with a multiplex immunoassay (Luminex). Human neocartilage pellets were exposed to 20% MS, 2% OAS-CM or 1 ng/mL Interleukin-1ß (IL-1ß). Effects on expression levels of Wnt signaling members were determined by reverse transcription-quantitative polymerase chain reaction. Additionally, the expression of these members in articular cartilage from human OA joints was analyzed in association with joint space narrowing (JSN) and osteophyte scores. Protein levels of IL-1ß, IL-6, IL-8, IL-10, tumor necrosis factor α, and granulocyte-macrophage colony-stimulating factor positively correlated with each other. MS increased noncanonical WNT5A and FOS expression. In contrast, these genes were downregulated upon stimulation with OAS-CM or IL-1ß. Furthermore, Wnt inhibitors DKK1 and FRZB decreased in response to OAS-CM or IL-1ß exposure. Finally, expression of WNT5A in OA articular cartilage was associated with increased JSN scores, but not osteophyte scores. Our results demonstrate that MS and inflammatory stimuli have opposite effects on canonical and noncanonical Wnt signaling in human neocartilage. Considering the extent to which MS and inflammation contribute to OA in individual patients, we hypothesize that targeting specific Wnt pathways offers a more effective, individualized approach.

Hyper-physiologic mechanical cues, as an osteoarthritis disease relevant environmental perturbation, cause a critical shift in set-points of methylation at transcriptionally active CpG sites in neo-cartilage organoids.

Background: Osteoarthritis (OA) is a complex, age-related multifactorial degenerative disease of diarthrodial joints marked by impaired mobility, joint stiffness, pain, and a significant decrease in quality of life. Among other risk factors, such as genetics and age, hyper-physiological mechanical cues are known to play a critical role in the onset and progression of the disease (1). It has been shown that post-mitotic cells, such as articular chondrocytes, heavily rely on methylation at CpG sites to adapt to environmental cues and maintain phenotypic plasticity. However, these long-lasting adaptations may eventually have a negative impact on cellular performance. We hypothesize that hyper-physiologic mechanical loading leads to the accumulation of altered epigenetic markers in articular chondrocytes, resulting in a loss of the tightly regulated balance of gene expression that leads to a dysregulated state characteristic of the OA disease state. Results: We showed that hyper-physiological loading evokes consistent changes in ML-tCpGs associated with expression changes in ITGA5, CAV1, and CD44, among other genes, which together act in pathways such as anatomical structure morphogenesis (GO:0009653) and response to wound healing (GO:0042060). Moreover, by comparing the ML-tCpGs and their associated pathways to tCpGs in OA pathophysiology, we observed a modest but particular interconnected overlap with notable genes such as CD44 and ITGA5. These genes could indeed represent lasting detrimental changes to the phenotypic state of chondrocytes due to mechanical perturbations that occurred earlier in life. The latter is further suggested by the association between methylation levels of ML-tCpGs mapped to CD44 and OA severity. Conclusion: Our findings confirm that hyper-physiological mechanical cues evoke changes to the methylome-wide landscape of chondrocytes, concomitant with detrimental changes in positional gene expression levels (ML-tCpGs). Since CAV1, ITGA5, and CD44 are subject to such changes and are central and overlapping with OA-tCPGs of primary chondrocytes, we propose that accumulation of hyper-physiological mechanical cues can evoke long-lasting, detrimental changes in set points of gene expression that influence the phenotypic healthy state of chondrocytes. Future studies are necessary to confirm this hypothesis.

Inhibition of KDM7A/B histone demethylases restores H3K79 methylation and protects against osteoarthritis.

OBJECTIVES: In 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. METHODS: We 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. RESULTS: We 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. CONCLUSION: Targeting KDM7A/B histone demethylases could be an attractive strategy to protect joints against osteoarthritis.

CCN4/WISP1 Promotes Migration of Human Primary Osteoarthritic Chondrocytes.

OBJECTIVES: Previously, we have shown the involvement of cellular communication network factor 4/Wnt-activated protein Wnt-1-induced signaling protein 1 (CCN4/WISP1) in osteoarthritic (OA) cartilage and its detrimental effects on cartilage. Here, we investigated characteristics of CCN4 in chondrocyte biology by exploring correlations of CCN4 with genes expressed in human OA cartilage with functional follow-up. DESIGN: Spearman correlation analysis was performed for genes correlating with CCN4 using our previously established RNA sequencing dataset of human preserved OA cartilage of the RAAK study, followed by a pathway enrichment analysis for genes with ρ ≥|0.6.| Chondrocyte migration in the absence or presence of CCN4 was determined in a scratch assay, measuring scratch size using a live cell imager for up to 36 h. Changes in expression levels of 12 genes, correlating with CCN4 and involved in migratory processes, were determined with reverse transcription-quantitative polymerase chain reaction (RT-qPCR). RESULTS: Correlation of CCN4 with ρ ≥|0.6| was found for 58 genes in preserved human OA cartilage. Pathway analysis revealed "neural crest cell migration" as most significant enriched pathway, containing among others CORO1C, SEMA3C, and SMO. Addition of CCN4 to primary chondrocytes significantly enhance chondrocyte migration as demonstrated by reduced scratch size over the course of 36 h, but at the timepoints measured no effect was observed on mRNA expression of the 12 genes. CONCLUSION: CCN4 increases cell migration of human primary OA chondrocytes. Since WISP1 expression is known to be increased in OA cartilage, this may serve to direct chondrocytes toward cartilage defects and orchestrate repair.

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.

Allelic expression imbalance in articular cartilage and subchondral bone refined genome-wide association signals in osteoarthritis.

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.

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.

Characterizing the secretome of licensed hiPSC-derived MSCs.

Although mesenchymal stromal cells (MSCs) from primary tissues have been successfully applied in the clinic, their expansion capabilities are limited and results are variable. MSCs derived from human-induced pluripotent stem cells (hiMSCs) are expected to overcome these limitations and serve as a reproducible and sustainable cell source. We have explored characteristics and therapeutic potential of hiMSCs in comparison to hBMSCs. RNA sequencing confirmed high resemblance, with average Pearson correlation of 0.88 and Jaccard similarity index of 0.99, and similar to hBMSCs the hiMSCs released extracellular vesicles with in vitro immunomodulatory properties. Potency assay with TNFα and IFNγ demonstrated an increase in well-known immunomodulatory genes such as IDO1, CXCL8/IL8, and HLA-DRA which was also highlighted by enhanced secretion in the media. Notably, expression of 125 genes increased more than 1000-fold. These genes were predicted to be regulated by NFΚB signaling, known to play a central role in immune response. Altogether, our data qualify hiMSCs as a promising source for cell therapy and/or cell-based therapeutic products. Additionally, the herewith generated database will add to our understanding of the mode of action of regenerative cell-based therapies and could be used to identify relevant potency markers.

An epigenome-wide view of osteoarthritis in primary tissues.

Osteoarthritis is a complex degenerative joint disease. Here, we investigate matched genotype and methylation profiles of primary chondrocytes from macroscopically intact (low-grade) and degraded (high-grade) osteoarthritis cartilage and from synoviocytes collected from 98 osteoarthritis-affected individuals undergoing knee replacement surgery. We perform an epigenome-wide association study of knee cartilage degeneration and report robustly replicating methylation markers, which reveal an etiologic mechanism linked to the migration of epithelial cells. Using machine learning, we derive methylation models of cartilage degeneration, which we validate with 82% accuracy in independent data. We report a genome-wide methylation quantitative trait locus (mQTL) map of articular cartilage and synovium and identify 18 disease-grade-specific mQTLs in osteoarthritis cartilage. We resolve osteoarthritis GWAS loci through causal inference and colocalization analyses and decipher the epigenetic mechanisms that mediate the effect of genotype on disease risk. Together, our findings provide enhanced insights into epigenetic mechanisms underlying osteoarthritis in primary tissues.

Inhibiting thyroid activation in aged human explants prevents mechanical induced detrimental signalling by mitigating metabolic processes.

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.

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.

A molecular map of long non-coding RNA expression, isoform switching and alternative splicing in osteoarthritis.

Osteoarthritis is a prevalent joint disease and a major cause of disability worldwide with no curative therapy. Development of disease-modifying therapies requires a better understanding of the molecular mechanisms underpinning disease. A hallmark of osteoarthritis is cartilage degradation. To define molecular events characterizing osteoarthritis at the whole transcriptome level, we performed deep RNA sequencing in paired samples of low- and high-osteoarthritis grade knee cartilage derived from 124 patients undergoing total joint replacement. We detected differential expression between low- and high-osteoarthritis grade articular cartilage for 365 genes and identified a 38-gene signature in osteoarthritis cartilage by replicating our findings in an independent dataset. We also found differential expression for 25 novel long non-coding RNA genes (lncRNAs) and identified potential lncRNA interactions with RNA-binding proteins in osteoarthritis. We assessed alterations in the relative usage of individual gene transcripts and identified differential transcript usage for 82 genes, including ABI3BP, coding for an extracellular matrix protein, AKT1S1, a negative regulator of the mTOR pathway and TPRM4, coding for a transient receptor potential channel. We further assessed genome-wide differential splicing, for the first time in osteoarthritis, and detected differential splicing for 209 genes, which were enriched for extracellular matrix, proteoglycans and integrin surface interactions terms. In the largest study of its kind in osteoarthritis, we find that isoform and splicing changes, in addition to extensive differences in both coding and non-coding sequence expression, are associated with disease and demonstrate a novel layer of genomic complexity to osteoarthritis pathogenesis.

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.

The role of TNFRSF11B in development of osteoarthritic cartilage.

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.

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.

Higher thyrotropin leads to unfavorable lipid profile and somewhat higher cardiovascular disease risk: evidence from multi-cohort Mendelian randomization and metabolomic profiling.

BACKGROUND: Observational studies suggest interconnections between thyroid status, metabolism, and risk of coronary artery disease (CAD), but causality remains to be proven. The present study aimed to investigate the potential causal relationship between thyroid status and cardiovascular disease and to characterize the metabolomic profile associated with thyroid status. METHODS: Multi-cohort two-sample Mendelian randomization (MR) was performed utilizing genome-wide significant variants as instruments for standardized thyrotropin (TSH) and free thyroxine (fT4) within the reference range. Associations between TSH and fT4 and metabolic profile were investigated in a two-stage manner: associations between TSH and fT4 and the full panel of 161 metabolomic markers were first assessed hypothesis-free, then directional consistency was assessed through Mendelian randomization, another metabolic profile platform, and in individuals with biochemically defined thyroid dysfunction. RESULTS: Circulating TSH was associated with 52/161 metabolomic markers, and fT4 levels were associated with 21/161 metabolomic markers among 9432 euthyroid individuals (median age varied from 23.0 to 75.4 years, 54.5% women). Positive associations between circulating TSH levels and concentrations of very low-density lipoprotein subclasses and components, triglycerides, and triglyceride content of lipoproteins were directionally consistent across the multivariable regression, MR, metabolomic platforms, and for individuals with hypo- and hyperthyroidism. Associations with fT4 levels inversely reflected those observed with TSH. Among 91,810 CAD cases and 656,091 controls of European ancestry, per 1-SD increase of genetically determined TSH concentration risk of CAD increased slightly, but not significantly, with an OR of 1.03 (95% CI 0.99-1.07; p value 0.16), whereas higher genetically determined fT4 levels were not associated with CAD risk (OR 1.00 per SD increase of fT4; 95% CI 0.96-1.04; p value 0.59). CONCLUSIONS: Lower thyroid status leads to an unfavorable lipid profile and a somewhat increased cardiovascular disease risk.

Circulating MicroRNAs Highly Correlate to Expression of Cartilage Genes Potentially Reflecting OA Susceptibility-Towards Identification of Applicable Early OA Biomarkers.

OBJECTIVE: To identify and validate circulating micro RNAs (miRNAs) that mark gene expression changes in articular cartilage early in osteoarthritis (OA) pathophysiology process. METHODS: Within the ongoing RAAK study, human preserved OA cartilage and plasma (N = 22 paired samples) was collected for RNA sequencing (respectively mRNA and miRNA). Spearman correlation was determined for 114 cartilage genes consistently and significantly differentially expressed early in osteoarthritis and 384 plasma miRNAs. Subsequently, the minimal number of circulating miRNAs serving to discriminate between progressors and non-progressors was assessed by regression analysis and area under receiver operating curves (AUC) was calculated with progression data and plasma miRNA sequencing from the GARP study (N = 71). RESULTS: We identified strong correlations (ρ ≥ |0.7|) among expression levels of 34 unique plasma miRNAs and 21 genes, including 4 genes that correlated with multiple miRNAs. The strongest correlation was between let-7d-5p and EGFLAM (ρ = -0.75, P = 6.9 × 10-5). Regression analysis of the 34 miRNAs resulted in a set of 7 miRNAs that, when applied to the GARP study, demonstrated clinically relevant predictive value with AUC > 0.8 for OA progression over 2 years and near-clinical value for progression over 5 years- (AUC = 0.8). CONCLUSIONS: We show that plasma miRNAs levels reflect gene expression levels in cartilage and can be exploited to represent ongoing pathophysiological processes in articular cartilage. We advocate that identified signature of 7 plasma miRNAs can contribute to direct further studies toward early biomarkers predictive for progression of osteoarthritis over 2 and 5 years.

Cripto favors chondrocyte hypertrophy via TGF-ß SMAD1/5 signaling during development of osteoarthritis.

Chondrocytes in mice developing osteoarthritis (OA) exhibit an aberrant response to the secreted cytokine transforming growth factor (TGF)-ß, consisting in a potentiation of intracellular signaling downstream of the transmembrane type I receptor kinase activin receptor-like kinase (ALK)1 against canonical TGF-ß receptor ALK5-mediated signaling. Unfortunately, the underlying mechanisms remain elusive. In order to identify novel druggable targets for OA, we aimed to investigate novel molecules regulating the ALK1/ALK5 balance in OA chondrocytes. We performed gene expression analysis of TGF-ß signaling modulators in joints from three different mouse models of OA and found an upregulated expression of the TGF-ß co-receptor Cripto (Tdgf1), which was validated in murine and human cartilage OA samples at the protein level. In vitro and ex vivo, elevated expression of Cripto favors the hypertrophic differentiation of chondrocytes, eventually contributing to tissue calcification. Furthermore, we found that Cripto participates in a TGF-ß-ALK1-Cripto receptor complex in the plasma membrane, thereby inducing catabolic SMAD1/5 signaling in chondrocytes. In conclusion, we demonstrate that Cripto is expressed in OA and plays a functional role promoting chondrocyte hypertrophy, thereby becoming a novel potential therapeutic target in OA, for which there is no efficient cure or validated biomarker. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.