Type I interferon response gene expression in established rheumatoid arthritis is not associated with clinical parameters
Arthritis Research & Therapy volume 18, Article number: 290 (2016)
A peripheral blood interferon (IFN) signature (i.e., elevated type I interferon response gene [IRG] expression) has been described in a subset of patients with rheumatoid arthritis (RA). In the present study, we systematically assessed the association between this IRG expression and clinical parameters.
Expression of 19 IRGs was determined in peripheral blood from 182 consecutive patients with RA and averaged into an IFN score per individual. Correlation and unpaired analyses were performed on the complete patient group. The analyses were internally validated by using an algorithm to randomize the patient group 1000 times into two equally sized sets, and then analyses were performed on both sets.
Associations were assessed between IFN score and disease duration, 28-joint Disease Activity Score and its components, the occurrence of erosions and nodules, autoantibody positivity, and immunosuppressive treatment. This analysis revealed lower IFN scores in patients using hydroxychloroquine, prednisone, and/or sulfasalazine, but it did not show significant associations between the other parameters and the IFN score. Selecting patients who were not treated with hydroxychloroquine, prednisone, and/or sulfasalazine (n = 95) did not reveal any significant associations either.
IRG expression in RA is affected by immunosuppressive treatment with prednisone, hydroxychloroquine, and/or sulfasalazine, but it is not evidently associated with other clinical parameters. Hence, the IFN signature appears to describe a subgroup of patients with RA but does not seem to reflect disease activity.
Rheumatoid arthritis (RA) is a systemic autoimmune disease characterized by chronic joint inflammation. It manifests as a heterogeneous disease that is partly reflected at the level of gene expression. Genome-wide gene expression analysis revealed evidence for molecular differences between patients with RA, in particular in the type I interferon (IFN) system. Approximately 50% of patients with RA display a peripheral blood IFN signature (i.e., relatively high expression of interferon response genes [IRGs]) .
Type I IFNs were initially known for their antiviral effects, but increasing insight into their activities revealed their role as pleiotropic cytokines with a critical role in modulating immune responses, such as cellular activation, major histocompatibility complex upregulation, induction of apoptosis, and inhibition of angiogenesis . It is thought that type I IFNs contribute to autoimmunity by initiating a break of tolerance (e.g., by the induction of dendritic cell maturation and inhibition of regulatory T cells) . The exact role of the IFN signature in RA is yet unknown, although it was shown to have potential clinical relevance. That is, (1) the presence of the IFN signature was shown to be a risk factor for arthritis development in preclinical disease , and (2) the presence of the IFN signature in established RA was found to be associated with the clinical response to treatment with rituximab  and tocilizumab .
Earlier studies have addressed whether the IFN signature in RA could be associated with clinical parameters, which inconclusively revealed a potential relationship of the IFN signature with anticitrullinated protein antibody (ACPA) titers [7, 8]. However, these study cohorts were rather small (35 subjects or less) and therefore highly subject to a lack of power. Hence, the relationships between the peripheral blood IFN signature and disease- and inflammation-related clinical parameters have never been thoroughly assessed. In the present study, we used a larger cohort of patients with established RA (n = 182) in combination with a random sampling algorithm to systematically investigate whether the peripheral blood IFN signature in RA could be associated with parameters such as disease activity, laboratory parameters, and the use of immunosuppressive treatment.
Patient recruitment and blood collection
Patients with RA (n = 182) were consecutively recruited from the Jan van Breemen Research Institute, Reade center, Amsterdam, The Netherlands. All patients were Caucasian and were diagnosed with RA according to the American College of Rheumatology 1987 criteria . Ninety-five percent of the patients displayed a 28-joint Disease Activity Score (DAS28) of ≥3.2 despite earlier treatment with at least two disease-modifying antirheumatic drugs (DMARDs). All patients provided written informed consent, and this study was approved by the medical ethics committee of Reade. The patients’ characteristics are shown in Table 1. From all patients, 2.5 ml of blood was drawn into PAXgene tubes (PreAnalytiX GmbH, Hombrechtikon, Switzerland) and stored at −20 °C until further processing.
RNA isolation, complementary DNA synthesis, and real-time polymerase chain reaction
Total RNA was isolated from the PAXgene tubes according to the manufacturer’s instructions. The quantity and purity of the RNA samples were checked using a NanoDrop spectrophotometer (Thermo Scientific, Wilmington, DE, USA). We reverse-transcribed 0.25 μg of RNA into complementary DNA (cDNA) using a RevertAid H Minus cDNA Synthesis Kit (Thermo Fisher, Waltham, MA, USA). A single aliquot of each cDNA sample was first subjected to 14 cycles of specific target amplification using a 0.2× mixture of all TaqMan gene expression assays in combination with the TaqMan PreAmp Master Mix (Applied Biosystems, Foster City, CA, USA). Following preamplification, the samples were diluted 1:5 (vol/vol) in Tris-ethylenediaminetetraacetic acid buffer, pH 8.0. Multiplex real-time quantitative polymerase chain reaction was performed using the 96.96 Biomark Dynamic Array systems (Fluidigm Corporation, South San Francisco, CA, USA) at ServiceXS (Leiden, The Netherlands) according to the manufacturer’s instructions. Quantities were calculated relative to GAPDH using the standard curve method. Expression levels were log2-transformed.
Calculation of the IFN score and statistical analyses
Nineteen IRGs described to be components of the IFN signature in RA  were measured (see Additional file 1: Table S1). Log2-transformed expression levels of the IRGs were highly correlated (r ≥ 0.7 for 90% of the combinations, p ≤ 0.002); therefore, an IFN score was calculated by averaging these values of all genes for each sample.
Data were analyzed using IBM SPSS Statistics version 22 (IBM, Armonk, NY, USA), R version 3.1.3 , and GraphPad Prism version 5.01 (GraphPad Software, La Jolla, CA, USA) software. For internal cross-validation, a 1000-times random sampling method was used to randomize the group of 182 patients into 2 equally sized sets and to execute Spearman correlation for continuous variables and Mann-Whitney U analysis for dichotomous variables on each set . p Values <0.05 were considered to be significant. Correction for multiple testing was performed using the method of Benjamini and Hochberg.
We studied the association between the peripheral blood IFN score and the following parameters: disease duration, DAS28 and its individual components, the occurrence of erosions and nodules, autoantibody positivity, and immunosuppressive treatment. As demonstrated in Table 2, prednisone use and dose, hydroxychloroquine (HCQ) use, and sulfasalazine (SSZ) use were the only variables that showed a significant result, of which only HCQ use remained significant after correction for multiple testing. Similar results were obtained in the cross-validation using the random sampling algorithm, which revealed significance for HCQ use only: a p value <0.05 was detected in both sets in 521 of the 1000 iterations, in 1 of the 2 sets for 479 of 1000 iterations, and never in none of the sets (median p value 0.015). A slight trend was also observed for prednisone use and dose and SSZ use (median p values 0.090–0.14, median coefficient for prednisone dose −0.18), although significance was never found in both sets for these variables. Each IRG was also analyzed individually, which revealed similar results (data not shown).
For prednisone, HCQ, and SSZ treatment, the IFN score was lower in the treated group than in the untreated group. Combining HCQ use, prednisone use, and SSZ use also revealed a significantly lower IFN score in patients using HCQ and/or prednisone and/or SSZ than in patients not treated with any of these agents, both in the analysis of the complete cohort and in the cross-validation (p = 0.0080 with Benjamini-Hochberg correction, median cross-validation p ≤ 0.012) (see Table 2 and Fig. 1a). Moreover, the suppressive effect appeared larger for patients treated with two or more of those agents than for patients treated with one agent (Fig. 1b). No association was found between IFN score and methotrexate (MTX) treatment or dose. The relation between IFN score and treatment did not appear to be confounded by DAS28 or disease duration.
Because the suppression of IFN score in HCQ-, prednisone-, and/or SSZ-treated patients could have a masking effect on other associations between IFN signature and clinical parameters, we also performed the analyses for the selection of patients who were not treated with prednisone, HCQ, and/or SSZ (n = 95) at the moment of blood collection. This did not result in any significant associations between the IFN score and other variables (uncorrected p value ≥0.099, median p value cross-validation ≥0.23) (see Additional file 1: Table S2).
The present study is the first use of a systematic approach in a relatively large cohort to study the relationship between the IFN signature in established RA and clinical parameters. We demonstrated that the IFN signature was suppressed in patients treated with HCQ, prednisone, and/or SSZ, but not with MTX. Furthermore, we did not observe any associations between the IFN signature and the other clinical parameters.
Van der Pouw Kraan et al. showed that a subgroup of patients with RA displays a common pathogen-response program, which was characterized by a higher incidence of the IFN signature as well as higher ACPA titers, suggesting that these parameters might be associated with one another . However, a causal relationship was not established, and our data indicate that this is not the case. The IFN signature was not significantly different between ACPA-negative and ACPA-positive patients, nor did it significantly correlate with ACPA titers. Possibly, the IFN signature and ACPA positivity are independently associated with activation of the common pathogen response program, because they are both implied to be induced via certain pathogens [12, 13].
Our cohort consisted mainly of established patients with RA with a DAS28 ≥ 3.2 despite treatment with at least two DMARDs. Although our data do not suggest any association between IFN score and DAS28 (p > 0.1 in all analyses), we cannot fully exclude the possibility that significant associations could have been found if the cohort had included more patients with low disease activity. Inclusion of patients with inactive and/or early disease could provide a more complete view of the IFN response in different states of the disease.
Remarkably, the IFN scores were decreased in HCQ-, prednisone-, and/or SSZ-treated patients, even though the beneficial effects of these treatments were supposedly diminished. Moreover, cotreatment with these agents appeared to have an additive suppressive effect. Interference of both prednisone and HCQ with type I IFN signaling has been described before [14, 15], but the influence of SSZ remains to be elucidated. It has been shown that SSZ reduces the levels of RA-related cytokines, such as interleukin-1β and tumor necrosis factor-α , suggesting that SSZ might function through overall suppression of inflammatory cytokines, including type I IFNs. Furthermore, it was demonstrated that SSZ is able to accelerate apoptosis of neutrophils , which we have recently shown to be major inducers of the type I IFN response in RA . Consequently, suppression of the type I IFN response by SSZ might be mediated via an increase in neutrophil apoptosis.
As previously described, suppression of the IFN score by certain treatment could affect the applicability of the IFN signature as a biomarker for therapy response, particularly to rituximab [5, 19]. That is, the treatment-related suppression of IFN score might impair the discriminative capacity of the biomarker, which would consequently lead to more false predictions. Because multiple studies have demonstrated that the extent of the IFN signature is highly variable between patients [1, 20, 21], we considered it important to assess the IFN score as a continuous variable rather than as a dichotomous variable. However, inclusion of a healthy control population would allow us to determine whether patients with prednisone, SSZ, and/or HCQ treatment would still display an IFN response above the levels of healthy control subjects. This could also give more insight into the applicability of the IFN response as a biomarker in these patients. Future studies should be done to elucidate the effect of each individual treatment, as well as combinatory therapy, on the IFN signature and the corresponding response prediction. Alternatively, presence of the IFN signature in individuals with arthralgia was shown to be associated with a higher risk for developing arthritis . It would be interesting to study whether early treatment with one of the implied suppressors of the IFN response could delay or even prevent disease onset.
Our data indicate that there are no evident associations between the peripheral blood IFN signature in established RA and clinical parameters. This suggests that the IFN signature is not an indication of disease activity per se, but its presence could indicate a potential difference in pathology or immune pathway activation compared with patients without this signature. Consequently, this could influence the response to therapy, particularly to biologics because these are specific modulators of these immune pathways.
Anticitrullinated protein antibodies
28-joint Disease Activity Score
Disease-modifying antirheumatic drug
Erythrocyte sedimentation rate
Immunoglobulin M rheumatoid factor
Interferon response gene
Swollen joint count in 28 joints
Tender joint count in 28 joints
Visual analogue scale
van der Pouw Kraan TC, Wijbrandts CA, van Baarsen LG, Voskuyl AE, Rustenburg F, Baggen JM, et al. Rheumatoid arthritis subtypes identified by genomic profiling of peripheral blood cells: assignment of a type I interferon signature in a subpopulation of patients. Ann Rheum Dis. 2007;66:1008–14.
Hervas-Stubbs S, Perez-Gracia JL, Rouzaut A, Sanmamed MF, Le BA, Melero I. Direct effects of type I interferons on cells of the immune system. Clin Cancer Res. 2011;17:2619–27.
Bacher N, Graulich E, Jonuleit H, Grabbe S, Steinbrink K. Interferon-alpha abrogates tolerance induction by human tolerogenic dendritic cells. PLoS One. 2011;6, e22763.
Lubbers J, Brink M, van de Stadt LA, Vosslamber S, Wesseling JG, van Schaardenburg D, et al. The type I IFN signature as a biomarker of preclinical rheumatoid arthritis. Ann Rheum Dis. 2013;72:776–80.
Raterman HG, Vosslamber S, de Ridder S, Nurmohamed MT, Lems WF, Boers M, et al. The interferon type I signature towards prediction of non-response to rituximab in rheumatoid arthritis patients. Arthritis Res Ther. 2012;14:R95.
Sanayama Y, Ikeda K, Saito Y, Kagami S, Yamagata M, Furuta S, et al. Prediction of therapeutic responses to tocilizumab in patients with rheumatoid arthritis: biomarkers identified by analysis of gene expression in peripheral blood mononuclear cells using genome-wide DNA microarray. Arthritis Rheumatol. 2014;66:1421–31.
van der Pouw Kraan TC, van Baarsen LG, Wijbrandts CA, Voskuyl AE, Rustenburg F, Baggen JM, et al. Expression of a pathogen-response program in peripheral blood cells defines a subgroup of rheumatoid arthritis patients. Genes Immun. 2008;9:16–22.
Cantaert T, van Baarsen LG, Wijbrandts CA, Thurlings RM, van de Sande MG, Bos C, et al. Type I interferons have no major influence on humoral autoimmunity in rheumatoid arthritis. Rheumatology (Oxford). 2010;49:156–66.
Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31:315–24.
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2014.
van de Wiel MA, Berkhof J, van Wieringen WN. Testing the prediction error difference between 2 predictors. Biostatistics. 2009;10:550–60.
Jenner RG, Young RA. Insights into host responses against pathogens from transcriptional profiling. Nat Rev Microbiol. 2005;3:281–94.
Valesini G, Gerardi MC, Iannuccelli C, Pacucci VA, Pendolino M, Shoenfeld Y. Citrullination and autoimmunity. Autoimmun Rev. 2015;14:490–7.
Flammer JR, Dobrovolna J, Kennedy MA, Chinenov Y, Glass CK, Ivashkiv LB, et al. The type I interferon signaling pathway is a target for glucocorticoid inhibition. Mol Cell Biol. 2010;30:4564–74.
Kuznik A, Bencina M, Svajger U, Jeras M, Rozman B, Jerala R. Mechanism of endosomal TLR inhibition by antimalarial drugs and imidazoquinolines. J Immunol. 2011;186:4794–804.
Barrera P, Boerbooms AM, van de Putte LB, van der Meer JW. Effects of antirheumatic agents on cytokines. Semin Arthritis Rheum. 1996;25:234–53.
Bertolotto M, Dallegri F, Dapino P, Quercioli A, Pende A, Ottonello L, et al. Sulphasalazine accelerates apoptosis in neutrophils exposed to immune complex: role of caspase pathway. Clin Exp Pharmacol Physiol. 2009;36:1132–5.
de Jong TD, Lubbers J, Turk S, Vosslamber S, Mantel E, Bontkes HJ, et al. The type I interferon signature in leukocyte subsets from peripheral blood of patients with early arthritis: a major contribution by granulocytes. Arthritis Res Ther. 2016;18:165.
de Jong TD, Vosslamber S, Blits M, Wolbink G, Nurmohamed MT, van der Laken CJ, et al. Effect of prednisone on type I interferon signature in rheumatoid arthritis: consequences for response prediction to rituximab. Arthritis Res Ther. 2015;17:78.
Higgs BW, Liu Z, White B, Zhu W, White WI, Morehouse C, et al. Patients with systemic lupus erythematosus, myositis, rheumatoid arthritis and scleroderma share activation of a common type I interferon pathway. Ann Rheum Dis. 2011;70:2029–36.
de Jong TD, Vosslamber S, Mantel E, de Ridder S, Wesseling JG, van der Pouw Kraan TCTM, et al. Physiological evidence for diversification of IFNα- and IFNβ-mediated response programs in different autoimmune diseases. Arthritis Res Ther. 2016;18:49.
This research was performed with support from the Dutch Arthritis Foundation (project number LLP20).
Availability of data and material
The datasets analyzed during the present study are available from the corresponding author on reasonable request.
TDdJ was involved in study conception and design, analysis and interpretation of the data, and drafting of the manuscript. MB was involved in the acquisition of patient material and data and participated in the data analysis and interpretation. SdR was involved in the data analysis and interpretation. SV participated in the acquisition, analysis, and interpretation of the data. MTN and GW participated in the acquisition of patient material and in data interpretation. CV was involved in study conception and design and in the drafting of the manuscript. TDdJ, MB, and CLV had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis, which were agreed to by all authors. All authors were involved in critical revision of the manuscript for important intellectual content, and all authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All patients provided written informed consent, and this study was approved by the medical ethics committee of Reade.
Cornelis L. Verweij is deceased.
List of the 19 IFN response genes measured. Table S2. Results of analysis for associations between IFN scores and clinical parameters after 1000 times of random sampling using only patients who were not treated with prednisone, HCQ, or SSZ. (PDF 82 kb)
About this article
Cite this article
de Jong, T.D., Blits, M., de Ridder, S. et al. Type I interferon response gene expression in established rheumatoid arthritis is not associated with clinical parameters. Arthritis Res Ther 18, 290 (2016). https://doi.org/10.1186/s13075-016-1191-y