Skip to main content

Type I interferon in patients with systemic autoimmune rheumatic disease is associated with haematological abnormalities and specific autoantibody profiles

Abstract

Objectives

To investigate the relationships between interferon alpha (IFNα) and the clinical and serological phenotype of patients with systemic autoimmune rheumatic disease (SARDs) in order to determine whether a distinct subpopulation of patients can be identified.

Methods

We recruited patients with at least 1 SARD clinical feature and at least 1 SARD-related autoantibody from two NHS Trusts in Greater Manchester. A 6-gene interferon-stimulated gene (ISG) score was calculated in all patients, and in a subgroup, a 30-gene ISG score was produced using NanoString. A digital Single Molecule Array (Simoa) was used to measure plasma IFNα protein. In an exploratory analysis, whole blood RNA sequencing was conducted in 12 patients followed by RT-qPCR confirmation of expression of 6 nucleic acid receptors (NARs) in the whole cohort.

Results

Sixty three of 164 (38%) patients had a positive ISG score. The 3 measures of IFNα all correlated strongly with each other (p < 0.0001). There were no differences in mucocutaneous or internal organ involvement between the ISG subgroups. The ISG-positive group had increased frequency of specific autoantibodies and haematological abnormalities which remained significant after adjusting for the SARD subtype. Expression of DDX58, MB21D1 and TLR7 was correlated with the ISG score whilst TLR3, TLR9 and MB21D1 were associated with neutrophil count.

Conclusion

In SARD patients, IFNα-positivity was associated with specific autoantibodies and haematological parameters but not with other clinical features. The variable NAR expression suggests that different pathways may drive IFNα production in individual patients. The identification of an IFNα-positive subgroup within a mixed SARD cohort supports a pathology-based approach to treatment.

Introduction

Systemic Autoimmune Rheumatic Diseases (SARDs) are a group of multisystem autoimmune disorders with overlapping clinical and serological manifestations. Patients may transition between clinical disease categories over time; some evolve from an undifferentiated connective tissue disease (UCTD) to a specific disease type, whilst others develop an overlap syndrome. The significant commonality in clinical features and autoantibody profiles between some conditions suggests a shared molecular aetiology, whilst within a single-disease group significant molecular heterogeneity has been described [1]. Studies seeking to understand disease pathogenesis, identify reliable biomarkers, or develop novel treatment strategies may therefore be improved if undertaken in a molecularly defined subgroup of SARD patients, rather than in traditionally classified individual conditions [2]. The molecular taxonomy of SARDs offers new opportunities in clinical trial design; therefore, it is essential to better understand whether true immunopathological subsets of patients exists within clinical SARD cohorts.

Previous single-disease studies have demonstrated an increase in type I interferon (IFN) in some SARDs (notably Systemic Lupus Erythematosus (SLE) and Sjogren’s syndrome (SS)) [3,4,5]. However, it is less clear whether those SARD patients expressing a high level of IFNα represent a homogenous molecular subgroup.

The type I IFNs are a family of 17 cytokines compromising predominantly of IFNα subtypes. Thirteen IFNα subtypes have been identified in humans, with the coding genes all located on chromosome 9 [6]. Studies to date have been challenging due to the relatively low abundance of IFNα in the circulation, and the lack of adequately specific antibodies to confidently distinguish between the IFNα subtypes and between IFNα and other type I and type III IFNs (IFNλ). The ‘gold standard’ has therefore been the measurement of interferon-stimulated genes (ISGs) either directly in human samples or via an ex vivo stimulation assay [7, 8]. More recently, a highly sensitive digital ELISA has been trialled in adult and paediatric patients with SLE and juvenile dermatomyositis [9, 10].

We aimed to investigate type I IFN pathway activity using multiple methodologies in an unselected adult SARD cohort and to determine how this related to the clinical and immunological phenotype of patients. We also used unbiased transcriptome analysis of a subgroup of patients and investigated relevant IFN pathway molecules in order to move beyond classical physician diagnoses and towards a molecular taxonomy of SARDs.

Methods

Patients and healthy volunteers

Patients were recruited into the Lupus Extended Autoimmune Phenotype (LEAP) cohort from Manchester University NHS Foundation Trust (3 sites) and Salford Royal NHS Foundation Trust (1 site). The LEAP cohort includes patients with SLE, SS, UCTD, mixed connective tissue disease (MCTD), systemic sclerosis (SSc) and idiopathic inflammatory myopathies (IIM). Ethical approval was obtained from the Greater Manchester East Research Ethics Committee (13/NW/0564), and the study was conducted in accordance with the Declaration of Helsinki. Patients with an established diagnosis and clinically stable disease were eligible for inclusion if they had ≥ 1 clinical feature of a connective tissue disease and ≥ 1 relevant autoantibody. Autoantibodies were measured according to physician discretion but all included the BioPlex 2200 ANA Screen with MDSS (comprising the following autoantigens: dsDNA, Chromatin/nucelosome, Ro/SS-A, La/SS-B, Smith [Sm], SmRNP, RNP, Scl-70, Jo-1 and centromere). The physician (rheumatologist) diagnosis was used as the primary classifier of patients. Detailed clinical and serological data were recorded. Whole blood was collected for subsequent RNA analysis. Routine clinical biochemical and haematological parameters were measured, and additional plasma samples were collected for IFNα protein measurement (see below).

Healthy volunteers were recruited to provide comparator cohorts. For measurement of the ISG score, 29 healthy controls (16/29, 55% female), with a median age 24 years, were recruited (this cohort has been as previously described [11]). For the NanoString analysis, 27 of the 29 volunteers above were sampled: 16/27 (59%) female, median age 34 years. For the nucleic acid receptors analysis, 16 healthy controls (14/16 (88%) female) with a median age of 42 years were used.

ISGs

Blood was collected into PAXgene tubes (PreAnalytix) kept at room temperature for between 2 and 72 h, frozen at − 20 °C, then stored at − 80 °C until extraction. Total RNA was extracted from whole blood using a PAXgene (PreAnalytix) RNA isolation kit. RNA concentration was assessed by Qubit 2.0 fluorometer (Thermo Fisher Scientific). After DNAse treatment, RNA was converted into cDNA with the high-capacity cDNA reverse transcriptase kit (Qiagen). Quantitative reverse transcription polymerase chain reaction (RT-qPCR) analysis was performed using the TaqMan Universal PCR Master Mix (Applied Biosystems), and cDNA derived from 40 ng total RNA. Using TaqMan probes for IFI27 (Hs01086370_m1), IFI44L (Hs00199115_m1), IFIT1 (Hs00356631_g1), ISG15 (Hs00192713_m1), RSAD2 (Hs01057264_m1) and SIGLEC1 (Hs00988063_m1), the relative abundance of each target transcript was normalised to the expression level of HPRT1 (Hs03929096_g1) and 18S (Hs999999001_s1) and assessed with the Applied Biosystems StepOne Software v2.1 and Data Assist Software v.3.01. For each of the six probes, individual (patient and control) data were expressed relative to a single calibrator. The median fold change of the six ISGs, when compared to the median of the combined data of healthy controls, was used to create an interferon score for each patient as previously described [11]. RQ is equal to 2−ΔΔCt, i.e. the normalised fold change relative to a control. The mean interferon score of the healthy controls plus two standard deviations above the mean (+ 2SD) was calculated. Scores above this value (> 2·466) were designated as positive.

Nucleic acid receptors (NARs)

RNA was extracted from whole blood, and RT-qPCR analysis was performed as described above using 18S and HPRT1 as reference probes. TaqMan probes for DDX58 (Hs01061436_m1), MB21D1 (Hs00403553_m1), TMEM173 (Hs00736956_m1), TLR3 (Hs01551079_g1), TLR7 (Hs01933259_s1) and TLR9 (Hs00370913_s1) were used. For each of the six probes, individual (patient and control) data were expressed relative to a single calibrator.

Quantification of interferon alpha (IFNα) protein levels in plasma using a single molecule array (Simoa)

As described by Rodero et al., the Simoa IFNα assay was developed using a Quanterix Homebrew Simoa assay according to the manufacturer’s instructions and utilising two autoantibodies specific for IFNα isolated and cloned from 2 APS1/APECED patients [9, 12]. The 8H1 antibody clone was used as a capture antibody after coating on paramagnetic beads (0.3 mg/mL), and the 12H5 was biotinylated (biotin/Ab ratio = 30/1) and used as the detector. Recombinant IFNα17/αI (PBL Assay Science) was used as a standard curve after cross-reactivity testing. The limits of detection (LOD) were calculated by the mean value of all blank runs + 3SDs and was 0.23 fg/ml.

NanoString ISG analysis

Gene expression of 30 interferon-stimulated genes and 3 housekeeping genes (HPRT1, NRDC, OTUD5, for details see Additional file 1) was measured. A 30 gene score was calculated based on the same methods for the RT-qPCR score (+ 2SD of healthy controls). Scores above 1.642 were considered positive. All ISG, NAR and NanoString measurements were conducted on single RNA samples.

Gene expression analysis by RNA sequencing (RNA-Seq)

Whole transcriptome expression analysis was performed using samples from 12 participants, selected according to the ISG score and presence/absence of anti-Smith antibodies (see Additional file 1: Supplementary methods for details). Gene lists were uploaded into Ingenuity Pathway Analysis (http://www.ingenuity.com) in order to determine differentially regulated canonical pathways in the patients. ISGs were identified within the RNA-Seq dataset by comparing differentially expressed genes with the Interferome v2.0 database [13] (accessed from http://www.interferome.org/interferome/home.jspx on 8/6/2018).

Availability of data

The RNA-Seq dataset supporting the conclusions of this article is available in the Array Express repository, (E-MTAB-7080, https://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-7080).

Statistical analysis

Statistics were calculated using STATA v.13 (STATcorp, USA), GraphPad Prism version 5.0d for Mac OS X and version 7.00 for Windows (San Diego, CA, USA) and R v3.4.1. Non-parametric statistical tests were used (Mann-Whitney U and Kruskal-Wallis for comparison of 2 or > 2 groups, respectively, Spearman r for correlations) unless specified otherwise. Univariate linear and logistic regression models were used. Multivariable models all included age, gender and ethnicity plus additional relevant variables as described.

Results

ISG expression in patients with SARDs

We recruited a total of 164 patients with a median (IQR) age of 48.5 (36.8, 57.3) years and disease duration of 7.11 (3.16, 15.8) years (Table 1). The majority (122/163, 74.8%) of patients were Caucasian reflecting a clinical population in North West England and 155 (94.5%) were female. The most prevalent conditions in our cohort were SLE (67 patients, 40.9%) and UCTD (43 patients, 26.2%) (Table 1).

Table 1 Demographic characteristics of the patient population

In the whole cohort the median (IQR) ISG score was 0.8115 (0.3155, 7.052). The log-transformed ISG score had a bimodal distribution (Fig. 1a). Using an ISG score cut-off of > 2.466 [11], 63 (38.4%) patients had a positive score. A finite mixture model was also used to group patients into IFN-high and IFN-low on the basis of the bimodal distribution and compared directly to our pre-validated cut-off. Using this approach, all IFN-low patients were classified as ISG score negative. Conversely, 9 patients were classified as IFN-high, but ISG score negative (see Additional file 1: Table S2).

Fig. 1
figure 1

Distribution of the 6-gene ISG score in the SARD cohort. a The log-transformed 6-gene ISG score shows a bimodal distribution. b The 6-gene ISG score correlates both with measurement of IFN protein and with a more extensive 30-gene score measured using NanoString. Spearman’s r is shown

In a subset of 92 patients, the IFNα plasma protein level was also measured using Simoa as previously described [14]. The median (IQR) IFNα protein concentration was 3.178 (0.361, 26.3) fg/ml. Using a cut-off of 10 fg/ml [9], 38/93 (40.9%) patients had a positive score. There was an excellent correlation between IFNα protein levels and the 6 gene RT-qPCR ISG score (0.8307, p < 0.0001). The 30 gene NanoString panel was measured in a subgroup of 30 patients, and there were excellent correlations between NanoString and the 6-gene RT-qPCR score (r = 0.9226, p < 0.0001) and between NanoString and IFNα protein (r = 0.7807, p < 0.0001) (Fig. 1b).

The number of patients with a positive ISG score varied between disease groups (Fisher’s exact test p < 0.001) (Fig. 2a and Table 1). Patients with SS and MCTD were most likely to have a positive ISG score (14/20, [70%] and 9/13 [69.2%], respectively) whilst positive ISG scores were less frequent in patients with SSc and UCTD (2/13 [15.4%] and 9/43 [20.9%], respectively). The distribution of ISG scores was similar when the relevant classification criteria were applied in lieu of the physician diagnosis (see Additional file 1: Figure S1).

Fig. 2
figure 2

The ISG score is increased in a subset of SARD patients. a The number of patients with a positive ISG score varies between disease groups. Each data point represents a single subject according to their clinical diagnosis. The horizontal line shows the 95th centile for healthy subjects (ISG score of 2.466). b The ISG score is increased in patients with a greater number of autoantibodies. The total number of autoantibodies (excluding ANA) detected in each patient is shown. Comparisons are made against patients with specific autoantibodies (Dunn’s multiple comparison’s test), **p < 0.01, ****p < 0.0001. c The graph shows the predicted probabilities of a positive ISG score according to whether the patients have SLE (blue) or another SARD (red). These were obtained using logistic regression models adjusted for age, gender, ethnicity and concomitant mediation. The points and error bars show the mean and standard deviation of predicted probabilities

In univariate logistic regression models, ISG score positivity was significantly associated with MCTD (OR [95% CI] 8.50 [2.12, 34.0]), SLE (2.55 [1.06, 6.16]) and SS (8.81 [2.64, 29.4]). In a multivariable model which also included medication (concomitant use (Y/N) of steroids, immunosuppressants and anti-malarials), these significant associations remained. We also noted an inverse association with age (0.97, [0.94, 0.998]), Caucasian ethnicity (0.38 [0.16, 0.91]) and anti-malarial use (0.41 [0.16, 1.04]).

Higher ISG expression is associated with haematological abnormalities in SARDs

In the whole cohort, we found a higher frequency of haematological abnormalities, as defined by the 1997 Modified ACR Classification Criteria for SLE [15] in patients with a positive ISG score (39/61 [63.9]% vs. 35/99 [35.4]%, p < 0.0001). In addition, patients with a positive ISG score demonstrated lower haemoglobin concentration and lower total white cell count (WCC), lymphocyte and neutrophil counts (Table 2). The platelet count was also slightly lower in ISG score-positive patients, although this was not statistically significant. There was no association between ISG score positivity and other common SARD clinical features.

Table 2 Association between ISG score and haematological parameters

In multivariable linear regression models which also included disease group, number of autoantibodies and medication use, the total WCC, lymphocyte count and neutrophil count all remained significantly inversely associated with the ISG score (see Additional file 1: Table S3).

Autoantibodies against RNA proteins are associated with a positive ISG score

The association between the autoantibody profiles of patients was determined using logistic regression models. In unadjusted univariate models, the odds of being in the ISG score-positive group were significantly increased in patients with anti-Ro, anti-La, anti-Sm, anti-RNP and anti-chromatin antibodies and rheumatoid factor (Table 3). There were significant differences in the frequency of anti-Ro, anti-La, anti-Smith, anti-RNP, anti-chromatin, anti-dsDNA and SSc-specific antibodies between the clinical disease groups. In a fully-adjusted model (including disease group) anti-Smith, anti-chromatin and rheumatoid factor all remained significantly associated with the ISG score.

Table 3 Association between autoantibodies and ISG score using logistic regression models

In patients with no specific autoantibodies (beyond ANA), the median (IQR) ISG score was 0.407 (0.303, 0.998) and there was a significant increase in ISG score as the total number of autoantibodies in each patient (defined as anti-dsDNA, anti-Smith, anti-RNP, anti-Ro, anti-La, anti-CCP, anti-chromatin/nucleosome or SSc-specific antibodies [anti-centromere, anti-Scl-70 or RNA polymerase III]) increased (Kruskal-Wallis, p = 0.0001) (Fig. 2b). In a multivariable model including disease group and medication use, each additional autoantibody increased the odd of having a positive ISG score by over 2-fold (unadjusted OR 2.44 [176, 3.36], adjusted OR 2.14 [1.49, 30.5]). The association between autoantibody number and ISG score was also determined for SLE patients alone compared to the rest of the SARD cohort. In both cases, a significant association was observed (OR 1.87 [1.18, 2.96] and OR 4.62 [2.35, 90.6] respectively, Fig. 2c).

ISG score-positive patients also had a higher median (IQR) ESR (erythrocyte sedimentation rate) (14.0 [8, 35.5] vs. 8.0 [5, 17.5] mm/h, p = 0.007) but not CRP (serum C-reactive protein) concentration (2.0 [1.0, 4.0] vs. 1.5 [1.0, 5.0] mg/dl, p = 0.7196). There was no association between a positive ISG score and history of low C3 or low C4 complement.

Whole blood RNA-Seq analysis

An exploratory RNA-Seq analysis was performed on samples from 12 SLE participants from 4 ‘cohorts’: 3 ISG positive/anti-Sm positive; 3 ISG positive/anti-Sm negative; 3 ISG-negative/anti-Sm-positive participants and 3 ISG-negative/anti-Sm-negative. One ISG-negative/anti-Sm-positive patient failed QC and was removed from the analysis. Anti-Sm was selected as it was most strongly related to the ISG score, allowing analysis of patients who were both ISG and Sm positive (‘double positive’), and patients with anti-Sm without a positive ISG. Hierarchical clustering of all coding transcripts was conducted (after removal of any globin/HLA genes). This demonstrated some similarities between patients belonging to the groupings described (Fig. 3a). Strikingly, of the 100 genes showing the greatest variance across all 11 samples, only a small number were ISGs. In addition, although the 3 ISG positive/anti-Sm positive patients shared high expression of the ISGs, 1 of the 3 patients (lane 1) showed markedly different expression levels of the other genes.

Fig. 3
figure 3

Whole blood transcriptome analysis of ISG-positive and ISG-negative SLE patients. a Heatmap of the 100 genes with the greatest variance. Upregulated genes are shown in red; downregulated genes are shown in blue. The patient samples are hierarchically clustered (Euclidean distance) over all coding genes. The heatmap is annotated to show known ISGs (dark green). Each patient sample is annotated according to the ISG score (positive or negative) and the anti-Sm antibody status (positive or negative). b Gene ontology analysis showing the canonical pathways that are over-represented in the ISG score-positive patients compared to the negative patients. c The expression of nucleic acid receptors (NARs) within the RNA-Seq dataset between the ISG score-positive (n = 6) and ISG score-negative (n = 5) patients

Gene ontology analysis of over-represented canonical biological pathways, between the ISG positive (n = 6) and negative (n = 5) participants, identified differential expression of nucleic acid receptors (NARs). The GO terms with the greatest enrichment were (i) interferon signalling, (ii) activation of IRF by cytosolic pattern recognition receptors and (iii) role of pattern recognition receptors in recognition of bacteria and viruses (Fig. 3b). The expression of a number of NARs was increased in ISG score-positive patients, and a group of these were then taken forward for further study (Fig. 3c).

Expression of nucleic acid receptors (NARs) differs between disease groups and is associated with the ISG score

The expression of 6 NARs (TLR3, TLR7, TLR9, DDX58 [RIG-I], MB21D1 [c-GAS] and TMEM173 [STING]) was measured in 155 patients and 16 HC (healthy controls were 14/16 [87.5]% female with a median age of 42 years). Expression of all NARs (with the exception of TLR3) was increased in SARDs compared to HC (Fig. 4a). There was a strong correlation between the ISG score and DDX58 expression (r = 0.7386, p < 0.0001) and modest correlations between the ISG score and TLR7 and MB21D1 (r = 0.3262, p < 0.0001 and r = 0.2266, p = 0.005, respectively) (Fig. 4b). Using logistic regression DDX58 and TLR7 were associated with a positive ISG score (OR 3.23 [2.25, 4.67] and 1.68 [1.14, 2.46] with DDX58 remaining significant in a multivariable model containing disease group (OR 3.69, [2.32, 5.85])).

Fig. 4
figure 4

Nucleic acid receptors are differentially expressed in patients with SARDs. a CTD patients have increased expression of all of the NARs except TLR3. Comparisons were made using Dunn’s multiple comparison’s test. *p < 0.05, ***p < 0.001, ****p < 0.0001. RQ = relative quantification of gene expression. b Heatmap showing the relative expression of each of the 6 NARs. Each row represents a single patient and is K-means clustered into 6 groups. The rows are annotated by diagnostic group and by ISG score (ISG score-positive patients in green). c Correlation between the NAR expression and ISG score in the combined CTD cohort. The graphs show the Spearman r for each NAR

Expression of all of the NARs (except TLR3) was higher in patients compared to healthy controls according to disease phenotype (see Additional file 1: Figure S2). In linear regression models of log-adjusted NAR expression, patients with anti-Ro, anti-RNP, anti-Smith and anti-chromatin antibodies had significantly increased DDX58 expression (see Additional file 1: Table S4). These associations remained significant after adjustment for age and gender, but were no longer significant when the ISG score was included in the model (data not shown). Patients with anti-Ro antibodies also had significantly lower expression of TLR9 and MB21D1 which remained significant in a multivariable model that included ISG score (β − 0.173 [− 0.233, − 0.015] and − 0.190 [− 0.346, − 0.034].

In an exploratory analysis, we identified clusters of patients based on the expression of the 6 NARs (Fig. 4c). Of those patients with a positive ISG score and increased DDX58 expression, a subgroup also had increased expression of TLR3. In contrast, patients with increased TMEM173 also had increased expression of either TLR7 or TLR9, but not both.

TLR3, TLR9 and MB21D1 (c-GAS) expression are independently associated with neutrophil counts

Significant correlations between a number of the NARs and haematological parameters (haemoglobin, lymphocyte count, neutrophil count and platelets) were identified (Fig. 5). Interestingly, both positive and negative correlations were observed for neutrophils. Linear regression models identified that TLR3 was negatively associated, whilst TLR9 and MB21D1 were positively associated, with neutrophil counts even after adjusted for prednisolone use and the ISG score (Table 4).

Fig. 5
figure 5

Correlation between the ISG score, NARs and haematological parameters. The figure shows a Spearman correlation matrix for the ISG score, NAR expression and haematological parameters. Only statistically significant (p < 0.05) parameters are shown. The size of the dot represents the strength of correlation (Spearman’s r); blue = positive correlation, red = negative correlation. Hb = haemoglobin

Table 4 Association between NAR expression and neutrophil count

Discussion

In this mixed SARD cohort, we used multiple methodologies to identify subgroups of patients with shared immunopathological phenotypes despite different clinical disease manifestations. Interestingly, in our multivariable regression models, type I IFN pathway activity was more strongly associated with specific autoantibodies and haematological parameters than clinical disease subtypes or individual clinical features.

Since type I interferon was first reported to be present in the serum of SLE patients in 1975 [16], a number of studies have examined type I IFN expression in patients with SARDs [17]. A panel of downstream ISGs, termed the ‘interferon signature’, is commonly used as a surrogate of IFNα expression. In this study, we used our established 6-probe RT-qPCR ISG assay [11] which has been extensively validated in patients with monogenic interferonopathies [18]. Although our patient and healthy control groups varied in terms of age, we have previously reported that our ISG score is not associated with age [11]. We also employed NanoString technology assessing 30 probes of interest. The significant correlation between these measures demonstrates that either method may be used to assess ISGs, with NanoString potentially offering a less operator-dependent technique, and thus greater potential for widespread clinical use. We also demonstrated a strong correlation with IFNα protein using a high-sensitivity digital ELISA, which could detect IFNα protein in the fg/ml range. This gives confidence that our six-gene panel is measuring the downstream effects of in vivo IFNα activity. In 6 (6.5%) patients, we observed a difference between the ISG score and digital ELISA and these patients are the subject of further study. Similarly to El-Sherbiny et al. [19], we also observed that the ISG score has a clear bimodal distribution. Although we utilised our pre-defined threshold for the ISG score, this approach and the results of the finite mixture model were highly concordant.

In our cohort, increased IFN expression was identified in 38% of patients, with a subpopulation in each clinical disease group demonstrating elevated levels. The overall number of ISG score-positive patients is lower than that observed in other autoimmune populations, perhaps reflecting the fact that we report an established mixed SARD cohort (including patients with UCTD), or that the patients had clinically stable disease, and that our cohort was predominantly Caucasian. We confirmed the findings of others that ISG score-positive patients were younger, but that there was no association with disease duration [10]. In a large study by Weckerle et al., IFNα activity was significantly increased in patients of non-European ancestry [20]. The relationship between IFNα and disease activity is less clear; in patients with SLE, some groups have reported an association [21, 22] whilst others have not [23, 24]. Most recently, a study of the whole blood ISG signature in SLE patients showed remarkable stability over time with no association with changes in disease activity [25]. This suggests that the IFN signature may be used to identify a phenotypic subpopulation of SARD patients in cross-sectional studies.

It has been proposed that increased interferon expression may pre-date the development of SARDs. In a study by Wither et al., increased IFN expression was identified in asymptomatic patients with a positive ANA [26]. Furthermore, a positive IFN score may predict those patients who go on to develop SLE or SS. In a study of 118 patients considered to be at risk for developing a SARDs, the IFN score was greater in those who progressed to a SARD than those who did not [27]. We also identified an ISG-positive subgroup in our UCTD patients with established disease.

In our mixed SARD cohort, we observed that this type I IFN signature was associated with haematological indices, but not with other autoimmune features. An inverse association between lymphocyte count and IFN score has been reported previously [19], but in our cohort, we also found associations with neutrophils, haemoglobin and a trend towards a lower platelet count. In the general population, viral infections with a corresponding IFN response are associated with lymphopenia [28]. Similarly, chronic IFNα exposure results in reduced haemoglobin, neutrophils and lymphocytes in patients with chronic hepatitis C, supporting a direct role for IFNα in the haematological abnormalities observed in our cohort [29].

The proportion of patients with haemolytic anaemia was very small, and we only observed a trend towards lower platelet counts; this is in contrast to patients with monogenic type I interferonopathies in whom thrombocytopenia is often the most common haematological manifestation [30, 31]. In addition, iatrogenic thrombocytopenia has been reported following interferon-β therapy [32]. The mechanism by which IFN may result in these haematological abnormalities is yet to be fully elucidated; however, it is likely that a number of mechanisms are involved including increased lymphocyte migration to tissue and reduced egress from lymph nodes [33].

In this larger clinical cohort, we corroborated our previous findings of a positive association between elevated type I interferon and the RNA-associated autoantibodies anti-Ro, anti-La, anti-Sm and anti-RNP [9]. Additionally, we observed an association between anti-chromatin/nucleosome antibodies and rheumatoid factor and a positive ISG score which appeared to be independent of the clinical diagnosis. Importantly, anti-chromatin antibodies were associated with increased IFNα whilst anti-dsDNA antibodies were not, suggesting that the DNA-histone complex may be more interferonogenic than DNA alone. This observation is in contrast with the study by Kirou et al. which found an association between IFNα and anti-dsDNA titres in SLE patients, although this may have been confounded by increased disease activity [21]. RNA-containing immune complexes appear to directly increase IFNα production in plasmacytoid dendritic cells (pDCs) in a mechanism that is dependent on both toll-like receptors (TLRs) and the Fc gamma receptor IIa [34], but may also trigger IFNα indirectly via effects on neutrophils [35].

We observed increased expression of all of the nuclear receptors, except TLR3, in whole blood of patients with SARDs compared to healthy controls. This supports the study by Chauhan et al. which identified increased PBMC TLR7 and TLR9 expression in SLE patients compared to healthy controls [36]. However, in our cohort, DDX58/RIG-I, MB21D1 and TLR7 expression were associated with the ISG score which may reflect our measurements in whole blood which does not allow adjustment for immune cell heterogeneity. Maria et al. reported that TLR7 was increased and TLR9 decreased in the pDCs and monocytes of IFN-positive SS patients [37]. In addition, DDX58/RIG-I was also increased in monocytes from IFN-positive patients. Similarly in childhood SLE, increased expression of TLR7 and a number of RIG-like receptors (including DDX58) has been reported in IFN score-positive patients [38].

Unsupervised hierarchical clustering of NAR expression identified groups of ISG-score-negative patients with high TLR3 expression, or high TLR9 and/or TMEM173/STING expression. The potential roles of these receptors in IFN-low subgroups warrant further investigation. Interestingly, although TLR3 was not increased in SARD patients compared to controls, there was a significant independent negative association with neutrophil count (in opposition to MB21D1 and TLR9) suggesting that non-IFN mechanisms might be involved. In contrast, TLR3 ligation in mice has been demonstrated to result in increased neutrophil numbers suggesting that increased TLR expression and increased TLR3 signalling may not explain this association [39].

Conclusions

In summary, transcriptomic and protein studies demonstrate that for subpopulations of the SARDs there are overlapping IFN-related biomarkers which may represent shared molecular drivers of disease and are more associated with specific disease features and ENA than clinical diagnoses. We also noted that NAR expression differs between ISG-positive and ISG-negative patients which may point to different immunopathological mechanisms between ISG “immunophenotypes”. Our observations suggest that the molecular taxonomy of a mixed SARD cohort is inconsistent with the clinical phenotype. A targeted approach to treatment of patients with a shared immunopathogenesis may therefore offer advantages over the current classification of disease.

References

  1. Banchereau R, Hong S, Cantarel B, Baldwin N, Baisch J, Edens M, et al. Personalized immunomonitoring uncovers molecular networks that stratify lupus patients. Cell. 2016;165(6):1548–50.

    Article  CAS  Google Scholar 

  2. Barturen G, Beretta L, Cervera R, Van VR, Alarcon-Riquelme ME. Moving towards a molecular taxonomy of autoimmune rheumatic diseases. Nat Rev Rheumatol. 2018;14(2):75–93.

    Article  CAS  Google Scholar 

  3. Ronnblom L, Alm GV, Eloranta ML. The type I interferon system in the development of lupus. Semin Immunol. 2011;23(2):113–21.

    Article  Google Scholar 

  4. 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(11):2029–36.

    Article  CAS  Google Scholar 

  5. Brkic Z, Maria NI, van Helden-Meeuwsen CG, van de Merwe JP, van Daele PL, Dalm VA, et al. Prevalence of interferon type I signature in CD14 monocytes of patients with Sjogren's syndrome and association with disease activity and BAFF gene expression. Ann Rheum Dis. 2013;72(5):728–35.

    Article  CAS  Google Scholar 

  6. Manry J, Laval G, Patin E, Fornarino S, Itan Y, Fumagalli M, et al. Evolutionary genetic dissection of human interferons. J Exp Med. 2011;208(13):2747–59.

    Article  CAS  Google Scholar 

  7. Rice GI, Del Toro DY, Jenkinson EM, Forte GM, Anderson BH, Ariaudo G, et al. Gain-of-function mutations in IFIH1 cause a spectrum of human disease phenotypes associated with upregulated type I interferon signaling. Nat Genet. 2014;46(5):503–9.

    Article  CAS  Google Scholar 

  8. Hua J, Kirou K, Lee C, Crow MK. Functional assay of type I interferon in systemic lupus erythematosus plasma and association with anti-RNA binding protein autoantibodies. Arthritis Rheum. 2006;54(6):1906–16.

    Article  CAS  Google Scholar 

  9. Rodero MP, Decalf J, Bondet V, Hunt D, Rice GI, Werneke S, et al. Detection of interferon alpha protein reveals differential levels and cellular sources in disease. J Exp Med. 2017;214(5):1547–55.

    Article  CAS  Google Scholar 

  10. Mathian A, Mouries-Martin S, Dorgham K, Devilliers H, Barnabei L, Ben SE, et al. Monitoring disease activity in systemic lupus erythematosus with single-molecule array digital ELISA quantification of serum interferon-alpha. Arthritis Rheumatol. 2019;71(5):756–765.

    Article  CAS  Google Scholar 

  11. Rice GI, Forte GM, Szynkiewicz M, Chase DS, Aeby A, Abdel-Hamid MS, et al. Assessment of interferon-related biomarkers in Aicardi-Goutieres syndrome associated with mutations in TREX1, RNASEH2A, RNASEH2B, RNASEH2C, SAMHD1, and ADAR: a case-control study. Lancet Neurol. 2013;12(12):1159–69.

    Article  CAS  Google Scholar 

  12. Meyer S, Woodward M, Hertel C, Vlaicu P, Haque Y, Karner J, et al. AIRE-deficient patients harbor unique high-affinity disease-ameliorating autoantibodies. Cell. 2016;166(3):582–95.

    Article  CAS  Google Scholar 

  13. Rusinova I, Forster S, Yu S, Kannan A, Masse M, Cumming H, et al. Interferome v2.0: an updated database of annotated interferon-regulated genes. Nucleic Acids Res. 2013;41(Database issue):D1040–6.

    CAS  PubMed  Google Scholar 

  14. Llibre A, Bondet V, Rodero MP, Hunt D, Crow YJ, Duffy D. Development and validation of an ultrasensitive single molecule array digital enzyme-linked immunosorbent assay for human interferon-alpha. J Vis Exp. 2018;136. https://doi.org/10.3791/57421.

  15. Hochberg MC. Updating the American college of rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum. 1997;40(9):1725.

    Article  CAS  Google Scholar 

  16. Skurkovich SV, Eremkina EI. The probable role of interferon in allergy. Ann Allergy. 1975;35(6):356–60.

    CAS  PubMed  Google Scholar 

  17. Ronnblom L, Eloranta ML. The interferon signature in autoimmune diseases. Curr Opin Rheumatol. 2013;25(2):248–53.

    Article  Google Scholar 

  18. Rice GI, Melki I, Fremond ML, Briggs TA, Rodero MP, Kitabayashi N, et al. Assessment of type I interferon signaling in pediatric inflammatory disease. J Clin Immunol. 2017;37(2):123–32.

    Article  CAS  Google Scholar 

  19. El-Sherbiny YM, Psarras A, Yusof MYM, Hensor EMA, Tooze R, Doody G, et al. A novel two-score system for interferon status segregates autoimmune diseases and correlates with clinical features. Sci Rep. 2018;8(1):5793.

    Article  CAS  Google Scholar 

  20. Weckerle CE, Franek BS, Kelly JA, Kumabe M, Mikolaitis RA, Green SL, et al. Network analysis of associations between serum interferon-alpha activity, autoantibodies, and clinical features in systemic lupus erythematosus. Arthritis Rheum. 2011;63(4):1044–53.

    Article  CAS  Google Scholar 

  21. Kirou KA, Lee C, George S, Louca K, Peterson MG, Crow MK. Activation of the interferon-alpha pathway identifies a subgroup of systemic lupus erythematosus patients with distinct serologic features and active disease. Arthritis Rheum. 2005;52(5):1491–503.

    Article  CAS  Google Scholar 

  22. Rose T, Grutzkau A, Hirseland H, Huscher D, Dahnrich C, Dzionek A, et al. IFNalpha and its response proteins, IP-10 and SIGLEC-1, are biomarkers of disease activity in systemic lupus erythematosus. Ann Rheum Dis. 2013;72(10):1639–45.

    Article  CAS  Google Scholar 

  23. Kennedy WP, Maciuca R, Wolslegel K, Tew W, Abbas AR, Chaivorapol C, et al. Association of the interferon signature metric with serological disease manifestations but not global activity scores in multiple cohorts of patients with SLE. Lupus Sci Med. 2015;2(1):e000080.

    Article  Google Scholar 

  24. Landolt-Marticorena C, Bonventi G, Lubovich A, Ferguson C, Unnithan T, Su J, et al. Lack of association between the interferon-alpha signature and longitudinal changes in disease activity in systemic lupus erythematosus. Ann Rheum Dis. 2009;68(9):1440–6.

    Article  CAS  Google Scholar 

  25. Petri M, Fu W, Ranger A, Allaire N, Cullen P, Magder LS, et al. Association between changes in gene signatures expression and disease activity among patients with systemic lupus erythematosus. BMC Med Genet. 2019;12(1):4.

    Google Scholar 

  26. Wither J, Johnson SR, Liu T, Noamani B, Bonilla D, Lisnevskaia L, et al. Presence of an interferon signature in individuals who are anti-nuclear antibody positive lacking a systemic autoimmune rheumatic disease diagnosis. Arthritis Res Ther. 2017;19(1):41.

    Article  Google Scholar 

  27. Md Yusof MY, Psarras A, El-Sherbiny YM, Hensor EMA, Dutton K, Ul-Hassan S, et al. Prediction of autoimmune connective tissue disease in an at-risk cohort: prognostic value of a novel two-score system for interferon status. Ann Rheum Dis. 2018;77:1432–9.

    PubMed  PubMed Central  Google Scholar 

  28. Schattner A, Meshorer A, Wallach D. Involvement of interferon in virus-induced lymphopenia. Cell Immunol. 1983;79(1):11–25.

    Article  CAS  Google Scholar 

  29. Schmid M, Kreil A, Jessner W, Homoncik M, Datz C, Gangl A, et al. Suppression of haematopoiesis during therapy of chronic hepatitis C with different interferon alpha mono and combination therapy regimens. Gut. 2005;54(7):1014–20.

    Article  CAS  Google Scholar 

  30. Crow YJ, Manel N. Aicardi-Goutieres syndrome and the type I interferonopathies. Nat Rev Immunol. 2015;15(7):429–40.

    Article  CAS  Google Scholar 

  31. Briggs TA, Rice GI, Adib N, Ades L, Barete S, Baskar K, et al. Spondyloenchondrodysplasia due to mutations in ACP5: a comprehensive survey. J Clin Immunol. 2016;36(3):220–34.

    Article  CAS  Google Scholar 

  32. Koudriavtseva T, Plantone D, Renna R, Mandoj C, Giannarelli D, Mainero C. Interferon-beta therapy and risk of thrombocytopenia in multiple sclerosis patients. Neurol Sci. 2015;36(12):2263–8.

    Article  Google Scholar 

  33. Kamphuis E, Junt T, Waibler Z, Forster R, Kalinke U. Type I interferons directly regulate lymphocyte recirculation and cause transient blood lymphopenia. Blood. 2006;108(10):3253–61.

    Article  CAS  Google Scholar 

  34. Lovgren T, Eloranta ML, Kastner B, Wahren-Herlenius M, Alm GV, Ronnblom L. Induction of interferon-alpha by immune complexes or liposomes containing systemic lupus erythematosus autoantigen- and Sjogren's syndrome autoantigen-associated RNA. Arthritis Rheum. 2006;54(6):1917–27.

    Article  Google Scholar 

  35. Lande R, Ganguly D, Facchinetti V, Frasca L, Conrad C, Gregorio J, et al. Neutrophils activate plasmacytoid dendritic cells by releasing self-DNA-peptide complexes in systemic lupus erythematosus. Sci Transl Med. 2011;3(73):73ra19.

    Article  Google Scholar 

  36. Chauhan SK, Singh VV, Rai R, Rai M, Rai G. Distinct autoantibody profiles in systemic lupus erythematosus patients are selectively associated with TLR7 and TLR9 upregulation. J Clin Immunol. 2013;33(5):954–64.

    Article  CAS  Google Scholar 

  37. Maria NI, Steenwijk EC, IJpma AS, van Helden-Meeuwsen CG, Vogelsang P, Beumer W, et al. Contrasting expression pattern of RNA-sensing receptors TLR7, RIG-I and MDA5 in interferon-positive and interferon-negative patients with primary Sjogren's syndrome. Ann Rheum Dis. 2017;76(4):721–30.

    Article  CAS  Google Scholar 

  38. Wahadat MJ, Bodewes ILA, Maria NI, van Helden-Meeuwsen CG, van Dijk-Hummelman A, Steenwijk EC, et al. Type I IFN signature in childhood-onset systemic lupus erythematosus: a conspiracy of DNA- and RNA-sensing receptors? Arthritis Res Ther. 2018;20(1):4.

    Article  Google Scholar 

  39. Downes JE, Marshall-Clarke S. Innate immune stimuli modulate bone marrow-derived dendritic cell production in vitro by toll-like receptor-dependent and -independent mechanisms. Immunology. 2010;131(4):513–24.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research was funded by the NIHR Manchester Biomedical Research Centre. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health. This work was also supported by the Centre for Epidemiology Versus Arthritis (grant number 21755). Professor Bruce is an NIHR Senior Investigator and is funded by Versus Arthritis, the National Institute for Health Research Manchester Biomedical Research Centre and the NIHR Manchester Clinical Research Facility. Funding is also acknowledged from the National Institute for Health Research (NIHR) (NIHR Transitional Research Fellowship, Dr. Tracy Briggs, TRF-2016-09-002). This study was also supported by the Medical Research Council (MR/N003322/1). DD acknowledges support from the ANR (CE17001002) and Immunoqure for provision of mAbs for Simoa.

Availability of data and materials

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

Author information

Authors and Affiliations

Authors

Contributions

The research project was conceived by IB, BP, JR, EMMc, CH, ALH and TB. Clinical data were obtained and reviewed by EB, SH, HC, ALH, BP and IB. Laboratory analyses were performed by DS, GR, DD VB and AH. Data analysis was principally performed by JR, TB, LZ, GI and DD although all authors contributed. All authors contributed to the writing of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ian N. Bruce.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was obtained from the Greater Manchester East Research Ethics Committee (13/NW/0564), and the study was conducted in accordance with the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional file

Additional file 1:

Supplementary methods (DOCX 1004 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Reynolds, J.A., Briggs, T.A., Rice, G.I. et al. Type I interferon in patients with systemic autoimmune rheumatic disease is associated with haematological abnormalities and specific autoantibody profiles. Arthritis Res Ther 21, 147 (2019). https://doi.org/10.1186/s13075-019-1929-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13075-019-1929-4

Keywords