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

Jean-Charles Rousseau, Marjorie Millet, Martine Croset, Elisabeth

Sornay-Rendu, Olivier Borel, Roland Chapurlat

► To cite this version:

Jean-Charles Rousseau, Marjorie Millet, Martine Croset, Elisabeth Sornay-Rendu, Olivier Borel, et al. Association of circulating microRNAs with prevalent and incident knee osteoarthritis in women: the OFELY study. Arthritis Research & Therapy, 2020, 22 (1), pp.2. 10.1186/s13075-019-2086-5. inserm-05022311

HAL Id: inserm-05022311

https://inserm.hal.science/inserm-05022311v1

Submitted on 6 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Association of circulating microRNAs with prevalent and incident knee osteoarthritis in women: the OFELY study

<figure><img src="image_2.png" /></figure>

Jean-Charles Rousseau<fnref n="1" />*, Marjorie Millet<fnref n="1" />, Martine Croset<fnref n="1" />, Elisabeth Sornay-Rendu<fnref n="1" />,², Olivier Borel<fnref n="1" />,² and Roland Chapurlat<fnref n="1" />,²,³

## Abstract

**Objectives:** 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.

**Methods:** 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.

**Results:** 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).

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

**Keywords:** Prevalent knee OA, Incident knee OA, Circulating miRNAs

## Key messages

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

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

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

## Introduction

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

<footnote>* Correspondence: jean-charles.rousseau@inserm.fr</footnote>

<footnote n="1">¹INSERM 1033, Pavillon F, Hôpital E. Herriot, 69437 Lyon Cedex 03, France
Full list of author information is available at the end of the article</footnote>
[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–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–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–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.

## Methods

### 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.
<figure><img src="image_3.png" /></figure>

### 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 10⁶. The differential expression of miRNA between the OA and control groups was expressed by the TPM ratio converted as the Log₂ (fold change) (Log₂FC).

### 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/m²) 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/m²)

and in 25 healthy subjects without incident knee OA matched for age and BMI (Table 1).

#### 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.

#### RNA isolation and miRNA RT-qPCR analysis

Total RNA was extracted from 200 µl of serum (miRCURY RNA isolation kit for biofluids, EXIQON,
<table><caption>**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</caption><thead><tr><th rowspan="3"></th><th rowspan="3"></th><th colspan="6">Quantitative real-time PCR</th></tr><tr><th colspan="3">Prevalent OA</th><th colspan="3">Incident OA</th></tr><tr><th>Control</th><th>Knee OA</th><th>p value</th><th>Control</th><th>Knee OA</th><th>p value</th></tr></thead><tbody><tr><td>Number of subjects</td><td></td><td>42</td><td>43</td><td></td><td>25</td><td>23</td><td></td></tr><tr><td>Age (years, mean ± SD)</td><td></td><td>68.3 ± 6.6</td><td>68.3 ± 6.6</td><td>0.77</td><td>67.7 ± 7.2</td><td>68.4 ± 8.0</td><td>0.72</td></tr><tr><td>BMI (kg/m²)</td><td></td><td>25.9 ± 4.4</td><td>26.6 ± 4.4</td><td>0.38</td><td>26.2 ± 5.2</td><td>25.2 ± 4.0</td><td>0.73</td></tr><tr><td rowspan="2">Knee OA</td><td>KL score at baseline</td><td>0</td><td>2/3</td><td></td><td>0</td><td>0/1</td><td></td></tr><tr><td>KL score at the end of the follow-up</td><td>0/1</td><td>2/3/4</td><td></td><td>0</td><td>2/3</td><td></td></tr><tr><td rowspan="3">OA at other sites</td><td>Hip (Lane score)</td><td>2</td><td>5</td><td>0.25</td><td>5</td><td>5</td><td>0.88</td></tr><tr><td>Hand (ACR criteria)</td><td>23</td><td>25</td><td>0.75</td><td>11</td><td>7</td><td>0.33</td></tr><tr><td>Spine (self-reported)</td><td>34</td><td>37</td><td>0.67</td><td>20</td><td>21</td><td>0.19</td></tr></tbody></table>

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) = Log₂(2⁻²⁻²ΔΔCT) (for protocol details, see Additional file 1: Data S1).

## 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).

## Results

### 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).

### 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).

#### 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).
<table><caption>**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</caption><thead><tr><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td>Murata K.<br>2010 [13]</td><td>Plasma and<br>synovial fluid</td><td>Plasma: 30 knee OA, 30 RA, 30 CTL<br>Synovial fluid: 30 OA, 30 RA</td><td>Plasma: OA: 75.1 yrs, 77 % ; RA: 60.1 yrs,<br>73%; CTL: 46.5 yrs, 57 %<br>Synovial fluid: RA: 63.1, 80%; OA: 75.3,<br>80%</td><td>RA and knee OA were diagnosed according to the ACR criteria<br>no indication of OA at other sites</td></tr><tr><td>Zhang L.<br>2012 [14]</td><td>Serum</td><td>Screening: 5 OA, 5 RA and 6 CTL<br>Validation: 102 ACL and 60 CTL</td><td>41 yrs or youngsters</td><td>one year after anterior cruciate ligament injury<br>no indication of OA at other sites</td></tr><tr><td>Okuhara A.<br>2012 [15]</td><td>Peripheral blood<br>mononuclear<br>cells</td><td>36 OA<br>36 CTL</td><td>OA: 68 yrs, 81 %<br>CTL: 32 yrs, 47 %</td><td>Knee OA ACR criteria<br>no indication of OA at other sites</td></tr><tr><td>Borgonio-<br>Cuadra VM,<br>2014 [16]</td><td>Plasma</td><td>Screening: 14 OA and 5 CTL<br>Validation: 27 OA and 27 CTL</td><td>Screening: OA: 55.7 yrs, 71.4 % CTL:<br>47.5 yrs, 100%<br>Validation: OA: 55.6 yrs, 88.9% CTL: 52.9<br>yrs, 81.5%</td><td>Knee OA, KL 2/3 and BMI &lt; 27<br>no indication of OA at other sites</td></tr><tr><td>Beyer C.<br>2015 [17]</td><td>Serum</td><td>Screening: pooled serum from 13 individuals with<br>knee/hip arthroplasty<br>Pooled serum from 13 individuals without knee/<br>hip arthroplasty<br>Validation: 749 OA and 67 CTL</td><td>OA: 65 yrs, 58.2% CTL: 62.7 yrs, 49.3%</td><td>Knee/hip arthroplasty (KL 3.4)<br>no indication of OA at other sites</td></tr><tr><td>Li YH, 2016<br>[18]</td><td>Synovial fluid</td><td>Screening: 4 early OA and 4 late OA<br>Validation: 22 early OA and 26 late OA</td><td>Screening: Early OA: 51 yrs, 100 % Late<br>OA: 64 yrs, 100%<br>Validation: Early OA: 56 yrs, 36.4% Late<br>OA: 60 yrs, 61.5%</td><td>Early stage OA: degenerative meniscal tears undergoing arthroscopic<br>surgery (KL grade 1.2)<br>Late stage OA: total knee replacement surgery (KL grade 3.4)<br>no indication of OA at other sites</td></tr><tr><td>Soyocak A.<br>2017 [19]</td><td>Peripheral blood<br>mononuclear<br>cell</td><td>100 patients with knee OA and 50 CTL</td><td>OA: from 47 to 70yrs, 84%<br>CTL: from 35 to 38 yrs, 84%</td><td>Knee OA ACR criteria<br>no indication of OA at other sites</td></tr><tr><td>Kong R.<br>2017 [20]</td><td>Plasma</td><td>Screening: 8 knee OA and 8 CTL<br>Validation: 100 OA and 100 CTL</td><td>Screening: OA: 51.13 yrs, 62.5% CTL:<br>50.75 yrs, 62.5%<br>Validation: OA: 51.69 yrs, 69% CTL:<br>51.09 yrs, 61%</td><td>Knee OA ACR criteria<br>no indication of OA at other sites</td></tr><tr><td>Ntoumou<br>E, 2017<br>[21]</td><td>serum</td><td>Screening: 12 primary OA and 12 CTL<br>Validation: 12 OA and 12 CTL</td><td>Screening: OA: 69.8 yrs, 75% CTL: 64.2<br>yrs, 50%<br>Validation: not indicated</td><td>Knee OA, KL ≥3<br>no indication of OA at other sites</td></tr><tr><td>Murata,<br>2010</td><td>High pure<br>miRNA isolation<br>kit (Roche)</td><td>Ncode VILO miRNA cDNA Synthesis kit<br>(Invitrogen)</td><td>Express SYBR GreenER qPCR Supermix<br>(Invitrogen)<br>control: cel-miR-39</td><td>Applied Biosystems<br>7300<br>SDS Relative<br>Quantification 1.3<br>(Applied Biosystems)</td></tr><tr><td>Zhang,<br>2012</td><td>miRNeasy kit<br>(Quiagen)</td><td>Validation: TaqMan miRNA reverse transcription kit<br>(Invitrogen) + a pulsed RT reaction with a<br>Eppendorf mastercycler (Eppendorf)</td><td>Screening: Megaplex RT human pool A<br>and B (Applied)<br>Validation: Preamp, TaqMan PreAmp<br>master mix; qPCR, TaqMan qPCR assays<br>control: U6 (Applied)</td><td>Plasma miRNAs had distinct pattern from SF<br>miRNAs<br>mir-132: potential diagnostic marker for patients<br>with OA or RA<br>U38 and U48 upregulated in patients developing<br>cartilage damage at one year after ACL injury</td></tr><tr><td>Okuhara,<br>2012</td><td>Trizol reagent</td><td>Thermocycler (BioRad)</td><td>TaqMan miRNA assay kit (Applied)</td><td>Screening: 7900HT<br>(Applied)<br>Validation: 7500<br>(Applied)<br>SDS Relative<br>Quantification 2.2.3<br>(Applied)</td></tr><tr><td></td><td></td><td></td><td>Mini Opticon Real-<br>146, 155, 181a, 223 upregulated in OA vs CTL</td><td></td></tr></tbody></table>
<table><caption>**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 (Continued)</caption><thead><tr><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td>2012</td><td>(Invitrogen)</td><td>Mini miRNeasy kit (Quiagen)</td><td>Screening: RT Megaplex Pool A on a GeneAmp PCR 9700 System (Applied) Validation: specific miRNA primer and TaqMan probes (Applied)</td><td>Screening: preamp with Megaplex PreAmp MasterMix, TLDA ver.2.0 plate A (Applied) Validation: preamp with Megaplex PreAmp MasterMix, qPCR control: U6</td><td>time PCR System (BioRad) Screening and Validation: 7900HT (Applied)</td><td>ealy stage: 146a and 223 higher than in late stage</td></tr><tr><td>Borgonio-Cuadra, 2014</td><td></td><td>Mini miRNeasy kit (Quiagen)</td><td>Megaplex Primer Pools (Human Pools A V.2.1)</td><td>Screening: Human TaqMan miRNA Array Card A V.2.1 (Applied) Validation: TaqMan miRNA assays (Applied) control: U6 or Ct average of all miRNA measurements for each sample</td><td>Screening and validation: 7900HT (Applied) SDS 2.2 software (Applied) and LinRegPCR software</td><td>12 miRNAs overexpressed in OA vs CTL: 16, 20b, 29c, 30b, 93, 126, 146a, 184, 186, 195, 345, 885-5p</td></tr><tr><td>Beyer, 2015</td><td></td><td>Mini miRNeasy kit (Quiagen)</td><td>Universal cDNA synthesis kit II (EXIQON)</td><td>miRNA ready-to-use PCR array (Human panel I + II, EXIQON) using ExiLENT SYBR Green master mix (EXIQON)</td><td>not indicated</td><td>let-7e, 454, 885-5p potential predictors for severe knee or hip OA</td></tr><tr><td>Li YH, 2016</td><td>miRCURY RNA isolation kit-biofluids (EXIQON)</td><td></td><td></td><td>TaqMan Small RNA Assays, TaqMan Gene Expression Master Mix control: U44 and 18S</td><td>qPCR system (Mx3000p, Stratagene)</td><td>23a-3p, 24-3p, 27b-3p, 29c-3p, 34a-5p, 186-5p upregulated and 27a-5p, 329, 655, 708-3p, 934 downregulated in late stage OA, severe early OA 27a-3p, 101-5p, 378-5p only expressed in late stage</td></tr><tr><td>Soyocak A, 2017</td><td>miRVana miRNA Isolation kit (Applied)</td><td>TaqMan MicroRNA Reverse transcription Kit (Applied)</td><td>Screening: microarray hybridation with the GeneChip miRNA 4.0 Array (Affymetrix) Validation: TaqMan microRNA Reverse Transcription kit (Life Technologies)</td><td>Validation: TaqMan miRNA assays (Applied) control: U6</td><td>miR-155: increased in OA miR-146a and miR-155 increased in the progressive stages</td><td></td></tr><tr><td>Kong R, 2017</td><td>LeukoLOCK kit (Ambion)</td><td></td><td></td><td></td><td>19b-3p, 92a-3p, 122-5p, 486-5p, 320b increased in OA</td><td></td></tr><tr><td>Ntroumou E, 2017</td><td>RiboEXTMLS kit (Geneall)</td><td></td><td>Screening: miRNA complete labeling and hybridization kit (Agilent) using SurePrint G3 Human miRNA 8X60K platform (Agilent) Validation: miScript II Reverse Transcription Kit (QIAGEN)</td><td>Validation: quantification with miScript SYBR Green PCR kit and miScript Primer Assays (QIAGEN) Control: Hsa-miR-25-1</td><td>19b-3p, 122-5p, 486-5p, great diagnostic value 19b-3p and 486-5p positively correlated with disease severity</td><td>3 miRNAs significantly downregulated in OA patients: hsa-miR-140-3p, hsa-miR-33b-3p, hsa-miR-671-3p</td></tr></tbody></table>

CTL control, OA osteoarthritis, yrs years, BMI Body Mass Index, KL Kellgren et Lawrence, ACR American College of Rheumatology, Preamp preamplification
<table><caption>**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)</caption><thead><tr><td>miR base ID</td><td>NCBI accession number</td><td>TaqMan™ Advanced miRNA Assay (ID)</td><td>Sequence of the mature miRNA (5'-3')</td></tr></thead><tbody><tr><td>hsa-miR-139-5p</td><td>MIMAT0000250</td><td>478312_mir</td><td>UCUACAGUGCACGUGUCUCCAGU</td></tr><tr><td>hsa-miR-200a-3p</td><td>MIMAT0000682</td><td>478490_mir</td><td>UAACACUGUCUGGUAAACGAUGU</td></tr><tr><td>hsa-miR-1299</td><td>MIMAT0005887</td><td>478696_mir</td><td>UUCUGGAAUUCUGUGUGAGGGA</td></tr><tr><td>hsa-let-7e-5p</td><td>MIMAT0000066</td><td>478579_mir</td><td>UGAGGUAGGAGGUUGUAUAGUU</td></tr><tr><td>hsa-miR-16-5p</td><td>MIMAT0000069</td><td>477860_mir</td><td>UAGCAGCACGUAAAUAUUGGCG</td></tr><tr><td>hsa-miR-29a-3p</td><td>MIMAT0000086</td><td>478587_mir</td><td>UAGCACCAUCUGAAAUCGGUUA</td></tr><tr><td>hsa-miR-29b-3p</td><td>MIMAT0000100</td><td>478369_mir</td><td>UAGCACCAUUUGAAAUCAGUGUU</td></tr><tr><td>hsa-miR-29c-3p</td><td>MIMAT0000681</td><td>479229_mir</td><td>UAGCACCAUUUGAAAUCGGUUA</td></tr><tr><td>hsa-miR-93-5p</td><td>MIMAT0000093</td><td>478210_mir</td><td>CAAAGUGCUGUUCGUGCAGGUAG</td></tr><tr><td>hsa-miR-126-3p</td><td>MIMAT0000445</td><td>477887_mir</td><td>UCGUACCGUGAGUAAAUAAUGCG</td></tr><tr><td>hsa-miR-132-3p</td><td>MIMAT0000426</td><td>477900_mir</td><td>UAACAGUCUACAGCCAUGGUCG</td></tr><tr><td>hsa-miR-146a-5p</td><td>MIMAT0000449</td><td>478399_mir</td><td>UGAGAACUGAAUUCCAUGGGUU</td></tr><tr><td>hsa-miR-184</td><td>MIMAT0000454</td><td>477938_mir</td><td>UGGACGGAGAACUGAUAAAGGU</td></tr><tr><td>hsa-miR-186-5p</td><td>MIMAT0000456</td><td>477940_mir</td><td>CAAAGAAUUCUCCUUUUGGCU</td></tr><tr><td>hsa-miR-195-5p</td><td>MIMAT0000461</td><td>477957_mir</td><td>UAGCAGCACAGAAAUUUGGC</td></tr><tr><td>hsa-miR-199a-5p</td><td>MIMAT0000231</td><td>478231_mir</td><td>CCCAGUGUUCAGACUACCUGUUC</td></tr><tr><td>hsa-miR-345-5p</td><td>MIMAT0000772</td><td>478366_mir</td><td>GCUGACUCCUAGUCCAGGGCUC</td></tr><tr><td>hsa-miR-375</td><td>MIMAT0000728</td><td>478074_mir</td><td>UUUGUUCGUUCGGUCGCGUGA</td></tr><tr><td>hsa-miR-885-5p</td><td>MIMAT0004947</td><td>478207_mir</td><td>UCCAUUACACUACCCUGCCUCU</td></tr><tr><td>hsa-miR-191-5p</td><td>MIMAT0000440</td><td>477952_mir</td><td>CAACGGAAUCCCAAAAGCAGCUG</td></tr><tr><td>hsa-miR-222-3p</td><td>MIMAT0000279</td><td>477982_mir</td><td>AGCUACAUCUGGCUACUGGGU</td></tr><tr><td>hsa-miR-361-5p</td><td>MI0000760</td><td>481127_mir</td><td>UUAUCAGAAUCUCCAGGGGUAC</td></tr><tr><td>cel-miR-39-3p</td><td>MI0000010</td><td>478293_mir</td><td>UCACCGGGUGUAAAUCAGCUUG</td></tr></tbody></table>

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

Discussion

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–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
<table><caption>**Table 4** Differential expression of the miRNA level in the serum from controls and OA patients, analyzed by next-generation sequencing</caption><thead><tr><th>Names</th><th>Log2 fold change</th><th>p value</th><th>FDR adjusted p value</th><th>Healthy average TPM</th><th>OA average TPM</th></tr></thead><tbody><tr><td>hsa-miR-139-5p</td><td>0.73</td><td>0.0001</td><td>0.0434</td><td>90.1</td><td>143.3</td></tr><tr><td>hsa-miR-1299</td><td>-3.38</td><td>0.0002</td><td>0.0434</td><td>12</td><td>0.8</td></tr><tr><td>hsa-miR-200a-3p</td><td>-1.88</td><td>0.0003</td><td>0.0473</td><td>77.2</td><td>29.4</td></tr><tr><td>hsa-miR-129-5p</td><td>0.81</td><td>0.0022</td><td>0.1640</td><td>3.6</td><td>6.5</td></tr><tr><td>hsa-miR-429</td><td>-1.43</td><td>0.0023</td><td>0.1640</td><td>5.5</td><td>2.5</td></tr><tr><td>hsa-miR-9-5p</td><td>0.83</td><td>0.0045</td><td>0.2410</td><td>3.2</td><td>6.1</td></tr><tr><td>hsa-miR-375</td><td>-0.93</td><td>0.0065</td><td>0.2795</td><td>737</td><td>411.3</td></tr><tr><td>hsa-miR-139-3p</td><td>0.61</td><td>0.0072</td><td>0.2820</td><td>48.8</td><td>69.7</td></tr><tr><td>hsa-miR-200b-3p</td><td>-1.18</td><td>0.0089</td><td>0.3193</td><td>28</td><td>15.7</td></tr><tr><td>hsa-miR-150-5p</td><td>0.72</td><td>0.0164</td><td>0.4619</td><td>173.6</td><td>252.6</td></tr><tr><td>hsa-miR-192-5p</td><td>-1.03</td><td>0.0171</td><td>0.4619</td><td>1837.8</td><td>1161.9</td></tr><tr><td>hsa-miR-155-5p</td><td>0.52</td><td>0.0193</td><td>0.4907</td><td>64.9</td><td>87.5</td></tr><tr><td>hsa-miR-15b-3p</td><td>-0.73</td><td>0.0247</td><td>0.5867</td><td>43.8</td><td>24.3</td></tr><tr><td>hsa-miR-1228-5p</td><td>0.78</td><td>0.0295</td><td>0.5918</td><td>4.9</td><td>8.7</td></tr><tr><td>hsa-miR-126-3p</td><td>0.38</td><td>0.0304</td><td>0.5918</td><td>3243.9</td><td>3976.7</td></tr><tr><td>hsa-miR-206</td><td>1.10</td><td>0.0338</td><td>0.6011</td><td>51.3</td><td>95.9</td></tr><tr><td>hsa-miR-199b-5p</td><td>-0.57</td><td>0.0373</td><td>0.6011</td><td>29.9</td><td>18.6</td></tr><tr><td>hsa-miR-1246</td><td>0.67</td><td>0.0375</td><td>0.6011</td><td>29.9</td><td>46.9</td></tr><tr><td>hsa-miR-10b-5p</td><td>0.46</td><td>0.0428</td><td>0.6390</td><td>3848.4</td><td>5039.8</td></tr><tr><td>hsa-miR-642a-3p</td><td>0.69</td><td>0.0444</td><td>0.6405</td><td>4.5</td><td>7.2</td></tr><tr><td>hsa-miR-186-5p</td><td>-0.45</td><td>0.0473</td><td>0.6405</td><td>777.7</td><td>534</td></tr><tr><td>hsa-miR-10a-5p</td><td>0.45</td><td>0.0493</td><td>0.6469</td><td>2613.1</td><td>3361.1</td></tr></tbody></table>

When we compared the OA vs healthy groups, 22 miRNAs showed differential expression (p < 0.05) between controls and OA patients. After Benjamini-Hochberg false discovery rate (FDR) correction, hsa-miR-139-5p, hsa-miR-1299, and hsa-miR-200a-3p remained significantly different between OA patients and controls (p < 0.05, FDR at 5%)

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
<table><caption>**Table 5** Differential expression of the miRNA level in the serum from controls and OA patients analyzed by quantitative real-time PCR</caption><thead><tr><th>miR</th><th>Prevalent OA (-), median (IQR 25-75%)</th><th>Prevalent OA (+), median (IQR 25-75%)</th><th>p value</th><th>Incident OA (-), median (IQR 25-75%)</th><th>Incident OA (+), median (IQR 25-75%)</th><th>p value</th></tr></thead><tbody><tr><td>126-3p</td><td>1.08 (0.85-1.36)</td><td>1.06 (0.89-1.30)</td><td>0.75</td><td>0.89 (0.80-1.04)</td><td>0.91 (0.73-1.23)</td><td>0.66</td></tr><tr><td>1299</td><td>0.28 (0.16-1.73)</td><td>0.33 (0.16-2.35)</td><td>0.80</td><td>0.27 (0.14-11.36)</td><td>0.34 (0.12-0.90)</td><td>0.82</td></tr><tr><td>132-3p</td><td>1.37 (0.71-2.12)</td><td>0.86 (0.43-1.85)</td><td>0.12</td><td>0.91 (0.63-1.19)</td><td>1.04 (0.59-1.97)</td><td>0.70</td></tr><tr><td>139-5p</td><td>1.50 (0.34-3.63)</td><td>2.09 (0.23-8.24)</td><td>0.27</td><td>0.28 (0.05-1.98)</td><td>0.81 (0.14-6.59)</td><td>0.16</td></tr><tr><td>146a-5p</td><td>0.85 (0.62-1.03)</td><td>1.12 (0.73-1.46)</td><td>0.015</td><td>0.83 (0.70-0.97)</td><td>0.88 (0.62-1.29)</td><td>0.94</td></tr><tr><td>16-5p</td><td>0.92 (0.72-1.82)</td><td>0.85 (0.61-1.50)</td><td>0.35</td><td>0.98 (0.59-1.24)</td><td>1.01 (0.73-1.44)</td><td>0.33</td></tr><tr><td>184</td><td>0.95 (0.64-2.18)</td><td>1.05 (0.58-2.38)</td><td>0.86</td><td>0.69 (0.41-1.76)</td><td>0.80 (0.31-1.24)</td><td>0.89</td></tr><tr><td>186-5p</td><td>1.09 (0.81-1.39)</td><td>1.03 (0.79-1.31)</td><td>0.47</td><td>0.82 (0.72-1.09)</td><td>1.05 (0.82-1.46)</td><td>0.09</td></tr><tr><td>195-5p</td><td>0.91 (0.75-1.42)</td><td>0.90 (0.52-1.60)</td><td>0.50</td><td>0.93 (0.65-1.24)</td><td>1.03 (0.69-1.47)</td><td>0.42</td></tr><tr><td>199a-3p</td><td>0.96 (0.77-1.12)</td><td>0.93 (0.82-1.47)</td><td>0.46</td><td>0.89 (0.76-0.98)</td><td>1.00 (0.71-1.39)</td><td>0.13</td></tr><tr><td>200a-3p</td><td>0.84 (0.56-1.93)</td><td>1.00 (0.58-1.57)</td><td>0.62</td><td>0.72 (0.49-1.06)</td><td>0.73 (0.48-1.43)</td><td>0.93</td></tr><tr><td>29a-3p</td><td>1.12 (0.76-1.44)</td><td>1.00 (0.68-1.38)</td><td>0.40</td><td>1.12 (0.63-1.35)</td><td>0.84 (0.61-1.50)</td><td>0.66</td></tr><tr><td>29b-3p</td><td>1.14 (0.81-1.74)</td><td>1.06 (0.81-1.51)</td><td>0.38</td><td>1.02 (0.71-1.42)</td><td>0.99 (0.67-1.53)</td><td>0.93</td></tr><tr><td>29c-3p</td><td>1.04 (0.66-1.59)</td><td>0.91 (0.64-1.59)</td><td>0.63</td><td>0.73 (0.44-1.16)</td><td>0.86 (0.45-1.56)</td><td>0.35</td></tr><tr><td>345-5p</td><td>3.10 (0.10-10.01)</td><td>2.81 (0.22-7.12)</td><td>0.99</td><td>0.53 (0.05-2.11)</td><td>1.14 (0.09-5.97)</td><td>0.24</td></tr><tr><td>375</td><td>1.57 (0.20-4.05)</td><td>1.40 (0.48-4.77)</td><td>0.56</td><td>0.98 (0.26-2.99)</td><td>0.86 (0.24-3.07)</td><td>0.91</td></tr><tr><td>885-5p</td><td>1.40 (0.60-2.16)</td><td>1.16 (0.28-2.40)</td><td>0.48</td><td>0.97 (0.49-2.85)</td><td>1.54 (0.70-2.38)</td><td>0.54</td></tr><tr><td>93-5p</td><td>1.11 (0.79-1.51)</td><td>0.91 (0.68-1.38)</td><td>0.21</td><td>1.07 (0.79-1.17)</td><td>1.19 (0.79-1.53)</td><td>0.25</td></tr><tr><td>let-7e-5p</td><td>1.15 (0.41-4.52)</td><td>2.28 (0.64-5.83)</td><td>0.25</td><td>1.36(0.78-2.72)</td><td>1.37 (0.36-2.36)</td><td>0.89</td></tr></tbody></table>

The miRNA levels were analyzed in the serum from individuals of the validation cohort, and differences in the serum miRNA levels are reported according to the prevalent or incident OA status. The relative miRNA level was used for the statistical comparisons. IQR = interquartile range [(+)] = women with prevalent or incident knee OA, (-) = controls

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-κβ activity and suppressing NF-κβ 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-
1A/ Prevalent OA, miR-146a-5p

% of prevalent OA

<figure><img src="image_4.png" /><figcaption>**Fig. 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</figcaption></figure>

1B/ Incident OA, miR-186-5p

% of incident OA

<figure><img src="image_5.png" /></figure>

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 chemo-attractant, 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.

**Supplementary information**

**Supplementary information** accompanies this paper at https://doi.org/10.1186/s13075-019-2086-5.

**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.

**Authors' contributions**

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.

**Funding**

This work was supported by INSERM.

**Availability of data and materials**

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

**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.

**Author details**

Received: 22 July 2019 Accepted: 9 December 2019

Published online: 02 January 2020

**References**

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

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

<footnote n="1">¹INSERM 1033, Pavillon F, Hôpital E. Herriot, 69437 Lyon Cedex 03, France.</footnote>

<footnote n="2">²Hôpital E. Herriot, Hospices Civils de Lyon, Lyon, France. ³Université de Lyon, Lyon, France.</footnote>
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.

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.

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

19. Li YH, Tavallaeae 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.

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.

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.

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.

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.

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

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.

26. Gamero 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.

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

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

## Publisher's Note

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

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_6.png" /></figure>