- Research article
- Open Access
Biomarker profiles in serum and saliva of experimental Sjögren's syndrome: associations with specific autoimmune manifestations
Arthritis Research & Therapyvolume 10, Article number: R22 (2008)
Sjögren's syndrome (SS) is a systemic autoimmune disease that mainly targets the exocrine glands. The aim of this study was to investigate the involvement of 87 proteins measured in serum and 75 proteins analyzed in saliva in spontaneous experimental SS. In addition, we intended to compute a model of the immunological situation representing the overt disease stage of SS.
Nondiabetic, nonobese diabetic (NOD) mice aged 21 weeks were evaluated for salivary gland function, salivary gland inflammation and extraglandular disease manifestations. The analytes, comprising chemokines, cytokines, growth factors, autoantibodies and other biomarkers, were quantified using multi-analyte profile technology and fluorescence-activated cell sorting. Age-matched and sex-matched Balb/c mice served as a reference.
We found NOD mice to exhibit impaired salivary flow, glandular inflammation and increased secretory SSB (anti-La) levels. Thirty-eight biomarkers in serum and 34 in saliva obtained from NOD mice were significantly different from those in Balb/c mice. Eighteen biomarkers in serum and three chemokines measured in saliva could predict strain membership with 80% to 100% accuracy. Factor analyses identified principal components mostly correlating with one clinical aspect of SS and having distinct associations with components extracted from other families of proteins.
Autoimmune manifestations of SS are greatly independent and associated with various immunological processes. However, CD40, CD40 ligand, IL-18, granulocyte chemotactic protein-2 and anti-muscarinic M3 receptor IgG3 may connect the different aspects of SS. Processes related to the adaptive immune system appear to promote SS with a strong involvement of T-helper-2 related proteins in hyposalivation. This approach further established saliva as an attractive biofluid for biomarker analyses in SS and provides a basis for the comparison and selection of potential drug targets and diagnostic markers.
Over recent decades the immune system has been subject to much investigation. Growing complexity has often been a major byproduct of the discoveries reported, and subsequently models were established to cope with such complexity. Regarding autoimmune diseases in general, and Sjögren's syndrome (SS) in particular, verifying and expanding such models is desirable, because it has proved difficult to extrapolate findings to existing models that were often developed in different contexts [1–3]. Recent technological advances have greatly increased the amount of information and the number of proteins that can be investigated in any given system and put into a scientific context simultaneously. These technologies, termed transcriptomics, proteomics, metabolomics and other '-omics', were followed by the an increase in systems-based thinking across different scientific disciplines . This trend has promoted systems biology from a technology-driven enterprise to an innovative tool in drug discovery, and it may lead to a more complete perspective on how specific components contribute at different system levels to the immune response. Contextualization and integration have been key drivers of such approaches.
SS (for review [5, 6]), a systemic autoimmune disease, is manifested by severe impairment of exocrine gland function and focal mononuclear cell infiltrates within the salivary and lacrimal glands. The identification of anti-M3 receptor (M3R) autoantibodies for the first time attributed a defined pathogenic role to an autoantibody in SS [7–9]. The roles of other autoantibodies in the pathogenesis of SS (especially SSA [anti-Ro] and SSB [anti-La], which are frequently present in patients with SS) remain to be determined. The disease can involve organs other than the exocrine glands, and the worst disease outcome – lymphoid malignancy – develops in up to 5% of patients with SS. Currently, applied treatments provide merely marginal symptomatic relief [10, 11].
The nonobese diabetic (NOD) mouse, which spontaneously develops both SS-like histopathology and hyposalivation, is the most widely accepted model for SS [12, 13]. Based on the findings of studies conducted in these mice , it is thought that the various SS-related manifestations develop according to a specific time course. However, similar to human SS, the immunological relationship between the two hallmarks of SS, namely salivary gland inflammation and hyposalivation, is far from being understood in NOD mice. Although some diabetes-related genetic loci might contribute to the SS-like disease in NOD mice, both autoimmune diseases can develop independently from each other . Onset of SS in NOD mice is not critically dependent on the diabetes-related H2g7 haplotype. NOD.B10.H2b congenic mice also exhibit an SS-like disease in the absence of overt diabetes , which prevents exclusion of diabetic animals from studies conducted in SS. However, similar to nondiabetic parental NOD mice, they exhibit lymphoid infiltration in the pancreas and, rarely, insulitis . Another model of SS, the C57BL/6.NOD-Aec1Aec2 strain, has not been screened for SS-unrelated autoimmune manifestations other than insulitis . However, the background strain C57BL/6 develops spontaneous organ-specific autoimmune lesions in salivary glands, pancreas, kidneys, lungs and liver, and produces a variety of autoantibodies .
The analytes investigated in this study, which were previously studied in humans or NOD mice within the context of SS, are listed in Additional file 1 (Supplementary table 1). Few studies have assessed interactions between several immune molecules and their association with disease parameters.
Summarizing findings regarding immune mediators such as cytokines in SS, it was concluded that T-helper (Th)2 cytokines are predominant in an early phase of SS, whereas Th1 cytokines are associated with a later stage of the disease . In opposition stands the proposed principle that decreased salivary flow, potentially associated with Th2 cytokines, follows the emergence of glandular inflammation, which was linked to a Th1 response [14, 20]. In addition, the transition between the preclinical and the overt disease state has been associated with shifts in cytokine profiles in NOD mice  and the IL-4/signal transducer and activator of transcription (STAT)6 pathway . Chemokines, small secreted proteins, have been implicated in leucocyte chemoattraction, angiogenesis, fibrosis and malignancy . Despite their uncontested potential as targets for therapeutic intervention, few studies have examined chemokines within the context of SS (Additional file 1 [Supplementary table 1]).
The purpose of the present study was to expand knowledge regarding 87 analytes in serum and 75 proteins in saliva. Thirty and 62 of these molecules have not yet been investigated in SS patients and SS-like disease in NOD mice, respectively. Thirty-six and 54 of the biomarkers have not yet been analyzed in serum and saliva obtained from patients with SS, and neither have 70 of the analytes in serum and 62 biomarkers in saliva from NOD mice been evaluated in a SS-specific context. Based on direct comparison with Balb/c mice, we intended to identify differentially expressed proteins and investigate their potential to discriminate between the disease model and the control strain. This pool of data should also allow computation of a correlation network, representing associations of biomarkers with relevant clinical features of SS in nondiabetic NOD mice, both systemically and locally.
Materials and methods
Animals and assessment of diabetes
Twenty-two female NOD/LtJ (stock #001976) and 19 female Balb/cJ (stock #000651) mice (The Jackson Laboratory, Bar Harbor, ME, USA) were housed in individually ventilated cages at the animal facility of the Department of Physiology, University of Bergen, Bergen, Norway. The study was approved by the Committee for Research on Animals/Forsøksdyrutvalget (project #12-05/BBB).
To serve as controls for subsequent immunostimulatory intervention studies, all mice were injected subcutaneously at 7 weeks of age with 25 μl incomplete Freund's adjuvant emulsified in phosphate-buffered saline. From 10 weeks onward NOD mice were screened weekly for diabetes. Two repeated measurements of glucosuria (>50 mg/dl; Keto-Diabur-Test strips, Roche, Mannheim, Germany) were considered to represent onset of diabetes. At weeks 20 and 21, all mice were screened for hyperglycaemia (>300 mg/dl; Ascensia-microfill, Bayer Healthcare, Mishawaka, IN, USA). At 21 weeks, 10 out of 22 mice (45.5%) were considered to be diabetic and were excluded from all subsequent analyses. We excluded these animals in order to eliminate from our findings any SS-unrelated impact of hyperglycaemia on the physiological process of saliva secretion and the anticipated distorting effect of hyperglycaemia on biomarker profiles.
Measurement of stimulated salivary flow
Mice were fasted but given water ad libitum and anaesthetized with an intramuscular injection of ketamine and medetomidine. Salivary secretion was induced by intraperitoneal injection of 0.5 μg pilocarpine/g body weight (#P6503; Sigma, St. Louis, MO, USA) and collected during 10 minutes. Pre-weighed tubes were weighed again after collection to determine the amount of saliva (1 μg = 1 μl). Protease inhibitor cocktail (#P8340; Sigma) was added at a concentration of 1:500 and samples were kept at -80°C until analysis.
Blood sampling and organ collection
Blood was collected from the saphenous vein from nonanaesthetized mice and by heart puncture on the day of euthanasia. The blood was allowed to clot and centrifuged for 10 minutes at 800 g to obtain serum. The organs were fixed in 4% formalin before embedding in paraffin, sectioning, and staining with haematoxylin and eosin. Sections obtained from the kidneys were also stained using the periodic acid-Schiff staining technique.
Evaluation of salivary gland inflammation and insulitis in the pancreas
After qualitative evaluation of three independent sections, the section with the highest degree of inflammation was recorded as a whole, creating a multiple image-composite picture. A graph tablet was used to select and morphometrically measure the total glandular area and the individual size of each focus. Subsequently, focus score (FS; number of foci of 50 or more mononuclear cells/mm2 glandular tissue) and ratio index (RI; area of inflammation/area of glandular tissue) were determined.
To determine the insulitis score (IS), at least five haematoxylin and eosin stained tissue sections of the pancreas were analyzed in a blinded manner. On average, 32 islets per mouse were scored, as described by Leiter .
Multi-analyte profiles from serum and saliva
A bead-based multiplex sandwich immunofluorescence assay was used to generate multi-analyte profiles (MAPs) from serum and saliva from the 12 nondiabetic NOD and 12 Balb/c mice, comprising 82 analytes for serum and 75 for saliva (Additional file 1 [Supplementary table 2]). Analyses were conducted at Rules Based Medicine Inc. (Austin, TX, USA) using a fully automated system. For each multiplex, eight-point calibrators and three-level controls were included on each microtitre plate. Antibodies used in the MAP to recognize and quantify the specific autoantibodies were directed against all isotypes.
Quantification of anti-M3R antibodies
Levels of anti-M3R autoantibodies were measured as described previously . In brief, aliquots of 2 × 105 Chinese hamster ovary cells, transfected with pcDNA5/FRT/V5-His MsM3R-Flp-In cells, were incubated for 1.5 hours at 4°C with 10 μl of serum before incubation with one of the following fluorescein isothiocyanate (FITC)-conjugated goat anti-mouse detection antibodies (purchased from Southern Biotech) diluted 1:50: isotype control, goat IgG (#0110-02); IgG (H+L; #1031-02); IgG1 F(ab')2 (#1072-02); IgG2b F(ab')2 (#1092-02); IgG2c F(ab')2 (#1079-02); and IgG3 F(ab')2 (#1102-02). The cells were analyzed using a FACSCalibur flow cytometer using Cell Quest software (BD Biosciences, San Jose, CA, USA) and FlowJo (Tree Star Inc., Ashland, OR, USA). The quantities of anti-M3R autoantibodies were analyzed by gating on the FITC-positive population situated above the threshold, set by the sample stained with the secondary antibody alone. The percentage of positive cells was calculated to represent the quantity of anti-M3R autoantibodies.
Means were compared using independent Student's t-test (two-tailed). Bivariate linear associations, used to generate the correlation matrixes, were computed using two-tailed Pearson correlation (r). Strain membership prediction was assessed by discriminant analyses (DA) and subsequent cross-validated (leave one out) group prediction. The quality of the DA function is expressed by its canonical correlation (R*).
Principal component analyses (PCAs) were computed from MAP obtained from NOD mice with the purpose being to uncover the latent structure within protein families. Protein family membership was defined based on the Protein ANalysis THrough Evolutionary Relationships (PANTHER) classification system . PCA seeks a linear combination of variables so that the maximum variance is extracted from the variables. It then removes this variance and seeks a second linear combination, and so forth. Loadings greater than 0.6 were considered defining parts of the component. For proper model specification, variables being either differentially expressed between the two strains (P < 0.05) and/or significantly correlated (r > 0.6; P < 0.05) with one of the disease parameters were included. As a rotation method, Varimax was chosen. The number of components was determined using the Kaiser criterion (Eigenvalue > 1.0). An explanatory criterion (>80%) was applied, in addition, for growth factors and cytokines in serum. Variables being defining parts of a component were also combined for DA and entered simultaneously. In serum no data were missing and in saliva missing values were excluded pair-wise from all analyses except PCA and DA. All analyses were computed using SPSS 13 (SPSS Inc., Chicago, IL, USA).
Autoimmune disease manifestations
At 21 weeks of age salivary secretion (expressed as μl/minute per g bodyweight) in NOD mice (n = 12; 0.367 ± 0.026) was decreased by 42% compared with Balb/c mice (n = 12; 0.637 ± 0.024; P < 0.001; Additional file 1 [Supplementary table 2]). Salivary secretion rate classified 100% of the mice according to their strain membership, confirming the onset of overt SS in all NOD mice.
FS averaged 1.007 ± 0.087 foci/mm2 and RI 0.045 ± 0.007 mm2/mm2 in NOD, whereas Balb/c mice were free from glandular inflammation. IS in nondiabetic NOD mice averaged 0.454 ± 0.052 (Additional file 1 [Supplementary table 2]), representing mild to intermediate insulitis expressed on a scale from 0 (no inflammation) to 1 (all islets are to a large extent invaded by lymphocytes). Importantly, among all variables only serum glutamic-oxaloacetic transaminase (SGOT; r = -0.615, P = 0.033), serum IgA (r = -0.581, P = 0.048) and anti-M3R IgG1 (r = 0.702, P = 0.011) correlated with IS.
Histopathological evaluation of the kidneys, thyroid gland, thymus, heart, lungs, liver, stomach, small and large intestines, appendix and skin revealed a subset of NOD mice exhibiting mononuclear cell infiltration in the kidneys (n = 5; Additional file 1 [Supplementary figure 1A]), accompanied by hyaline casts in two cases (Additional file 1 [Supplementary figure 1B]). In one case, hyaline casts were found in the absence of lymphoid infiltration. In addition, some kidneys exhibited glomeruli with increased numbers of mesangial cells. In one kidney hyaline material was found in glomerular capillaries (Additional file 1 [Supplementary figure 1C]). Necrosis or crescent formation in the glomeruli, however, was not observed. NOD mice exhibiting signs of kidney pathology (n = 6) had significantly lower β2-microglobulin and lower anti-proteinase 3 antibody levels compared with mice free from such alterations (n = 6). Sections of lungs from 11 NOD mice presented foamy cells in the alveoli (Additioinal file 1 [Supplementary figure 1D, E]), which in three cases were accompanied by focal lymphoid infiltrates in the lungs (Additional file 1 [Supplementary figure 1D, F]). However, convincing histological patterns of interstitial lung disease were absent. Using light microscopic screening, no signs of other extraglandular disease manifestations or other independent autoimmune diseases were found.
In serum the levels of 87 biomarkers, including 14 autoantibodies and four anti-M3R antibody subclasses, were measured (Additional file 1 [Supplementary table 2]). Ten proteins, mostly cytokines, were undetectable in serum of NOD and Balb/c. Interferon-γ was detectable in one NOD mouse and circulating fibroblast growth factor-9 in one Balb/c mouse only. In saliva 75 biomarkers were analyzed, of which IL-3 and leptin were undetectable in Balb/c mice; IL-3, leptin and glutathione S-transferase-α were detectable in only one, four and four NOD mice, respectively (Additional file 1 [Supplementary table 2]). Thirty-eight of the analytes assessed in serum were found at significantly different concentrations in NOD mice compared with Balb/c mice, whereas in saliva 34 analytes were significantly different (Table 1).
DA were subsequently computed to investigate each analyte's individual potential and relative importance in accurately predicting mouse strain. Cross-validated classification revealed 18 biomarkers in serum that predicted strain membership with 80% to 100% accuracy (hit rate) and with specificity and sensitivity up to 100%, whereas such capacity was identified for three chemokines measured in saliva (Table 2). Compared with nonvalidated group prediction, cross-validated prediction is based on all cases except the case being classified and it is thought to give a better estimate of the hit rate in the population. For each analyte shown in Table 2 the specificity (percentage of correct predictions in the NOD group) and sensitivity (percentage of correct predictions in the Balb/c group) were calculated.
An immune response is an orchestrated process that involves several protein families and several molecules with similar molecular function. The PANTHER classification system was used to classify the analytes into families of proteins with shared function based on scientific experimental evidence and evolutionary relationships. A correlation matrix comprising the variables identified to be suitable for modelling purposes exhibited profound differences in protein associations within and among protein families (Figure 1). To gain further understanding of these interrelationships in NOD mice, PCA was computed to uncover the underlying dimensions of the immune response (Additional file 1 [Supplementary table 3]). PCA identifies patterns in data and uses a correlation matrix to find the linear combination of original variables, which accounts for most of the variance. It then represents those objects in terms of these new linear combinations, which are called principal components. The first component will be defined to account for the maximum of variation possible (Additional file 1 [Supplementary figure 3]). The second component will subsequently account for the maximum of the remaining variation, and so forth, until a certain criterion set by the researcher is met. In our dataset such patterns were clearly recognizable (Figure 1). Consequently, a large number of the original variables, entered according to protein family membership, could be combined into 27 components that still accounted for at least 80% of the original variance (Additional file 1 [Supplementary table 3]).
Subsequently, linear interrelationships were analyzed using Pearson correlation. The extraction of principal components reduced the correlation matrix considerably, from 9'604 (Figure 1) to 1'994 coefficients (Figure 2). Components associated with disease manifestations are presented in Figures 2 and 3, whereas correlations between components are shown in Figures 2 and 4. Variables for which no reasonable component structure could be computed were included as original variables, if a significant correlation with an autoimmune manifestation was detected (Figures 2 and 3). Results form multivariate DA, based on the simultaneous entry of the variables combined in the different components, are listed in Additional file 1 (Supplementary table 4). They represent the relative importance of the collinear variables, combined in an individual component, in predicting strain membership.
Correlation analyses indicated no linear association between salivary flow and the parameters of glandular inflammation (Figure 2 and 3). PCA identified three components in serum (Se-C) and two in saliva (Sa-C) that correlated with salivary flow (Figure 3). The defining variables of the components are specified in parentheses. Chemokine Se-C-2 (macrophage-inflammatory protein [MIP]-1α, MIP-1γ and monocyte chemoattractant protein [MCP]-5), cytokine Se-C-4 (IL-10) and growth factor Se-C-2 (macrophage-colony stimulating factor and growth hormone) were all positively correlated with salivary flow (Figure 3 and Additional file 1 [Supplementary table 3]). The same applied for the following serum analytes, for which no PCA-based solution could be computed: C-reactive protein, SGOT, vascular cell adhesion molecule (VCAM)-1 and IgA.
In saliva, C-C chemokine ligand/C chemokine ligand Sa-C-2 (eotaxin and macrophage-derived chemokine [MDC]) and cytokine Sa-C-3 (CD40 and CD40L [IL-18 borderline]) correlated with decreased salivary flow. Not reflected by their corresponding component, salivary IL-5 (r = -0.708, P = 0.010), thrombopoietin (r = -0.766, P = 0.004) and factor III (r = -0.614, P = 0.034) also correlated with salivary flow.
Glandular inflammation was negatively correlated with cytokine Se-C-2 (negative loading for IL-1α [-0.94] and positive loading for CD40 ligand [CD40L; 0.68]). Interestingly, this component was correlated with SSB in serum.
In saliva multiple factors exhibited a linear interrelationship with measures of glandular inflammation. Acute phase reactants Sa-C-1 (C-reactive protein, SGOT and serum amyloid P) and coagulation factor Sa-C-2 (von Willebrand factor [vWF] and fibrinogen) correlated negatively with FS and RI. vWF was negatively loaded on coagulation factor Sa-C-2, consistent with its initial positive correlation with increased IR (r = 0.792, P = 0.002). In addition, secretory IgA was inversely correlated with FS. Components extracted from protein families involved in specific immune reactions revealed C-X-C chemokine ligand (CXCL) Sa-C-2 (granulocyte chemotactic protein [GCP]-2) to be positively correlated with IR (r = 0.664, P = 0.018), and cytokine Sa-C-1, combining eight cytokines, exhibited a negative association with IR. In contrast, cytokine Sa-C-2 (leukaemia inhibitory factor, IL-10 and IL-1β) showed an almost significant positive correlation with FS and IR (both r = 0.581, P = 0.061). Correlation patterns compared with other components further supports its connection with glandular inflammation (Figure 2).
We found all isotypes of M3R autoantibodies that we investigated to be significantly increased in NOD mice compared with Balb/c mice. Circulating SSB was the only autoantibody that significantly correlated with any SS disease manifestation. PCA failed to extract components from circulating autoantibody levels, and therefore individual associations are reported in Figures 2 and 3. Interestingly, SSB levels were associated with cytokine Se-C-2 (r = -0.590, P = 0.043) and VCAM-1 (r = 0.587, P = 0.045) and anti-M3R IgG1 correlated with cytokine Se-C-5 (IL-18; r = 0.649, P = 0.023). Levels of anti-M3R IgG3 correlated negatively with MIP-1α (r = -0.626, P = 0.029) and MCP-5 (r = -0.580, P = 0.048), which in turn were associated with salivary flow (r = 0.735, P = 0.006 and r = 0.742, P = 0.006, respectively). Consistent with this observation, levels of anti-M3R IgG3 correlated with cytokine Sa-C-3 (CD40, CD40L and borderline IL-18; r = 0.724, P = 0.012) on its part associated with decreased salivary flow. In addition, the positive association between M3R IgG3 and CXCL Sa-C-2 (GCP-2; r = 0.789, P = 0.002) in its turn correlating with RI, suggests that anti-M3R IgG3 has a role in connecting different SS-related disease manifestations. In addition, salivary matrix metalloproteinase (MMP)-9 also correlated positively with increased anti-M3R IgG3.
Autoantibodies in saliva were grouped into autoantibody Sa-C-1 (anti-RNP, anti-insulin, anti-mitochondrial, SSB and anti-scleroderma-70 antibodies) and autoantibody Sa-C-2 (anti-β2 glycoprotein). Components positively correlated with autoantibody Sa-C-1 were chemokine Se-C-3 (MIP-3β, lymphotactin and MDC), C-C chemokine ligand/XCL Sa-C-1 (lymphotactin, MCP-3, MIP-1β, MCP-1 and RANTES [regulated upon activation, normal T-cell expressed and secreted]), CXCL Sa-C-1 (melanoma growth stimulatory-activity protein, MIP-2, inducible protein [IP]-10 and GCP-2) and growth factor Sa-C-1 (stem-cell factor, fibroblast growth factor-9 and thrombopoietin). In contrast, chemokine Se-C-2 and cytokine Se-C-4, in turn also associated with salivary flow, exhibited a significant negative correlation with autoantibody Sa-C-1. Growth factor Se-C-3 (epidermal growth factor) and coagulation factor Se-C-1 (factors VII and III), together with endothelin-1, correlated negatively with autoantibody Sa-C-1 as well. Autoantibody Sa-C-2 correlated positively with chemokine Se-C-1 (MCP-1, MCP-3, MIP-2 and IP-10) and circulating insulin.
In summary, proteins were mostly associated exclusively with either salivary flow, parameters of glandular inflammation, or autoantibody levels. Dual associations were only apparent for the following: chemokine Se-C-2, cytokine Se-C-4, cytokine Se-C-2, cytokine Sa-C-3 and CXCL Sa-C-2. From variables not included in PCA, only VCAM-1 in serum correlated with two autoimmune parameters.
Associations among components
Correlations among components associated with autoimmune manifestations are presented in Figures 2 and 4. Between components associated with salivary flow, an antagonistic interplay between cytokine Sa-C-3 (associated with decreased salivary flow) and all three components extracted from serum (positively associated with salivary flow) was apparent. Components associated with FS and/or RI correlated strongly with each other, indicating a process in which cytokine Sa-C-2 and CXCL Sa-C-2 oppose acute-phase reactant Sa-C-1 and secreted IgA. In addition, cytokine Se-C-2 and coagulation factor Sa-C-2 played a dual role by combining the positively associated IL-1α and vWF with the negatively associated CD40L and fibrinogen, respectively. Regarding components associated with autoantibodies, the components related to anti-M3R IgG1 correlated negatively with autoantibody Sa-C-1, whereas all components that were positively associated with autoantibody Sa-C-1 correlated with each other.
We observed a complete absence of original variables or components that correlated with salivary flow and either FS or RI. This absence was as absolute when analyzing inter-component correlations (Figure 2). In contrast, some components were associated with either salivary flow or glandular inflammation, having positive or negative inter-relationships with components related to autoantibody levels, such as cytokine Se-C-4 (IL-10) and cytokine Sa-C-1 (interferon-γ, IL-7, oncostatin M, granulocyte macrophage-colony stimulating factor, IL-17, IL-11, IL-5 and IL-18). CXCL Sa-C-2 (GCP-2) exhibited negative associations with the proteins related to increased anti-M3R IgG1 and positive associations with anti-M3R IgG3 and MMP-9. Interestingly, MMP-9 expression correlated positively with components associated with increased glandular inflammation.
By studying the immune system through the application of reductionist principles, its mediators have been thoroughly analyzed over recent decades, which has yielded tremendous scientific advances. However, studying the properties of its isolated components is limited in terms of elucidating how system properties emerge, because they may strongly rely on and arise from interactions between its numerous constituents. In this study we present a substantial amount of data on immunologically relevant proteins, for which often only scant or no information within the context with SS has been published. Multiplexing of analytes thus significantly diminished the sample volume required and enabled us to investigate the analytes' connectivity through correlation networks. Such a study cannot be as conclusive in defining the role of a single protein as a component-focused experimental study design. However, it represents a novel way to analyze the implications of multiple molecules in a specific condition, by providing insight into the inter-relationships that define a specific system state.
Analyses of protein levels instead of mRNA expression data can allow exclusion of certain factors that cause uncertainty, such as RNA stability and correlations between mRNA levels and corresponding protein levels. A recent study combined global gene expression analyses with quantitative proteomics based on two-dimensional gel electrophoresis and mass spectrometry in saliva obtained from patients with SS . Indeed, the correlation between mRNA levels and proteins was proven to be poor. Using mass spectrometry, 42 proteins were identified that were significantly altered when comparing pooled saliva from SS patients with healthy control individuals . None of the proteins identified by Hu and coworkers  was included in our MAP. Two-dimensional gel electrophoresis is limited in terms of its sensitivity in identifying proteins with concentrations in the pg/ml range, especially if they may be masked by abundant proteins present in the biofluid of interest. However, methods such as stable-isotope protein tagging or subtractive proteomics may improve the number of immune system related proteins identified by mass spectrometry. Nevertheless, the two approaches of antibody-based biomarker identification and mass spectrometry-based global proteome profiling may well complement one another in delineating SS-specific disease signatures. Nevertheless, further technological advances are required in both fields to achieve more complete coverage of crucial immune mediators. Unfortunately, molecules of proven importance in SS, such as B-cell activating factor  and type I interferons  were not part of our MAP. Including such crucial molecules in future studies and the present correlation network would indeed greatly increase completeness and enahce the plausibility of an integrated model of SS.
By accounting for SS-related extraglandular autoimmune manifestations and insulitis, we are confident that analyses in serum reflect to a large extent systemic aspects of SS. Nevertheless, although we could not associate the presence of extraglandular disease manifestations with the situation in the salivary glands, it would be of interest to investigate further the involvement of kidneys and lungs in both experimental and human SS [28–30]. The potential presence of SS-unrelated histopathology or of other subclinical autoimmune diseases in mouse models for SS represents a common problem of all mouse models used in SS research [13, 31]. This divergence of these models from clinical SS may lead one to draw erroneous conclusions from measurements in serum. Indeed, biomarker profiling solely focusing on serum or plasma may require considerable validation efforts to prove their specificity for a specific autoimmune disease such as SS.
In contrast, because it is directly collected from the site of inflammation, saliva may largely reflect the immunological situation in the salivary glands. The presence and local production of some inflammatory mediators and autoantibodies in the salivary glands were previously described (Additional file 1 [Supplementary table 1]); these findings were largely confirmed by our study. Importantly, we did not observe a concentration effect related to lesser fluid secretion. In a noninflammatory situation (represented by Balb/c mice), the average correlation coefficient of analytes included in the model did not show a negative association with salivary flow (r = -0.076). The noninvasive collection method and the lack of extraction procedures predispose biomarkers, measured in saliva, as potential surrogate markers of disease and disease activity [25, 32]. In addition, they may reflect the disease independent from other inflammatory conditions in patients. The identification of a set of biomarkers in human saliva, similar to the three chemokines we identified in mice, that predicts the presence of glandular inflammation with high accuracy would represent a major advance in the field.
Traditionally, loss of secretory capacity, degree of lymphoid infiltration and production of specific autoantibodies have been anticipated to correlate with each other and to indicate disease state and severity [5, 6]. However, the correctness of this assumption was difficult to prove. Our findings strongly argue against such close interrelationships and suggest that there is much independence of the various hallmarks of SS. Only RI and SSB correlated directly with each other, and the separation between proteins in terms of whether they were associated with hyposalivation or with glandular inflammation was absolute.
Circulating MIP-1α and MCP-5, which we found to be associated with higher salivary flow, are Th1-related chemokines that are negatively regulated by STAT6 . In opposition, eotaxin and MDC, correlating with decreased salivary flow, are dependent on STAT6  and are considered to be Th2-related chemokines . In accordance with these findings, STAT6-deficient NOD mice did not develop hyposalivation . Both eotaxin and MDC are produced by Th2-promoting dendritic cell types upon engagement of CD40/CD40L ; in our study, these two proteins were associated with low salivary secretion capacity as well. CD40 and CD40L are expressed on salivary gland epithelial cells and infiltrating lymphocytes in biopsies obtained from SS patients . The observed correlation between CD40/CD40L and anti-M3R IgG3 levels may therefore be related to the primary role played by CD40 and CD40L in B-cell survival, B-cell proliferation, antibody production and antibody isotype switching . We found all of the antibody isotypes binding M3R that we measured to be upregulated, with (perhaps most importantly) positive associations between M3R IgG3 and secretory CD40, CD40L, IL-18, GCP-2 and MMP-9. Previous studies identified IgG1 as the crucial isotype in anti-M3R antibody mediated inhibition of salivary flow in STAT6-/-  and IL4-/-  deficient NOD mice. Nevertheless, IgG2 subclasses and IgG3 are generally considered to be significantly more potent in mediating pathogenic effects . In contrast to other isotypes, the effect of IgG3 is FcR independent and strictly related to complement activation; the latter component of the immune system was recently related to SS pathology . In addition, IgG3 can form complexes through self-association and generate cryoglobulins , a feature observed in SS . With regard to the quality of the antibody response, we found total circulating and secretory IgA to be related to lower disease activity.
Apart from B-cell fate, CD40/CD40L ligation plays a central role in converting a tolerogenic antigen presentation into pathogenic immune activation . In conditions with chronic inflammation, CD40 and other inflammatory mediators can have adverse effects on tissue renewal and repair processes . Indeed, we found specific growth factors to be negatively correlated with CD40 and CD40L and secreted autoantibodies. Furthermore, we found circulating IL-10 to correlate negatively with CD40/CD40L; the anti-inflammatory properties of IL-10 may apply to our model, because increased circulating IL-10 did correlate with higher salivary flow and lower autoantibody concentrations in saliva. IL-10 gene transfer in NOD mice partially suppressed the appearance of SS-like features . However, we found IL-10 in saliva to be associated with glandular inflammation, which corroborates reports indicating that IL-10 transgenic mice develop progressive histopathology and hyposalivation, evocative of SS .
With the exception of vWF, increased concentrations of proteins in saliva related to acute tissue injury were all associated with a lower degree of glandular inflammation. These events may be followed later by inflammatory cell invasion into the glandular tissue and/or mirror an imbalance between processes that restore homeostasis and factors that promote chronic inflammation. Higher levels of IL-1 family members did correlate with both worsening of hyposalivation and increased glandular inflammation. IL-1β is also a major inducer of GCP-2  and was, similar to IL-1β, associated with glandular inflammation. Leukaemia inhibitory factor, inducible through IL-1β  and loaded on the same component together with IL-1β and IL-10, has been shown to have parallels with IL-1, tumour necrosis factor-α and IL-6 within the context of promoting inflammation in rheumatoid arthritis .
Among the biomarkers analyzed in saliva, GCP-2, IP-10 and RANTES exhibited the greatest potential in predicting strain membership. GCP-2 is expressed at neutrophil and macrophage dominated inflammatory sites, IP-10 is related to Th1 immune responses, and RANTES can be involved in Th1 and Th2 immune responses . GCP-2 and IP-10 have opposite roles in angiogenesis ; this issue has not yet been addressed in SS, despite the recognized importance of neovascularization in promoting the influx of inflammatory cells . Proteins such as GCP-2, vascular endothelial growth factor, epidermal growth factor, IL-1, VCAM-1 and MMP-9, which we found to be significantly increased in NOD mice, are involved in angiogenesis . In addition, angiostatic molecules such as IP-10, tissue inhibitor of metalloproteinase-1 and endostatin-1 were also altered . Based on our findings, pathogenic neovascularization deserves to receive more attention in SS research.
Increased concentrations of molecules that govern adaptive immune responses were associated with an aggravation of SS coupled with a strong association of Th2-related proteins with hyposalivation. In addition, we observed great independence among the different disease manifestations of SS, which must be considered when one is evaluating potential drug targets. Nevertheless, CD40, CD40L, IL-18, GCP-2 and anti-M3R IgG3 may, however, represent a pivotal point that intersects the various aspects of SS pathology. As multiplex technology allowed us to draw a comprehensive picture of the overt disease stage, we suggest that treatment outcome may similarly be monitored . Hypothesis-driven verification of our findings, together with validation of our findings in human SS specimens, are important goals for the future. We also believe that discovery-driven studies, such as the one presented here, can widen the horizon and inspire new component-focused basic research. Nevertheless, bridging the gap between these two approaches represents a major challenge for the years to come, but the reward will be a more integrated perspective on immunology and autoimmune disease.
= CD40 ligand
= C-X-C chemokine ligand
= discriminant analyses
= fluorescein isothiocyanate
= focus score
= granulocyte chemotactic protein
= inducible protein
= insulitis score
= multi-analyte profile
= monocyte chemoattractant protein
= macrophage-derived chemokine
= macrophage-inflammatory protein
= matrix metalloproteinase
= M3 receptor
= nonobese diabetic
= Protein ANalysis THrough Evolutionary Relationships
= principal component analysis
= regulated upon activation, normal T-cell expressed and secreted
= ratio index
= serum glutamic-oxaloacetic transaminase
= Sjögren's syndrome
= signal transducer and activator of transcription
= vascular cell adhesion molecule
= von Willebrand factor.
Zinkernagel RM, Hengartner H: On immunity against infections and vaccines: credo 2004. Scand J Immunol. 2004, 60: 9-13. 10.1111/j.0300-9475.2004.01460.x.
Anderson CC: Time, space and contextual models of the immunity tolerance decision: bridging the geographical divide of Zinkernagel and Hengartner's 'Credo 2004'. Scand J Immunol. 2006, 63: 249-256. 10.1111/j.1365-3083.2006.01742.x.
Cohn M: A commentary on the Zinkernagel-Hengartner 'Credo 2004'. Scand J Immunol. 2005, 62: 187-194. 10.1111/j.1365-3083.2005.001561.x.
van der Greef J, Martin S, Juhasz P, Adourian A, Plasterer T, Verheij ER, McBurney RN: The art and practice of systems biology in medicine: mapping patterns of relationships. J Proteome Res. 2007, 6: 1540-1559. 10.1021/pr0606530.
Delaleu N, Jonsson R, Koller MM: Sjögren's syndrome. Eur J Oral Sci. 2005, 113: 101-113. 10.1111/j.1600-0722.2004.00183.x.
Jonsson R, Bowman SJ, Gordon T: Sjögren's Syndrome. Arthritis and Allied Conditions: A Textbook in Rheumatology. Edited by: Koopman W. 2005, Philadelphia, PA: Lippencott William & Wilkins, 2: 1681-1706. 15
Robinson CP, Brayer J, Yamachika S, Esch TR, Peck AB, Stewart CA, Peen E, Jonsson R, Humphreys-Beher MG: Transfer of human serum IgG to nonobese diabetic Igmu null mice reveals a role for autoantibodies in the loss of secretory function of exocrine tissues in Sjögren's syndrome. Proc Natl Acad Sci USA. 1998, 95: 7538-7543. 10.1073/pnas.95.13.7538.
Dawson LJ, Stanbury J, Venn N, Hasdimir B, Rogers SN, Smith PM: Antimuscarinic antibodies in primary Sjögren's syndrome reversibly inhibit the mechanism of fluid secretion by human submandibular salivary acinar cells. Arthritis Rheum. 2006, 54: 1165-1173. 10.1002/art.21764.
Gao J, Cha S, Jonsson R, Opalko J, Peck AB: Detection of anti-type 3 muscarinic acetylcholine receptor autoantibodies in the sera of Sjögren's syndrome patients by use of a transfected cell line assay. Arthritis Rheum. 2004, 50: 2615-2621. 10.1002/art.20371.
Delaleu N, Jonsson MV, Jonsson R: Disease mechanisms of Sjögren's syndrome. Drug Discovery Today. 2004, 1: 329-336. 10.1016/j.ddmec.2004.11.014.
Fox RI: Sjögren's syndrome. Lancet. 2005, 366: 321-331. 10.1016/S0140-6736(05)66990-5.
Cha S, Peck AB, Humphreys-Beher MG: Progress in understanding autoimmune exocrinopathy using the non-obese diabetic mouse: an update. Crit Rev Oral Biol Med. 2002, 13: 5-16.
Jonsson MV, Delaleu N, Jonsson R: Animal models of Sjögren's syndrome. Clin Rev Allergy Immunol. 2007, 32: 215-224. 10.1007/s12016-007-8012-7.
Jonsson MV, Delaleu N, Brokstad KA, Berggreen E, Skarstein K: Impaired salivary gland function in NOD mice: association with changes in cytokine profile but not with histopathologic changes in the salivary gland. Arthritis Rheum. 2006, 54: 2300-2305. 10.1002/art.21945.
Cha S, Nagashima H, Brown VB, Peck AB, Humphreys-Beher MG: Two NOD Idd-associated intervals contribute synergistically to the development of autoimmune exocrinopathy (Sjögren's syndrome) on a healthy murine background. Arthritis Rheum. 2002, 46: 1390-1398. 10.1002/art.10258.
Robinson CP, Yamachika S, Bounous DI, Brayer J, Jonsson R, Holmdahl R, Peck AB, Humphreys-Beher MG: A novel NOD-derived murine model of primary Sjögren's syndrome. Arthritis Rheum. 1998, 41: 150-156. 10.1002/1529-0131(199801)41:1<150::AID-ART18>3.0.CO;2-T.
Wicker LS, Appel MC, Dotta F, Pressey A, Miller BJ, DeLarato NH, Fischer PA, Boltz RC, Peterson LB: Autoimmune syndromes in major histocompatibility complex (MHC) congenic strains of nonobese diabetic (NOD) mice. The NOD MHC is dominant for insulitis and cyclophosphamide-induced diabetes. J Exp Med. 1992, 176: 67-77. 10.1084/jem.176.1.67.
Hayashi Y, Utsuyama M, Kurashima C, Hirokawa K: Spontaneous development of organ-specific autoimmune lesions in aged C57BL/6 mice. Clin Exp Immunol. 1989, 78: 120-126.
Mitsias DI, Tzioufas AG, Veiopoulou C, Zintzaras E, Tassios IK, Kogopoulou O, Moutsopoulos HM, Thyphronitis G: The Th1/Th2 cytokine balance changes with the progress of the immunopathological lesion of Sjögren's syndrome. Clin Exp Immunol. 2002, 128: 562-568. 10.1046/j.1365-2249.2002.01869.x.
Nguyen CQ, Gao JH, Kim H, Saban DR, Cornelius JG, Peck AB: IL-4-STAT6 signal transduction-dependent induction of the clinical phase of Sjögren's syndrome-like disease of the nonobese diabetic mouse. J Immunol. 2007, 179: 382-390.
Mehrad B, Keane MP, Strieter RM: Chemokines as mediators of angiogenesis. Thromb Haemost. 2007, 97: 755-762.
Leiter EH: The NOD mouse: as model for insulin-dependent diabetes mellitus. Current Protocols in Immunology. Edited by: Coligan JE, Kruisbeek AM, Margulies DH, Shevach EM, Strober W. 1997, John Wiley and Sons, Inc. Hoboken, NJ, USA, 15.19.11-15.19.23.
Gao J, Killedar S, Cornelius JG, Nguyen C, Cha S, Peck AB: Sjögren's syndrome in the NOD mouse model is an interleukin-4 time-dependent, antibody isotype-specific autoimmune disease. J Autoimmun. 2006, 26: 90-103. 10.1016/j.jaut.2005.11.004.
PANTHER Classification System. [http://www.pantherdb.org]
Hu S, Wang J, Meijer J, Ieong S, Xie Y, Yu T, Zhou H, Henry S, Vissink A, Pijpe J, Kallenberg C, Elashoff D, Loo JA, Wong DT: Salivary proteomic and genomic biomarkers for primary Sjögren's syndrome. Arthritis Rheum. 2007, 56: 3588-3600. 10.1002/art.22954.
Groom J, Kalled SL, Cutler AH, Olson C, Woodcock SA, Schneider P, Tschopp J, Cachero TG, Batten M, Wheway J, et al: Association of BAFF/BLyS overexpression and altered B cell differentiation with Sjögren's syndrome. J Clin Invest. 2002, 109: 59-68.
Båve U, Nordmark G, Lövgren T, Ronnelid J, Cajander S, Eloranta ML, Alm GV, Rönnblom L: Activation of the type I interferon system in primary Sjögren's syndrome: a possible etiopathogenic mechanism. Arthritis Rheum. 2005, 52: 1185-1195. 10.1002/art.20998.
Parambil JG, Myers JL, Lindell RM, Matteson EL, Ryu JH: Interstitial lung disease in primary Sjögren syndrome. Chest. 2006, 130: 1489-1495. 10.1378/chest.130.5.1489.
Tu WH, Shearn MA, Lee JC, Hopper J: Interstitial nephritis in Sjögren's syndrome. Ann Intern Med. 1968, 69: 1163-1170.
Moutsopoulos HM, Balow JE, Lawley TJ, Stahl NI, Antonovych TT, Chused TM: Immune complex glomerulonephritis in sicca syndrome. Am J Med. 1978, 64: 955-960. 10.1016/0002-9343(78)90449-7.
Hjelmervik TO, Lindqvist AK, Petersen K, Johannesson M, Stavrum AK, Johansson A, Jonsson R, Holmdahl R, Bolstad AI: The influence of the NOD Nss1/Idd5 loci on sialadenitis and gene expression in salivary glands of congenic mice. Arthritis Res Ther. 2007, 9: R99-10.1186/ar2300.
Larsen MR, Jensen SS, Jakobsen LA, Heegaard NH: Exploring the sialiome using titanium dioxide chromatography and mass spectrometry. Mol Cell Proteomics. 2007, 6: 1778-1787. 10.1074/mcp.M700086-MCP200.
Fulkerson PC, Zimmermann N, Hassman LM, Finkelman FD, Rothenberg ME: Pulmonary chemokine expression is coordinately regulated by STAT1, STAT6, and IFN-gamma. J Immunol. 2004, 173: 7565-7574.
Lloyd CM, Delaney T, Nguyen T, Tian J, Martinez AC, Coyle AJ, Gutierrez-Ramos JC: CC chemokine receptor (CCR)3/eotaxin is followed by CCR4/monocyte-derived chemokine in mediating pulmonary T helper lymphocyte type 2 recruitment after serial antigen challenge in vivo. J Exp Med. 2000, 191: 265-274. 10.1084/jem.191.2.265.
Lebre MC, Burwell T, Vieira PL, Lora J, Coyle AJ, Kapsenberg ML, Clausen BE, De Jong EC: Differential expression of inflammatory chemokines by Th1- and Th2-cell promoting dendritic cells: a role for different mature dendritic cell populations in attracting appropriate effector cells to peripheral sites of inflammation. Immunol Cell Biol. 2005, 83: 525-535. 10.1111/j.1440-1711.2005.01365.x.
Dimitriou ID, Kapsogeorgou EK, Moutsopoulos HM, Manoussakis MN: CD40 on salivary gland epithelial cells: high constitutive expression by cultured cells from Sjögren's syndrome patients indicating their intrinsic activation. Clin Exp Immunol. 2002, 127: 386-392. 10.1046/j.1365-2249.2002.01752.x.
Clark LB, Foy TM, Noelle RJ: CD40 and its ligand. Adv Immunol. 1996, 63: 43-78.
Baudino L, Azeredo da Silveira S, Nakata M, Izui S: Molecular and cellular basis for pathogenicity of autoantibodies: lessons from murine monoclonal autoantibodies. Springer Semin Immunopathol. 2006, 28: 175-184. 10.1007/s00281-006-0037-0.
Nguyen CQ, Kim H, Cornelius JG, Peck AB: Development of Sjögren's syndrome in nonobese diabetic-derived autoimmune-prone C57BL/6.NOD-Aec1Aec2 mice is dependent on complement component-3. J Immunol. 2007, 179: 2318-2329.
Quezada SA, Bennett K, Blazar BR, Rudensky AY, Sakaguchi S, Noelle RJ: Analysis of the underlying cellular mechanisms of anti-CD154-induced graft tolerance: the interplay of clonal anergy and immune regulation. J Immunol. 2005, 175: 771-779.
Peguet-Navarro J, Dalbiez-Gauthier C, Moulon C, Berthier O, Reano A, Gaucherand M, Banchereau J, Rousset F, Schmitt D: CD40 ligation of human keratinocytes inhibits their proliferation and induces their differentiation. J Immunol. 1997, 158: 144-152.
Kok MR, Yamano S, Lodde BM, Wang J, Couwenhoven RI, Yakar S, Voutetakis A, Leroith D, Schmidt M, Afione S, et al: Local adeno-associated virus-mediated interleukin 10 gene transfer has disease-modifying effects in a murine model of Sjögren's syndrome. Hum Gene Ther. 2003, 14: 1605-1618. 10.1089/104303403322542257.
Saito I, Haruta K, Shimuta M, Inoue H, Sakurai H, Yamada K, Ishimaru N, Higashiyama H, Sumida T, Ishida H, et al: Fas ligand-mediated exocrinopathy resembling Sjögren's syndrome in mice transgenic for IL-10. J Immunol. 1999, 162: 2488-2494.
Zachlederova M, Jarolim P: Gene expression profiles of microvascular endothelial cells after stimuli implicated in the pathogenesis of vasoocclusion. Blood Cells Mol Dis. 2003, 30: 70-81. 10.1016/S1079-9796(03)00011-1.
Okamoto H, Yamamura M, Morita Y, Harada S, Makino H, Ota Z: The synovial expression and serum levels of interleukin-6, interleukin-11, leukemia inhibitory factor, and oncostatin M in rheumatoid arthritis. Arthritis Rheum. 1997, 40: 1096-1105. 10.1002/art.1780400614.
Olson TS, Ley K: Chemokines and chemokine receptors in leukocyte trafficking. Am J Physiol Regul Integr Comp Physiol. 2002, 283: R7-R28.
Strieter RM, Gomperts BN, Keane MP: The role of CXC chemokines in pulmonary fibrosis. J Clin Invest. 2007, 117: 549-556. 10.1172/JCI30562.
Bodolay E, Koch AE, Kim J, Szegedi G, Szekanecz Z: Angiogenesis and chemokines in rheumatoid arthritis and other systemic inflammatory rheumatic diseases. J Cell Mol Med. 2002, 6: 357-376. 10.1111/j.1582-4934.2002.tb00514.x.
Delaleu N, Madureira AC, Immervoll H, Jonsson R: Inhibition of experimental Sjögren's syndrome through immunization with heat-shock protein 60 kDa and its peptide aa437–460. Arthritis Rheum.
Funding for this study was provided by The Meltzer Fond, Olaf and Gullborg Johannessens Legacy, the Broegelmann Foundation and the Strategic Research Program, and Western Norway Regional Health Authority. We thank Ana Carina Madureira, Gry Bernes and the staff at the Animal Facility, Department of Physiology, University of Bergen and at the Department of Pathology, Haukeland University Hospital, Bergen, Norway for excellent technical assistance. In addition, we thank Inge Jonassen from the Bergen Center of Computational Biology for helpful advice, Dr B Delaleu-Justitz for careful revision of the manuscript, and Dr S Küster for her support.
The authors declare that they have no competing interests.
ND and RJ designed the study. ND carried out the animal experiments, salivary flow measurements and collection of organs. HI participated in preparation and carried out the qualitative evaluation of the organs. ND carried out the morphometrical analyses. JC and ND carried out the fluorescence-activated cell sorting analyses for anti-M3R autoanatibodies. ND carried out the data analyses. ND and RJ wrote the manuscript. RJ supervised the study. All authors read and approved the final manuscript.