- Research article
- Open access
- Published:
Correlation analyses of clinical and molecular findings identify candidate biological pathways in systemic juvenile idiopathic arthritis
BMC Medicine volume 10, Article number: 125 (2012)
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.
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).
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).
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.
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.
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.
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.
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.
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
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.
Woo P: Systemic juvenile idiopathic arthritis: diagnosis, management, and outcome. Nat Clin Pract Rheumatol. 2006, 2: 28-34. 10.1038/ncprheum0084.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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
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.
Ravelli A, Martini A: Juvenile idiopathic arthritis. Lancet. 2007, 369: 767-778. 10.1016/S0140-6736(07)60363-8.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Anandarajah A: Role of RANKL in bone diseases. Trends Endocrinol Metab. 2009, 20: 88-94. 10.1016/j.tem.2008.10.007.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Pre-publication history
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1741-7015/10/125/prepub
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
Corresponding author
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)
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)
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
Below are the links to the 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.
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
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1741-7015-10-125