Skip to main content
  • Research article
  • Open access
  • Published:

Association of circulating microRNAs with prevalent and incident knee osteoarthritis in women: the OFELY study



In the context of the scarcity of biomarkers for knee osteoarthritis (OA), we examined the associations of prevalent and incident OA with the expression levels of serum miRNAs in subjects with and without OA.


With a next-generation sequencing approach, we compared the miRome expression of 10 women with knee OA and 10 age-matched healthy subjects. By real-time qPCR, we analyzed the expression levels of 19 miRNAs at baseline selecting 43 women with prevalent knee OA (Kellgren Lawrence score of 2/3), 23 women with incident knee OA over a 4-year follow-up and 67 healthy subjects without prevalent or incident OA matched for age and body mass index.


Serum miR-146a-5p was significantly increased in the group of prevalent knee OA compared with controls (relative quantification (RQ); median [Interquartile range] 1.12 [0.73; 1.46] vs 0.85 [0.62; 1.03], p = 0.015). The likelihood of prevalent knee OA was significantly increased (odds ratio [95% confidence interval (CI)] 1.83 [1.21–2.77], p = 0.004) for each quartile increase in serum miR-146a-5p. The women with miR-146a-5p levels above the median (0.851) had a higher risk of prevalent knee OA compared to those below the median [95% CI] 4.62 [1.85–11.5], p = 0.001. Moreover, we found a significant association between the baseline level of serum miR-186-5p and the risk of incident knee OA (Q4 vs Q1–3; odds ratio [95% CI] 6.13 [1.14–32.9], p = 0.034).


We showed for the first time that miR-146a-5p and miR-186-5p are significantly associated with prevalent and incident knee OA, respectively.

Key messages

  1. 1.

    Circulating miR-146a-5p is significantly associated with prevalent knee OA in women.

  2. 2.

    Circulating miR-186-5p is significantly associated with incident knee OA in women.

  3. 3.

    Circulating miR-186-5p is a potential non-invasive biomarker of early-stage OA in women.


Osteoarthritis (OA) is the most frequent chronic musculoskeletal disease affecting approximately 40% of adults aged 70 years and over [1]. It is considered as a slowly progressive disease degrading all tissues of the affected joint, appearing as a destruction of the cartilage, the hallmark of OA, but also a mild-to-moderate synovial inflammation and an alteration of the subchondral bone structure. So far, the gold standard to estimate the extent of the disease has been plain radiography, but the poor sensitivity and the relatively large precision error lead to disease detection only when significant cartilage degradation has already occurred. Measuring circulating biomarkers appears as a potential non-invasive approach to diagnose the disease and prognosticate its evolution. Research has focused on the structural components of the cartilage matrix, especially type II collagen degradation markers [2]. Because of the large overlap between biomarker concentrations in controls and OA patients, the diagnosis of OA in an individual patient cannot be made by currently available biomarkers. Given the limitations of the tools that are currently available for OA assessment, there is considerable interest in the identification of specific biological markers that reflect quantitative and dynamic variations in joint tissue remodeling. For this purpose, microRNAs (miRNAs) are new potential targets. These small non-coding RNAs of 22–28 nucleotides in length can silence gene expression by binding to target messenger RNA (mRNA) repressing the translation. This regulation process is redundant because 1 miRNA can inhibit a large number of mRNA, and 1 gene can be targeted by multiple miRNAs. The miRNAs are important players in maintaining the health of all joint tissues including the cartilage, ligament, tendon, and muscle. Some of them have a potential role in regulating homeostasis of joint tissues in the context of OA [3]. Several in vitro and in vivo studies have reported the miRNA involvement in the OA onset and progression, by targeting cartilage-associated genes [4]. They regulate the expression of genes involved in pathways altered in OA chondrocytes such as apoptosis [5], expression levels of MMPs and ADAMTS [6,7,8], and chondrocyte signaling [9, 10]. To date, it has been reported that about 80 miRNAs are involved in the pathology of OA and that their circulating level in patients might reflect the underlying disease state [4]. Indeed, circulating miRNAs are considered to be attractive molecules to serve as biomarkers because they are easily accessible, stable in biofluids, and quantifiable with a high degree of sensitivity and specificity. Circulating miRNAs are protected from degradation because they circulate bound to proteins such as Argonaute 2 and nucleoplasmin or inside microvesicles such as exosomes [11,12,13]. Finally, although the expression of miRNAs has been well documented in OA cartilage, it remains difficult to understand to which extent circulating miRNAs reflect the tissue alteration. Moreover, the aberrant expression levels previously reported in circulating miRNAs differ considerably between clinical studies due to the different characteristics of these settings [14,15,16,17,18,19,20,21,22].

In this study, a first screening phase was conducted by comparing the serum miRome in subjects with (n = 10) and without (n = 10) knee OA in the OFELY cohort by next-generation sequencing (NGS) analysis. Then, on the basis of our NGS data and on the previously reported function of miRNAs in the OA process, both in basic and/or clinical research, we have selected 19 candidate miRNAs for quantitative real-time polymerase chain reaction (RT-qPCR) analysis. By conducting the validation phase of these candidate miRNAs in the OFELY cohort, we report that miR-146a-5p and miR-186-5p are significantly associated with prevalent and incident knee OA, respectively.


Patients selected from the OFELY cohort

The study group for both steps, NGS and RT-qPCR, included French women belonging to a population-based cohort. These women were participants in a prospective investigation of the determinants of bone loss, the Os des FEmmes de LYon (OFELY) study. This cohort has previously been described in detail elsewhere [23, 24]. In the present analysis, we studied the miRNA expression in women with knee OA at baseline for this analysis, 8 years after the recruitment of the cohort (ninth follow-up visit). The primary or post-traumatic etiology of knee OA was not known, but the selected women had no other skeletal diseases. The evaluation of OA disease was performed by radiography for knee OA and spine disc degeneration, by clinical examination for hand OA, and by a questionnaire for the hip OA (for details on the OA evaluation, see [8]). After baseline OA assessment, women were followed prospectively for 4 years.

Assessment of knee OA

Radiographs of both knees were obtained in all women. Radiologic evaluation of the knees consisted of bilateral posteroanterior weight-bearing knee radiographs with fixed flexion using the SynaFlex X-ray Positioning Frame (Synarc, San Francisco, CA), as previously described [25]. Radiographs were obtained in a single radiography unit by the same staff of two technicians using a previously described standardized protocol [26]. The severity of OA was performed and graded according to Kellgren and Lawrence (KL) classification at baseline and 4 years later [27]. Prevalent knee OA was defined by a KL score higher or equal to 2 at baseline. Incident OA was defined by a KL score higher or equal to 2 at year 4 and a KL score < 2 at baseline. All knee radiographs were scored by a single trained rheumatologist (ES-R). Measurements were made paired but not blinded to sequence, which has been shown not to modify the sensitivity to change [28].

The miRNA analysis by next-generation sequencing

Selection of patients

The expression levels of miRNAs in the serum were measured in ten women with knee OA (KL score of 2 and 3) and in ten healthy women without OA at any site (for details, see [24]). Both groups were matched for age (healthy 61.9 ± 3.03 years and OA 63.9 ± 3.4 years, p = 0.17). The flow chart for the design of the experimental study is reported in Fig. 1.

Fig. 1
figure 1

The flow chart of the study

NGS analysis

Total RNA was extracted from 400 μl of serum with the miRCURY Biofluids extraction kit (Exiqon®, Denmark) and analyzed by small RNA sequencing on a NextSeq500 sequencing instrument (Illumina platform). After RNA conversion into miRNA NGS libraries using the NEBNEXT library generation kit, cDNA was pre-amplified prior to library purification and quantification (for detailed procedures of extraction, library preparation and purification, normalization and quality controls, see Additional file 1: Data S1). Measurements were expressed as tags per million (TPM) in which the number of reads for a particular miRNA is divided by the total number of mapped reads in a sample and multiplied by 106. The differential expression of miRNA between the OA and control groups was expressed by the TPM ratio converted as the Log2 (fold change) (Log2FC).

The miRNA analysis by real-time quantitative polymerase chain reaction

Selection of patients

All women were postmenopausal and belonged to the OFELY cohort. The expression levels of selected miRNAs were measured, for the prevalent OA, in 43 women with prevalent knee OA (KL score of 2 and 3; early and intermediate knee OA, age 68.3 ± 6.6 years, body mass index (BMI) 26.6 ± 4.4 kg/m2) and in 42 healthy women without knee OA matched for age and BMI, and for the incident OA, in 23 women with incident knee OA over the next 4 years (age 68.4 ± 8 years, BMI 25.2 ± 4 kg/m2) and in 25 healthy subjects without incident knee OA matched for age and BMI (Table 1).

Table 1 Baseline characteristics of postmenopausal women included in the knee OA study with the localization of OA, if present at other sites than the knee

Preselection of miRNAs

The first selection of candidate miRNAs was based on the data obtained from the initial screening phase. We have selected miR-139-5p, miR-200a-3p, and miR-1299 in regard to their significant differential expression between the groups after Benjamini-Hochberg false discovery rate (FDR) correction. We have also added miR-375 for its high level of expression and a tendency to be differentially expressed, despite the absence of significant difference after FDR correction. Then, a second miRNA selection was performed by a critical assessment of literature data until December 2016 to highlight miRNAs that are differentially expressed in the biofluids of patients with knee OA vs controls. We have considered the criteria for patient selection, the size of the groups studied, the technical approaches used for screening and validation of miRNAs in biofluids, and the involvement of these differentially expressed miRNAs in the regulation of metabolic process contributing to OA pathology (see Table 2). Thus, we have further selected 15 candidate miRNAs: let-7e-5p, 16-5p, 29a-3p, 29b-3p, 29c-3p, 93-5p, 126-3p, 132-3p, 146a-5p, 184, 186-5p, 195-5p, 199a-3p, 345-5p, and 885-5p.

Table 2 A summary of studies on human miRNAs with dysregulated expression in patients suffering from OA, including the biological fluid tested, the number of patients with their baseline characteristics, the criteria for patient selection, the quantification methods, and the miRNAs significantly dysregulated

RNA isolation and miRNA RT-qPCR analysis

Total RNA was extracted from 200 μl of serum (miRCURY RNA isolation kit for biofluids, EXIQON, Danemark) by adding a lysis solution containing cel-miR-39-3p as exogenous control, followed by protein precipitation and miRNA purification on a silica column. miRNAs quantification was performed by the TaqMan® Advanced miRNA technology (Applied Biosystems, Thermo Fisher Scientific) in which miRNA was reverse-transcribed and cDNA pre-amplified with the TaqMan Advanced miRNA-cDNA synthesis kit, followed by analysis on pre-designed miRNA TaqMan arrays (TLDA) according to the manufacturer’s protocols (Applied Biosystems, CA, USA). Amplified cDNAs (15 μl) were mixed with 75 μl of TaqMan Fast Advanced Mastermix and water (60 μl), and 100 μl of the sample was added to each array tank of the card. The TLDA wells were pre-spotted with TaqMan®Advanced miRNA assays that allow the quantification in duplicate of 19 miRNAs (Table 3) by RT-qPCR reaction on a QuantStudio7 flex (Applied Biosystems) according to the manufacturer’s protocol. The Ct values recorded by the software Expression Suite were normalized by the Ct mean of 3 endogenous miRNAs for the relative quantification (RQ) of miRNA levels as fold change (FC) = Log2(2−ΔΔCT) (for protocol details, see Additional file 1: Data S1).

Table 3 Identification of the miRNAs analyzed in the validation phase and of the endogenous and exogenous miRNA normalizers. Each miRNA is identified by its NCBI accession number and sequence, according to miRBase v20. The pre-designed TaqMan Advanced miRNA Assays used to quantify the serum miRNA level are specified by their identification number (Applied Biosystems, Thermo Fisher Scientific)

Statistical analysis

A miRNA with a p value < 0.05 and a false discovery rate of 5% [Benjamini-Hochberg false discovery rate (FDR) correction for NGS approach] were considered as differentially expressed. Wilcoxon tests were used to compare the miRNA levels between women with and without OA because of the skewness of the data. We have examined the likelihood of OA, expressed as odds ratios (ORs) and 95% confidence intervals (CIs), per quartile increase in miRNA levels in a logistic regression model. All statistical analyses were performed using Stata 12 (Stata Corp LP, College Station, TX, USA).


Screening: serum miRNA profiling of patients with knee OA and controls

We identified 421 miRNAs with an expression level ≥ 1 TPM and 241 with an expression level ≥ 10 TPM. When we compared both groups, 22 miRNAs showed differential expression (p < 0.05) between controls and OA patients, 13 upregulated and 9 downregulated. After the Benjamini-Hochberg FDR correction, hsa-miR-139-5p, hsa-miR-1299, and hsa-miR-200a-3p remained significantly differently expressed between OA patients and controls (p < 0.05, FDR at 5%) (Table 4).

Table 4 Differential expression of the miRNA level in the serum from controls and OA patients, analyzed by next-generation sequencing

Validation: differential expression of candidate miRNAs in serum of patients with knee OA and controls

Prevalent knee OA and miRNA-146a-5p

When considered as a continuous variable, miR-146a-5p was significantly increased in the group of prevalent knee OA compared with controls (relative quantification (RQ); median [interquartile range] 1.12 [0.73; 1.46] vs 0.85 [0.62; 1.03], p = 0.015). The likelihood of prevalent knee OA was significantly increased (odds ratio [95% confidence interval (CI)] 1.83 [1.21–2.77], p = 0.004) for each quartile increase in serum miR-146a-5p. Moreover, the women with miR-146a-5p levels above the median (0.851) had a greater risk of prevalent knee OA compared to those below the median [95% CI] 4.62 [1.85–11.5], p = 0.001 (Table 5 and Fig. 2a and Additional file 2: Data S2 for the 19 miRs tested, mean and standard deviation (SD) for each quartile).

Table 5 Differential expression of the miRNA level in the serum from controls and OA patients analyzed by quantitative real-time PCR
Fig. 2
figure 2

Risk of prevalent and incident OA according to the median of miR-146a-5p and the quartiles of miR-186-5p, respectively. Multiple logistic regression analyses to determine the ability of serum miR-146a-5p and miR-186-5p levels to predict prevalent or incident OA, respectively. a Prevalent OA: percentage of women with osteoarthritis under and upper the median of miR-146a-5p (number of women with prevalent OA − total number of women). b Incident OA: percentage of women with osteoarthritis in the first three quartiles vs the upper quartile (number of women with incident OA − total number of women). CI, confidence interval; OD, odds ratio

Incident knee OA and miRNA-186-5p

When considered as a continuous variable, serum miR-186-5p levels were not significantly associated with an increased risk of incident knee OA (p = 0.09). However, in an analysis across quartiles, we found a significant association between baseline miR-186-5p levels and the likelihood of incident knee OA for each quartile increase (odds ratio [95% CI] 1.71 [1.00–2.95], p = 0.049). Women in the upper quartile were 6.13 times more likely to develop incident knee OA over the next 4 years compared to those in the three others quartiles (Q4 vs Q1–3; odds ratio [95% CI] 6.13 [1.14–32.9], p = 0.034) (Table 5 and Fig. 2b).


In the present work, we have studied the serum miRNA profiling of women with knee OA to identify potential non-invasive biomarkers of the disease. We first screened the serum miRNAs with significant differences in the expression between women with knee OA and healthy controls by a NGS approach. Then, we have conducted a validation study in the OFELY cohort of women having an incident or prevalent knee OA with a long-time follow-up. In this, we analyzed by RT-qPCR whether the 3 differentially expressed miRNAs (139-5p, 200a-3p, 1299) in patients with knee OA in the screening step were also aberrantly expressed in the women with incident and prevalent knee OA of the validation cohort. In addition, miR-375 has been chosen because of its high expression level and a tendency to be differentially expressed in our NGS approach despite the absence of significant difference after FDR correction. Moreover, we have measured the expression levels of 15 additional miRNAs previously highlighted as having a dysregulated expression in patients with knee OA. The miRs-29 family (a, b, c) and miR-199a-3p were chosen for their activity as key regulators of chondrocyte gene expression with aberrant expression in OA cartilage [29,30,31]. The miRNAs let-7e-5p, 16-5p, 93-5p, 126-3p, 132-3p, 146a-5p, 184, 186-5p, 195-5p, 345-5p, and 885-5p were selected because previous studies analyzing circulating miRNAs in the serum and plasma [14, 17, 18] have reported their dysregulated expression in patients with knee OA (see Table 4).

The NGS approach revealed that miR-139-5p, miR-1299, and miR-200a-3p had levels of expression significantly different between OA patients and controls. However, we did not validate these 3 miRNAs in the second step of the analysis when we measured their level in the largest validation cohort. This lack of validation might come from the small number of samples in the discovery step (10 OA vs 10 non-OA) even if this number is comparable to those used in previous studies (Beyer: 13 OA vs 13 non-OA [18]; Borgonio-Cuadra: 14 OA vs 5 non-OA [17]; Kong: 8 OA vs 8 non-OA [21]) and/or differences in the specificity and sensibility of the platforms used (Illumina for this study; Borgonio-Cuadra [17] and Beyer [18]: Applied Biosystems; Kong [21]: Affimetrix).

It appears that the previous studies performed with biological fluids are difficult to compare each other because of the heterogeneity of the clinical situations (see Table 2). Nevertheless, our study can be compared to that of Borgonio-Cuadra et al. [17], particularly the patient selection is based on the same criteria, radiologic KL scores of 2 and 3 and a BMI smaller than 27. However, among the ten miRNAs in common, only miR-146a-5p (miR-146) showed significant overexpression, emphasizing the challenges faced with patient heterogeneity. The reasons for this discrepancy could be the ethnic origin, Mexican vs European, enhancing the differences in exposure to environmental factors; the sex of the participants, men and women vs women only; the biological fluid used-plasma vs serum; or the mean age of participants, 55 vs 68 years.

Nevertheless, we found that miR-146a-5p expression is significantly increased in the OA group of postmenopausal women compared to healthy controls independently of age and BMI. Moreover, we reported for the first time that the risk of prevalent knee OA was significantly increased for each quartile increase in the serum level of miR-146a-5p. However, the interpretation of this increase remains challenging in regard to the dual role of miR-146a-5p in the pathophysiology of osteoarthritis. MiR-146 is expressed in chondrocytes which begin to undergo degenerative changes [32]; its expression is stimulated by pro-inflammatory cytokines in an NF-κB-dependent manner [33]. It may contribute to OA by impairing the TGF-β signaling pathway targeting Smad4 and increasing apoptosis [5, 34]. Recently, Zhang et al. showed that miR-146 facilitates OA progression by targeting Camk2d and Ppp3r2 that are required to maintain the phenotype of mature chondrocytes [33].

In contrast, miR-146 could play a protective role in OA. It is induced during hypoxia by HIF-1 and promotes autophagy by decreasing the Bcl-2 expression, an autophagy inhibitor. Autophagy has a protective function for cartilage during OA [35]. Moreover, miR-146 exerts negative control on inflammatory responses by suppressing IRAK1 and TRAF6, two molecular targets of miR-146 impairing NF-kβ activity and suppressing NF-kβ target genes such as IL-6 and TNF-α [36, 37]. Recently, Zhong et al. [38] showed that miR-146 could increase the proliferation and inhibit apoptosis of OA chondrocytes by inhibiting TRAF6 expression and suppressing the activation of the NF-κB signaling pathway. Finally, Guan et al. [39] showed that miR-146 is activated by cyclic loads in the physiological range but suppressed by mechanical overload in human articular chondrocytes. They proposed that miR-146 protects joint cartilage from degeneration through inhibition of the feedback loop of Notch-1/IL-6 during aging or trauma.

It must be noted that circulating miR-146 has been highlighted in two other studies analyzing the peripheral blood mononuclear cells. Okuhara et al. [16] found it overexpressed in patients suffering from knee OA compared to healthy controls and highly expressed in OA peripheral blood mononuclear cells from patients with early OA. Soyocak et al. [20] observed that, in patients classified according to the KL score, miR146a expression increased in the progressive stages (grades 3 and 4).

Finally, miR-146 is one of the four miRNAs with miR-223, miR-16, and miR-30b that were differently expressed in the cartilage and in the circulation with a role in the cartilage homeostasis [4]. Moreover, miR-146 is also involved in the inhibition of osteoclastogenesis [40] suggesting that its measurement in the serum could reflect the deregulations affecting the cartilage and bone tissues. Although the activity of this miRNA seems dependent on the disease phases, it appears as a key element in the OA pathophysiology, and our results suggest that miR-146a-5p could be a new biological marker for knee OA.

We report for the first time a significant association, independent of age and BMI, of serum miR-186-5p and incident knee OA risk in a population of postmenopausal women followed prospectively for 4 years. This result is in accordance with Borgonio-Cuadra et al. [17] showing that miR-186-5p was overexpressed in the plasma of early-stage OA patients compared to controls. In silico analysis has suggested that potential pathways regulated by miR-186 could be signaling by PDGF, developmental biology, membrane trafficking, and collagen formation. However, Li et al. [19] showed that miR-186 originated from the synovial membrane and belongs to a panel of seven miRNAs that were significantly upregulated in late-stage OA synovial fluid compared to early-stage irrespective of age, gender, and BMI. The reasons for this discrepancy are not clear but could be the geographic origin of patients, North America vs Europe, or the biological fluids analyzed, synovial fluid vs serum/plasma. Moreover, patients from Li’s study were selected from a group with degenerative meniscal tears undergoing arthroscopic surgery (early OA) or total knee replacement (late-stage OA) while Borgonio-Cuadra et al. [17] and us have selected patients on the radiographic score only, without any consideration regarding the OA etiology.

Together with the serum association of miR-146a-5p and prevalent OA, these data reinforce the notion that the miRNAs involved in OA pathophysiology vary according to the stage of the disease [19]. Finally, miR-186 has been essentially studied in cancer, Alzheimer’s disease, neuropathic disorder, and acute myocardial infarction [41]. It can be noticed that the last 3 diseases have been shown to be linked to OA. Recently, Weber et al. [42] demonstrated an association between OA and the risk of dementia. The miR-186 belonged to a set of 27 differentially expressed miRNAs between patients suffering from Alzheimer’s disease and controls [43]. The pathways potentially downregulated are linked with neuronal synaptic functions. Moreover, neuropathic pain is at least partially responsible for pain in OA [44]. Interestingly, miR-186 negatively regulated CXCL13, also known as B lymphocyte chemoattractant, and its downregulation in spinal neurons after spinal nerve ligation causes upregulation of CLCX13 to drive neuropathic pain [45]. Finally, several cross-sectional studies have shown a strong association between OA and cardiovascular diseases [46], and the overexpression of miR-186 in plasma have the potential to be used for the early detection of myocardial infarction [47]. However, for the present analysis, no woman suffered from dementia. Moreover, we had incomplete information regarding cardiovascular events, and we did not have data concerning pain in the OFELY study. So, the potential links between miR-186-5p, cardiovascular diseases, dementia, and OA could not be evaluated.

Our study has strengths and limitations. We investigated the association between the serum miRNAs expression and prevalence and for the first time incidence of knee OA in a well-characterized population of women followed prospectively over a long period of time of 4 years. We have chosen to carry out the validation step of our study by analysis of duplicate samples on TaqMan microRNA arrays, in order to reduce erratic manipulations and improve the reproducibility between the samples. Moreover, our study combines the validation of our NGS results and the first attempt to replicate previously published results, a process that had been lacking in the miRNA research field. Finally, we have indicated that women selected for the analyses of prevalent and incident knee OA had also OA at other sites (Table 1). This is a key point because we showed in the OFELY cohort that at this mean age, 68 ± 8 years, the prevalence of OA was high with 75% and 88% of women having osteophytosis at the lumbar spine and thoracic spine, respectively [48, 49] in agreement with previous prevalence data [50, 51]. Consequently, in both groups, knee OA and healthy women, we have verified that the presence of OA at other sites did not cause any bias between the groups (see Table 1) in order to decrease as much as possible the effects of this unavoidable situation of generalized OA on the detection of miRNAs associated with knee OA. However, we studied knee OA in women only, and the number of participants is limited. It has been demonstrated that gender may influence the level of the expression of miRNAs [52]. Consequently, our results need to be confirmed by additional larger studies and particularly in males. It can be noticed that miR-146 and miR-186 have been quantified in serum patients because the miR-3p forms of both miRNAs were almost undetectable in the 20 tested sera in the screening phase (data not shown). This lack of detection in the serum is probably related to thermodynamic instability of the 3p arms and to their subsequent degradation in the cellular p-bodies. Finally, the molecular mechanisms associating serum miR-146 and miR-186 to OA remain to be further clarified, for example, the role of the single-nucleotide polymorphisms in the promoter region of miR-146 gene [53].

Collectively, our results show that miR-146a-5p is increased in women with mild to moderate knee OA compared to healthy women. Importantly, miR-186-5p is also increased in those women who will develop radiographic knee OA over the next 4 years; therefore, this miRNA has the potential to detect preclinical knee OA.

Availability of data and materials

All data generated and analyzed during this study are included in this published article.


  1. Martel-Pelletier J, Barr AJ, Cicuttini FM, Conaghan PG, Cooper C, Goldring MB, et al. Osteoarthritis. Nat Rev Dis Primers. 2016;2:16072.

    Article  PubMed  Google Scholar 

  2. Rousseau J, Garnero P. Biological markers in osteoarthritis. Bone. 2012;51(2):265–77.

    Article  CAS  PubMed  Google Scholar 

  3. Shorter E, Sannicandro AJ, Poulet B, Goljanek-Whysall K. Skeletal muscle wasting and its relationship with osteoarthritis: a mini-review of mechanisms and current interventions. Curr Rheumatol Rep. 2019;21(8):40.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Trachana V, Ntoumou E, Anastasopoulou L, Tsezou A. Studying microRNAs in osteoarthritis: critical overview of different analytical approaches. Mech Ageing Dev. 2018;171:15–23.

    Article  CAS  PubMed  Google Scholar 

  5. Jin L, Zhao J, Jing W, Yan S, Wang X, Xiao C, et al. Role of miR-146a in human chondrocyte apoptosis in response to mechanical pressure injury in vitro. Int J Mol Med. 2014;34(2):451–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Meng F, Zhang Z, Chen W, Huang G, He A, Hou C, et al. MicroRNA-320 regulates matrix metalloproteinase-13 expression in chondrogenesis and interleukin-1beta-induced chondrocyte responses. Osteoarthr Cartil. 2016;24(5):932–41.

    Article  CAS  Google Scholar 

  7. Ji Q, Xu X, Xu Y, Fan Z, Kang L, Li L, et al. miR-105/Runx2 axis mediates FGF2-induced ADAMTS expression in osteoarthritis cartilage. J Mol Med (Berl). 2016;94(6):681–94.

    Article  CAS  Google Scholar 

  8. Mao G, Wu P, Zhang Z, Zhang Z, Liao W, Li Y, et al. MicroRNA-92a-3p regulates Aggrecanase-1 and Aggrecanase-2 expression in chondrogenesis and IL-1beta-induced catabolism in human articular chondrocytes. Cell Physiol Biochem. 2017;44(1):38–52.

    Article  PubMed  Google Scholar 

  9. Zhang D, Cao X, Li J, Zhao G. MiR-210 inhibits NF-kappaB signaling pathway by targeting DR6 in osteoarthritis. Sci Rep. 2015;5:12775.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Hu W, Zhang W, Li F, Guo F, Chen A. miR-139 is up-regulated in osteoarthritis and inhibits chondrocyte proliferation and migration possibly via suppressing EIF4G2 and IGF1R. Biochem Biophys Res Commun. 2016;474(2):296–302.

    Article  CAS  PubMed  Google Scholar 

  11. Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, et al. Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci U S A. 2011;108(12):5003–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Turchinovich A, Weiz L, Langheinz A, Burwinkel B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. 2011;39(16):7223–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. 2011;13(4):423–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Murata K, Yoshitomi H, Tanida S, Ishikawa M, Nishitani K, Ito H, et al. Plasma and synovial fluid microRNAs as potential biomarkers of rheumatoid arthritis and osteoarthritis. Arthritis Res Ther. 2010;12(3):R86.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Zhang L, Yang M, Marks P, White LM, Hurtig M, Mi QS, et al. Serum non-coding RNAs as biomarkers for osteoarthritis progression after ACL injury. Osteoarthr Cartil. 2012;20(12):1631–7.

    Article  CAS  Google Scholar 

  16. Okuhara A, Nakasa T, Shibuya H, Niimoto T, Adachi N, Deie M, et al. Changes in microRNA expression in peripheral mononuclear cells according to the progression of osteoarthritis. Mod Rheumatol. 2012;22(3):446–57.

    Article  CAS  PubMed  Google Scholar 

  17. Borgonio Cuadra VM, Gonzalez-Huerta NC, Romero-Cordoba S, Hidalgo-Miranda A, Miranda-Duarte A. Altered expression of circulating microRNA in plasma of patients with primary osteoarthritis and in silico analysis of their pathways. PLoS One. 2014;9(6):e97690.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Beyer C, Zampetaki A, Lin NY, Kleyer A, Perricone C, Iagnocco A, et al. Signature of circulating microRNAs in osteoarthritis. Ann Rheum Dis. 2015;74(3):e18.

    Article  PubMed  CAS  Google Scholar 

  19. Li YH, Tavallaee G, Tokar T, Nakamura A, Sundararajan K, Weston A, et al. Identification of synovial fluid microRNA signature in knee osteoarthritis: differentiating early- and late-stage knee osteoarthritis. Osteoarthr Cartil. 2016;24(9):1577–86.

    Article  CAS  Google Scholar 

  20. Soyocak A, Kurt H, Ozgen M, Turgut Cosan D, Colak E, Gunes HV. miRNA-146a, miRNA-155 and JNK expression levels in peripheral blood mononuclear cells according to grade of knee osteoarthritis. Gene. 2017;627:207–11.

    Article  CAS  PubMed  Google Scholar 

  21. Kong R, Gao J, Si Y, Zhao D. Combination of circulating miR-19b-3p, miR-122-5p and miR-486-5p expressions correlates with risk and disease severity of knee osteoarthritis. Am J Transl Res. 2017;9(6):2852–64.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Ntoumou E, Tzetis M, Braoudaki M, Lambrou G, Poulou M, Malizos K, et al. Serum microRNA array analysis identifies miR-140-3p, miR-33b-3p and miR-671-3p as potential osteoarthritis biomarkers involved in metabolic processes. Clin Epigenetics. 2017;9:127.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Arlot ME, Sornay-Rendu E, Garnero P, Vey-Marty B, Delmas PD. Apparent pre- and postmenopausal bone loss evaluated by DXA at different skeletal sites in women: the OFELY cohort. J Bone Miner Res. 1997;12(4):683–90.

    Article  CAS  PubMed  Google Scholar 

  24. Rousseau JC, Sornay-Rendu E, Bertholon C, Garnero P, Chapurlat R. Serum periostin is associated with prevalent knee osteoarthritis and disease incidence/progression in women: the OFELY study. Osteoarthr Cartil. 2015;23(10):1736–42.

    Article  CAS  Google Scholar 

  25. Peterfy C, Li J, Zaim S, Duryea J, Lynch J, Miaux Y, et al. Comparison of fixed-flexion positioning with fluoroscopic semi-flexed positioning for quantifying radiographic joint-space width in the knee: test-retest reproducibility. Skelet Radiol. 2003;32(3):128–32.

    Article  CAS  Google Scholar 

  26. Garnero P, Sornay-Rendu E, Arlot M, Christiansen C, Delmas PD. Association between spine disc degeneration and type II collagen degradation in postmenopausal women: the OFELY study. Arthritis Rheum. 2004;50(10):3137–44.

    Article  PubMed  Google Scholar 

  27. Kellgren JH, Lawrence JS. Radiological assessment of osteo-arthrosis. Ann Rheum Dis. 1957;16(4):494–502.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Gensburger D, Roux JP, Arlot M, Sornay-Rendu E, Ravaud P, Chapurlat R. Influence of blinding sequence of radiographs on the reproducibility and sensitivity to change of joint space width measurement in knee osteoarthritis. Arthritis Care Res (Hoboken). 2010;62(12):1699–705.

    Article  Google Scholar 

  29. Akhtar N, Haqqi TM. MicroRNA-199a* regulates the expression of cyclooxygenase-2 in human chondrocytes. Ann Rheum Dis. 2012;71(6):1073–80.

    Article  CAS  PubMed  Google Scholar 

  30. Ukai T, Sato M, Akutsu H, Umezawa A, Mochida J. MicroRNA-199a-3p, microRNA-193b, and microRNA-320c are correlated to aging and regulate human cartilage metabolism. J Orthop Res. 2012;30(12):1915–22.

    Article  CAS  PubMed  Google Scholar 

  31. Le LT, Swingler TE, Crowe N, Vincent TL, Barter MJ, Donell ST, et al. The microRNA-29 family in cartilage homeostasis and osteoarthritis. J Mol Med (Berl). 2016;94(5):583–96.

    Article  CAS  Google Scholar 

  32. Yamasaki K, Nakasa T, Miyaki S, Ishikawa M, Deie M, Adachi N, et al. Expression of MicroRNA-146a in osteoarthritis cartilage. Arthritis Rheum. 2009;60(4):1035–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Zhang X, Wang C, Zhao J, Xu J, Geng Y, Dai L, et al. miR-146a facilitates osteoarthritis by regulating cartilage homeostasis via targeting Camk2d and Ppp3r2. Cell Death Dis. 2017;8(4):e2734.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Li J, Huang J, Dai L, Yu D, Chen Q, Zhang X, et al. miR-146a, an IL-1beta responsive miRNA, induces vascular endothelial growth factor and chondrocyte apoptosis by targeting Smad4. Arthritis Res Ther. 2012;14(2):R75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Zhang F, Wang J, Chu J, Yang C, Xiao H, Zhao C, et al. MicroRNA-146a induced by hypoxia promotes chondrocyte autophagy through Bcl-2. Cell Physiol Biochem. 2015;37(4):1442–53.

    Article  CAS  PubMed  Google Scholar 

  36. Wang JH, Shih KS, Wu YW, Wang AW, Yang CR. Histone deacetylase inhibitors increase microRNA-146a expression and enhance negative regulation of interleukin-1beta signaling in osteoarthritis fibroblast-like synoviocytes. Osteoarthr Cartil. 2013;21(12):1987–96.

    Article  CAS  Google Scholar 

  37. Gu SX, Li X, Hamilton JL, Chee A, Kc R, Chen D, et al. MicroRNA-146a reduces IL-1 dependent inflammatory responses in the intervertebral disc. Gene. 2015;555(2):80–7.

    Article  CAS  PubMed  Google Scholar 

  38. Zhong JH, Li J, Liu CF, Liu N, Bian RX, Zhao SM, et al. Effects of microRNA-146a on the proliferation and apoptosis of human osteoarthritis chondrocytes by targeting TRAF6 through the NF-kappaB signalling pathway. Biosci Rep. 2017;37(2).

  39. Guan YJ, Li J, Yang X, Du S, Ding J, Gao Y, et al. Evidence that miR-146a attenuates aging- and trauma-induced osteoarthritis by inhibiting Notch1, IL-6, and IL-1 mediated catabolism. Aging Cell. 2018;17(3):e12752.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Nakasa T, Shibuya H, Nagata Y, Niimoto T, Ochi M. The inhibitory effect of microRNA-146a expression on bone destruction in collagen-induced arthritis. Arthritis Rheum. 2011;63(6):1582–90.

    Article  CAS  PubMed  Google Scholar 

  41. Islam F, Gopalan V, Vider J, Wahab R, Ebrahimi F, Lu CT, et al. MicroRNA-186-5p overexpression modulates colon cancer growth by repressing the expression of the FAM134B tumour inhibitor. Exp Cell Res. 2017;357(2):260–70.

    Article  CAS  PubMed  Google Scholar 

  42. Weber A, Mak SH, Berenbaum F, Sellam J, Zheng YP, Han Y, et al. Association between osteoarthritis and increased risk of dementia: a systemic review and meta-analysis. Medicine (Baltimore). 2019;98(10):e14355.

    Article  Google Scholar 

  43. Satoh J, Kino Y, Niida S. MicroRNA-Seq data analysis pipeline to identify blood biomarkers for Alzheimer’s disease from public data. Biomark Insights. 2015;10:21–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Dimitroulas T, Duarte RV, Behura A, Kitas GD, Raphael JH. Neuropathic pain in osteoarthritis: a review of pathophysiological mechanisms and implications for treatment. Semin Arthritis Rheum. 2014;44(2):145–54.

    Article  PubMed  Google Scholar 

  45. Jiang BC, Cao DL, Zhang X, Zhang ZJ, He LN, Li CH, et al. CXCL13 drives spinal astrocyte activation and neuropathic pain via CXCR5. J Clin Invest. 2016;126(2):745–61.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Nielen MM, van Sijl AM, Peters MJ, Verheij RA, Schellevis FG, Nurmohamed MT. Cardiovascular disease prevalence in patients with inflammatory arthritis, diabetes mellitus and osteoarthritis: a cross-sectional study in primary care. BMC Musculoskelet Disord. 2012;13:150.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Wang KJ, Zhao X, Liu YZ, Zeng QT, Mao XB, Li SN, et al. Circulating MiR-19b-3p, MiR-134-5p and MiR-186-5p are promising novel biomarkers for early diagnosis of acute myocardial infarction. Cell Physiol Biochem. 2016;38(3):1015–29.

    Article  CAS  PubMed  Google Scholar 

  48. Sornay-Rendu E, Munoz F, Duboeuf F, Delmas PD, Study O. Disc space narrowing is associated with an increased vertebral fracture risk in postmenopausal women: the OFELY Study. J Bone Miner Res. 2004;19(12):1994–9.

    Article  PubMed  Google Scholar 

  49. Sornay-Rendu E, Allard C, Munoz F, Duboeuf F, Delmas PD. Disc space narrowing as a new risk factor for vertebral fracture: the OFELY study. Arthritis Rheum. 2006;54(4):1262–9.

    Article  PubMed  Google Scholar 

  50. O’Neill TW, McCloskey EV, Kanis JA, Bhalla AK, Reeve J, Reid DM, et al. The distribution, determinants, and clinical correlates of vertebral osteophytosis: a population based survey. J Rheumatol. 1999;26(4):842–8.

    PubMed  Google Scholar 

  51. Liu G, Peacock M, Eilam O, Dorulla G, Braunstein E, Johnston CC. Effect of osteoarthritis in the lumbar spine and hip on bone mineral density and diagnosis of osteoporosis in elderly men and women. Osteoporos Int. 1997;7(6):564–9.

    Article  CAS  PubMed  Google Scholar 

  52. Wu J, Cai H, Xiang YB, Matthews CE, Ye F, Zheng W, et al. Intra-individual variation of miRNA expression levels in human plasma samples. Biomarkers. 2018;23(4):339–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Dole NS, Delany AM. MicroRNA variants as genetic determinants of bone mass. Bone. 2016;84:57–68.

    Article  CAS  PubMed  Google Scholar 

Download references


This work was supported by INSERM.

Author information

Authors and Affiliations



JCR and RC contributed to the conception and design. MC was responsible for the miRNA expertise. MM and OB contributed to the collection and assembly of the data. JCR, MC, ESR, and RC contributed to the analysis and interpretation of the data. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jean-Charles Rousseau.

Ethics declarations

Ethics approval and consent to participate

All subjects gave their informed consent, and the OFELY cohort has been approved by the local ethics committee (Comité de Protection des personnes Lyon Sud-Est 2).

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Data S1.

Protocol details for NGS analysis, RNA isolation and RT-qPCR analysis.

Additional file 2: Data S2.

Mean and standard deviation (SD) for each quartile of the 19 miRs tested in the validation phase.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rousseau, JC., Millet, M., Croset, M. et al. Association of circulating microRNAs with prevalent and incident knee osteoarthritis in women: the OFELY study. Arthritis Res Ther 22, 2 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: