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

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) \u3c 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

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 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 inter-donor variation in their in vitro chondrogenic differentiation potential. Underlying molecular mechanisms responsible for culture-expanded MSC heterogeneity remains poorly understood. In this study, we sought to identify variations in miRNA signatures associated with cultured equine MSC chondrogenic differentiation potential from different donors. Neocartilage tissue generated from equine cord blood-derived MSCs (CB-MSCs) was categorized as having either high or low chondrogenic potential based on their histological appearance and quantification of glycosaminoglycan deposition. Using next-generation sequencing, we identified 30 differentially expressed miRNAs amongst undifferentiated MSC cultures that corresponded with their chondrogenic potential. Of note, MSCs with low chondrogenic potential upregulated miR-146a and miR-487b-3p, which was also observed by qRT-PCR. Our findings suggest that miRNA profiling of equine MSC cultures may have prognostic value in selecting MSC donors with regards to their chondrogenic differentiation potential

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

Statistical data of differentially expressed miRNAs equine plasma of OA vs OCD with a 1.5-log2FC or greater change in expression (p < 0.05).

Statistical data of differentially expressed miRNAs equine plasma of OA vs OCD with a 1.5-log2FC or greater change in expression (p < 0.05).</p

ddPCR count data.

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.</div

Summary of miRNAs selected for validation and their expression profiles compared in NGS vs ddPCR in equine blood plasma.

Summary of miRNAs selected for validation and their expression profiles compared in NGS vs ddPCR in equine blood plasma.</p

RNA Seq count data.

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.</div