# Deregulation of microRNA expression in monocytes and CD4+ T lymphocytes from patients with axial spondyloarthritis

Olivier Fogel, Andreas Bugge Tinggaard, Maud Fagny, Nelly Sigrist, Elodie Roche, Laurence Leclere, Jean-François Deleuze, Frederic Batteux, Maxime Dougados, Corinne Miceli-Richard, et al.

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

## To cite this version:

Olivier Fogel, Andreas Bugge Tinggaard, Maud Fagny, Nelly Sigrist, Elodie Roche, et al. Deregulation of microRNA expression in monocytes and CD4+ T lymphocytes from patients with axial spondyloarthritis. Arthritis Research & Therapy, 2019, 21 (1), pp.1-14. 10.1186/s13075-019-1829-7. hal-03642295

HAL Id: hal-03642295

https://hal.inrae.fr/hal-03642295

Submitted on 19 Apr 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>
## RESEARCH ARTICLE

## Open Access

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

CrossMark

# Deregulation of microRNA expression in monocytes and CD4⁺ T lymphocytes from patients with axial spondyloarthritis

Olivier Fogel<fnref n="1" />,², Andreas Bugge Tinggaard<fnref n="1" />,³, Maud Fagny<fnref n="1" />, Nelly Sigrist<fnref n="1" />, Elodie Roche<fnref n="1" />, Laurence Leclere<fnref n="1" />, Jean-François Deleuze<fnref n="1" />, Frederic Batteux⁴, Maxime Dougados²,⁵,⁶, Corinne Miceli-Richard²,⁵† and Jörg Tost<fnref n="1" />*‡

### Abstract

**Background:** MicroRNAs (MiRs) play an important role in the pathogenesis of chronic inflammatory diseases. This study is the first to investigate miR expression profiles in purified CD4⁺ T lymphocytes and CD14⁺ monocytes from patients with axial spondyloarthritis (axSpA) using a high-throughput qPCR approach.

**Methods:** A total of 81 axSpA patients fulfilling the 2009 ASAS classification criteria, and 55 controls were recruited from October 2014 to July 2017. CD14⁺ monocytes and CD4⁺ T lymphocytes were isolated from peripheral blood mononuclear cells. MiR expression was investigated by qPCR using the Exiqon Human MiRnome panel I analyzing 372 miRNAs. Differentially expressed miRNAs identified in the discovery cohort were validated in the replication cohort.

**Results:** We found a major difference in miR expression patterns between T lymphocytes and monocytes regardless of the patient or control status. Comparing disease-specific differentially expressed miRs, 13 miRs were found consistently deregulated in CD14⁺ cells in both cohorts with miR-361-3p, miR-223-3p, miR-484, and miR-16-5p being the most differentially expressed. In CD4⁺ T cells, 11 miRs were differentially expressed between patients and controls with miR-16-1-3p, miR-28-5p, miR-199a-5p, and miR-126-3p were the most strongly upregulated miRs among patients. These miRs are involved in disease relevant pathways such as inflammation, intestinal permeability or bone formation. Mir-146a-5p levels correlated inversely with the degree of inflammation in axSpA patients.

**Conclusions:** We demonstrate a consistent deregulation of miRs in both monocytes and CD4⁺ T cells from axSpA patients, which could contribute to the pathophysiology of the disease with potential interest from a therapeutic perspective.

**Keywords:** Spondyloarthritis, MiRNA, Epigenetics, CD4⁺ T lymphocytes, Monocytes

### Background

Axial spondyloarthritis (axSpA) is a chronic inflammatory rheumatic disease that mainly affects the spine and sacroiliac joints in young adults. AxSpA includes both radiographic axSpA, also known as ankylosing spondylitis (AS), characterized by radiographic sacroiliitis according to the modified New York criteria [1], and non-radiographic axSpA (nr-axSpA), characterized by

the absence of structural damage of the sacroiliac joints on X-rays. The disease is frequently associated with extra-articular manifestations including psoriasis, uveitis, or inflammatory bowel disease [2, 3]. The long-term outcome of AxSpA is determined by structural damages such as the formation of syndesmophytes, excessive new bone formation of the spine [4].

Strong evidence suggests that disease pathogenesis is driven by innate immune cells (monocytes/macrophages, dendritic cells) and CD4⁺ T lymphocytes [3, 5-8]. Genetic variation associated with AxSpA was found enriched in open chromatin regions in monocytes and CD4⁺ T lymphocytes, suggesting these cell populations might

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

<footnote>* Correspondence: tost@cng.fr</footnote>

<footnote>†Corinne Miceli-Richard and Jörg Tost contributed equally to this work.</footnote>

<footnote n="1">¹Laboratory for Epigenetics and Environment, Centre National de Recherche en Génomique Humaine, CEA - Institut de Biologie François Jacob, 2 rue Gaston Crémieux, Evry, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
play an important role in the disease [9]. Genetic analyses identified the HLA-B27 and IL23R loci, whose proteins are present at the cell membrane of monocytes and T lymphocytes respectively, as the genetic factors most strongly associated with the disease [10, 11]. Monocytes/macrophages are crucially involved in the disease pathogenesis notably through the HLA-B27-induced unfolded protein response (UPR) stress leading to the release of pro-inflammatory cytokines (IL-1α, TNF, IL-6, and IL-23) [12]. Furthermore, T cells and macrophages have been found in the cellular infiltrates of tissue biopsies from AS patients [13–18]. Adoptive transfer studies from disease-prone B27 transgenic rats to B27 nude rats demonstrated that CD4⁺ T cells were able to induce the disease [7].

MicroRNAs (miRs) are small non-coding RNAs composed of 18–25 nucleotides that play an important regulatory role at the post-transcriptional level in diverse biological processes including cell differentiation, proliferation, apoptosis, or cellular function [19]. MiRs are involved in immune functions during granulopoiesis, T and B cell ontogenesis, TLR signaling, and cytokine production [20, 21]. We have previously shown that miR expression profiles are highly cell type specific with different expression profiles observed between CD4⁺ T and CD19⁺ B lymphocytes [22], reflecting their specific roles in a broad range of cellular functions. While miR deregulation has been comprehensively investigated in several autoimmune or rheumatic diseases including rheumatoid arthritis [23–25] and Sjögren's syndrome [22], similar studies on radiographic or non-radiographic axSpA are scarce.

Given the cell specificity of miR expression profiles and considering the biological importance of CD4⁺ T cells and monocytes in the pathophysiology of the disease, we investigated and validated miR expression profiles in these two cell types from axSpA patients and controls in two independent cohorts, totaling 81 patients and 55 controls.

# Methods

## Patients

Two independent cohorts of 22 and 59 patients with axSpA were recruited from October 2014 to July 2017 in the Department of Rheumatology at Cochin Hospital in Paris, France. All patients fulfilled the 2009 ASAS classification criteria for axSpA [2]. Demographic and phenotypic data, HLA-B27 status, imaging (X-rays and/or MRI assessment of sacroiliac joints), and treatment history were collected. Disease activity was assessed using both the Bath Ankylosing Spondylitis Disease Index (BASDAI) and the Ankylosing Spondylitis Disease Activity Score (ASDAS) [26, 27]. All patients were TNF-blocker

naïve or free of any biological treatment for more than 3 months.

The exploratory cohort included 22 patients and 17 age- and sex-matched controls. The replication cohort included 59 patients and 38 age and sex-matched controls (Fig. 1). The main characteristics of patients and controls for both cohorts are shown in Table 1. Briefly, the mean age was 40 ± 13 years at enrolment, 73% were male, 77% were HLA-B27 positive, and 70% had radiographic sacroiliitis. The median disease duration from the onset of symptoms was 4 (0–46) years. Patients exhibited high disease activity with a mean BASDAI score of 49 ± 19 and a mean CRP of 12.5 ± 17 mg/l. Most patients were under NSAIDs at inclusion (80%), but few were treated with csDMARDs (15%).

Controls were followed in the Department of Rheumatology for peripheral osteoarthritis without any history of inflammatory, infectious or tumoral disease; none had any immunosuppressive treatment. Written informed consent was obtained from all participants prior to enrolment, and the study was approved by the ethics committee Ile De France VII (number PP 14-039).

## Cell isolation and RNA isolation

PBMCs were isolated by density-gradient centrifugation using Unisep maxi+ tube (Eurobio, Courtaboeuf, France). Monocytes and CD4⁺ T lymphocytes were sorted from PBMCs by direct magnetic labeling with CD14⁺ and CD4⁺ microbeads according to the manufacturer's instructions (Miltenyi Biotec, Paris, France). MicroRNAs were extracted using a 2-step protocol using both the RNeasy minElute Cleanup kit and the miRNeasy mini kit (Qiagen, Hilden, Germany) to obtain a smallRNA-enriched fraction according to the manufacturer's protocol. Thirty nanograms of smallRNA were reverse transcribed into complementary DNA using the universal cDNA synthesis kit II following manufacturer's protocol (Exiqon, Vedbaek, Denmark).

## Analysis of microRNA expression by real-time qPCR

The expression levels of 372 microRNAs were screened using the Human miRnome panel I (Version 4R, Exiqon, Vedbaek, Denmark) as previously described [22]. After quality controls and data processing steps (see Additional file 1), we retained only miRs that were expressed in at least 60% of the samples. Of the 372 miR assessed by real-time qPCR, 193 and 223 miR passed this threshold when analyzing monocytes and T lymphocytes together, 204 and 243 miRs passed this threshold in monocytes alone and 202 and 218 in CD4⁺ T lymphocytes in the exploratory and replication cohort, respectively. Only the differentially expressed microRNAs identified in the exploratory cohort were further analyzed in the 59 patients and 38 controls from the
<figure><img src="image_5.png" /><figcaption>**Fig. 1** Description of the population studied and main results. The number of patients and controls with sufficient material used for the different analyses are shown in the graph. Seventeen parameters were investigated for association with the phenotype on the 81 patients (sex, age, disease duration, smoking status, HLA-B27, peripheral arthritis, enthesitis, uveitis, psoriasis, familial history of spondyloarthritis, radiographic sacro-iliitis, magnetic sacro-iliitis, CRP, BASDAI score, ASDAS, csDMARD intake, or NSAID intake). Differentially expressed miR expression were identified using a two-step process: Discovery in an exploratory cohort and confirmed in a replication cohort. MicroRNAs written in red are upregulated in the disease and those written in blue are downregulated in patients. A nominal *p* value < 0.05 was considered as significant</figcaption></figure>

validation cohort (Fig. 1). MiR-155 and miR-146a were added to the analysis of the second cohort because of their relevant role in key inflammatory pathways and inflammatory diseases [20, 28–31].

## Comparison with results from previous work

We compared the identified panel of differentially expressed miR with data previously published in the literature. The Pubmed database was searched for articles in English analysing miR expression in spondyloarthritis with the following key words “microrna” or “miR” and “ankylosing spondylitis” or “spondyloarthritis” (until December 31, 2017). Only the miRs validated by qPCR in an independent cohort within the same publication and replicated by an independent group were included in the analysis.

## Statistical analysis

Analyses of differential expression and hierarchical clustering were performed using GenEx software (v.6, Exiqon). MiRs differentially expressed between axSpA patients and controls were identified using an unpaired *t* test with a nominal *p* value < 0.05 considered as significant. The study of miR expression in function of clinical

parameters was performed on the 81 axSpA patients by combining the data from both cohorts using the R stat package (R Core Team 2017). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria). A Spearman's test was used to study the correlation between miR expression levels and quantitative clinical parameters such as ASDAS, BASDAI, or CRP. The Mann-Whitney *U* test or Kruskall-Wallis test were used to test for differences in the distribution of miR expression between classes for categorical clinical parameters. Results were corrected for multiple testing using the Benjamini and Hochberg (BH) stratified false discovery approach. We considered results with a FDR-corrected *p* value < 0.05 as significant. The level of expression of miRs differentially expressed in CD14⁺ monocytes of SpA patients was compared between patients not treated with csDMARDs and patients treated with sulfasalazine or methotrexate patients using a three-way Kruskall-Wallis test.

KEGG pathway enrichment analysis was performed among target genes of each significantly differentially expressed miR in monocytes or CD4⁺ T lymphocytes separately using the DIANA tools [32]. Pathways were significantly enriched if the *p* value was lower or equal
<table><caption>**Table 1** Demographic and disease-related characteristics of the population study</caption><thead><tr><th rowspan="2">Variables</th><th colspan="2">Overall population</th><th colspan="2">Exploratory cohort</th><th colspan="2">Replication cohort</th></tr><tr><th>Patients</th><th>Controls</th><th>Patients</th><th>Controls</th><th>Patients</th><th>Controls</th></tr></thead><tbody><tr><td>No. of cases/controls</td><td>81</td><td>55</td><td>22</td><td>17</td><td>59</td><td>38</td></tr><tr><td>Male, n (%)</td><td>59 (73)</td><td>33 (62)</td><td>16 (72)</td><td>10 (59)</td><td>43 (73)</td><td>23 (64)</td></tr><tr><td>Age (years)<sup>a</sup></td><td>40 ± 13</td><td>38 ± 8</td><td>41 ± 12</td><td>39 ± 9</td><td>39 ± 13</td><td>37 ± 8</td></tr><tr><td>Disease duration (years)<sup>b</sup></td><td>4 (0–46)</td><td></td><td>6.5 (0.2–46)</td><td></td><td>4 (0–44)</td><td></td></tr><tr><td>HLA-B27 positive, n (%)</td><td>58 (77)</td><td></td><td>16 (80)</td><td></td><td>42 (76)</td><td></td></tr><tr><td>Current smoker, n (%)</td><td>38 (47.5)</td><td></td><td>11 (50)</td><td></td><td>27 (46)</td><td></td></tr><tr><td>Clinical characteristics</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>AxSpA, n (%)</td><td>78 (97.5)</td><td></td><td>22 (100%)</td><td></td><td>56 (96)</td><td></td></tr><tr><td>AxSpA X-ray positive, n (%)</td><td>56 (70)</td><td></td><td>18 (82)</td><td></td><td>38 (65)</td><td></td></tr><tr><td>AxSpA MRI positive, n (%)</td><td>63 (86)</td><td></td><td>18 (86)</td><td></td><td>45 (87)</td><td></td></tr><tr><td>Peripheral arthritis</td><td>17 (21)</td><td></td><td>6 (27)</td><td></td><td>11 (19)</td><td></td></tr><tr><td>Enthesitis</td><td>40 (50)</td><td></td><td>11 (50)</td><td></td><td>29 (50)</td><td></td></tr><tr><td>IBD</td><td>4 (5)</td><td></td><td>0 (0)</td><td></td><td>4 (7)</td><td></td></tr><tr><td>Uveitis</td><td>28 (35)</td><td></td><td>8 (36)</td><td></td><td>20 (34)</td><td></td></tr><tr><td>Psoriasis</td><td>11 (14)</td><td></td><td>2 (9)</td><td></td><td>9 (16)</td><td></td></tr><tr><td>Baseline CRP, mg/l<sup>a</sup></td><td>12.5 ± 17</td><td></td><td>15 ± 24</td><td></td><td>12 ± 14</td><td></td></tr><tr><td>BASDAI M0<sup>a</sup></td><td>49 ± 19</td><td></td><td>42 ± 18</td><td></td><td>52 ± 19</td><td></td></tr><tr><td>ASDAS M0<sup>a</sup></td><td>3 ± 0.9</td><td></td><td>2.9 ± 0.8</td><td></td><td>3 ± 1</td><td></td></tr><tr><td>Treatment</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NSAIDs, last 3 months, n (%)</td><td>64 (80)</td><td></td><td>19 (86)</td><td></td><td>45 (78)</td><td></td></tr><tr><td>Steroids, last 3 months, n (%)</td><td>3 (3.75)</td><td></td><td>0 (0)</td><td></td><td>3 (5)</td><td></td></tr><tr><td>csDMARDS, last 3 months, n (%)</td><td>12 (15)</td><td></td><td>7 (32)</td><td></td><td>5 (9)</td><td></td></tr><tr><td>Methotrexate, n (%)</td><td>5 (42)</td><td></td><td>2 (29)</td><td></td><td>3 (60)</td><td></td></tr><tr><td>Sulfasalazine, n (%)</td><td>7 (58)</td><td></td><td>5 (71)</td><td></td><td>2 (40)</td><td></td></tr></tbody></table>

IBD inflammatory bowel disease, CRP C-reactive protein, BASDAI Bath Ankylosing Spondylitis Disease Activity Index, ASDAS Ankylosing Spondylitis Disease

Activity Score

ªMean ± SD

ᵇMedian (max-min)

to 0.05 after correction for multiple testing, and the number of target genes belonging to this pathway was greater or equal to 5.

## Results

### MiR expression in patients and controls is cell type specific

One hundred sixty-three of the 193 miRs and 203 of the 223 miRs that could be reliably detected were significantly differentially expressed between CD4⁺ T lymphocytes and monocytes in the exploratory and replication cohort, respectively, after correction for multiple testing (data not shown). Principal component analysis of miR expression showed two main clusters corresponding to monocytes and T lymphocytes from both patients and controls. Hierarchical clustering revealed a lack of segregation between patients and controls within each cell type at this level of analysis (Additional file 1: Figure S1).

### MiR expression signatures in CD14⁺ monocytes and CD4⁺ T lymphocytes from axSpA patients

In monocytes, 51 of the 204 detectable miRs were significantly differentially expressed in the exploratory cohort between patients and controls. In the replication cohort, we confirmed that 13 miRs were differentially expressed: 5 miRs were downregulated whereas 8 miRs were upregulated in monocytes from axSpA patients compared to controls (Table 2, Fig. 1). Among these miRs, miR-361-3p and miR-223-3p were the most significantly associated with the disease in the replication cohort (log2 fold change (log₂FC): −1.35, *p* = 5.80 × 10⁻⁵ and log₂FC: 1.27, *p* = 4.90 × 10⁻⁴, respectively) (Fig. 2a). The expression level of miR-146a-5p was significantly lower in monocytes from axSpA patients compared to controls (log₂FC: −1.31, *p* = 2.00 × 10⁻³).

In CD4⁺ T lymphocytes, 70 of the 202 detectable miRs were differentially expressed in the exploratory cohort. The analysis in the replication cohort confirmed
<table><caption>**Table 2** Differentially expressed microRNA in CD4⁺ or CD14⁺ cells between axSpA patients and controls</caption><thead><tr><th>MicroRNA</th><th>Cell type</th><th>Log2 fold change<sup>1</sup></th><th>exploratory cohort P value</th><th>replication cohort P value</th><th>Previously validated target genes<sup>2</sup></th><th>Pathways<sup>3</sup></th><th>Diseases</th></tr></thead><tbody><tr><td>miR-16-1-3p</td><td>CD4<sup>+</sup></td><td>1.85</td><td>0.04</td><td>1.30E-05</td><td>IGF1R, CREBBP, EXT1/2</td><td>Autophagy</td><td></td></tr><tr><td>miR-28-5p</td><td>CD4<sup>+</sup></td><td>1.15</td><td>3.06E-03</td><td>2.50E-04</td><td>MAPK1, B7-H3</td><td>Osteoblastic differentiation, JAK/STAT signaling pathway</td><td></td></tr><tr><td>miR-199a-5p</td><td>CD4<sup>+</sup></td><td>2.11</td><td>3.20E-06</td><td>8.90E-04</td><td>CEBP/α, HIF-1α, IKKβ, NFkB1, EXT1/2, SMAD4, WNT2</td><td>TGF-β, monocyte to macrophage differentiation, osteoclast differentiation, WNT signaling</td><td>AS [43]</td></tr><tr><td>miR-126-3p</td><td>CD4<sup>+</sup></td><td>1.74</td><td>6.57E-07</td><td>2.40E-03</td><td>IKB, VCAM-1, LRP6, DNMT1, CXCL12</td><td>NFkB, FcγR-mediated phagocytosis, HIF-1α signaling, PI3K/AKT signaling pathway</td><td>IBD [47], AS [40]</td></tr><tr><td>let-7d-3p</td><td>CD4<sup>+</sup></td><td>1.24</td><td>0.03</td><td>0.04</td><td></td><td></td><td></td></tr><tr><td>let-7b-5p</td><td>CD14<sup>+</sup></td><td>1.24</td><td>0.02</td><td>3.10E-03</td><td></td><td></td><td></td></tr><tr><td>miR-361-3p</td><td>CD4<sup>+</sup></td><td>−1.26</td><td>1.11E-03</td><td>4.80E-03</td><td>IFNβ1, TLR4, SOX9, IL17RC, ERAP2</td><td>PI3K/AKT signaling pathway, WNT signaling, Apoptosis</td><td>CD [72], AS [41]</td></tr><tr><td>miR-484</td><td>CD4<sup>+</sup></td><td>1.2</td><td>5.24E-04</td><td>0.02</td><td>IFNβ1, IL1RA</td><td>WNT signaling</td><td>CD [72]</td></tr><tr><td>miR-16-5p</td><td>CD4<sup>+</sup></td><td>1.11</td><td>9.89E-04</td><td>1.10E-03</td><td>IFNAR2, IL2, RUNX2</td><td>Osteoblast differentiation</td><td></td></tr><tr><td>miR-197-3p</td><td>CD4<sup>+</sup></td><td>1.52</td><td>3.77E-04</td><td>0.03</td><td>SMAD7, IL1R1, PDCD4, IL18</td><td>M1 macrophages polarization, T lymphocyte activation, TGF-β, PI3K/AKT<sup>4</sup> and NOD pathways</td><td>AS [41, 73], RA [74]</td></tr><tr><td>miR-181c-5p</td><td>CD4<sup>+</sup></td><td>−1.23</td><td>7.41E-03</td><td>0.04</td><td>IL2, CCL22, IL1α, BCL2</td><td>TGF-β, Osteoclast differentiation, Apoptosis</td><td></td></tr><tr><td>miR-874-3p</td><td>CD4<sup>+</sup></td><td>−1.18</td><td>5.27E-04</td><td>0.05</td><td>STAT3, IKKβ</td><td>Apoptosis, FcγR mediated phagocytosis, PI3K/AKT signaling pathway</td><td></td></tr><tr><td>miR-223-3p</td><td>CD14<sup>+</sup></td><td>1.27</td><td>0.05</td><td>4.90E-04</td><td>NLRP3, CLDN8, FOXO3, IL6</td><td>Osteoblastic differentiation, FcγR mediated phagocytosis</td><td>SpA [40, 76]</td></tr><tr><td>miR-146a-5p</td><td>CD14<sup>+</sup></td><td>−1.31</td><td>0.08</td><td>2.01E-03</td><td>NFkB, TRAF6, IRAK1, TLR2/4, STAT1, RORα, MAP3K7</td><td>Monocyte survival, macrophage differentiation, NFkB pathway, Epithelial cell adhesion</td><td></td></tr><tr><td>miR-15a-5p</td><td>CD14<sup>+</sup></td><td>−1.15</td><td>2.75E-04</td><td>4.60E-03</td><td>BCL2, NFkB1, CARD8, TNFSF9, IKKα, SMAD7, IFNγ</td><td>NFkB and TNFα pathways, cytokines production, myeloid lineage differentiation</td><td></td></tr><tr><td>miR-376a-3p</td><td>CD14<sup>+</sup></td><td>1.62</td><td>3.00E-06</td><td>6.05E-03</td><td>CDK2, SMAD7, ATG4C, IL7</td><td>PI3K/AKT and NOD-like signaling pathways, BCR and TCR signaling, NFkB, TLR signaling, osteoclast differenciation, Apoptosis</td><td>IBD [62, 75], psoriasis [37], RA [74]</td></tr><tr><td>miR-885-5p</td><td>CD14<sup>+</sup></td><td>1.42</td><td>5.57E-03</td><td>0.01</td><td>CDK2, CASP16, EGF1R, CLDN1</td><td>Cell cycle, Bacterial invasion of epithelial cells, TGF-β pathway, Autophagy</td><td></td></tr><tr><td>miR-574-3p</td><td>CD14<sup>+</sup></td><td>1.35</td><td>1.78 E-03</td><td>1.89E-03</td><td>RXRα, IRF4, EGFR</td><td>Cell cycle, Epithelial cell adhesion, Bacterial invasion of epithelial cells</td><td>FMF [77]</td></tr></tbody></table>
<table><caption>**Table 2** Differentially expressed microRNA in CD4⁺ or CD14⁺ cells between axSpA patients and controls (Continued)</caption><thead><tr><th>MicroRNA</th><th>Cell type</th><th>Log2 fold change<sup>1</sup></th><th>exploratory cohort P value</th><th>replication cohort P value</th><th>Previously validated target genes*</th><th>Pathways**</th><th>Diseases</th></tr></thead><tbody><tr><td>miR-92b-3p</td><td>CD14<sup>+</sup></td><td>1.34</td><td>0.02</td><td>0.04</td><td>SMAD2/5, BMP7, MAPK1</td><td>mediated proteolysis</td><td></td></tr><tr><td>miR-30e-5p</td><td>CD14<sup>+</sup></td><td>-1.18</td><td>9.06E-04</td><td>0.05</td><td>BCL2, ATG12</td><td>Osteoblastic differentiation, TGF-B pathway</td><td></td></tr><tr><td colspan="8"><sup>1</sup> Ankylosing spondylitis, CD Crohn disease, <sup>IBD</sup> inflammatory bowel disease, <sup>RA</sup> rheumatoid arthritis, <sup>FMF</sup> familial Mediterranean fever</td></tr><tr><td colspan="8"><sup>2</sup> In the replication cohort</td></tr><tr><td colspan="8"><sup>3</sup> Experimentally validated target genes on miRWalk 2.0 website (http://zmf.umm.uni-heidelberg.de/apps/zmf/mirwalk2/) [78]</td></tr><tr><td colspan="8"><sup>4</sup> Pathway analysis using DIANA-mirpath [32] or search on PubMed</td></tr></tbody></table>
<figure><img src="image_6.png" /><figcaption>**Fig. 2** Distribution of the expression of the most differentially expressed microRNAs between axSpA patients and controls in CD14⁺ monocytes (a) and in CD4⁺ T lymphocytes (b) in the exploratory and replication cohorts. In monocytes, miR-484, miR-16-5p, and miR-223-3p were significantly upregulated and miR-361-3p was downregulated in patients (grey) compared to controls (white). In CD4⁺ T lymphocytes, the four most significantly deregulated miRs were upregulated in patients (grey) compared to controls (white). \*p < 0.05, \*\*p < 0.01, \*\*\*p < 0.001</figcaption></figure>
11 differentially expressed miRs, of which 3 were down-regulated and 8 upregulated in CD4⁺ T lymphocytes from axSpA patients compared with controls (Table 2, Fig. 1). MiR-16-1-3p and miR-28-5p were the miRs most significantly associated with the disease (log₂FC: 1.85, *p* = 1.30 × 10⁻⁵ and log₂FC: 1.15, *p* = 2.50 × 10⁻⁴, respectively) (Fig. 2b). An upregulation of miR-484, miR-126-3p, and miR-16-5p and a downregulation of miR-361-3p were found in the two cell types in the same direction suggesting a broad impact on the immune response.

## Key pathways

The bioinformatic analysis using DIANA tools identified an enrichment of deregulated miRs in the extracellular membrane-receptor interaction pathway and fatty acid biosynthesis and metabolism (Additional file 1: Figure S2). Some miRs may be of particular interest because of their role in key pathways of the disease. Mir-199a-5p, mir-126-3p, mir-16-5p, mir-197-3p, mir-223-3p, mir-146a-5p, mir-15a-5p, mir-376a-3p, and miR-92b-3p are involved in the immune response through the NF-κB, TGFβ, and NOD pathways or by modulating monocyte activation or macrophage polarization (Table 2). Bone remodeling may be also altered by miRs regulating the osteoblastic or osteoclast differentiation or through the WNT signaling pathway. Finally, mir-223-3p, mir-885, and mir-376a-3p could act on the interplay between the microbiota and the immune response by regulating the permeability of the intestinal epithelium during bacterial invasion (Table 2).

## Analysis of relationships between miRNA expression and clinical characteristics of the patients

We further explored whether the differentially expressed validated miRs between patients and controls were associated with some specific patient's characteristics. Twenty parameters were independently analysed covering clinical and biological parameters of the patients from both cohorts (Table 1). Inflammatory bowel disease (IBD) or corticosteroid intake was not included in the analyses as the effectives were not well-balanced for these two parameters. MiRs differentially expressed according to the disease's phenotype with a nominal *p* < 0.05 are reported in Table 3. After Benjamini-Hochberg correction for multiple testing, only age, the ASDAS score, and the CRP value in CD14⁺ cells remained significantly associated with the expression of specific miRs. The level of miR-146a-5p was negatively associated with patients who had a positive CRP at least once in their disease history (log₂FC: −0.3; FDR = 0.02) and inversely correlated with the CRP levels (*r* = −0.3; FDR = 0.04) and ASDAS score (*r* = −0.40; FDR = 0.007) (Fig. 3). However, no association was found between any miR expression and the BASDAI score or the

presence of a radiographic or magnetic sacro-iliitis. MiR-92b-3p was positively correlated with CRP value (*r* = 0.36, *p* = 0.01) and miR-126-3p and miR-574-3p were correlated with patients' age (*r* = −0.31, *p* = 0.05 and *r* = 0.29, *p* = 0.05 respectively). The use of NSAIDs did not impact miR expression while there was a trend toward a potential effect of csDMARD on three miRs (Table 3). When comparing patients not treated with csDMARDs to those treated with methotrexate (MTX) or sulfasalazine (SLZ), only miR-223-3p in CD14⁺ cells was found differentially expressed, with a decreased expression level in patients treated with sulfasalazine.

## Comparison of the differentially expressed miRs between our results and previously published literature

To provide further evidence for a specific miR signature in SpA, we compared the validated miRs found deregulated in our study with previous reports on miRNA changes in axSpA. Of the initially identified 43 studies, 16 studies were finally retained (see Additional file 1: Figure S3 for details). When analysing results on monocytes and on T lymphocytes separately, 53 miRs were found deregulated in the disease in at least one study, but only 12 miRs were found in at least two studies. Three miRs were found consistently and significantly upregulated in patients compared to controls across different tissue samples: miR-126-3p, miR-484, and miR-92b-3p (Additional file 1: Figure S4). MiR-16-5p was found to be upregulated in two studies analysing T cells, including ours.

## Discussion

In this large-scale miR profiling study analysing two independent cohorts, we showed that patients with axSpA present solid evidence for a deregulation of specific miRs in monocytes (*n* = 13) and CD4⁺ T lymphocytes (*n* = 11). Only four miRs were consistently differentially expressed between patients and controls in both CD4⁺ and CD14⁺ cells, demonstrating two distinct specific signatures of miR expression profiles in monocytes and CD4⁺ T lymphocytes. The major strength of this study is the large number of miRs screened in sorted cells from a large number of clinically well-characterized subjects.

To date, miR expression profiles have been poorly studied in SpA. Two different approaches have been pursued: (i) the first strategy consisting in analysing miR expression in sorted cells and (ii) the second strategy consisting in analysing circulating miRs isolated from serum or plasma. The first approach is mainly performed to advance our understanding of cell type-specific mechanisms involved in disease pathophysiology; the second aims mainly at identifying biomarkers of the disease in easily accessible biological material for clinical routine. However, miR expression profiles obtained from one strategy cannot be
<table><caption>**Table 3** Association between miRs expression and clinical or biological parameters in CD14⁺ and CD4⁺ cells from axSpA patients</caption><thead><tr><th rowspan="2">Parameters</th><th colspan="3">CD14<sup>+</sup></th><th colspan="3">CD4<sup>+</sup></th></tr><tr><th>MicroRNA</th><th>p value</th><th>FDR corrected p value</th><th>MicroRNA</th><th>p value</th><th>FDR corrected p value</th></tr></thead><tbody><tr><td rowspan="3">Age<sup>2</sup></td><td>miR-126-3p</td><td>4.45E-03</td><td>0.05</td><td rowspan="3">miR-16-5p</td><td rowspan="3">0.05</td><td rowspan="3">0.49</td></tr><tr><td>miR-574-3p</td><td>8.45E-03</td><td>0.05</td></tr><tr><td>miR-146a-5p</td><td>0.04</td><td>0.18</td></tr><tr><td>B27</td><td>miR-885-5p</td><td>0.03</td><td>0.46</td><td></td><td></td><td></td></tr><tr><td>Male</td><td>miR-223-3p</td><td>0.03</td><td>0.37</td><td></td><td></td><td></td></tr><tr><td>Disease duration<sup>2</sup></td><td>miR-92b-3p</td><td>0.03</td><td>0.44</td><td></td><td></td><td></td></tr><tr><td>Uveitis</td><td>miR-92-3p</td><td>0.02</td><td>0.42</td><td></td><td></td><td></td></tr><tr><td rowspan="2">Psoriasis</td><td>miR-15a-5p</td><td>0.01</td><td>0.10</td><td rowspan="2">miR-16-1-3p</td><td rowspan="2">0.04</td><td rowspan="2">0.39</td></tr><tr><td>miR-484</td><td>0.01</td><td>0.10</td></tr><tr><td rowspan="2">Familial history of SpA</td><td>miR-126-3p</td><td>0.02</td><td>0.18</td><td rowspan="2">miR-361-3p</td><td rowspan="2">0.04</td><td rowspan="2">0.28</td></tr><tr><td>miR-574-3p</td><td>0.03</td><td>0.18</td></tr><tr><td>Axial involvement</td><td></td><td></td><td></td><td>miR-874-3p</td><td>0.04</td><td>0.45</td></tr><tr><td>Peripheral arthritis</td><td>miR-16-5p</td><td>0.03</td><td>0.43</td><td>miR-361-3p</td><td>5.76E-03</td><td>0.06</td></tr><tr><td rowspan="2">Elevated CRP<sup>1</sup></td><td>miR-146a-5p</td><td>1.82E-03</td><td>0.02</td><td rowspan="2">miR-484</td><td rowspan="2">0.04</td><td rowspan="2">0.34</td></tr><tr><td>miR-361-3p</td><td>0.04</td><td>0.29</td></tr><tr><td rowspan="4">CRP value<sup>2</sup></td><td>miR-92b-3p</td><td>1.06E-03</td><td>0.01</td><td rowspan="4">miR-874-3p</td><td rowspan="4">0.04</td><td rowspan="4">0.39</td></tr><tr><td>miR-146a-5p</td><td>6.64E-03</td><td>0.04</td></tr><tr><td>miR-484</td><td>0.02</td><td>0.09</td></tr><tr><td>miR-574-3p</td><td>0.03</td><td>0.09</td></tr><tr><td rowspan="2">BASDAI > 40</td><td>miR-219a-5p</td><td>0.03</td><td>0.21</td><td rowspan="2">miR-28-5p</td><td rowspan="2">0.01</td><td rowspan="2">0.11</td></tr><tr><td>miR-361-3p</td><td>0.03</td><td>0.21</td></tr><tr><td rowspan="2">BASDAI score<sup>2</sup></td><td>miR-361-3p</td><td>0.01</td><td>0.14</td><td></td><td></td><td></td></tr><tr><td>miR-146a-5p</td><td>5.98E-04</td><td>7.78E-03</td><td></td><td></td><td></td></tr><tr><td rowspan="2">ASDAS score<sup>2</sup></td><td>miR-92b-3p</td><td>0.01</td><td>0.11</td><td></td><td></td><td></td></tr><tr><td>let-7b-5p</td><td>0.05</td><td>0.30</td><td></td><td></td><td></td></tr><tr><td rowspan="3">Current smoker</td><td>miR-223-3p</td><td>0.05</td><td>0.30</td><td></td><td></td><td></td></tr><tr><td>miR-376a-3p</td><td>0.02</td><td>0.09</td><td></td><td></td><td></td></tr><tr><td>miR-126-3p</td><td>0.02</td><td>0.09</td><td></td><td></td><td></td></tr><tr><td>csDMARDs</td><td>miR-223-3p</td><td>0.01</td><td>0.09</td><td></td><td></td><td></td></tr></tbody></table>

Only associations with a nominal *p* value < 0.05 are reported

MicroRNAs with a *p* value < 0.05 after Benjamin-Hochberg correction are written in italic

extrapolated to the other, as multiple cell types and mechanisms contribute to the prevalence and abundance of miRs present as cell-free circulating miRs.

We decided to investigate the miR-expression profile of CD4⁺ T lymphocytes and monocytes because of their relevance for the pathophysiology of axSpA [3, 5, 9]. We found that the main factor influencing the miR expression profile is the cell type: almost 90% of the measured miRs were differentially expressed between CD4⁺ and CD14⁺ cells in our study, regardless of the disease status. Thus, the assessment of miR expression profiles in PBMCs or whole blood may be challenging as any

observed changes could be related to variable cell subset composition, requiring thus the assessment of miR expression in sorted cells.

Some studies have addressed the impact of csDMARDs on miR expression in psoriasis or rheumatoid arthritis (RA). However, current data is conflicting. MiR-10a, miR-132-3p, miR-146a-5p, and miR-155-5p measured in whole blood or serum have been shown to be increased after MTX treatment in RA [33, 34] whereas miR-16 in plasma was not altered by the use of MTX in juvenile idiopathic arthritis [35]. In contrast, miR-146a measured in whole blood from patients with

<footnote n="1">¹At least one positive CRP value according to the laboratory reference during the disease course due to the disease's activity</footnote>

<footnote n="2">²Spearman test was used for quantitative values</footnote>
<figure><img src="image_7.png" /><figcaption>**Fig. 3** Association between the level of miR-146a-5p and the CRP level or ASDAS score in monocytes from axSpA patients. Negative correlation between the level of miR-146a-5p and **a.** the ASDAS score and **b.** the CRP levels, respectively. **c.** The miR-146a-5p level was significantly higher in patients who have never had a positive CRP value according to the laboratory threshold during the course of their disease. \*p < 0.05</figcaption></figure>

psoriasis was decreased after 3 months of methotrexate [36]. Also, in psoriasis, miR-223-3p was found downregulated in PBMC, but not altered in serum after the use of MTX and decreased in serum of patients treated by Etanercept [37, 38]. The authors suggested that the downregulation of miR-223 was not related to the treatment but rather to the immunosuppressive effect of the treatment on cells. We found only a single study reporting the impact of sulfasalazine on miR expression. Mir-21 was found to be negatively associated with *PDCD4* expression in whole blood from SpA patients not treated with csDMARDs, while positively associated with the expression of the same gene in patients treated with sulfasalazine [39]. In our study, we did not find miR-21 to be differentially expressed. Considering the presence of multiple cell types in most analyses, it is currently difficult to conclude if the use of csDMARDs has a direct influence on the miR expression profile or leads to a shift in cell populations. In our sorted cell populations, we found little evidence for the modification of the miRNA expression profiles by csDMARDs. However, our absolute number of patients treated with DMARDs were low, and larger studies will be required to fully address this question.

Several miRs found deregulated in CD4⁺ T cells or in monocytes in our study have been previously shown to be deregulated in other cell types in AS. Overexpression of mir-92b-3p, mir-16-5p, and mir-126-3p were also found overexpressed in T cells, fibroblasts, or PBMCs by others [40–42], reinforcing their likely involvement in the disease. No other study has focused on miR expression in monocytes from AS patients. Contrary to our result, Wang et al. found a lower expression of miR-199a-5p in T cells from 41 AS patients compared to healthy controls [43]. The main difference with our study is that miR expression was assessed on CD3⁺ T cells comprising both CD4⁺ and CD8⁺ T cells.

Mir-16-5p was found upregulated in both monocytes and T lymphocytes in our study and in CD3⁺ T lymphocytes from AS patients in a previous study [41]. It was assigned an anti-inflammatory function as it binds to the 3'-untranslated region (UTR) of *Tnfa* and *Il12b*, repressing their expression [44], or to the 3'UTR of *Pdcd4* (programmed cell death), resulting in a decreased level of Il-6 or Tnfa in mouse macrophages [45]. Nevertheless, miR-16 has also been shown to be able to create a pro-inflammatory environment [46].
The functional role of miR-126-3p has been well studied in IBD, where it was found upregulated in the colon tissue of patients with active ulcerative colitis [47, 48]. Furthermore, IκB, an inhibitor of NFκB, is a downstream target of miR-126-3p and the transfection of a miR-126-3p mimic decreases the expression of IκB in a dose-dependent manner [48]. MiR-126 also impairs intestinal mucosal barrier function under LPS stimulation by targeting *S1PR2* through the PI3k/AKT signaling pathway [49]. Functional experiments are required to further clarify its role in axSpA.

Macrophages play a major role in SpA as demonstrated by their abundant infiltration on biopsies from sacroiliac joints of AS patients, within the synovial tissue in chronic inflammatory arthritis, and in the gut mucosa [13, 50]. The macrophage infiltrate was shown to correlate with disease activity and to decrease after anti-TNF treatment [51, 52]. Macrophages represent a heterogeneous population of classically activated macrophages (M1), the main producers of inflammatory mediators such as TNFα, which are specialized in the clearance of pathogens and alternative M2 macrophages, which are involved in tissue repair and immunosuppression [53]. Whether M1 or M2 are causally implicated in SpA is still a matter of debate. CD163, a marker of M2 macrophages was found increased in the intimal lining layer of the synovial membrane and in gut mucosa from SpA patients [16, 50, 54]. In AS, the IFN-γ signature classically seen in M1 macrophages was impaired and the level of M1-derived cytokines was decreased in synovial fluid compared to patients with rheumatoid arthritis [17, 55]. However, these studies were focused on peripheral SpA and little is known regarding the role of M1 or M2 macrophages in axSpA. In PBMCs, a M2-like predominant polarization of monocytes has been reported, together with a M2/M1 ratio negatively correlated with CRP levels or BASDAI score [56]. A shift toward a M2 macrophage profile could facilitate bacterial survival and growth, especially for intracellular pathogens, such as Chlamydia or Yersinia, that have been involved in reactive arthritis and axial and peripheral arthritis in SKG mice [57, 58].

MiRs play a critical role in controlling macrophage differentiation [59]. MiR-16 was able to induce differentiation of mouse peritoneal macrophages into M1 macrophages from either the basal M0 or M2 polarized state and to activate purified CD4⁺ T lymphocytes, thus creating a pro-inflammatory environment [46]. In contrast, miR-124 and miR-223-3p were reported to contribute to M2 polarization or to limit M1 activation [60, 61]. The increased expression of miR-223-3p in our analysis suggests that circulating monocytes should be more prone to evolve into M2 macrophages in peripheral tissue. The function of miR-223-3p in SpA

has never been addressed before, but it was found to be up-regulated in psoriasis and in inflammatory bowel disease [37, 62–64]. A dual role has been attributed to miR-223-3p in IBD: a protective role through the decreased in the release of inflammatory mediators and a facilitating role by increasing intestinal permeability. The overexpression of miR-223-3p in RAW 264.7 macrophages inhibited LPS-stimulated secretion of Il-6 or Il-1 by targeting *STAT3* [60]. In H. pylori-infected THP-1 monocytes, increased expression of miR-223-3p was correlated with decreased expression of NLRP3 [65]. MiR-223⁻/⁻ mice exhibit exacerbated DSS-induced colitis and an increased level of IL-1β through the loss of repressive function on *NLRP3* targeted by miR-223-3p [64]. MiR-223-3p may also increase the intestinal permeability by targeting *CLDN8* in response to IL-23 [62]. Therapeutic approaches in mice model have been elaborated targeting this dual role of miR-223. A systemic treatment with an antagomiR-223 was able to attenuate the phenotype in TNBS-induced colitis in mice [62]. By contrast, an intravenous nanoparticle delivery of a miR-223 mimic attenuated DSS-colitis in mice by targeting *Nlrp3* [64]. The upregulation of miR-223-3p among SpA patients in our study could both reflect the activated state of the monocyte/macrophages cells and an induced anti-inflammatory regulatory loop.

MiR-146a-5p is the most studied miR in autoimmune or inflammatory diseases. Also called anti-inflammoMir, miR-146a can bind to many targets involved in inflammatory signaling such as IRAK1, TLR2/4 or NFκB. It acts as an effective regulator to prevent an overstimulated inflammatory state [31, 66]. IRAK1⁻/⁻ CD4⁺ T cells are less prone to produce IL17 under TGF-B and IL-6 stimuli and expressed a higher level of Foxp3 under TGF-B stimulation compared to WT CD4⁺ T cells [67]. Thus, IRAK1 can control the Th17/regulatory T-cells balance under specific stimuli, promoting Th17 differentiation. Keratinocytes from miR-146a⁻/⁻ mice were more sensitive to IL-17A stimulation that keratinocytes from WT mice, suggesting that mir-146a is a potent suppressor of IL17-driven disease [68]. MiR-146a-5p expression was found to be increased in skin biopsies and PBMCs from patients with psoriasis [69]. In SpA, discordant results have been published with an overexpression in fibroblasts and a decreased expression in PBMCs [70, 71]. We found a downregulation of miR-146a-5p in monocytes and its level was negatively correlated with ASDAS and CRP. Given its target genes, miR-146a-5p deregulation in monocytes might facilitate the polarization of Th17 cells from naïve CD4⁺ T cells through the release of pro-inflammatory cytokines.
## Conclusions

We have demonstrated a cell type-specific miRNA expression profile, whose deregulation could be related to the pathophysiology of the disease. More specifically, this study has highlighted some pathways that could be crucial in the onset or the maintenance of the disease such as the role of macrophage polarization, the impact of intestinal permeability or the modulation of the expression of pro-inflammatory cytokines.

## Additional files

**Additional file 1:** MicroRNAs expression by real-time qPCR. **Figure S1.** Principal component analysis of microRNAs expression and hierarchical clustering analysis in the exploratory cohort (A, C) and in the replication cohort (B, D). **Figure S2.** Pathways analysis of significant differentially expressed miRs in monocytes (A) and CD4⁺ T lymphocytes (B). **Figure S3.** Flow chart for the selection of published study for the comparison analysis. **Figure S4.** Cluster analysis of differentially expressed miRs in published studies from literature. (DOCX 2104 kb)

### Abbreviations

AS: Ankylosing spondylitis; ASDAS: Ankylosing Spondylitis Disease Activity Score; axSpA: Axial spondyloarthritis; BASDAI: Bath Ankylosing Spondylitis Disease Index; BH: Benjamini and Hochberg; CRP: C-reactive protein; csDMARD: Conventional synthetic disease modifying anti rheumatic drug; IBD: Inflammatory bowel disease; miR: MicroRNA; MTX: Methotrexate; NSAID: Non steroid anti-inflammatory drug; SLZ: Sulfasalazine

### Acknowledgements

The authors thank Nathalie Menage (Department of Rheumatology, Cochin hospital) for help in collecting clinical data and Sandrine Chouzenoux for help with the cell sorting.

### Funding

The research described in this manuscript was mainly funded by a Sirius research award (UCB Pharma S.A.). Additional funding was obtained through a Passerelle research award (Pfizer), iCARE (MSD Avenir) and the institutional budget of the CNRGH. O.F. is the recipient of a PhD fellowship from Assistance Publique - Hopitaux de Paris / CEA (poste accueil AP-HP).

### Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

### Authors' contributions

CMR and JT designed the study and provided funding; OF, ABT, LL, FB, MD, and CMR collected samples; OF, ABT, NS, ER, LL, FB, JFD processed samples and performed experiments; OF, ABT, MF, and JT analysed the data. OF wrote the first draft of the manuscript; JFD, MD, CMR, and JT revised the manuscript. All authors read and approved the final manuscript.

### Ethics approval and consent to participate

Written informed consent was obtained from all participants prior to enrolment and the study was approved by the ethics committee Ile De France VII (number PP 14–039).

### Consent for publication

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

## Publisher's Note

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

## Author details

¹Laboratory for Epigenetics and Environment, Centre National de Recherche en Génomique Humaine, CEA - Institut de Biologie François Jacob, 2 rue Gaston Crémieux, Evry, France. ²Department of Rheumatology - Hôpital Cochin. Assistance Publique - Hôpitaux de Paris, Paris Descartes University, Paris, France. ³Department of Biomedicine, Aarhus University, Aarhus, Denmark. ⁴Department of Immunology, Cochin Hospital, Paris, France. ⁵Unité Mixte AP-HP/ Institut Pasteur, Institut Pasteur, Immunoregulation Unit, Paris, France. ⁶INSERM (U1153) : Clinical Epidemiology and Biostatistics, PRES Sorbonne Paris-Cité, Paris, France.

Received: 9 November 2018 Accepted: 18 January 2019

Published online: 12 February 2019

## References

1. van der Linden S, Valkenburg HA, Cats A. Evaluation of diagnostic criteria for ankylosing spondylitis. A proposal for modification of the New York criteria. Arthritis Rheum. 1984;27:361–8.

2. Rudwaleit M, van der Heijde D, Landewé R, Listing J, Akkoc N, Brandt J, et al. The development of Assessment of SpondyloArthritis international Society classification criteria for axial spondyloarthritis (part II): validation and final selection. Ann Rheum Dis. 2009;68:777–83.

3. Sieper J, Braun J, Dougados M, Baeten D. Axial spondyloarthritis. Nat Rev Dis Primer. 2015;1:15013.

4. Van Mechelen M, Gulino GR, de Vlam K, Lories R. Bone disease in axial spondyloarthritis. Calcif Tissue Int. 2018;102:547–58.

5. Ranganathan V, Gracey E, Brown MA, Inman RD, Haroon N. Pathogenesis of ankylosing spondylitis - recent advances and future directions. Nat Rev Rheumatol. 2017;13:359–67.

6. Wendling D, Guillot X, Prati C. The IL-23/Th 17 pathway in spondyloarthritis: the royal road? Jt Bone Spine Rev Rhum. 2015;82:1–4.

7. Breban M, Fernández-Sueiro JL, Richardson JA, Hadavand RR, Maika SD, Hammer RE, et al. T cells, but not thymic exposure to HLA-B27, are required for the inflammatory disease of HLA-B27 transgenic rats. J Immunol Baltim Md 1950. 1996;156:794–803.

8. Conrad K, Wu P, Sieper J, Syrbe U. In vivo pre-activation of monocytes in patients with axial spondyloarthritis. Arthritis Res Ther. 2015;17:179.

9. Farh KK-H, Marson A, Zhu J, Kleinewietfeld M, Housley WJ, Beik S, et al. Genetic and epigenetic fine mapping of causal autoimmune disease variants. Nature. 2015;518:337–43.

10. Brown MA, Kenna T, Wordsworth BP. Genetics of ankylosing spondylitis—insights into pathogenesis. Nat Rev Rheumatol. 2016;12:81–91.

11. Ellinghaus D, Jostins L, Spain SL, Cortes A, Bethune J, Han B, et al. Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci. Nat Genet. 2016;48:510–8.

12. Colbert RA, Tran TM, Layh-Schmitt G. HLA-B27 misfolding and ankylosing spondylitis. Mol Immunol. 2014;57:44–51.

13. Braun J, Bollow M, Neure L, Seipelt E, Seyrekbasan F, Herbst H, et al. Use of immunohistologic and in situ hybridization techniques in the examination of sacroiliac joint biopsy specimens from patients with ankylosing spondylitis. Arthritis Rheum. 1995;38:499–505.

14. Bollow M, Fischer T, Reisshauer H, Backhaus M, Sieper J, Hamm B, et al. Quantitative analyses of sacroiliac biopsies in spondyloarthropathies: T cells and macrophages predominate in early and active sacroiliitis - cellularity correlates with the degree of enhancement detected by magnetic resonance imaging. Ann Rheum Dis. 2000;59:135–40.

15. Appel H, Kuhne M, Spiekermann S, Ebhardt H, Grozdanovic Z, Köhler D, et al. Immunohistologic analysis of zygapophyseal joints in patients with ankylosing spondylitis. Arthritis Rheum. 2006;54:2845–51.

16. Baeten D, Demetter P, Cuvelier CA, Kruithof E, Van Damme N, De Vos M, et al. Macrophages expressing the scavenger receptor CD163: a link between immune alterations of the gut and synovial inflammation in spondyloarthropathy. J Pathol. 2002;196:343–50.

17. Vandooren B, Noordenbos T, Ambarus C, Krausz S, Cantaert T, Yeremenko N, et al. Absence of a classically activated macrophage cytokine signature in peripheral spondylarthritis, including psoriatic arthritis. Arthritis Rheum. 2009;60:966–75.

18. Ciccia F, Alessandro R, Rizzo A, Accardo-Palumbo A, Raimondo S, Raiata F, et al. Macrophage phenotype in the subclinical gut inflammation of patients with ankylosing spondylitis. Rheumatol Oxf Engl. 2014;53:104–13.
19. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281-97.

20. Pauley KM, Cha S, Chan EKL. MicroRNA in autoimmunity and autoimmune diseases. J Autoimmun. 2009;32:189-94.

21. Duroux-Richard I, Presumey J, Courts G, Gay S, Gordeladze J, Jorgensen C, et al. MicroRNAs as new player in rheumatoid arthritis. Jt Bone Spine Rev Rhum. 2011;78:17-22.

22. Wang-Renault S-F, Boudaoud S, Nocturne G, Roche E, Sigrist N, Daviaud C, et al. Deregulation of microRNA expression in purified T and B lymphocytes from patients with primary Sjögren's syndrome. Ann Rheum Dis. 2018;77:133-40.

23. Pauley KM, Satoh M, Chan AL, Bubb MR, Reeves WH, Chan EK. Upregulated miR-146a expression in peripheral blood mononuclear cells from rheumatoid arthritis patients. Arthritis Res Ther. 2008;10:R101.

24. Stanczyk J, Pedrioli DML, Brentano F, Sanchez-Pernaute O, Kolling C, Gay RE, et al. Altered expression of MicroRNA in synovial fibroblasts and synovial tissue in rheumatoid arthritis. Arthritis Rheum. 2008;58:1001-9.

25. Fulci V, Scappucci G, Sebastiani GD, Giannitti C, Franceschini D, Meloni F, et al. miR-223 is overexpressed in T-lymphocytes of patients affected by rheumatoid arthritis. Hum Immunol. 2010;71:206-11.

26. Garrett S, Jenkinson T, Kennedy LG, Whitelock H, Gaisford P, Calin A. A new approach to defining disease status in ankylosing spondylitis: the Bath Ankylosing Spondylitis Disease Activity Index. J Rheumatol. 1994;21:2286-91.

27. Lukas C, Landewé R, Sieper J, Dougados M, Davis J, Braun J, et al. Development of an ASAS-endorsed disease activity score (ASDAS) in patients with ankylosing spondylitis. Ann Rheum Dis. 2009;68:18-24.

28. Li L, Chen X-P, Li Y-J. MicroRNA-146a and human disease. Scand J Immunol. 2010;71:227-31.

29. Testa U, Pelosi E, Castelli G, Labbaye C. miR-146 and miR-155: two key modulators of immune response and tumor development. Non-coding RNA. 2017;3(3):22.

30. Kurowska-Stolarska M, Alivernini S, Ballantine LE, Asquith DL, Millar NL, Gilchrist DS, et al. MicroRNA-155 as a proinflammatory regulator in clinical and experimental arthritis. Proc Natl Acad Sci U S A. 2011;108:11193-8.

31. Taganov KD, Boldin MP, Chang K-J, Baltimore D. NF-kappaB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. Proc Natl Acad Sci U S A. 2006;103:12481-6.

32. Vlachos IS, Kostoulas N, Vergoulis T, Georgakilas G, Reczko M, Maragkakis M, et al. DIANA miRPath v.2.0: investigating the combinatorial effect of microRNAs in pathways. Nucleic Acids Res. 2012;40:W498-504.

33. Singh A, Patro PS, Aggarwal A. MicroRNA-132, miR-146a, and miR-155 as potential biomarkers of methotrexate response in patients with rheumatoid arthritis. Clin Rheumatol. Clin Rheumatol. 2018. https://doi.org/10.1007/s10067-018-4380-z.

34. Hong H, Yang H, Xia Y. Circulating miR-10a as predictor of therapy response in rheumatoid arthritis patients treated with methotrexate. Curr Pharm Biotechnol. 2018;19:79-86.

35. Demir F, Çebi AH, Kalyoncu M. Evaluation of plasma microRNA expressions in patients with juvenile idiopathic arthritis. Clin Rheumatol. 2018;37:3255-62.

36. Ele-Refaei AM, El-Esawy FM. Effect of narrow-band ultraviolet B phototherapy and methotrexate on microRNA (146a) levels in blood of psoriatic patients. Dermatol Res Pract. 2015;2015:145769.

37. Løvendorf MB, Zibert JR, Gyldenløve M, Røpke MA, Skov L. MicroRNA-223 and miR-143 are important systemic biomarkers for disease activity in psoriasis. J Dermatol Sci. 2014;75:133-9.

38. Pivarci A, Meisgen F, Xu N, Ståhle M, Sonkoly E. Changes in the level of serum microRNAs in patients with psoriasis after antitumour necrosis factor-a therapy. Br J Dermatol. 2013;169:563-70.

39. Huang C-H, Wei JC-C, Chang W-C, Chiou S-Y, Chou C-H, Lin Y-J, et al. Higher expression of whole blood microRNA-21 in patients with ankylosing spondylitis associated with programmed cell death 4 mRNA expression and collagen cross-linked C-telopeptide concentration. J Rheumatol. 2014;41:1104-11.

40. Lv Q, Li Q, Zhang P, Jiang Y, Wang X, Wei Q, et al. Disorders of MicroRNAs in peripheral blood mononuclear cells: as novel biomarkers of ankylosing spondylitis and provocative therapeutic targets. Biomed Res Int. 2015;2015:504208.

41. Lai N-S, Yu H-C, Chen H-C, Yu C-L, Huang H-B, Lu M-C. Aberrant expression of microRNAs in T cells from patients with ankylosing spondylitis contributes to the immunopathogenesis. Clin Exp Immunol. 2013;173:47-57.

42. Zhang C, Wang C, Jia Z, Tong W, Liu D, He C, et al. Differentially expressed mRNAs, lncRNAs, and miRNAs with associated co-expression and ceRNA networks in ankylosing spondylitis. Oncotarget. 2017;8:113543-57.

43. Wang Y, Luo J, Wang X, Yang B, Cui L. MicroRNA-199a-5p induced autophagy and inhibits the pathogenesis of ankylosing spondylitis by modulating the mTOR signaling via directly targeting Ras homolog enriched in brain (Rheb). Cell Physiol Biochem Int J Exp Cell Physiol Biochem Pharmacol. 2017;42:2481-91.

44. Huang Z, Ma J, Chen M, Jiang H, Fu Y, Gan J, et al. Dual TNF-α/IL-12p40 interference as a strategy to protect against colitis based on miR-16 precursors with macrophage targeting vectors. Mol Ther J Am Soc Gene Ther. 2015;23:1611-21.

45. Liang X, Xu Z, Yuan M, Zhang Y, Zhao B, Wang J, et al. MicroRNA-16 suppresses the activation of inflammatory macrophages in atherosclerosis by targeting PDCD4. Int J Mol Med. 2016;37:967-75.

46. Jia X, Li X, Shen Y, Miao J, Liu H, Li G, et al. MiR-16 regulates mouse peritoneal macrophage polarization and affects T-cell activation. J Cell Mol Med. 2016;20:1898-907.

47. Wu F, Zikusoka M, Trindade A, Dassopoulos T, Harris ML, Bayless TM, et al. MicroRNAs are differentially expressed in ulcerative colitis and alter expression of macrophage inflammatory peptide-2 alpha. Gastroenterology. 2008;135:1624-1635.e24.

48. Feng X, Wang H, Ye S, Guan J, Tan W, Cheng S, et al. Up-regulation of microRNA-126 may contribute to pathogenesis of ulcerative colitis via regulating NF-kappaB inhibitor IkBa. PLoS One. 2012;7:e52782.

49. Chen T, Xue H, Lin R, Huang Z. MiR-126 impairs the intestinal barrier function via inhibiting S1PR2 mediated activation of PI3K/AKT signaling pathway. Biochem Biophys Res Commun. 2017;494:427-32.

50. Demetter P, De Vos M, Van Huysse JA, Baeten D, Ferdinande L, Peeters H, et al. Colon mucosa of patients both with spondyloarthritis and Crohn's disease is enriched with macrophages expressing the scavenger receptor CD163. Ann Rheum Dis. 2005;64:321-4.

51. Baeten D, Kruithof E, De Rycke L, Boots AM, Mielants H, Veys EM, et al. Infiltration of the synovial membrane with macrophage subsets and polymorphonuclear cells reflects global disease activity in spondyloarthropathy. Arthritis Res Ther. 2005;7:R359-69.

52. Kruithof E, Baeten D, Van den Bosch F, Mielants H, Veys EM, De Keyser F. Histological evidence that infliximab treatment leads to downregulation of inflammation and tissue remodelling of the synovial membrane in spondyloarthropathy. Ann Rheum Dis. 2005;64:529-36.

53. Sica A, Mantovani A. Macrophage plasticity and polarization: in vivo veritas. J Clin Invest. 2012;122:787-95.

54. Ambarus CA, Noordenbos T, de Hair MHJ, Tak PP, Baeten DLP. Intimal lining layer macrophages but not synovial sublining macrophages display an IL-10 polarized-like phenotype in chronic synovitis. Arthritis Res Ther. 2012;14:R74.

55. Smith JA, Barnes MD, Hong D, DeLay ML, Inman RD, Colbert RA. Gene expression analysis of macrophages derived from ankylosing spondylitis patients reveals interferon-gamma dysregulation. Arthritis Rheum. 2008;58:1640-9.

56. Zhao J, Yuan W, Tao C, Sun P, Yang Z, Xu W. M2 polarization of monocytes in ankylosing spondylitis and relationship with inflammation and structural damage. APMIS Acta Pathol Microbiol Immunol Scand. 2017;125:1070-5.

57. Gracey E, Lin A, Akram A, Chiu B, Inman RD. Intracellular survival and persistence of Chlamydia muridarum is determined by macrophage polarization. PLoS One. 2013;8:e69421.

58. Baillet AC, Rehaume LM, Benham H, O'Meara CP, Armitage CW, Ruscher R, et al. High Chlamydia burden promotes tumor necrosis factor-dependent reactive arthritis in SKG mice. Arthritis Rheumatol Hoboken NJ. 2015;67:1535-47.

59. Essandoh K, Li Y, Huo J, Fan G-C. MiRNA-mediated macrophage polarization and its potential role in the regulation of inflammatory response. Shock. 2016;46:122-31.

60. Chen Q, Wang H, Liu Y, Song Y, Lai L, Han Q, et al. Inducible microRNA-223 down-regulation promotes TLR-triggered IL-6 and IL-1β production in macrophages by targeting STAT3. PLoS One. 2012;7:e42971.

61. Veremeyko T, Siddiqui S, Sotnikov I, Yung A, Ponomarev ED. IL-4/IL-13-dependent and independent expression of miR-124 and its contribution to M2 phenotype of monocytic cells in normal conditions and during allergic inflammation. PLoS One. 2013;8:e81774.

62. Wang H, Chao K, Ng SC, Bai AH, Yu Q, Yu J, et al. Pro-inflammatory miR-223 mediates the cross-talk between the IL23 pathway and the intestinal barrier in inflammatory bowel disease. Genome Biol. 2016;17:58.
63. Polytarchou C, Oikonomopoulos A, Mahurkar S, Touroutoglou A, Koukos G, Hommes DW, et al. Assessment of circulating microRNAs for the diagnosis and disease activity evaluation in patients with ulcerative colitis by using the nanostring technology. Inflamm Bowel Dis. 2015;21:2533-9.

64. Neudecker V, Haneklaus M, Jensen O, Khailova L, Masterson JC, Tye H, et al. Myeloid-derived miR-223 regulates intestinal inflammation via repression of the NLRP3 inflammasome. J Exp Med. 2017;214:1737-52.

65. Pachathundikandi SK, Backert S. Helicobacter pylori controls NLRP3 expression by regulating hsa-miR-223-3p and IL-10 in cultured and primary human immune cells. Innate Immun. 2018;24:11-23.

66. Hou J, Wang P, Lin L, Liu X, Ma F, An H, et al. MicroRNA-146a feedback inhibits RIG-I-dependent Type I IFN production in macrophages by targeting TRAF6, IRAK1, and IRAK2. J Immunol Baltim Md 1950. 2009; 183:2150-8.

67. Maitra U, Davis S, Reilly CM, Li L. Differential regulation of Foxp3 and IL-17 expression in CD4 T helper cells by IRAK-1. J Immunol Baltim Md 1950. 2009;182:5763-9.

68. Srivastava A, Nikamo P, Lohcharoenkal W, Li D, Meisgen F, Xu Landén N, et al. MicroRNA-146a suppresses IL-17-mediated skin inflammation and is genetically associated with psoriasis. J Allergy Clin Immunol. 2017;139:550-61.

69. Xia P, Fang X, Zhang Z, Huang Q, Yan K, Kang K, et al. Dysregulation of miRNA146a versus IRAK1 induces IL-17 persistence in the psoriatic skin lesions. Immunol Lett. 2012;148:151-62.

70. Wei C, Zhang H, Wei C, Mao Y. Correlation of the expression of miR-146a in peripheral blood mononuclear cells of patients with ankylosing spondylitis and inflammatory factors. Exp Ther Med. 2017;14:5027-31.

71. Di G, Kong L, Zhao Q, Ding T. MicroRNA-146a knockdown suppresses the progression of ankylosing spondylitis by targeting dickkopf 1. Biomed Pharmacother Biomedecine Pharmacother. 2018;97:1243-9.

72. Guo Z, Wu R, Gong J, Zhu W, Li Y, Wang Z, et al. Altered microRNA expression in inflamed and non-inflamed terminal ileal mucosa of adult patients with active Crohn's disease. J Gastroenterol Hepatol. 2015;30:109-16.

73. Magrey MN, Haqqi T, Haseeb A. Identification of plasma microRNA expression profile in radiographic axial spondyloarthritis-a pilot study. Clin Rheumatol. 2016;35:1323-7.

74. Dunaeva M, Blom J, Thurlings R, Pruijn GJM. Circulating serum miR-223-3p and miR-16-5p as possible biomarkers of early rheumatoid arthritis. Clin Exp Immunol. 2018;193(3):376-85.

75. Kim H-Y, Kwon HY, Ha Thi HT, Lee HJ, Kim GI, Hahm K-B, et al. MicroRNA-132 and microRNA-223 control positive feedback circuit by regulating FOXO3a in inflammatory bowel disease. J Gastroenterol Hepatol. 2016;31:1727-35.

76. Prajzlerová K, Grobelná K, Hušáková M, Forejtová Š, Jüngel A, Gay S, et al. Association between circulating miRNAs and spinal involvement in patients with axial spondyloarthritis. PLoS One. 2017;12:e0185323.

77. Akkaya-Ulum YZ, Balci-Peynircioglu B, Karadag O, Eroglu FK, Kalyoncu U, Kiraz S, et al. Alteration of the microRNA expression profile in familial Mediterranean fever patients. Clin Exp Rheumatol. 2017;35(Suppl 108):90-4.

78. Dweep H, Sticht C, Pandey P, Gretz N. miRWalk--database: prediction of possible miRNA binding sites by "walking" the genes of three genomes. J Biomed Inform. 2011;44:839-47.

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_8.png" /></figure>