Exploring the Extracellular Vesicle MicroRNA Expression Repertoire in Patients with Rheumatoid Arthritis and Ankylosing Spondylitis Treated with TNF Inhibitors

Wielińska, Joanna; Crossland, Rachel E.; Łacina, Piotr; Świerkot, Jerzy; Bugaj, Bartosz; Dickinson, Anne M.; Bogunia-Kubik, Katarzyna

doi:https://doi.org/10.1155/2021/2924935

Disease Markers

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Circulating Non-Coding RNAs as Promising Candidates for Biomarker-Guided Therapeutics

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 2924935 | https://doi.org/10.1155/2021/2924935

Exploring the Extracellular Vesicle MicroRNA Expression Repertoire in Patients with Rheumatoid Arthritis and Ankylosing Spondylitis Treated with TNF Inhibitors

Joanna Wielińska,¹Rachel E. Crossland,²Piotr Łacina,¹Jerzy Świerkot,³Bartosz Bugaj,³Anne M. Dickinson,²and Katarzyna Bogunia-Kubik¹

Academic Editor: Anum Naz

Received12 Jul 2021

Accepted03 Sept 2021

Published30 Sept 2021

Abstract

Rheumatoid arthritis (RA) and ankylosing spondylitis (AS) belong to the most common inflammatory rheumatic diseases. MicroRNAs (miRNAs) are small 18–22 RNA molecules that function as posttranscriptional regulators. They are abundantly present within extracellular vesicles (EVs), small intercellular communication vesicles that can be found in bodily fluids and that have key functions in pathological and physiological pathways. Recently, EVs have gained much interest because of their diagnostic and therapeutic potential. Using NanoString profiling technology, the miRNA repertoire of serum EVs was determined and compared in RA and AS patients before and after anti-TNF therapy to assess its potential use as a diagnostic and prognostic biomarker. Furthermore, possible functional effects of those miRNAs that were characterized by the most significant expression changes were evaluated using in silico prediction algorithms. The analysis revealed a unique profile of differentially expressed miRNAs in RA and AS patient serum EVs. We identified 12 miRNAs whose expression profiles enabled differentiation between RA and AS patients before induction of anti-TNF treatment, as well as 4 and 14 miRNAs whose repertoires were significantly changed during the treatment in RA and AS patients, respectively. In conclusion, our findings suggest that extracellular vesicle miRNAs could be used as potential biomarkers associated with RA and AS response to biological treatment.

1. Introduction

MicroRNAs (miRNAs, miRs) are a family of single-stranded, noncoding endogenous regulatory RNAs derived from double-stranded precursors, typically composed of 21-23 nucleotides. They are involved in the regulation of gene expression at the posttranscriptional level [1, 2]. They function by the way of complementary binding to the 3UTR regions of target mRNA via the RNA-induced silencing complex (RISC), resulting in inhibition of the translation process [3]. One miRNA molecule can regulate the expression of several genes (transcripts), and a transcript may have a 3-UTR region that is recognized by many miRNAs [4, 5]. It is estimated that encoding miRNA genes constitute 1-5% of all the genes in humans and animals [6]. Regulatory miRNAs and exosomes are becoming increasingly important in identification of molecular markers related to pathogenesis and prognosis of disease [7].

Extracellular vesicles (EVs) are lipid membrane-enclosed vesicles, released by cells into the extracellular space. EVs are a heterogeneous collection of exosomes, microvesicles, and apoptotic bodies, ranging in size from 40 nm to 4000 nm. Classification is based on their size, origin, and biological function. EVs are present in blood, saliva, urine, milk, and amniotic fluid and are secreted by all mammalian cell types [8]. Their main function is the transport of lipids, proteins, miRNAs, and mRNAs. Extracellular vesicles are also mediators in intercellular communication and immune-regulatory processes such as bone remodelling, which implicates them in pathogenesis of rheumatic diseases [9]. Although the unique mechanism of immune complex formation remains unclear, it has been observed that synovial exosomes contain citrullinated peptides, which are well-known autoantigens in RA [10]. It has been reported that EVs widely participate in RA development, including antigen presentation and immune complex formation, inflammation, delivery of miRNA, and destruction of extracellular matrix [11]. Moreover, EVs are considered to be promising biomarkers in joint diseases such as RA [12]. Recent studies have indicated a significantly higher level of plasma EVs in RA patients compared to healthy individuals [13]. Other studies established a role of EVs as biomarkers in arthritis by showing an association of exosomal amyloid A and lymphatic vessel endothelial hyaluronic acid receptor-1 with disease activity in RA [14] and identified higher expression of long noncoding RNA, HOX Transcript Antisense RNA (HOTAIR), in serum exosomes of RA cases [15].

Pathogenesis of rheumatoid arthritis (RA) and ankylosing spondylitis (AS) is considered to be multifactorial, with disease susceptibility being associated with genetic, environmental, and stochastic factors. RA and AS are the most common inflammatory rheumatic diseases [16, 17]. They are characterized by different clinical, laboratory, and imaging hallmarks [18, 19]. RA is an autoimmune disease characterized by symmetric and erosive arthritis typically affecting small- and medium-sized joints. It is marked by inflammation of joint synovial tissue, resulting in progression of cartilage and bone tissue damage, ultimately leading to disability [20]. AS is a systemic inflammatory disorder that affects the sacroiliac joints and the spine and can also affect peripheral joints, causing characteristic inflammatory back pain, which can lead to structural and functional impairments [21].

A number of studies have reported that alteration of miRNA profiles may play an important role in the pathogenesis of rheumatic diseases [22, 23], and thus, these profiles may constitute potential biomarkers [24]. Increased expression of miRNAs has been detected in various cell types of RA patients, and miR-146a was shown to mediate negative feedback in the immune response of RA [25]. This specific miRNA was found to be overexpressed in synovial fibroblasts, synovial tissue, synovial fluid monocytes, peripheral blood mononuclear cells, and serum plasma of RA patients [26]. miR-146a regulates gene expression of TNF receptor-associated factor 6 (TRAF6) and IL-1 receptor-associated kinase 1 (IRAK1) in inflammation and participates in a negative feedback loop [27]. Moreover, single-nucleotide polymorphisms (SNPs) in miR-146a may alter its expression. The rs2910164 SNP was found to reduce the expression level of miR-146a, which led to less efficient inhibition of target genes, including two molecules important for RA development, TRAF6 and IRAK1, suggesting that miR-146a rs2910164 could contribute to RA development [28]. Moreover, our previous study in a group of Polish patients with RA showed an association between rs2910164 C variant and higher expression of miRNA-146a in serum after three months of therapy with TNF inhibitors [29].

Despite the fact that anti-TNF treatment constitutes a breakthrough in management of RA and other rheumatic diseases, approximately 30% of patients do not achieve any improvement. Because patients with RA or AS have a variable response to treatment, identification of biomarkers capable of predicting therapeutic response is imperative.

Both RA and AS would benefit from discovery of biomarkers that could be detected when disease is present, distinguish between the two disorders, that are associated with disease progression and outcome, and help to predict the response to treatment.

This study is aimed at investigating whether analysis and comparison of miRNA profiles in patients with RA and AS could be used (i) for detection of diagnostic miRNA markers to distinguish RA from AS and (ii) before and after TNF-α inhibitor treatment to predict the outcome and effectiveness of this biological therapy on modulation of proinflammatory response.

2. Materials and Methods

2.1. Patients

AS patients were classified according to the 1984 modified New York criteria [30]. Adult Caucasians () included in the study were characterized with high-disease activity (defined as and ) before initiation of anti-TNF therapy and failed to respond to at least two nonsteroidal anti-inflammatory drugs (NSAIDs) for at least four weeks at maximum doses (if there were no contraindications). Effectiveness of current drug therapy was assessed using Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), which is based on a 0-10 scale measuring discomfort, pain, and fatigue.

RA patients enrolled in this project were classified according to the 2010 American College of Rheumatology (ACR)/European League Against Rheumatism (EULAR) criteria, as well as the presence of active disease represented by the disease activity score in 28 joints (). Patients qualified for the study were also before biologic agents therapy initiation and failed to respond to at least two disease-modifying antirheumatic drugs (DMARDs). They were ≥18 years, of Caucasian origin, and with complete medical records.

Exclusion criteria for both RA and AS patients comprised clinically significant impairment of hepatic and renal function, coexistence of connective tissue diseases, infections with hepatotropic viruses, infections resistant to therapy, an ongoing history of cancer or uncontrolled diabetes, alcohol abuse, pregnancy or breastfeeding, and insufficient clinical records.

The ASAS/EULAR criteria were used to assess the clinical outcome of anti-TNF treatment.

All RA patients responded to therapy. Good response was a reduction of the DAS28 score to a posttreatment value of , and moderate response was interpreted as and a posttreatment or and [31].

For AS, a good response was determined as a reduction of from baseline and at the endpoint. A moderate response was defined as a reduction of from baseline and at the endpoint [32].

The study was approved by the Wrocław Medical University Ethics Committee (identification code KB-625/2016), and written informed consent was obtained from all participants.

Clinical characteristics of the patients are presented in Table 1 as mean and standard deviation (± SD).

2.2. Sample Preparation

Serum samples were collected from three patients (two men and one woman, aged: 23-46, mean age: 35) with RA and three patients (two men and one woman, aged: 36-70, mean age: 48) with AS at two-time points before and after three months of anti-TNF treatment initialization. All RA patients were treated with Etanercept. In the AS group, two patients were treated with Etanercept and one with Adalimumab. Further characteristics of patients enrolled in this study are detailed in Table 1. Sera were prepared from whole blood samples, which had been collected in 8.5 ml BD Vacutainer® SST™ II advance blood tubes (Becton, Dickinson and Company, Franklin Lakes, NJ, USA) with clot activator and acrylic gel separation. Samples were centrifuged at 1500 g for 10 min and stored at −70°C.

2.3. EV and RNA Isolation

Extracellular vesicles (EVs) were isolated from 3 ml of serum using Total Exosome Isolation Reagent (from Serum) (ThermoFisher Scientific), following supplier’s guidelines. EV pellets were stored at -80°C.

The EV pellets were resuspended in particle-free PBS, and RNA was isolated using Total Exosome RNA and Protein Isolation Kit (ThermoFisher Scientific), following supplier’s protocol. RNA was concentrated incorporating Amicon Ultra-0.5 centrifugal filter unit with Ultracel-3 membrane (Merck Millipore), according to NanoString recommendations.

RNA concentration was determined using the 2100 Bioanalyzer and the RNA 6000 Pico Kit (Agilent Technologies). While accurate RIN numbers could not be generated, due to the lack of ribosomal RNA in EVs, RNA concentration and peak area were assessed. All RNA was stored at −80°C.

2.4. EV Morphology and Size Measurement

To find out whether EVs were properly isolated, EV morphology and size distribution were assessed. For EV morphology, transmission electron microscopy (TEM) was used. 300-mesh grids were filmed with Pioloform® resin (SPI Supplies), carbon coated, and plasma etched before use. EVs were directly applied on the grid, stained with 10 μl of 2% aqueous uranyl, and air-dried. Examination was conducted on a Hitachi HT7800 transmission electron microscope equipped with an Emsis Xarosa camera with Radius software, in cooperation with the Electron Microscopy Research Services, Newcastle University.

For EV size distribution analysis, nanoparticle tracking analysis (NTA) was performed employing a NanoSight LM10 microscope supplied with NTA software version 3.0 (NanoSight Ltd., UK). Background extraction with blur settings and maximum jump distance was applied automatically, and second recordings were taken for each sample.

2.5. miRNA Profiling

Serum EV microRNA expression profiling was performed using the nCounter® Human v3 miRNA Expression Assay kit (NanoString Technologies) as previously described [33]. This code set comprises 98% of microRNA sequences found in miRbase v22 and includes 798 mature microRNAs, six positive and eight negative controls, six ligation controls, and five reference controls. The procedure was performed according to manufacturer’s guidelines. Data normalization was performed using nSolver Analysis Software v4.0 (NanoString Technologies), with positive control normalization using the geometric mean and normalization flagging outside the normalization factor range 0.3-3.0. Codeset content normalization was performed using the top 100 microRNAs for normalization, based on geometric mean and flagging outside the normalization factor range 0.1-10.0.

2.6. Statistical Analysis

Data normalization and fold change (FC) expression differences between groups were conducted using nSolver v4.0 software (NanoString Technologies). Further analyses were performed using the R software (version 3.6.1) with RStudio 1.2.5001 (RStudio, Inc., USA) applying an analysis pipeline designed by Newcastle University, Haematological Sciences Department. “ggplot2” (v2.1.0) package function was used to construct Volcano plots. Heatmaps with unsupervised clustering were generated using “gplots” (v2.17.0) and “RColorBrewer” (v1.1-2). values between the two groups were calculated using a two-tailed -test. Significance was set at .

2.7. Target Prediction and Pathway Analysis

To perform gene prediction and pathway analyses for miRNAs obtained from NanoString, we incorporated an approach described by Lou et al. [34]. Potential microRNA gene targets were identified by the miRNet database (http://www.mirnet.ca/) [35]. The STRING database (http://string-db.org) was used to construct networks of protein-protein interactions based on target genes obtained from miRNet [36]. Hub genes, or genes from the protein-protein interaction network that have the highest degree of connectivity, were determined using the Cytoscape software (version 3.7.2) [37]. They were subsequently used as input in the KEGG pathway enrichment analysis performed in the Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.abcc.ncifc http://rf.gov/) [38].

3. Results

3.1. NanoString Experimental Setup and MicroRNA Expression Counts

In this study, serum samples from three RA patients and three AS patients before and after administration of biological treatment were analysed. Samples were labelled according to diagnosis (RA or AS), time points of sample collection: before anti-TNF treatment induction (indicated by odd numbers, for example: RA1 and RA3) and after three month of drug administration (indicated by respective even numbers, RA2 and RA4). This labelling method resulted in the following groups: six RA samples, including three before treatment (RA1, RA3, and RA5) and three after therapy (RA2, RA4, and RA6), and similarly six AS samples, three before treatment (AS1, AS3, and AS5) and additional three samples after three months after treatment initialization (AS2, AS4, and AS6).

All NanoString samples passed quality control parameters according to their microRNA expression profiles. A total of 159 microRNAs were positively expressed in >2 RA and AS samples and were included in the final analysis.

3.2. MicroRNA Expression Analysis in RA EV Samples before vs. after Therapy

Unsupervised hierarchical clustering analysis separated RA serum EV samples before and after three months of anti-TNF treatment as shown in Figure 1(a). EVs from serum collected after three months of anti-TNF therapy were characterized by a unique miRNA signature consisting of four miRNA molecules that were overexpressed (), compared to samples collected before treatment: miR-520 h (4.77-fold, ), miR-498 (3.28-fold, ), miR-548n (1.66-fold, ), and miR-19b-3p (1.35-fold, ) (Figures 1(b) and 1(c), Table 2).

Figure 1

Serum EV microRNA expression in RA and AS patients before vs. after anti-TNF treatment. (a, d) Unsupervised hierarchical clustering analysis. Samples before (1, 3, and 5) vs. after (2, 4, and 6) anti-TNF therapy. (b, e) Volcano plots showing the relationship between fold change and significance for RA and AS patients before and after anti-TNF therapy. The horizontal dashed line indicates cutoff for significance (−log10 value > 1.3) and the vertical lines for /≤−1.5. Significantly different miRNAs are highlighted in blue. (c) Heatmaps showing unsupervised hierarchical clustering of differentially expressed miRNAs in serum EVs of patients before () vs. after Etanercept treatment. The colour scale indicates relative fold change (red: high; blue: low).

3.3. MicroRNA Expression Analysis in AS EV Samples before vs. after Therapy

AS samples were similarly evaluated and unsupervised clustering analysis separated serum EV samples taken before and after three months of anti-TNF treatment (Figure 1(d)). Expression profiling analysis identified 14 miRNAs significantly () differentially expressed in AS after therapy, of which eight were overexpressed: miR-130a-3p, miR-146a-5p, miR-21-5p, miR-22-3p, miR-23a-3p, miR-30a-5p, miR-362-3p, miR-548ah-5p (with FC range 1.49-2.82, value range: 0.002-0.048), while 6 were downregulated: let-7c-5p, let-7f-5p, miR-125a-5p, miR-18a-5p, miR-374b-5p, and miR-98-5p (FC range: 1.48-10.99, value range: 0.023-0.044). The greatest fold change () was detected for hsa-miR-98-5p (Figures 1(e) and 1(f), Table 2).

3.4. Comparison of RA vs. AS EV miRNA Expression before Anti-TNF Therapy Implementation

To assess whether RA and AS patients can be distinguished based on their serum EV expression profiles, we compared miRNA repertoires at the baseline, before induction of biological treatment (Figure 2(a)). The comparison of RA and AS serum EVs before anti-TNF treatment implementation identified 12 miRNAs with significantly different expression (FC range: 1.08-11.58, value range: 0.001-0.048) (Figures 2(b) and 2(c)). Among those miRNA molecules, higher expression of miR-125a-5p, miR-130b-3p, miR-151a-5p, miR-301a-3p, and miR-324-5p characterized patients with RA (FC range: 4.31-11.58, value range: 0.012-0.031), while miR-376c-3p, miR-378h, miR-411-5p, miR-548a-5p, miR-548n, miR-548q, and miR-579-3p were overexpressed in AS patients, but with a smaller FC range (1.08-3.50, value range: 0.001-0.048) (Figures 2(b) and 2(c), Table 2).

Figure 2

Serum EV microRNA expression in RA vs. AS patients before and three months after anti-TNF treatment. (a, d) Unsupervised hierarchical clustering analysis. RA vs. AS patients before biological treatment initialization and after agent administration. (b, e) Volcano plots showing the relationship between fold change and significance for RA vs. AS patients before and after anti-TNF therapy. The horizontal dashed line indicates cutoff for significance (−log10 value > 1.3) and the vertical lines for /≤−1.5. Significantly different miRNAs are highlighted in blue. (c) Heatmaps showing unsupervised hierarchical clustering of differentially expressed miRNAs in serum EVs of RA vs. AS patients before () and after treatment. The colour scale indicates relative fold change (red: high; blue: low).

3.5. Comparison of RA vs. AS EV miRNA Expression after Therapy

A similar comparison of miRNA repertoires between patients with RA vs. AS was performed on the RA and AS samples collected three months after initialization of anti-TNF treatment. A total of 18 microRNAs were significantly differentially expressed after therapy in RA compared to AS patients (FC range: 1.77-12.76, value range: 0.004-0.050), of which eight were upregulated in AS (miR-7i-5p, miR-1915-3p, miR-30a-5p, miR-3158-3p, miR-379-5p, miR-496, miR-612, and miR-649; FC range: 1.77-3.08, value range: 0.004-0.050) (Figures 2(e) and 2(f), Table 2).

Furthermore, miR-30a was upregulated in AS after therapy in comparison to its expression before treatment and its expression in patients with RA after anti-TNF drug administration. On the other hand, miR-98 was downregulated after initialization of anti-TNF therapy in both cases. Interestingly, miR-151a and miR-125a were characterized with lower expression levels in AS patients compared to those with RA after, as well as before biological therapy.

3.6. Identification of Target Genes and Potential Pathways

The potential target genes were analysed separately for upregulated and downregulated microRNAs. 10 hub genes were identified for each of the analysed groups. In RA patients, all upregulated microRNAs after therapy were linked to TP53, EP300, PTEN, MAPK1, STAT3, ESR1, EZH2, CCNB1, BRCA1, and CASP3. In AS patients, both upregulated and downregulated microRNAs after treatment were associated with TP53, AKT1, MYC, UBC, EGFR, and IL6. Additionally, downregulated microRNAs after therapy were linked to UBA52, CCND1, PTEN, and STAT3, whereas the upregulated ones were associated with RPS27A, MAPK1, UBB, and VEGFA. In comparison between RA and AS patients before treatment, downregulated microRNAs were connected to TP53, AKT1, MYC, UBC, EGFR, RPS27A, MAPK1, UBB, HSPA8, and PTEN, while upregulated microRNAs were linked to HSP90AA1, CCNB1, STAT3, CDC5L, MDM2, CASP3, ATM, SKP1, ACTB, and IGF1R. The hub genes detected in comparison between RA and AS patients after therapy linked to downregulated, as well as upregulated microRNAs were TP53, AKT1, MYC, EGFR, and CTNNB1, although downregulated microRNAs were also related to UBC, PTEN, CCND1, HSPA8, and VEGFA and upregulated to MAPK1, JUN, HDAC1, NOTCH1, and MAPK8.

The KEGG pathway enrichment analysis showed that targets of miRNAs identified as upregulated/downregulated in the treatment-related analyses (RA before vs. after, AS before vs. after) were enriched for pathways associated with response to drugs and cellular response to drugs. Additionally, targets of miRNAs upregulated in AS vs. RA (after therapy) and in AS patients after treatment were enriched in pathways associated with the TNF signalling pathway.

4. Discussion

Anti-TNF treatment constitutes a breakthrough in management of rheumatic diseases, although many patients do not achieve significant or any improvement. This makes it essential to find new biomarkers capable of predicting therapeutic response. In the present study, we examined whether serum EV miRNA profiles could be used to distinguish rheumatoid arthritis from ankylosing spondylitis, as well as to predict outcome of TNF inhibitor treatment in both of these rheumatic diseases. We found that in RA patients, four miRNAs (miR-520h, miR-498, miR-548n, and miR-19b-3p) were differentially expressed after anti-TNF treatment, while fourteen were differentially expressed in AS patients (miR-130a-3p, miR-146a-5p, miR-21-5p, miR-22-3p, miR-23a-3p, miR-30a-5p, miR-362-3p, miR-548ah-5p, let-7c-5p, let-7f-5p, miR-125a-5p, miR-18a-5p, miR-374b-5p, and miR-98-5p). Additionally, we identified twelve miRNAs that distinguished the two diseases (miR-125a-5p, miR-130b-3p, miR-151a-5p, miR-301a-3p, miR-324-5p, miR-376c-3p, miR-378h, miR-411-5p, miR-548a-5p, miR-548n, miR-548q, and miR-579-3p).

Recently, there have been many studies investigating EV miRNAs as markers for evaluating RA progression. To date, the association between miRNAs and parameters of disease activity were reported in East Asian populations, specifically in Korean [39] and Chinese [40] RA patients. Wang et al. showed downregulation of exosome-delivered miR-548a-3p in serum of RA patients in contrast with healthy controls, as well as a lower level of miR-548a-3p correlating with higher levels of CRP, RF, and ESR [40]. Similarly, serum exosomal miR-6089 was significantly decreased in RA Chinese patients compared to controls [41]. Regarding plasma exosomes, secretion of miR-17, miR-19b, and miR-121 was significantly higher in RA patients compared to healthy individuals. Moreover, it was found that transport of miR-17 into T cells represses Treg induction and differentiation [42].

MiRNAs delivered by exosomes were also described in the context of their therapeutic potential in RA. miR-150-5p reduced joint destruction by inhibiting angiogenesis mediated by downregulation of matrix metalloproteinase 14 (MMP14) and vascular endothelial growth factor (VEGF) [43]. However, upregulation of exosomal miR-92 boosted bone destruction by blocking apoptosis of fibroblast-like synoviocytes (FLSs) and inflammatory cytokine release [44]. Another microRNA, miR-let-7b, and its ligation to TLR-7 were able to induce joint inflammation through M1 macrophage differentiation [45]. Furthermore, overexpression of miR-221-3p in exosomes isolated from inflamed FLSs might suppress bone formation at erosion sites [46]. These findings provide evidence that exosomal miRNAs can modulate inflammatory responses during RA pathogenesis.

To date, researchers have focused mostly on microRNAs extracted from whole blood, serum/plasma, peripheral blood mononuclear cells (PBMC), and rheumatoid arthritis synovial fibroblasts (RASF) to investigate their role in RA progression [23, 47, 48]. Herein, we present isolated serum exosomes as a novel source of microRNAs. Various techniques can be used to study the miRNA profile. The most common practice is to employ microarray profiling to analyse a wide range of microRNAs on a qRT-PCR instrument [49]. However, a novel technique, NanoString, described previously in a study on serum miRNAs in patients with graft-versus-host disease (GvHD) [33, 50], allows analysis of around 800 microRNAs using digital barcoding technology. This removes the need for reverse transcription or preamplification of RNA, as the technology is able to directly count isolated microRNA molecules. During our study, NanoString was used to assess the repertoire of serum EV microRNAs in RA and AS cases subjected to biological treatment with anti-TNF agents. In the present study, we established for the first time that miR-498, miR-520h, and miR-548n are present in serum EVs of RA patients. Furthermore, expression of miR-7c-5p, miR-7f-5p, miR-362-3p, miR-3746-5p, and miR-548ah-5p has been characterized by us for the first time in AS EVs; in fact, these miRNAs have never been described in any rheumatoid disease at all.

As described in the literature (summary presented in Table 3), miR-155 [51] and miR-146a [52] are among the most widely studied microRNAs in RA and AS patients. Although our analysis did not identify miR-155 expression as significantly different between RA and AS patients, our results regarding miR-146a are consistent with those previously described. Higher expression of miR-146a in RA is well-characterized in multiple components such as synovial fibroblasts [53], CD4+ T cells [54], serum [55], and peripheral blood cells [26]. Moreover, upregulation of miR-146a in various samples from RA patients can distinguish them from osteoarthritis (OA) patients [56]. Interestingly, miR-146a is an NF-κB-dependent gene and controls inflammatory responses through inhibition of IL-1 receptor-associated kinase 1 and TNF receptor-associated factor 6 proteins [25] and downregulation of TLR4/NF-κB pathway [57]. By studying exosomes from rheumatoid arthritis synovial fibroblast cell line, Takamura et al. observed upregulation of miR-146a caused by TNF-α stimulation [58]. A positive correlation between the miR-146a level and IL-1β, IL-6, and TNF-α expression was also observed in AS [59].

There is some evidence to suggest a role for miR-146a in the course of the disease. Filková et al. characterized a decreased level of circulating miR-146a at an early stage of RA, compared to established disease [60]. Furthermore, the expression of miR-146a correlated positively with disease activity in RA [61, 62]. A baseline level of miR-146a was lower in methotrexate responders than in nonresponders [63]. Regarding anti-TNF therapy, previous studies reported increased levels of serum miR-146a in patients under treatment [29, 64]. Besides, Liu et al. demonstrated the potential predictive value of miR-146a measurement for biological agent therapy outcome in RA patients [65].

Recent studies also highlight a role of miR-146a in AS. Significantly higher expression levels have been found in patient serum samples compared to healthy controls [66]. miR-146a overexpression can also cause inhibition of fibroblast proliferation and osteogenic potential while its knockdown blocked disease progression by regulating Dickkopf Wnt Signalling Pathway Inhibitor 1 (DKK1) expression [67]. Other results suggested a positive correlation between miR-146a expression in peripheral blood mononuclear cells and duration of morning stiffness, ESR, CRP, and BASDAI [68]. Prajzlerová et al. also established the role of miR-146a in the pathogenesis of Axial Spondyloarthritis (AxSpA) [69]. Our current study does not show any significant differences in miR-146a expression between RA patients in two time points (before and after three months of the treatment) or RA group compared to AS. However, we detected upregulation of miR-146a-5p in AS cases after anti-TNF therapy. This observation is consistent with previous results obtained by our group for RA patients [29] where the expression level of miR-146a was lower in serum of patients before therapy compared to that after three months of treatment. Taken together, those conclusions implicate miR-146a as a potential noninvasive biomarker in RA and AS prognosis.

Another microRNA, miR-125b, can also serve as a predictor of anti-TNF therapy response [70]. We identified downregulation of miR-125a-5p in AS patient serum EVs after therapy. In contrast, Castro-Villegas et al. observed that patients responding to anti-TNF/DMARD combination therapy exhibited overexpression of miR125b-5p after treatment and lower RF, CRP, TNFα, IL-17, and IL-6 levels [64]. Previous studies confirmed that the upregulation of miR-125a and miR-125b directly affects NF-κB [71, 72]. Likewise, decreased miR-125a-5p expression enhanced the protein level of TNFRSF1B and reduced osteoclast activity [73].

On the other hand, Perez-Sanchez et al. found that expression levels of miR-146a-5p, miR125a-5p, miR-151-3p, miR-22-3p, and miR-451a were higher in AS patients compared to psoriatic arthritis patients. Furthermore, miR-146a-5p, miR-125a-5p, and miR-22-3p can distinguish the active and nonactive stage of the disease. Expression of miR-125a-5p, miR-151a-3p, miR-150-5p, and miR-451a was also related to the presence of syndesmophytes in AS patients [74]. Although we demonstrated that miR-146a-5p and miR-22-3p were upregulated in AS patient EVs after therapy, miR-125a-5p was downregulated. Moreover, miR-125a-5p was characterized by lower expression in AS patient EVs compared to RA EVs before and after treatment. Similarly, miR-151-5p was decreased in AS EV samples before and after anti-TNF treatment compared to RA samples. These results implicate miR-125a-5p as a potential diagnostic biomarker in AS.

Little is known about miR-130a and miR-301a-3p. In AS patients after biological therapy, we observed an upregulation of miR-130a. A molecular mechanism involving histone deacetylase 3 (HDAC3) was published by Jiang and Wang [75]. A study on miR-130a in human chondrocytes identified its crucial role in regulating TNF-α expression [76]. miR-301a-3p was found to be overexpressed in the PBMCs and associated with Th17 cell frequency in RA patients [77]. We hypothesize that miR-130b-3p and miR301a-3p EV expression may distinguish RA patients from AS patients based on our observation that those microRNAs are higher in RA cases.

To date, publications describing miR-18a referred to osteoarthritis [78, 79], primary Sjögren’s syndrome [80], and rheumatoid arthritis [81]. Herein, we demonstrate for the first time downregulation of miR-18a-5p in AS patient EVs after anti-TNF therapy.

We also observed that EVs isolated from sera of RA patients showed overexpression of miR-19b-3p, miR-498, miR-520h, and miR-548n. miR-19b-3p was previously investigated in knee osteoarthritis patients and associated with disease severity [82]. Duan et al. revealed miR-19b-3p involvement in OA through the GRK6-NF-κB pathway [83]. Gantier et al. suggested that miR-19 regulates NF-κB signalling [84]; however, data collected from an experiment with RA FLSs underline the role of miR-19a/b in the stimulation of TLR2 expression [85].

To the best of our knowledge, we showed for the first time that miR-548n and miR-548ah-5p are overexpressed after therapy in RA and AS patients, respectively. Moreover, we identified miR-548n, miR-548ah-5p, and miR-548q as possible biomarkers of AS. Even though no results describing miR-548n, miR-548ah-5p, or miR-548q in rheumatic diseases have been published before, our results are consistent with those reporting about miR-548a-3p, which belongs to the hsa-miR-548 family. miR-548a-3p has previously been shown to be downregulated in serum EVs of RA compared to healthy individuals. Furthermore, expression was negatively correlated with ESR, CRP, and RF levels in patients [40]. In our study, expression of miR-548a-5p was also decreased in RA patients compared to AS patients, before initialization of anti-TNF therapy.

Several studies reported the differential expression of miR-29a and miR-21, making them potentially involved in AS pathogenesis [86–89]. Additionally, overexpression of miR-29a was described in PBMCs from AS patients after etanercept treatment [90]. However, we did not observe differential expression of miR-29a in either RA or AS samples. The analysis showed a significant increase in expression of miR-21 in AS patients after anti-TNF treatment. In contrast, Huang et al. reported greater miR-21 expression in whole blood of patients with AS compared to controls [91]. Recently, Zou et al. showed that upregulation of miR-21 was associated with radiographic severity of AS [92]. The studies regarding RA reported decreased miR-21 level in PBMCs, CD4+ T cells [93], and plasma [94]. Balzano et al. hypothesized that low expression of miR-21-5p was a result of corticosteroids that inhibit NF-κB [94]. Further analyses revealed that miR-21 may be implicated in several signalling pathways such as the IL-34/STAT3/miR-21 pathway, essential for synovial fibroblast survival in RA [95], as well as a mediator between inflammation and bone formation through the JAK2/STAT3 pathway in AS [96]. Additionally, miR-21 plays a role in mediation of RANKL-induced osteoclastogenesis and downregulation of programmed cell death 4 (PDCD4) protein levels [97].

Other microRNA considered to be involved in AS are miR-30a and let-7i. MiR-30a was downregulated in serum/plasma of radiographic axial spondyloarthritis patients compared to healthy individuals [98]. In turn, let-7i was upregulated in T cells [99] and plasma [100] of AS patients compared to controls. In our analysis, both miR-30a-5p and let-7i-5p levels were higher in AS patient EVs after therapy compared to that before treatment and RA cases, respectively. In the study conducted by Lai et al., expression of let-7i in T cells was also positively correlated with the Bath Ankylosing Spondylitis Radiology Index (BASRI) of the lumbar spine, a scoring system which we did not include in our analysis. The authors concluded that overexpression of these microRNAs suppressed TLR-4 expression, which led to downregulation of TLR-4 [99].

The present work is in line with the previous finding identifying miR-23a as a potential predictor of etanercept response [101]. miR-23b can repress IL-17-associated autoimmune inflammation in human fibroblast-like synoviocytes [102]. This observation was confirmed for miR-23a in articular cartilage tissues from RA patients [103].

In summary, to the best of our knowledge, the present study identified for the first time 12 microRNAs differently expressed in serum EVs between RA and AS patients before biological agent administration. miR-125a-5p, miR-130b-3p, miR-151a-5p, miR-301a-3p, and miR-324-5p were upregulated in RA EVs, and miR-376c-3p, miR-378h, miR-411-5p, miR-548a-5p, miR-548n, miR-548q, and miR-579-3p in AS EVs. We believe that these microRNAs have the potential to distinguish RA pathogenesis from AS. The pathway prediction analysis performed using mirPATH v.3 DIANA tools (data not shown) found that all the microRNAs are involved in the Wnt signalling pathway.

5. Conclusions

Our analysis has revealed a unique profile of differentially expressed miRNAs in RA and AS patient serum EVs, both before and three months after anti-TNF treatment administration. These results suggest that EV miRNA profiling of RA and AS patients can be used for detection of diagnostic and predictive biomarkers. They confirm a potential role of these miRNAs in distinguishing the two diseases and in the prediction of response to treatment in RA and AS patients. Nevertheless, the results reported herein should be considered as a pilot study that was conducted on a limited number of patients; therefore, validation in larger verification cohorts is required.

Data Availability

All relevant data are within the manuscript. More detailed information cannot be shared publicly because it allows identification of the patients. It can be accessed from the Hirszfeld Institute of Immunology and Experimental Therapy (contact via [email protected]) for researchers who meet the criteria for access to confidential data.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

This research was funded by the National Science Centre (Poland) (grant number 2016/21/B/NZ5/01901).

References

W. Filipowicz, S. N. Bhattacharyya, and N. Sonenberg, “Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?” Nature Reviews. Genetics, vol. 9, no. 2, pp. 102–114, 2008.
View at: Publisher Site | Google Scholar
J. I. Satoh and H. Tabunoki, “Comprehensive analysis of human microRNA target networks,” BioData Mining, vol. 4, 2011.
View at: Publisher Site | Google Scholar
D. P. Bartel, “MicroRNAs: genomics, biogenesis, mechanism, and function,” Cell, vol. 116, pp. 281–297, 2004.
View at: Publisher Site | Google Scholar
L. F. Gulyaeva and N. E. Kushlinskiy, “Regulatory mechanisms of microRNA expression,” Journal of Translational Medicine, vol. 14, p. 143, 2016.
View at: Publisher Site | Google Scholar
J. V. Fritz, A. Heintz-Buschart, A. Ghosal et al., “Sources and functions of extracellular small RNAs in human circulation,” Annual Review of Nutrition, vol. 36, pp. 301–336, 2016.
View at: Publisher Site | Google Scholar
H. Liang, L. Huang, J. Cao, K. Zen, X. Chen, and C.-Y. Zhang, “Regulation of mammalian gene expression by exogenous microRNAs,” WIREs RNA, vol. 2012, pp. 733–742, 2012.
View at: Publisher Site | Google Scholar
Z. Zakeri, A. Salmaninejad, N. Hosseini et al., “MicroRNA and exosome: key players in rheumatoid arthritis,” Journal of Cellular Biochemistry, vol. 120, no. 7, pp. 10930–10944, 2019.
View at: Publisher Site | Google Scholar
R. Pourakbari, M. Khodadadi, A. Aghebati-Maleki, L. Aghebati-Maleki, and M. Yousefi, “The potential of exosomes in the therapy of the cartilage and bone complications; emphasis on osteoarthritis,” Life Sciences, vol. 236, p. 116861, 2019.
View at: Publisher Site | Google Scholar
J. Behera and N. Tyagi, “Exosomes: mediators of bone diseases, protection, and therapeutics potential,” Oncoscience, vol. 5, pp. 181–195, 2018.
View at: Publisher Site | Google Scholar
K. Skriner, K. Adolph, P. R. Jungblut, and G. R. Burmester, “Association of citrullinated proteins with synovial exosomes,” Arthritis and Rheumatism, vol. 54, pp. 3809–3814, 2006.
View at: Publisher Site | Google Scholar
J. Withrow, C. Murphy, Y. Liu, M. Hunter, S. Fulzele, and M. W. Hamrick, “Extracellular vesicles in the pathogenesis of rheumatoid arthritis and osteoarthritis,” Arthritis Research & Therapy, vol. 18, p. 286, 2016.
View at: Publisher Site | Google Scholar
J. Malda, J. Boere, C. H. A. Van De Lest, P. R. Van Weeren, and M. H. M. Wauben, “Extracellular vesicles - new tool for joint repair and regeneration,” Nature Reviews Rheumatology, vol. 12, pp. 243–249, 2016.
View at: Publisher Site | Google Scholar
J. Sellam, V. Proulle, A. Jüngel et al., “Increased levels of circulating microparticles in primary Sjögren’s syndrome, systemic lupus erythematosus and rheumatoid arthritis and relation with disease activity,” Arthritis Research & Therapy, vol. 11, p. R156, 2009.
View at: Publisher Site | Google Scholar
J. Yoo, S. K. Lee, M. Lim, D. Sheen, E. H. Choi, and S. A. Kim, “Exosomal amyloid A and lymphatic vessel endothelial hyaluronic acid receptor-1 proteins are associated with disease activity in rheumatoid arthritis,” Arthritis Research & Therapy, vol. 19, 2017.
View at: Publisher Site | Google Scholar
J. Song, D. Kim, J. Han, Y. Kim, M. Lee, and E. J. Jin, “PBMC and exosome-derived Hotair is a critical regulator and potent marker for rheumatoid arthritis,” Clinical and Experimental Medicine, vol. 15, pp. 121–126, 2014.
View at: Publisher Site | Google Scholar
J. Braun and J. Sieper, “Ankylosing spondylitis,” Lancet, vol. 369, pp. 1379–1390, 2007.
View at: Publisher Site | Google Scholar
D. L. Scott, F. Wolfe, and T. W. J. Huizinga, “Rheumatoid arthritis,” Lancet, vol. 376, pp. 1094–1108, 2010.
View at: Publisher Site | Google Scholar
B. Michelsen, R. Fiane, A. P. Diamantopoulos et al., “A comparison of disease burden in rheumatoid arthritis, psoriatic arthritis and axial spondyloarthritis,” PLoS One, vol. 10, article e0123582, 2015.
View at: Publisher Site | Google Scholar
G. Schett, L. C. Coates, Z. R. Ash, S. Finzel, and P. G. Conaghan, “Structural damage in rheumatoid arthritis, psoriatic arthritis, and ankylosing spondylitis: traditional views, novel insights gained from TNF blockade, and concepts for the future,” Arthritis Research & Therapy, vol. 13, p. S4, 2011.
View at: Publisher Site | Google Scholar
S. Jung, K. Kim, C. Yang, S. Park, and J. Ju, “Cytokine-mediated bone destruction in rheumatoid arthritis,” Journal of Immunology Research, vol. 2014, 2014.
View at: Publisher Site | Google Scholar
M. C. Turina, J. J. de Winter, J. E. Paramarta et al., “Clinical and imaging signs of spondyloarthritis in first-degree relatives of HLA–B27–positive ankylosing spondylitis patients: the pre-spondyloarthritis (pre-SpA) cohort study,” Arthritis & Rhematology, vol. 68, pp. 2444–2455, 2016.
View at: Publisher Site | Google Scholar
I. Alevizos and G. G. Illei, “MicroRNAs as biomarkers in rheumatic diseases,” Nature Reviews Rheumatology, vol. 6, pp. 391–398, 2010.
View at: Publisher Site | Google Scholar
G. Evangelatos, G. E. Fragoulis, V. Koulouri, and G. I. Lambrou, “MicroRNAs in rheumatoid arthritis: from pathogenesis to clinical impact,” Autoimmunity Reviews, vol. 18, p. 102391, 2019.
View at: Publisher Site | Google Scholar
J. Wielinska and K. Bogunia-Kubik, “miRNAs as potential biomarkers of treatment outcome in rheumatoid arthritis and ankylosing spondylitis,” Pharmacogenomics, vol. 22, pp. 291–301, 2021.
View at: Publisher Site | Google Scholar
K. D. Taganov, M. P. Boldin, K.-J. Chang, and D. Baltimore, “NF-B-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses,” PNAS, vol. 103, pp. 12481–12486, 2006.
View at: Publisher Site | Google Scholar
K. M. Pauley, M. Satoh, A. L. Chan, M. R. Bubb, W. H. Reeves, and E. K. L. Chan, “Upregulated miR-146a expression in peripheral blood mononuclear cells from rheumatoid arthritis patients,” Arthritis Research & Therapy, vol. 10, pp. 1–10, 2008.
View at: Publisher Site | Google Scholar
A. Chatzikyriakidou, P. V. Voulgari, I. Georgiou, and A. A. Drosos, “The role of microRNA-146a (miR-146a) and its target IL-1R-associated kinase (IRAK1) in psoriatic arthritis susceptibility,” Scandinavian Journal of Immunology, vol. 71, pp. 382–385, 2010.
View at: Publisher Site | Google Scholar
K. Jazdzewski, E. L. Murray, K. Franssila, B. Jarzab, D. R. Schoenberg, and A. De La Chapelle, “Common SNP in pre-miR-146a decreases mature miR expression and predisposes to papillary thyroid carcinoma,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, pp. 7269–7274, 2008.
View at: Publisher Site | Google Scholar
K. Bogunia-Kubik, B. Wysoczańska, D. Piątek, M. Iwaszko, M. Ciechomska, and J. Świerkot, “Significance of polymorphism and expression of miR-146a and NFkB1 genetic variants in patients with rheumatoid arthritis,” Archivum Immunologiae et Therapiae Experimentalis (Warsz), vol. 64, pp. 131–136, 2016.
View at: Publisher Site | Google Scholar
S. van der Linden, H. A. Valkenburg, and A. Cats, “Evaluation of diagnostic criteria for ankylosing spondylitis. A proposal for modification of the New York criteria,” Arthritis and Rheumatism, vol. 27, pp. 361–368, 1984.
View at: Publisher Site | Google Scholar
J. Fransen and P. L. C. M. van Riel, “The disease activity score and the EULAR response criteria,” Clinical and Experimental Rheumatology, vol. 23, pp. S93–S99, 2005.
View at: Publisher Site | Google Scholar
J. Zochling and J. Braun, “Assessment of ankylosing spondylitis,” Clinical and Experimental Rheumatology, vol. 23, 5 Suppl 39, pp. S133–S141, 2005.
View at: Google Scholar
R. E. Crossland, J. Norden, M. K. Juric et al., “Expression of serum microRNAs is altered during acute graft-versus-host disease,” Frontiers in Immunology, vol. 8, pp. 1–14, 2017.
View at: Publisher Site | Google Scholar
W. Lou, J. Liu, B. Ding et al., “Identification of potential miRNA–mRNA regulatory network contributing to pathogenesis of HBV-related HCC,” Journal of Translational Medicine, vol. 17, 2019.
View at: Publisher Site | Google Scholar
Y. Fan, K. Siklenka, S. K. Arora, P. Ribeiro, S. Kimmins, and J. Xia, “miRNet-dissecting miRNA-target interactions and functional associations through network-based visual analysis,” Nucleic Acids Research, vol. 44, pp. 135–141, 2016.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. Franceschini, S. Wyder et al., “STRING v10: protein–protein interaction networks, integrated over the tree of life,” Nucleic Acids Research, vol. 43, pp. D447–D452, 2015.
View at: Publisher Site | Google Scholar
P. Shannon, “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
W. Huang Da, B. Sherman, and R. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.
View at: Publisher Site | Google Scholar
J. Yoo, M.-K. Lim, and D.-H. Sheen, “MicroRNA-1915-3P in serum exosome is associated with disease activity of rheumatoid arthritis in Korea. Osteoarthritis and Cartilage,” Elsevier BV, vol. 34, p. S100, 2019.
View at: Publisher Site | Google Scholar
Y. Wang, F. Zheng, G. Gao et al., “MiR-548a-3p regulates inflammatory response via TLR4/NF-κB signaling pathway in rheumatoid arthritis,” Journal of Cellular Biochemistry, vol. 120, pp. 1133–1140, 2019.
View at: Publisher Site | Google Scholar
D. Xu, M. Song, C. Chai et al., “Exosome-encapsulated miR-6089 regulates inflammatory response via targeting TLR4,” Journal of Cellular Physiology, vol. 234, pp. 1502–1511, 2019.
View at: Publisher Site | Google Scholar
L. Wang, C. Wang, X. Jia, and J. Yu, “Circulating exosomal miR-17 inhibits the induction of regulatory T cells via suppressing TGFBR II expression in rheumatoid arthritis,” Cellular Physiology and Biochemistry, vol. 50, pp. 1754–1763, 2018.
View at: Publisher Site | Google Scholar
Z. Chen, H. Wang, Y. Xia, F. Yan, and Y. Lu, “Therapeutic potential of mesenchymal cell–derived miRNA-150-5p–expressing exosomes in rheumatoid arthritis mediated by the modulation of MMP14 and VEGF,” Journal of Immunology, vol. 201, pp. 2472–2482, 2018.
View at: Publisher Site | Google Scholar
Y. Xin, Z. S. Yang, X. Y. Fei et al., “Plasma exosomal MIR-92A are involved in the occurrence and development of bone destruction in ra patients by inhibiting apoptosis of fibroblast-like synoviocytes,” Ann Rheum Dis, vol. 77, 2018.
View at: Publisher Site | Google Scholar
S.-j. Kim, Z. Chen, A. B. Essani et al., “Identification of a novel toll-like receptor 7 endogenous ligand in rheumatoid arthritis synovial fluid that can provoke arthritic joint inflammation,” Arthritis & Rhematology, vol. 68, pp. 1099–1110, 2016.
View at: Publisher Site | Google Scholar
Y. Maeda, N. H. Farina, M. M. Matzelle, P. J. Fanning, J. B. Lian, and E. M. Gravallese, “Synovium-derived microRNAs regulate bone pathways in rheumatoid arthritis,” Journal of Bone and Mineral Research, vol. 32, pp. 461–472, 2017.
View at: Publisher Site | Google Scholar
A. V. Churov, E. K. Oleinik, and M. Knip, “MicroRNAs in rheumatoid arthritis: altered expression and diagnostic potential,” Autoimmunity Reviews, vol. 14, pp. 1029–1037, 2015.
View at: Publisher Site | Google Scholar
I. K. Y. Lam, J. X. Chow, C. S. Lau, and V. S. F. Chan, “MicroRNA-mediated immune regulation in rheumatic diseases,” Cancer Letters, vol. 431, pp. 201–212, 2018.
View at: Publisher Site | Google Scholar
M. F. Romo-García, Y. Bastian, M. Zapata-Zuñiga et al., “Identification of putative miRNA biomarkers in early rheumatoid arthritis by genome-wide microarray profiling: a pilot study,” Gene, vol. 720, p. 144081, 2019.
View at: Publisher Site | Google Scholar
R. E. Crossland, F. Perutelli, K. Bogunia-Kubik et al., “Potential novel biomarkers in chronic graft-versus-host disease,” Frontiers in Immunology, vol. 11, p. 3374, 2020.
View at: Publisher Site | Google Scholar
M. Kurowska-Stolarska, S. Alivernini, L. E. Ballantine et al., “MicroRNA-155 as a proinflammatory regulator in clinical and experimental arthritis,” Proceedings of the National Academy of Sciences of the United States of America, vol. 108, pp. 11193–11198, 2011.
View at: Publisher Site | Google Scholar
K. Smigielska-Czepiel, A. van den Berg, P. Jellema et al., “Comprehensive analysis of miRNA expression in T-cell subsets of rheumatoid arthritis patients reveals defined signatures of naive and memory Tregs,” Genes and Immunity, vol. 15, pp. 115–125, 2014.
View at: Publisher Site | Google Scholar
J. Stanczyk, D. M. Leslie Pedrioli, F. Brentano et al., “Altered expression of microRNA in synovial fibroblasts and synovial tissue in rheumatoid arthritis,” Arthritis and Rheumatism, vol. 58, pp. 1001–1009, 2008.
View at: Publisher Site | Google Scholar
J. Li, Y. Wan, Q. Guo et al., “Altered microRNA expression profile with miR-146a upregulation in CD4+T cells from patients with rheumatoid arthritis,” Arthritis Research & Therapy, vol. 12, p. R81, 2010.
View at: Publisher Site | Google Scholar
G. Ayeldeen, Y. Nassar, H. Ahmed, O. Shaker, and T. Gheita, “Possible use of miRNAs-146a and -499 expression and their polymorphisms as diagnostic markers for rheumatoid arthritis,” Molecular and Cellular Biochemistry, vol. 449, no. 1-2, pp. 145–156, 2018.
View at: Publisher Site | Google Scholar
M. Kriegsmann, T. M. Randau, S. Gravius et al., “Expression of miR-146a, miR-155, and miR-223 in formalin-fixed paraffin-embedded synovial tissues of patients with rheumatoid arthritis and osteoarthritis,” Virchows Archiv, vol. 469, pp. 93–100, 2016.
View at: Publisher Site | Google Scholar
W. Liu, Y.-H. Wu, L. Zhang et al., “MicroRNA-146a suppresses rheumatoid arthritis fibroblast-like synoviocytes proliferation and inflammatory responses by inhibiting the TLR4/NF-kB signaling,” Oncotarget, vol. 9, no. 35, pp. 23944–23959, 2018.
View at: Publisher Site | Google Scholar
Y. Takamura, W. Aoki, A. Satomura, S. Shibasaki, and M. Ueda, “Small RNAs detected in exosomes derived from the MH7A synovial fibroblast cell line with TNF-α stimulation,” PLoS One, vol. 13, pp. 1–18, 2018.
View at: Publisher Site | Google Scholar
Y. Li, S. Zhang, C. Zhang, and M. Wang, “LncRNA MEG3 inhibits the inflammatory response of ankylosing spondylitis by targeting miR-146a,” Molecular and Cellular Biochemistry, vol. 466, pp. 17–24, 2020.
View at: Publisher Site | Google Scholar
M. Filková, B. Aradi, L. Šenolt et al., “Association of circulating miR-223 and miR-16 with disease activity in patients with early rheumatoid arthritis,” Annals of the Rheumatic Diseases, vol. 73, pp. 1898–1904, 2014.
View at: Publisher Site | Google Scholar
K. Murata, H. Yoshitomi, S. Tanida et al., “Plasma and synovial fluid microRNAs as potential biomarkers of rheumatoid arthritis and osteoarthritis,” Arthritis Research & Therapy, vol. 12, p. R86, 2010.
View at: Publisher Site | Google Scholar
T. Niimoto, T. Nakasa, M. Ishikawa et al., “MicroRNA-146a expresses in interleukin-17 producing T cells in rheumatoid arthritis patients,” BMC Musculoskeletal Disorders, vol. 11, p. 209, 2010.
View at: Publisher Site | Google Scholar
A. Singh, P. S. Patro, and A. Aggarwal, “MicroRNA-132, miR-146a, and miR-155 as potential biomarkers of methotrexate response in patients with rheumatoid arthritis,” Clinical Rheumatology, vol. 38, pp. 877–884, 2019.
View at: Publisher Site | Google Scholar
C. Castro-Villegas, C. Pérez-Sánchez, A. Escudero et al., “Circulating miRNAs as potential biomarkers of therapy effectiveness in rheumatoid arthritis patients treated with anti-TNFα,” Arthritis Research & Therapy, vol. 17, p. 49, 2015.
View at: Publisher Site | Google Scholar
Y. Liu, Y. Han, H. Qu, J. Fang, M. Ye, and W. Yin, “Correlation of microRNA expression profile with clinical response to tumor necrosis factor inhibitor in treating rheumatoid arthritis patients: a prospective cohort study,” Journal of Clinical Laboratory Analysis, vol. 33, 2019.
View at: Publisher Site | Google Scholar
B.-p. Qian, M.-l. Ji, Y. Qiu et al., “Identification of serum miR-146a and miR-155 as novel noninvasive complementary biomarkers for ankylosing spondylitis,” Spine (Phila Pa 1976), vol. 41, pp. 735–742, 2016.
View at: Publisher Site | Google Scholar
G. Di, L. Kong, Q. Zhao, and T. Ding, “MicroRNA-146a knockdown suppresses the progression of ankylosing spondylitis by targeting dickkopf 1,” Biomedicine & Pharmacotherapy, vol. 97, pp. 1243–1249, 2018.
View at: Publisher Site | Google Scholar
C. Wei, H. Zhang, C. Wei, and Y. Mao, “Correlation of the expression of miR-146a in peripheral blood mononuclear cells of patients with ankylosing spondylitis and inflammatory factors,” Experimental and Therapeutic Medicine, vol. 14, pp. 5027–5031, 2017.
View at: Publisher Site | Google Scholar
K. Prajzlerová, K. Grobelná, M. Hušáková et al., “Association between circulating miRNAs and spinal involvement in patients with axial spondyloarthritis,” PLoS One, vol. 12, article e0185323, 2017.
View at: Publisher Site | Google Scholar
I. Duroux-Richard, Y.-M. Pers, S. Fabre et al., “Circulating miRNA-125b is a potential biomarker predicting response to rituximab in rheumatoid arthritis,” Mediators of Inflammation, vol. 2014, 9 pages, 2014.
View at: Publisher Site | Google Scholar
L. de la Rica, A. García-Gómez, N. R. Comet et al., “NF-κB-direct activation of microRNAs with repressive effects on monocyte-specific genes is critical for osteoclast differentiation,” Genome Biology, vol. 16, p. 2, 2015.
View at: Publisher Site | Google Scholar
B. Zhang, L. S. Wang, and Y. H. Zhou, “Elevated microRNA-125b promotes inflammation in rheumatoid arthritis by activation of NF-κB pathway,” Biomedicine & Pharmacotherapy, vol. 93, pp. 1151–1157, 2017.
View at: Publisher Site | Google Scholar
L. Sun, J. X. Lian, and S. Meng, “MiR-125a-5p promotes osteoclastogenesis by targeting TNFRSF1B,” Cellular & Molecular Biology Letters, vol. 24, 2019.
View at: Publisher Site | Google Scholar
C. Perez-Sanchez, P. Font-Ugalde, P. Ruiz-Limon, C. Lopez-Pedrera, M. C. Castro-Villegas, and M. C. Abalos-Aguilera, “Circulating microRNAs as potential biomarkers of disease activity and structural damage in ankylosing spondylitis patients,” Human Molecular Genetics, vol. 27, pp. 875–890, 2018.
View at: Publisher Site | Google Scholar
Y. Jiang and L. Wang, “Role of histone deacetylase 3 in ankylosing spondylitis via negative feedback loop with microRNA-130a and enhancement of tumor necrosis factor-1α expression in peripheral blood mononuclear cells,” Molecular Medicine Reports, vol. 13, pp. 35–40, 2016.
View at: Publisher Site | Google Scholar
Z.-C. Li, N. Han, X. Li, G. Li, Y.-Z. Liu, and G.-X. Sun, “Decreased expression of microRNA-130a correlates with TNF-α in the development of osteoarthritis,” International Journal of Clinical and Experimental Pathology, vol. 8, no. 3, pp. 2555–2564, 2015.
View at: Google Scholar
X. Tang, K. Yin, H. Zhu et al., “Correlation between the expression of microRNA-301a-3p and the proportion of Th17 cells in patients with rheumatoid arthritis,” Inflammation, vol. 39, pp. 759–767, 2016.
View at: Publisher Site | Google Scholar
B. Ding, S. Xu, X. Sun, J. Gao, W. Nie, and H. Xu, “miR-18a-3p encourages apoptosis of chondrocyte in osteoarthritis via HOXA1 pathway,” Curr Molecular Pharmacology, vol. 13, 2020.
View at: Publisher Site | Google Scholar
B. Li, L. Bai, P. Shen, Y. Sun, Z. Chen, and Y. Wen, “Identification of differentially expressed microRNAs in knee anterior cruciate ligament tissues surgically removed from patients with osteoarthritis,” International Journal of Molecular Medicine, vol. 40, pp. 1105–1113, 2017.
View at: Publisher Site | Google Scholar
W. Yang, X. Yan, Q. Xia et al., “Predisposition of six well-characterized microRNAs to syndesmophytes among Chinese patients with ankylosing spondylitis,” Modern Rheumatology, vol. 29, pp. 173–180, 2019.
View at: Publisher Site | Google Scholar
M. Trenkmann, M. Brock, R. E. Gay, B. A. Michel, S. Gay, and L. C. Huber, “Tumor necrosis factor α-induced microRNA-18a activates rheumatoid arthritis synovial fibroblasts through a feedback loop in NF-κB signaling,” Arthritis and Rheumatism, vol. 65, pp. 916–927, 2013.
View at: Publisher Site | Google Scholar
R. Kong, J. Gao, Y. Si, and D. Zhao, “Combination of circulating miR-19b-3p, miR-122-5p and miR-486-5p expressions correlates with risk and disease severity of knee osteoarthritis,” American Journal of Translational Research, vol. 9, no. 6, pp. 2852–2864, 2017.
View at: Google Scholar
L. Duan, D. Duan, W. Wei et al., “MiR-19b-3p attenuates IL-1β induced extracellular matrix degradation and inflammatory injury in chondrocytes by targeting GRK6,” Molecular and Cellular Biochemistry, vol. 459, pp. 205–214, 2019.
View at: Publisher Site | Google Scholar
M. P. Gantier, H. J. Stunden, C. E. McCoy et al., “A miR-19 regulon that controls NF-κB signaling,” Nucleic Acids Research, vol. 40, pp. 8048–8058, 2012.
View at: Publisher Site | Google Scholar
L. Philippe, G. Alsaleh, G. Suffert et al., “TLR2 expression is regulated by microRNA miR-19 in rheumatoid fibroblast-like synoviocytes,” Journal of Immunology, vol. 188, pp. 454–461, 2012.
View at: Publisher Site | Google Scholar
J. Huang, G. Song, Z. Yin, X. Luo, and Z. Ye, “Elevated miR-29a expression is not correlated with disease activity index in PBMCs of patients with ankylosing spondylitis,” Modern Rheumatology, vol. 24, pp. 331–334, 2014.
View at: Publisher Site | Google Scholar
J. Liu, D. Fei, J. Xing, and J. Du, “MicroRNA-29a inhibits proliferation and induces apoptosis in rheumatoid arthritis fibroblast-like synoviocytes by repressing STAT3,” Biomedicine & Pharmacotherapy, vol. 96, pp. 173–181, 2017.
View at: Publisher Site | Google Scholar
J. Huang, G. Song, Z. Yin, Z. Fu, and L. Zhang, “Altered expression of microRNAs targeting Dkk-1 in peripheral blood mononuclear cells of patients with ankylosing spondylitis,” Cent Eur J Immunol., vol. 44, pp. 59–64, 2019.
View at: Publisher Site | Google Scholar
X. Li, Q. Lv, L. Tu et al., “Aberrant expression of microRNAs in peripheral blood mononuclear cells as candidate biomarkers in patients with axial spondyloarthritis,” International Journal of Rheumatic Diseases, vol. 22, pp. 1188–1195, 2019.
View at: Publisher Site | Google Scholar
Q. Lv, Q. Li, P. Zhang et al., “Disorders of microRNAs in peripheral blood mononuclear cells: as novel biomarkers of ankylosing spondylitis and provocative therapeutic targets,” BioMed Research International, vol. 2015, 7 pages, 2015.
View at: Publisher Site | Google Scholar
C.-H. Huang, J. C.-C. Wei, W.-C. Chang 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,” The Journal of Rheumatology, vol. 41, pp. 1104–1111, 2014.
View at: Publisher Site | Google Scholar
Y. C. Zou, Y. P. Gao, H. D. Yin, and G. Liu, “Serum miR-21 expression correlates with radiographic progression but also low bone mineral density in patients with ankylosing spondylitis: a cross-sectional study,” Innate Immunity, vol. 25, pp. 314–321, 2019.
View at: Publisher Site | Google Scholar
L. Dong, X. Wang, J. Tan et al., “Decreased expression of microRNA-21 correlates with the imbalance of Th17 and Treg cells in patients with rheumatoid arthritis,” Journal of Cellular and Molecular Medicine, vol. 18, pp. 2213–2224, 2014.
View at: Publisher Site | Google Scholar
F. Balzano, M. Deiana, S. D. Giudici et al., “MicroRNA expression analysis of centenarians and rheumatoid arthritis patients reveals a common expression pattern,” International Journal of Medical Sciences, vol. 14, pp. 622–628, 2017.
View at: Publisher Site | Google Scholar
S. Yang, S. Jiang, Y. Wang, S. Tu, Z. Wang, and Z. Chen, “Interleukin 34 upregulation contributes to the increment of microRNA 21 expression through STAT3 activation associated with disease activity in rheumatoid arthritis,” The Journal of Rheumatology, vol. 43, no. 7, pp. 1312–1319, 2016.
View at: Publisher Site | Google Scholar
Y. C. Zou, L. M. Yan, Y. P. Gao, Z. Y. Wang, and G. Liu, “miR-21 may act as a potential mediator between inflammation and abnormal bone formation in ankylosing spondylitis based on TNF-α concentration-dependent manner through the JAK2/STAT3 pathway,” Dose-Response, vol. 18, 2020.
View at: Publisher Site | Google Scholar
T. Sugatani, J. Vacher, and K. A. Hruska, “A microRNA expression signature of osteoclastogenesis,” Blood, vol. 117, pp. 3648–3657, 2011.
View at: Publisher Site | Google Scholar
M. N. Magrey, T. Haqqi, and A. Haseeb, “Identification of plasma microRNA expression profile in radiographic axial spondyloarthritis—a pilot study,” Clinical Rheumatology, vol. 35, pp. 1323–1327, 2016.
View at: Publisher Site | Google Scholar
N. S. Lai, H. C. Yu, H. C. Chen, C. L. Yu, H. B. Huang, and M. C. Lu, “Aberrant expression of microRNAs in T cells from patients with ankylosing spondylitis contributes to the immunopathogenesis,” Clinical and Experimental Immunology, vol. 173, pp. 47–57, 2013.
View at: Publisher Site | Google Scholar
P. Reyes-Loyola, P. Rodríguez-Henríquez, M. A. Ballinas-Verdugo et al., “Plasma let-7i, miR-16, and miR-221 levels as candidate biomarkers for the assessment of ankylosing spondylitis in Mexican patients naïve to anti-TNF therapy,” Clinical Rheumatology, vol. 38, pp. 1367–1373, 2019.
View at: Publisher Site | Google Scholar
B. V. J. Cuppen, on behalf of all SRU investigators, M. Rossato et al., “Can baseline serum microRNAs predict response to TNF-alpha inhibitors in rheumatoid arthritis?” Arthritis Research & Therapy, vol. 18, p. 189, 2016.
View at: Publisher Site | Google Scholar
S. Zhu, W. Pan, X. Song et al., “The microRNA miR-23b suppresses IL-17-associated autoimmune inflammation by targeting TAB2, TAB3 and IKK-α,” Nature Medicine, vol. 18, pp. 1077–1086, 2012.
View at: Publisher Site | Google Scholar
J. Hu, C. Zhai, J. Hu et al., “MiR-23a inhibited IL-17-mediated proinflammatory mediators expression via targeting IKKα in articular chondrocytes,” International Immunopharmacology, vol. 43, pp. 1–6, 2017.
View at: Publisher Site | Google Scholar
H. M. A. Elsayed, W. S. Khater, A. A. Ibrahim, M. S. E. Hamdy, and N. A. Morshedy, “MicroRNA-146a expression as a potential biomarker for rheumatoid arthritis in Egypt,” Egypt J Med Hum Genet., vol. 18, pp. 173–179, 2017.
View at: Publisher Site | Google Scholar
Z. Z. Chen, X. D. Zhang, Y. Chen, and Y. B. Wu, “The role of circulating miR-146a in patients with rheumatoid arthritis treated by Tripterygium wilfordii Hook F,” Med (United States), vol. 96, 2017.
View at: Publisher Site | Google Scholar
A. Abou-Zeid, M. Saad, and E. Soliman, “MicroRNA 146a expression in rheumatoid arthritis: association with tumor necrosis factor-alpha and disease activity,” Genetic Testing and Molecular Biomarkers, vol. 15, pp. 807–812, 2011.
View at: Publisher Site | Google Scholar
N. Mookherjee and H. S. El-Gabalawy, “High degree of correlation between whole blood and PBMC expression levels of miR-155 and miR-146a in healthy controls and rheumatoid arthritis patients,” Journal of Immunological Methods, vol. 400–401, pp. 106–110, 2013.
View at: Publisher Site | Google Scholar
M. J. Ormseth, J. F. Solus, K. C. Vickers, A. M. Oeser, P. Raggi, and C. M. Stein, “Utility of select plasma microRNA for disease and cardiovascular risk assessment in patients with rheumatoid arthritis,” The Journal of Rheumatology, vol. 42, pp. 1746–1751, 2015.
View at: Publisher Site | Google Scholar
H. Wang, W. Peng, X. Ouyang, W. Li, and Y. Dai, “Circulating microRNAs as candidate biomarkers in patients with systemic lupus erythematosus,” Translational Research, vol. 160, pp. 198–206, 2012.
View at: Publisher Site | Google Scholar
Q. Zhou, S. Haupt, J. T. Kreuzer et al., “Decreased expression of miR-146a and miR-155 contributes to an abnormal Treg phenotype in patients with rheumatoid arthritis,” Annals of the Rheumatic Diseases, vol. 74, pp. 1265–1274, 2015.
View at: Publisher Site | Google Scholar
V. Anaparti, I. Smolik, X. Meng, V. Spicer, N. Mookherjee, and H. El-Gabalawy, “Whole blood microRNA expression pattern differentiates patients with rheumatoid arthritis, their seropositive first-degree relatives, and healthy unrelated control subjects,” Arthritis Research & Therapy, vol. 19, p. 249, 2017.
View at: Publisher Site | Google Scholar
X. Li, F. Tian, and F. Wang, “Rheumatoid arthritis-associated microrna-155 targets socs1 and upregulates TNF-α and IL-1β in PBMCs,” International Journal of Molecular Sciences, vol. 14, pp. 23910–23921, 2013.
View at: Publisher Site | Google Scholar
L. Long, P. Yu, Y. Liu et al., “Upregulated microRNA-155 expression in peripheral blood mononuclear cells and fibroblast-like synoviocytes in rheumatoid arthritis,” Clinical & Developmental Immunology, vol. 2013, pp. 1–10, 2013.
View at: Publisher Site | Google Scholar
R. S. Abdul-Maksoud, A. M. Sediq, A. A. A. Kattaia et al., “Serum miR-210 and miR-155 expression levels as novel biomarkers for rheumatoid arthritis diagnosis,” British Journal of Biomedical Science, vol. 74, pp. 209–213, 2017.
View at: Publisher Site | Google Scholar
Y. Wang, T. Feng, S. Duan et al., “miR-155 promotes fibroblast-like synoviocyte proliferation and inflammatory cytokine secretion in rheumatoid arthritis by targeting FOXO3a,” Experimental and Therapeutic Medicine, vol. 19, pp. 1288–1296, 2020.
View at: Publisher Site | Google Scholar
M. Rajasekhar, A. M. Olsson, K. J. A. Steel et al., “MicroRNA-155 contributes to enhanced resistance to apoptosis in monocytes from patients with rheumatoid arthritis,” Journal of Autoimmunity, vol. 79, pp. 53–62, 2017.
View at: Publisher Site | Google Scholar
A. Elmesmari, A. R. Fraser, C. Wood et al., “MicroRNA-155 regulates monocyte chemokine and chemokine receptor expression in rheumatoid arthritis,” Rheumatol (United Kingdom), vol. 55, pp. 2056–2065, 2016.
View at: Publisher Site | Google Scholar
C. Zhang, L. Fang, X. Liu et al., “miR-22 inhibits synovial fibroblasts proliferation and proinflammatory cytokine production in RASF via targeting SIRT1,” Gene, vol. 724, 2020.
View at: Publisher Site | Google Scholar
L. Ouboussad, L. Hunt, E. M. A. Hensor et al., “Profiling microRNAs in individuals at risk of progression to rheumatoid arthritis,” Arthritis Research & Therapy, vol. 19, p. 288, 2017.
View at: Publisher Site | Google Scholar
J. Lin, R. Huo, L. Xiao et al., “A novel p53/microRNA-22/Cyr61 axis in synovial cells regulates inflammation in rheumatoid arthritis,” Arthritis & Rhematology, vol. 66, pp. 49–59, 2014.
View at: Publisher Site | Google Scholar
K. Xu, P. Xu, J. F. Yao, Y. G. Zhang, W. K. Hou, and S. M. Lu, “Reduced apoptosis correlates with enhanced autophagy in synovial tissues of rheumatoid arthritis,” Inflammation Research, vol. 62, pp. 229–237, 2013.
View at: Publisher Site | Google Scholar
K. Murata, M. Furu, H. Yoshitomi et al., “Comprehensive microRNA analysis identifies miR-24 and miR-125a-5p as plasma biomarkers for rheumatoid arthritis,” PLoS One, vol. 8, article e69118, 2013.
View at: Publisher Site | Google Scholar
A. S. Abo ElAtta, Y. B. M. Ali, I. H. Bassyouni, and R. M. Talaat, “Upregulation of miR-221/222 expression in rheumatoid arthritis (RA) patients: correlation with disease activity,” Clinical and Experimental Medicine, vol. 19, pp. 47–53, 2019.
View at: Publisher Site | Google Scholar
S. Yang and Y. Yang, “Downregulation of microRNA-221 decreases migration and invasion in fibroblast-like synoviocytes in rheumatoid arthritis,” Molecular Medicine Reports, vol. 12, pp. 2395–2401, 2015.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Joanna Wielińska et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

628

Downloads

600

Citations