Can baseline serum microRNAs predict response to TNF-alpha inhibitors in rheumatoid arthritis?
Arthritis Research & Therapy volume 18, Article number: 189 (2016)
In rheumatoid arthritis, prediction of response to TNF-alpha inhibitor (TNFi) treatment would be of clinical value. This study aims to discover miRNAs that predict response and aims to replicate results of two previous studies addressing this topic.
From the observational BiOCURA cohort, 40 adalimumab- (ADA) and 40 etanercept- (ETN) treated patients were selected to enter the discovery cohort and baseline serum profiling on 758 miRNAs was performed. The added value of univariately selected miRNAs (p < 0.05) over clinical parameters in prediction of response was determined by means of the area under the receiver operating characteristic curve (AUC-ROC). Validation was performed by TaqMan single qPCR assays in 40 new patients.
Expression of miR-99a and miR-143 predicted response to ADA, and miR-23a and miR-197 predicted response to ETN. The addition of miRNAs increased the AUC-ROC of a model containing only clinical parameters for ADA (0.75 to 0.97) and ETN (0.68 to 0.78). In validation, none of the selected miRNAs significantly predicted response. miR-23a was the only overlapping miRNA compared to the two previous studies, however inversely related with response in one of these studies. The reasons for the inability to replicate previously proposed miRNAs predicting response to TNFi and replicate those from the discovery cohort were investigated and discussed.
To date, no miRNA consistently predicting response to TNFi therapy in RA has been identified. Future studies on this topic should meet a minimum of standards in design that are addressed in this study, in order to increase the reproducibility.
Rheumatoid arthritis (RA) is a chronic, disabling disease that mainly affects the synovial joints, with a prevalence of 0.5–1.0 % in Western countries [1, 2]. The introduction of TNF-α-inhibiting therapy (TNFi), such as adalimumab (ADA) and etanercept (ETN), has dramatically improved the outlook for RA patients. Nevertheless, a substantial proportion of patients (approximately 30–40 %) fail to respond to TNFi therapy [3, 4]. As we cannot predict before initiation of therapy which patients will be non-responders , TNFi treatment is administered in a trial and error approach. However in the time frame from initiation of therapy until response can be judged, which is usually 3–6 months later, non-responding patients suffer from uncontrolled disease with possible joint damage and the potential harmful side effects from treatment. The challenge is therefore to identify responders and non-responders to TNFi beforehand, so that TNFi use or considering alternatives can be encouraged.
microRNAs (miRNAs) are a large family of highly conserved noncoding genes that play a fundamental role in biological processes by controlling protein expression [6–8]. miRNAs execute these actions by binding to protein-coding messenger RNAs (mRNAs), resulting in translational repression or mRNA degradation . Besides intracellulary, miRNAs are also found in several biological fluids, including saliva, plasma, serum and urine, either circulating in conjunction with specific carrier proteins or enclosed in extracellular vesicles [9, 10]. Exploring the use of circulating miRNAs as biomarkers for diseases has gained momentum in recent years because of the easy accessibility, the associations with specific disease conditions and their good stability [10, 11]. In RA, a systemic inflammatory disease primarily manifesting in the joints, biomarkers in the circulation would intuitively not be the most relevant compartment. However, the levels of miRNAs are frequently higher in the circulation than in the synovial fluid and correlate with disease activity in RA, indicating that the systemic compartment provides a useful compartment for studying the ongoing pathophysiological processes  In addition, abnormal expression of both synovial and systemic miRNAs have been linked to disease activity and pathogenesis, even though their direct targets are not always known [13–19]. Three recent studies focused on the prediction of response to therapy in RA by using circulating miRNAs, of which two investigated response to therapy with TNFi [20, 21] and one to rituximab . According to these studies, promising predictors for TNFi therapy were miR-22 , miR-23a , miR-223 , and miR-886 . Circulating miR-23a seems of particular interest, since it was the only identified candidate biomarker that was overlapping among both studies in univariate analyses. However, upregulation of miR23a was found in whole blood , whereas a downregulation was found in serum  of future responders.
In this study we explored the serum miRNAs associated with good and bad response to TNFi therapy, in order to replicate the results that have been published before. In addition to the previous studies performed, we involved clinical parameters in the prediction and attempted to validate the miRNAs and prediction models in a separate cohort.
Clinical data collection
Patients initiating ADA or ETN therapy were selected from the “Biologicals and Outcome Compared and predicted Utrecht region in Rheumatoid Arthritis” (BiOCURA) study. BiOCURA is an observational cohort, in which RA patients eligible for biological treatment according to regular clinical practice were enrolled and followed after start of treatment, in one academic hospital and seven regional hospitals in the Netherlands (see Acknowledgements). Re-inclusion after switching to a different biological treatment was possible, at which patients re-entered baseline again. The study was approved by the local ethics committee of the University Medical Center Utrecht and the institutional review boards of the participating centers, and was performed in accordance with the Declaration of Helsinki. Informed consent was obtained from each patient.
Trained nurses gathered all data during a dedicated visit, which included all clinical parameters, joint counts and collection of blood. Visits were scheduled at baseline (before initiation) and after 3, 6, and 12 months of treatment. Disease activity was assessed using the disease activity score based on a 28-joint count (DAS28)  and subsequently the European League Against Rheumatism (EULAR) response was calculated . This study design, allowed the determination of a clinical response of each patient, on the basis of three EULAR responses over the course of 1 year.
Two separate cohorts were composed from the available patients in the BiOCURA study: a discovery cohort was used to screen the expression of a full panel of 758 miRNAs, while a validation cohort was used to test if the results found in the discovery phase were reproducible. The discovery cohort was formed by selecting the most extreme patients regarding clinical response, among all ADA- and ETN-treated patients included between June 2009 (start of BiOCURA) and October 2012 (n = 74 ADA and n = 68 ETN). The top responding patients (n = 20 for both ADA and ETN, from now on called “responders”), were identified by the selection of patients with the best three EULAR responses over the course of 1 year. The selection of bad responders (n = 20 for both ADA and ETN, from now on called “non-responders”), was based on the most negative EULAR responses over the course of 1 year and/or (early) discontinuation of TNFi treatment due to inefficacy. Patients with a baseline DAS28 < 2.6 were excluded from the analysis in order to reduce the chance that limited improvement in DAS28 resulted in misclassification as EULAR non-responders. For validation of results, responders (n = 10 for both ADA and ETN) and non-responders (n = 10 for both ADA and ETN) were selected using the same criteria as in the discovery cohort, among patients included from October 2012 until June 2015 (n = 25 ADA and n = 40 ETN). Since the validation cohort was smaller, relatively more patients were selected and the differences in clinical outcome between responders and non-responders were less extreme. The baseline characteristics for responders and non-responders are shown in Table 1 and for responders and non-responders split per cohort in Additional file 1. Additionally, the baseline characteristics for the discovery and validation cohorts are shown in Additional file 2.
Blood processing and RNA extraction
Blood was collected in Vacutainer® SSTII tubes (BD, Franklin Lakes, NJ, USA) and processed immediately after clotting. Samples were centrifuged for 10 min at 1500 g at room temperature and serum was aliquoted and stored at −80 °C until use. RNA was extracted from 240 μl of serum using the miRcury RNA Isolation kit for Biofluids (Exiqon), according to the manufacturer’s instructions. During extraction, 300 pg of a synthetic miRNA (Arabidopsis thaliana ath-miR-159a) was added to each sample as a spike-in to monitor technical variability along the isolation procedure and for later normalization.
miRNA profiling was performed by TaqMan RT-qPCR on the OpenArray platform (Life Technologies, Carlsbad, CA, USA). This method allows the simultaneous analysis of 758 miRNAs, split into two equal pools (A and B). Manufacturer’s instructions were followed with minor adjustments. Briefly, 2.5ul of isolated serum RNA was reverse-transcribed by using the miRNA multiplex RT primers pools, either v2.1 for pool A or v3.0 for pool B, and the TaqMan miRNA reverse transcription kit (Life Technologies). RT products were pre-amplified using the Megaplex PreAmp Primers pools A and B in the presence of the TaqMan PreAmp Master Mix (Life Technologies), by using the following thermal cycler conditions: 10 min, 95 °C; 2 min, 55 °C; 2 min, 72 °C and 16 cycles of 15 sec, 95 °C and 4 min, 60 °C and one single cycle of 10 min, 96 °C. The miRNA OpenArray profiling was performed on the amplified cDNA, diluted to 1:40, with 0.1 × TE buffer pH 8.0 and subsequently 1:2 by using the TaqMan OpenArray Master Mix on the QuantStudio 12 K Flex Real-Time PCR System (Life Technologies).
miRNA profiling data was analyzed using the Relative Quantification application in the online accessible Thermo Fisher Cloud (https://apps.thermofisher.com/apps/dashboard/), using the relative threshold cycle (Crt) and the comparative threshold cycle method . Briefly, miRNA expression was calculated after normalization by exogenous ath-miR-159a spike-in (ΔCrt = Crt mean target – Crt mean miR-159a). The relative fold change (FC) of each sample was determined by setting the FC of a random ADA or ETN non-responder sample at 1, and calculating the FC compared to this reference (FC = 2–ΔΔCrt, where –ΔΔCrt = ΔCrt reference – ΔCrt sample). Low expressed miRNAs, i.e., having Crt higher than 27 were set to 27, and samples with a low amplification quality (i.e., amplification score < 1.24) were excluded from the analysis.
Individual miRNA analysis
miRNA-specific TaqMan Real-Time quantitative PCR (RT-qPCR) assays were purchased from Life Technologies for hsa-miR-23a-3p (ID 000399), hsa-miR-99a-5p (ID 000435), hsa-miR-143-3p (ID 002249), hsa-miR-197-3p (ID 000497), and for the exogenous control ath-miR-159a (ID 000338). From 2.5 μl baseline serum RNA, cDNA was synthesized by using individual miRNA-specific RT primers contained in the miRNA assay in the presence of 3.3 U/μl MultiScribe RT enzyme (Life Technologies), by using the following thermal cycler conditions: 10 min, 4 °C; 30 min,16 °C; 30 min, 42 °C; and 5 min, 85 °C. Circulating miRNA levels were quantified in duplicate from 3 μl cDNA, with TaqMan Fast Advance Master Mix and specific primers of the miRNA assay, using the following amplification condition on the Quantstudio 12 K Flex Real-Time PCR system: 2 min, 50 °C; 20 sec, 95 °C; 40 cycles of 1 sec, 95 °C; and 20 sec, 60 °C. RT-qPCR data were calculated as described above, with the difference that baseline threshold cycles (Ct) were used.
Differential expression of miRNAs between responders and non-responders was calculated separately for ADA and ETN by means of an independent sample t tests on the –ΔΔCrt/–ΔΔCt, with a threshold for significance of 0.05 (uncorrected p value). The levels of differentially expressed miRNAs were plotted in GraphPad Prism (GraphPad, La Jolla, CA, USA) as FC of responders versus non-responders. Validation was considered successful when both the t test was significant and plots of the FC showed the same direction (i.e., up/downregulation).
In order to determine the added value of the miRNAs over clinical parameters, we built two prediction models for each treatment, using multivariable logistic regression. The first model consisted of all baseline clinical parameters that were significantly different between responders and non-responders (the “clinical model”). The second model included the clinical parameters and the selected miRNAs (–ΔΔCrt values) (the “combined model”). Per model, the area under the receiver operating characteristic curve (AUC-ROC) was calculated as an indicator of the predictive ability. We considered an AUC-ROC of < 0.7 limitedly, 0.7–0.8 moderately and > 0.8 highly predictive of response. The sensitivity and specificity were calculated for the best cutoff value per model, according to Youden’s index . Evaluation of the added value of miRNAs was based on the increase of predictive abilities when switching from the clinical to the combined model. In order to validate the findings from multivariable analysis, the prediction rules of the clinical and combined models were applied in the validation cohort, thereby freezing the regression coefficients of the individual parameters from the original model. Again, the AUC-ROC, the sensitivity and specificity were calculated to interpret the added value of miRNAs over clinical parameters alone.
Identification of miRNAs as predictor of TNFi response
We analyzed the profile of miRNAs present in the circulation of responders versus non-responders with a broad panel of 758 miRNAs. In the discovery cohort (n = 80), four miRNAs were significantly differentially expressed between responders and non-responders: high and low baseline levels of respectively miR-99a and miR-143 predicted response to ADA, while patients with high levels of miR-23a and miR-197 more frequently responded to ETN (Fig. 1). Expression values of patients in the discovery cohort were also plotted for miRNAs proposed by the previous studies as predictors for response (Additional file 3). miR-23a in ETN-treated patients was the only miRNA overlapping between this study and previous ones published.
Since the measurement of miRNAs can be costly when incorporated in clinical practice, we wanted to rule out the possibility that the miRNAs identified do not increase the magnitude of prediction that is already possible based on clinical parameters. We therefore compared the predictive abilities of models based on clinical parameters alone, and clinical parameters together with the miRNA expression levels. The clinical characteristics that were used, were those that presented a significant difference between responders and non-responders at baseline, namely the DAS28 (p = 0.01), swollen joint count (SJC, p = 0.03) and amount of glucocorticoid (GC) users (p = 0.01) for patients treated with ADA, and C-reactive protein (CRP, p = 0.03) for those treated with ETN (Additional file 1). The predictive properties of these models without and with miRNAs are shown in Table 2. The clinical model for ADA showed a moderate predictive value (AUR-ROC 0.75), that was increased by the addition of miR-99a and miR-143 in the combined model (AUR-ROC 0.97). For ETN, the CRP alone was only limitedly able to predict response (AUC-ROC 0.68), however, the predictive value increased by the addition of miR-23a and miR-197 in the combined model (AUC-ROC 0.78).
Since replication in (prognostic) research is key to prove validity, we tried to confirm our results in an additional cohort of 40 patients. The differentially abundant miRNAs from the discovery cohort were analyzed in the validation cohort by using single RT-qPCR assays (Fig. 2). None of the miRNAs could significantly predict the response to TNFi in the validation cohort (p > 0.05). For miR-99a and miR-143 in ADA users, inverse directions were seen compared to the results in the discovery cohort.
Multivariable analyses did not confirm the predictive abilities of the combined models found in the discovery phase (Additional file 4). Contrary to what was observed in the discovery phase, application of the prediction model for ADA including only clinical parameters showed better predictive abilities than the corresponding combined model (AUC-ROC from 0.93 to 0.57). This is most likely explained by the inverse relationship of the miRNAs and response in the validation compared to the discovery cohort. For ETN, the addition of miRNA added to the prediction of response (AUC-ROC from 0.59 to 0.66), generating a model that was only limitedly predictive.
Factors that can contribute to the inability to validate findings in multiple cohorts
It is of importance to investigate which factors are involved in the inability to validate findings in multiple cohorts, since these issues could also be applicable to other studies aiming to identify miRNAs predicting the response to TNFi therapy. Possible (combinations of) factors could be, but are not limited to, the usage of different miRNA detection methods, selection of false positive results in the discovery phase, or clinical parameters influencing the relationship between miRNA levels and response.
Despite that miRNA expression analysis in the discovery and validation phase were both based on the same method of detection, i.e., miRNA-specific retrotranscription combined with TaqMan-based RT-qPCR, the scale (high throughput versus single assay) of the techniques sufficiently varies. To evaluate whether these differences could impact the final result, we performed a technical replication in all 40 ADA or 40 ETN samples from the discovery cohort using single assays for the four selected miRNAs (as described in "Methods" - miRNA analyses - Individual miRNA analysis).
Correlations of the results obtained by the profiling versus those measured by single assay were assessed by calculating the Spearman correlation (r) between the normalized detection levels (ΔCrt and ΔCt respectively) without excluding samples based on amplification scores (Fig. 3). The correlations of test-retest values ranged from 0.45 to 0.88 (all p < 0.0001), which can be considered ranging from reasonable to good, thus demonstrating that the two analyses are concordant. Comparison of the single assay miRNA expression levels between responders versus non-responders confirmed that miR-143 was significantly lower in ADA-responders, whereas the other miRNAs showed the same direction as in the profiling, though did not reach significance (Additional file 5).
To verify whether these results could be related to a false discovery rate, we recalculated the differential expression for all miRNAs in the discovery cohort while applying the Benjamini and Hochberg false discovery rate (B&H FDR), which showed corrected p values of 1.00 for all miRNAs. Considering none of the miRNAs was significantly different after correction, there is the possibility that only false positive results were selected in the discovery phase.
Another possible explanation why we were unable to replicate the results from previous studies and our discovery cohort is that clinical parameters interact with miRNA levels and these clinical parameters were not equally distributed between the cohorts. Differences in case-mix between cohorts that were unaccounted for (see Additional file 2) would then lead to different estimations of each miRNA and response. Despite this, an adjustment for these clinical parameters would then give comparable estimations for the miRNAs involved. We investigated this theory by running a crude model of response including the specific miRNA only, and an adjusted model considering both the miRNA and clinical parameters, and run these models for the two cohorts analyzed (Table 3). Despite the adjustment, the odds ratios (OR) of these miRNAs for response were still (very) different between the discovery and validation cohort. This indicates that the clinical parameters do not explain why results could not be validated. On the other hand, these analyses showed that clinical parameters have a strong effect on the association between miRNA levels and the response to therapy, as indicated by the (relatively large) differences between crude ORs and adjusted ORs. Therefore, adjustment for clinical parameters will contribute to externalization of results to cohorts with a different case-mix, as is a common occurrence in a heterogeneous disease such as (established) RA. Considering that clinical parameters could affect the expression of miRNAs, miRNA levels may be, to a certain extent, a representation of patient’s clinical characteristics. Therefore, we evaluated the correlation of clinical parameters and miRNAs levels, as measured by single assays, irrespective of response The analysis revealed that all miRNAs associated with either CRP or erythrocyte sedimentation rate (ESR) (Additional file 6). However, since the explained variance of each miRNA by clinical parameters was less than 35 %, miRNA levels are not a complete reflection of clinical characteristics and can thus contain informative additional information.
Since the exact origin of circulating of miRNAs is unknown and blood cells have been proposed as a key source , we wanted to rule out the possibility that the serum levels of miRNAs are a reflection of the composition of circulating leukocytes. We therefore correlated the levels of serum miRNAs with the percentage of peripheral blood leukocyte subsets, as measured by flow cytometry in 20 randomly selected patients (Additional file 7). Out of all comparisons that were considering the surface markers CD3, CD4, CD8, CD14, CD19, CD45, and CD16 + 56 and the FC of all four miRNAs, one significant correlation was observed between the levels of miR-197, as measured in the profiling, and the percentage of natural killer-like T cells (r = 0.587, p = 0.008). However, this association was not reproduced when considering the miRNA levels measured by the single assay (r = 0.083, p = 0.831). Even though we cannot exclude the contribution of other rare cell subsets that were not identified, these results suggested that the levels of serum miRNAs in responders versus non-responders are not dependent on the composition of circulating leukocytes, thus making a different leukocyte composition an unlikely cause for the inconsistency between the two cohorts analyzed and the other studies.
Prediction of TNFi response is needed for a more personalized approach in RA treatment. Since two previous studies have addressed this question and identified candidate miRNAs [20, 21], we aimed at verifying whether these could be validated in an independent cohort, and eventually, whether we could find new predictors. High values of circulating miR23a were univariately predictive of response in our study and one previous work , however, miR-23a was inversely related to response in a third study  and is thus not a consistent predictor. On the basis of the profiling results, four miRNAs showed an added value to clinical parameters in predicting response to TNFi. However, these miRNAs could not be validated in a separate cohort of consecutively included patients.
Several factors could have contributed to the inability to replicate findings from previous works and our own selection of miRNAs. A first possible contributing factor is technical variability of the techniques employed for the miRNA analysis. Within our study and the previous studies, pre-analytical and analytical protocols were standardized, as is considered to contribute to more reliable measurements in miRNA studies . Correlations between the two protocols used to analyze miRNA expression were reasonably good, though could explain (some) difference in the outcome. Indeed, a significant differential expression between responders and non-responders could only be replicated for miR-143, which was the miRNA with the best amplification quality in the array among those selected (i.e., amplification score ≥ 1.24 in all samples). The previous studies used comparable techniques for biomarkers discovery: TaqMan single miRNA assays  and the same platform as in our discovery step (TaqMan OpenArray, Pool A) . Furthermore, the OpenArray appears one of the most reliable high-throughput technique for miRNA analysis , and it was successfully used for profiling of serum miRNA in multiple studies by us [Chouri E, et al, manuscript in preparation] and others [30, 31]. Interestingly, miR-23a was positively related to response both in our discovery cohort and in a previous study that was also using the OpenArray as profiling platform . On the contrary, the absence of a relation or inverse relation with response was found in all cases employing TaqMan single assays . Altogether, the analytical techniques and their technical variability are unlikely to be explanatory for all the differences within and between the studies, though they might be a contributing factor.
The lack of correction for multiple testing in the discovery phase might have led to a subsequent (wrongful) selection of candidate miRNAs. However, a correction for multiple testing was not applied, because this could be too strict when trying to replicate findings already identified by others and would also have increased the chance of excluding potentially useful miRNAs (false negative results). In addition, a separate validation step will reveal which of the (less strictly) selected miRNAs has a true biological meaning, and thus compensates for the absence of correction for false discovery rate. Yet, the fact that the application of B&H FDR resulted in a p value of 1.00 for all analyzed miRNAs is suggestive for the possibility that false positive miRNAs were selected. The identification of miRNAs in the previous studies might also be based on false positive results, since no correction for multiple testing was performed and results were not validated in a separate cohort. Not unimportant, if all proposed predictors so far are indeed false positives, this might entail an absence of biologically relevant miRNAs for the prediction of TNFi response in RA.
A third contributing factor for the inability to validate our findings could be related to differences in baseline clinical parameters, which could influence the relation between miRNAs and response. In additional analyses, however, we demonstrated that the clinical parameters were not able to explain why results could not be replicated between the discovery and the validation cohort. However, we did observe that the predictive values of miRNAs were dependent on clinical parameters (large difference in crude OR and adjusted OR) and that, to a certain extent, circulating miRNA levels are a reflection of clinical parameters. Since the other studies did not correct for clinical parameters, heterogeneity in baseline characteristics of included patients might explain why miRNAs were not reproduced across studies. This is especially likely to have occurred if the heterogeneity across studies involves clinical characteristics relating to response, such as baseline DAS28, SJC, GC use and CRP, which were indeed the parameters that differed across the studies considered. In particular, patients included in the work by Castro-Villegas et al.  were more frequently treated with GCs (64.7 % versus 30.0 %), whereas those in the cohort used by Krintel et al. , showed a higher median CRP (15 versus 6 mg/l), TJC (15 versus 7), SJC (10 versus 1) and VAS-GH (70 versus 60) as compared to our cohort. In both previous studies the prediction model with miRNAs was not compared to or combined with baseline clinical parameters, which would have made the predictive estimations more generalizable. Concluding, heterogeneity cannot explain why results could not be validated within our study, although it might explain to a certain extent why the predictors identified across the studies are different.
Another possible contributing factor is represented by the chosen study design, in terms of inclusion criteria, measurement and time point of response, statistical analyses, etc. Krintel et al. , analyzed a cohort of TNFi-naïve patients treated with MTX and intra-articular triamcinolone, which were additionally randomized to ADA treatment (n = 90) or placebo treatment (n = 90). To identify miRNAs specifically predictive for EULAR response to ADA combination therapy, an interaction term for each miRNA with the received treatment was added to the prediction model. Castro-Villegas et al.  used a cohort of patients treated with ADA, ETN or infliximab and investigated how serum miRNAs changed over time in ten patients after TNFi initiation using miScript miRNA PCR array (Qiagen, Hilden, Germany). In the following step, the ten most relevant miRNAs were measured by single-miRNA assays in 85 additional patients, and univariate and multivariable tests were applied to predict response to any TNFi. In case the true predictive ability of miRNAs is weak, the design might make the difference in inclusion or exclusion of each miRNA, which would explain the differences in identified miRNAs across studies. In addition, discrepancy in found miRNAs across studies might arise from the fact that each TNFi treatment is analyzed independently (our study and Krintel et al. ) or in combination with others to find universal miRNAs for TNFi response (Castro-Villegas et al. ). Indeed, despite that all registered TNF-alpha- inhibiting therapies target TNF-alpha, they have small chemical differences and etanercept, in particular, also targets lymphotoxin-alpha . It is therefore possible that biomarkers predictive of response to TNFi therapy are to some extent TNFi-specific.
Despite the advantage in terms of stability, the identification of circulating miRNAs with a concrete potential of application in clinical practice is very limited. In other inflammatory diseases, such as inflammatory bowel disease, the usage of miRNAs as potential biomarkers is still being explored, though has so far not revealed usable predictors of response to therapy . In this line, the usefulness and robustness of miRNAs as biomarkers has been questioned; e.g., in non-neoplastic diseases only 33 % (139/416) of the reported miRNAs were considered either biologically plausible, specific for the disease or interpretable with the current knowledge . Another study showed that up to 58 % of the reported circulating miRNAs related to cancer subtypes were not disease-specific and most likely derived from blood cells . These studies indicate that false positive results in studies exploring circulating miRNA are lurking, and warrant additional carefulness when proposing a miRNA as a marker for a specific disease or disease state. For RA specifically, biomarkers are frequently identified in the circulation [13–19]. However, because the disease primarily affects the joints, the synovial compartment might constitute an alternative good source of biomarkers for RA, as has been demonstrated for cell-derived microparticles in the synovium compared to the circulation [36–38]. Since in RA the synovial miRNAs do not necessarily correlate with plasma miRNAs , it is possible that synovial miRNAs have better predictive abilities than the circulating ones. If so, incorporation of any predictive test on the synovial fluid or synovial tissue instead of the circulation, will affect the clinical feasibility negatively. In the future, the discovery of potential biomarker can be boosted by the implementation of novel high-throughput techniques. One of the most promising at this regard is next-generation sequencing (NGS) that has the potential to also identify novel and not previously annotated miRNAs (currently 1882 known (http://www.mirbase.org, accessed July 4, 2016 ), as it is not restricted to a predefined selection of miRNAs such as the multiplex-based techniques used in this study . NGS might identify previously unknown targets and discover novel miRNAs for the prediction of response to RA.
So far, there are no miRNAs that can be used in the prediction of response to TNFi therapy. We believe that a combination of differences in study design, technical variability, lack of multiple testing corrections, and heterogeneity between studies could contribute to these discrepancies. However, it is also conceivable that the irreproducibility of results is caused by the absence of truly biologically relevant miRNAs in the prediction of response to TNFi. Overall, our study demonstrated that in order to increase reproducibility of the results, future studies addressing this topic should (1) standardize detection methods, (2) investigate the added value of miRNAs over clinical parameters, (3) technically replicate findings using a method suitable in case of implementation in clinical practice (i.e., single assay), (4) validate findings in a separate cohort, especially when correction for multiple testing in the discovery phase is not performed, and (5) since heterogeneity influences the ORs of miRNAs, prediction models should preferably be validated in a cohort close to the target population.
Area under the receiver operating characteristic curve
- B&H FDR:
Benjamini and Hochberg false discovery rate
Biologicals and Outcome Compared and predicted Utrecht region in Rheumatoid Arthritis
relative threshold cycle
Disease activity score based on a 28-joint count
Erythrocyte sedimentation rate
European League Against Rheumatism
real-time quantitative PCR
swollen joint count
Wolfe AM, Kellgren JH, Masi AT. The epidemiology of rheumatoid arthritis: a review. II. Incidence and diagnostic criteria. Bull Rheum Dis. 1968;19:524–9.
Silman AJ, Pearson JE. Epidemiology and genetics of rheumatoid arthritis. Arthritis Res. 2002;4 Suppl 3:S265–72.
Keystone EC, Kavanaugh AF, Sharp JT, Tannenbaum H, Hua Y, Teoh LS, et al. Radiographic, clinical, and functional outcomes of treatment with adalimumab (a human anti-tumor necrosis factor monoclonal antibody) in patients with active rheumatoid arthritis receiving concomitant methotrexate therapy: a randomized, placebo-controlled, 52-week trial. Arthritis Rheum. 2004;50:1400–11.
Weinblatt ME, Kremer JM, Bankhurst AD, Bulpitt KJ, Fleischmann RM, Fox RI, et al. A trial of etanercept, a recombinant tumor necrosis factor receptor:Fc fusion protein, in patients with rheumatoid arthritis receiving methotrexate. N Engl J Med. 1999;340:253–9.
Cuppen BV, Welsing PM, Sprengers JJ, Bijlsma JW, Marijnissen AC, van Laar JM, et al. Personalized biological treatment for rheumatoid arthritis: a systematic review with a focus on clinical applicability. Rheumatology (Oxford). 2016;55:826–39.
Mookherjee N, El-Gabalawy HS. High degree of correlation between whole blood and PBMC expression levels of miR-155 and miR-146a in healthy controls and rheumatoid arthritis patients. J Immunol Methods. 2013;400–401:106–10.
Sayed D, Abdellatif M. MicroRNAs in development and disease. Physiol Rev. 2011;91:827–87.
Long L, Yu P, Liu Y, Wang S, Li R, Shi J, et al. Upregulated microRNA-155 expression in peripheral blood mononuclear cells and fibroblast-like synoviocytes in rheumatoid arthritis. Clin Dev Immunol. 2013;2013:296139.
Wang J, Chen J, Sen S. MicroRNA as biomarkers and diagnostics. J Cell Physiol. 2016;231:25–30.
Cortez MA, Bueso-Ramos C, Ferdin J, Lopez-Berestein G, Sood AK, Calin GA. MicroRNAs in body fluids--the mix of hormones and biomarkers. Nat Rev Clin Oncol. 2011;8:467–77.
Pauley KM, Cha S, Chan EK. MicroRNA in autoimmunity and autoimmune diseases. J Autoimmun. 2009;32:189–94.
Murata K, Yoshitomi H, Tanida S, Ishikawa M, Nishitani K, Ito H, et al. Plasma and synovial fluid microRNAs as potential biomarkers of rheumatoid arthritis and osteoarthritis. Arthritis Res Ther. 2010;12:R86.
Duroux-Richard I, Jorgensen C, Apparailly F. What do microRNAs mean for rheumatoid arthritis? Arthritis Rheum. 2012;64:11–20.
Luo X, Ranade K, Talker R, Jallal B, Shen N, Yao Y. microRNA-mediated regulation of innate immune response in rheumatic diseases. Arthritis Res Ther. 2013;15:210.
Wittmann J, Jack HM. microRNAs in rheumatoid arthritis: midget RNAs with a giant impact. Ann Rheum Dis. 2011;70 Suppl 1:i92–6.
Bottini N, Firestein GS. Epigenetics in rheumatoid arthritis: a primer for rheumatologists. Curr Rheumatol Rep. 2013;15:372.
Ceribelli A, Yao B, Dominguez-Gutierrez PR, Nahid MA, Satoh M, Chan EK. MicroRNAs in systemic rheumatic diseases. Arthritis Res Ther. 2011;13:229.
Salehi E, Eftekhari R, Oraei M, Gharib A, Bidad K. MicroRNAs in rheumatoid arthritis. Clin Rheumatol. 2015;34:615–28.
Churov AV, Oleinik EK, Knip M. MicroRNAs in rheumatoid arthritis: altered expression and diagnostic potential. Autoimmun Rev. 2015;14:1029–37.
Krintel SB, Dehlendorff C, Hetland ML, Horslev-Petersen K, Andersen KK, Junker P, et al. Prediction of treatment response to adalimumab: a double-blind placebo-controlled study of circulating microRNA in patients with early rheumatoid arthritis. Pharmacogenomics J. 2015;16:141–6.
Castro-Villegas C, Perez-Sanchez C, Escudero A, Filipescu I, Verdu M, Ruiz-Limon P, et al. Circulating miRNAs as potential biomarkers of therapy effectiveness in rheumatoid arthritis patients treated with anti-TNFalpha. Arthritis Res Ther. 2015;17:49.
Duroux-Richard I, Pers YM, Fabre S, Ammari M, Baeten D, Cartron G, et al. Circulating miRNA-125b is a potential biomarker predicting response to rituximab in rheumatoid arthritis. Mediators Inflamm. 2014;2014:342524.
Prevoo ML, van 't Hof MA, Kuper HH, van Leeuwen MA, van de Putte LB, van Riel PL. Modified disease activity scores that include twenty-eight-joint counts. Development and validation in a prospective longitudinal study of patients with rheumatoid arthritis. Arthritis Rheum. 1995;38:44–8.
van Gestel AM, Haagsma CJ, van Riel PL. Validation of rheumatoid arthritis improvement criteria that include simplified joint counts. Arthritis Rheum. 1998;41:1845–50.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–8.
Youden W. Index for rating diagnostic tests. Cancer. 1950;3:32–5.
Ma R, Jiang T, Kang X. Circulating microRNAs in cancer: origin, function and application. J Exp Clin Cancer Res. 2012;31:38.
Witwer KW. Circulating microRNA biomarker studies: pitfalls and potential solutions. Clin Chem. 2015;61:56–63.
Farr RJ, Januszewski AS, Joglekar MV, Liang H, McAulley AK, Hewitt AW, et al. A comparative analysis of high-throughput platforms for validation of a circulating microRNA signature in diabetic retinopathy. Sci Rep. 2015;5:10375.
Fujioka S, Nakamichi I, Esaki M, Asano K, Matsumoto T, Kitazono T. Serum microRNA levels in patients with Crohn’s disease during induction therapy by infliximab. J Gastroenterol Hepatol. 2014;29:1207–14.
Nadal E, Truini A, Nakata A, Lin J, Reddy RM, Chang AC, et al. A novel serum 4-microRNA signature for lung cancer detection. Sci Rep. 2015;5:12464.
Tansey MG, Szymkowski DE. The TNF superfamily in. 2009: new pathways, new indications, and new drugs. Drug Discov Today. 2009;14:1082–8.
Kalla R, Ventham NT, Kennedy NA. MicroRNAs: new players in inflammatory bowel disease. Gut. 2015;64:1008.
Haider BA, Baras AS, McCall MN, Hertel JA, Cornish TC, Halushka MK. A critical evaluation of microRNA biomarkers in non-neoplastic disease. PLoS One. 2014;9, e89565.
Pritchard CC, Kroh E, Wood B, Arroyo JD, Dougherty KJ, Miyaji MM, et al. Blood cell origin of circulating microRNAs: a cautionary note for cancer biomarker studies. Cancer Prev Res (Phila). 2012;5:492–7.
Knijff-Dutmer EA, Koerts J, Nieuwland R, Kalsbeek-Batenburg EM, van de Laar MA. Elevated levels of platelet microparticles are associated with disease activity in rheumatoid arthritis. Arthritis Rheum. 2002;46:1498–503.
Vinuela-Berni V, Doniz-Padilla L, Figueroa-Vega N, Portillo-Salazar H, Abud-Mendoza C, Baranda L, et al. Proportions of several types of plasma and urine microparticles are increased in patients with rheumatoid arthritis with active disease. Clin Exp Immunol. 2015;180:442–51.
Berckmans RJ, Nieuwland R, Kraan MC, Schaap MC, Pots D, Smeets TJ, et al. Synovial microparticles from arthritic patients modulate chemokine and cytokine release by synoviocytes. Arthritis Res Ther. 2005;7:R536–44.
Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42:D68–73.
Pimentel F, Bonilla P, Ravishankar YG, Contag A, Gopal N, LaCour S, et al. Technology in microRNA profiling: circulating microRNAs as noninvasive cancer biomarkers in breast cancer. J Lab Autom. 2015;20:574–88.
We thank Katja Coeleveld, Kim van der Wurff-Jacobs and Dorien van de Berg, for bio-banking, Karin A.L. Schrijvers-te-Brake, Annemiek Sloeserwij, Joke Nijdeken and Marieke Vianen who facilitated in the acquisition of data, Cornelis Bekker for his technical assistance, Paco Welsing for his aid in statistical analyses, and the Society for Rheumatology Research Utrecht (SRU) for including patients. We are grateful to all of the patients who have participated in this study.
Contributing SRU hospitals and their contact persons: University Medical Center Utrecht, Dr. J. Tekstra; Antonius Hospital Nieuwegein/Utrecht, Dr. E.J. ter Borg; Diakonessen Hospital Utrecht, Dr. Y. Schenk; Meander Medical Center Amersfoort, Dr. S. Linn-Rasker; Sint Maartenskliniek Woerden, Dr. W.H. van der Laan; Hospital St. Jansdal Harderwijk, Dr. D.G. Kuiper-Geertsma; Tergooi Hospital Hilversum, Dr. M. Nabibux; Flevo Hospital Almere, Dr. C.M. Verhoef.
This research was supported ZonMW: The Netherlands Association for Health Research and Development (project no. 152001008), Foundation Friends of the UMCU based on a private gift and unrestricted grants from Pfizer, BMS, UCB, Roche, and Schering-Plough. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The datasets supporting the conclusions of this article are not publically available though can be sent on request.
BC collected clinical data, performed experiments, performed statistical analyses, and drafted the manuscript. MR coordinated the experiments, aided in statistical analyses and interpretation of the results and helped to draft and revise the manuscript. RF-S included patients, coordinated the study, and helped to draft and revise the manuscript. ANC performed single assay experiments and helped to draft the manuscript. YS included patients and helped to draft the manuscript. JWJB conceived and coordinated the study and helped to draft the manuscript. TRDJR conceived the study and helped to draft the manuscript. FPJGL conceived and coordinated the study and helped to draft the manuscript. All authors read and approved the final manuscript.
JWJB received research grants (to his department) and consultancy fees from AbbVie, BMS, Crescendo, MSD, Mundipharma, Pfizer, Roche, Sun, and UCB. JvL received grants from Roche, MSD, Eli Lilly and Pfizer. The other authors report no conflict of interests.
Consent for publication
Ethics approval and consent to participate
The study was approved by the local ethics committee of the University Medical Center Utrecht and the institutional review boards of the participating centers (see Acknowledgements). Informed consent was obtained from each patient.
Baseline characteristics of responders and non-responders, split for cohort and treatment. (DOCX 22 kb)
Baseline characteristics of patients in the discovery and validation cohort, split for received treatment. (DOCX 19 kb)
Expression of miRNAs predicting TNFi response identified in other studies. In this analysis, the predictive ability of miR-22, miR-23a, miR-223 and miR-886-3p is explored. (DOCX 217 kb)
Validation of multivariable models for prediction of response to TNFi. The prediction models from the discovery cohort were applied in the validation cohort using the miRNA levels as measured by single assays (DOCX 15 kb)
Technical replication of selected miRNAs from the discovery cohort. Using single miRNA assays, the four selected miRNAs were retested in the same patients for their predictive abilities. (DOCX 104 kb)
Clinical parameters at baseline associated with current miRNA levels. Serum values of the four selected miRNAs were tested for their relationship with clinical parameters, independent from response. (DOCX 16 kb)
Cell subsets related to miRNA expression level, as measured by OpenArray or qPCR assay. In this additional analysis, miRNA expression levels are compared with FACS results on PBMC cell subsets. (XLSX 63 kb)
About this article
Cite this article
Cuppen, B.V.J., Rossato, M., Fritsch-Stork, R.D.E. et al. Can baseline serum microRNAs predict response to TNF-alpha inhibitors in rheumatoid arthritis?. Arthritis Res Ther 18, 189 (2016). https://doi.org/10.1186/s13075-016-1085-z
- Rheumatoid arthritis
- TNF-alpha inhibitor