Cell and transcriptomic diversity of infrapatellar fat pad during knee osteoarthritis.

OBJECTIVES: In this study, we employ a multiomic approach to identify major cell types and subsets, and their transcriptomic profiles within the infrapatellar fat pad (IFP), and to determine differences in the IFP based on knee osteoarthritis (KOA), sex and obesity status. METHODS: Single-nucleus RNA sequencing of 82 924 nuclei from 21 IFPs (n=6 healthy control and n=15 KOA donors), spatial transcriptomics and bioinformatic analyses were used to identify contributions of the IFP to KOA. We mapped cell subclusters from other white adipose tissues using publicly available literature. The diversity of fibroblasts within the IFP was investigated by bioinformatic analyses, comparing by KOA, sex and obesity status. Metabolomics was used to further explore differences in fibroblasts by obesity status. RESULTS: We identified multiple subclusters of fibroblasts, macrophages, adipocytes and endothelial cells with unique transcriptomic profiles. Using spatial transcriptomics, we resolved distributions of cell types and their transcriptomic profiles and computationally identified putative cell-cell communication networks. Furthermore, we identified transcriptomic differences in fibroblasts from KOA versus healthy control donor IFPs, female versus male KOA-IFPs and obese versus normal body mass index (BMI) KOA-IFPs. Finally, using metabolomics, we defined differences in metabolite levels in supernatants of naïve, profibrotic stimuli-treated and proinflammatory stimuli-treated fibroblasts from obese compared to normal BMI KOA-IFPs. CONCLUSIONS: Overall, by employing a multiomic approach, this study provides the first comprehensive map of the cellular and transcriptomic diversity of human IFP and identifies IFP fibroblasts as key cells contributing to transcriptomic and metabolic differences related to KOA disease, sex or obesity.

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.

Transcriptomics and metabolomics: Challenges of studying obesity in osteoarthritis.

Objective: Obesity is a leading risk factor for both the incidence and progression of osteoarthritis (OA). Omic technologies, including transcriptomics and metabolomics are capable of identifying RNA and metabolite profiles in tissues and biofluids of OA patients. The objective of this review is to highlight studies using transcriptomics and metabolomics that contribute to our understanding of OA pathology in relation to obesity. Design: We conducted a targeted search of PUBMED for articles, and GEO for datasets, published up to February 13, 2024, screening for those using high-throughput transcriptomic and metabolomic techniques to study human or pre-clinical animal model tissues or biofluids related to obesity-associated OA. We describe relevant studies and discuss challenges studying obesity as a disease-related factor in OA. Results: Of the 107 publications identified by our search criteria, only 15 specifically used transcriptomics or metabolomics to study joint tissues or biofluids in obesity-related OA. Specific transcriptomic and metabolomic signatures associated with obesity-related OA have been defined in select local joint tissues, biofluids and other biological material. However, considerable challenges exist in understanding contributions of obesity-associated modifications of transcriptomes and metabolomes related to OA, including sociodemographic, anthropometric, dietary and molecular redundancy-related factors. Conclusions: A number of additional transcriptomic and metabolomic studies are needed to comprehensively understand how obesity affects OA incidence, progression and outcomes. Integration of transcriptome and metabolome signatures from multiple tissues and biofluids, using network-based approaches will likely help to better define putative therapeutic targets that could enable precision medicine approaches to obese OA patients.

microRNAs are differentially expressed in equine plasma of horses with osteoarthritis and osteochondritis dissecans versus control horses.

Osteoarthritis (OA) is a leading cause of lameness in horses with no effective disease-modifying treatment and challenging early diagnosis. OA is considered a disease of the joint involving the articular cartilage, subchondral bone, synovial membrane, and ligaments. Osteochondritis dissecans (OCD) is a joint disease consisting of focal defects in the osteochondral unit which may progress to OA later in life. MicroRNAs (miRNAs) have been recognized as small non-coding RNAs that regulate a variety of biological processes and have been detected in biological fluids. MiRNAs are currently investigated for their utility as biomarkers and druggable targets for a variety of diseases. The current study hypothesizes that miRNA profiles can be used to actively monitor joint health and differences in miRNA profiles will be found in healthy vs diseased joints and that differences will be detectable in blood plasma of tested horses. Five horses with OA, OCD, and 4 controls (C) had blood plasma and synovial fluid collected. Total RNA, including miRNA was isolated before generating miRNA libraries from the plasma of the horses. Libraries were sequenced at the Schroeder Arthritis Institute (Toronto). Differential expression analysis was done using DESeq2 and validated using ddPCR. KEGG pathway analysis was done using mirPath v.3 (Diana Tools). 57 differentially expressed miRNAs were identified in OA vs C plasma, 45 differentially expressed miRNAs in OC vs C plasma, and 21 differentially expressed miRNAs in OA vs OCD plasma. Notably, miR-140-5p expression was observed to be elevated in OA synovial fluid suggesting that miR-140-5p may serve as a protective marker early on to attenuate OA progression. KEGG pathway analysis of differentially expressed plasma miRNAs showed relationships with glycan degradation, glycosaminoglycan degradation, and hippo signaling pathway. Interestingly, ddPCR was unable to validate the NGS data suggesting that isomiRs may play an integral role in miRNA expression when assessed using NGS technologies.

MicroRNAs as Prognostic Markers for Chondrogenic Differentiation Potential of Equine Mesenchymal Stromal Cells.

Mesenchymal stromal cells (MSCs) are a promising cell source for cartilage tissue regeneration in animals and humans but with large interdonor variation in their in vitro chondrogenic differentiation potential. Underlying molecular mechanisms responsible for culture-expanded MSC heterogeneity remain poorly understood. In this study, we sought to identify variations in microRNA (miRNA) signatures associated with cultured equine MSC chondrogenic differentiation potential from different donors. Neocartilage tissue generated from equine cord blood-derived MSCs was categorized as having either high or low chondrogenic potential (LCP) based on their histological appearance and quantification of glycosaminoglycan deposition. Using next-generation sequencing, we identified 30 differentially expressed miRNAs among undifferentiated MSC cultures that corresponded with their chondrogenic potential. Of note, MSCs with LCP upregulated miR-146a and miR-487b-3p, which was also observed by quantitative real-time polymerase chain reaction. Our findings suggest that miRNA profiling of equine MSC cultures may have prognostic value in selecting MSC donors with regard to their chondrogenic differentiation potential.

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.

Contribution of MicroRNA-27b-3p to Synovial Fibrotic Responses in Knee Osteoarthritis.

OBJECTIVE: Synovial fibrosis contributes to osteoarthritis (OA) pathology, but the underlying mechanisms remain unknown. We have observed increased microRNA-27b-3p (miR-27b-3p) levels in synovial fluid of patients with late-stage radiographic knee OA. Here, we investigated the contribution of miR-27b-3p to synovial fibrosis in patients with severe knee OA and in a mouse model of knee OA. METHODS: We stained synovium sections obtained from patients with radiographic knee OA scored according to the Kellgren/Lawrence scale and mice that underwent destabilization of the medial meniscus (DMM) for miR-27b-3p using in situ hybridization. We examined the effects of intraarticular injection of miR-27b-3p mimic into naive mouse knee joints and intraarticular injection of a miR-27b-3p inhibitor into mouse knee joints after DMM. We performed transfection with miR-27b-3p mimic and miR-27b-3p inhibitor in human OA fibroblast-like synoviocytes (FLS) using reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) array, RNA sequencing, RT-qPCR, Western blotting, immunofluorescence, and migration assays. RESULTS: We observed increased miR-27b-3p expression in the synovium from patients with knee OA and in mice with DMM-induced arthritis. Injection of the miR-27b-3p mimic in mouse knee joints induced a synovial fibrosis-like phenotype, increased synovitis scores, and increased COL1A1 and α-smooth muscle actin (α-SMA) expression. In the mouse model of DMM-induced arthritis, injection of the miR-27b-3p inhibitor decreased α-SMA but did not change COL1A1 expression levels or synovitis scores. Transfection with the miR-27b-3p mimic in human OA FLS induced profibrotic responses, including increased migration and expression of key extracellular matrix (ECM) genes, but transfection with the miR-27b-3p inhibitor had the opposite effects. RNA sequencing identified a PPARG/ADAMTS8 signaling axis regulated by miR-27b-3p in OA FLS. Human OA FLS transfected with miR-27b-3p mimic and then treated with the PPARG agonist rosiglitazone or with ADAMTS8 small interfering RNA exhibited altered expression of select ECM genes. CONCLUSION: Our findings demonstrate that miR-27b-3p has a key role in ECM regulation associated with synovial fibrosis during 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.

A bioinformatics approach to microRNA-sequencing analysis.

The rapid expansion of Next Generation Sequencing (NGS) data availability has made exploration of appropriate bioinformatics analysis pipelines a timely issue. Since there are multiple tools and combinations thereof to analyze any dataset, there can be uncertainty in how to best perform an analysis in a robust and reproducible manner. This is especially true for newer omics applications, such as miRNomics, or microRNA-sequencing (miRNA-sequencing). As compared to transcriptomics, there have been far fewer miRNA-sequencing studies performed to date, and those that are reported seldom provide detailed description of the bioinformatics analysis, including aspects such as Unique Molecular Identifiers (UMIs). In this article, we attempt to fill the gap and help researchers understand their miRNA-sequencing data and its analysis. This article will specifically discuss a customizable miRNA bioinformatics pipeline that was developed using miRNA-sequencing datasets generated from human osteoarthritis plasma samples. We describe quality assessment of raw sequencing data files, reference-based alignment, counts generation for miRNA expression levels, and novel miRNA discovery. This report is expected to improve clarity and reproducibility of the bioinformatics portion of miRNA-sequencing analysis, applicable across any sample type, to promote sharing of detailed protocols in the NGS field.

MicroRNA-34a-5p Promotes Joint Destruction During Osteoarthritis.

OBJECTIVE: MicroRNA-34a-5p (miR-34a-5p) expression is elevated in the synovial fluid of patients with late-stage knee osteoarthritis (OA); however, its exact role and therapeutic potential in OA remain to be fully elucidated. This study was undertaken to examine the role of miR-34a-5p in OA pathogenesis. METHODS: Expression of miR-34a-5p was determined in joint tissues and human plasma (n = 71). Experiments using miR-34a-5p mimic or antisense oligonucleotide (ASO) treatment were performed in human OA chondrocytes, fibroblast-like synoviocytes (FLS) (n = 7-9), and mouse OA models, including destabilization of the medial meniscus (DMM; n = 22) and the accelerated, more severe model of mice fed a high-fat diet and subjected to DMM (n = 11). Wild-type (WT) mice (n = 9) and miR-34a-knockout (KO) mice (n = 11) were subjected to DMM. Results were expressed as the mean ± SEM and analyzed by t-test or analysis of variance, with appropriate post hoc tests. P values less than 0.05 were considered significant. RNA sequencing was performed on WT and KO mouse chondrocytes. RESULTS: Expression of miR-34a-5p was significantly increased in the plasma, cartilage, and synovium of patients with late-stage OA and in the cartilage and synovium of mice subjected to DMM. Plasma miR-34a-5p expression was significantly increased in obese patients with late-stage OA, and in the plasma and knee joints of mice fed a high-fat diet. In human OA chondrocytes and FLS, miR-34a-5p mimic increased key OA pathology markers, while miR-34a-5p ASO improved cellular gene expression. Intraarticular miR-34a-5p mimic injection induced an OA-like phenotype. Conversely, miR-34a-5p ASO injection imparted cartilage-protective effects in the DMM and high-fat diet/DMM models. The miR-34a-KO mice exhibited protection against DMM-induced cartilage damage. RNA sequencing of WT and KO chondrocytes revealed a putative miR-34a-5p signaling network. CONCLUSION: Our findings provide comprehensive evidence of the role and therapeutic potential of miR-34a-5p in OA.

Effect of autotaxin inhibition in a surgically-induced mouse model of osteoarthritis.

Objective: Osteoarthritis (OA) is a degenerative joint disease with no approved disease modifying therapy. The enzyme autotaxin (ATX) converts lysophoshatidylcholine analogues to lysophosphatidic acid. Systemic inhibition of ATX reduces pain in animal models of OA; however, OA disease-modifying effects associated with ATX inhibition remain unknown. Here, we sought to determine whether local (knee joint) injection of an ATX inhibitor attenuates surgically-induced OA in mice. Methods: ATX expression was evaluated in human knee OA cartilage. Ten-week-old mice were subjected to surgically-induced OA. ATX inhibitor (PF-8380, 2.5ng/joint) was injected intra-articularly either at three and five weeks post-surgery or at two, four, six and eight weeks post-surgery and knee joints were evaluated by histopathology and immunohistochemistry to study the expression of catabolic and cell death markers. mRNA sequencing of human OA chondrocytes treated with/without ATX inhibitor was performed to identify differentially expressed transcripts, followed by pathway enrichment analysis. Results: ATX expression was elevated in severely degenerated cartilage compared to less degenerated human OA cartilage. In surgically-induced OA mice, intra-articular injection of ATX inhibitor at three and five weeks post-surgery partially protected knee joints from cartilage degeneration. Interestingly, earlier and more frequent ATX inhibitor injections did not confer significant protection. Immunohistochemical analysis showed decreased expression of catabolic and apoptotic markers with two ATX inhibitor injections. mRNA sequencing followed by pathway analysis identified pathways related to cholesterol analogue metabolism and cell-cycle that could be modulated by ATX inhibition in human OA chondrocytes. Conclusion: Local delivery of ATX inhibitor partially attenuates surgically-induced OA in mice.