- Open Access
Proteomic analysis of synovial fluid from rheumatic arthritis and spondyloarthritis patients
Clinical Proteomics volume 17, Article number: 29 (2020)
The aetiologies and pathogeneses of the joint diseases rheumatoid arthritis (RA) and spondyloarthritis (SpA) are still not fully elucidated. To increase our understanding of the molecular pathogenesis, we analysed the protein composition of synovial fluid (SF) from rheumatoid arthritis (RA) and spondyloarthritis (SpA) patients.
Fifty-six synovial fluid samples (RA, n = 32; SpA, n = 24) were digested with trypsin, and the resulting peptides were separated by liquid chromatography and analysed by tandem mass spectrometry. Additionally, the concentration of cell-free DNA (cfDNA) in the synovial fluid was measured, and plasma C-reactive protein (CRP) was determined.
Three hundred thirty five proteins were identified within the SF. The more abundant proteins seen in RA SF were inflammatory proteins, including proteins originating from neutrophil granulocytes, while SpA SF had less inflammatory proteins and a higher concentration of haptoglobin. The concentration of cell-free DNA in the SF increased together with proteins that may have originated from neutrophils. Plasma CRP levels in both RA and SpA, correlated to other acute phase reactants.
The proteomic results underline that neutrophils are central in the RA pathology but not in SpA, and even though inhibitors of neutrophils (migration, proteinase inhibitors) were present in the SF it was not sufficient to interrupt the disease process.
The rheumatic diseases constitute a group of diseases that affects joints, ligaments, tendons, bones and can also show systemic manifestations. Rheumatoid arthritis (RA) and spondyloarthritis (SpA) are common inflammatory systemic joint diseases, with a prevalence of 0.5–1% and 0.1–0.3% respectively [1, 2]. RA is characterized by autoantibodies, including antibodies to citrullinated proteins (adaptive immune system) and neutrophil infiltration (innate immune system) of the synovial fluid (SF), whereas SpA is an autoinflammatory disorder of the innate immune system . RA is twice as common in women, while SpA is twice as common in males [1, 4]. RA typically affects the small joints of the extremities and as the disease progress, cartilage and bone destruction can occur. In SpA, arthritis often affects the small joints of the spine, sacroiliac joints, and large joints of the extremities. As both diseases progress, cartilage and bone destruction often occur. Extra articular manifestations can include enthesitis, psoriasis, uveitis, and inflammatory bowel disease.
Diagnosis of RA and SpA is based on the clinical manifestations, genetic- and biochemical markers, accompanied by imaging techniques (radiographs and magnetic resonance imaging)[5, 6]. The two main serological tests for the RA diagnosis are rheumatoid factor (RF) , and anti-citrullinated protein antibodies (ACPA) . Citrullination is a deamination of the side chain of the amino acid arginine catalysed by citrullinating enzymes, peptidyl arginine deiminases (PAD), in particular PAD4 and PAD2, of which single nucleotide polymorphisms in PAD4 are associated with RA susceptibility . In RA both RF and ACPA are positively associated with the development of a more severe disease progression .
In SpA, MHC class I type HLA-B27 is present in up to 90% of the patients  and with less than 5% of the patients being RF or ACPA positive . The two diseases can thereby be differentiated, but in rare cases the two diseases can co-exist . In both diseases C-reactive protein (CRP) can be increased during active disease. Additionally, in spite the known differences, understanding of the RA and SpA aetiologies remains incomplete.
The protein composition (Proteome) of SF has previously been investigated. By use of 2D-gel electrophoresis Noh et al.  analysed the proteome of SF from RA patients early and late in the disease development and compared it with serum from healthy donors. They found several low molecular weight proteins in the SF compared with serum. In addition, they showed that tumour necrosis-alpha-induced Adipose-Related Protein and Zinc Finger Protein, ZNF658, could be detected as possible markers for RA in serum. Using LC-MS/MS based proteomics Mateos et al.  pooled SF from 20 RA patients and compared the results to a pool of 20 SF from osteoarthrosis (OA) patients. They found proteins related to inflammation to be dominant in the RA group, and proteins involved in the formation and remodelling of the extracellular matrix in the OA group. Similar results were obtained by Balakrishnan et al.  who pooled 5 RA and 5 OA samples, respectively, and also found increased amounts of inflammatory markers of the S100 protein family and matrix metalloproteinases (MMP) in the RA group.
Using quantitative proteomics, Mahendran et al.  analysed individual SF samples from 10 RA, 10 psoriatic arthritis patients (PsA) and 10 controls. PsA is a subgroup of SpA and, compared to the controls, MMP3 was highly increased in both disease groups. Overall, the SF proteomics of PsA and SpA were highly similar.
To determine the molecular differences in SF from RA and SpA patients, we characterized the proteome of 56 individual SF samples from RA (n=32) and SpA (n=24) patients. This allowed us to avoid pooling samples and thereby reveal minor differences between the disease groups. To compare proteomic results with biochemical parameters, C-reactive protein (CRP) in plasma as acute phase marker, and cell-free DNA (cfDNA) as SF pseudo-marker for neutrophil extracellular traps (NETs) were determined .
Human subjects and biobank samples
The clinical samples were collected under the project “INflamation in ARThritis (INART), approved by The Central Denmark Region Committees on Health Research Ethics (1-10-72-291-12). All patients were >18 years of age, fulfilling the ACR/EULAR criteria for RA and ASAS criteria for SpA, respectively [6, 19]. The diagnosis was supported by determination of RF (38% positive) and ACPA (42% positive) for RA patients, and by determination of HLA-B27 for SpA patients (84% positive) determined by Department of Biochemistry at Aarhus University Hospital, Denmark.
Sodium dodecyl sulphate-polyacrylamide gel electrophoresis of SF
SF was cleared for cellular debris by centrifugation at 600×g for 15 min at 20 °C, before storage at -80 °C. After thawing, the samples were centrifuged at 20,000×g for 10 min at 4 °C. The protein concentration was measured with the Bicinchonic Acid (BCA) Protein Assay (Thermo Fisher Scientific, Waltham, MA) according the manufactures instruction. Five µg protein in SDS-samples buffer (Expedeon, San Diego, CA) were separated by 12% SDS-polyacrylamide gel electrophoresis (Expedeon). As molecular weight standard, 2.5 µl Mark12® (Thermo Fisher Scientific) was used. The proteins were stained with Krypton™ Fluorescent Protein Stain (Thermo Fisher Scientific) according the manufactures instruction and scanned on an Amersham Typhoon Biomolecular Imager (GE Healthcare, Chicago, IL).
Sample preparation for proteomics
For sample preparation “filter-aided sample preparation” (FASP) was used [20, 21]. Briefly, 100 µg SF-protein was dissolved in 5% (w/v) sodium deoxycholate (SDC) in 50 mM triethylammonium bicarbonate (TEAB). The samples were heated to 90 °C for 5 min. Molecular weight cut-off Spinfilters 10 kDa (YM10; Millipore, Sigma-Aldrich, St. Louis, MO, USA) were used for buffer exchange between the different steps. The samples were reduced with 12 mM Tris(2-carboxyethyl) phosphine hydrochloride (TCEP), alkylated with 40 mM iodoacetamide (IAA) and digested with 0.4 μg sequencing grade modified trypsin (Promega, Fitchburg, Wisconsin, USA) resuspended in 0.5% SDC, 50 mM TEAB. After digestion, the peptides were collected, and acidified with 0.1% trifluoroacetic acid (TFA). The peptide product was purified using ethyl acetate extraction and the final product was dried down in a vacuum centrifuge and stored at − 80 °C. Prior to analysis, the samples were resuspended in 2% acetonitrile (ACN) and 0.1% TFA.
Mass spectrometry-based proteomics analysis
The mass spectrometry-based analysis was performed according to Bennike et al.  in a randomized patient order. The protein solution was analysed on an automated LC-electrospray ionization (ESI) MS/MS system using an Ultimate 3000 UPLC system with a nanopump module. The system was coupled online to a Thermo-Electron Q Exactive Plus mass spectrometer (Thermo Scientific, Waltham, USA) with an emitter for nanospray ionization. Triplicate runs of each sample (5% of digested material) were loaded onto the C18 reversed phase column (Dionex; Acclaim PepMap100 C18, 5 μm precolumn and 50 cm Acclaim Pepmap RSLC, 75 μm ID main column, Thermo Scientific) and eluted with a linear gradient of 96% solvent A (1% formic acid) and 4% solvent B (acetonitrile) which was increased to 35% solvent B on a 90 min ramp gradient. The MS was operated in data dependent acquisition (DDA) mode, selecting the 12 precursor-ions with the highest intensity for higher energy collisional dissociation (HCD) fragmentation. The raw- and processed data have been made available via ProteomeXchange with identifier PXD010723 .
Proteomics data analysis
A label-free analysis of the proteomics data was performed in MaxQuant v22.214.171.124. The fragment scans were searched against a Uniprot database containing all reviewed Homo sapiens proteins.
(Uniprot reference proteome UPID5640; downloaded 08.2017). The following abundant peptide modifications were included in the analysis: carbamidomethylated cysteine residues (fixed), acetylation of N peptides from the N-terminal of proteins (variable), and oxidation of methionine (variable). The build-in MaxQuant target-decoy search strategy was applied and used to adjust the false discovery rate (FDR) of identified peptides and proteins to max 1%. The MaxQuant MaxLFQ feature, which estimates peptide and protein abundances based on normalized summed peptide precursor intensities, was applied. The resulting label free protein abundance (LFQ) data was processed in Perseus v126.96.36.199 . All protein abundances were log2-transformed. Only for the unsupervised principle component analysis (PCA), did we replace (imputed) missing values with values drawn from a normal distribution to circumvent the problem that PCA cannot handle missing values [24, 25]. This was done by using standard parameters in Perseus for label-free proteomics data, to simulate signals from low-abundant proteins (width = 0.3, downshift = 1.8). Technical replicates were combined by the median, and differentially expressed proteins were identified by t-tests, corrected for multiple hypothesis testing using permutation-based false-positive control with standard parameters in Perseus (s0 = 0.1, FDR < 0.05). Protein function was analysed using Gene Ontology (GO) nomenclature from UniProt protein knowledgebase (UniProtKB) (http://www.uniprot.org) annotation and the software tool “Software tool for researching annotations of proteins” (STRAP) . Unsupervised hierarchical clustering with Euclidean distance calculations was performed on z-score normalized data using standard parameters in Perseus (300 clusters, 10 iterations). Pearson’s correlation analyses of the LFQ values were performed as previously described . Finally, for exploratory analyses we performed linear mixed effect models and random forest modelling using R.
CRP was measured at the Department of Biochemistry at Aarhus University Hospital as part of routine care using a Cobas 6000 (Chemistry XPT).
Cell-free DNA measurement
SF was thawed and centrifuged at 15,000×g for 15 min, diluted 1:25 in 10 mM Tris, pH 8.0 with 1 mM EDTA (TE-buffer). The Quant-iT™ PicoGreen™ dsDNA Assay Kit (ThermoFisher Scientific) was used according to the manufacture’s instruction using 96 well Microplates, PP, F-Bottom black chimney well design (Sigma Aldrich). Fourfold dilution series of DNA were included on all plates (1 µg/ml, 250 ng/ml, 62 ng/ml, 15.6 ng/ml, 3.9 ng/ml, 970 pg/ml, 243 pg/ml, 0 pg/ml). Samples and standards were prepared and measured in duplicates. Plates were measured on an Enspire Multimode Plate Reader (Perkin Elmer, Waltham, MA) with excitation 480 nm and emission 520 nm.
SF was obtained from patients visiting the outpatient clinic at Aarhus University Hospital at the time when therapeutic arthrocentesis was performed. Of the 32 RA patients, five were in treatment with TNF inhibitors (Adalimumab®, Certolizumab® or Etanercept®), four with the IL-6 receptor antagonist Tocilizumab® and three with T-cell activation inhibitor Abatacept®. Nineteen of the RA patients were in methotrexate treatment. Of the 24 SpA patients 6 were in treatment with TNF inhibitors and one with Abatacept®. In the RA group 38% was positive for RF and 42% positive for ACPA. In the group of SpA patients 84% was of the HLA-B27 tissue type.
SDS-PAGE analysis of SF
The protein concentration in the 56 SF-samples varied from 15.7 to 55.4 mg/ml with a mean of 38.1 mg/ml, which is slightly higher than the ~ 25 mg/ml reported in SF from healthy persons. The higher protein-concentration is likely caused by the inflammation . To analyse for major differences in the protein composition, all SF-samples were visualized by 12% SDS-PAGE (Fig. 1a, patient number 12–22). Human albumin bands with a size of 66.5 kDa were seen in all samples with a similar intensity, showing that the adjusted protein load was identical for all samples. Variation was only seen in band patterns below 20 kDa. Comparing the samples, there was no observable correlation between the low molecular weight band patterns and the diagnosis. SF is a filtrate of plasma plus proteins produced locally in the joint. Therefore, we compared a SF sample from a SpA patient to plasma from a healthy participant using SDS-PAGE (Fig. 1b), and only minor differences in low molecular weight bands were seen between normal plasma and SF from SpA patient 19. This indicates that the majority of proteins in an inflamed joint is plasma-derived proteins.
Quantitative proteomics analysis of SF proteins in RA and SpA
Because of the high similarity between the overall protein composition of RA and SpA samples, a label-free proteomics analysis was performed to get a deeper proteome coverage and obtain relative quantitative information to reveal the differences in SF of the 56 samples (RA n = 32, SpA n =24). Each SF sample was digested in solution with trypsin, separated by nanoUPLC and peptides were sequenced using tandem MS by HCD fragmentation. All samples were analysed in triplicates, resulting in 168 MS-runs. Cumulated, we identified 335 proteins (false discovery rate (FDR) < 1%). Following stringent filtering to ensure high-quality quantitative data, 266 proteins were quantifiable cumulated in the SF samples (Additional file 1: Table S1).
Overall SF proteome segregation by gene ontology characterization
The quantifiable proteins were classified according to their functional Gene ontology (GO) information obtained from UniProtKB and GO annotation and visualized with the STRAP software according to “biological process” and “cellular component” (Fig. 2). As shown in Fig. 2a (biological processes), characterization of the biological process (Fig. 2a) identified 85 proteins associated with the immune system, and in addition the groups “regulation” and “cellular process” were dominant. The GO cellular component (Fig. 2b) showed that 252 proteins were classified as “extracellular” in agreement with SF being an extracellular fluid. In addition, many proteins have the annotation “Other intracellular organelles”. This GO nomenclature covers the secretory pathway group. Furthermore, few proteins are seen in the groups: “nucleus”, “cytoplasm” and “cytoskeleton”, supporting a high degree of cellular infiltration in SF from these patients.
Comparison of RA and SpA proteomes by Principal component analysis (PCA)
The overall similarities between the 168 LC–MS runs were investigated by an unsupervised PCA plot on the complete data set before merging the technical triplicates. The scores plot revealed that the technical repeats mostly cluster together (Fig. 2c, encircled). This demonstrates a smaller technical variance than interpersonal difference between the samples, as expected for a sensitive and robust analysis method. The first principal component (the largest possible variation 21%, X-axis) does not discriminate between RA and SpA patients (Fig 2c), whereas the second component (highest variance to first component 11%, Y-axis) is discriminative for some of the RA patients (Fig 2c).
Proteins that contribute the most to the variation in component one and two were analysed by a factor loading plot (Fig. 2d). Principle component 1 was largely separated based on haptoglobin (HP), Fibrinogen beta chain (FGB) and several immunoglobulin proteins encoded from variable (V) gene segments of both heavy chain and light chain. Principle component 2 became separated on basis of inflammatory proteins, including Myeloperoxidase, S100-A12 protein, Lysozyme C, Cathepsin G and Neutrophil gelatinase-associated lipocalin. (Fig. 2d). This indicates that the PC2 separation is based on the inflammatory joint status of the patients. Accordingly, the four RA patients with the highest PC2-value (pt. 22, 25, 34, and 64) had significantly increased SF IgM (mean increase 6.34, p-value 0.04). In these patients CRP was also increased, although not significantly.
Analysis of SF protein differences between RA and SpA
Comparing all identified proteins of RA to all of SpA using t-tests, none of the proteins passed multiple hypothesis correction (q-value < 0.05). However, applying a less strict cut-off without multiple hypothesis correction (p-value < 0.05, log2 (fold change) > 0.5 or < − 0.5), 25 more abundant proteins in RA compared to SpA were identified, in addition to 4 less abundant proteins (Fig. 2e, Table 1). The majority of the more abundant RA-proteins are present in neutrophil granulocytes and monocytes (Table 1 marked a) or proteins involved in the glycolysis (Table 1 marked b). This is in agreement with the preliminary conclusions from the PCA, and indicates the presence of infiltrating neutrophils in the SF from RA patients, in agreement with previous reports . In addition, the data demonstrates the significantly lower/lack of infiltrating neutrophils in SpA, highlighting a central difference in the pathogenesis of RA and SpA.
Correlation of protein changes to concentration of cell-free DNA in SF
Cell-free DNA (cfDNA) in plasma is an inflammatory marker for RA and has been proposed as a predictive marker for biological disease-modifying anti-rheumatic drugs (DMARD) treatment [29, 30]. cfDNA in the SF was measured by a fluorometric method. The SF cfDNA concentration was found to vary between 0.5 to 42.2 µg/ml. These values are an order of tree magnitude higher than can be observed in plasma from RA patients [29, 30], and it reflects that cfDNA originates from cells present in SF. SF cfDNA could be measured in both patient groups, but no statistically significant difference between the intra articular cfDNA concentrations was observed between the RA and SpA patients (RAmean = 7.3 µg/ml, SpAmean = 5.1 µg/ml, t-test p-value = 0.30) (Fig. 2f). The SF samples were centrifuged to remove cells before freezing and to avoid the release of cellular DNA. However, it cannot be excluded that cell lysis occurred after collection and prior to centrifugation. With this reservation in mind, the SF levels of cfDNA was correlated to the SF protein abundance levels to identify functional proteins, which could elucidate the biological interpretation of the cfDNA measurements. Sixty-eight proteins correlated significantly (p-value < 0.05) with the amount of cfDNA in the SF (Table 2).
Correlation of proteins found changed/increased by proteomic in SF to regulatory pathways
To identify underlying biological themes and pathways of the correlating proteins, a GO enrichment analysis using the Reactome pathway database, calculated with all identified SF-proteins as background was performed . The list of synovial fluid proteins with a positive correlation with cfDNA, was significantly enriched (FDR < 0.05) in five biological pathways (Table 3). Thirty-nine of the proteins were tagged as “Immune System” (R-HSA-168256), 31 as “Innate Immune System” (R-HSA-168249) and 26 of the 68 proteins (38%) were categorized as “Neutrophil degranulation” (R-HAS-6798695). The Glycolysis (R-HAS-70171) pathway was also significant enriched, in agreement with the anaerobic metabolism of neutrophils and liberation of the cytoplasmic enzymes upon cell disruption .
NETs are formed in the process of NETosis where neutrophils ejects DNA with histones, in addition to cytoplasmic and secretory granules to form a web-like structure [33, 34]. Of the 11 most correlated proteins to cfDNA, 9 were known NETs proteins (Table 2, labelled “+”). Therefore, the presence of cfDNA correlates with proteins predominantly found in neutrophil granulocytes in both RA and SpA SF, and in vivo liberation of DNA prior to sample collection is therefore likely and not a sampling artefact.
Twenty-nine proteins known to be associated with NETs formation have been identified by MS [34, 35]. Of these 29 proteins, 21 were identified in the SF. To visualize the cfDNA correlation to NETs proteins a hierarchical clustering with the 21 proteins and SF cfDNA was performed (Fig. 3). The analysis showed that high concentration of cfDNA was correlated to presence of NETs proteins. NETs formation is enhanced by Resistin, that can be produced by synoviocytes in joints of RA patients [36, 37]. The receptor for Resistin is the Adenylyl cyclase-associated protein 1 (CAP1) . Accordingly, CAP1 was found to be more abundant in the RA group (Table 1) and correlated significantly and positively with cfDNA (R = 0.5698, p-value = 1.53*10−4) (Table 2). The findings indicate that NETs formation is more abundant in RA-joints compared to SpA, again pointing to the involvement of NETs in RA which recently has been demonstrated to sustain inflammation in other inflammatory diseases.
Correlation of plasma CRP to SF proteins
A common inflammatory marker used in SpA and RA is plasma-CRP, which has a half-life of 19 h . Therefore, in addition to cfDNA, the measured plasma-CRP was correlated to protein abundances in SF, with the aim to identify inflammation-associated proteins. Forty-one proteins were found to correlate significantly with plasma-CRP (p-value < 0.05), 22 of which correlated positively and 19 negatively. Reassuringly, SF-CRP correlated positively with plasma-CRP (R = 0.7148, p-value = 6.12*10−10), verifying the validity of the MS-based proteomics data (Table 4). The other proteins correlating positively to plasma-CRP were also acute phase reactants, e.g. alpha-1-antichymotrypsin, serum amyloid A-1 protein, and, further, other protein groups influenced by plasma-CRP such as complement system proteins (Table 4). Proteins correlating negatively included known “Negative” acute-phase proteins such as albumin, transferrin, transthyretin and retinol-binding protein 4 [39, 40].
Correlation of SF proteins to plasma-CRP and SF cfDNA
The two inflammatory markers, plasma-CRP and SF cfDNA correlated, but not strongly (R = 0.4765, p-value = 2.01*10−4). However, the proteins correlating positively or negatively to each of the two markers were considerably different, and of the 109 statistically significantly correlating proteins, only 12 (11%) were common (Tables 2, 4). Hierarchical clustering with positive correlated proteins to both cfDNA (Fig. 3) and positive as well as negative correlated proteins to plasma CRP, showed a distinct clustering (Fig. 3) with a group of 11, predominantly RA patients with NETs markers (Fig. 3, left). Three of these patients were in biological treatment, showing that the treatment was not responsible for the generation of cfDNA. A group of 10, predominantly SpA patients, correlated with acute phase reactants without NETs proteins (Fig. 3, right), demonstrating the involvement of NETs in RA and the likely subgrouping of the SpA patients. On average a stronger correlating with cfDNA to granulocyte proteins (mean absolute difference 0.13, p-value 1.15*10−6) were seen than to plasma CRP.
The findings presented in the present study is, to our knowledge, the first unbiased proteomic approach of examining and comparing RA and SpA SF protein composition. The findings were further correlated with plasma CRP levels and SF cfDNA for discriminating factors. Among the RA patients, proteins from neutrophils were more dominant in SF. Haptoglobin was the only protein found reduced in SF from RA compared with SpA patients, in accordance with higher levels of the scavenging receptor for haptoglobin–haemoglobin CD163 in RA patients . Proteins from neutrophils and glycolytic enzymes correlated strongly to cfDNA in SF predominantly from RA patients. CRP and other acute phase reactants were seen in both RA and SpA patients, but high amounts of acute phase reactants were also detected in SpA patients without neutrophil granulocyte markers.
NETs is formed by neutrophils undergoing NETosis, where DNA is expelled from the cells within the tissue lining and in biofluids to form web-like structures together with histones, S100-proteins, and proteins stored in secretory granules (Fig. 4) . NETs markers, as circulating cfDNA, can be used to monitor treatment efficiency of biological DMARD in RA patients .
Correlation of SF cfDNA to known NETs proteins in SF was found, indicating that the origin of cfDNA was from neutrophils and thereby that NETosis likely had occurred in the joints in vivo. Due to the presence of both the cytoplasmic S100-proteins and the nuclear histone proteins, degranulation of neutrophils without NETosis is unlikely. The DNA-binding protein, Histone H4, is the most correlated protein with cfDNA. Furthermore, MMP9, NGAL (Neutrophil gelatinase-associated lipocalin) and other NETs proteins were also strongly correlated to the presence of cfDNA in the samples (Table 2), indicating that the cfDNA originated from NETosis process. Cytoplasmic proteins as Glucose-6-phosphate isomerase and Glyceraldehyde-3-phosphate dehydrogenase correlated also to cfDNA indicating cell rupture (Fig. 4).
Neutrophils contain many proteases with tissue destructive effects. In the SF with high cfDNA, MMP9, MMP8, neutrophil elastase, Myeloblastin and Cathepsin G, that can contribute to destruction, were found (Table 2). MMP9 forms a high molecular weight complex with NGAL, that protects MMP9 from degradation and thereby prolongs its proteolytic activity and tissue destruction . NGAL is upregulated in human neutrophils by granulocyte-macrophage colony-stimulating (GM-CSF), and in SF of RA compared to OA patients, higher concentrations of NGAL was measured . This is in agreement with reports of higher plasma MMP9 in RA patients compared with OA, and that a higher enzymatic activity was found in SF of RA patients [44, 45]. Furthermore, neutrophils from RA patients compared to healthy controls have a tendency to undergo NETosis easier when treated with PMA . MMP9 could also be produced by synoviocytes in SF from RA patients, but its high correlation to cfDNA shows that the MMP9 most likely originated from neutrophils .
Metalloproteinase inhibitor 1 and Leukocyte elastase inhibitor both correlated to cfDNA, but the cell type producing these inhibitors are unclear (Fig. 4). Metalloproteinase inhibitor 1 can form a complex with MMP9 and NGAL. Therefore, the in vivo MMP9 activity is difficult to determine. Ahrens et al.  have shown active MMP9 gelatinase activity in SF from RA patients by use of gel zymography, and thus, MMP9 was separated from Metalloproteinase inhibitor 1. Annexin 1, released by neutrophils, is an inhibitor of migration of neutrophils from the blood to the inflamed site, and even though Annexin 1 was present in the SF with cfDNA, it was unable to prevent disease progression (Fig 4).
Katano et al.  proposed GM-CSF as a target for treatment of RA and this is now supported by human trials. Patients with high SF cfDNA and NGAL could therefore be a distinct clinical endotype that may benefit from such a treatment . Sato et al.  proposed that Resistin produced by synovial tissue could be important for the pathogenesis of RA. The receptor for Resistin, CAP1, was present in SF and correlated to cfDNA. CAP1 is a cytoplasmic protein, but is translocated to the cell membrane when Resistin is present . Resistin enhances NETosis , but whether CAP1 in SF can function as a receptor antagonist in its free form is presently unknown.
Plasma CRP is a central marker used in combination with 28-joint Disease Activity Score (DAS28CRP) and Ankylosing Spondylitis Disease Activity Score (ASDAS) to determine disease progression and treatment success, that correlates to disease progression . When plasma CRP was correlated to SF proteomic data, other known positively and negatively regulated acute phase reactants were found (Tables 3 and 4, Fig. 3). Proteins correlated to plasma CRP and cfDNA, respectively, had only 12% proteins in common. Leucine-rich alpha-2-glycoprotein (LRG1) is correlated to both CRP and cfDNA. It is an acute phase protein induced primarily by IL-6 in the liver [50, 51], but is also a neutrophil secondary granule protein released together with lactoferrin, and thus its correlation to both CRP and cfDNA is explainable . LRG1 inhibits the anti-proliferative effect of transforming growth factor ß1 (TGFß1) on myeloid cells . TGFß1 regulates the anti-inflammatory process in the synovial membrane in RA patients, and is important for the self-regulation that can result in remission periods . In RA patients treated with anti-interleukin-6 receptor antibody (RoActemra) LRG1 was a better marker for remission than CRP, matrix metalloproteinase 3 level and erythrocyte sedimentation rate . As shown in the present study, this can be due to LRG1 being a maker for both acute phase reactants and for neutrophil degranulation/NETosis. As shown in Fig. 3, a group of manly SpA patients had high CRP and acute phase reactants without granulocyte markers, confirming the central role of neutrophil granulocytes in the pathogenesis of RA.
The differences between RA and SpA identified here are in line with results from randomized clinical trials with already approved drugs. Thus, T-cell targeted therapies such as inhibitors of IL-17 has shown efficacy in SpA but not in RA. In contrast, therapies targeting myeloid derived cytokines such as inhibitors of IL-1 and IL-6 are effective in RA but not in SpA . Our findings also support that RA is a very heterogenous disease as proposed by others . Thus, a subgroup of RA patients with high neutrophil activation was found. High degree of neutrophil priming and NETosis is also present in other diseases such as systemic lupus erythematosus and granulomatosis with polyangiitis [33, 57]. Some of the drugs used for treatment of these diseases such as the C5a receptor inhibitor Avacopan® could, therefore, also be effective in RA patients with high neutrophil activation. In this way, the study could help guide future drug development to treat immune mediated inflammatory arthritis.
The proteomics of SF from SpA and RA patients showed a marked difference in the amounts of proteins from the innate immune system, primarily originating from neutrophil granulocytes. The presence of these proteins was more pronounced in the RA patient group. These proteins were also correlated to SF cfDNA indicating NETosis. Neutrophils produce IL-6 that induces acute phase reactants, but surprisingly, little correlation between NETs proteins/cfDNA and acute phase reactants/CRP was seen, indicating that two different inflammatory mechanisms are used for increase in CRP and cfDNA. This is in agreement with the recent finding that mature neutrophils are unresponsive to IL-6 due to the absence of gp130 in the IL-6 receptor complex . Some of the patients with high cfDNA were in treatment with anti-interleukin-6 receptor (RoActemra®) or TNF-α antagonists. This may have influenced the amounts of acute phase reactants. However, in a patient material as presented in this study, SF cfDNA is an indicator of intraarticular NETosis. Therefore, SF cfDNA measurement may be used as an indicator of severe arthritis, and our findings demonstrate the involvement of NETs in RA but less in SpA pathogenesis.
Availability of data and materials
The raw- and processed data have been made available via ProteomeXchange with identifier PXD010723
anti-cyclic citrullinated protein antibodies
false discovery rate
higher energy collisional dissociation
label free quantification
major histocompatibility complex
peptidyl arginine deiminases
principle component analysis
tandem mass spectrometry
neutrophil extracellular traps
UniProt protein knowledgebase
ultra-performance liquid chromatography
Linos A, Worthington JW, O’Fallon WM, Kurland LT. The epidemiology of rheumatoid arthritis in Rochester, Minnesota: a study of incidence, prevalence, and mortality. Am J Epidemiol. 1980;111:87–98.
Wang R, Ward MM. Epidemiology of axial spondyloarthritis: an update. Curr Opin Rheumatol. 2018;30:137–43.
Ottonello L, Cutolo M, Frumento G, Arduino N, Bertolotto M, Mancini M, et al. Synovial fluid from patients with rheumatoid arthritis inhibits neutrophil apoptosis: role of adenosine and proinflammatory cytokines. Rheumatol Narnia. 2002;41:1249–60. https://doi.org/10.1093/rheumatology/41.11.1249.
Will R, Edmunds L, Elswood J, Calin A. Is there sexual inequality in ankylosing spondylitis? A study of 498 women and 1202 men. J Rheumatol. 1990;17:1649–52.
Allen A, Carville S, McKenna F. Diagnosis and management of rheumatoid arthritis in adults: summary of updated NICE guidance. BMJ. 2018;362:k3015. https://doi.org/10.1136/bmj.k3015.
Sieper J, Rudwaleit M, Baraliakos X, Brandt J, Braun J, Burgos-Vargas R, et al. The Assessment of SpondyloArthritis international Society (ASAS) handbook: A guide to assess spondyloarthritis. Ann Rheum Dis. 2009;6:8.
Edelman GM, Kunkel HG, Franklin EC. Interaction of the rheumatoid factor with antigen-antibody complexes and aggregated gamma globulin. J Exp Med. 1958;108:105–20.
Bizzaro N, Mazzanti G, Tonutti E, Villalta D, Tozzoli R. Diagnostic accuracy of the anti-citrulline antibody assay for rheumatoid arthritis. Clin Chem. 2001;47:1089–93.
Suzuki A, Yamada R, Chang X, Tokuhiro S, Sawada T, Suzuki M, et al. Functional haplotypes of PADI4, encoding citrullinating enzyme peptidylarginine deiminase 4, are associated with rheumatoid arthritis. Nat Genet. 2003;34:395–402.
van de Stadt LA, van der Horst AR, de Koning MHMT, Bos WH, Wolbink GJ, van de Stadt RJ, et al. The extent of the anti-citrullinated protein antibody repertoire is associated with arthritis development in patients with seropositive arthralgia. Ann Rheum Dis. 2011;70:128–33. https://doi.org/10.1136/ard.2010.132662.
Harper BE, Reveille JD. Spondyloarthritis: clinical suspicion, diagnosis, and sports. 2009. p. 29–34.
Swart A, Burlingame RW, Gürtler I, Mahler M. Third generation anti-citrullinated peptide antibody assay is a sensitive marker in rheumatoid factor negative rheumatoid arthritis. Clin Chim Acta. 2012;414:266–72.
Barczyńska TA, Węgierska M, Żuchowski P, Dura M, Zalewska J, Waszczak M, et al. Coexistence of rheumatoid arthritis and ankylosing spondylitis. Reumatologia. 2015;53:279–85.
Noh R, Park SG, Ju JH, Chi SW, Kim S, Lee CK, et al. Comparative proteomic analyses of synovial fluids and serums from rheumatoid arthritis patients. J Microbiol Biotechnol. 2014;24:119–26.
Mateos J, Lourido L, Fernández-Puente P, Calamia V, Fernández-López C, Oreiro N, et al. Differential protein profiling of synovial fluid from rheumatoid arthritis and osteoarthritis patients using LC-MALDI TOF/TOF. J Proteomics. 2012;75:2869–78. https://doi.org/10.1016/j.jprot.2011.12.042.
Balakrishnan L, Bhattacharjee M, Ahmad S, Nirujogi R, Renuse S, Subbannayya Y, et al. Differential proteomic analysis of synovial fluid from rheumatoid arthritis and osteoarthritis patients. Clin Proteomics. 2014;11:1. https://doi.org/10.1186/1559-0275-11-1.
Mahendran SM, Keystone EC, Krawetz RJ, Liang K, Diamandis EP, Chandran V. Elucidating the endogenous synovial fluid proteome and peptidome of inflammatory arthritis using label-free mass spectrometry. Clin Proteomics BioMed Central. 2019;16:1–13. https://doi.org/10.1186/s12014-019-9243-3.
Sur Chowdhury C, Giaglis S, Walker UA, Buser A, Hahn S, Hasler P. Enhanced neutrophil extracellular trap generation in rheumatoid arthritis: analysis of underlying signal transduction pathways and potential diagnostic utility. Arthritis Res. Ther. 2014;16:R122. https://doi.org/10.1186/ar4579.
Dejaco C, Singh YP, Perel P, Hutchings A, Camellino D, Mackie S, et al. 2015 recommendations for the management of polymyalgia rheumatica: a European League Against Rheumatism/American College of Rheumatology collaborative initiative. Ann Rheum Dis. 2015;74:1799–807.
Bennike TB, Carlsen TG, Ellingsen T, Bonderup OK, Glerup H, Bøgsted M, et al. Neutrophil extracellular traps in ulcerative colitis. Inflamm Bowel Dis. 2015;21:2052–67.
León IR, Schwämmle V, Jensen ON, Sprenger RR. Quantitative assessment of in-solution digestion efficiency identifies optimal protocols for unbiased protein analysis. Mol Cell Proteomics. 2013;12:2992–3005.
Vizcaíno J, Deutsch EEW, Wang R, Vizcaino JA, Deutsch EEW, Wang R, et al. ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat Biotechnol. 2014;32:223–6. https://doi.org/10.1038/nbt.2839%5cn.
Cox J, Hein MY, Luber CA, Paron I, Nagaraj N, Mann M. Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ. Mol Cell Proteomics. 2014;13:2513–26. https://doi.org/10.1074/mcp.m113.031591.
Grung B, Manne R. Missing values in principal component analysis. Chemom Intell Lab Syst. 1998;42:125–39.
Bennike TB, Carlsen TG, Ellingsen T, Bonderup OK, Glerup H, Bøgsted M, et al. Neutrophil extracellular traps in ulcerative colitis. Inflamm Bowel Dis. 2015;1:1.
Bhatia VN, Perlman DH, Costello CE, McComb ME. Software tool for researching annotations of proteins: open-source protein annotation software with data visualization. Anal Chem. 2009;81:9819–23. https://doi.org/10.1021/ac901335x.
Bennike TB, Ellingsen T, Glerup H, Bonderup OK, Carlsen TG, Meyer MK, et al. Proteome analysis of rheumatoid arthritis gut mucosa. J Proteome Res. 2017;16:356–64.
Bennike T, Ayturk U, Haslauer CM, Froehlich JW, Proffen BL, Barnaby O, et al. A normative study of the synovial fluid proteome from healthy porcine knee joints. J Proteome Res. 2014;13:4377–87. https://doi.org/10.1021/pr500587x.
Hashimoto T, Yoshida K, Hashimoto N, Nakai A, Kaneshiro K, Suzuki K, et al. Circulating cell free DNA: a marker to predict the therapeutic response for biological DMARDs in rheumatoid arthritis. Int J Rheum Dis. 2017;20:722–30.
Pérez-Sánchez C, Ruiz-Limón P, Aguirre MA, Jiménez-Gómez Y, Arias-delaRosa I, Ábalos-Aguilera MC, et al. Diagnostic potential of NETosis-derived products for disease activity, atherosclerosis and therapeutic effectiveness in Rheumatoid Arthritis patients. J Autoimmun. 2017;82:31–40.
Haw R, Hermjakob H, D’Eustachio P, Stein L. Reactome pathway analysis to enrich biological discovery in proteomics data sets. Proteomics. 2011;11:3598–613.
Kramer PA, Ravi S, Chacko B, Johnson MS, Darley-Usmar VM. A review of the mitochondrial and glycolytic metabolism in human platelets and leukocytes: implications for their use as bioenergetic biomarkers. Redox Biol. 2014;2:206–10. https://doi.org/10.1016/j.redox.2013.12.026.
Kaplan MJ, Radic M. Neutrophil extracellular traps: double-edged swords of innate immunity. J Immunol. 2012;189:2689–95. https://doi.org/10.4049/jimmunol.1201719.
Urban CF, Ermert D, Schmid M, Abu-Abed U, Goosmann C, Nacken W, et al. Neutrophil extracellular traps contain calprotectin, a cytosolic protein complex involved in host defense against Candida albicans. PLoS Pathog. 2009;5:8.
Carmona-Rivera C, Zhao W, Yalavarthi S, Kaplan MJ. Neutrophil extracellular traps induce endothelial dysfunction in systemic lupus erythematosus through the activation of matrix metalloproteinase-2. Ann Rheum Dis. 2015;74:1417–24.
Jiang S, Park DW, Tadie J-M, Gregoire M, Deshane J, Pittet JF, et al. Human resistin promotes neutrophil proinflammatory activation and neutrophil extracellular trap formation and increases severity of acute lung injury. J. Immunol. 2014;192:4795–803. https://doi.org/10.4049/jimmunol.1302764.
Sato H, Muraoka S, Kusunoki N, Masuoka S, Yamada S, Ogasawara H, et al. Resistin upregulates chemokine production by fibroblast-like synoviocytes from patients with rheumatoid arthritis. Arthritis Res. Ther. 2017;19:6–11.
Vigushin DM, Pepys MB, Hawkins PN. Metabolic and scintigraphic studies of radioiodinated human C-reactive protein in health and disease. J Clin Invest. 1993;91:1351–7.
Larson LM, Namaste SM, Williams AM, Engle-Stone R, Addo OY, Suchdev PS, et al. Adjusting retinol-binding protein concentrations for inflammation: Biomarkers Reflecting Inflammation and Nutritional Determinants of Anemia (BRINDA) project. Am J Clin Nutr. 2017;106:390S–401S.
Ritchie RF, Palomaki GE, Neveux LM, Navolotskaia O, Ledue TB, Craig WY. Reference distributions for alpha2-macroglobulin: a practical, simple and clinically relevant approach in a large cohort. J Clin Lab Anal. 2004;18:139–47.
Greisen SR, Moller HJ, Stengaard-Pedersen K, Hetland ML, Hørslev-Petersen K, Jørgensen A, et al. Soluble macrophage-derived CD163 is a marker of disease activity and progression in early rheumatoid arthritis. Clin Exp Rheumatol. 2011;29:689–92.
Yan L, Borregaard N, Kjeldsen L, Moses MA. The high molecular weight urinary matrix metalloproteinase (MMP) activity is a complex of gelatinase B/MMP-9 and neutrophil gelatinase-associated lipocalin (NGAL): modulation of MMP-9 activity by NGAL. J Biol Chem. 2001;276:37258–65.
Katano M, Okamoto K, Arito M, Kawakami Y, Kurokawa MS, Suematsu N, et al. Implication of granulocyte-macrophage colony-stimulating factor induced neutrophil gelatinase-associated lipocalin in pathogenesis of rheumatoid arthritis revealed by proteome analysis. Arthritis Res Ther. 2009;11:R3.
Sopata I, Wize J, Filipowicz-Sosnowska A, Stanisławska-Biernat E, Brzezińska B, Maślinński S. Neutrophil gelatinase levels in plasma and synovial fluid of patients with rheumatic diseases. Rheumatol Int. 1995;15:9–14.
Ahrens D, Koch AE, Pope RM, Stein-Picarella M, Niedbala MJ. Expression of matrix metalloproteinase 9 (96-kd gelatinase B) in human rheumatoid arthritis. Arthritis Rheum. 1996;39:1576–87. https://doi.org/10.1002/art.1780390919.
Zhou M, Qin S, Chu Y, Wang F, Chen L, Lu Y. Immunolocalization of MMP-2 and MMP-9 in human rheumatoid synovium. Int J Clin Exp Pathol. 2014;7:3048–56.
Burmester GR, McInnes IB, Kremer JM, Miranda P, Vencovský J, Godwood A, et al. Mavrilimumab, a fully human granulocyte-macrophage colony-stimulating factor receptor α monoclonal antibody. Arthritis Rheumatol. 2018;70:679–89. https://doi.org/10.1002/art.40420.
Lee S, Lee HC, Kwon YW, Lee SE, Cho Y, Kim J, et al. Adenylyl cyclase-associated protein 1 is a receptor for human resistin and mediates inflammatory actions of human monocytes. Cell Metab. 2014;19:484–97. https://doi.org/10.1016/j.cmet.2014.01.013.
Wells G, Becker J-C, Teng J, Dougados M, Schiff M, Smolen J, et al. Validation of the 28-joint Disease Activity Score (DAS28) and European League Against Rheumatism response criteria based on C-reactive protein against disease progression in patients with rheumatoid arthritis, and comparison with the DAS28 based on erythr. Ann Rheum Dis. 2009;68:954–60. https://doi.org/10.1136/ard.2007.084459.
Bini L, Magi B, Marzocchi B, Cellesi C, Berti B, Raggiaschi R, et al. Two-dimensional electrophoretic patterns of acute-phase human serum proteins in the course of bacterial and viral diseases. Electrophoresis. 1996;17:612–6.
Shirai R, Hirano F, Ohkura N, Ikeda K, Inoue S. Up-regulation of the expression of leucine-rich α2-glycoprotein in hepatocytes by the mediators of acute-phase response. Biochem Biophys Res Commun. 2009;382:776–9. https://doi.org/10.1016/j.bbrc.2009.03.104.
Druhan LJ, Lance A, Li S, Price AE, Emerson JT, Baxter SA, et al. Leucine rich α-2 glycoprotein: a novel neutrophil granule protein and modulator of myelopoiesis. PLoS ONE. 2017;12:e0170261. https://doi.org/10.1371/journal.pone.0170261.
Taketazu F, Kato M, Gobl A, Ichijo H, ten Dijke P, Itoh J, et al. Enhanced expression of transforming growth factor-beta s and transforming growth factor-beta type II receptor in the synovial tissues of patients with rheumatoid arthritis. Lab Invest. 1994;70:620–30.
Fujimoto M, Serada S, Suzuki K, Nishikawa A, Ogata A, Nanki T, et al. Brief report: leucine-rich α 2 -glycoprotein as a potential biomarker for joint inflammation during anti-interleukin-6 biologic therapy in rheumatoid arthritis. Arthritis Rheumatol. 2015;67:2056–60. https://doi.org/10.1002/art.39164.
Schett G, Elewaut D, McInnes IB, Dayer JM, Neurath MF. Toward a cytokine-based disease taxonomy. Nat Med Nat Med. 2013;19:822–4.
McGonagle D, Watad A, Savic S. Mechanistic immunological based classification of rheumatoid arthritis. Autoimmun. Rev. 2018;17:1115–23. https://doi.org/10.1016/j.autrev.2018.06.001.
Nakazawa D, Masuda S, Tomaru U, Ishizu A. Pathogenesis and therapeutic interventions for ANCA-associated vasculitis. Nat Rev Rheumatol. 2019;15:91–101.
Wilkinson AN, Gartlan KH, Kelly G, Samson LD, Olver SD, Avery J, et al. Granulocytes are unresponsive to IL-6 due to an absence of gp130. J Immunol. 2018;200:3547–55.
We thank laboratory assistant Ditte Bech Laursen for technical assistance with the project.
This study was supported by Aase and Ejner Danielsens Foundation (Grant 10-001785), Danish National Mass Spectrometry Platform for Functional Proteomics (PRO-MS), Danish Rheumatism Association (Grant R116-Rp4652), The Obel Family Foundation (Grant 25508), The Beckett Foundation, The Hertha Christensen Foundation, The Lundbeck Foundation (R181-2014-3372), and The Carlsberg Foundation (CF14-0561).
Ethics approval and consent to participate
The clinical samples were collected under the project “INflamation in ARThritis (INART), approved by The Central Denmark Region Committees on Health Research Ethics (1-10-72-291-12). All patients were > 18 years of age, fulfilling the ACR/EULAR criteria for RA and ASAS criteria for SpA, respectively.
Consent for publication
We agree for publication upon acceptance and we agree that all copyright ownership for the article is transferred to Clinical Proteomics. The material submitted is new, original and has not been submitted to another journal for concurrent consideration.
The authors declared that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Birkelund, S., Bennike, T.B., Kastaniegaard, K. et al. Proteomic analysis of synovial fluid from rheumatic arthritis and spondyloarthritis patients. Clin Proteom 17, 29 (2020). https://doi.org/10.1186/s12014-020-09292-9
- Synovial fluid
- Rheumatoid arthritis
- Cell-free DNA
- Neutrophil extracellular traps