|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Proteomics and Protein Markers |
1 Laboratory of Rheumatology, GIGA Research, CHU, University of Liège; 2 Laboratory of Clinical Chemistry, GIGA Research, University of Liège; 3 GIGA Bioinformatics Platform, University of Liège; 4 Bioinformatics and Modeling Unit, Department of Electrical Engineering & Computer Science, GIGA Research, University of Liège;5 Laboratory of Hepato-Gastroenterology, CHU, University of Liège, Liège, Belgium.
aAddress correspondence to this author at: Laboratory of Rheumatology, Tour GIGA +2, CHU, 4000 Liège, Belgium. Fax (32) 4 366 45 34; e-mail ddeseny{at}chu.ulg.ac.be
| Abstract |
|---|
|
|
|---|
Methods: We used SELDI-TOF MS to analyze serum samples from patients with various forms of inflammatory arthritis. Several protein profiles were collected on different Bio-Rad Laboratories ProteinChip arrays (CM10 and IMAC-Cu2+) and were evaluated statistically to select potential biomarkers.
Results: SELDI-TOF MS analyses identified several calgranulin proteins [S100A8 (calgranulin A), S100A9 (calgranulin B), S100A9*, and S100A12 (calgranulin C)], serum amyloid A (SAA), SAA des-Arg (SAA-R), and SAA des-Arg/des-Ser (SAA-RS) as biomarkers and confirmed the results with other techniques, such as western blotting, immunoprecipitation, and nano-LC-MS/MS. The S100 proteins were all able to significantly differentiate samples from patients with rheumatoid arthritis (RA), psoriatic arthritis (PsA), and ankylosing spondylitis (AS) from those of patients with inflammatory bowel diseases used as an inflammatory control (IC) group, whereas the SAA, SAA-R, and SAA-RS proteins were not, with the exception of AS. The 4 S100 proteins were coproduced in all of the pathologies and were significantly correlated with the plasma calprotectin concentration; however, these S100 proteins were correlated with the SAA peak intensities only in the RA and IC patient groups. In RA, these S100 proteins (except for S100A12) were significantly correlated with the serum concentrations of C-reactive protein, matrix metalloproteinase 3, and anti–cyclic citrullinated peptide and with the Disease Activity Score (DAS28).
Conclusions: The SELDI-TOF MS technology is a powerful approach for analyzing the status of monomeric, truncated, or posttranslationally modified forms of arthritis biomarkers, such as the S100A8, S100A9, S100A12, and SAA proteins. The fact that the SELDI-TOF MS data were correlated with results obtained with the classic calprotectin ELISA test supports the reliability of this new proteomic technique.
| Introduction |
|---|
|
|
|---|
S100A12 is produced mainly by granulocytes upon inflammatory activation and acts independently of S100A8 and S100A9(17). Its interaction with the receptor for advanced glycation end products (RAGE) induces proinflammatory signals in endothelium and in cells of the immune system(18). Increased S100A12 concentrations have been found in the serum, synovium, and synovial fluid of patients with RA(11)(19), PsA(19), and AS(19), as well as in the serum of IBD patients(20). Serum amyloid A (SAA), a key promoter of inflammatory events in RA and shown to be produced by inflamed synovial tissue, also induces RAGE activation(21).
The identification of the S100 and SAA groups of proteins as biomarkers is therefore challenging. In the absence of ELISAs specific for the monomeric forms of S100A8 and S100A9 and of a commercially available test for S100A12, we have hypothesized that mass spectrometry (MS) might be a possible approach for detecting these proteins(22). S100A8 has already been identified by 2-dimensional gel electrophoresis to be present in synovium(13), but not in serum(12). Two-dimensional liquid chromatography–coupled tandem MS has also been used successfully to identify several S100 proteins, including S100A8, in serum samples from a few RA patients(11); however, the labor-intensive nature of these 2 proteomics technologies allow the investigation of only small numbers of biological samples. Accordingly, we postulated SELDI-TOF MS technology to be a more appropriate proteomic approach for detecting low molecular weight proteins (<20 kDa), because such an approach permits rapid analysis of hundreds of serum samples at a time(23). The potential of this technology for discovering biomarkers has been demonstrated for such chronic inflammatory conditions as RA and PsA(24), as well as for IBD(25). We therefore investigated S100A8, S100A9, S100A12, and SAA status in RA, PsA, and AS patients, with IBD patients used as a positive inflammatory control (IC) group. We then evaluated the data generated for these markers with respect to the plasma calprotectin concentration and several clinical and biological variables associated with arthritis activity.
| Patients and Methods |
|---|
|
|
|---|
The epidemiologic characteristics of the patient groups are summarized in Table 1
. RA patients fulfilling the 1987 American College of Rheumatology criteria(26) had a median DAS28 of 6.3 (range, 3.5–8.8), with 86% of the scores >5.1 (high disease activity). PsA patients had active disease with at least 3 tender and swollen joints. The AS patients [modified New York criteria(27)] had active disease, as indicated by a median Bath Ankylosing Spondylitis Disease Activity Index (BASDAI) of 6/10 (range, 4/10–10/10). Diagnosis of IBD patients was made according to validated criteria(28). Active Crohn disease was defined by a Harvey–Bradshaw index
7, and active ulcerative colitis was defined by both clinical and endoscopic signs of activity.
|
The preparation and reproducibility of ProteinChip arrays (Vermillion/Ciphergen Biosystems) are described in the supplemental methods in the Data Supplement that accompanies the online version of this article at http://www.clinchem.org/content/vol54/issue6.
analysis of seldi-tof ms data
Peaks were detected with ProteinChip Biomarker Wizard software (version 3.0; Bio-Rad Laboratories). We analyzed the data by 2 approaches, nonparametric Mann–Whitney U tests and a machine-learning algorithm called random forests(29), after we had completed various preprocessing steps(24). Random forests is a decision-tree multivariate analysis that estimates the relevance or relative contribution of each peak to the classification of 2 groups(29). The latter approach allows m/z values to be ranked according to their relevance for differentiating the 2 groups based on quantitative estimates of the percentage of information (% of info) supplied. P values <0.05 were considered to be statistically significant.
western blotting
S100A8, S100A9, S100A12, and SAA were assessed by Western blot (WB) analysis. In brief, 2 µL of serum was run on 12% NuPAGE Bis-Tris polyacrylamide gels (Invitrogen), transferred, and incubated with anti-S100A8 monoclonal antibody (20 µL diluted in 10 mL; BMA Biomedicals), anti-S100A9 polyclonal antibody (10 µL diluted in 10 mL; Santa Cruz Biotechnology), anti-S100A12 polyclonal antibody (20 µL diluted in 10 mL: Santa Cruz Biotechnology), or anti-SAA monoclonal antibody (5 µL diluted in 10 mL; Abcam). We then incubated with a mouse secondary antibody (2 µL diluted in 10 mL; GE Healthcare/Amersham Biosciences) to specifically detect S100A8 and SAA, with a rabbit secondary antibody (2 µL diluted in 10 mL; Dako) to detect S100A9, or with a goat secondary antibody (2 µL diluted in 10 mL; Santa Cruz Biotechnology) to detect S100A12. Proteins were revealed with an enhanced chemiluminescence detection method according to the manufacturers instructions (GE Healthcare/Amersham Biosciences).
The supplemental methods in the online Data Supplement describe the immunoprecipitation and 1-dimensional gel electrophoresis methods as well as LC-MS/MS identification of SAA.
elisa tests
We carried out ELISAs for anti–cyclic citrullinated peptide 2 (anti-CCP2) antibody (cutoff, 5 000 relative units/L; Euroimmun), calprotectin (range, 1.6–100 µg/L; Hycult biotechnology), and matrix metalloproteinase 3 (MMP-3) (range, 1.25–20 µg/L; BioSource)(30) as recommended by the respective manufacturers. We used serum samples or, in the case of calprotectin assays, the corresponding plasma samples.
correlation analysis
We used ranked data to calculate Spearman correlation coefficients and used the
2 test to compare qualitative data. P values <0.05 were considered statistically significant.
| Results |
|---|
|
|
|---|
|
We also loaded the 139 serum samples in triplicate onto immobilized metal affinity capture ProteinChip (IMAC-Cu2+) arrays and collected 417 spectra. Both statistical approaches detected 3 peaks as significantly increased in RA patients compared with the NIC group. These peaks had m/z values of 11 438 (P < 10–6; % of info, 1.2%; rank, 20), 11 528 (P < 10–9; % of info, 3.7%; rank, 5), and 11 680 (P < 10–9; % of info, 6.5%; rank, 1). These peaks were thought to be related to 3 SAA variant proteins: SAA (calculated mass, 11 682 Da), SAA des-Arg (SAA-R) (calculated mass, 11 526 Da) and SAA des-Arg/des-Ser (SAA-RS) (calculated mass, 11 439 Da). Fig. 1
shows these potential biomarkers on spectra collected on CM10 and IMAC-Cu2+ arrays of serum samples from patients with various forms of arthritis and from NIC group individuals.
|
identification of s100 family proteins
We confirmed the identities of the S100A8, S100A9, and S100A12 proteins by comparing WB and SELDI-TOF MS results for the serum samples (Fig. 2
). The WB analysis was performed with 8 serum samples (positions 1–8) from each of the 5 sets of patients (NIC, RA, PsA, AS, and IC). The same RA sample was run in position 9 on each gel as a positive control. SELDI-TOF MS spectra for the same serum samples were monitored at m/z values of 10 835, 13 272, and 10 444 and compared with the corresponding WB results. We obtained similar profiles in the WB and SELDI analyses for each protein tested (Fig. 2
, A–C). These data suggest that the peaks at 10 835, 13 272, and 10 444 m/z values are the S100A8, S100A9, and S100A12 proteins, respectively. According to these 2 semiquantitative approaches, some RA, PsA, AS, and, to a lesser extent, IC serum samples were positive for the 3 S100 proteins. We detected 2 S100A9 variants. The m/z value of 13 272 agreed well with the calculated mass for an oxidized form of S100A9 (13 242 Da plus 32 Da), and the 12 688 m/z value likely represents an S100A9 variant known as S100A9*, which has previously been characterized by ultraviolet MALDI MS(31). S100A9* results from translation beginning at amino acid residue 5 and acetylation of amino acid residue 6 (a Ser residue), yielding a calculated mass of 12 691 Da.
|
We also confirmed the identities of the S100A8 and S100A9 proteins by eluting the 2 proteins from their immunocomplexes on NP20 arrays after immunodepletion of a serum sample from an RA patient (see Fig. 1
in the online Data Supplement). We also detected a slight cross-reactivity with S100A12, S100A9*, and S100A9 on the NP20 spectra after eluting S100A8 (the sequence homologies of the S100A8, S100A9, and S100A12 proteins are around 40%)(32). We were unable to immunoprecipitate the S100A12 protein.
identification of saa proteins
To identify the proteins responsible for the 11 438, 11 528, and 11 680 m/z peak values in the RA spectra on IMAC-Cu2+ arrays, we collected serum samples from a healthy control individual and an RA patient, depleted the samples of albumin and IgG, and ran them on a 1-dimensional SDS-PAGE gel (see Fig. 2 in the online Data Supplement). After silver staining, we excised a band from the RA sample on the gel with an apparent molecular weight of 11 kDa and subjected it to LC-MS/MS analysis (see Fig. 2 in the online Data Supplement). We also excised another band from the gel at the same position for the serum sample from the healthy control individual. We digested these bands with trypsin and analyzed the resulting peptides by tandem MS (see Fig. 2 in the online Data Supplement). Sequencing analysis of 4 major tryptic fragments (sequence coverage, 51%; total score, 186) revealed the excised protein to be SAA.
We confirmed the identity of SAA by comparing the WB results with the SELDI-TOF MS results for the serum samples (Fig. 2D
). Each set of patient samples (but not those of the NIC group) was positive for SAA in the WB analysis. This result confirmed the identity of the peak at m/z 11 680 to be the SAA protein (calculated mass, 11 682 Da). The peaks at m/z values of 11 438, 11 528, and 11 680 were clustered on each spectrum. We therefore hypothesized that the 11 438 and 11 528 m/z peak values represented 2 variants of the original SAA protein. The 11 528 m/z value corresponds to the calculated mass for the SAA protein without its first N-terminal Arg residue (–156 Da); this peak represents the SAA-R protein (calculated mass, 11 526 Da). Similarly, the 11 438 m/z value corresponds to the calculated mass of the SAA protein truncated at the N-terminal end by 2 residues, Arg and Ser (–243 Da); this peak represents the SAA-RS protein (calculated mass, 11 439 Da). We also confirmed the identities of these proteins by eluting them from their immunocomplexes on NP20 after immunodepleting a serum sample from an RA patient (see Fig. 1 in the online Data Supplement).
the distribution of peak intensities among patient groups
The discriminatory power of the peaks at m/z values of 10 835 (S100A8), 12 688 (S100A9*), 13 272 (S100A9), 10 444 (S100A12), 11 680 (SAA), 11 528 (SAA-R), and 11 438 (SAA-RS) detected in the spectra of samples from each of the arthritis groups was assessed with the Mann–Whitney U test (Table 2
). In brief, S100 proteins were found to be highly significantly effective (P < 10–9) for distinguishing RA, PsA, and AS patients, but not IC patients, from NIC individuals. Similarly, the SAA, SAA-R, and SAA-RS proteins were effective for distinguishing RA, PsA, and IC patients (but not AS patients) from NIC individuals. All of the S100 proteins significantly distinguished RA, PsA, and AS patients from individuals in the IC group, whereas the SAA, SAA-R, and SAA-RS proteins did not, with the exception of AS patients. Lastly, a comparison within the rheumatic inflammatory disease groups revealed that the S100 proteins weakly distinguished PsA and AS patients, whereas the S100A12 protein was the only variable that significantly distinguished RA patients from PsA patients. SAA, SAA-R, and SAA-RS proteins uniformly and weakly distinguished RA from AS patients.
We next analyzed the percentages of serum samples that were positive for the 7 biomarkers. These percentages were calculated according to the peak intensities measured at m/z values of 10 444 (S100A12), 10 835 (S100A8), 13 272 (S100A9), and 12 688 (S100A9*) on CM10 arrays, and at m/z values of 11 680 (SAA), 11 528 (SAA-R), and 11 438 (SAA-RS) on IMAC-Cu2+ arrays. Peak intensity values were averaged. The cutoff value was defined as the highest peak intensity of the corresponding m/z values in the spectra for the NIC serum samples. Similarly, the percentages of plasma samples positive for calprotectin were calculated with the cutoff defined as the highest calprotectin concentration observed in the NIC group. Increased S100 protein peak intensities and increased calprotectin concentrations were detected in 43%–89% and 33%–56% of arthritis serum samples, respectively, and in 15%–34% and 33% of IC patients, respectively (Table 3
). The positivity rates among the 4 S100 proteins were significantly linked [
2 (9) = 84; P < 0.0001] and were significantly linked with the positivity rates for calprotectin [
2 (7) = 38; P < 0.0001]. Increased intensities for SAA proteins were detected in 29%–44% of RA patients, 22% of PsA patients, 16% of AS patients, and 33%–42% of IC patients. No significant differences were found between SAA, SAA-R, and SAA-RS intensities in any of the sample sets. The positivities of these variants were also significantly linked [
2 (3) = 113; P < 0.0001].
|
Statistically significant correlations were found between the various S100 proteins, SAA, calprotectin, and other evaluated variables in the RA, PsA, AS, and IC patient groups (Table 4
). In brief, values for the 4 S100 proteins were both highly intercorrelated and correlated with plasma calprotectin concentration in all of the patient groups. S100 protein peaks were also correlated with SAA peaks, but only in the RA and IC groups. In the RA group, S100 protein peaks and the calprotectin concentration were correlated with CRP, log MMP-3, and anti-CCP2 antibody (except for S100A12) serum concentrations and with the DAS28, but not with the number of tender or swollen joints (data not shown). RA patients who produced S100 proteins or calprotectin had longer disease durations than S100-negative and calprotectin-negative patients (132 months vs 17 months and 117 months vs 16 months, respectively; P < 0.05). S100 protein peaks were not correlated in PsA or AS patient groups with the CRP, SAA, or log MMP-3 (except for S100A12 in the PsA group) serum concentration, or with clinical variables, including the number of swollen or tender joints and the BASDAI (data not shown). S100 protein peaks were also correlated with log MMP-3 serum concentration in the IC patient group.
|
SAA was correlated with the serum concentration of CRP in all patient groups and with DAS28 and logarithm MMP-3 serum concentration in the RA patient group. The plasma calprotectin concentration was correlated with log MMP-3 serum concentration in the RA, PsA, and IC patient groups and with the SAA protein in the RA, AS, and IC patient groups.
| Discussion |
|---|
|
|
|---|
The highly significant linear correlations between the peak intensities for the S100 proteins and the calprotectin plasma concentration as measured by ELISA confirm the reliability of this new proteomic approach for investigating these inflammation-related proteins. We conclude that these proteins are up-regulated in each of the diseases we studied because the S100 proteins and calprotectin are correlated both qualitatively, as shown by the strong concordance between WB positivity and SELDI-TOF MS positivity, and quantitatively, as shown by the significant linear correlations in peak intensities. Furthermore, the peak intensities of S100A8, S100A9, S100A9*, and S100A12, as well as the calprotectin concentration, are correlated with variables that reflect the biological and clinical activity of RA, such as the DAS28 and serum concentrations of CRP, MMP-3, and anti-CCP2 antibody. These results provide strong support for the clinical relevance of proteomic detection of calgranulins. The anti-CCP2 antibody may also be related to radiologically observed structural damage in RA because anti-CCP antibodies have been demonstrated to be independent predictors of joint damage(35), as is the serum concentration of calprotectin(16). We could not address this question directly because we conducted no x-ray–imaging studies; however, we did find that patients with increased S100 protein or calprotectin concentrations had significantly longer disease durations, which are expected to be associated with greater radiographically detectable damage.
We observed a discrepancy between RA patients and PsA and AS patients in that the 4 S100 proteins were significantly correlated with CRP and MMP-3 concentrations in RA patients, but not in the other 2 groups of patients. We attribute this finding to the fact that AS and PsA patients exhibit abnormal CRP and MMP-3 concentrations less frequently than RA patients(30). The strong correlation between the 4 S100 proteins and the serum concentration of MMP-3 in the IC group agrees with findings of abundant production of both types of proteins in inflamed intestinal tissue(36).
We also observed that the 4 S100 proteins distinguished arthritis conditions from nonarthritis conditions (NIC and IC groups), whereas SAA patterns distinguished inflammatory diseases (RA, PsA, and IC groups) from noninflammatory conditions (the NIC group). It is also interesting that within the arthritis groups the peak intensities for the 4 S100 proteins were correlated with the SAA peak intensities in the RA group but not in the AS and PsA groups. These results mimic what we had already observed with the CRP and MMP-3 variables. We therefore conclude that the arthritis process is correlated with the inflammatory process in RA, in which the acute-phase response is well developed, but not in AS or PsA, in which it is weaker. This conclusion thus suggests that the regulatory pathways at the sites of inflammation and the patterns of local production of S100 and SAA proteins in RA patients are different from those in AS and PsA patients(37).
SELDI-TOF MS is a unique method for distinguishing S100 monomers from multimeric forms and for detecting SAA variants, which are not possible with current ELISAs. The identification of the various SAA forms may be important because they may have different pathophysiological roles. The recent identification via SELDI-TOF analysis of several truncated forms of S100A8 and S100A12 in cystic fibrosis patients suggests that C-terminal truncations affect protein function(34). The presence of S100A8 and S100A12, but not calprotectin, have been found to be characteristic of intra-amniotic inflammation(22), and SAA variants with different properties, such as differential susceptibility to matrix metalloproteinase digestion, have been described(38). These findings demonstrate the relevance of developing new proteomic approaches sufficiently powerful for investigating proteins in their modified or monomeric forms.
In conclusion, we have used SELDI-TOF MS technology to identify several relevant arthritis biomarkers that are correlated with several biological or clinical variables associated with disease activity. We could not address the functional role of these biomarkers in the pathophysiology of arthritis, because we selected serum samples from individuals with the respective disease characteristics and not from individuals at different stages of each disease. Studies that take advantage of the SELDI-TOF MS technology to evaluate the effects of disease stage on these biomarkers may shed light on this question.
| Acknowledgments |
|---|
Financial Disclosures: None declared.
Acknowledgments: The authors thank Aline Desoroux and Gaël Cobraiville for their expert technical assistance. M.F. is a Research Associate, and E.L. and M.P.M. are Senior Research Associates at FNRS (National Fund for Scientific Research).
| Footnotes |
|---|
| References |
|---|
|
|
|---|
treatment of infiltrating versus resident tissue macrophage subsets in autoimmune arthritis. J Pathol 2005;206:17-27.[CrossRef][Web of Science][Medline]
[Order article via Infotrieve]
-helices can determine specific association of two EF-hand proteins. J Mol Biol 2007;370:887-898.[CrossRef][Web of Science][Medline]
[Order article via Infotrieve]
B-dependent signal transduction pathway. Arthritis Rheum 2006;54:105-114.[CrossRef][Web of Science][Medline]
[Order article via Infotrieve]The following articles in journals at HighWire Press have cited this article:
![]() |
A. Baillet, C. Trocme, S. Berthier, M. Arlotto, L. Grange, J. Chenau, S. Quetant, M. Seve, F. Berger, R. Juvin, et al. Synovial fluid proteomic fingerprint: S100A8, S100A9 and S100A12 proteins discriminate rheumatoid arthritis from other inflammatory joint diseases Rheumatology, January 25, 2010; (2010): kep452v1 - kep452. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Chung, A. E. Nelson, K. K. Y. Ho, and R. C. Baxter Proteomic Profiling of Growth Hormone-Responsive Proteins in Human Peripheral Blood Leukocytes J. Clin. Endocrinol. Metab., August 1, 2009; 94(8): 3038 - 3043. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |