RNA sequencing data integration reveals an miRNA interactome of osteoarthritis cartilage

Widespread regulation of gene expression by glucocorticoids in chondrocytes from patients with osteoarthritis as determined by RNA-Seq

Background

Intra-articular glucocorticoid (GC) injections are widely used as a symptomatic treatment for osteoarthritis (OA). However, there are also concerns about their potentially harmful effects, and their detailed effects on chondrocyte phenotype remain poorly understood.

Methods

We studied the effects of dexamethasone on gene expression in OA chondrocytes with RNA-Seq. Chondrocytes were isolated from the cartilage from OA patients undergoing knee replacement surgery and cultured with or without dexamethasone for 24 h. Total RNA was isolated and sequenced, and functional analysis was performed against the Gene Ontology (GO) database. Results for selected genes were confirmed with RT-PCR. We also investigated genes linked to OA in recent genome-wide expression analysis (GWEA) studies.

Results

Dexamethasone increased the expression of 480 and reduced that of 755 genes with a fold change (FC) 2.0 or greater. Several genes associated with inflammation and cartilage anabolism/catabolism as well as lipid and carbohydrate metabolism were among the most strongly affected genes. In the GO analysis, genes involved in the extracellular matrix organization, cell proliferation and adhesion, inflammation, and collagen synthesis were enriched among the significantly affected genes. In network analysis, NGF, PI3KR1, and VCAM1 were identified as central genes among those most strongly affected by dexamethasone.

Conclusions

This is the first study investigating the genome-wide effects of GCs on the gene expression in OA chondrocytes. In addition to clear anti-inflammatory and anticatabolic effects, GCs affect lipid and glucose metabolism in chondrocytes, an observation that might be particularly important in the metabolic phenotype of OA.

Inhibition of miR-490-5p Promotes Human Adipose-Derived Stem Cells Chondrogenesis and Protects Chondrocytes via the PITPNM1/PI3K/AKT Axis

MicroRNAs (miRNAs) play a pivotal role in cartilage development and homeostasis in osteoarthritis (OA). While the fundamental roles of miRNAs in cartilage degeneration have been extensively studied, their effects on chondrogenic differentiation induced by human adipose-derived stem cells (hADSCs) and the underlying mechanisms remain largely elusive. Here, we investigated the roles and mechanisms of miRNAs in hADSC chondrogenic differentiation and chondrocyte homeostasis. Using microarray analysis, we screened miRNAs expressed in the chondrogenic differentiated hADSCs and identified miR-490-5p as the most significantly down-regulated miRNA. We analyzed its expression patterns during chondrogenesis in vivo and in vitro. Our study showed that miR-490-5p overexpression promoted the transition of hADSCs from chondrogenesis to osteogenesis. In addition, based on miRNA-mRNA prediction analysis and dual-luciferase reporter assay, we proposed and proved that miR-490-5p targeted PITPNM1 by binding to its 3'-UTR and inhibiting its translation. Moreover, loss- and gain-of-function experiments identified the involvement of the PI3K/AKT signaling pathway, and a rescue experiment determined the effect and specific mechanism of the miR-490-5p/PITPNM1/PI3K/AKT axis in hADSC chondrogenic differentiation and chondrocyte homeostasis. Inhibition of miR-490-5p alleviated cartilage injury in vivo as demonstrated using the destabilization of the medial meniscus (DMM) OA model. We identified miR-490-5p as a novel modulator of hADSC-mediated chondrogenesis and chondrocyte phenotype. This study highlighted that miR-490-5p attenuated hADSC chondrogenesis and accelerated cartilage degradation through activation of the PI3K/AKT signaling pathway by targeting PITPNM1. Inhibition of miR-490-5p facilitated hADSC chondrogenic differentiation and protected chondrocyte phenotype via the PITPNM1/PI3K/AKT axis, thus providing a novel stem cell potential therapeutic target for OA treatment.

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

OBJECTIVE

MicroRNAs act locally and systemically to impact osteoarthritis (OA) pathophysiology, but comprehensive profiling of the circulating miRNome in early vs late stages of OA has yet to be conducted. Sequencing has emerged as the preferred method for microRNA profiling since it offers high sensitivity and specificity. Our objective was to sequence the miRNome in plasma from 91 patients with early [Kellgren-Lawrence (KL) grade 0 or 1 (n = 41)] or late [KL grade 3 or 4 (n = 50)] symptomatic radiographic knee OA to identify unique microRNA signatures in each disease state.

DESIGN

MicroRNA libraries were prepared using the QIAseq miRNA Library Kit and sequenced on the Illumina NextSeq 550. Counts were produced for microRNAs captured in miRBase and for novel microRNAs. Statistical, bioinformatics, and computational biology approaches were used to refine and interpret the final list of microRNAs.

RESULTS

From 215 differentially expressed microRNAs (FDR < 0.01), 97 microRNAs showed an increase or decrease in expression in ≥85% of samples in the early OA group as compared to the median expression in the late OA group. Increasing this threshold to ≥95%, seven microRNAs were identified: hsa-miR-335-3p, hsa-miR-199a-5p, hsa-miR-671-3p, hsa-miR-1260b, hsa-miR-191-3p, hsa-miR-335-5p, and hsa-miR-543. Four novel microRNAs were present in ≥50% of early OA samples and had 27 predicted gene targets in common with the prioritized set of predicted gene targets from the 97 microRNAs, suggesting common underlying mechanisms.

CONCLUSION

Sequencing of well-characterized patient cohorts produced unbiased profiling of the circulating miRNome and identified a unique panel of 11 microRNAs in early radiographic knee OA.

OATargets: a knowledge base of genes associated with osteoarthritis joint damage in animals

Objectives

To collate the genes experimentally modulated in animal models of osteoarthritis (OA) and compare these data with OA transcriptomics data to identify potential therapeutic targets.

Methods

PubMed searches were conducted to identify publications describing gene modulations in animal models. Analysed gene expression data were retrieved from the SkeletalVis database of analysed skeletal microarray and RNA-Seq expression data. A network diffusion approach was used to predict new genes associated with OA joint damage.

Results

A total of 459 genes were identified as having been modulated in animal models of OA, with ageing and post-traumatic (surgical) models the most prominent. Ninety-eight of the 143 genes (69%) genetically modulated more than once had a consistent effect on OA joint damage severity. Several discrepancies between different studies were identified, providing lessons on interpretation of these data. We used the data collected along with OA gene expression data to expand existing annotations and prioritise the most promising therapeutic targets, which we validated using the latest reported associations. We constructed an online database OATargets to allow researchers to explore the collated data and integrate it with existing OA and skeletal gene expression data.

Conclusions

We present a comprehensive survey and online resource for understanding gene regulation of animal model OA pathogenesis.

Elucidating Epigenetic Regulation by Identifying Functional cis-Acting Long Noncoding RNAs and Their Targets in Osteoarthritic Articular Cartilage

Objective

To identify robustly differentially expressed long noncoding RNAs (lncRNAs) with osteoarthritis (OA) pathophysiology in cartilage and to explore potential target messenger RNA (mRNA) by establishing coexpression networks, followed by functional validation.

Methods

RNA sequencing was performed on macroscopically lesioned and preserved OA cartilage from patients who underwent joint replacement surgery due to OA (n = 98). Differential expression analysis was performed on lncRNAs that were annotated in GENCODE and Ensembl databases. To identify potential interactions, correlations were calculated between the identified differentially expressed lncRNAs and the previously reported differentially expressed protein‐coding genes in the same samples. Modulation of chondrocyte lncRNA expression was achieved using locked nucleic acid GapmeRs.

Results

By applying our in‐house pipeline, we identified 5,053 lncRNAs that were robustly expressed, of which 191 were significantly differentially expressed (according to false discovery rate) between lesioned and preserved OA cartilage. Upon integrating mRNA sequencing data, we showed that intergenic and antisense differentially expressed lncRNAs demonstrate high, positive correlations with their respective flanking sense genes. To functionally validate this observation, we selected P3H2‐AS1, which was down‐regulated in primary chondrocytes, resulting in the down‐regulation of P3H2 gene expression levels. As such, we can confirm that P3H2‐AS1 regulates its sense gene P3H2.

Conclusion

By applying an improved detection strategy, robustly differentially expressed lncRNAs in OA cartilage were detected. Integration of these lncRNAs with differential mRNA expression levels in the same samples provided insight into their regulatory networks. Our data indicate that intergenic and antisense lncRNAs play an important role in regulating the pathophysiology of OA.

Comprehensive analysis of differentially expressed microRNAs and mRNAs involved in diabetic corneal neuropathy

AIMS

Corneal nerve fibers are derived from the ophthalmic division of the trigeminal ganglion (TG). Here, by sequencing of microRNAs (miRNAs) and messenger RNAs (mRNAs) from diabetic and normal TG tissues, we aimed to uncover potential miRNAs, mRNAs, and the network of their interactions involved in the pathogenesis of diabetic corneal neuropathy.

MAIN METHODS

We performed RNA sequencing to systematically screen out differentially expressed miRNAs and mRNAs in TG tissues from diabetic and normal mice. Functional enrichment analyses were performed to illustrate the biological functions of differentially expressed mRNAs (DEmRNAs). Following this, miRNA-mRNA regulatory networks were built by means of bioinformatics methods to suggest regulatory role for miRNAs in the pathogenesis of diabetic corneal neuropathy. Finally, the credibility of the sequencing-based results was validated using qRT-PCR.

KEY FINDINGS

Sequencing analyses disclosed that 68 miRNAs and 114 mRNAs were differentially expressed in diabetic TG tissues compared with normal TG samples. The functional analyses showed that DEmRNAs participated in diabetes-related biological processes. After applying an optimized approach to predict miRNA-mRNA pairs, a miRNA-mRNA interacting network was inferred. Subsequently, the expression and correlation of miR-350-5p and Mup20, miR-592-5p and Angptl7 as well as miR-351-5p and Elovl6 were preliminarily validated.

SIGNIFICANCE

Our study provides a systematic characterization of miRNA and mRNA expression in the TG during diabetic corneal neuropathy and will contribute to the development of clinical diagnostic and therapeutic strategies for diabetic corneal neuropathy.

The miRNA-mRNA interactome of murine induced pluripotent stem cell-derived chondrocytes in response to inflammatory cytokines

Osteoarthritis (OA) is a degenerative joint disease, and inflammation within an arthritic joint plays a critical role in disease progression. Pro-inflammatory cytokines, specifically IL-1 and TNF-α, induce aberrant expression of catabolic and degradative enzymes and inflammatory cytokines in OA and result in a challenging environment for cartilage repair and regeneration. MicroRNAs (miRNAS) are small noncoding RNAs and are important regulatory molecules that act by binding to target messenger RNAs (mRNAs) to reduce protein synthesis and have been implicated in many diseases, including OA. The goal of this study was to understand the mechanisms of miRNA regulation of the transcriptome of tissue-engineered cartilage in response to IL-1β and TNF-α using an in vitro murine induced pluripotent stem cell (miPSC) model system. We performed miRNA and mRNA sequencing to determine the temporal and dynamic responses of genes to specific inflammatory cytokines as well as miRNAs that are differentially expressed (DE) in response to both cytokines or exclusively to IL-1β or TNF-α. Through integration of mRNA and miRNA sequencing data, we created networks of miRNA-mRNA interactions which may be controlling the response to inflammatory cytokines. Within the networks, hub miRNAs, miR-29b-3p, miR-17-5p, and miR-20a-5p, were identified. As validation of these findings, we found that delivery of miR-17-5p and miR-20a-5p mimics significantly decreased degradative enzyme activity levels while also decreasing expression of inflammation-related genes in cytokine-treated cells. This study utilized an integrative approach to determine the miRNA interactome controlling the response to inflammatory cytokines and novel mediators of inflammation-driven degradation in tissue-engineered cartilage.

Genetics in Cartilage Lesions: Basic Science and Therapy Approaches

Cartilage lesions have a multifactorial nature, and genetic factors are their strongest determinants. As biochemical and genetic studies have dramatically progressed over the past decade, the molecular basis of cartilage pathologies has become clearer. Several homeostasis abnormalities within cartilaginous tissue have been found, including various structural changes, differential gene expression patterns, as well as altered epigenetic regulation. However, the efficient treatment of cartilage pathologies represents a substantial challenge. Understanding the complex genetic background pertaining to cartilage pathologies is useful primarily in the context of seeking new pathways leading to disease progression as well as in developing new targeted therapies. A technology utilizing gene transfer to deliver therapeutic genes to the site of injury is quickly becoming an emerging approach in cartilage renewal. The goal of this work is to provide an overview of the genetic basis of chondral lesions and the different approaches of the most recent systems exploiting therapeutic gene transfer in cartilage repair. The integration of tissue engineering with viral gene vectors is a novel and active area of research. However, despite promising preclinical data, this therapeutic concept needs to be supported by the growing body of clinical trials.

Knee Osteoarthritis: A Review of Pathogenesis and State-Of-The-Art Non-Operative Therapeutic Considerations

Being the most common musculoskeletal progressive condition, osteoarthritis is an interesting target for research. It is estimated that the prevalence of knee osteoarthritis (OA) among adults 60 years of age or older is approximately 10% in men and 13% in women, making knee OA one of the leading causes of disability in elderly population. Today, we know that osteoarthritis is not a disease characterized by loss of cartilage due to mechanical loading only, but a condition that affects all of the tissues in the joint, causing detectable changes in tissue architecture, its metabolism and function. All of these changes are mediated by a complex and not yet fully researched interplay of proinflammatory and anti-inflammatory cytokines, chemokines, growth factors and adipokines, all of which can be measured in the serum, synovium and histological samples, potentially serving as biomarkers of disease stage and progression. Another key aspect of disease progression is the epigenome that regulates all the genetic expression through DNA methylation, histone modifications, and mRNA interference. A lot of work has been put into developing non-surgical treatment options to slow down the natural course of osteoarthritis to postpone, or maybe even replace extensive surgeries such as total knee arthroplasty. At the moment, biological treatments such as platelet-rich plasma, bone marrow mesenchymal stem cells and autologous microfragmented adipose tissue containing stromal vascular fraction are ordinarily used. Furthermore, the latter two mentioned cell-based treatment options seem to be the only methods so far that increase the quality of cartilage in osteoarthritis patients. Yet, in the future, gene therapy could potentially become an option for orthopedic patients. In the following review, we summarized all of the latest and most important research in basic sciences, pathogenesis, and non-operative treatment.

Knee Osteoarthritis: A Review of Pathogenesis and State-Of-The-Art Non-Operative Therapeutic Considerations

Being the most common musculoskeletal progressive condition, osteoarthritis is an interesting target for research. It is estimated that the prevalence of knee osteoarthritis (OA) among adults 60 years of age or older is approximately 10% in men and 13% in women, making knee OA one of the leading causes of disability in elderly population. Today, we know that osteoarthritis is not a disease characterized by loss of cartilage due to mechanical loading only, but a condition that affects all of the tissues in the joint, causing detectable changes in tissue architecture, its metabolism and function. All of these changes are mediated by a complex and not yet fully researched interplay of proinflammatory and anti-inflammatory cytokines, chemokines, growth factors and adipokines, all of which can be measured in the serum, synovium and histological samples, potentially serving as biomarkers of disease stage and progression. Another key aspect of disease progression is the epigenome that regulates all the genetic expression through DNA methylation, histone modifications, and mRNA interference. A lot of work has been put into developing non-surgical treatment options to slow down the natural course of osteoarthritis to postpone, or maybe even replace extensive surgeries such as total knee arthroplasty. At the moment, biological treatments such as platelet-rich plasma, bone marrow mesenchymal stem cells and autologous microfragmented adipose tissue containing stromal vascular fraction are ordinarily used. Furthermore, the latter two mentioned cell-based treatment options seem to be the only methods so far that increase the quality of cartilage in osteoarthritis patients. Yet, in the future, gene therapy could potentially become an option for orthopedic patients. In the following review, we summarized all of the latest and most important research in basic sciences, pathogenesis, and non-operative treatment.

LDHA-mediated ROS generation in chondrocytes is a potential therapeutic target for osteoarthritis

The contribution of inflammation to the chronic joint disease osteoarthritis (OA) is unclear, and this lack of clarity is detrimental to efforts to identify therapeutic targets. Here we show that chondrocytes under inflammatory conditions undergo a metabolic shift that is regulated by NF-κB activation, leading to reprogramming of cell metabolism towards glycolysis and lactate dehydrogenase A (LDHA). Inflammation and metabolism can reciprocally modulate each other to regulate cartilage degradation. LDHA binds to NADH and promotes reactive oxygen species (ROS) to induce catabolic changes through stabilization of IκB-ζ, a critical pro-inflammatory mediator in chondrocytes. IκB-ζ is regulated bi-modally at the stages of transcription and protein degradation. Overall, this work highlights the function of NF-κB activity in the OA joint as well as a ROS promoting function for LDHA and identifies LDHA as a potential therapeutic target for OA treatment.

Exosomes from patients with major depression cause depressive-like behaviors in mice with involvement of miR-139-5p-regulated neurogenesis

Exosomal microRNAs (miRNAs) have been suggested to participate in the pathogenesis of neuropsychiatric diseases, but their role in major depressive disorder (MDD) is unknown. We performed a genome-wide miRNA expression profiling of blood-derived exosomes from MDD patients and control subjects and revealed the top differentially expressed exosomal miRNA, i.e. hsa-miR-139-5p (upregulation), had good performance to differentiate between MDD patients and controls. Tail vein injection of blood exosomes isolated from MDD patients into normal mice caused their depressive-like behaviors as determined by the forced swimming, tail suspension, and novelty suppressed feeding tests, and injection of blood exosomes isolated from healthy volunteers into unpredictable mild stress (CUMS)-treated mice alleviated their depressive-like behaviors. CUMS mice also showed significantly increased blood and brain levels of exosomal miR-139-5p. Furthermore, the depressive-like behaviors in CUMS-treated mice were rescued by intranasal injection of miR-139-5p antagomir, suggesting that increased exosomal miR-139-5p levels may mediate stress-induced depression-like behavior in mice. Both exosome treatment and miR-139-5p antagomir treatment increased hippocampal neurogenesis in the CUMS-treated mice, and treatment of exosome from MDD patients decreased hippocampal neurogenesis in the normal mice. The role of miR-139-5p in neurogenesis was validated by in vitro experiments, demonstrating that miR-139-5p is a negative regulator for neural stem cell proliferation and neuronal differentiation. Our findings together suggest that exosomes from patients with major depression caused depressive-like behaviors in mice with involvement of miR-139-5p-regulated neurogenesis. Therefore, exosomal miRNAs are promising targets for the diagnosis and treatment of MDD.

Overexpression of miR-10a-5p facilitates the progression of osteoarthritis

The current study was aimed at exploring the potential roles and possible mechanisms of miR-10a-5p in osteoarthritis (OA). We performed RT-qPCR, Western blot, CCK8, EdU Assay, and flow cytometry assay to clarify the roles of miR-10a-5p in OA. Furthermore, the whole transcriptome sequencing together with integrated bioinformatics analyses were conducted to elucidate the underlying mechanisms of miR-10a-5p involving in OA. Our results demonstrated that miR-10a-5p was upregulated in OA and acted as a significant contributing factor for OA. A large number of circRNAs, lncRNAs, miRNAs, and mRNAs were identified by overexpressing miR-10a-5p. Functional enrichment analyses indicated that these differentially-expressed genes were enriched in some important terms including PPAR signaling pathway, PI3K-Akt signaling pathway, and p53 signaling pathway. A total of 42 hub genes were identified in the protein-protein interaction network including SERPINA1, TTR, APOA1, and A2M. Also, we constructed the network regulatory interactions across coding and noncoding RNAs triggered by miR-10a-5p, which revealed the powerful regulating effects of miR-10a-5p. Moreover, we found that HOXA3 acted as the targeted genes of miR-10a-5p and miR-10a-5p contributed to the progression of OA by suppressing HOXA3 expression. Our findings shed insight on regulatory mechanisms of miR-10a-5p, which might provide novel therapeutic targets for OA.

Osteoarthritis year in review 2019: genetics, genomics and epigenetics

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