Osteoarthritis year in review: genetics, genomics, epigenetics.

OBJECTIVE: In this review, we have highlighted the advances over the past year in genetics, genomics and epigenetics in the field of osteoarthritis (OA). METHODS: A literature search of PubMed was performed using the criteria: "osteoarthritis" and one of the following terms "genetic(s), genomic(s), epigenetic(s), polymorphism, noncoding ribonucleic acid (RNA), microRNA, long noncoding RNA, lncRNA, circular RNA, RNA sequencing (RNA-seq), single cell sequencing, transcriptomics, or deoxyribonucleic acid (DNA) methylation between April 01, 2020 and April 30, 2021. RESULTS: In total we identified 765 unique publications, which eventually reduced to 380 of relevance to the field as judged by two assessors. Many of these studies included multiple search terms. We summarised advances relating to genetics, functional genetics, genomics and epigenetics, focusing on our personal key papers during the year. CONCLUSIONS: This year few studies have identified new genetic variants contributing to OA susceptibility, but a focus has been on refining risk loci or their functional validation. The use of new technologies together with investigating the cross-talk between multiple tissue types, greater sample sizes and/or better patient classification (OA subtypes) will continue to increase our knowledge of disease mechanisms and progress towards understanding and treating OA.

Osteoarthritis year in review 2019: genetics, genomics and epigenetics.

Although osteoarthritis (OA) aetiology is complex, genetic, genomic and epigenetic studies published within the last decade have advanced our understanding of the molecular processes underlying this common musculoskeletal disease. The purpose of this narrative review is to highlight the key research articles within the OA genetics, genomics and epigenetics fields that were published between April 2018 and April 2019. The review focuses on the identification of new OA genetic risk loci, genomics techniques that have been used for the first time in human cartilage and new publicly available databases, and datasets that will aid OA functional studies. Fifty-six new OA susceptibility loci were identified by two large scale genome wide association study meta-analyses, increasing the number of genome-wide significant risk loci to 90. OA risk variants are enriched near genes involved in skeletal development and morphology, and show genetic overlap with height, hip shape, bone area and developmental dysplasia of the hip. Several functional studies of OA loci were published, including a genome-wide analysis of genetic variation on cartilage gene expression. A specialised data portal for exploring cross-species skeletal transcriptomic datasets has been developed, and the first use of cartilage single cell RNAseq analysis reported. This year also saw the systematic identification of all microRNAs, long non-coding RNAs and circular RNAs expressed in human OA cartilage. Putative transcriptional regulatory regions have been mapped in human chondrocytes genome-wide, providing a dataset that will facilitate the prioritisation and characterisation of OA genetic and epigenetic loci.

Identification of long non-coding RNAs expressed in knee and hip osteoarthritic cartilage.

OBJECTIVE: Long intergenic non-coding RNAs (lincRNAs) are emerging as key regulators in gene expression; however, little is known about the lincRNA expression changes that occur in osteoarthritis (OA). Here we aimed to define a transcriptome of lncRNAs in OA cartilage, specifically comparing the lincRNA transcriptome of knee and hip cartilage. METHOD: RNA-seq was performed on nucleic acid extracted from hip cartilage from patients undergoing joint replacement surgery because of either OA (n = 10) or because of a neck of femur fracture (NOF; n = 6). After transcript alignment, counts were performed using Salmon and differential expression for ENSEMBL lincRNAs determined using DESeq2. Hip RNA-seq lincRNA expression was compared to a knee dataset (ArrayExpress; E-MTAB-4304). ChIP-seq data from ENCODE was used to determine whether lincRNAs were associated with promoters (plncRNA) or unidirectional enhancer-like regulatory elements (elncRNAs). RESULTS: Our analysis of the hip transcriptome identified 1692 expressed Transcripts Per Million (TPM ≥1) Ensembl lincRNAs, of which 198 were significantly (FDR ≤0.05) differentially expressed in OA vs normal (NOF) cartilage. Similar analysis of knee cartilage transcriptome identified 648 Emsembl lincRNAs with 93 significantly (FDR ≤0.05) differentially expressed in intact vs damaged cartilage. In total, 1834 lincRNAs were expressed in both hip and knee cartilage, with a highly significant correlation in expression between the two cartilages. CONCLUSION: This is the first study to use RNA-seq to map and compare the lincRNA transcriptomes of hip and knee cartilage. We propose that lincRNAs expressed selectively in cartilage, or showing differential expression in OA, will play a role in cartilage homoeostasis.

Dickkopf-3 is upregulated in osteoarthritis and has a chondroprotective role.

OBJECTIVE: Dickkopf-3 (Dkk3) is a non-canonical member of the Dkk family of Wnt antagonists and its upregulation has been reported in microarray analysis of cartilage from mouse models of osteoarthritis (OA). In this study we assessed Dkk3 expression in human OA cartilage to ascertain its potential role in chondrocyte signaling and cartilage maintenance. METHODS: Dkk3 expression was analysed in human adult OA cartilage and synovial tissues and during chondrogenesis of ATDC5 and human mesenchymal stem cells. The role of Dkk3 in cartilage maintenance was analysed by incubation of bovine and human cartilage explants with interleukin-1ß (IL1ß) and oncostatin-M (OSM). Dkk3 gene expression was measured in cartilage following murine hip avulsion. Whether Dkk3 influenced Wnt, TGFß and activin cell signaling was assessed in primary human chondrocytes and SW1353 chondrosarcoma cells using qRT-PCR and luminescence assays. RESULTS: Increased gene and protein levels of Dkk3 were detected in human OA cartilage, synovial tissue and synovial fluid. DKK3 gene expression was decreased during chondrogenesis of both ATDC5 cells and humans MSCs. Dkk3 inhibited IL1ß and OSM-mediated proteoglycan loss from human and bovine cartilage explants and collagen loss from bovine cartilage explants. Cartilage DKK3 expression was decreased following hip avulsion injury. TGFß signaling was enhanced by Dkk3 whilst Wnt3a and activin signaling were inhibited. CONCLUSIONS: We provide evidence that Dkk3 is upregulated in OA and may have a protective effect on cartilage integrity by preventing proteoglycan loss and helping to restore OA-relevant signaling pathway activity. Targeting Dkk3 may be a novel approach in the treatment of OA.

Detecting new microRNAs in human osteoarthritic chondrocytes identifies miR-3085 as a human, chondrocyte-selective, microRNA.

OBJECTIVE: To use deep sequencing to identify novel microRNAs (miRNAs) in human osteoarthritic cartilage which have a functional role in chondrocyte phenotype or function. DESIGN: A small RNA library was prepared from human osteoarthritic primary chondrocytes using in-house adaptors and analysed by Illumina sequencing. Novel candidate miRNAs were validated by northern blot and qRT-PCR. Expression was measured in cartilage models. Targets of novel candidates were identified by microarray and computational analysis, validated using 3'-UTR-luciferase reporter plasmids. Protein levels were assessed by western blot and functional analysis by cell adhesion. RESULTS: We identified 990 known miRNAs and 1621 potential novel miRNAs in human osteoarthritic chondrocytes, 60 of the latter were expressed in all samples assayed. MicroRNA-140-3p was the most highly expressed microRNA in osteoarthritic cartilage. Sixteen novel candidate miRNAs were analysed further, of which six remained after northern blot analysis. Three novel miRNAs were regulated across models of chondrogenesis, chondrocyte differentiation or cartilage injury. One sequence (novel #11), annotated in rodents as microRNA-3085-3p, was preferentially expressed in cartilage, dependent on chondrocyte differentiation and, in man, is located in an intron of the cartilage-expressed gene CRTAC-1. This microRNA was shown to target the ITGA5 gene directly (which encodes integrin alpha5) and inhibited adhesion to fibronectin (dependent on alpha5beta1 integrin). CONCLUSION: Deep sequencing has uncovered many potential microRNA candidates expressed in human cartilage. At least three of these show potential functional interest in cartilage homeostasis and osteoarthritis (OA). Particularly, novel #11 (microRNA-3085-3p) which has been identified for the first time in man.

Identification of the pathogenic pathways in osteoarthritic hip cartilage: commonality and discord between hip and knee OA.

OBJECTIVE: To define for the first time the transcriptomes of normal and end-stage osteoarthritis (OA) hip cartilage. MATERIALS AND METHODS: RNA was isolated from cartilage within 2h of joint replacement surgery. Gene expression was analyzed using Agilent GeneSpring GX 11 following hybridization to Illumina Human HT-12 V3 microarrays. Real-time reverse-transcription polymerase chain reaction (RT-PCR) was used to validate the expression of six genes identified by microarray as differentially expressed. Gene Set Enrichment Analysis (GSEA) and Ingenuity Pathway Analysis (IPA) were used to investigate enriched functions or canonical pathways amongst differentially expressed genes respectively. RESULTS: In total we identified 998 differentially expressed genes (fold change ≥ ±1.5, P-value ≤ 0.01) between neck of femur fracture (NOF) (n = 10) and OA hip (n = 9) patient cartilage. These differentially expressed genes were enriched within 71 canonical pathways. A comparison between a comparable knee dataset(20) only identified 229 genes similarly differentially expressed although remarkably 34 canonical pathways overlapped between experiments. CONCLUSIONS: This study is the first to report a comprehensive gene expression analysis of human hip OA cartilage compared to control (NOF) cartilage at the whole-genome level. Our differential gene expression dataset shows excellent correlation with similar defined studies using comparable tissue but reveals discord between hip and knee OA at the individual gene status but with commonality with regards the molecular pathways involved.

Epigenetic mechanisms in cartilage and osteoarthritis: DNA methylation, histone modifications and microRNAs.

Osteoarthritis (OA) is a complex multifactorial disease with a strong genetic component. Several studies have suggested or identified epigenetic events that may play a role in OA progression and the gene expression changes observed in diseased cartilage. The aim of this review is to inform about current research in epigenetics and epigenetics in OA. Epigenetic mechanisms include DNA methylation, histone modifications, and microRNAs. Collectively, these enable the cell to respond quickly to environmental changes and can be inherited during cell division. However, aberrant epigenetic modifications are associated with a number of pathological conditions, including OA. Advancements in epigenetic research suggests that global analysis of such modifications in OA are now possible, however, with the exception of microRNAs, it will be a significant challenge to demonstrate how such events impact on the disease.

Differential Toll-like receptor-dependent collagenase expression in chondrocytes.

OBJECTIVES: To characterise the catabolic response of osteoarthritic chondrocytes to Toll-like receptor (TLR) ligands. METHODS: Induction of the collagenases, matrix metalloproteinase (MMP)1 and MMP13, by TLR ligands was assessed in chondrocytes by real-time reverse transcriptase (RT)-PCR. TLR signalling pathway activation and their involvement in collagenase induction were confirmed by immunoblotting and use of pathway inhibitors and siRNA. TLR expression was compared in the femoral head cartilage of normal controls and patients with osteoarthritis (OA) by real-time RT-PCR. RESULTS: Ligands for TLR6/2 and TLR3 showed the greatest upregulation of MMP1 and MMP13 respectively, although all TLR ligands upregulated these MMPs. MMP1 and MMP13 induction by TLR3 and TLR1/2 or TLR6/2 ligands were dependent on Trif and MyD88, respectively. These inductions were dependent upon the nuclear factor (NF)kappaB pathway, but were differentially inhibited by various mitogen-activated protein kinase inhibitors, with MMP13 induction most reliant on the extracellular signal-regulated kinase pathway. In addition, ligands for TLR1/2 and TLR6/2, but not TLR3, induced significant collagenolysis in a cartilage resorption assay. Finally, TLR2 was significantly downregulated and TLR3 upregulated in OA, compared to normal, cartilage. CONCLUSIONS: Activation of chondrocyte TLRs leads to differential collagenase gene activation. Treatment of chondrocytes with TLR1/2 or TLR6/2 ligands resulted in collagen resorption. The modulated expression of chondrocyte TLR2 and TLR3 in OA cartilage, compared to normal, may reflect a response to repair cartilage or prevent further extracellular matrix destruction. These data suggest modulation of TLR-mediated signalling as a potential therapeutic strategy for the treatment of OA.