Participants
The analysis was performed in two steps (additional file 1). In the first step, the differential salivary protein composition between patients with SLE and healthy controls (HCs) were analyzed by 2-DE with MS. In the next step, the differentially expressed proteins identified in the 2-DE with MS analysis were validated by western blotting and enzyme-linked immunosorbent assay (ELISA).
There were two participant groups in this study: The first group included 11 patients with SLE and 11 HCs, whose salivary samples (the samples were pooled for each group) were subjected to 2-DE with MS proteomic analysis. The second participant group included 94 patients with SLE, 57 patients with RA, and 62 HCs (Table 1). The concentration of proteins in the saliva samples of HCs and patients with SLE or RA was validated by western blotting and ELISA.
All enrolled patients with SLE met the revised American College of Rheumatology classification criteria [19]. Additionally, the age-matched and sex-matched patients with RA and HCs who had no history of autoimmune or inflammatory disorders were enrolled in the study. This study was conducted according to the Declaration of Helsinki and Good Clinical Practice guidelines. All subjects provided their informed consents for participating in the study. The study protocol was approved by the institutional review of board of our hospital (BMR-SMP-13-199).
The information on medical history and clinical manifestations was collected from a chart review and blood test results, such as complete blood count, erythrocyte sedimentation rate (ESR) and the levels of anti-nuclear antibody, complement 3 and 4, and anti-dsDNA antibody. The patients with RA were enrolled as a disease control to analyze the differential expression of specific proteins between SLE and RA, which are both chronic autoimmune diseases.
Basic characteristics of second participant group is here: the mean age of patients with SLE, patients with RA, and HCs was 39.8 ± 9.8, 41 ± 7.9, and 39.5 ± 6.9 years, respectively. The mean age was not different between the three groups. Among the patients with SLE, 41 patients (43.6%) tested positive for anti-dsDNA antibody, 28 patients (29.8%) had mucocutaneous symptoms, 31 patients (33.0%) had arthritis, and 29 patients (30.9%) had nephritis. The mean SLE disease activity index (SLEDAI) was 3.8 ± 4.2. Among the 57 patients with RA, 45 patients (77.6%) tested positive for rheumatoid factors and their mean disease activity score including 28 joints (DAS28) was 3.3 ± 1.15.
Table 1 Clinical characteristics of the subjects
|
SLE
|
RA
|
HC
|
Number
|
94
|
57
|
62
|
Age, years
|
39.8 ± 9.8
|
41 ± 7.9
|
39.5 ± 6.9
|
Sex (F/M)
|
87/7
|
50/7
|
58/4
|
Leukocyte, /µL
|
5165 ± 2365
|
|
|
Hemoglobin, /µL
|
12.2 ± 2.4
|
|
|
Platelet, x 103/µL
|
221.8 ± 76.1
|
|
|
Lymphocyte, /µL
|
1454 ± 654
|
|
|
ESR, mm/h
|
16.4 ± 18
|
|
|
Complement 3, mg/dL
|
85.2 ± 26.9
|
|
|
Complement 4, mg/dL
|
18.6 ± 8.9
|
|
|
Anti-dsDNA Ab (+), n (%)
|
41 (43.6)
|
|
|
Rheumatoid factor (+), n (%)
|
|
45 (77.6)
|
|
Mucocutaneous involvement, n (%)
|
28 (29.8)
|
|
|
Arthritis, n (%)
|
31 (33)
|
|
|
Nephritis, n (%)
|
29 (30.9)
|
|
|
Serositis, n (%)
|
4 (3.8)
|
|
|
Hematologic involvement, n (%)
|
35 (37.6)
|
|
|
SLEDAI
|
3.8 ± 4.2
|
|
|
DAS-28
|
|
3.3 ± 1.2
|
|
All values presented as number (%) or mean ± standard deviation
Ab, antibody; DAS28, disease activity score including 28 joints; dsDNA, double-strand deoxyribonucleic acid; ESR, erythrocyte sedimentation rate; HC, healthy control; RA, rheumatoid arthritis; SLE, systemic lupus erythematosus; SLEDAI, systemic lupus erythematosus disease activity index
Saliva sample collection
As salivary proteins exhibit diurnal variations, the saliva samples were collected from all participants between 9:00 and 11:00 am. The subjects were not allowed to eat, drink, smoke, or perform oral hygiene procedures for at least 1 h prior to the sample collection. The saliva samples were collected for 5 min after the subjects rinsed their mouth with water [20]. The saliva secretion was not stimulated in the study subjects. The subjects were asked to keep their mouths closed and expectorate the saliva into a tube once per minute. Each saliva sample was immediately treated with the protease inhibitors to preserve the integrity of the protein constituents. The saliva samples were centrifuged at 3,000 rpm for 15 min at 4°C. After removing the clear supernatant, the samples were aliquoted and stored at -20°C until further use.
Two-dimensional gel electrophoresis (2-DE)
The samples from 11 patients with SLE or 11 HCs were pooled equally to avoid intra-class variations that were detected between the patients in 2-DE analyses. A 1 mL aliquot of the sample was concentrated 10 times using the Amicon-3K centrifugal filters at 14,000 g and 4°C for 20 min. The proteins in the salivary samples were precipitated using 500 μL of trichloroacetic acid/acetone (90%; v/v)-dithiothreitol mixture overnight at -20°C. The samples were centrifuged at 10,000 rpm and 10°C for 10 min. The supernatant was collected and the samples were pretreated with 250 μL of rehydration buffer. Next, the samples were centrifuged at 10,000 rpm and 10°C for 10 min to remove any insoluble material. The protein concentration of the samples was estimated by Bradford protein assay (Bio-Rad, Hercules, CA, USA).
Liquid chromatography tandem mass spectrometry (LC-MS)
With 2-DE proteomic analysis of the saliva samples of 11 patients with SLE and 11 HCs, the proteins separated into numerous spots with different concentrations. The proteins in 20 spots were subjected to liquid chromatography tandem-mass spectrometry (LC-MS) to analyze the proteins with high specificity [21]. The gel pieces containing the protein spots were destained, reduced, alkylated, and digested with modified sequencing grade trypsin (Sigma, MO, US A), as previously described [22]. Peptide mixtures were lyophilized and stored at -80°C for further LC-MS analysis.
Each sample was resuspended in 0.1% trifluoroacetic acid and injected into a Zorbox 300SB-C18 75 μm × 15 cm column (Agilent, Santa Clara, CA, USA) via the trap column. The peptides were separated in an acetonitrile gradient at a flow rate of 200 nL/min in an UltiMate 3000 nano HPLC system (Dionex, Sunnyvale, CA, USA). The peptides were then applied on-line to an LTQ (Thermo Fisher, Waltham, MA, USA) ion-trap mass spectrometer. The mobile phase gradient was initiated with an increase from 5 to 40% buffer within 110 min. Next, the gradient was increased to 80% buffer in 1 min. The gradient was maintained at isocratic conditions of 80% buffer for 15 min. The main working liquid-junction electrospray ionization source parameters were as follows: ion spray voltage, 1.6 kV; capillary voltage, 24 V; and capillary temperature, 200°C. The spectra were obtained in full scan mode using the dynamic exclusion criteria. LC-MS runs were analyzed using the DeCyder MS (version 2.0; GE Healthcare, Uppsala, Sweden) software.23 The peptide peaks were detected with an average peak width of 1 min and matched with a mass accuracy of at least 0.6 Da and a maximum time window of 4 min. The abundance of individual peptides in the respective gradient fraction was calculated by peak integration.
The data were manually examined, and the overlapping peaks were discarded (additional file 2). The threshold level for differentially expressed proteins was defined as at least 2-fold increase or decrease in spot intensity that was statistically significant. The MS spectra of the peptide peaks were searched against the Uniprot Human database using the Mascot version 2.3 (Matrix Science, London, UK). For quantitative protein profiling, only the proteins identified by multiple peptides with significant MASCOT score (p < 0.05) were considered.
Western blotting analysis
The immunoglobulin gamma-3 chain C region (IGHG3), immunoglobulin alpha-1 chain C region (IGHA1), protein S100-A8 (S100A8), lactoferrin, and 8-oxoguanine DNA glycosylase (OGG1) were analyzed by western blotting using the rabbit anti-human IGHG3 polyclonal (MBS248789, MyBiosource, San Diego, CA, USA), rabbit anti-human IGHA1 polyclonal (MBS9206028, MyBiosource), rabbit anti-human rat S100A8 polyclonal (MBS127619, MyBiosource), mouse anti-human lactoferrin monoclonal (ab10110, Abcam, Cambridge, UK), and rabbit anti-human OGG1 polyclonal (NB100-106, Novusbio, Centennial, CO, USA) antibodies, respectively. The proteins were subjected polyacrylamide gel electrophoresis using 10% (for IGHG3) or 15% (for IGHA1, S100/A8, lactoferrin, and OGG1) gel. The resolved proteins were transferred to a polyvinylidene fluoride membrane. The membranes were incubated with secondary antibody (goat anti-rabbit antibody A120-101P for IGHG3, IGHA1, S100A8, and OGG1, and goat anti-mouse antibody for lactoferrin, Bethyl Laboratories, Mongomery, TX, USA) diluted at 1:10,000 (IGHG3 and S100/A8), and 1:2,000 (IGHA1, lactoferrin, and OGG). All analyses were performed in triplicates. The protein concentration was determined by the optical density of specific immunoreactive bands, and the optical density of each bands was measured using image J software (NIH, Bethesda, MD, USA).
Enzyme-linked immunosorbent assay (ELISA)
The levels of salivary IGHG3 and lactoferrin were measured in patients with SLE or RA and HC by ELISA using the human IGHG3 ELISA kit (ab137981, Abcam) and human lactoferrin ELISA kit (ab108882, Abcam), respectively, following the manufacturer’s instructions. All measurements were performed in duplicates.
Statistical analysis
The difference in the expression of salivary IGHG3, IGHA1, S100A8, lactoferrin, and OGG1 determined by western blotting and the concentrations of salivary IGHG3 and lactoferrin measured by ELISA in patients with SLE or RA and HCs were compared using the two-sample Wilcoxon rank-sum (Mann-Whitney) test. As the data were not normally distributed and the variance was not homogeneous among the groups, one-way analysis of variance (ANOVA) was not used. Hence, post-hoc test with Bonferroni correction was used. The correlations between the levels of salivary IGHG3 or lactoferrin and clinical features in patients with SLE were determined using the Spearman’s rank correlation coefficient. On the receiver operating characteristic (ROC) curve analysis of salivary proteins, area under curve (AUC), sensitivity, and specificity were calculated. The difference was considered statistically significant when the p value was less than 0.05. All statistical analyses were performed in the Statistical Package for the Social Sciences version 22.0 (IBS Corp, Armonk, NY, USA) and SAS9.4 (SAS institute Inc, Cary, NC, USA).