Skip to main content

Correlation analyses of clinical and molecular findings identify candidate biological pathways in systemic juvenile idiopathic arthritis

Abstract

Background

Clinicians have long appreciated the distinct phenotype of systemic juvenile idiopathic arthritis (SJIA) compared to polyarticular juvenile idiopathic arthritis (POLY). We hypothesized that gene expression profiles of peripheral blood mononuclear cells (PBMC) from children with each disease would reveal distinct biological pathways when analyzed for significant associations with elevations in two markers of JIA activity, erythrocyte sedimentation rate (ESR) and number of affected joints (joint count, JC).

Methods

PBMC RNA from SJIA and POLY patients was profiled by kinetic PCR to analyze expression of 181 genes, selected for relevance to immune response pathways. Pearson correlation and Student's t-test analyses were performed to identify transcripts significantly associated with clinical parameters (ESR and JC) in SJIA or POLY samples. These transcripts were used to find related biological pathways.

Results

Combining Pearson and t-test analyses, we found 91 ESR-related and 92 JC-related genes in SJIA. For POLY, 20 ESR-related and 0 JC-related genes were found. Using Ingenuity Systems Pathways Analysis, we identified SJIA ESR-related and JC-related pathways. The two sets of pathways are strongly correlated. In contrast, there is a weaker correlation between SJIA and POLY ESR-related pathways. Notably, distinct biological processes were found to correlate with JC in samples from the earlier systemic plus arthritic phase (SAF) of SJIA compared to samples from the later arthritis-predominant phase (AF). Within the SJIA SAF group, IL-10 expression was related to JC, whereas lack of IL-4 appeared to characterize the chronic arthritis (AF) subgroup.

Conclusions

The strong correlation between pathways implicated in elevations of both ESR and JC in SJIA argues that the systemic and arthritic components of the disease are related mechanistically. Inflammatory pathways in SJIA are distinct from those in POLY course JIA, consistent with differences in clinically appreciated target organs. The limited number of ESR-related SJIA genes that also are associated with elevations of ESR in POLY implies that the SJIA associations are specific for SJIA, at least to some degree. The distinct pathways associated with arthritis in early and late SJIA raise the possibility that different immunobiology underlies arthritis over the course of SJIA.

Peer Review reports

Background

Systemic juvenile idiopathic arthritis (SJIA) is currently classified as a subtype of juvenile idiopathic arthritis [1], and is characterized by a combination of arthritis and systemic inflammation, including fever, rash and serositis. SJIA has distinct demographic characteristics compared to other JIA subtypes, including onset throughout childhood and lack of gender preference. At clinical presentation, SJIA may resemble other diseases in children, including viral infection and Kawasaki disease [2–4]. The outcome in SJIA is variable, with close to half of children having a monocyclic course, less than 10% having an intermittent course, and over half having a persistent course [5, 6], the latter often dominated by chronic arthritis. An adult form of SJIA is called Adult Onset Still Disease (AOSD) and occurs rarely [7].

There are also unique immunophenotypic features in SJIA compared to other JIA subtypes, such as the lack of human leukocyte antigen (HLA) class II allele association, low or absent autoantibodies (specifically, antinuclear antibodies, rheumatoid factor or anti-CCP antibodies [8]), a tendency toward monocytosis [9, 10], high levels of IL-18 [11, 12] and natural killer cell abnormalities in at least a subset of patients [13]. These immunologic features, together with the therapeutic efficacy of inhibitors of IL-1 or IL-6 in SJIA and AOSD, suggest that these diseases might be best classified as autoinflammatory rather than autoimmune [14–17].

Despite our knowledge of some important immunological characteristics of active SJIA, the pathogenesis of SJIA remains unknown. One of the unanswered questions is whether independent biological processes underlie the systemic symptoms and the arthritis. Evidence from clinical studies shows that earlier in the disease, IL-1 inhibitors (and perhaps also IL-6 blockade) are efficacious, especially against systemic symptoms, but at a later stage, where arthritis may predominate, patients may develop resistance to these therapies [18–20]. These findings suggest that distinct biological processes may be associated with different manifestations and/or different stages of the disease.

Transcriptional profiling of peripheral blood cells has been a useful approach for identifying biological pathways involved in SJIA and other complex diseases, such as polyarticular JIA (POLY), rheumatoid arthritis (RA), systemic lupus erythematosus and Kawasaki disease [21–24]. Previous studies of SJIA using microarray analyses have revealed transcriptional signatures in peripheral blood associated with active disease and with patient subsets [25–29].

We hypothesized that distinct gene expression patterns may be associated with individual clinical parameters used as measures of the systemic inflammation and the arthritis. We analyzed expression in peripheral blood mononuclear cells (PBMC) of a panel of inflammation-associated genes to determine patterns associated with elevations in two markers of disease activity in JIA, erythrocyte sedimentation rate (ESR) and number of active joints (joint count, JC). ESR is a marker of inflammation that is elevated in association with systemic as well as organ-specific inflammation, including arthritis [30]. Active joints are defined as joints with non-bony swelling or limited range of motion, with either tenderness or pain on motion; we chose active joint count as a marker of arthritis.

We asked if common or unique expression profiles are associated with ESR and JC in SJIA. In order to assess the specificity of our results for SJIA, we also asked whether the expression of the panel of tested genes differed in SJIA patients compared to patients with polyarticular course JIA (POLY), which is characterized by chronic polyarthritis. We then analyzed if JC associated genes differ during the early and late phase of SJIA. Based on the gene expression patterns, we identified candidate biological pathways associated with the systemic and arthritis components of SJIA.

Methods

Subject population and clinical data collection

The study was approved by the Stanford University Administrative Panel on Human Subjects in Medical Research (protocol ID 13932). Informed consent was obtained from patients or parents or guardians before blood sample collection. Venous blood samples from all subjects were treated anonymously throughout the analysis. All JIA patients were followed at the Pediatric Rheumatology Clinic at Lucile Packard Children's Hospital. SJIA and POLY patients met amended ILAR criteria for diagnosis [1]. Thirty-one SJIA and 18 POLY individual patients participated in this study. A total of 46 SJIA samples (22 Flare and 24 Quiescence samples), and 25 POLY samples (17 Flare and 8 Quiescence samples) were analyzed. Some patients, (SJIA n = 15; POLY n = 7) contributed samples during both flare and quiescent disease states. Twelve POLY patients were rheumatoid factor (RF) negative, and six were RF positive. All samples were classified as flare (F) or quiescence (Q) based on a scheme we developed for this and other studies of JIA [10, 31, 32] (Tables 1, 2 and 3 and [33, 34]). SJIA flare samples had a systemic score of ≥ 1 and/or an arthritis score of ≥ B (≥ 5 active joints). POLY flare samples had an arthritis score of ≥ 1 (≥ 1 to 10 active joints). Arthritis severity is scored differently for SJIA and POLY patients, because the patterns of joint involvement generally are different between the two groups [34], with the exception that some SJIA patients develop POLY-like arthritis with symmetric, small joint involvement. The arthritis scoring system is based on frequency analyses of numbers of active joints in early active SJIA [34] and in active POLY [Sandborg C, frequency data not shown]. Comprehensive clinical information was collected at each patient visit, including history, physical exam and clinical laboratory values [10]. As shown in Table 4, and consistent with the known demographics of JIA [35], SJIA patients are younger than POLY patients and are gender-balanced, whereas there are more female than male POLY patients. As expected, flare (F) patients from both SJIA and POLY cohorts differ significantly from quiescent (Q) patients for variables reflecting active inflammation: erythrocyte sedimentation rate (ESR), white blood cell count (WBC), platelets (PLT) and joint count (JC, number of affected joints).

Table 1 Systemic scoring system for SJIA patients.
Table 2 Arthritis scoring system for SJIA patients.
Table 3 Arthritis scoring system for POLY JIA patients*
Table 4 Subjects demographic and clinical characteristics

Sample processing

Blood samples were obtained only when there was a clinical need for blood tests. A total of 3 to 4 ml of blood was collected directly in vacutainer cell preparation tubes (CPT) with sodium citrate (Becton Dickinson, Franklin Lakes, NJ, USA). Peripheral blood mononuclear cells (PBMCs) were isolated within three hours of collection by centrifugation of CPT tubes, per the manufacturer's instructions.

RNA preparation

Purified PBMCs were lysed in RLT reagent (Qiagen, Valencia, CA, USA) and lysate was stored at -80°C until RNA extraction. RNA was isolated using the RNeasy mini kit (Qiagen), per the manufacturer's instructions with an additional on-column DNase I (Qiagen) treatment for 40 minutes. The RNA concentration was measured by the Ribogreen assay (Molecular Probes, Grand Island, NY, USA) or by absorbance at 260 nm. The purity of RNA was assessed by the ratio of the absorbance readings at 260 and 280 nm. The integrity of the RNA samples was also checked by either agarose gel electrophoresis or with the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).

Gene panel selection

In a pilot study, paired flare/remission PBMC samples from 14 SJIA patients were processed for RNA as described [36] and analyzed using Lymphochip cDNA microarrays (Patrick Brown, Stanford University, Stanford, CA, USA) [37, 38]. A large number of genes were identified as differentially expressed in flare versus remission samples by Significance Analysis of Microarrays (SAM) [39]. Hierarchical clustering was performed with the Cluster program and visualized using TreeView [40] (Eisen Lab, University of California, Berkeley, CA), as illustrated for a subset of genes in Additional file 1, and also in [36, 41]. The full data set, GSE37388, is released to the public on the Gene Expression Omnibus (GEO) database. From the large set, we selected genes (n = 131) representing various ontologic categories, such as signaling, transcription, inflammation and immune function. We then added other immune-related genes (n = 50) that are expressed in PBMC and implicated in JIA or RA by published reports. The genes were selected prior to analysis of any blood samples for this study, and the samples used for the microarray experiment were not re-used here. The 181 selected genes are shown on Additional file 2; we confirmed that many are immune-related using the program PANTHER 7.0 (Protein ANalysis THrough Evolutionary Relationships) Classification System (Thomas Lab, University of Southern California, Los Angeles, CA, USA), which classifies proteins by their functions, using published experimental evidence and evolutionary relationships [42] (http://www.pantherdb.org/) to categorize their biological functions. This analysis showed that the largest functional category is inflammatory chemokine and cytokine signaling pathways (14.6% of the genes), followed by interleukin signaling pathways (10.8%), apoptosis signaling pathways (9.9%) and toll receptor signaling pathways (6.2%). A full list of categories covered is shown on Additional file 3.

Gene expression detection by kinetic PCR

The kinetic RT-PCR assay was performed as described [43]. Briefly, all reactions were carried out in duplicate as a single-step RT-PCR reaction, using SYBR green chemistry. Data from duplicate reactions for each gene were averaged and normalized based on levels of expression of four housekeeping genes: eukaryotic translation elongation factor 1 alpha1 (EEF1A1), protein phosphatase 1, catalytic subunit, gamma isoform (PPP1CC), ribosomal protein L12 (RPL12), and ribosomal protein L41 (RPL41). The normalized expression level, housekeeping normalized units, of each gene was used to determine the fold change among samples. In a preliminary experiment, we found that a subset (n = 75) of our gene panel showed very limited variation in level (± 2-fold difference from the mean value) in five healthy individuals (two females and three males) over a four-month period (data not shown).

Identification of ESR or JC significantly associated genes in SJIA and POLY

Genes significantly associated with SJIA and POLY were determined using Pearson's correlation and Student's t-test, as explained in the Results section.

Significance analysis of the canonical biological pathways

The biological pathways indicated by the group of genes associated with each clinical parameter/patient cohort subset were determined by pathway analysis with Ingenuity IPA system (Ingenuity Systems, Redwood City, CA, USA; http://www.ingenuity.com). The significance of either ESR or JC related pathways was analyzed using sparse linear discriminant analysis method, as previously described [44]. Correlation between SJIA ESR-related and JC-related pathways was analyzed by Pearson correlation. To determine a threshold to extract pathways that significantly differentiate ESR and JC in SJIA, 500 simulated SJIA ESR-related and 500 simulated JC-related pathway data sets were created by permutation of canonical pathway identifications and their associated pathway P-values for SJIA ESR or JC. For each canonical pathway, the absolute P-value difference in logarithm form between SJIA ESR and JC was computed using one of the 500 simulated SJIA ESR and one of the 500 simulated JC pathway P-value data sets. This led to 500 absolute log P-value differences for each canonical pathway between SJIA ESR and JC, which later were sorted and 20%, 50% and 80% values were computed. Densities of the absolute differences between SJIA ESR and JC-related pathways for the original and the simulated data sets (20%, 50%, and 80%) were plotted using the R package. Comparison of the original data set and the 80th percentile simulated data set determined the threshold to select significantly different pathways between SJIA ESR and JC. A similar approach was applied to the analysis of significantly different pathways between SJIA ESR and POLY ESR.

Results

ESR and JC-associated gene expression in JIA

ESR was chosen as a quantitative measure of systemic inflammation for our analysis, as it typically rises in association with flares of systemic symptoms and was assessed in the largest number of samples. We also considered another measure of systemic inflammation, C-reactive protein (CRP), but few samples were assessed for CRP, precluding the use of this parameter in our analysis. The number of affected counts (joint count, JC), as defined above (Introduction), was used as a quantitative measure of arthritis.

Our samples were initially classified as flare or quiescence based on criteria that we have developed for analysis of JIA (Tables 1, 2 and 3), as previously published [10], and ESR and JC are part of these criteria. We performed a distribution analysis of ESR or JC values by disease states (flare/quiescence) using R Epicalc package (http://cran.r-project.org/web/packages/epicalc/) to investigate if additional subgroups would be revealed. Visual inspection of the results show that the SJIA and the POLY flare patients could be partitioned into two groups related to their ESR values (Figure 1A): F1, with ESR values below 20, and F2, with ESR values above 20. All patients in the F1 subgroup had mild flares by our other criteria (not shown). Quiescence samples all had ESR below 20, clustering together with the F1 flare group. This analysis also showed that, in our samples, JC values in the flare and quiescence disease states are generally non-overlapping in both SJIA and POLY patients, with quiescence samples with 0 or 1 joint count, and all flare samples above zero (Figure 1B).

Figure 1
figure 1

Distribution of ESR or JC values by Flare (F) or Quiescence (Q) disease states. (A) The SJIA/POLY F patients are partitioned into two groups (F1 and F2) by ESR value < 20. (B) SJIA/POLY JC values are generally non-overlapping, in accordance with the F and Q disease states.

We analyzed the association of the 181 gene panel with ESR and JC in both SJIA and POLY samples, using the strategies delineated in Figure 2. Genes whose expression was significantly associated with ESR or JC in SJIA and POLY cohorts were identified in two ways. As described in Figure 2A, Pearson correlation analyses were performed to correlate ESR or JC values with patient expression data sets. To assess the significance of these findings, we calculated the global false discovery rate (gFDR) by 100-fold permutation of normalized kPCR data. After determining the gFDR, local FDR (lFDR) analysis can compute and assign significance measures to all features [45]. A cut-off value of lFDR ≤ 0.05 was used to select significant genes for downstream pathway analysis.

Figure 2
figure 2

Schematic of the experimental design and result summary. (A) Pearson correlation analysis and (B) Student's t-test were used to identify ESR or JC associated genes in SJIA and POLY.

We also analyzed gene expression association using Student's t-test, as shown in Figure 2B. For ESR, based on the analysis from Figure 1A, we initially divided our samples into three groups: flare samples with ESR < 20 (F1), flare samples with ESR > 20 (F2), and quiescence (to ensure that differences between F1 and Quiescence were not overlooked). We identified genes whose mean expression value differed significantly between the F1 (ESR < 20) and F2 (ESR > 20) patient groups, but no differences in genes expressed by the F1 and the quiescence group were found. Subsequently, we grouped the flare F1 and the quiescence groups into one group for ESR analysis. For JC, no other partitioning was necessary, as shown in Figure 1B, and samples were grouped into flare and quiescence groups. As we did previously for Pearson analysis, we calculated local FDR and a value of < 0.05 was considered significant (Figure 2B). This second analysis found genes missed by correlation analysis, as the latter requires a linear relationship and captures genes with more tightly regulated expression (small differences between F and Q samples).

Pearson correlation analysis of expression data from SJIA subjects found 79 genes from our panel to be ESR-correlated and 36 genes to be JC-correlated. Student's t-test found 66 ESR-associated and 79 JC-associated genes in SJIA. This pattern differed from relationships of the expression levels of the same genes with these clinical parameters in POLY-course JIA patients: 20 ESR-correlated and no JC-correlated genes were found in POLY, and none of the genes were ESR-associated or JC-associated by Student's t-test in POLY. Combining both analyses, we found 91 ESR-related and 92 JC-related genes in SJIA, and 20 ESR-related and no JC-related genes in POLY. A list of significantly associated genes is on Additional file 4. Additional file 5 (Supplementary Figure 2A-F) diagrams the fold changes in expression of the selected genes between groups (for example, F2 versus F1 + Q) and between quartiles of ESR or JC. The probability density analysis graphically represents the normalized frequency distribution of the fold ratio of the selected two groups. This result indicates that our selected genes have significant variation between groups (limited variation = fold ratio close to 1) while showing strong association. This analysis further supports our approach (Figure 2) for the identification of significant associations. The reduced number of genes associated with these clinical parameters in the POLY cohort was not surprising given that the gene list was chosen in large part using expression data from SJIA PBMC. Indeed, this finding implies a degree of specificity of the associated genes for SJIA (see discussion). Another likely contribution to this difference might be the extent that disease-related processes are reflected in peripheral blood in the two disease types.

Comparative analysis of SJIA ESR and JC related pathways

Using the lists of associated genes, we determined biological pathways associated with each clinical parameter/patient cohort by pathway analysis with the Ingenuity IPA system. ESR-related and JC-related pathways were then compared, to investigate whether the same biological pathways are involved in ESR and JC elevations in SJIA. As shown in Figure 3A, there is strong correlation (Pearson correlation coefficient, 0.91) between SJIA ESR and JC-related pathways (n = 189), implying that some of the same pathways play roles in the systemic and arthritic components of the disease. Shown in Figure 3B, densities of the absolute log P-value differences of all pathways between SJIA ESR and JC, for the original and the 20, 50 and 80 percentile of the simulated random data sets, were computed and plotted.

Figure 3
figure 3

Analyses of biological pathways associated with SJIA ESR, SJIA JC and POLY ESR. Comparative analyses determine pathways that are significantly different in their association with SJIA ESR vs JC (top panels) or that differentiate SJIA/ESR and POLY/ESR cohorts (bottom panels). (A) Pearson correlation analysis between SJIA ESR and JC related pathways identified by Ingenuity Analysis. (B) Simulation analyses to determine the threshold for significantly differentiating pathways between ESR and JC in SJIA. Densities of the absolute differences between SJIA ESR- and JC-related pathways, for the original and the simulation data sets (20%, 50% and 80%) are shown. See Methods for more details. (C) Pearson correlation between SJIA and POLY ESR-related pathways identified by Ingenuity Analysis. (D) Simulation analysis to determine the threshold for significantly differentiating pathways between SJIA and POLY ESR, performed as in (B).

Significantly differentiating pathways between ESR and JC in SJIA revealed by this analysis are in Table 5 (top two pathways), as were pathways that were correlated comparably with both ESR and JC (Table 5). The only pathway more significantly related to SJIA ESR than to SJIA JC was the glucocorticoid (GC) receptor signaling pathway. The expression of most of the genes in this pathway was higher in samples with higher ESR compared to samples with lower ESR. In contrast, the PI3K/Akt signaling pathway was more significantly related to SJIA JC. Though the significance of the association favored JC, the expression levels of most genes in this cell survival pathway were higher in samples with higher ESR or higher JC. However, as might be expected, TP53, which encodes p53, a pro-apoptotic, negative regulator of the Akt pathway [46], was down-regulated in association with JC and ESR elevations. Consistent with these results, we have previously reported that purified monocytes have lower TP53 transcript levels and increased cellular resistance to apoptotic stimuli during SJIA flare compared to quiescence [36]. Overall, the identification of some pathways that are differentially correlated with ESR and JC raises the possibility of differences in aspects of the immunobiology of arthritis compared to systemic inflammation in SJIA, as discussed below.

Table 5 Biologic pathways that distinguish or are shared between SJIA ESR and SJIA JC

A number of pathways were significantly related to SJIA ESR and JC to the same degree (Table 5). For several of these pathways, the expression of most of the associated genes was higher in samples with higher ESR or JC. These pathways include, among others, protein kinase receptor (PKR, a pattern-recognition receptor) signaling in interferon induction, T cell and B cell signaling in the pattern of rheumatoid arthritis (RA), and (macrophage) migration inhibition factor (MIF) regulation of innate immunity. Other genes in pathways associated with activating innate responses, such as lipopolysaccharide (LPS) signaling and triggering receptor expressed on myeloid cells (TREM1) signaling are also higher samples with either higher ESR or JC. Genes in other pathways showed lower expression in samples with higher ESR or JC, such as T helper cell differentiation, iCOS-iCOSL (inducible T-cell co-stimulator/ligand) signaling in T helper cells and CD40 (co-stimulatory molecule on antigen presenting cells) signaling. Notably, these down-regulated pathways are associated with adaptive immune responses. Also down-regulated in association with elevations of both SJIA ESR and JC is the pathway for crosstalk between dendritic cells and natural killer cells, which can be involved in restriction of innate responses [47].

Two genes, the DNA repair enzyme ATM and the transcription factor NFATC2 (also known as NFAT1), are in the pathway for RANK signaling in osteoclasts and are both down-regulated in association with systemic (ESR) and arthritic (JC) disease activity. The Rank/RankL pathway is an important regulator of bone remodeling [48]. An ATM deficiency has been described in CD4+ T cells from rheumatoid arthritis (RA) patients [49], associated with premature immunosenescence. However, ATM may also be involved in bone formation, and ATM deficient animals show increased numbers of osteoclasts [50]. The transcription factor NFATC2 has been identified as a negative regulator of cartilage cell growth [51]. It is also important in T cell effector function, translocating to the nucleus following T cell receptor activation and regulating expression several cytokines in CD4 T cells (reviewed in [52]). Thus, its inverse correlation with ESR and JC may be similar to the other T cell-related pathways described above. Interestingly, hyperactivation of NFATC2 in T cells is associated with decreased susceptibility to experimental autoimmune encephalomyelitis, indicating that increased NFATC2 activity may have immunomodulatory effects that down-regulate autoaggressive reactions [53].

Comparative analysis of SJIA and POLY ESR-related pathways

We next asked whether some biological pathways involved in SJIA ESR elevation are also involved in POLY ESR elevation, by comparing SJIA and POLY ESR-related genes. As shown in Figure 3C, there is reduced correlation (correlation coefficient, 0.59) between SJIA and POLY ESR-related pathways (n = 119), compared to the correlation we observed between SJIA ESR- and SJIA JC-related pathways. Shown in Figure 3D, densities of the absolute logarithm P-value differences of all pathways between SJIA and POLY ESR, for the original and the 20, 50 and 80 percentile of the simulated random data sets, were computed and plotted.

Several pathways differ significantly between SJIA and POLY ESR, as quantified by absolute difference between the SJIA and POLY pathways (Table 6). These include: the role of macrophages, fibroblasts and endothelial cells in RA, IL-10 signaling, glucocorticoid receptor signaling, among others. These data suggest a greater role for these pathways in SJIA, compared to POLY. However, very few genes were associated with POLY ESR in most pathways, resulting in low significance of association with POLY (not shown). In these differentiating pathways, the (few) genes correlating with ESR were higher in POLY samples with higher ESR values, suggesting that these genes, perhaps in the context of other pathways or in the context of the identified pathways but within the joint, contribute to inflammation in polyarticular course JIA. Indeed, evidence from RA, the adult disease most similar to polyarticular JIA, implicates monocyte and macrophage activation [54] and endothelial cell dysfunction [55], both in joints and in the periphery.

Table 6 Biologic pathways that distinguish or are shared between SJIA ESR and POLY ESR

As observed in the previous analysis, pathways associated with T cell responses are significantly associated with ESR in SJIA but the genes in these pathways show lower expression in samples with higher ESR in comparison to samples with lower ESR. In addition, this analysis showed that genes in B cell activating factor (BAFF) signaling, April (A proliferation-inducing ligand, TNFSF13)-mediated signaling and IL-15 signaling pathways show lower expression samples with in elevated ESR in SJIA (Table 6).

Comparative analysis of joint count (JC) correlated genes in systemic and arthritis phase (SAF) and arthritis phase (AF) SJIA patients

Using the previously identified JC-associated genes (Additional file 4), we then investigated whether arthritis-related gene pathways change when the disease phenotype changes from the earlier systemic and arthritic activity/flare (SAF) phase to arthritis-only activity/flare (AF) phase. Shown in Figure 4A, SJIA patients were distributed according to values of JC and systemic scores (Tables 1 and 2-) to identify SAF and AF subgroups. Figure 4B shows the JC associated genes that are significantly correlated with JC in SAF and AF subgroups. Within the SJIA SAF group, only IL-10 was identified to positively correlate with JC (P-value 0.026). In contrast, 12 genes were found to significantly correlate (negatively) with JC in AF subgroup (Figure 4B, listed in order of decreasing significance): TRAP1, IL2RG, CD40LG, PARP1, TP53, ATM, NFATC2, GZMA, CASP10, PFKFB3, IRF3 and IRF4. Canonical pathway analysis (Figure 4C) mapped 8 of the 12 JC-correlated genes to a single network with the Th2 cytokine, IL-4, at its center. These functional relationships suggest that lack of IL-4 may contribute to arthritis in the AF subgroup. The difference in JC-correlated genes/pathways between AF and SAF supports the hypothesis that different biological pathways are engaged in the chronic arthritis stage versus the more acute (or systemic symptom-associated) arthritis of SJIA.

Figure 4
figure 4

Analysis of JC-associated genes in SJIA patients with high JC counts but different systemic scores. (A) Distribution analysis of SJIA patients, according to both joint count (JC) and systemic scores, established systemic plus arthritis flare (SAF) and arthritis flare (AF) subgroups. (B) Pearson correlation analysis relating JC values to the associated gene expression in both the SAF and AF subgroups. (C) Pathway analysis of the JC correlated genes in AF subgroup. Differentially JC-correlated genes were analyzed using data mining software (Ingenuity Systems) to identify relevant biological pathways associated with each clinical subgroup. The node color indicates the degree of positive (red) or negative (green) correlation to SJIA AF. Nodes are displayed using shapes that represent the functional class of the gene product and different relationships are represented by line type (see legend). Relationships are primarily due to co-expression, but can also include phosphorylation/dephosphorylation, proteolysis, activation/deactivation, transcription, binding, inhibition, biochemical modification.

Discussion

In this study, we sought to identify molecular pathways involved in the systemic and arthritic components of SJIA by investigating the gene expression pathways associated with increases in ESR and active joint count. We chose ESR as a marker of systemic inflammation, but we note that SJIA flares associated with elevated ESR may also include arthritis. Further, SJIA flares with macrophage activation syndrome (MAS) may actually lower ESR from fibrinogen consumption as part of the coagulopathy [56]. The latter issue does not confound our analysis, as the three flare samples with low ESR were from patients with mild flares without MAS. Strictly speaking, our approach delineated gene associations with ESR; however, in our group of SJIA samples, ESR typically correlated closely with other evidence of systemic disease.

Several variables that influence transcriptional profiles should be considered in relation to our results. It is possible that some of the observed differences in gene expression are due to differences in cell type composition of PBMC between SJIA and POLY, or between flare and quiescence [10]. Changes in abundance of cell types may be relevant to disease mechanisms. For monocyte-related genes, we [41] and others [9] have shown that the differences in transcript abundance are not explained by differences in monocyte numbers alone, but reflect activation state. The use of medication and disease duration at the time of sampling may influence the pattern of gene expression. A larger, likely multi-center, study will be needed to rigorously control for these important variables.

Our analysis revealed overlap in molecular pathways involved in increased ESR and elevated JC in SJIA. This result was not unexpected, given reported correlations between these two parameters [30, 57]. However, the glucocorticoid receptor (GCR) signaling pathway was more significantly related to ESR than JC. Systemic symptoms of SJIA respond to exogenous steroids, suggesting the elevation of GCR signaling may represent an endogenous effort to dampen systemic inflammation. The comparable doses of exogenous steroids in the F and Q groups make it less likely that steroid therapy is inducing this pathway. Notably, polymorphism in the GCR gene is associated with the level of inflammatory activity in JIA [58]. Involvement of GCR signaling in systemic inflammation in SJIA and stronger association of this pathway with inflammation in SJIA versus POLY (at least as reflected in blood cells) is consistent with reduced responses in SJIA patients to non-glucocorticoid drugs that are efficacious in subsets of POLY patients (for example, methotrexate and anti-TNFα [59, 60]).

We also found that the PI3K/Akt signaling pathway is more significantly related to SJIA JC than ESR. This pathway, which is activated by a variety of stimuli, including IL-1β, TNFα and IL-6, is potentially involved in IL-17 production [61]. IL-17 could be an important factor in SJIA arthritis [19], particularly in the later phase. We did not assess expression of IL-17 in this study, but our preliminary data suggest that CD4+ T cells from SJIA patients secrete higher levels of IL-17 than control cells when cultured in TH17-polarizing conditions [Wong M, Mellins E, unpublished results]. Recently, enrichment of Th17 (and Th1) cells in blood of SJIA patients has been described [62].

Our findings are consistent with the hypothesis that dysregulation of the innate immune system makes a more prominent contribution to SJIA immunopathology than alterations of the adaptive immune system [17, 63], whereas adaptive responses are thought to drive oligoarticular and polyarticular JIA [64, 65]. However, our results implicate deficiencies in genes associated with T cell-related responses in SJIA pathology, similar to observations in other studies [29]. For example, reduced cytolytic cell activity [66] and diminished function of T regulatory cells may play roles in SJIA etiology [67]. Down-regulated genes associated with cytolytic function also participate in dendritic cell/NK cell and monocyte/NK cell interaction. Some cytolysis genes are part of the IL-15 signaling pathway, and IL-15 is involved in the development of NK cells [68].

In the systemic plus arthritic stage of SJIA, we found that expression of IL-10 in PBMC was positively associated with arthritis. In in vitro studies of SJIA monocytes, we and others [26, 29] observe that IL-10 is expressed after TLR stimulation, and IL-10 signaling is intact [Macaubas et al., unpublished]. Given the immunosuppressive effect of IL-10, association of this gene with arthritis in SJIA may represent an attempt by the immune system to reduce inflammation. The level of IL-10 may be inadequate to deal with the inflammatory challenge, as the frequency of a promoter allele associated with low IL-10 expression is increased in SJIA patients [69, 70]. We found that LPS-induced production of IL-10 protein in SJIA monocytes is comparable to controls [41].

A striking finding of this study is that deficiency in IL-4-related pathways correlates with JC in the arthritic phase of SJIA. IL-4 has been implicated in protection against arthritis. Polymorphism in the IL4Rα gene that confers reduced responsiveness to IL-4 is associated with worse outcome in RA [71]. Low levels of circulating IL-4 are observed in patients with active POLY [72]. IL-4 has been shown to suppress growth factor-induced proliferation of cultured rheumatoid synovial cells by interfering with the cell cycle and by decreasing cell survival [73]. In the murine model of collagen-induced arthritis, IL-4 is protective against cartilage and bone destruction [74], and neutralization of IL-4 in the same model results in reversal of arthritis suppression [75]. IL-4 is also protective in the model of proteoglycan-induced arthritis [76]. Interestingly, in proteoglycan-induced arthritis, mice deficient in IL-4Rα showed higher IL-1β, IL-6 and MIP1a, whereas levels of IFNγ and autoantibodies were less affected. These results imply that IL-4 suppresses innate immune activity more than the adaptive system in this arthritis model [77]. This might model the arthritis of late stage SJIA.

IL-4 inhibits expression of pro-inflammatory cytokines, such as IL-1β, TNFα [78, 79] and IL-17 [80]. As mentioned, IL-17 is an attractive candidate for a driver of inflammation in the later arthritic phase of SJIA. Th17 cells may become IL-1 independent in SJIA, as seen in an animal model [81]. The IL-1β independence of IL-17 action would be consistent with the decreased efficiency of anti-IL1 therapy in the later arthritic phase of SJIA [19]. The ability of IL-4 to suppress reactivation of committed Th17 cells [82] may be another mechanism by with IL-4 deficiency could contribute to arthritis in SJIA. Finally, in a small, open label study, oral histone deacetylase inhibitors in patients with mean SJIA duration of five years showed significant therapeutic benefit, specifically for arthritis [83]. This finding is consistent with the idea that distinct biology may be involved in later phase arthritis in SJIA.

We found no gene association or correlation linked with POLY joint count, and a limited number of somewhat different genes were associated with elevated ESR in POLY-JIA subjects. Our gene panel was largely derived from a SJIA-based microarray study and, as such, it has a significant bias towards SJIA-related genes. Further, the systemic nature of SJIA predicts more changes in peripheral blood than for POLY, where pathology is more localized. Our POLY cohort was itself heterogeneous, including RF+ and RF- patients, which were analyzed as one group. Most gene expression studies have analyzed RF- patients only [29, 84, 85]; some have not determined the RF status [86]. Griffin et al. 2009 showed that RF+ and RF- patients can share a similar gene signature [21].

It will be of interest to determine the cell type within PBMC that is responsible for particular transcripts. Based on correlated expression patterns with more lineage specific genes, it is most likely that IL-4 transcripts derive from CD4 T cells; the IL-4 message expression is correlated with expression of CD40LG and IL2RG (not shown). In contrast, the IL-10 expression correlates with expression of IL-1, IL-1-related genes and IL-6 (not shown), suggesting IL-10 transcripts are expressed in monocytes. Further studies are also needed to determine the specificity of the SJIA gene signature in relation to other acute inflammatory diseases, such as bacterial and viral infections and other rheumatologic pediatric diseases [27]. Nonetheless, our current results add to the growing evidence that different molecular mechanisms distinguish SJIA from other JIA subtypes [10, 11, 26, 29, 32].

Conclusions

This study demonstrates that analysis of individual clinical parameters in a complex disease like SJIA may reveal unique and informative molecular associations. In addition to elucidating disease immunopathology, this approach may help identify therapeutic targets and strategies tailored to the different phases of SJIA.

Abbreviations

AF:

arthritis-predominant phase

AOSD:

Adult Onset Still Disease

BAFF:

B cell activating factor

CPT:

cell preparation tubes

CRP:

C-reactive protein

EEF1A1:

eukaryotic translation elongation factor 1 alpha1

ESR:

erythrocyte sedimentation rate

F:

SJIA flare

FDR:

false discovery rate

GEO:

Gene Expression Omnibus

GC:

glucocorticoid

GCR:

glucocorticoid receptor

IL:

interleukin

JC:

joint count

LPS:

lipopolysaccharide

MAS:

macrophage activation syndrome

PANTHER:

Protein ANalysis THrough Evolutionary Relationships) Classification System

PBMC:

peripheral blood mononuclear cells

PCR:

polymerase chain reaction

PKR:

protein kinase receptor

PLT:

platelets

POLY:

polyarticular juvenile idiopathic arthritis

PPP1CC:

protein phosphatase 1, catalytic subunit, gamma isoform

Q:

SJIA quiescence

qPCR:

quantitative PCR

RA:

rheumatoid arthritis

RF:

rheumatoid factor

RPL12:

ribosomal protein L12

RPL41:

ribosomal protein L41

SAF:

systemic plus arthritic phase

SAM:

Significance Analysis of Microarrays

SJIA:

systemic onset juvenile idiopathic arthritis

WBC:

white blood count.

References

  1. Petty RE, Southwood TR, Manners P, Baum J, Glass DN, Goldenberg J, He X, Maldonado-Cocco J, Orozco-Alcala J, Prieur A-M, Suarez-Almazor ME, Woo P: International League of Associations for Rheumatology classification of juvenile idiopathic arthritis: second revision, Edmonton, 2001. J Rheumatol. 2004, 31: 390-392.

    PubMed  Google Scholar 

  2. Woo P: Systemic juvenile idiopathic arthritis: diagnosis, management, and outcome. Nat Clin Pract Rheumatol. 2006, 2: 28-34. 10.1038/ncprheum0084.

    Article  PubMed  Google Scholar 

  3. Kumar S, Vaidyanathan B, Gayathri S, Rajam L: Systemic onset juvenile idiopathic arthritis with macrophage activation syndrome misdiagnosed as Kawasaki disease: case report and literature review. Rheumatol Int. 2010.

    Google Scholar 

  4. Wittkowski H, Frosch M, Wulffraat N, Goldbach Mansky R, Kallinich T, Kuemmerle Deschner J, Frühwald M, Dassmann S, Pham T-H, Roth J, Foell D: S100A12 is a novel molecular marker differentiating systemic-onset juvenile idiopathic arthritis from other causes of fever of unknown origin. Arthritis Rheum. 2008, 58: 3924-3931. 10.1002/art.24137.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Lomater C, Gerloni V, Gattinara M, Mazzotti J, Cimaz R, Fantini F: Systemic onset juvenile idiopathic arthritis: a retrospective study of 80 consecutive patients followed for 10 years. J Rheumatol. 2000, 27: 491-496.

    CAS  PubMed  Google Scholar 

  6. Singh-Grewal D, Schneider R, Bayer N, Feldman BM: Predictors of disease course and remission in systemic juvenile idiopathic arthritis: significance of early clinical and laboratory features. Arthritis Rheum. 2006, 54: 1595-1601. 10.1002/art.21774.

    Article  CAS  PubMed  Google Scholar 

  7. Efthimiou P, Kontzias A, Ward C, Ogden N: Adult-onset Still's disease: can recent advances in our understanding of its pathogenesis lead to targeted therapy?. Nat Clin Pract Rheumatol. 2007, 3: 328-335. 10.1038/ncprheum0510.

    Article  CAS  PubMed  Google Scholar 

  8. Syed R, Gilliam B, Moore T: Rheumatoid factors and anticyclic citrullinated peptide antibodies in pediatric rheumatology. Curr Rheumatol Rep. 2008, 10: 156-163. 10.1007/s11926-008-0027-4.

    Article  CAS  PubMed  Google Scholar 

  9. Fall N, Barnes M, Thornton S, Luyrink L, Olson J, Ilowite NT, Gottlieb BS, Griffin T, Sherry DD, Thompson S, Glass DN, Colbert RA, Grom AA: Gene expression profiling of peripheral blood from patients with untreated new-onset systemic juvenile idiopathic arthritis reveals molecular heterogeneity that may predict macrophage activation syndrome. Arthritis Rheum. 2007, 56: 3793-3804. 10.1002/art.22981.

    Article  CAS  PubMed  Google Scholar 

  10. Macaubas C, Nguyen K, Deshpande C, Phillips C, Peck A, Lee T, Park J, Sandborg C, Mellins E: Distribution of circulating cells in systemic juvenile idiopathic arthritis across disease activity states. Clin Immunol. 2010, 134: 206-216. 10.1016/j.clim.2009.09.010.

    Article  CAS  PubMed  Google Scholar 

  11. Maeno N, Takei S, Nomura Y, Imanaka H, Hokonohara M, Miyata K: Highly elevated serum levels of interleukin-18 in systemic juvenile idiopathic arthritis but not in other juvenile idiopathic arthritis subtypes or in Kawasaki disease: comment on the article by Kawashima et al. Arthritis Rheum. 2002, 46: 2539-2541. 10.1002/art.10389.

    Article  PubMed  Google Scholar 

  12. Lotito APN, Campa A, Silva CAA, Kiss MHB, Mello SBV: Interleukin 18 as a marker of disease activity and severity in patients with juvenile idiopathic arthritis. J Rheumatol. 2007, 34: 823-830.

    CAS  PubMed  Google Scholar 

  13. Grom AA, Villanueva J, Lee S, Goldmuntz EA, Passo MH, Filipovich A: Natural killer cell dysfunction in patients with systemic-onset juvenile rheumatoid arthritis and macrophage activation syndrome. J Pediatr. 2003, 142: 292-296. 10.1067/mpd.2003.110.

    Article  CAS  PubMed  Google Scholar 

  14. Lequerre T, Quartier P, Rosellini D, Alaoui F, De Bandt M, Mejjad O, Kone-Paut I, Michel M, Dernis E, Khellaf M, Limal N, Job-Deslandre C, Fautrel B, Le Loet X, Sibilia J: Interleukin-1 receptor antagonist (anakinra) treatment in patients with systemic-onset juvenile idiopathic arthritis or adult onset Still disease: preliminary experience in France. Ann Rheum Dis. 2008, 67: 302-308.

    Article  CAS  PubMed  Google Scholar 

  15. Iwamoto M, Nara H, Hirata D, Minota S, Nishimoto N, Yoshizaki K: Humanized monoclonal anti-interleukin-6 receptor antibody for treatment of intractable adult-onset Still's disease. Arthritis Rheumat. 2002, 46: 3388-3389. 10.1002/art.10620.

    Article  CAS  PubMed  Google Scholar 

  16. Masters S, Simon A, Aksentijevich I, Kastner D: Horror autoinflammaticus: the molecular pathophysiology of autoinflammatory disease (*). Ann Rev Immunol. 2009, 27: 621-668. 10.1146/annurev.immunol.25.022106.141627.

    Article  CAS  Google Scholar 

  17. Mellins E, Macaubas C, Grom A: Pathogenesis of systemic juvenile idiopathic arthritis: some answers, more questions. Nat Rev Rheumatol. 2011, 7: 416-426. 10.1038/nrrheum.2011.68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Gattorno M, Piccini A, Lasiglie D, Tassi S, Brisca G, Carta S, Delfino L, Ferlito F, Pelagatti MA, Caroli F, Buoncompagni A, Viola S, Loy A, Sironi M, Vecchi A, Ravelli A, Martini A, Rubartelli A: The pattern of response to anti-interleukin-1 treatment distinguishes two subsets of patients with systemic-onset juvenile idiopathic arthritis. Arthritis Rheum. 2008, 58: 1505-1515. 10.1002/art.23437.

    Article  CAS  PubMed  Google Scholar 

  19. Nigrovic PA, Mannion M, Prince FH, Zeft A, Rabinovich CE, van Rossum MA, Cortis E, Pardeo M, Miettunen PM, Janow G, Birmingham J, Eggebeen A, Janssen E, Shulman AI, Son MB, Hong S, Jones K, Ilowite NT, Cron RQ, Higgins GC: Anakinra as first-line disease-modifying therapy in systemic juvenile idiopathic arthritis: report of forty-six patients from an international multicenter series. Arthritis Rheum. 2011, 63: 545-555. 10.1002/art.30128.

    Article  CAS  PubMed  Google Scholar 

  20. Quartier P, Allantaz F, Cimaz R, Pillet P, Messiaen C, Bardin C, Bossuyt X, Boutten A, Bienvenu J, Duquesne A, Richer O, Chaussabel D, Mogenet A, Banchereau J, Treluyer J-M, Landais P, Pascual V: A multicentre, randomised, double-blind, placebo-controlled trial with the interleukin-1 receptor antagonist anakinra in patients with systemic-onset juvenile idiopathic arthritis (ANAJIS trial). Ann Rheum Dis. 2011, 70: 747-754. 10.1136/ard.2010.134254.

    Article  CAS  PubMed  Google Scholar 

  21. Griffin T, Barnes M, Ilowite N, Olson J, Sherry D, Gottlieb B, Aronow B, Pavlidis P, Hinze C, Thornton S, Thompson S, Grom A, Colbert R, Glass D: Gene expression signatures in polyarticular juvenile idiopathic arthritis demonstrate disease heterogeneity and offer a molecular classification of disease subsets. Arthritis Rheum. 2009, 60: 2113-2123. 10.1002/art.24534.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. van der Pouw Kraan TC, Wijbrandts CA, van Baarsen LG, Voskuyl AE, Rustenburg F, Baggen JM, Ibrahim SM, Fero M, Dijkmans BA, Tak PP, Verweij CL: 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-1014. 10.1136/ard.2006.063412.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Bennett L, Palucka AK, Arce E, Cantrell V, Borvak J, Banchereau J, Pascual V: Interferon and granulopoiesis signatures in systemic lupus erythematosus blood. J Exp Med. 2003, 197: 711-723. 10.1084/jem.20021553.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Popper S, Shimizu C, Shike H, Kanegaye J, Newburger J, Sundel R, Brown P, Burns J, Relman D: Gene-expression patterns reveal underlying biological processes in Kawasaki disease. Genome Biol. 2007, 8: R261-10.1186/gb-2007-8-12-r261.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Pascual V, Allantaz F, Arce E, Punaro M, Banchereau J: Role of interleukin-1 (IL-1) in the pathogenesis of systemic onset juvenile idiopathic arthritis and clinical response to IL-1 blockade. J Exp Med. 2005, 201: 1479-1486. 10.1084/jem.20050473.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Ogilvie EM, Khan A, Hubank M, Kellam P, Woo P: Specific gene expression profiles in systemic juvenile idiopathic arthritis. Arthritis Rheum. 2007, 56: 1954-1965. 10.1002/art.22644.

    Article  CAS  PubMed  Google Scholar 

  27. Allantaz F, Chaussabel D, Stichweh D, Bennett L, Allman W, Mejias A, Ardura M, Chung W, Wise C, Palucka K, Ramilo O, Punaro M, Banchereau J, Pascual V: Blood leukocyte microarrays to diagnose systemic onset juvenile idiopathic arthritis and follow the response to IL-1 blockade. J Exp Med. 2007, 204: 2131-2144. 10.1084/jem.20070070.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Ishikawa S, Mima T, Aoki C, Yoshio-Hoshino N, Adachi Y, Imagawa T, Mori M, Tomiita M, Iwata N, Murata T, Miyoshi M, Takei S, Aihara Y, Yokota S, Matsubara K, Nishimoto N: Abnormal expression of the genes involved in cytokine networks and mitochondrial function in systemic juvenile idiopathic arthritis identified by DNA microarray analysis. Ann Rheum Dis. 2009, 68: 264-272. 10.1136/ard.2007.079533.

    Article  CAS  PubMed  Google Scholar 

  29. Barnes M, Grom A, Thompson S, Griffin T, Pavlidis P, Itert L, Fall N, Sowders D, Hinze C, Aronow B, Luyrink L, Srivastava S, Ilowite N, Gottlieb B, Olson J, Sherry D, Glass D, Colbert R: Subtype-specific peripheral blood gene expression profiles in recent-onset juvenile idiopathic arthritis. Arthritis Rheum. 2009, 60: 2102-2112. 10.1002/art.24601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Palmisani E, Solari N, Magni-Manzoni S, Pistorio A, Lab E, Panigada S, Martini A, Ravelli A: Correlation between juvenile idiopathic arthritis activity and damage measures in early, advanced, and longstanding disease. Arthritis Rheum. 2006, 55: 843-849. 10.1002/art.22357.

    Article  PubMed  Google Scholar 

  31. Ling XB, Lau K, Deshpande C, Park JL, Milojevic D, Macaubas C, Xiao C, Lopez-Avila V, Kanegaye J, Burns JC, Cohen H, Schilling J, Mellins ED: Urine peptidomic and targeted plasma protein analyses in the diagnosis and monitoring of systemic juvenile idiopathic arthritis. Clin Proteomics. 2010, 6: 175-193. 10.1007/s12014-010-9058-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Ling X, Park J, Carroll T, Nguyen K, Lau K, Macaubas C, Chen E, Lee T, Sandborg C, Milojevic D, Kanegaye J, Gao S, Burns J, Schilling J, Mellins E: Plasma profiles in active systemic juvenile idiopathic arthritis: Biomarkers and biological implications. Proteomics. 2010, 10: 4415-4430. 10.1002/pmic.201000298.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Sandborg CI, Nepom BS, Mellins E: Juvenile arthritis. Clinical Immunology. Edited by: Rich RR, Fleisher T, Shearer W, Kotzin B, Schroeder H. 2001, London: Mosby, 62.1-62.13. 2

    Google Scholar 

  34. Sandborg C, Holmes TH, Lee T, Biederman K, Bloch DA, Emery H, McCurdy D, Mellins ED: Candidate early predictors for progression to joint damage in systemic juvenile idiopathic arthritis. J Rheumatol. 2006, 33: 2322-2329.

    PubMed  Google Scholar 

  35. Ravelli A, Martini A: Juvenile idiopathic arthritis. Lancet. 2007, 369: 767-778. 10.1016/S0140-6736(07)60363-8.

    Article  CAS  PubMed  Google Scholar 

  36. Srivastava S, Macaubas C, Deshpande C, Alexander H, Chang S-Y, Sun Y, Park J, Lee T, Begovich A, Mellins E: Monocytes are resistant to apoptosis in systemic juvenile idiopathic arthritis. Clin Immunol. 2010, 136: 257-268. 10.1016/j.clim.2010.04.003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Alizadeh A, Eisen M, Davis R, Ma C, Sabet H, Tran T, Powell J, Yang L, Marti G, Moore D, Hudson JJ, Chan W, Greiner T, Weisenburger D, Armitage J, Lossos I, Levy R, Botstein D, Brown P, Staudt L: The lymphochip: a specialized cDNA microarray for the genomic-scale analysis of gene expression in normal and malignant lymphocytes. Cold Spring Harb Symp Quant Biol. 1999, 64: 71-78. 10.1101/sqb.1999.64.71.

    Article  CAS  PubMed  Google Scholar 

  38. Alizadeh AA, Eisen MB, Davis RE, Ma C, Lossos IS, Rosenwald A, Boldrick JC, Sabet H, Tran T, Yu X, Powell JI, Yang L, Marti GE, Moore T, Hudson J, Lu L, Lewis DB, Tibshirani R, Sherlock G, Chan WC, Greiner TC, Weisenburger DD, Armitage JO, Warnke R, Levy R, Wilson W, Grever MR, Byrd JC, Botstein D, Brown PO, et al: Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature. 2000, 403: 503-511. 10.1038/35000501.

    Article  CAS  PubMed  Google Scholar 

  39. Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Macaubas C, Nguyen K, Peck A, Buckingham J, Deshpande C, Wong E, Alexander H, Chang S-Y, Begovich A, Sun Y, Park J, Pan K-H, Lin R, Lih C-J, Augustine E, Phillips C, Hadjinicolaou A, Lee T, Mellins E: Alternative activation in systemic juvenile idiopathic arthritis monocytes. Clin Immunol. 2012, 142: 362-372. 10.1016/j.clim.2011.12.008.

    Article  CAS  PubMed  Google Scholar 

  42. Mi H, Lazareva-Ulitsky B, Loo R, Kejariwal A, Vandergriff J, Rabkin S, Guo N, Muruganujan A, Doremieux O, Campbell MJ, Kitano H, Thomas PD: The PANTHER database of protein families, subfamilies, functions and pathways. Nucleic Acids Res. 2005, 33: D284-D288. 10.1093/nar/gki418.

    Article  CAS  PubMed  Google Scholar 

  43. Baechler EC, Batliwalla FM, Karypis G, Gaffney PM, Moser K, Ortmann WA, Espe KJ, Balasubramanian S, Hughes KM, Chan JP, Begovich A, Chang SY, Gregersen PK, Behrens TW: Expression levels for many genes in human peripheral blood cells are highly sensitive to ex vivo incubation. Genes Immun. 2004, 5: 347-353. 10.1038/sj.gene.6364098.

    Article  CAS  PubMed  Google Scholar 

  44. Wu MC, Zhang L, Wang Z, Christiani DC, Lin X: Sparse linear discriminant analysis for simultaneous testing for the significance of a gene set/pathway and gene selection. Bioinformatics. 2009, 25: 1145-1151. 10.1093/bioinformatics/btp019.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Ling X, Cohen H, Jin J, Lau I, Schilling J: FDR made easy in differential feature discovery and correlation analyses. Bioinformatics. 2009, 25: 1461-1462. 10.1093/bioinformatics/btp176.

    Article  CAS  PubMed  Google Scholar 

  46. Gottlieb T, Leal JFM, Seger R, Taya Y, Oren M: Cross-talk between Akt, p53 and Mdm2: possible implications for the regulation of apoptosis. Oncogene. 2002, 21: 1299-1303. 10.1038/sj.onc.1205181.

    Article  CAS  PubMed  Google Scholar 

  47. Walzer T, Dalod M, Robbins SH, Zitvogel L, Vivier E: Natural-killer cells and dendritic cells: "l'union fait la force". Blood. 2005, 106: 2252-2258. 10.1182/blood-2005-03-1154.

    Article  CAS  PubMed  Google Scholar 

  48. Anandarajah A: Role of RANKL in bone diseases. Trends Endocrinol Metab. 2009, 20: 88-94. 10.1016/j.tem.2008.10.007.

    Article  CAS  PubMed  Google Scholar 

  49. Shao L, Fujii H, Colmegna I, Oishi H, Goronzy JJ, Weyand CM: Deficiency of the DNA repair enzyme ATM in rheumatoid arthritis. J Exp Med. 2009, 206: 1435-1449. 10.1084/jem.20082251.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Rasheed N, Wang X, Niu Q-T, Yeh J, Li B: Atm-deficient mice: an osteoporosis model with defective osteoblast differentiation and increased osteoclastogenesis. Hum Mol Genet. 2006, 15: 1938-1948. 10.1093/hmg/ddl116.

    Article  CAS  PubMed  Google Scholar 

  51. Ranger AM, Gerstenfeld LC, Wang J, Kon T, Bae H, Gravallese EM, Glimcher MJ, Glimcher LH: The nuclear factor of activated T cells (Nfat) transcription factor Nfatp (Nfatc2) Is a repressor of chondrogenesis. J Exp Med. 2000, 191: 9-22. 10.1084/jem.191.1.9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Hogan PG, Chen L, Nardone J, Rao A: Transcriptional regulation by calcium, calcineurin, and NFAT. Genes Dev. 2003, 17: 2205-2232. 10.1101/gad.1102703.

    Article  CAS  PubMed  Google Scholar 

  53. Ghosh S, Koralov SB, Stevanovic I, Sundrud MS, Sasaki Y, Rajewsky K, Rao A, Müller MR: Hyperactivation of nuclear factor of activated T cells 1 (NFAT1) in T cells attenuates severity of murine autoimmune encephalomyelitis. Proc Natl Acad Sci USA. 2010, 107: 15169-15174. 10.1073/pnas.1009193107.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Gierut A, Perlman H, Pope RM: Innate immunity and rheumatoid arthritis. Rheum Dis Clin North Am. 2010, 36: 271-296. 10.1016/j.rdc.2010.03.004.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Khan F, Galarraga B, Belch JJF: The role of endothelial function and its assessment in rheumatoid arthritis. Nat Rev Rheumatol. 2010, 6: 253-261. 10.1038/nrrheum.2010.44.

    Article  PubMed  Google Scholar 

  56. Stephan JL, Zeller J, Hubert P, Herbelin C, Dayer JM, Prieur AM: Macrophage activation syndrome and rheumatic disease in childhood: a report of four new cases. Clin Exp Rheumatol. 1993, 11: 451-456.

    CAS  PubMed  Google Scholar 

  57. Crowson CS, Rahman MU, Matteson EL: Which measure of inflammation to use? A comparison of erythrocyte sedimentation rate and C-reactive protein measurements from randomized clinical trials of golimumab in rheumatoid arthritis. J Rheumatol. 2009, 36: 1606-1610. 10.3899/jrheum.081188.

    Article  CAS  PubMed  Google Scholar 

  58. Kostik MM, Klyushina AA, Moskalenko MV, Scheplyagina LA, Larionova VI: Glucocorticoid receptor gene polymorphism and juvenile idiopathic arthritis. Pediatr Rheumatol Online J. 2011, 9: 2-10.1186/1546-0096-9-2.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Russo RA, Katsicas MM: Clinical remission in patients with systemic juvenile idiopathic arthritis treated with anti-tumor necrosis factor agents. J Rheumatol. 2009, 36: 1078-1082. 10.3899/jrheum.090952.

    Article  CAS  PubMed  Google Scholar 

  60. Otten MH, Prince FH, Armbrust W, ten Cate R, Hoppenreijs EP, Twilt M, Koopman-Keemink Y, Gorter SL, Dolman KM, Swart JF, van den Berg JM, Wulffraat NM, van Rossum MA, van Suijlekom-Smit LW: Factors associated with treatment response to etanercept in juvenile idiopathic arthritis. JAMA. 2011, 306: 2340-2347. 10.1001/jama.2011.1671.

    Article  CAS  PubMed  Google Scholar 

  61. Kim K-W, Cho M-L, Park M-K, Yoon C-H, Park S-H, Lee S-H, Kim H-Y: Increased interleukin-17 production via a phosphoinositide 3-kinase/Akt and nuclear factor kappaB-dependent pathway in patients with rheumatoid arthritis. Arthritis Res Ther. 2005, 7: R139-R148. 10.1186/ar1470.

    Article  CAS  PubMed  Google Scholar 

  62. Omoyinmi E, Hamaoui R, Pesenacker A, Nistala K, Moncrieffe H, Ursu S, Wedderburn LR, Woo P: Th1 and Th17 cell subpopulations are enriched in the peripheral blood of patients with systemic juvenile idiopathic arthritis. Rheumatology (Oxford). 2012.

    Google Scholar 

  63. Thompson S, Barnes M, Griffin T, Grom A, Glass D: Heterogeneity in juvenile idiopathic arthritis: impact of molecular profiling based on DNA polymorphism and gene expression patterns. Arthritis Rheum. 2010, 62: 2611-2615. 10.1002/art.27561.

    Article  PubMed  Google Scholar 

  64. Macaubas C, Nguyen K, Milojevic D, Park J, Mellins E: Oligoarticular and polyarticular JIA: epidemiology and pathogenesis. Nat Rev Rheumatol. 2009, 5: 616-626. 10.1038/nrrheum.2009.209.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Lin Y-T, Wang C-T, Gershwin ME, Chiang B-L: The pathogenesis of oligoarticular/polyarticular vs systemic juvenile idiopathic arthritis. Autoimmun Rev. 2011, 10: 482-489. 10.1016/j.autrev.2011.02.001.

    Article  CAS  PubMed  Google Scholar 

  66. Wulffraat NM, Rijkers GT, Elst E, Brooimans R, Kuis W: Reduced perforin expression in systemic juvenile idiopathic arthritis is restored by autologous stem-cell transplantation. Rheumatology. 2003, 42: 375-379. 10.1093/rheumatology/keg074.

    Article  CAS  PubMed  Google Scholar 

  67. de Kleer I, Vastert B, Klein M, Teklenburg G, Arkesteijn G, Yung GP, Albani S, Kuis W, Wulffraat N, Prakken B: Autologous stem cell transplantation for autoimmunity induces immunologic self-tolerance by reprogramming autoreactive T cells and restoring the CD4+CD25+ immune regulatory network. Blood. 2006, 107: 1696-1702. 10.1182/blood-2005-07-2800.

    Article  CAS  PubMed  Google Scholar 

  68. Kennedy MK, Glaccum M, Brown SN, Butz EA, Viney JL, Embers M, Matsuki N, Charrier K, Sedger L, Willis CR, Brasel K, Morrissey PJ, Stocking K, Schuh JCL, Joyce S, Peschon JJ: Reversible defects in natural killer and memory Cd8 T cell lineages in interleukin 15-deficient mice. J Exp Med. 2000, 191: 771-780. 10.1084/jem.191.5.771.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Fife MS, Gutierrez A, Ogilvie EM, Stock CJW, Samuel JM, Thomson W, Mack LF, Lewis CM, Woo P: Novel IL10 gene family associations with systemic juvenile idiopathic arthritis. Arthritis Res Ther. 2006, 8: R148-10.1186/ar2041.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Moller J, Paul D, Ganser G, Range U, Gahr M, Kelsch R, Rosen-Wolff A, Hedrich C: IL10 promoter polymorphisms are associated with systemic onset juvenile idiopathic arthritis (SoJIA). Clin Exp Rheumatol. 2010, 28: 912-918.

    CAS  PubMed  Google Scholar 

  71. Prots I, Skapenko A, Wendler J, Mattyasovszky S, Yon C, Spriewald B, Burkhardt H, Rau R, Kalden J, Lipsky P, Schulze Koops H: Association of the IL4R single-nucleotide polymorphism I50V with rapidly erosive rheumatoid arthritis. Arthritis Rheum. 2006, 54: 1491-1500. 10.1002/art.21832.

    Article  CAS  PubMed  Google Scholar 

  72. Spelling P, Bonfa E, Caparbo VF, Pereira RM: Osteoprotegerin/RANKL system imbalance in active polyarticular-onset juvenile idiopathic arthritis: a bone damage biomarker?. Scand J Rheumatol. 2008, 37: 439-444. 10.1080/03009740802116224.

    Article  CAS  PubMed  Google Scholar 

  73. Dechanet J, Briolay J, Rissoan M, Chomarat P, Galizzi J, Banchereau J, Miossec P: IL-4 inhibits growth factor-stimulated rheumatoid synoviocyte proliferation by blocking the early phases of the cell cycle. J Immunol. 1993, 151: 4908-4917.

    CAS  PubMed  Google Scholar 

  74. Joosten LA, Lubberts E, Helsen MM, Saxne T, Coenen-de Roo CJ, Heinegrd D, van den Berg WB: Protection against cartilage and bone destruction by systemic interleukin-4 treatment in established murine type II collagen-induced arthritis. Arthritis Res. 1999, 1: 81-91. 10.1186/ar14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Myers LK, Tang B, Stuart JM, Kang AH: The role of IL-4 in regulation of murine collagen-induced arthritis. Clin Immunol. 2002, 102: 185-191. 10.1006/clim.2001.5162.

    Article  CAS  PubMed  Google Scholar 

  76. Horsfall A, Butler D, Marinova L, Warden P, Williams R, Maini R, Feldmann M: Suppression of collagen-induced arthritis by continuous administration of IL-4. J Immunol. 1997, 159: 5687-5696.

    CAS  PubMed  Google Scholar 

  77. Cao Y, Brombacher F, Tunyogi Csapo M, Glant T, Finnegan A: Interleukin-4 regulates proteoglycan-induced arthritis by specifically suppressing the innate immune response. Arthritis Rheum. 2007, 56: 861-870. 10.1002/art.22422.

    Article  CAS  PubMed  Google Scholar 

  78. Essner R, Rhoades K, McBride WH, Morton DL, Economou JS: IL-4 down-regulates IL-1 and TNF gene expression in human monocytes. J Immunol. 1989, 142: 3857-3861.

    CAS  PubMed  Google Scholar 

  79. Hart PH, Vitti GF, Burgess DR, Whitty GA, Piccoli DS, Hamilton JA: Potential antiinflammatory effects of interleukin 4: suppression of human monocyte tumor necrosis factor alpha, interleukin 1, and prostaglandin E2. Proc Natl Acad Sci USA. 1989, 86: 3803-3807. 10.1073/pnas.86.10.3803.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Lubberts E, Joosten LA, Chabaud M, van Den Bersselaar L, Oppers B, Coenen-De Roo CJ, Richards CD, Miossec P, van Den Berg WB: IL-4 gene therapy for collagen arthritis suppresses synovial IL-17 and osteoprotegerin ligand and prevents bone erosion. J Clin Invest. 2000, 105: 1697-1710. 10.1172/JCI7739.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Koenders M, Lubberts E, Oppers Walgreen B, van den Bersselaar L, Helsen M, Kolls J, Joosten LAB, van den Berg WB: Induction of cartilage damage by overexpression of T cell interleukin-17A in experimental arthritis in mice deficient in interleukin-1. Arthritis Rheum. 2005, 52: 975-983. 10.1002/art.20885.

    Article  CAS  PubMed  Google Scholar 

  82. Cooney LA, Towery K, Endres J, Fox DA: Sensitivity and Resistance to Regulation by IL-4 during Th17 Maturation. J Immunol. 2011, 187: 4440-4450. 10.4049/jimmunol.1002860.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Vojinovic J, Damjanov N, D'Urzo C, Furlan A, Susic G, Pasic S, Iagaru N, Stefan M, Dinarello C: Safety and efficacy of an oral histone deacetylase inhibitor in systemic-onset juvenile idiopathic arthritis. Arthritis Rheum. 2011, 63: 1452-1458. 10.1002/art.30238.

    Article  CAS  PubMed  Google Scholar 

  84. Barnes M, Grom A, Thompson S, Griffin T, Luyrink L, Colbert R, Glass D: Biologic similarities based on age at onset in oligoarticular and polyarticular subtypes of juvenile idiopathic arthritis. Arthritis Rheum. 2010, 62: 3249-3258. 10.1002/art.27657.

    Article  PubMed  PubMed Central  Google Scholar 

  85. Jarvis J, Jiang K, Frank M, Knowlton N, Aggarwal A, Wallace C, McKee R, Chaser B, Tung C, Smith L, McGhee J, Chen Y, Osban J, O'Neil K, Centola M: Gene expression profiling in neutrophils from children with polyarticular juvenile idiopathic arthritis. Arthritis Rheum. 2009, 60: 1488-1495. 10.1002/art.24450.

    Article  PubMed  PubMed Central  Google Scholar 

  86. Barnes MG, Aronow BJ, Luyrink LK, Moroldo MB, Pavlidis P, Passo MH, Grom AA, Hirsch R, Giannini EH, Colbert RA, Glass DN, Thompson SD: Gene expression in juvenile arthritis and spondyloarthropathy: pro-angiogenic ELR+ chemokine genes relate to course of arthritis. Rheumatology. 2004, 43: 973-979. 10.1093/rheumatology/keh224.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

The authors thank colleague scientists at the Stanford University Pediatric Proteomics group for critical discussions, and the Stanford University IT group for excellence in Linux cluster support. XBL was supported by the Stanford NIH/NCRR CTSA award number UL1 RR025744 and by the Lucile Packard Foundation for Children's Health. EDM, CM, YS and CD were supported by funding from the Lucile Packard Foundation for Children's Health (LPFCH) and the Stanford NIH CTSA (UL1 RR025744), The Wasie Foundation, Dana Foundation and the National Institutes of Health (NIH); these funds also supported the manuscript preparation. K-HP, RL and C-JL were supported by the Defense Advanced Research Projects Agency. SNC is the Kwoh-Ting Li Professor in the School of Medicine. Funding bodies played no role in the collection, analysis and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Elizabeth D Mellins.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

XBL performed the statistical analysis and wrote the paper. CM performed statistical analysis, interpreted the data and wrote the paper. HCA, S-YPC and ABB designed and performed the kinetic PCR. QW and EC performed statistical analysis. YS and CD processed patients' samples, prepared RNA, performed and analyzed the microarrays. K-HP, RL, C-JL and SHP performed microarray and initial data analysis. TL and CS provided patients' samples and clinical information. SNC helped design the study and the initial strategy for data analysis. EM contributed to study design, interpreted data and wrote the manuscript. All authors read and approved the final manuscript.

Xuefeng B Ling, Claudia Macaubas contributed equally to this work.

Electronic supplementary material

12916_2012_623_MOESM1_ESM.DOC

Additional file 1: Supplementary Figure 1. Unsupervised hierarchical clustering analysis of a subset of differentially expressed genes from SJIA flare and quiescence samples studied by microarray. Paired samples from 14 subjects at flare (red) and quiescence (green) were studied. Independent samples from the same individual are indicated by *. Each column represents a separate sample; each row represents a separate gene. (DOC 670 KB)

Additional file 2: Supplementary Table 1. The 181 immune-related genes analyzed in this study. (DOC 61 KB)

12916_2012_623_MOESM3_ESM.DOC

Additional file 3: Supplementary Table 2. Annotation of the 181 analyzed genes into different functional categories using PANTHER software. (DOC 32 KB)

Additional file 4: Supplementary Table 3. SJIA ESR, SJIA JC and POLY ESR related gene lists. (DOC 36 KB)

12916_2012_623_MOESM5_ESM.DOC

Additional file 5: Supplementary Figure 2. The fold changes in expression of the genes that differ between groups and between quartiles of ESR or JC. (A) The probability density of fold changes with the selected ESR-related genes between F1 + Q and F2 of SJIA; (B) The probability density of fold changes with the selected ESR-related genes between quartiles (the first quartile to the third quartile) of SJIA; (C) The probability density of fold changes with the selected JC-related genes between Q and F of SJIA; (D) The probability density of fold changes with the selected JC- related genes between quartiles (the first quartile to the third quartile) of SJIA; (E) The probability density of fold changes with the selected ESR-related genes between F2 and F1 + Q of POLY; (F) The probability density of fold changes with the selected ESR-related genes between quartiles (the first quartile to the third quartile) of POLY. (DOC 202 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Ling, X.B., Macaubas, C., Alexander, H.C. et al. Correlation analyses of clinical and molecular findings identify candidate biological pathways in systemic juvenile idiopathic arthritis. BMC Med 10, 125 (2012). https://doi.org/10.1186/1741-7015-10-125

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1741-7015-10-125

Keywords