Sequencing identifies a distinct signature of circulating microRNAs in early radiographic knee osteoarthritis

MicroRNA signature for early prediction of knee osteoarthritis structural progression using integrated machine and deep learning approaches

OBJECTIVE

Conventional methodologies are ineffective in predicting the rapid progression of knee osteoarthritis (OA). MicroRNAs (miRNAs) show promise as biomarkers for patient stratification. We aimed to develop a miRNA prognosis model for identifying knee OA structural progressors/non-progressors using integrated machine/deep learning tools.

METHODS

Baseline serum miRNAs from Osteoarthritis Initiative (OAI) participants were isolated and sequenced. Participants were categorized based on their likelihood of knee structural progression/non-progression using magnetic resonance imaging and X-ray data. For prediction model development, 152 OAI participants (91 progressors, 61 non-progressors) were used. MiRNA features were reduced through VarClusHi clustering. Key miRNAs and OA determinants (age, sex, body mass index, race) were identified using seven machine learning tools. The final prediction model was developed using advanced machine/deep learning techniques. Model performance was assessed with area under the curve (AUC) (95% confidence intervals) and accuracy. Monte Carlo cross-validation ensured robustness. Model validation used 30 OAI baseline plasma samples from an independent set of participants (14 progressors, 16 non-progressors).

RESULTS

Feature clustering selected 107 miRNAs. Elastic Net was chosen for feature selection. An optimized prediction model based on an Artificial Neural Network comprising age and four miRNAs (hsa-miR-556-3p, hsa-miR-3157-5p, hsa-miR-200a-5p, hsa-miR-141-3p) exhibited excellent performance (AUC, 0.94 [0.89, 0.97]; accuracy, 0.84 [0.77, 0.89]). Model validation performance (AUC, 0.81 [0.63, 0.92]; accuracy, 0.83 [0.66, 0.93]) demonstrated the potential for generalization.

CONCLUSION

This study introduces a novel miRNA prognosis model for knee OA patients at risk of structural progression. It requires five baseline features, demonstrates excellent performance, is validated with an independent set, and holds promise for future personalized therapeutic monitoring.

Circulating miR-126-3p is a mechanistic biomarker for knee osteoarthritis

Osteoarthritis is a major contributor to pain and disability worldwide, yet there are currently no validated soluble biomarkers or disease-modifying treatments. Given that microRNAs are promising mechanistic biomarkers that can be therapeutically targeted, in this study, we aimed to identify and prioritize reproducible circulating microRNAs associated with radiographic knee osteoarthritis. Across four independent cohorts, we find circulating miR-126-3p is elevated in knee osteoarthritis versus controls. Across six primary human knee osteoarthritis tissues, miR-126-3p is highest in subchondral bone, fat pad and synovium, and lowest in cartilage. Following both intravenous and intra-articular miR-126-3p mimic treatment in a surgical mouse model of knee osteoarthritis, we show reduced disease severity in males. In human knee osteoarthritis biospecimens, miR-126-3p mimic treatment reduces genes and markers associated with angiogenesis, as well as genes linked to osteogenesis, adipogenesis, and synovitis-processes secondary to angiogenesis. Our findings indicate that miR-126-3p is elevated in knee osteoarthritis and mitigates disease severity, supporting its potential as a biomarker and therapeutic target.

Progress in multi-omics studies of osteoarthritis

Osteoarthritis (OA), a ubiquitous degenerative joint disorder, is marked by pain and disability, profoundly impacting patients' quality of life. As the population ages, the global prevalence of OA is escalating. Omics technologies have become instrumental in investigating complex diseases like OA, offering comprehensive insights into its pathogenesis and progression by uncovering disease-specific alterations across genomics, transcriptomics, proteomics, and metabolomics levels. In this review, we systematically analyzed and summarized the application and recent achievements of omics technologies in OA research by scouring relevant literature in databases such as PubMed. These studies have shed light on new potential therapeutic targets and biomarkers, charting fresh avenues for OA diagnosis and treatment. Furthermore, in our discussion, we highlighted the immense potential of spatial omics technologies in unraveling the molecular mechanisms of OA and in the development of novel therapeutic strategies, proposing future research directions and challenges. Collectively, this study encapsulates the pivotal advances in current OA research and prospects for future investigation, providing invaluable references for a deeper understanding and treatment of OA. This review aims to synthesize the recent progress of omics technologies in the realm of OA, aspiring to furnish theoretical foundations and research orientations for more profound studies of OA in the future.

Identification and characterisation of temporal abundance of microRNAs in synovial fluid from an experimental equine model of osteoarthritis

BACKGROUND

MicroRNAs, a class of small noncoding RNAs, serve as post-transcriptional regulators of gene expression and are present in a stable and quantifiable form in biological fluids. MicroRNAs may influence intra-articular responses and the course of disease, but very little is known about their temporal changes in osteoarthritis.

OBJECTIVES

To identify miRNAs and characterise the temporal changes in their abundance in SF from horses with experimentally induced osteoarthritis. We hypothesised that the abundance of miRNA would change during disease progression.

STUDY DESIGN

In vivo experiments.

METHODS

RNA extracted from synovial fluid obtained sequentially (Day 0, 28 and 70) from nine horses with experimentally induced osteoarthritis was subjected to small RNA sequencing using the Illumina Hiseq 4000 sequencing platform. Differentially abundant miRNAs underwent further validation and mapping of temporal abundance (Day 0, 14, 17, 21, 28, 35, 42, 49, 56, 63 and 70 days after osteoarthritis induction) by microfluidic reverse transcription quantitative real-time PCR. Bioinformatic analyses were performed to predict potential biological associations and target genes of the differentially abundant microRNAs.

RESULTS

Small RNA sequencing revealed 61 differentially abundant microRNAs at an early osteoarthritis stage (Day 28), and subsequent reverse transcription quantitative real-time PCR analysis validated 20 of these. Significant biological functions of the differentially abundant microRNAs were apoptosis, necrosis, cell proliferation and cell invasion. Following validation, four microRNAs (miRNA-199b-3p, miRNA-139-5p, miRNA-1839 and miRNA-151-5p) were detected in more than 50% of the synovial fluid samples and had higher abundance in osteoarthritic than in control joints.

MAIN LIMITATIONS

Limited sample size.

CONCLUSION

This is the first study to determine longitudinal changes in synovial fluid microRNA abundance in an equine model of osteoarthritis. Larger studies are needed in naturally occurring osteoarthritis to interrogate putative changes identified by this study.

The Study of Small RNA Sequencing from Biological Samples (miRNA-seq)

Small RNA sequencing (sRNA-seq) has greatly transformed the study of molecular biology by allowing thorough analysis of several types of short RNA molecules, such as microRNAs (miRNAs), small interfering RNAs (siRNAs), and piwi-interacting RNAs (piRNAs). This chapter presents a comprehensive methodology for conducting miRNA -seq on biological materials. The text introduces the principles and uses of miRNA-seq, emphasizing its importance in gene regulation, disease processes, and therapeutic interventions. The chapter subsequently explores the essential stages of miRNA-seq, encompassing sample acquisition, RNA isolation, library construction, sequencing, and data interpretation. Every part contains pragmatic suggestions, guidance for resolving issues, and optimal methods to guarantee excellent outcomes. Upon completion of this chapter, readers will possess a comprehensive comprehension of the miRNA-seq workflow and will be adequately prepared to apply this potent approach in their own research endeavors.

Investigating the Impact of Circulating MicroRNAs on Knee and Hip Osteoarthritis: Causal Links, Biological Mechanisms, and Drug Interactions

Osteoarthritis (OA), particularly in the knee and hip, poses a significant global health challenge due to limited therapeutic options. To elucidate the molecular mechanisms of OA and identify potential biomarkers and therapeutic targets, we utilized genome-wide association studies (GWAS) and cis-miRNA expression quantitative trait loci (cis-miR-eQTL) datasets to identify miRNAs associated with OA, revealing 16 that were linked to knee OA and 21 to hip OA. Among these, hsa-miR-1303 was significantly upregulated in both knee and hip OA (IVW: p = 6.8164×10-36 and 4.7919×10-2 respectively, OR > 1) and identified as a key factor in disease progression. Hsa-miR-1303 potentially regulates 30 genes involved in critical signaling pathways, such as the neurotrophin signaling pathway, and interacts with competing endogenous RNAs (ceRNAs) like circ_0041843 and LINC01338, thereby influencing key regulatory proteins such as SUMO2 and PARP1. Pharmacologically, hsa-miR-1303 targets nine druggable genes, including NRAS, H2AZ1, and RPS3, which have implications for drugs like cantharidin and diindolylmethane, potentially critical for developing novel OA treatments. Conversely, hsa-miR-125a-5p and hsa-miR-125b-5p, which are downregulated in both knee and hip OA, are associated with pathways such as HIF-1 and JAK-STAT, which modulate apoptotic signaling and transcriptional regulation. These miRNAs also interact with ceRNAs such as circ_0000254 and SPACA6P-AS, impacting proteins like STAT3, MCL1, and TRAF6. A drug interaction analysis identified 47 potential treatments, including Resveratrol and Acetaminophen, suggesting new therapeutic possibilities for OA management. This study not only highlights the role of miRNAs like hsa-miR-1303 and hsa-miR-125 in OA but also opens avenues for miRNA-based therapeutic development.

A scoping review of statistical methods to investigate colocalization between genetic associations and microRNA expression in osteoarthritis

Background

Genetic colocalization analysis is a statistical method that evaluates whether two traits (e.g., osteoarthritis [OA] risk and microRNA [miRNA] expression levels) share the same or distinct genetic association signals in a locus typically identified in genome-wide association studies (GWAS). This method is useful for providing insights into the biological relevance of genetic association signals, particularly in intergenic regions, which can help to elucidate disease mechanisms in OA and other complex traits.

Objectives

To review the existing literature on genetic colocalization methods, assess their suitability for studying OA, and investigate their capacity to integrate miRNA data, while bearing in view their statistical assumptions.

Design

We followed scoping review methodology and used Covidence software for data management. Search terms for colocalization, GWAS, and genetic or statistical models were used in the databases MEDLINE and EMBASE, searched till March 4, 2024.

Results

Our search returned 546 peer-reviewed papers, of which 96 were included following title/abstract and full-text screening. Based on both cumulative and annual publication counts, the most cited method for colocalization analysis was coloc. Four papers examined OA-related phenotypes, and none examined miRNA. An approach to colocalization analysis using miRNA was postulated based on further hand-searching.

Conclusions

Colocalization analysis is a largely unexplored method in OA. Many of the approaches to colocalization analysis identified in this review, including the integration of GWAS and miRNA data, may help to elucidate genetic and epigenetic factors implicated in OA and other complex traits.

Circulating miRNA-122 is associated with knee osteoarthritis progression: A 6-year longitudinal cohort study in the Yakumo study

OBJECTIVE

The association between knee osteoarthritis (OA) and miRNAs has been widely reported. However, the utility of miRNAs as predictors of knee osteoarthritis (KOA) progression in longitudinal studies has not been reported. We aimed to identify circulating miRNAs (c-miRNAs) associated with KOA progression in the general population and to examine their potential use as predictors of KOA progression.

METHODS

In 2012 and 2018, 66 participants (128 knees) took part in a resident health check-up in the Yakumo study. If the KL classification progressed two or more levels, the patient was classified as having progressive OA. Quantitative real-time polymerase chain reaction was used to screen 21 c-miRNAs. The expression levels of those c-miRNAs were compared between the progressive OA group and non-progressive OA group using student-t-test. Logistic analysis was performed in c-miRNAs less than p < 0.10 in univariate analysis.

RESULTS

The progressive OA group consisted of 78 knees. The results of the comparison between the progressive OA group and the non-progressive OA group showed that six c-miRNAs as follows; let7d (p = 0.030), c-miRNA-122 (p < 0.001), 150 (p = 0.070), 199 (p = 0.078), 21 (p = 0.016) and 320 (p = 0.093) were extracted as factors related to the progression of knee OA. In addition, logistic regression analysis identified c-miRNA-122 as an independent factor involved in the progression of knee osteoarthritis (odds ratio: 1.510, 95% confidence interval: 1.060-2.140, p = 0.023). The ROC curve showed by c-miRNA-122 for the progression of OA risk had an area under the curve of 0.702 (95% CI: 0.609-0.795). The threshold of c-miRNA-122 was -4.609.

CONCLUSION

The expression level of c-miRNA-122 was associated with the risk of KOA progression in community dwelling Japanese people.

Construction of an interactome network among circRNA-miRNA-mRNA reveals new biomarkers in hBMSCs osteogenic differentiation

Human bone marrow mesenchymal stem cells (hBMSCs) are adult stem cells residing in the bone marrow, characterized by their capacity for multi-directional differentiation, self-renewal, migration, and engraftment. Serving as seed cells, BMSCs play a pivotal role in the regeneration of bone defects. Hence, investigating the transcription factors and signaling pathways involved in the regulation of osteogenic differentiation in BMSCs holds significant importance. Recent research has unveiled that certain circular RNAs (circRNAs) can function as molecular sponges, influencing the osteogenic differentiation process of mesenchymal stem cells. However, many circRNAs remain undiscovered, and their precise mechanisms remain elusive. Therefore, the objective of this study is to construct an osteogenic differentiation-related circRNA-miRNA-mRNA network in hBMSCs. Subsequently, through bioinformatics analysis, we constructed a ceRNA network related to the osteogenic differentiation ability of hBMSCs, comprising 22 circRNAs, 17 miRNAs, and 15 mRNAs. The potential circRNA-miRNA-mRNA axes, including the role of hsa_circ_0001600 in promoting the osteogenic differentiation of hBMSCs through the targeted regulation of hsa-miR-542-3p, were validated through in vitro experiments.

Early microRNA and metabolite changes after anterior cruciate ligament reconstruction surgery

OBJECTIVES

Anterior cruciate ligament (ACL) reconstruction after injury does not prevent post-traumatic osteoarthritis (PTOA). Circulating microRNA (miRNA) and metabolite changes emerging shortly after ACL injury and reconstruction remain insufficiently defined, potentially harbouring early cues contributing to PTOA evolution. Moreover, their differential expression between females and males also may influence PTOA's natural trajectory. This study aims to determine alterations in plasma miRNA and metabolite levels in the early stages following ACL reconstruction and between females and males.

METHODS

A cohort of 43 ACL reconstruction patients was examined. Plasma was obtained at baseline, 2 weeks, and 6 weeks post-surgery (129 biospecimens in total). High-throughput miRNA sequencing and metabolomics were conducted. Differentially expressed miRNAs and metabolites were identified using negative binomial and linear regression models, respectively. Associations between miRNAs and metabolites were explored using time and sex as co-variants, (pre-surgery versus 2 and 6 weeks post-surgery). Using computational biology, miRNA-metabolite-gene interaction and pathway analyses were performed.

RESULTS

Levels of 46 miRNAs were increased at 2 weeks post-surgery compared to pre-surgery (baseline) using miRNA sequencing. Levels of 13 metabolites were significantly increased while levels of 6 metabolites were significantly decreased at 2 weeks compared to baseline using metabolomics. Hsa-miR-145-5p levels were increased in female subjects at both 2 weeks (log2-fold-change 0.71, 95%CI 0.22,1.20) and 6 weeks (log2-fold-change 0.75, 95%CI 0.07,1.43) post-surgery compared to males. In addition, hsa-miR-497-5p showed increased levels in females at 2 weeks (log2-fold-change 0.77, 95%CI 0.06,1.48) and hsa-miR-143-5p at 6 weeks (log2-fold-change 0.83, 95%CI 0.07,1.59). Five metabolites were decreased at 2 weeks post-surgery in females compared to males: L-leucine (-1.44, 95%CI -1.75,-1.13), g-guanidinobutyrate (-1.27, 95%CI 1.54,-0.99), creatinine (-1.17, 95%CI -1.44,-0.90), 2-methylbutyrylcarnitine (-1.76, 95%CI -2.17,-1.35), and leu-pro (-1.13, 95%CI -1.44,-0.83). MiRNA-metabolite-gene interaction analysis revealed key signalling pathways based on post-surgical time-point and in females versus males.

CONCLUSION

MiRNA and metabolite profiles were modified by time and by sex early after ACL reconstruction surgery, which could influence surgical response and ultimately risk of developing PTOA.

Evolution and advancements in genomics and epigenomics in OA research: How far we have come

OBJECTIVE

Osteoarthritis (OA) is the most prevalent musculoskeletal disease affecting articulating joint tissues, resulting in local and systemic changes that contribute to increased pain and reduced function. Diverse technological advancements have culminated in the advent of high throughput "omic" technologies, enabling identification of comprehensive changes in molecular mediators associated with the disease. Amongst these technologies, genomics and epigenomics - including methylomics and miRNomics, have emerged as important tools to aid our biological understanding of disease.

DESIGN

In this narrative review, we selected articles discussing advancements and applications of these technologies to OA biology and pathology. We discuss how genomics, deoxyribonucleic acid (DNA) methylomics, and miRNomics have uncovered disease-related molecular markers in the local and systemic tissues or fluids of OA patients.

RESULTS

Genomics investigations into the genetic links of OA, including using genome-wide association studies, have evolved to identify 100+ genetic susceptibility markers of OA. Epigenomic investigations of gene methylation status have identified the importance of methylation to OA-related catabolic gene expression. Furthermore, miRNomic studies have identified key microRNA signatures in various tissues and fluids related to OA disease.

CONCLUSIONS

Sharing of standardized, well-annotated omic datasets in curated repositories will be key to enhancing statistical power to detect smaller and targetable changes in the biological signatures underlying OA pathogenesis. Additionally, continued technological developments and analysis methods, including using computational molecular and regulatory networks, are likely to facilitate improved detection of disease-relevant targets, in-turn, supporting precision medicine approaches and new treatment strategies for OA.

Neuronal small extracellular vesicles carrying miR-181c-5p contribute to the pathogenesis of epilepsy by regulating the protein kinase C-δ/glutamate transporter-1 axis in astrocytes

Information exchange between neurons and astrocytes mediated by extracellular vesicles (EVs) is known to play a key role in the pathogenesis of central nervous system diseases. A key driver of epilepsy is the dysregulation of intersynaptic excitatory neurotransmitters mediated by astrocytes. Thus, we investigated the potential association between neuronal EV microRNAs (miRNAs) and astrocyte glutamate uptake ability in epilepsy. Here, we showed that astrocytes were able to engulf epileptogenic neuronal EVs, inducing a significant increase in the glutamate concentration in the extracellular fluid of astrocytes, which was linked to a decrease in glutamate transporter-1 (GLT-1) protein expression. Using sequencing and gene ontology (GO) functional analysis, miR-181c-5p was found to be the most significantly upregulated miRNA in epileptogenic neuronal EVs and was linked to glutamate metabolism. Moreover, we found that neuronal EV-derived miR-181c-5p interacted with protein kinase C-delta (PKCδ), downregulated PKCδ and GLT-1 protein expression and increased glutamate concentrations in astrocytes both in vitro and in vivo. Our findings demonstrated that epileptogenic neuronal EVs carrying miR-181c-5p decrease the glutamate uptake ability of astrocytes, thus promoting susceptibility to epilepsy.

Three decades of advancements in osteoarthritis research: insights from transcriptomic, proteomic, and metabolomic studies

OBJECTIVE

Osteoarthritis (OA) is a complex disease involving contributions from both local joint tissues and systemic sources. Patient characteristics, encompassing sociodemographic and clinical variables, are intricately linked with OA rendering its understanding challenging. Technological advancements have allowed for a comprehensive analysis of transcripts, proteomes and metabolomes in OA tissues/fluids through omic analyses. The objective of this review is to highlight the advancements achieved by omic studies in enhancing our understanding of OA pathogenesis over the last three decades.

DESIGN

We conducted an extensive literature search focusing on transcriptomics, proteomics and metabolomics within the context of OA. Specifically, we explore how these technologies have identified individual transcripts, proteins, and metabolites, as well as distinctive endotype signatures from various body tissues or fluids of OA patients, including insights at the single-cell level, to advance our understanding of this highly complex disease.

RESULTS

Omic studies reveal the description of numerous individual molecules and molecular patterns within OA-associated tissues and fluids. This includes the identification of specific cell (sub)types and associated pathways that contribute to disease mechanisms. However, there remains a necessity to further advance these technologies to delineate the spatial organization of cellular subtypes and molecular patterns within OA-afflicted tissues.

CONCLUSIONS

Leveraging a multi-omics approach that integrates datasets from diverse molecular detection technologies, combined with patients' clinical and sociodemographic features, and molecular and regulatory networks, holds promise for identifying unique patient endophenotypes. This holistic approach can illuminate the heterogeneity among OA patients and, in turn, facilitate the development of tailored therapeutic interventions.

Association of synovial fluid and urinary C2C-HUSA levels with surgical outcomes post-total knee arthroplasty

OBJECTIVES

After total knee arthroplasty (TKA), ∼30% of knee osteoarthritis (KOA) patients show little symptomatic improvement. Earlier studies have correlated urinary (u) type 2 collagen C terminal cleavage peptide assay (C2C-HUSA), which detects a fragment of cartilage collagen breakdown, with KOA progression. This study determines whether C2C levels in urine, synovial fluid, or their ratio, are associated with post-surgical outcomes.

METHODS

From a large sample of 489 subjects, diagnosed with primary KOA undergoing TKA, Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) pain and function scores were collected at baseline (time of surgery) and one-year post-TKA. Baseline urine (u) and synovial fluid (sf) were analysed using the IBEX-C2C-HUSA assay, with higher values indicating higher amounts of cartilage degradation. For urine, results were normalised to creatinine. Furthermore, subjects' changes in WOMAC scores were categorised based on percent reduction in pain or improvement in function, compared to baseline, such that >66.7%, >33.3 to ≤66.7%, and ≤33.3% denoted "strong", "moderate" and "mild/worse" responses, respectively. Associations of individual biofluid C2C-HUSA levels, or their ratio, with change in WOMAC pain and function scores up to one-year post-TKA, or category of change, were analysed by linear, logistic, or cumulative odds models.

RESULTS

Higher baseline uC2C-HUSA levels or a lower ratio of baseline sfC2C-HUSA to uC2C-HUSA were associated with improvements in WOMAC pain by linear multivariable modelling [odds ratio -0.40 (95% confidence interval -0.76, -0.05) p = 0.03; 0.36 (0.01, 0.71), p = 0.04, respectively], while sfC2C-HUSA alone was not. However, lower ratios of sfC2C-HUSA to uC2C-HUSA were associated with improvements in WOMAC function [1.37 (0.18, 2.55), p = 0.02], while sfC2C-HUSA and uC2C-HUSA alone were not. Lower ratios of sfC2C-HUSA to uC2C-HUSA were also associated with an increased likelihood of a subject being categorised in a group where TKA was beneficial in both univariable [pain, 0.81 (0.68, 0.96), p = 0.02; function, 0.92 (0.85, 0.99), p = 0.035] and multivariable [pain, 0.81 (0.68, 0.97) p = 0.02; function, 0.92 (0.85, 1.00), p = 0.043] ordinal modelling, while sfC2C-HUSA and uC2C-HUSA alone were not.

CONCLUSIONS

Overall, ratios of baseline sfC2C-HUSA to uC2C-HUSA, and baseline uC2C-HUSA, may play an important role in studying post-TKA surgical outcomes.

ECM1-associated miR-1260b promotes osteogenic differentiation by targeting GDI1

Osteoporosis (OP) is a common disease among older adults. The promotion of osteoblast differentiation plays a crucial role in alleviating OP symptoms. Extracellular matrix protein 1 (ECM1) has been reported to be closely associated with osteogenic differentiation. In this study, we constructed U2OS cell lines with ECM1 knockdown and ECM1a overexpression based on knockdown, and identified the target miRNA (miR-1260b) by sequencing. Overexpression of miR-1260b promoted the osteogenic differentiation of U2OS and MG63 cells, as demonstrated by increased alkaline phosphatase (ALP) activity, matrix mineralization, and Runt-Related Transcription Factor 2 (RUNX2), Osteopontin (OPN), Collagen I (COL1A1), and Osteocalcin (OCN) protein expressions, whereas low expression of miR-1260b had the opposite effect. In addition, miR-1260b expression was decreased in OP patients than in non-OP patients. Next, we predicted the target gene of miRNA through TargetScan and miRDB and found that miR-1260b negatively regulated GDP dissociation inhibitor 1 (GDI1) by directly binding to its 3'-untranslated region. Subsequent experiments revealed that GDI1 overexpression decreased ALP activity and calcium deposit, reduced RUNX2, OPN, COL1A1, and OCN expression levels, and reversed the effects of miR-1260b on osteogenic differentiation. In conclusion, ECM1-related miR-1260b promotes osteogenic differentiation by targeting GDI1 in U2OS and MG63 cells. Thus, this study has significant implication for osteoporosis treatment.

Proteomic and genomic profiling of plasma exosomes from patients with ankylosing spondylitis

INTRODUCTION

Recent advances in understanding the biology of ankylosing spondylitis (AS) using innovative genomic and proteomic approaches offer the opportunity to address current challenges in AS diagnosis and management. Altered expression of genes, microRNAs (miRNAs) or proteins may contribute to immune dysregulation and may play a significant role in the onset and persistence of inflammation in AS. The ability of exosomes to transport miRNAs across cells and alter the phenotype of recipient cells has implicated exosomes in perpetuating inflammation in AS. This study reports the first proteomic and miRNA profiling of plasma-derived exosomes in AS using comprehensive computational biology analysis.

METHODS

Plasma samples from patients with AS and healthy controls (HC) were isolated via ultracentrifugation and subjected to extracellular vesicle flow cytometry analysis to characterise exosome surface markers by a multiplex immunocapture assay. Cytokine profiling of plasma-derived exosomes and cell culture supernatants was performed. Next-generation sequencing was used to identify miRNA populations in exosomes enriched from plasma fractions. CD4+ T cells were sorted, and the frequency and proliferation of CD4+ T-cell subsets were analysed after treatment with AS-exosomes using flow cytometry.

RESULTS

The expression of exosome marker proteins CD63 and CD81 was elevated in the patients with AS compared with HC (q<0.05). Cytokine profiling in plasma-derived AS-exosomes demonstrated downregulation of interleukin (IL)-8 and IL-10 (q<0.05). AS-exosomes cocultured with HC CD4+ T cells induced significant upregulation of IFNα2 and IL-33 (q<0.05). Exosomes from patients with AS inhibited the proliferation of regulatory T cells (Treg), suggesting a mechanism for chronically activated T cells in this disease. Culture of CD4+ T cells from healthy individuals in the presence of AS-exosomes reduced the proliferation of FOXP3+ Treg cells and decreased the frequency of FOXP3+IRF4+ Treg cells. miRNA sequencing identified 24 differentially expressed miRNAs found in circulating exosomes of patients with AS compared with HC; 22 of which were upregulated and 2 were downregulated.

CONCLUSIONS

Individuals with AS have different immunological and genetic profiles, as determined by evaluating the exosomes of these patients. The inhibitory effect of exosomes on Treg in AS suggests a mechanism contributing to chronically activated T cells in this disease.

MicroRNA as Possible Mediators of the Synergistic Effect of Celecoxib and Glucosamine Sulfate in Human Osteoarthritic Chondrocyte Exposed to IL-1β

This study investigated the role of a pattern of microRNA (miRNA) as possible mediators of celecoxib and prescription-grade glucosamine sulfate (GS) effects in human osteoarthritis (OA) chondrocytes. Chondrocytes were treated with celecoxib (1.85 µM) and GS (9 µM), alone or in combination, for 24 h, with or without interleukin (IL)-1β (10 ng/mL). Cell viability was determined using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay, apoptosis and reactive oxygen species (ROS) by cytometry, nitric oxide (NO) by Griess method. Gene levels of miRNA, antioxidant enzymes, nuclear factor erythroid (NRF)2, and B-cell lymphoma (BCL)2 expressions were analyzed by quantitative real time polymerase chain reaction (real time PCR). Protein expression of NRF2 and BCL2 was also detected at immunofluorescence and western blot. Celecoxib and GS, alone or in combination, significantly increased viability, reduced apoptosis, ROS and NO production and the gene expression of miR-34a, -146a, -181a, -210, in comparison to baseline and to IL-1β. The transfection with miRNA specific inhibitors significantly counteracted the IL-1β activity and potentiated the properties of celecoxib and GS on viability, apoptosis and oxidant system, through nuclear factor (NF)-κB regulation. The observed effects were enhanced when the drugs were tested in combination. Our data confirmed the synergistic anti-inflammatory and chondroprotective properties of celecoxib and GS, suggesting microRNA as possible mediators.

Emerging molecular biomarkers in osteoarthritis pathology

Osteoarthritis (OA) is the most common form of arthritis resulting in joint discomfort and disability, culminating in decline in life quality. Attention has been drawn in recent years to disease-associated molecular biomarkers found in readily accessible biofluids due to low invasiveness of acquisition and their potential to detect early pathological molecular changes not observed with traditional imaging methodology. These biochemical markers of OA have been found in synovial fluid, blood, and urine. They include emerging molecular classes, such as metabolites and noncoding RNAs, as well as classical biomarkers, like inflammatory mediators and by-products of degradative processes involving articular cartilage. Although blood-based biomarkers tend to be most studied, the use of synovial fluid, a more isolated biofluid in the synovial joint, and urine as an excreted fluid containing OA biomarkers can offer valuable information on local and overall disease activity, respectively. Furthermore, larger clinical studies are required to determine relationships between biomarkers in different biofluids, and their impacts on patient measures of OA. This narrative review provides a concise overview of recent studies of OA using these four classes of biomarkers as potential biomarker for measuring disease incidence, staging, prognosis, and therapeutic intervention efficacy.

Characterization of miR-335-5p and miR-335-3p in human osteoarthritic tissues

OBJECTIVE

We aimed to characterize the expression patterns, gene targets, and functional effects of miR-335-5p and miR-335-3p among seven primary human knee and hip osteoarthritic tissue types.

METHODS

We collected synovial fluid, subchondral bone, articular cartilage, synovium, meniscus/labrum, infrapatellar/acetabular fat, anterior cruciate ligament/ligamentum teres, and vastus medialis oblique/quadratus femoris muscle (n = 7-20) from surgical patients with early- or late-stage osteoarthritis (OA) and quantified miR-335-5p and miR-335-3p expression by real-time PCR. Predicted gene targets were measured in knee OA infrapatellar fat following miRNA inhibitor transfection (n = 3), and prioritized gene targets were validated following miRNA inhibitor and mimic transfection (n = 6). Following pathway analyses, we performed Oil-Red-O staining to assess changes in total lipid content in infrapatellar fat.

RESULTS

Showing a 227-fold increase in knee OA infrapatellar fat (the highest expressing tissue) versus meniscus (the lowest expressing tissue), miR-335-5p was more abundant than miR-335-3p (92-fold increase). MiR-335-5p showed higher expression across knee tissues versus hip tissues, and in late-stage versus early-stage knee OA fat. Exploring candidate genes, VCAM1 and MMP13 were identified as putative direct targets of miR-335-5p and miR-335-3p, respectively, showing downregulation with miRNA mimic transfection. Exploring candidate pathways, predicted miR-335-5p gene targets were enriched in a canonical adipogenesis network (p = 2.1e - 5). Modulation of miR-335-5p in late-stage knee OA fat showed an inverse relationship to total lipid content.

CONCLUSION

Our data suggest both miR-335-5p and miR-335-3p regulate gene targets in late-stage knee OA infrapatellar fat, though miR-335-5p appears to be more prominent, with tissue-, joint-, and stage-specific effects.

CircRNAs in osteoarthritis: research status and prospect

Osteoarthritis (OA) is the most common joint disease globally, and its progression is irreversible. The mechanism of osteoarthritis is not fully understood. Research on the molecular biological mechanism of OA is deepening, among which epigenetics, especially noncoding RNA, is an emerging hotspot. CircRNA is a unique circular noncoding RNA not degraded by RNase R, so it is a possible clinical target and biomarker. Many studies have found that circRNAs play an essential role in the progression of OA, including extracellular matrix metabolism, autophagy, apoptosis, the proliferation of chondrocytes, inflammation, oxidative stress, cartilage development, and chondrogenic differentiation. Differential expression of circRNAs was also observed in the synovium and subchondral bone in the OA joint. In terms of mechanism, existing studies have mainly found that circRNA adsorbs miRNA through the ceRNA mechanism, and a few studies have found that circRNA can serve as a scaffold for protein reactions. In terms of clinical transformation, circRNAs are considered promising biomarkers, but no large cohort has tested their diagnostic value. Meanwhile, some studies have used circRNAs loaded in extracellular vesicles for OA precision medicine. However, there are still many problems to be solved in the research, such as the role of circRNA in different OA stages or OA subtypes, the construction of animal models of circRNA knockout, and more research on the mechanism of circRNA. In general, circRNAs have a regulatory role in OA and have particular clinical potential, but further studies are needed in the future.

DraculR: A Web-Based Application for In Silico Haemolysis Detection in High-Throughput microRNA Sequencing Data

The search for novel microRNA (miRNA) biomarkers in plasma is hampered by haemolysis, the lysis and subsequent release of red blood cell contents, including miRNAs, into surrounding fluid. The biomarker potential of miRNAs comes in part from their multicompartment origin and the long-lived nature of miRNA transcripts in plasma, giving researchers a functional window for tissues that are otherwise difficult or disadvantageous to sample. The inclusion of red-blood-cell-derived miRNA transcripts in downstream analysis introduces a source of error that is difficult to identify posthoc and may lead to spurious results. Where access to a physical specimen is not possible, our tool will provide an in silico approach to haemolysis prediction. We present DraculR, an interactive Shiny/R application that enables a user to upload miRNA expression data from a short-read sequencing of human plasma as a raw read counts table and interactively calculate a metric that indicates the degree of haemolysis contamination. The code, DraculR web tool and its tutorial are freely available as detailed herein.

miR-199a-5p Reduces Chondrocyte Hypertrophy and Attenuates Osteoarthritis Progression via the Indian Hedgehog Signal Pathway

Osteoarthritis (OA), the most common type of arthritis, is an age-associated disease, characterized by the progressive degradation of articular cartilage, synovial inflammation, and degeneration of subchondral bone. Chondrocyte proliferation is regulated by the Indian hedgehog (IHH in humans, Ihh in animals) signaling molecule, which regulates hypertrophy and endochondral ossification in the development of the skeletal system. microRNAs (miRNAs, miRs) are a family of about 22-nucleotide endogenous non-coding RNAs, which negatively regulate gene expression. In this study, the expression level of IHH was upregulated in the damaged articular cartilage tissues among OA patients and OA cell cultures, while that of miR-199a-5p was the opposite. Further investigations demonstrated that miR-199a-5p could directly regulate IHH expression and reduce chondrocyte hypertrophy and matrix degradation via the IHH signal pathway in the primary human chondrocytes. The intra-articular injection of synthetic miR-199a-5p agomir attenuated OA symptoms in rats, including the alleviation of articular cartilage destruction, subchondral bone degradation, and synovial inflammation. The miR-199a-5p agomir could also inhibit the Ihh signaling pathway in vivo. This study might help in understanding the role of miR-199a-5p in the pathophysiology and molecular mechanisms of OA and indicate a potential novel therapeutic strategy for OA patients.

Research progress on the role and mechanism of miR-671 in bone metabolism and bone-related diseases

Bone metabolism consists of bone formation and resorption and maintains a dynamic balance in vivo. When bone homeostasis is broken, it can manifest as osteoarthritis (OA), rheumatoid arthritis (RA), osteosarcoma (OS), etc. MiR-671, an important class of non-coding nucleotide sequences in vivo, is regulated by lncRNA and regulates bone metabolism balance by regulating downstream target proteins and activating various signaling pathways. Based on the structure and primary function of miR-671, this paper summarizes the effect and mechanism of miR-671 in bone-related inflammation and cancer diseases, and prospects the application possibility of miR-671, providing reference information for targeted therapy of bone-related disorders.

A review on the role of miR-671 in human disorders

miR-671 is encoded by a gene on 7q36.1 and contributes to the pathogenesis of a variety of disorders, including diverse types of cancers, atherosclerosis, ischemic stroke, liver fibrosis, osteoarthritis, Parkinson's disease, rheumatoid arthritis, acute myocardial infarction and Crohn's disease. In the context of cancer, different studies have revealed opposite roles for this miRNA. In brief, it has been shown to be down-regulated in pancreatic ductal carcinoma, ovarian cancer, gastric cancer, osteosarcoma, esophageal squamous cell carcinoma and myelodysplastic syndromes. Yet, miR-671 has been up-regulated in glioma, colorectal cancer, prostate cancer and hepatocellular carcinoma. Studies in breast, lung and renal cell carcinoma have reported inconsistent results. The current review aims at summarization of the role of miR-671 in these disorders focusing on its target mRNA in each context and dysregulated signaling pathways. We also provide a summary of the role of this miRNA as a prognostic factor in malignancies.

MirDIP 5.2: tissue context annotation and novel microRNA curation

MirDIP is a well-established database that aggregates microRNA-gene human interactions from multiple databases to increase coverage, reduce bias, and improve usability by providing an integrated score proportional to the probability of the interaction occurring. In version 5.2, we removed eight outdated resources, added a new resource (miRNATIP), and ran five prediction algorithms for miRBase and mirGeneDB. In total, mirDIP 5.2 includes 46 364 047 predictions for 27 936 genes and 2734 microRNAs, making it the first database to provide interactions using data from mirGeneDB. Moreover, we curated and integrated 32 497 novel microRNAs from 14 publications to accelerate the use of these novel data. In this release, we also extend the content and functionality of mirDIP by associating contexts with microRNAs, genes, and microRNA-gene interactions. We collected and processed microRNA and gene expression data from 20 resources and acquired information on 330 tissue and disease contexts for 2657 microRNAs, 27 576 genes and 123 651 910 gene-microRNA-tissue interactions. Finally, we improved the usability of mirDIP by enabling the user to search the database using precursor IDs, and we integrated miRAnno, a network-based tool for identifying pathways linked to specific microRNAs. We also provide a mirDIP API to facilitate access to its integrated predictions. Updated mirDIP is available at https://ophid.utoronto.ca/mirDIP.

A pilot study of microRNA assessment as a means to identify novel biomarkers of spontaneous osteoarthritis in dogs

MicroRNAs (miRNAs) are important regulators of intercellular signaling and are promising biomarkers in osteoarthritis (OA). In this study, comprehensive analysis was performed to identify miRNAs involved in the pathogenesis of spontaneous OA in dogs. Dogs diagnosed with OA based on radiography and arthroscopy of the stifle joint were included in the OA group. Dogs without any evidence of orthopedic disease were included in the unaffected group. To investigate miRNA expression levels, RNA sequencing analysis (RNA-seq) was performed in synovial tissue (OA group: n = 3, Unaffected group: n = 3) and RT-qPCR was performed in synovial tissue, synovial fluid and serum (OA group: n = 17, Unaffected group: n = 6), and compared between the two groups. The RNA-seq results showed that 57 miRNAs were significantly upregulated and 42 were significantly downregulated in the OA group. Specifically, miR-542 and miR-543 expression levels in the synovial tissue, synovial fluid, and serum were consistently higher in the OA group than in the unaffected group, suggesting that these miRNAs may be used as biomarkers for detecting canine OA. This is the first report to comprehensively analyze the expression patterns of miRNAs in the synovial tissue of dogs with spontaneous OA.

Association of presurgical circulating MicroRNAs with 1-year postsurgical pain reduction in spine facet osteoarthritis patients with lumbar spinal stenosis

Purpose

Up to 30% of spine facet osteoarthritis patients with lumbar spinal stenosis (SF-OA ​+ ​LSS) have little to no improvement in their pain after surgery. Lack of meaningful improvement in pain following surgery provides a unique opportunity to identify specific predictive biomarker signatures that might be associated with the outcomes of surgical treatment. The objective of the present study was to determine whether a microRNA (miRNA) biomarker signature could be identified in presurgical blood plasma that corresponded with levels of SF-OA ​+ ​LSS patient post-surgical pain intensity one year later.

Methods

RNA was extracted from baseline plasma of SF-OA ​+ ​LSS patients and prepared for miRNA sequencing. Statistical approaches were performed to identify differentially expressed miRNAs associated with reduced 1-year postsurgical pain (n ​= ​56). Using an integrated computational approach, we further created predicted gene and pathway networks for each identified miRNA.

Results

We identified a panel of 4 circulating candidate miRNAs (hsa-miR-155-5p, hsa-let-7e-5p, hsa-miR-125a-5p, hsa-miR-99b-5p) with higher levels at presurgical baseline that were associated with greater changes in % NPRS20Δ, reflecting reduced pain intensity levels at one year. Genes encoding hsa-let-7e-5p, hsa-miR-125a-5p, and hsa-miR-99b-5p are part of an evolutionarily conserved miRNA cluster. Using integrated computational analyses, we showed that mammalian target of rapamycin, transforming growth factor-β1 receptor, Wnt signaling, epithelial-mesenchymal transition regulators, and cholecystokinin signaling were enriched pathways of predicted gene targets.

Conclusions

Taken together, our findings suggest that 4 presurgical baseline circulating miRNAs correlate with 1-year postsurgical SF-OA ​+ ​LSS patient pain intensity and represent possible candidate biomarker signature of surgical pain response.

Small non-coding RNA landscape of extracellular vesicles from a post-traumatic model of equine osteoarthritis

Extracellular vesicles comprise an as yet inadequately investigated intercellular communication pathway in the field of early osteoarthritis. We hypothesised that the small non-coding RNA expression pattern in synovial fluid and plasma would change during progression of experimental osteoarthritis. In this study, we conducted small RNA sequencing to provide a comprehensive overview of the temporal expression profiles of small non-coding transcripts carried by extracellular vesicles derived from plasma and synovial fluid for the first time in a posttraumatic model of equine osteoarthritis. Additionally, we characterised synovial fluid and plasma-derived extracellular vesicles with respect to quantity, size, and surface markers. The different temporal expressions of seven microRNAs in plasma and synovial fluid-derived extracellular vesicles, eca-miR-451, eca-miR-25, eca-miR-215, eca-miR-92a, eca-miR-let-7c, eca-miR-486-5p, and eca-miR-23a, and four snoRNAs, U3, snord15, snord46, and snord58, represent potential biomarkers for early osteoarthritis. Bioinformatics analysis of the differentially expressed microRNAs in synovial fluid highlighted that in early osteoarthritis these related to the inhibition of cell cycle, cell cycle progression, DNA damage and cell proliferation as well as increased cell viability and differentiation of stem cells. Plasma and synovial fluid-derived extracellular vesicle small non-coding signatures have been established for the first time in a temporal model of osteoarthritis. These could serve as novel biomarkers for evaluation of osteoarthritis progression or act as potential therapeutic targets.

Exosomes rewire the cartilage microenvironment in osteoarthritis: from intercellular communication to therapeutic strategies

Osteoarthritis (OA) is a prevalent degenerative joint disease characterized by cartilage loss and accounts for a major source of pain and disability worldwide. However, effective strategies for cartilage repair are lacking, and patients with advanced OA usually need joint replacement. Better comprehending OA pathogenesis may lead to transformative therapeutics. Recently studies have reported that exosomes act as a new means of cell-to-cell communication by delivering multiple bioactive molecules to create a particular microenvironment that tunes cartilage behavior. Specifically, exosome cargos, such as noncoding RNAs (ncRNAs) and proteins, play a crucial role in OA progression by regulating the proliferation, apoptosis, autophagy, and inflammatory response of joint cells, rendering them promising candidates for OA monitoring and treatment. This review systematically summarizes the current insight regarding the biogenesis and function of exosomes and their potential as therapeutic tools targeting cell-to-cell communication in OA, suggesting new realms to improve OA management.

Haemolysis Detection in MicroRNA-Seq from Clinical Plasma Samples

The abundance of cell-free microRNA (miRNA) has been measured in blood plasma and proposed as a source of novel, minimally invasive biomarkers for several diseases. Despite improvements in quantification methods, there is no consensus regarding how haemolysis affects plasma miRNA content. We propose a method for haemolysis detection in miRNA high-throughput sequencing (HTS) data from libraries prepared using human plasma. To establish a miRNA haemolysis signature we tested differential miRNA abundance between plasma samples with known haemolysis status. Using these miRNAs with statistically significant higher abundance in our haemolysed group, we further refined the set to reveal high-confidence haemolysis association. Given our specific context, i.e., women of reproductive age, we also tested for significant differences between pregnant and non-pregnant groups. We report a novel 20-miRNA signature used to identify the presence of haemolysis in silico in HTS miRNA-sequencing data. Further, we validated the signature set using firstly an all-male cohort (prostate cancer) and secondly a mixed male and female cohort (radiographic knee osteoarthritis). Conclusion: Given the potential for haemolysis contamination, we recommend that assays for haemolysis detection become standard pre-analytical practice and provide here a simple method for haemolysis detection.

Pentraxin 3 regulated by miR-224-5p modulates macrophage reprogramming and exacerbates osteoarthritis associated synovitis by targeting CD32

Emerging evidence has shown an imbalance in M1/M2 macrophage polarization to play an essential role in osteoarthritis (OA) progression. However, the underlying mechanistic basis for this polarization is unknown. RNA sequencing of OA M1-polarized macrophages found highly expressed levels of pentraxin 3 (PTX3), suggesting a role for PTX3 in OA occurrence and development. Herein, PTX3 was found to be increased in the synovium and articular cartilage of OA patients and OA mice. Intra-articular injection of PTX3 aggravated, while PTX3 neutralization reversed synovitis and cartilage degeneration. No metabolic disorder or proteoglycan loss were observed in cartilage explants when treated with PTX3 alone. However, cartilage explants exhibited an OA phenotype when treated with culture supernatants of macrophages stimulated with PTX3, suggesting that PTX3 did not have a direct effect on chondrocytes. Therefore, the OA anti-chondrogenic effects of PTX3 are primarily mediated through macrophages. Mechanistically, PTX3 was upregulated by miR-224-5p deficiency, which activated the p65/NF-κB pathway to promote M1 macrophage polarization by targeting CD32. CD32 was expressed by macrophages, that when stimulated with PTX3, secreted abundant pro-inflammation cytokines that induced severe articular cartilage damage. The paracrine interaction between macrophages and chondrocytes produced a feedback loop that enhanced synovitis and cartilage damage. The findings of this study identified a functional pathway important to OA development. Blockade of this pathway and PTX3 may prevent and treat OA.

Circulating microRNAs differentiate fast-progressing from slow-progressing and non-progressing knee osteoarthritis in the Osteoarthritis Initiative cohort

Introduction

The objective of this study is to identify circulating microRNAs that distinguish fast-progressing radiographic knee osteoarthritis (OA) in the Osteoarthritis Initiative cohort by applying microRNA-sequencing.

Methods

Participants with Kellgren-Lawrence (KL) grade 0/1 at baseline were included (N = 106). Fast-progressors were defined by an increase to KL 3/4 by 4-year follow-up (N = 20), whereas slow-progressors showed an increase to KL 2/3/4 only at 8-year follow-up (N = 35). Non-progressors remained at KL 0/1 by 8-year follow-up (N = 51). MicroRNA-sequencing was performed on plasma collected at baseline and 4-year follow-up from the same participants. Negative binomial models were fitted to identify differentially expressed (DE) microRNAs. Penalized logistic regression (PLR) analyses were performed to select combinations of DE microRNAs that distinguished fast-progressors. Area under the receiver operating characteristic curves (AUC) were constructed to evaluate predictive ability.

Results

DE analyses revealed 48 microRNAs at baseline and 2 microRNAs at 4-year follow-up [false discovery rate (FDR) < 0.05] comparing fast-progressors with both slow-progressors and non-progressors. Among these were hsa-miR-320b, hsa-miR-320c, hsa-miR-320d, and hsa-miR-320e, which were predicted to target gene families, including members of the 14-3-3 gene family, involved in signal transduction. PLR models included miR-320 members as top predictors of fast-progressors and yielded AUC ranging from 82.6 to 91.9, representing good accuracy.

Conclusion

The miR-320 family is associated with fast-progressing radiographic knee OA and merits further investigation as potential biomarkers and mechanistic drivers of knee OA.

MiR-1260b protects against LPS-induced degenerative changes in nucleus pulposus cells through targeting TCF7L2

Nucleus pulposus (NP) cells play a critical role in maintaining intervertebral disc integrity through producing the components of extracellular matrix (ECM). NP cell dysfunction, including senescence and hyper-apoptosis, has been regarded as critical events during intervertebral disc degeneration development. In the present study, we found that Transcription Factor 7-Like 2 (TCF7L2) was overexpressed within degenerative intervertebral disc tissue samples, and TCF7L2 silencing improved lipopolysaccharide (LPS)-induced repression on NP cell proliferation, ECM synthesis, and LPS-induced NP cell senescence. miR-1260b directly targeted TCF7L2 and inhibited TCF7L2 expression. miR-1260b overexpression improved LPS-induced degenerative changes in NP cells; more importantly, TCF7L2 overexpression significantly reversed the effects of miR-1260b overexpression on LPS-stimulated degenerative changes within NP cells. For the first time, we demonstrated the function of the miR-1260b/TCF7L2 axis on the phenotypic maintenance of chondrocyte-like NP cells and ECM synthesis by NP cells under LPS stimulation.

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.

Genetic biomarkers in osteoarthritis: a quick overview

Osteoarthritis (OA) is a chronic musculoskeletal disease with a polygenic and heterogeneous nature. In addition, when clinical manifestations appear, the evolution of the disease is usually already irreversible. Therefore, the efforts on OA research are focused mainly on the discovery of therapeutic targets and reliable biomarkers that permit the early identification of different OA-related parameters such as diagnosis, prognosis, or phenotype identification. To date, potential candidate protein biomarkers have been associated with different aspects of the disease; however, there is currently no gold standard. In this sense, genomic data could act as complementary biomarkers of diagnosis and prognosis or even help to identify therapeutic targets of the disease. In this review, we will describe the most recent advances in genetic biomarkers in OA over the past three years.

Mouse microRNA signatures in joint ageing and post-traumatic osteoarthritis

OBJECTIVE

This study investigated mice serum and joint microRNA expression profiles in ageing and osteoarthritis to elucidate the role of microRNAs in the development and progression of disease, and provide biomarkers for ageing and osteoarthritis.

DESIGN

Whole joints and serum samples were collected from C57BL6/J male mice and subjected to small RNA sequencing. Groups used included; surgically-induced post-traumatic osteoarthritis, (DMM; 24 months-old); sham surgery (24 months-old); old mice (18 months-old); and young mice (8 months-old). Differentially expressed microRNAs between the four groups were identified and validated using real-time quantitative PCR. MicroRNA differential expression data was used for target prediction and pathway analysis.

RESULTS

In joint tissues, miR-140-5p, miR-205-5p, miR-682, miR-208b-3p, miR-499-5p, miR-455-3p and miR-6238 were differentially expressed between young and old groups; miR-146a-5p, miR-3474, miR-615-3p and miR-151-5p were differentially expressed between DMM and Sham groups; and miR-652-3p, miR-23b-3p, miR-708-5p, miR-5099, miR-23a-3p, miR-214-3p, miR-6238 and miR-148-3p between the old and DMM groups. The number of differentially expressed microRNAs in serum was higher, some in common with joint tissues including miR-140-5p and miR-455-3p between young and old groups; and miR-23b-3p, miR-5099 and miR-6238 between old and DMM groups.We confirmed miR-140-5p, miR-499-5p and miR-455-3p expression to be decreased in old mouse joints compared to young, suggesting their potential use as biomarkers of joint ageing in mice.

CONCLUSIONS

MiR-140-5p, miR-499-5p and miR-455-3p could be used as joint ageing biomarkers in mice. Further research into these specific molecules in human tissues is now warranted to check their potential suitability as human biomarkers of ageing.

DNA methylation and noncoding RNA in OA: Recent findings and methodological advances

Introduction

Osteoarthritis (OA) is a chronic musculoskeletal disease characterized by progressive loss of joint function. Historically, it has been characterized as a disease caused by mechanical trauma, so-called 'wear and tear'. Over the past two decades, it has come to be understood as a complex systemic disorder involving gene-environmental interactions. Epigenetic changes have been increasingly implicated. Recent improvements in microarray and next-generation sequencing (NGS) technologies have allowed for ever more complex evaluations of epigenetic aberrations associated with the development and progression of OA.

Methods

A systematic review was conducted in the Pubmed database. We curated studies that presented the results of DNA methylation and noncoding RNA research in human OA and OA animal models since 1985.

Results

Herein, we discuss recent findings and methodological advancements in OA epigenetics, including a discussion of DNA methylation, including microarray and NGS studies, and noncoding RNAs. Beyond cartilage, we also highlight studies in subchondral bone and peripheral blood mononuclear cells, which highlight widespread and potentially clinically important alterations in epigenetic patterns seen in OA patients. Finally, we discuss epigenetic editing approaches in the context of OA.

Conclusions

Although a substantial body of literature has already been published in OA, much is still unknown. Future OA epigenetics studies will no doubt continue to broaden our understanding of underlying pathophysiology and perhaps offer novel diagnostics and/or treatments for human OA.

The non-coding RNA interactome in joint health and disease

Non-coding RNAs have distinct regulatory roles in the pathogenesis of joint diseases including osteoarthritis (OA) and rheumatoid arthritis (RA). As the amount of high-throughput profiling studies and mechanistic investigations of microRNAs, long non-coding RNAs and circular RNAs in joint tissues and biofluids has increased, data have emerged that suggest complex interactions among non-coding RNAs that are often overlooked as critical regulators of gene expression. Identifying these non-coding RNAs and their interactions is useful for understanding both joint health and disease. Non-coding RNAs regulate signalling pathways and biological processes that are important for normal joint development but, when dysregulated, can contribute to disease. The specific expression profiles of non-coding RNAs in various disease states support their roles as promising candidate biomarkers, mediators of pathogenic mechanisms and potential therapeutic targets. This Review synthesizes literature published in the past 2 years on the role of non-coding RNAs in OA and RA with a focus on inflammation, cell death, cell proliferation and extracellular matrix dysregulation. Research to date makes it apparent that 'non-coding' does not mean 'non-essential' and that non-coding RNAs are important parts of a complex interactome that underlies OA and RA.

Characterization and miRNA Profiling of Extracellular Vesicles from Human Osteoarthritic Subchondral Bone Multipotential Stromal Cells (MSCs)

Osteoarthritis (OA) is a heterogeneous disease in which the cross-talk between the cells from different tissues within the joint is affected as the disease progresses. Extracellular vesicles (EVs) are known to have a crucial role in cell-cell communication by means of cargo transfer. Subchondral bone (SB) resident cells and its microenvironment are increasingly recognised to have a major role in OA pathogenesis. The aim of this study was to investigate the EV production from OA SB mesenchymal stromal cells (MSCs) and their possible influence on OA chondrocytes. Small EVs were isolated from OA-MSCs, characterized and cocultured with chondrocytes for viability and gene expression analysis, and compared to small EVs from MSCs of healthy donors (H-EVs). OA-EVs enhanced viability of chondrocytes and the expression of chondrogenesis-related genes, although the effect was marginally lower compared to that of the H-EVs. miRNA profiling followed by unsupervised hierarchical clustering analysis revealed distinct microRNA sets in OA-EVs as compared to their parental MSCs or H-EVs. Pathway analysis of OA-EV miRNAs showed the enrichment of miRNAs implicated in chondrogenesis, stem cells, or other pathways related to cartilage and OA. In conclusion, OA SB MSCs were capable of producing EVs that could support chondrocyte viability and chondrogenic gene expression and contained microRNAs implicated in chondrogenesis support. These EVs could therefore mediate the cross-talk between the SB and cartilage in OA potentially modulating chondrocyte viability and endogenous cartilage regeneration.

Asporin regulated by miR-26b-5p mediates chondrocyte senescence and exacerbates osteoarthritis progression via TGF-β1/Smad2 pathway

OBJECTIVES

This study aimed to investigate the role and mechanism of asporin in modulating chondrocyte senescence in osteoarthritis (OA) pathology.

METHODS

Asporin and senescence-related hallmark expression were examined in human and experimental OA mouse cartilage samples. Twelve-week-old male C57 mice were administered with recombinant protein (rm-asporin)- or asporin-siRNA-expressing lentiviruses via intra-articular injection once a week after destabilization of the medial meniscus (DMM) surgery to induce OA. Cartilage damage was measured using the Osteoarthritis Research Society International score. Senescence-associated β-galactosidase (SA-βGal) staining, γH2AX, p21, and p16INK4a were analyzed by immunofluorescence staining and western blot to assess the specific role of asporin in chondrocyte senescence. The TGF-β1/Smad2 signaling pathway and miR-26b-5p were further evaluated to explore the mechanism of asporin in OA.

RESULTS

Asporin was upregulated in articular chondrocytes of OA patients and DMM mice and accompanied by accumulation of senescent cells. Asporin overexpression exaggerated OA progression, whereas silencing asporin restored chondrocyte homeostasis and deferred chondrocyte senescence, leading to markedly attenuated DMM-induced OA. Cellular and molecular analyses showed that asporin can be inhibited by miR-26b-5p, which was significantly downregulated in OA cartilage, leading to exacerbation of experimental OA partially through inhibition of TGF-β1/Smad2 signaling in chondrocytes.

CONCLUSIONS

Our findings indicate that asporin plays an essential role in chondrocyte senescence and OA pathogenesis. Upregulated by miR-26b-5p, asporin inhibits the TGF-β1/Smad2 pathway to accelerate chondrocyte senescence and exacerbate cartilage degeneration. Targeting the miR-26b-5p/asporin/Smad2 axis may serve as a practical therapeutic strategy to delay chondrocyte senescence and OA development.

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⁻⁵). 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.

A Network Biology Approach to Understanding the Tissue-Specific Roles of Non-Coding RNAs in Arthritis

Discovery of non-coding RNAs continues to provide new insights into some of the key molecular drivers of musculoskeletal diseases. Among these, microRNAs have received widespread attention for their roles in osteoarthritis and rheumatoid arthritis. With evidence to suggest that long non-coding RNAs and circular RNAs function as competing endogenous RNAs to sponge microRNAs, the net effect on gene expression in specific disease contexts can be elusive. Studies to date have focused on elucidating individual long non-coding-microRNA-gene target axes and circular RNA-microRNA-gene target axes, with a paucity of data integrating experimentally validated effects of non-coding RNAs. To address this gap, we curated recent studies reporting non-coding RNA axes in chondrocytes from human osteoarthritis and in fibroblast-like synoviocytes from human rheumatoid arthritis. Using an integrative computational biology approach, we then combined the findings into cell- and disease-specific networks for in-depth interpretation. We highlight some challenges to data integration, including non-existent naming conventions and out-of-date databases for non-coding RNAs, and some successes exemplified by the International Molecular Exchange Consortium for protein interactions. In this perspective article, we suggest that data integration is a useful in silico approach for creating non-coding RNA networks in arthritis and prioritizing interactions for further in vitro and in vivo experimentation in translational research.