Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Identification of Urinary Peptide Biomarkers Associated with Rheumatoid Arthritis

  • Angelique Stalmach ,

    Contributed equally to this work with: Angelique Stalmach, Hanna Johnsson, Duncan Porter

    Angelique.Stalmach@glasgow.ac.uk

    Affiliation University of Glasgow, Institute of Cardiovascular and Medical Sciences, Glasgow, United Kingdom

  • Hanna Johnsson ,

    Contributed equally to this work with: Angelique Stalmach, Hanna Johnsson, Duncan Porter

    Affiliation Western General Infirmary, Edinburgh, United Kingdom

  • Iain B. McInnes,

    Affiliation University of Glasgow, Institute of Infection, Immunity and Inflammation, Glasgow, United Kingdom

  • Holger Husi,

    Affiliation University of Glasgow, Institute of Cardiovascular and Medical Sciences, Glasgow, United Kingdom

  • Julie Klein,

    Affiliation Mosaiques-Diagnostics GmbH, Hannover, Germany

  • Mohammed Dakna,

    Affiliation Mosaiques-Diagnostics GmbH, Hannover, Germany

  • William Mullen,

    Affiliation University of Glasgow, Institute of Cardiovascular and Medical Sciences, Glasgow, United Kingdom

  • Harald Mischak,

    Affiliations University of Glasgow, Institute of Cardiovascular and Medical Sciences, Glasgow, United Kingdom, Mosaiques-Diagnostics GmbH, Hannover, Germany

  • Duncan Porter

    Contributed equally to this work with: Angelique Stalmach, Hanna Johnsson, Duncan Porter

    Affiliation University of Glasgow, Institute of Infection, Immunity and Inflammation, Glasgow, United Kingdom

Abstract

Early diagnosis and treatment of rheumatoid arthritis are associated with improved outcomes but current diagnostic tools such as rheumatoid factor or anti-citrullinated protein antibodies have shown limited sensitivity. In this pilot study we set out to establish a panel of urinary biomarkers associated with rheumatoid arthritis using capillary electrophoresis coupled to mass spectrometry. We compared the urinary proteome of 33 participants of the Scottish Early Rheumatoid Arthritis inception cohort study with 30 healthy controls and identified 292 potential rheumatoid arthritis-specific peptides. Amongst them, 39 were used to create a classifier model using support vector machine algorithms. Specific peptidic fragments were differentially excreted between groups; fragments of protein S100-A9 and gelsolin were less abundant in rheumatoid arthritis while fragments of uromodulin, complement C3 and fibrinogen were all increasingly excreted. The model generated was subsequently tested in an independent test-set of 31 samples. The classifier demonstrated a sensitivity of 88% and a specificity of 93% in diagnosing the condition, with an area under the receiver operating characteristic curve of 0.93 (p<0.0001). These preliminary results suggest that urinary biomarkers could be useful in the early diagnosis of rheumatoid arthritis. Further studies are currently being undertaken in larger cohorts of patients with rheumatoid arthritis and other athridities to assess the potential of the urinary peptide based classifier in the early detection of rheumatoid arthritis.

Introduction

Rheumatoid arthritis (RA) is a systemic autoimmune condition that primarily affects the joints and can lead to joint damage, disability and premature mortality. Early diagnosis and treatment of RA are associated with better patient response to treatment [1], [2], reduced co-morbidity and lower mortality [3]. In the initial stages of the disease, accurate diagnosis can be challenging. In 2010, the American College of Rheumatology and European League Against Rheumatism (ACR/EULAR) developed a new approach to classifying RA based on scoring criteria [4]. This classification system improves sensitivity for the early detection of the disease compared to the former ACR 1987 classification criteria. However it has been shown to lead to both over- and under-diagnosis [5]. The biomarkers rheumatoid factor (RF) and anti-citrullinated protein antibodies (ACPA) form one of the current classification criteria, but in isolation the markers lack sensitivity [6]. Novel biomarkers which could assist in accurate, early diagnosis would facilitate more effective early intervention whilst limiting exposure to disease modifying therapy in patients otherwise destined to remit spontaneously. Recent studies have sought novel RA biomarkers in peripheral blood and synovial fluid [7][10]. Despite early promises, none of these approaches have yet yielded combinations of biomarkers with better specificity and sensitivity than ACPA used alone.

A novel approach for identification of diagnostic biomarkers in RA uses capillary electrophoresis coupled to mass spectrometry (CE-MS). This method has already enabled the identification of urinary biomarker classifiers for the diagnosis of diseases like chronic kidney disease [11], acute kidney injury [12], stroke [13], and cardiovascular diseases [14]. It allows classification of case versus control groups with good accuracy [15]. The use of urine rather than blood for the identification of biomarkers has several advantages, including non-invasive sample collection, a high stability due to absence of proteolytic agents and a low dynamic range of analyte concentration which facilitates the detection and quantification of peptides [16]. Furthermore, discovery of proteomic biomarkers may be useful in understanding the molecular mechanisms involved in the onset and progression of disease [17]. In this study, we aimed to identify potential biomarkers for the early diagnosis of RA. We hypothesized that RA-specific peptides would be measured in urine samples of patients and produce a unique fingerprint of peptides compared to healthy controls. This was achieved by comparing the urinary peptide profile of patients with RA with that of healthy controls. Our secondary objective was to identify the proteases involved in generating the RA-specific peptide fragments, using an in silico approach.

Materials and Methods

RA population and data collection

Patients are recruited to the Scottish Early RA (SERA) inception cohort if they have newly diagnosed undifferentiated arthritis or rheumatoid arthritis. Detailed demographic and clinical data are recorded; samples of blood, urine and synovial fluid are taken every 6 months, and stored for future analyses.

Patients enrolled in SERA who fulfilled the ACR/EULAR diagnostic criteria for RA at diagnosis were identified, and 25 ACPA positive (>20 units) and 24 ACPA negative (<7 units) patients were selected at random. Clinical information on gender, age, disease duration from onset of symptoms, 28 joint count disease activity score (DAS28) [18], health assessment questionnaire (HAQ), C-reactive protein (CRP), ACPA and RF values were collected. Laboratory tests had been analysed in routine National Health Service (NHS) laboratories across Scotland.

Control population

Controls were volunteers who agreed to participate in nutritional-based interventions, were at least 18 years of age, non-smokers and in general good health. They were enrolled in other on-going proteomic studies and had given written informed consent prior to starting the studies.

Ethics statement

Cases provided enduring and generic written consent for their samples to be used in analyses, and the SERA study was approved by the West of Scotland Research Ethics Committee. Controls gave written informed consent prior to participating in on-going proteomic studies approved by the University of Glasgow Faculty of Medicine Ethics Committee.

Urine collection

For both cases and controls, spot urine samples were collected in sterile containers at the time of the study visit. One millilitre aliquots of unprocessed samples were stored at −80°C prior to being processed as recommended by the European Kidney and Urine Proteomics and Human Kidney and Urine Proteome Project and described previously [16]. Most samples were frozen within 4 hours of collection but samples from remote parts of Scotland were couriered in chilled containers overnight. Once frozen, samples were left to defrost at room temperature, occasionally mixing by gentle inversion, and samples were defrosted only once.

Urinary proteomic analysis

The urine samples were prepared as previously described by removing large proteins (>20 kDa), urea, electrolytes and salts, and by enriching polypeptides [11]. Briefly, 700 µL of urine were defrosted with the addition of 0.1% PMSF saturated in ethanol and diluted with 700 µL of a solution containing 2 M urea, 0.1 M NaCl, 10 mM NH4OH and 0.02% SDS. The mixture was then filtered through a 20 kDa MW cut-off ultra-centrifugation filter device (Sartorius Stedim UK Ltd, United Kingdom) at 2,600 × g for one hour at 4°C. A volume of 1.1 mL of the filtrate was then loaded onto a pre-equilibrated PD-10 desalting column (GE Healthcare, Sweden) and eluted using 0.01% aqueous NH4OH. The eluate was subsequently freeze-dried and stored at 4°C prior to being resuspended in HPLC-grade water to a final protein concentration of 2 mg/mL for capillary electrophoresis-mass spectrometry analysis.

Capillary Electrophoresis-Mass Spectrometry (CE-MS) analysis and Data processing

CE-MS analysis was performed as previously described using a P/ACE MDQ capillary electrophoresis system (Beckman Coulter, Fullerton, USA) on line coupled to a MicroTOF MS (BrukerDaltonic, Bremen, Germany) [19]. Samples were injected hydrodynamically at 2.0 psi for 99 sec (ca. 250 nL) and separation of peptides was achieved by reverse polarity at 25 kV for the first 30 min, and with increasing pressure (up to 0.5 psi) for another 34 min. The cartridge temperature was maintained at 25°C. Running buffer contained 79∶20∶1 (v/v) deionised filtered (0.2 µm) water, acetonitrile and formic acid. Sheath liquid consisted of 30% 2-propanol and 0.4% formic acid in deionised filtered (0.2 µm) water. The ESI sprayer (Agilent Technologies, Palo Alto, CA, USA) was grounded, and the ion spray interface potential was set between −4 and −4.5 kV. Spectra were accumulated every 3 seconds over a range of mass-to-charge ratios from 50 to 3000. Details on accuracy, precision, selectivity, sensitivity, reproducibility, and stability of the CE-MS method have been previously described [11], [20]. MosaiquesVisu was used to analyse the CE-MS data [21]. Peptides are initially characterised by their molecular mass, CE-migration time, and ion signal intensity (amplitude) value. Internal standard peptides were used for calibration, as previously described [22]. All detected peptides were deposited, matched, and annotated in a MicrosoftSQL database, allowing for further analysis and comparison between case and control groups.

Statistical analysis and classifier development

Of the initial cohort of RA patients (n = 49) and healthy controls (n = 45), 33 cases and 30 controls were selected at random to establish a panel of RA-specific urinary peptides. After testing for normal distribution, continuous data were compared by the Mann-Whitney test, as this test has proven to be of superior statistical power in proteomics datasets [23]. This test is particularly suited for proteomics data as such data suffer from missing values leading to non-normal skewed distributions even after log-transformation. A p-value of <0.05 was considered to be statistically significant, after correction for multiple testing. Only peptides with a frequency of at least 40% in either group were considered for further analysis. In order to control for the false discovery rate at 0.05, the p-values were adjusted by the Benjamini and Hochberg method [24] implemented in the Bioconductor package multtest [25].

Correlative association between significantly excreted peptides and potential confounding factors such as age and gender was assessed in both case and control groups using Spearman’s rank coefficient correlation (age) and Mann-Whitney test (gender) followed by p-value adjustment using the Benjamini and Hochberg method as previously described. Due to the small sample size of the study, we further analyzed the correlation of the potential biomarkers with age and gender in a cohort of 500 healthy patients from our database.

Significantly differentially excreted peptides were subsequently identified using MS-MS and used to develop a biomarker model to classify between RA and non-RA patients. The remaining 31 samples (16 cases and 15 controls) were used as an independent validation subset in the support-vector-machine (SVM) based MosaCluster software [14], [19]. The sensitivity, specificity and area under the receiver operating characteristic (ROC) curve of the resulting model were calculated using MedCalc version 12.1.3.0 (MedCalc Software bvba, Belgium).

The SVM classifier uses the log transformed intensities of n peptides as coordinates in an n dimensional space. It then builds an n-1 hyper plane that spans this space by performing a quadratic programming optimization of a Lagrangian using the training labels only while allowing for samples to lie on the wrong side of the plane. For such misclassification, the SVM introduces a cost parameter C. Due to the fact that non separable problems in low dimensions may be separable in higher dimensions, the SVM uses the so-called Kernel-trick to transform the samples to a higher dimensional space. Mosacluster uses the standard radial basis functions as kernel. These functions are Gaussians with the parameter gamma controlling for the width. The optimal parameters C and gamma are found via a leave-one-out cross validation error estimation. Implantation of SMV is popular in data mining software, and the Kernel-based Machine Learning Lab (kernlab) package in R in particular is used as a versatile tool for building SVM-based classifiers [26].

Proteases prediction

In order to link urinary fragments to the proteases involved in their generation, a predictive analysis was carried out using Proteasix. This is an open-source tool used to predict the proteases involved in naturally occurring peptide generation in silico, as previously described [27]. Briefly, Proteasix is a cleavage site database that can associate proteases with their corresponding cleavage site sequences based on octopeptides (P4P3P2P1-P1′P2′P3′P4′). Each peptide is described by the respective substrate SWISS-PROT identifier or name, peptide start and peptide end. The search predicts protease association with 0 and up to 3 mismatches in the cleavage site sequence. After entering the peptide list, the tool aligns each peptide sequence with the full-length SWISS-PROT sequence to identify N- and C-term cleavage sites. Each cleavage site is searched in the database to retrieve all predicted protease/cleavage site combinations. This required the generation of a discovery matrix of unique peptide versus protease, and every matrix point is either 0 (no cut), 1 (cut) or 2 (cuts both N-and C-terminus). The matrix was subsequently divided between peptides that were up- or down-regulated based on the absolute fold-change values, and the occurrences were computed by summing the matrix points. Data were assessed using a mathematical approach of frequency analysis. Frequency analysis uses the calculated frequencies of protease association with peptides per protease in either up- or down-regulated groups. The analysis was done using the normal/equal distribution as a reference. Frequency scores were calculated, per protease, using the absolute ratios of the difference over the sum of the frequencies of peptide occurrences that were found up- and down- regulated, multiplied by the difference of occurrences in the up- and down-regulated groups. This mathematical model provides an indication of the distribution of peptide occurrences from the normal distribution, where clustering of peptides around the normal distribution line being interpreted as a lack of specificity in the protease activity.

Results

Descriptive data

Median and interquartile ranges are shown in Table 1 for age, duration of symptoms, ACPA, RF, DAS28, HAQ and CRP values for the training and test sets. The patients all fulfilled the 2010 ACR/EULAR diagnostic criteria for RA [4], and had a median duration of symptoms of 113 days in the training set and 91 days in the test set. Most patients had moderate to severe disease activity in both groups, with a DAS28 of greater than 3.2 [28]. In the training set, 48% were positive for ACPA compared to 56% in the test set. A preliminary analysis investigating the correlation between DAS28, HAQ and CRP revealed a poor correlation between CRP and DAS28, and CRP and HAQ (coefficient of determination between DAS28 and HAQ of 0.310 with p-value = 0.001; coefficient of determination between CRP and DAS28 of 0.053 with p-value = 0.124; coefficient of determination between CRP and HAQ of 0.126 with p-value = 0.015).

thumbnail
Table 1. Baseline characteristics of the training and test set populations (case and control)1.

https://doi.org/10.1371/journal.pone.0104625.t001

Urinary biomarkers associated with RA

The work flow used for establishing urinary biomarkers associated with RA is shown in Figure 1. In order to establish potential urinary biomarkers associated with RA, urines samples from 33 RA patients and 30 healthy volunteers were run using CE-MS and analysed for their peptidomic profile.

thumbnail
Figure 1. Work flow used for the determination of urinary biomarkers associated with RA.

https://doi.org/10.1371/journal.pone.0104625.g001

A total of 292 peptides were significantly differentially regulated in the RA patients compared to controls, after adjustment for multiple testing using the Benjamini-Hochberg procedure for false discovery rate. Of the 292 potential biomarkers, 70 peptides were identified by MS/MS sequencing.

Of these, 34 were found in lower (Table 2) and 36 in higher concentrations (Table 3) in RA patients. Peptide fragments found in significantly lower abundance were identified as protein S100-A9, gelsolin, alpha-1-microglobulin, C-X-C motif chemokine 16, mucin-1 subunit alpha, carboxypeptidase A1 and T calcium channel alpha 1G subunit variant 249. Peptide fragments with higher abundance in RA urine were myosin light chain 3, uromodulin, vesicular integral-membrane protein VIP36, complement 3, fibrinogen alpha chain, clusterin and CD99 antigen. There were also proteins from which peptide fragments were found in both higher and lower concentrations in patients with RA and these were fragments of collagen 1A1 and 1A3, sodium/potassium transporting-ATPase subunit gamma and haemoglobin subunits.

thumbnail
Table 2. Urinary peptides which were significantly less abundant in patients with RA.

https://doi.org/10.1371/journal.pone.0104625.t002

thumbnail
Table 3. Urinary peptides which were found in significantly higher concentration in patients with RA.

https://doi.org/10.1371/journal.pone.0104625.t003

Model of RA-specific classifier biomarkers

As the RA population and the controls showed systematic differences in age and gender, their confounding aspect was checked by correlating each of the 70 biomarkers with age and gender in a cohort of 500 healthy patients from our database. This resulted in 31 biomarkers showing a correlation with age, gender or both (see Table S1). These biomarkers were therefore excluded from the RA classifier biomarker model. The biomarker model was established using the remaining 39 sequenced peptides that were significantly different between controls and cases. Accuracy of the model in the training set was 100% when tested employing complete take-one-out crossvalidation.

When tested on the blinded test set of 31 samples (16 cases and 15 controls), the proteomic signal was significantly different between groups, and the polypeptidic profiles obtained are shown in Figure 2. Accuracy of the model in the independent test set was 91%, and an area under the curve (AUC) of 0.93 on ROC analysis. Median value of the 39 biomarker classifier model was 0.955 (−2.563; 2.394) with sensitivity of 88% and specificity of 93% for identification of RA (p<0.0001). Median values (range) of the classifier model obtained for the case and control groups were 1.542 (0.175; 2.394) and −0.087 (−2.563; 1.809) respectively.

thumbnail
Figure 2. Urinary polypeptide signatures in cases and controls from the validation set based on 39 significantly different sequenced peptides.

Normalized molecular weight (500–15000 Da) in logarithmic scale is plotted against normalized migration time (18–45 minutes). The mean signal intensity of the polypeptide peak is given in 3-dimensional depiction (n = 15 controls and 16 cases).

https://doi.org/10.1371/journal.pone.0104625.g002

Prediction of changes in protease activity based on RA-associated urinary biomarkers

Predicted changes in protease activity, based on the association between a protease and its corresponding cleavage site sequence in the octapeptide form, was carried out on peptides which were significantly differentially regulated on the entire cohort of cases (n = 49) and controls (n = 45). Proteasix search and subsequent analysis resulted in the identification of 131 peptides associated with 32 proteases, resulting in 1271 protease-peptide pairs. The frequency distribution analysis presented here (computed as frequency scores) provides a mathematical modelling of the data taking into account the frequency of peptide occurrences per protease weighed by the difference of occurrences. This analysis, as opposed to a statistical approach based on chi-squared frequency distribution or test of independency, takes into account not only the frequency of distribution but also the difference of occurrences, therefore discarding proteases resulting in a similar number of peptides found equally up- and down- regulated. As plotted in Figure 3, proteases found the furthest away from the normal distribution line (demonstrating an increased specificity in cleavage activity) and with a greater difference in the number of occurrences between the up- and down- regulated groups (as indicated by the higher frequency score in Table 4, with negative scores emphasising a down-regulation) are potential candidates involved in the pathophysiology of RA. Amongst the 32 proteases associated with urinary fragments of RA-related peptides, the potential predicted activities of kallikrein 6 (KLK6) and plasminogen (PLG) were up-regulated, whereas the associated activities of matrix metalloproteinase 3 (MMP3), MMP8, MMP9 and MMP13 and cathepsin B were down-regulated in cases compared to controls (Table 4).

thumbnail
Figure 3. Graphical representation of the frequency distribution of proteases with modified activity associated with RA.

Percentage frequency of peptide occurrences in the down-regulation group is plotted on the x-axis, whereas the percentage frequency of occurrences in the up-regulated group is plotted on the y-axis. Circled data points represent the proteases which activity is the most affected in RA compared to that of healthy controls (see Table 4).

https://doi.org/10.1371/journal.pone.0104625.g003

thumbnail
Table 4. Predictive analysis of changes in protease activity associated with peptides differentially regulated in RA 1.

https://doi.org/10.1371/journal.pone.0104625.t004

Discussion

Based on the comparison of 33 newly diagnosed patients with RA and 30 healthy controls, we were able to identify 292 potential urinary biomarkers associated with the diagnosis of RA. Of these, 70 were sequenced and identified by MS/MS, and 39 used to develop a biomarker model for RA after adjusting for age and gender. This is the first study to show that a urinary biomarker model has potential to assist with the early diagnosis of RA, with good sensitivity and specificity (88% and 93% respectively). Although repeat analysis with better matched controls is required to evaluate if the model is truly specific for RA, these initial results compare favourably to ACPA [6]. Furthermore, a panel of biomarkers like this usually performs better as a diagnostic tool than single or few biomarkers alone [29].

Urinary peptides and protein fragments are the end products of upstream proteolytic processes so the differential urinary excretion of peptides between controls and RA patients may indicate their role in the pathophysiology of the disease. Some of the protein fragments identified originated from proteins known to be up/down regulated in RA, including collagens [30], [31], gelsolin [32], and fibrinogen alpha [7], [33]. To identify which proteases might be responsible for the urinary biomarkers identified, we performed a Proteasix prediction analysis. This suggested an increased activity of KLK-6 and PLG and a reduced activity of cathepsin B, MMP3, MMP8, MMP9 and MMP13.

The kallikrein-kinin system with activation plasma (KLKB1) as well as tissue kallikrein (KLK1) has been implicated in inflammation, and raised levels of both have been found in plasma and synovial fluid of patients with RA [34][38]. Less is known about the role of KLK6 in the pathogenesis of RA, but KLK6 has been found in synovial fluid in patients with psoriatic arthritis [39], and KLK6 promotes survival of murine lymphocytes with actions on proteinase-activated receptor 1 [40].

PLG has been implicated in early RA and although the mechanism remains to be fully established, it may cleave components of the complement system and activate protease-activated receptors and MMPs [41]. Both MMPs and cathepsin B contribute to joint destruction in RA [42], [43] and elevated levels of serum MMP1 and MMP3 correlate well with the progression of erosive disease in early disease [44]. It is therefore unexpected that the protease prediction analysis consistently suggested down-regulation of MMPs and cathepsin B in our RA cohort. Future studies are required to evaluate and explain these findings which may reflect technical elements of our approach or allude to as yet unexplained renal biology in RA patients.

Further validation is also needed to confirm the diagnostic value of urinary biomarkers in early RA. In particular, comparison should be made between the urinary proteome of patients with RA and aged and gender matched patients with other chronic inflammatory conditions and arthritidies. This will establish if the peptides and peptidases identified are specific to RA or reflect chronic inflammation and joint degradation in general. It would also be informative to take into account disease characteristics such as the 2010 classification score and the presence of erosions. This was not possible in the current study due to the small sample size, which also has implications for the power of the study as a whole. It will therefore be necessary to use a bigger and better characterised cohort when validating our findings. This will allow for estimations of sensitivity and specificity of our biomarker profile for diagnostic purposes and will also allow more powerful correlations between individual peptides and disease markers and characteristics. The current study is based on a cross-sectional design but a longitudinal design would be required to monitor disease progression, to investigate if the urinary proteome at baseline is predictive of outcome. A longitudinal design would also allow monitoring changes in the urinary proteome over time in the same patients with correlation to clinical assessments and response to treatment. In cardiovascular disease, for example, it has been shown that the urinary biomarker pattern becomes healthier in patients who had received treatment [14].

Supporting Information

Table S1.

31 biomarkers showed a correlation with age, gender or both.

https://doi.org/10.1371/journal.pone.0104625.s001

(XLS)

Author Contributions

Conceived and designed the experiments: DP HM IM. Performed the experiments: AS. Analyzed the data: AS HH HJ MD. Contributed reagents/materials/analysis tools: JK WM. Wrote the paper: AS DP HJ. Obtained SERA samples: HJ.

References

  1. 1. Anderson JJ, Wells G, Verhoeven AC, Felson DT (2000) Factors predicting response to treatment in rheumatoid arthritis: the importance of disease duration. Arthritis Rheum 43: 22–29.
  2. 2. O’Dell JR (2002) Treating rheumatoid arthritis early: a window of opportunity? Arthritis Rheum 46: 283–285.
  3. 3. Symmons DP, Jones MA, Scott DL, Prior P (1998) Longterm mortality outcome in patients with rheumatoid arthritis: early presenters continue to do well. The Journal of Rheumatology 25: 1072–1077.
  4. 4. Aletaha D, Neogi T, Silman AJ, Funovits J, Felson DT, et al. (2010) 2010 Rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Annals of the Rheumatic Diseases 69: 1580–1588.
  5. 5. Cader MZ, Filer A, Hazlehurst J, de Pablo P, Buckley CD, et al. (2011) Performance of the 2010 ACR/EULAR criteria for rheumatoid arthritis: comparison with 1987 ACR criteria in a very early synovitis cohort. Ann Rheum Dis 70: 949–955.
  6. 6. Whiting PF, Smidt N, Sterne JA, Harbord R, Burton A, et al. (2010) Systematic review: accuracy of anti-citrullinated Peptide antibodies for diagnosing rheumatoid arthritis. Ann Intern Med 152: W155–W166.
  7. 7. Chandra PE, Sokolove J, Hipp BG, Lindstrom TM, Elder JT, et al. (2011) Novel multiplex technology for diagnostic characterization of rheumatoid arthritis. Arthritis Res Ther 13: R102.
  8. 8. De Franceschi L, Bosello S, Scambi C, Biasi D, De Santis M, et al. (2011) Proteome analysis of biological fluids from autoimmune-rheumatological disorders. Proteomics Clin Appl 5: 78–89.
  9. 9. Kriegsmann M, Seeley EH, Schwarting A, Kriegsmann J, Otto M, et al. (2012) MALDI MS imaging as a powerful tool for investigating synovial tissue. Scand J Rheumatol 41: 305–309.
  10. 10. Wild N, Karl J, Grunert VP, Schmitt RI, Garczarek U, et al. (2008) Diagnosis of rheumatoid arthritis: multivariate analysis of biomarkers. Biomarkers 13: 88–105.
  11. 11. Good DM, Zurbig P, Argiles A, Bauer HW, Behrens G, et al. (2010) Naturally occurring human urinary peptides for use in diagnosis of chronic kidney disease. Mol Cell Proteomics 9: 2424–2437.
  12. 12. Metzger J, Kirsch T, Schiffer E, Ulger P, Mentes E, et al. (2010) Urinary excretion of twenty peptides forms an early and accurate diagnostic pattern of acute kidney injury. Kidney Int 78: 1252–1262.
  13. 13. Dawson J, Walters M, Delles C, Mischak H, Mullen W (2012) Urinary proteomics to support diagnosis of stroke. PLoS One 7: e35879.
  14. 14. Delles C, Schiffer E, von Zur Muhlen C, Peter K, Rossing P, et al. (2010) Urinary proteomic diagnosis of coronary artery disease: identification and clinical validation in 623 individuals. J Hypertens 28: 2316–2322.
  15. 15. Mischak H, Schanstra JP (2011) CE-MS in biomarker discovery, validation, and clinical application. Proteomics – Clinical Applications 5: 9–23.
  16. 16. Mischak H, Kolch W, Aivaliotis M, Bouyssié D, Court M, et al. (2010) Comprehensive human urine standards for comparability and standardization in clinical proteome analysis. Proteomics – Clinical Applications 4: 464–478.
  17. 17. Mischak H, Rossing P (2010) Proteomic biomarkers in diabetic nephropathy–reality or future promise? Nephrology Dialysis Transplantation 25: 2843–2845.
  18. 18. Prevoo ML, van ‘t Hof MA, Kuper HH, van Leeuwen MA, van de Putte LB, et al. (1995) Modified disease activity scores that include twenty-eight-joint counts. Development and validation in a prospective longitudinal study of patients with rheumatoid arthritis. Arthritis Rheum 38: 44–48.
  19. 19. Carty DM, Siwy J, Brennand JE, Zurbig P, Mullen W, et al. (2011) Urinary proteomics for prediction of preeclampsia. Hypertension 57: 561–569.
  20. 20. Mischak H, Vlahou A, Ioannidis JP (2013) Technical aspects and inter-laboratory variability in native peptide profiling: the CE-MS experience. Clin Biochem 46: 432–443.
  21. 21. Wittke S, Fliser D, Haubitz M, Bartel S, Krebs R, et al. (2003) Determination of peptides and proteins in human urine with capillary electrophoresis–mass spectrometry, a suitable tool for the establishment of new diagnostic markers. Journal of Chromatography A 1013: 173–181.
  22. 22. Jantos-Siwy J, Schiffer E, Brand K, Schumann G, Rossing K, et al. (2009) Quantitative urinary proteome analysis for biomarker evaluation in chronic kidney disease. J Proteome Res 8: 268–281.
  23. 23. Dakna M, Harris K, Kalousis A, Carpentier S, Kolch W et al.. (2010) Addressing the challenge of defining valid proteomic biomarkers and classifiers. BMC Bioinformatics 11.
  24. 24. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B (Methodological) 57: 289–300.
  25. 25. Gentleman R, Carey V, Bates D, Bolstad B, Dettling M, et al. (2004) Bioconductor: open software development for computational biology and bioinformatics. Genome Biology 5: R80 10.1186/gb-2004-5-10-r80.
  26. 26. Karatzoglou A, Smola A, Hornik K, Zeileis A (2004) kernlab – An S4 Package for Kernel Methods in R. Journal of Statistical Software. 11: 1–20.
  27. 27. Klein J, Eales J, Zürbig P, Vlahou A, Mischak H, et al. (2013) Proteasix: A tool for automated and large-scale prediction of proteases involved in naturally occurring peptide generation. Proteomics 13: 1077–1082.
  28. 28. Fransen J, van Riel PLCM (2005) The Disease Activity Score and the EULAR response criteria. Clin Exp Rheumatol 23: S93–S99.
  29. 29. Schiffer E, Bick C, Grizelj B, Pietzker S, SchÂfer W (2012) Urinary proteome analysis for prostate cancer diagnosis: Cost-effective application in routine clinical practice in Germany. International Journal of Urology 19: 118–125.
  30. 30. Mullan RH, Matthews C, Bresnihan B, FitzGerald O, King L, et al. (2007) Early changes in serum type ii collagen biomarkers predict radiographic progression at one year in inflammatory arthritis patients after biologic therapy. Arthritis & Rheumatism 56: 2919–2928.
  31. 31. Verstappen SMM, Poole AR, Ionescu M, King LE, Abrahamowicz M, et al. (2006) Radiographic joint damage in rheumatoid arthritis is associated with differences in cartilage turnover and can be predicted by serum biomarkers: an evaluation from 1 to 4 years after diagnosis. Arthritis Research & Therapy 8: R31 10.1186/ar1882.
  32. 32. Osborn T, Verdrengh M, Stossel T, Tarkowski A, Bokarewa M (2008) Decreased levels of the gelsolin plasma isoform in patients with rheumatoid arthritis. Arthritis Research & Therapy 10: R117 10.1186/ar2520.
  33. 33. Zhao X, Okeke N, Sharpe O, Batliwalla F, Lee A, et al. (2008) Circulating immune complexes contain citrullinated fibrinogen in rheumatoid arthritis. Arthritis Research & Therapy 10: R94 10.1186/ar2478.
  34. 34. Isordia-Salas I, Pixley RA, Sainz IM, Martinez-Murillo C, Colman RW (2004) The role of plasma high molecular weight kininogen in experimental intestinal and systemic inflammation. Archives of Medical Research 35: 369–377.
  35. 35. Rahman MM, Bhoola KD, Elson CJ, Lemon M, Dieppe PA (1995) Identification and functional importance of plasma kallikrein in the synovial fluids of patients with rheumatoid, psoriatic, and osteoarthritis. Annals of the Rheumatic Diseases 54: 345–350.
  36. 36. Williams RJ, Henderson LM, Naidoo Y, Cassim B, Elson CJ, et al. (1997) Immunocytochemical analysis of tissue kallikrein and the kinin moiety in rheumatoid synovial fluid neutrophils. Rheumatology 36: 420–425.
  37. 37. Cassim B, Shaw OM, Mazur M, Misso NL, Naran A, et al. (2009) Kallikreins, kininogens and kinin receptors on circulating and synovial fluid neutrophils: role in kinin generation in rheumatoid arthritis. Rheumatology 48: 490–496.
  38. 38. Rahman MM, Lemon MJC, Elson CJ, Dieppe PA, Bhoola KD (1994) Proinflammatory role of tissue kallikrein in modulating pain in inflamed joints. British Journal of Rheumatology 34: 88–90.
  39. 39. Eissa A, Cretu D, Soosaipillai A, Thavaneswaran A, Pellett F, et al. (2013) Serum kallikrein-8 correlates with skin activity, but not psoriatic arthritis, in patients with psoriatic disease. Clin Chem Lab Med 51: 317–325.
  40. 40. Scarisbrick IA, Epstein B, Cloud BA, Yoon H, Wu J, et al. (2011) Functional role of kallikrein 6 in regulating immune cell survival. PLoS One 6: e18376
  41. 41. Judex MO, Mueller BM (2005) Plasminogen activation/plasmin in rheumatoid arthritis: matrix degradation and more. Am J Pathol 166: 645–647.
  42. 42. Vincenti M, Brinckerhoff C (2002) Transcriptional regulation of collagenase (MMP-1, MMP-13) genes in arthritis: integration of complex signaling pathways for the recruitment of gene-specific transcription factors. Arthritis Res 4: 157–164 10.1186/ar401.
  43. 43. Trabandt A, Gay RE, Fassbender HG, Gay S (1991) Cathepsin B in synovial cells at the site of joint destruction in rheumatoid arthritis. Arthritis Rheum 34: 1444–1451.
  44. 44. Green MJ, Gough AKS, Devlin J, Smith J, Astin P, et al. (2003) Serum MMP-3 and MMP-1 and progression of joint damage in early rheumatoid arthritis. Rheumatology 42: 83–88.