- Open Access
Elucidating the endogenous synovial fluid proteome and peptidome of inflammatory arthritis using label-free mass spectrometry
Clinical Proteomics volume 16, Article number: 23 (2019)
Inflammatory arthritis (IA) is an immunological disorder in which loss of immune tolerance to endogenous self-antigens perpetuates synovitis and eventual destruction of the underlying cartilage and bone. Pathological changes in the joint are expected to be represented by synovial fluid (SF) proteins and peptides. In the present study, a mass spectrometry-based approach was utilized for the identification of key protein and peptide mediators of IA.
Age-matched SF samples from 10 rheumatoid arthritis patients, 10 psoriatic arthritis patients and 10 cadaveric controls were subjected to an integrated proteomic and peptidomic protocol using liquid chromatography tandem mass spectrometry. Significant differentially abundant proteins and peptides were identified between cohorts according to the results of a Mann–Whitney U test coupled to the Benjamini–Hochberg correction for multiple hypothesis testing. Fold change ratios were computed for each protein and peptide according to their log-transformed extracted ion current. Pathway analysis and antimicrobial peptide (AMP) prediction were conducted to clarify the pathophysiological relevance of identified proteins and peptides to IA.
We determined that 144 proteins showed significant differential abundance between the IA and control SF proteomes, of which 11 protein candidates were selected for future follow-up studies. Similar analyses applied to our peptidomic data identified 15 peptide sequences, originating from 4 protein precursors, to have significant differential abundance in IA compared to the control SF peptidome. Pathway enrichment analysis of the IA SF peptidome along with AMP prediction suggests a possible mechanistic role of microbes in eliciting an immune response which drives the development of IA.
The discovery-phase data generated herein has provided a basis for the identification of candidates with the greatest potential to serve as novel serum biomarkers specific to inflammatory arthritides. Moreover, these findings facilitate the understanding of possible disease mechanisms specific to each subtype.
Inflammatory arthritis (IA) is characterized by synovial hyperplasia leading to degradation of adjacent articular cartilage and bone . The term encompasses several forms of inflammatory joint diseases that when taken together, have an annual incidence ranging from 115 to 271 per 100,000 adults . IA is a multifactorial disease driven by the complex interplay of both genetics and the environment. Rheumatoid arthritis (RA), the most common and potentially destructive IA, has a well-established association with class II major histocompatibility complex (MHC) alleles while the spondyloarthritides, such as psoriatic arthritis (PsA), are more frequently associated with class I MHC alleles . Susceptibility to IA increases when genetic predisposition is complemented by environmental risk factors such as smoking, obesity and more recently, microbial infection and intestinal dysbiosis [4,5,6]. The exact etiology of IA is still poorly understood with studies aimed at delineating the molecular pathways driving loss of immunological tolerance to the body’s self-antigens. Alterations to the adaptive and innate immune system perpetuate systemic inflammation and lead to an elevated risk of developing comorbid conditions such as cardiovascular disease, metabolic syndrome, diabetes and depression [7, 8]. Naturally, there is a compelling need to identify markers of aberrant immune pathways relevant to IA which may advance current insights into the molecular mechanisms of the disease and serve as clinical markers for disease monitoring and treatment responses.
The rise in high-throughput technologies, such as next-generation gene sequencing and mass spectrometry (MS), facilitate the discovery of key modulators of disease. Specifically, MS-based approaches provide an essential analytical platform for the identification, quantification and characterization of candidate biomarkers. Biomarkers may come in the form of a molecular signature, a clinical feature or even as an imaging parameter. Molecular biomarkers may be further subtyped into the domains of genomics, transcriptomics, proteomics, metabolomics or peptidomics. Due to the importance of proteins in pathophysiological processes, there is increased interest in resolving the proteomic profile of biospecimens related to IA. Similarly, peptides play a seminal role in mediating physiological functions by serving as neurotransmitters, hormones, antibiotics and immune regulators . During IA, joint pain and inflammation are driven by aberrant proteolysis resulting in the production of inflammatory peptides and the destruction of joint cartilage and bone. Synovial fluid (SF), a proximal fluid which bathes the intrinsic joint structures, is an important reservoir of putative protein and peptide biomarkers whose abundance levels fluctuate in response to pathological changes due to disease .
In the current study, we performed MS-based proteomic and peptidomic analyses of SF from RA and PsA patients to identify and quantify significant proteins and peptides related to the aetiopathogenesis of IA. Differential abundance analyses highlighted the capacity for dysregulated SF proteins and peptides to reflect disease activity while pathway analysis and antimicrobial peptide (AMP) prediction alluded to a larger role of microbes in the initiation and progression of IA. These findings provide the means for discovering novel candidates to serve as possible biomarkers of IA while simultaneously, highlighting possible mechanistic networks responsible for the disease progression of RA and PsA.
Materials and methods
Patients and SF collection
Research ethics board approval was received for the study from the University Health Network, Mount Sinai Hospital and the University of Calgary. Informed consent was obtained from all patients.
SF samples for the study were obtained, retrospectively, from 10 cases with RA, 10 cases with PsA and 10 cadaveric controls. RA patients were classified according to the 1987 American College of Rheumatology (ACR) classification criteria . PsA patients satisfied the Classification Criteria for Psoriatic Arthritis (CASPAR) .
Cadaveric control SF were obtained from joints through the Southern Alberta Tissue Donation Program. Inclusion criteria consisted of an age of 18 years or older with no medical history of arthritis, joint injury or joint surgery (including visual inspection of cartilage surfaces during recovery), no prescription anti-inflammatory medications and availability within 4 h of death. Exclusion criteria for all disease cohorts included patients receiving therapeutic biological drugs and the presence of other causes of inflammation (e.g. infection and/or crystal disease) or co-morbidities (e.g. cancer).
SF sample preparation
IA SF samples (both RA and PsA) were obtained through needle aspiration of knee joints and kept on ice. Samples were transferred to centrifuge tubes and spun at 160 RCF for 10 min at 4 °C. The supernatant was transferred to a sterile 1.5 mL centrifuge tube and spun at 2000 RCF for another 10 min at 4 °C. Samples were immediately stored at − 80 °C until further processing. SF samples from cadavers were collected without the use of lavage. Samples were centrifuged at 3000 RCF for 15 min and stored at − 80 °C until further processing.
At the time of analysis, samples were blinded, randomized, thawed on ice and their respective total protein concentrations were measured with a Pierce Coomassie (Bradford) total protein assay.
SF sample preparation for proteomic analysis
For proteomic investigations, SF samples were first adjusted to 300 µg total protein in 50 mM ammonium bicarbonate (ABC). Protein concentration was conducted using Amicon Ultra-0.5 centrifugal filter units (10 kDa molecular weight cut-off; MilliporeSigma) which were pre-equilibrated with 400 uL of 50 mM ABC. Samples were loaded and spun at 10,000 RPM for 35 min at 4 °C and transferred to a new tube by spinning upside at 5000 RPM for 2 min.
Concentrates were collected and brought to a total volume of 100 µL with 50 mM ABC. Proteins were denatured with powdered urea to a final concentration of 8 M. Dithiothreitol (DTT) (Sigma-Aldrich) was added to each concentrate sample to a final concentration of 5 mM and incubated at 60 °C for 45 min. This was followed by alkylation with 15 mM iodoacetamide (IAM) (Sigma-Aldrich) at room temperature in the dark for 45 min. Samples were diluted fivefold with 50 mM ABC to prevent inhibition of trypsin activity by high concentrations of urea. Concentrate samples were digested with trypsin (Sigma-Aldrich) in a 1:50 (trypsin to total protein) ratio for 20 h at 37 °C and then dropwise acidified to a pH of 2 with formic acid (FA) to inhibit trypsin activity. Samples were reduced to 300 µL via speed vacuum concentration and stored at − 20 °C until subjected to liquid chromatography-tandem mass spectrometry (LC–MS/MS).
SF sample preparation for peptidomics analysis
Peptides were isolated based on a protocol described by Kamphorst et al. . Fifty microliters of SF were diluted in 235 µL of 50 mM ABC and 15 µL dimethyl sulfoxide (DMSO) for peptidomic analysis. Peptide concentration was conducted using Amicon Ultra-0.5 centrifugal filter units (10 kDa MWCO; MilliporeSigma) which were pre-equilibrated with 250 µL of 50 mM ABC. SF samples were spun at 10 000 RPM for 60 min at 4 °C then washed with 100 µL of 50 mM of ABC and spun for another 10 min. Filtrates were acidified with 5 µL of FA.
Peptides were desalted using one hydrophilic-lipophilic-balanced reverse-phase cartridge per sample (Oasis HLB). Each cartridge [1 mL (30 mg); Waters cat# WAT094225] was first pre-equilibrated with 1 mL of 90% acetonitrile (ACN), 0.1% FA and 0.02% trifluoroacetic acid (TFA) and then washed with 3 mL of buffer A (5% ACN, 0.1% FA, 0.02% TFA). The SF sample was then passed through the cartridge and washed a second time with 3 mL of buffer A. Peptides were eluted with 700 µL of 60% ACN, 0.1% FA and 0.02% TFA and each eluate was reduced to a volume of less than 300 µL and stored at − 20 °C until subjected to LC–MS/MS.
Processed samples were desalted using C-18 OMIX Pipette Tips (Agilent Technologies, USA) and eluted in 3 µL of MS buffer B (65% ACN, 0.1% FA in H2O). The eluates were then diluted with 57 µL of MS buffer A (0.1% FA in H2O) and 28 µL were injected onto a 2 cm C18 trap column, packed with Varian Pursuit (5 µm C18), with an 8 µm tip (New Objective). The LC setup was coupled online to a Q Exactive (Thermo Fisher Scientific, USA) mass spectrometer with a nanoelectrospray ionization source (Proxeon Biosystems). Samples for direct proteomic analysis as well as samples for direct peptidomics analysis underwent a 60-min linear gradient using MS buffer A and MS buffer B. Eluted peptides were subjected to tandem mass spectrometry in positive ion mode. Data acquisition was conducted via Thermo XCalibur v.3.0.63 (Thermo Fisher Scientific, USA).
Protein identification and quantification
The resulting proteomic and peptidomic raw data files were uploaded into MaxQuant v.126.96.36.199 (www.coxdocs.org)  with the integrated Andromeda search engine. MS and MS/MS spectra were searched against a reverted version of the SwissProt human protein database (version July 2017) for protein identification and a randomized version of the SwissProt human protein database for peptide identification. Search parameters for proteomic analysis included carbamidomethylation of cysteines as a fixed modification and oxidized methionine and N-terminal acetylation as variable modifications. Trypsin was the chosen digestion enzyme and a maximum of two missed cleavages were allowed. Search parameters for peptidomic analysis included oxidized methionine and oxidized proline as variable modifications. An unspecific enzyme search was the chosen digestion method. Both proteomic and peptidomic data were initially searched against a smaller “human first search” database with a peptide tolerance of 20 ppm for mass recalibration. The main search was performed using the Swissprot human protein database (version July 2017) with a peptide tolerance of 4.5 ppm. Data was analyzed using label-free quantification (LFQ) with a minimum ratio count of 1 and the “Match between runs” interval set to 2 min. The peptide-spectrum match and protein false discovery rate were set to 1%.
Pathway analysis of dysregulated proteins identified by LC–MS/MS was conducted using the functional-analysis tool Ingenuity Pathway Analysis (IPA; http://www.ingenuity.com) . To determine the specificity of identified proteins at the tissue and biological fluid level, proteomic datasets were searched against ProteomicsDB (http://www.ProteomicsDB.org), a web-based database of mass spectrometry-generated proteomics data . Pathway analysis of SF peptides was conducted through the Database for Annotation, Visualization and Integrated Discovery (DAVID) 6.8 with reference to the Kyoto Encyclopedia of Genes and Genomes (KEGG) . Annotations with q-values of less than 0.05 were considered statistically significant. Identification of known AMPs in the SF peptidome was determined by comparison with experimentally validated human AMPs taken from the Collection of Anti-Microbial Peptides (CAMPR3) (http://www.camp.bicnirrh.res.in/) database . AMP prediction of the identified peptides was performed using the support vector machine (SVM) learning algorithm developed for the CAMPR3 database. Peptides with an SVM score of 0.8 or higher were predicted to be antimicrobial.
Statistical analyses and data visualizations were completed with R (R Foundation for Statistical Computing). A linear model was fitted to examine the effects of age and sex on the protein and peptide expression data using the LIMMA package in R . Due to the nature of data generated by LC–MS/MS, protein quantification is often skewed and imposes limits on statistical inference. To circumvent assumptions of normality, the Mann–Whitney U test coupled to the Benjamini–Hochberg correction to control for multiple hypothesis testing was performed for comparisons between two independent groups. Adjusted p values of less than 0.05 were considered statistically significant. Differential abundance of proteins and peptides were computed with the myTAI package in R, generating a ratio of log-transformed extracted ion currents in one group against the second group, considered to be the reference group . A volcano plot was used to visualize the results of the Mann–Whitney U test.
Clinical characteristics of recruited patients
Demographics, disease characteristics and concomitant therapies of recruited patients are summarized in Table 1.
Holistic protein and peptide mining
Collectively, 389 unique proteins were identified across all IA SF proteomic samples. When assessing each cohort individually, 377 unique proteins were identified in RA patient samples, 369 unique proteins were identified in PsA patient samples and 399 proteins were identified in control patient samples. A review of the overlap between proteomes of each cohort revealed 347 proteins to be common to all three patient groups.
A total of 226 unique peptide sequences were identified across all IA SF samples originating from a total of 48 unique proteins. Inter-cohort comparisons identified 184 unique peptides in RA patient samples, 175 unique peptides in PsA patient samples and 192 unique peptides in control patient samples. Comparisons between the SF peptidomes of arthritic and control conditions revealed 95 peptides to be common to all three groups.
Next, we investigated the overlap between the proteins identified through our peptidomic approach and those identified through our proteomic approach by comparing the IA-associated proteins originating from both experiments. Of the 48 precursor proteins from our peptidomic study, 25 proteins were also found in the IA SF proteome (Fig. 1). Taken together, they have yielded the combined identification of 412 proteins in IA SF. A complete list of identified proteins and peptides are reported in Additional file 1: Tables S1, S2 and S3.
Dysregulated proteins in IA SF
Differential abundance analyses were conducted to detect dysregulated proteins in the SF of: (1) IA compared to control and (2) RA compared to PsA. Using non-parametric statistical tests, 144 proteins were determined to have statistically significant differential abundance in IA SF with 54 proteins showing significant upregulation and 90 proteins showing significant downregulation (Fig. 2). When comparing RA and PsA proteomes, no proteins showed significant differences in abundance after correcting for multiple hypothesis testing. However, with respect to an unadjusted p value, 22 proteins were differentially abundant between the two groups with 13 proteins demonstrating significant upregulation in RA relative to PsA and 9 proteins showing significant upregulation in PsA relative to RA. Significantly dysregulated proteins in IA compared to control and significantly dysregulated proteins in RA compared to PsA are summarized in Additional file 1: Tables S4 and S5, respectively.
Dysregulated functional pathways likely to be associated with the significantly upregulated and downregulated proteins of IA SF were detected with IPA. Core analyses determined the top 5 canonical pathways associated with upregulated proteins to be: (1) LXR/RXR activation, (2) FXR/RXR activation, (3) acute phase response signaling, (4) atherosclerosis signaling and (5) IL-12 signaling and production in macrophages, several of which have been previously associated with IA. Details regarding the top diseases and disorders as well as molecular and cellular functions associated with both groups of dysregulated proteins can be found in Additional file 1: Table S6. As the data suggests, upregulated proteins show more distinct relations to inflammatory and immunological processes while downregulated proteins demonstrate stronger relations to metabolic processes. Ultimately, to identify the strongest candidate biomarkers to be validated in IA patient serum, we focused on upregulated proteins in the SF.
Tissue and fluid specificity of upregulated proteins were used to narrow down the list of candidates deemed likely to be associated with IA, RA and PsA. We concentrated on proteins which displayed strong abundance in SF, bone, bone marrow or immune regulatory cells according to ProteomicsDB. Immunoglobulins were excluded from further analysis. The resulting list of upregulated proteins compared to the reference group consisted of 38 IA-specific, 8 RA-specific and 9 PsA-specific unique protein candidates. High abundance proteins in serum, as identified according to the literature [21, 22], were excluded due to the likelihood that they were serum contaminants at the time of joint aspiration. Moreover, this ensured protein candidates were unlikely to be overexpressed in the serum of non-diseased patients. Following this filtering step, the final list of candidate biomarkers consisted of 5, 4 and 2 upregulated proteins which we deemed likely to be associated with IA, RA and PsA, respectively (Table 2).
Dysregulated peptides in IA SF
Differential abundance analyses were conducted to detect strongly dysregulated peptides in the SF of: (1) IA compared to control and (2) RA compared to PsA. For both comparisons, no peptides showed statistically significant differences in abundance after correcting for multiple hypothesis testing, with the exception of the peptide sequence DSGEGDFLAEGGGV when comparing IA to the control. Alternatively, with respect to the unadjusted p value, 11 peptides were determined to be significantly differentially abundant in IA SF with 10 peptides showing significant upregulation and 1 peptide showing significant downregulation (Table 3). A complete list of dysregulated peptides in IA compared to control and dysregulated peptides in RA compared to PsA are summarized in Additional file 1: Tables S7 and S8, respectively.
When comparing RA and PsA peptidomes, 5 peptides showed differential abundance between the two groups with all 5 peptides demonstrating significant upregulation in PsA SF relative to RA SF (Table 4).
Pathway enrichment analysis of the SF peptidome
KEGG analysis revealed significantly enriched pathways (fold enrichment in brackets) related to the IA SF peptidome. Figure 3 illustrates the top KEGG pathways among which complement and coagulation cascades , Staphylococcus aureus infection , protein digestion and absorption  and extracellular matrix (ECM)-receptor interaction  were significantly enriched.
Antimicrobial peptides in IA SF
Accumulating evidence suggests a crucial role of intestinal resident flora in chronic activation of innate and adaptive immune responses leading to inflammatory disorders. Microorganisms residing in the intestine play an important role in maintaining systemic homeostasis through the delicate balance of the immune system response. Perturbations in the composition of the intestinal microbiota have been shown to elicit inappropriate immune cell activation leading to an inflammatory cascade and eventually, clinical disease . Specifically, perturbations of the gut epithelial cell layer and/or increased exposure to microbial metabolites may be primary triggers of an inflammatory cascade resulting in joint pathology . Protective mechanisms, such as the expression of AMPs, have naturally developed to oppose microbial dysbiosis. AMPs are fundamental effectors of the innate immune response with a broad spectrum of microbicidal activity. Under inflammatory conditions, the synovial membrane has demonstrated an altered pattern of expression of AMPs relative to healthy controls and suggests a valuable role of these proteins in the differential diagnosis of inflammatory joint disease .
Putative AMPs in the SF peptidome of IA were predicted with the assistance of a SVM learning algorithm (Additional file 1: Table S9). Overall, 26 peptide sequences originating from 8 proteins (complement C4-A, fibrinogen beta chain, fibrinogen alpha chain, annexin A1, collagen type III alpha 1 chain, collagen type I alpha 1 chain, gliomedin and EMI domain-containing protein (1) were predicted to have antimicrobial activity with an SVM score of 0.8 or higher (Table 5).
In the current study, a comparative MS-based approach coupled to statistical and bioinformatics analyses was performed on IA SF relative to control SF, and RA SF relative to PsA SF, to detect notable differences in both the proteomic and peptidomic data. Studies using an MS-based approach to evaluate the proteome of similar inflammatory diseases, including psoriasis , systemic lupus erythematosus , and ankylosing spondylitis , corroborate the robustness of such analytical methodologies. The investigation of a proximal joint fluid, such as SF, was preferred since its protein and peptide expression patterns are expected to be reflective of the pathophysiological state of the joint. As such, elucidating the SF proteome and peptidome during the progression of IA can provide novel insights into molecular drivers of the disease.
The molecular pathways involved in the pathogenesis of IA are also overrepresented in the current study based on functional network analysis of IA SF proteins and peptides. Prominent mechanisms related to the identification of upregulated proteins include: (1) acute phase response signaling, (2) antimicrobial response, (3) inflammatory response, (4) IL-12 signaling and production in macrophages and (5) cell-to-cell signaling and interaction. Similarly, interaction networks were established through pathway enrichment analysis of IA SF peptides. Of interest was the enrichment of Staphylococcus aureus infection. As previously highlighted, correlative studies are beginning to recognize a fundamental interplay between the microbiome and immune system response in the etiology of IA [29, 30]. Although the role of S. aureus in the progression of IA has yet to be clarified, the enrichment of this pathway, as reflected by the peptides identified in our study, reinforces this developing hypothesis.
Overall, our analyses resulted in the identification of 144 differentially expressed proteins in the IA SF proteome. Comparison of RA SF to PsA SF identified 22 differentially expressed proteins. Since we are interested in identifying putative markers which can be further validated in patient serum, we decided to focus solely on upregulated proteins in each arthritic condition. High-potential candidate biomarkers were selected on the basis of several molecular features including: differential abundance, fluid and tissue specificity, immunoglobulin status and abundance in the plasma proteome. Our list of dysregulated proteins in IA was reduced to a total of 5 promising protein candidates representative of intrinsic joint structures including the articular cartilage, synovial membrane and synoviocytes. The re-discovery of several upregulated proteins which have been previously implicated in IA, such as CD5 molecule-like (CD5L), matrix metalloproteinase (MMP)-3, defensin alpha 3 (DEFA3), S100 calcium-binding protein (S100) A8, and A9, provided an internal validation of our analytical proteomic approach [31, 32]. The application of similar, stringent filtering criteria on protein candidates of RA and PsA resulted in 4 RA-specific and 2 PsA-specific promising protein candidates.
Our analytical approach also yielded the discovery of novel putative biomarkers which, to our knowledge, have yet to be described in the context of IA. This includes the identification of alpha-ketoglutarate-dependent dioxygenase (FTO), family with sequence similarity 21 member C (FAM21C; more commonly known as WASH complex subunit 2C, WASHC2C) and T-box transcription factor (TBX3). Of these candidates, only TBX3 has been previously observed in IA at the genetic level . A genome-wide association study (GWAS) identified the single nucleotide polymorphism (SNP), rs12579024, located nearest the TBX3 gene, to be strongly associated with RA in a Korean population (p value < 0.0001). The functional roles of TBX3 have, thus far, been primarily described in relation to the morphogenesis of limbs and organs  as well as oncogenic processes . A recent study by Willmer et al.  attempted to delineate the molecular mechanisms driven by TBX3 and identified cyclin-dependent kinase inhibitor p21WAF1 (p21), a key mediator of cell cycle arrest, to be a primary repressed target of TBX3. Interestingly, p21 has also been implicated in the regulation of proinflammatory cytokines and MMP production in synovial fibroblasts, both of which greatly promote inflammation and joint destruction during the pathogenesis of RA . Isolated RA synovial fibroblasts have shown reduced expression of p21 relative to osteoarthritis (OA) synovial fibroblasts and adenovirus-mediated delivery of p21 suppresses the spontaneous production of IL-6 and MMP1 in RA synovial fibroblasts. In support of this, p21−/− mice maintain an enhanced experimental IA with markedly increased numbers of macrophages and articular destruction . This phenotype is resolved, however, with the administration of a p21-peptide mimetic. When taken with our own findings, it is conceivable that the upregulation of TBX3 in the synovial joint may lead to reduced p21 expression in synovial fibroblasts and promotes the proinflammatory state distinctive of IA pathogenesis. These findings corroborate with our hypothesis that delineating the IA proteome may highlight underlying mechanisms related to the progression of inflammatory arthritic disease and serve as novel targets for screening and therapeutic purposes.
Comparisons of RA and PsA revealed high-priority protein candidates specific to each disease. In RA SF, coagulation factor XII, SPARC-like protein 1, Rab GDP dissociation inhibitor beta and immunoglobulin gamma Fc region receptor III-A (FCGR3A) were notably upregulated; of which, activating FCGR3A has demonstrated important roles in sustaining the inflammatory response through the secretion of cytokines and proteases from the immune cell on which it is expressed . Likewise, allelic studies have demonstrated SNPs that may serve as susceptibility markers for RA . Taken together, the therapeutic targeting of FCGR3A may facilitate future management of RA.
Of the two PsA-specific protein candidates we identified, periostin (POSTN) has been previously investigated in our studies of the PsA tissue proteome as a potential serum marker of PsA . Although serum validation of POSTN did not reveal statistically significant differences between PsA and control serum, its elevated levels in both PsA lesional skin as well as SF alludes to an important role of the protein in the pathobiology of PsA and may serve as part of a panel of biomarkers to differentiate between the onset of PsA and RA.
Differential abundance analyses of peptide sequences identified 11 peptides to be significantly dysregulated in IA SF compared to the control group. Upregulated peptide sequences were primarily derived from FGA while single sequences originated from collagen type I alpha 1 (COL1A1) and coiled-coil serine rich protein 2 (CCSER2). All significant FGA-derived peptide fragments were representative of the 16-amino acid residue (ADSGEGDFLAEGGGVR) of fibrinopeptide A (FpA) located at the NH2-terminal end of FGA. The lack of detection of the full-length FpA peptide sequence in IA SF can be rationalized by the peptide’s short half-life of 3–5 min in the blood plasma . FpA is a cleavage product of thrombin-induced conversion of fibrinogen into a fibrin clot. Fibrin deposition in the SF or on the synovial membrane is a consistent feature of IA and is believed to perpetuate inflammation and joint tissue destruction through synovial cell activation [43, 44]. Liu et al. demonstrated that stimulation of synovial fibroblasts with fibrin(ogen) resulted in the upregulated expression of IL-8 and intercellular adhesion molecule 1 (ICAM-1) for the recruitment and retention, respectively, of lymphocytes within the arthritic joint . Elevated abundance of FGA and FpA in serum has been observed in patients with inflammation-associated diseases including systemic lupus erythematosus, Crohn’s disease, ischemic heart disease and gastric cancer [45,46,47,48]. These findings highlight the non-specific indication of inflammation by FpA and its associated peptide fragments, and advocates for its utility as a sensitivity index of disease activity in patients with IA. Moreover, targeting FGA in the synovial joint may be a necessary therapeutic intervention to modulate the inflammatory response. Comparisons of peptide abundance between RA and PsA identified FGA and FGB-related peptide sequences to be consistently downregulated in RA relative to PsA. Although this may be indicative of a discriminatory ability for FGA and FGB peptide fragments to differentiate between the onset of RA and PsA, this outcome does not corroborate with the finding that RA patients are at a greater increased risk of venous thromboembolism relative to PsA patients . Targeted quantification in a second set of SF samples is necessary to verify this finding.
The advent of high-throughput microbial DNA sequencing has marked a renewed interest in the complex interplay of the intestinal microbiome and inflammatory diseases. Studies suggest that the induction of autoimmunity is closely linked to intestinal dysbiosis and leads to distal synovitis and joint pathology . There exist several protective mechanisms to prevent changes in the gut microbiota including the physicochemical barrier of antimicrobial proteins and peptides. AMPs are a collective of naturally-occurring, cationic peptides released by lymphocytes of the innate immune system. Of the 26 peptides predicted to have antimicrobial activity, 13 of them originated from FGA or FGB precursor proteins. Despite the pro-inflammatory impression associated with the accumulation of FGA and FGB in the SF, their presence may be critical to the activation of microbicidal activity. Soluble fibrinogen and fibrin matrices have demonstrated antimicrobial host defense through their ability to physically entrap bacteria in addition to the recruitment and engagement of host immune cells which in turn, facilitate the removal of invading pathogens . Taken together, the deposition of fibrin during the progression of IA may initially serve the favourable purpose of limiting bacterial infection through the activation of antimicrobial host defense mechanisms. However, its added role in the recruitment and activation of leukocytes may exacerbate synovial joint inflammation thereby fueling joint disease.
Though these findings are limited by lack of verification in a subsequent set of SF samples, the identification of IA-specific candidates using a label-free, MS-based approach has shown biological relevance and prospective utility for clinical applications. Future follow-up studies will address verification and validation efforts of selected protein and peptide candidates in a new set of SF and serum samples, respectively. We do acknowledge the limitation of sex discrepancy amongst the IA SF samples in our study which may have influenced the proteins and peptides identified. However, to compensate for this discrepancy between each subtype of IA, our control group consisted of an equal number of male and female SF samples. Moreover, we tested the influence of both sex and age on our data using a linear model and found there to be no effect by either predictor.
A technical limitation of this study includes the lack of fractionation of digested proteins and peptides which likely contributed to the low fold change ratios of our candidate biomarkers. Pre-fractionation methods are important for reducing the complexity of biological fluids and tissues. The proteomic profile of IA SF is markedly shifted compared to healthy SF with a greater concentration of pro-inflammatory cytokines, immunoglobulins, matrix-degrading enzymes and acute-phase markers. The dynamic range between proteins in diseased SF can vary by a factor of 1010  and the likelihood, therefore, of masking potentially clinically-relevant proteins within the low-abundance proteome increases and may be exacerbated by analysis of unfractionated biological samples. However, improving accessibility to low-concentration proteins comes at the cost of longer analysis times and lower reproducibility . Moreover, fractionation technologies have previously failed to significantly extend the sampling of the proteome relative to the unfractionated proteome .
Chronic inflammation in IA is orchestrated by a complex network of signaling pathways which are expected to be represented in the protein and peptide expression patterns of SF. Therefore, proteomic and peptidomic analysis of SF can reflect the molecular underpinnings of IA and enhance our understanding of principal drivers at the apex of this disease. Overall, through the application of high-throughput, label-free MS, this discovery-phase study has generated a comprehensive proteomic dataset representative of IA SF and its specific subtypes. We discovered 5 protein candidates and 10 peptide candidates upregulated in IA SF, of which 3 proteins have yet to be described in IA. Moreover, subtype-specific analyses identified 4 RA-specific protein candidates, 2 PsA-specific protein candidates and 5 PsA-specific peptide candidates. Several of these candidates have been associated with inflammatory pathways at the genetic level but have not been investigated at the protein level and therefore, require functional experimentation to elucidate their role in the pathogenesis of IA. The data presented herein underscores the potential for proteins and peptides to elucidate mechanistic pathways related to the onset of arthritic disease in addition to their capacity to serve as informative clinical biomarkers.
Availability of data and materials
The mass spectrometry proteomics and peptidomics datasets supporting the conclusions of this article are available in the PRIDE Archive via the PRIDE partner repository with the data set identifier PXD011872; http://www.ebi.ac.uk/pride/archive/ (username: firstname.lastname@example.org and password: 3hXihB2 s).
major histocompatibility complex
American College of Rheumatology
classification criteria for psoriatic arthritis
liquid chromatography-tandem mass spectrometry
ingenuity pathway analysis
database for annotation, visualization and integrated discovery
Kyoto Encyclopedia of Genes and Genomes
- CAMPR3 :
collection of anti-microbial peptides
support vector machine
fibrinogen alpha chain
fibrinogen beta chain
tumor necrosis factor alpha
S100 calcium-binding protein
defensin alpha 3
family with sequence similarity 21 member C
T-box transcription factor
genome-wide association study
single nucleotide polymorphism
cyclin-dependent kinase inhibitor p21WAF1
immunoglobulin gamma Fc region receptor III-A
phosphoglycerate kinase 1
collagen type I alpha 1
coiled-coil serine rich protein 2
intercellular adhesion molecule 1
Biniecka M, Connolly M, Gao W, Ng CT, Balogh E, Gogarty M, et al. Redox-mediated angiogenesis in the hypoxic joint of inflammatory arthritis. Arthritis Rheumatol. 2014;66(12):3300–10.
Hazes JMW, Luime JJ. The epidemiology of early inflammatory arthritis. Nat Rev Rheumatol. 2011;7:381.
Ritchlin CT, Colbert RA, Gladman DD. Psoriatic arthritis. N Engl J Med. 2017;376(10):957–70.
Symmons DP, Bankhead CR, Harrison BJ, Brennan P, Barrett EM, Scott DG, et al. Blood transfusion, smoking, and obesity as risk factors for the development of rheumatoid arthritis: results from a primary care-based incident case-control study in Norfolk, England. Arthritis Rheum. 1997;40(11):1955–61.
Zhang X, Zhang D, Jia H, Feng Q, Wang D, Liang D, et al. The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment. Nat Med. 2015;21(8):895–905.
Scher JU, Ubeda C, Artacho A, Attur M, Isaac S, Reddy SM, et al. Decreased bacterial diversity characterizes the altered gut microbiota in patients with psoriatic arthritis, resembling dysbiosis in inflammatory bowel disease. Arthritis Rheum (Hoboken, NJ). 2015;67(1):128–39.
Ogdie A, Schwartzman S, Husni ME. Recognizing and managing comorbidities in psoriatic arthritis. Curr Opin Rheumatol. 2015;27(2):118–26.
Solomon DH, Karlson EW, Rimm EB, Cannuscio CC, Mandl LA, Manson JE, et al. Cardiovascular morbidity and mortality in women diagnosed with rheumatoid arthritis. Circulation. 2003;107(9):1303–7.
Greening DW, Kapp EA, Simpson RJ. The peptidome comes of age: mass spectrometry-based characterization of the circulating cancer peptidome. The Enzymes. 2017;42:27–64.
Gibson DS, Rooney ME. The human synovial fluid proteome: a key factor in the pathology of joint disease. Proteom Clin Appl. 2007;1(8):889–99.
Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31(3):315–24.
Taylor W, Gladman D, Helliwell P, Marchesoni A, Mease P, Mielants H. Classification criteria for psoriatic arthritis: development of new criteria from a large international study. Arthritis Rheum. 2006;54(8):2665–73.
Kamphorst JJ, van der Heijden R, DeGroot J, Lafeber FP, Reijmers TH, van El B, et al. Profiling of endogenous peptides in human synovial fluid by NanoLC-MS: method validation and peptide identification. J Proteome Res. 2007;6(11):4388–96.
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotech. 2008;26(12):1367–72.
Kramer A, Green J, Pollard J Jr, Tugendreich S. Causal analysis approaches in ingenuity pathway analysis. Bioinformatics (Oxford, England). 2014;30(4):523–30.
Schmidt T, Samaras P, Frejno M, Gessulat S, Barnert M, Kienegger H, et al. ProteomicsDB. Nucleic Acids Res. 2018;46(D1):D1271–81.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Waghu FH, Barai RS, Gurung P, Idicula-Thomas S. CAMPR3: a database on sequences, structures and signatures of antimicrobial peptides. Nucleic Acids Res. 2016;44(D1):D1094–7.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Drost HG, Gabel A, Liu J, Quint M, Grosse I. myTAI: Evolutionary Transcriptomics with R. Bioinformatics (Oxford, England). 2017.
Hortin GL, Sviridov D, Anderson NL. High-abundance polypeptides of the human plasma proteome comprising the top 4 logs of polypeptide abundance. Clin Chem. 2008;54(10):1608–16.
Fountoulakis M, Juranville JF, Jiang L, Avila D, Röder D, Jakob P, et al. Depletion of the high-abundance plasma proteins. Amino Acids. 2005;27:249–59.
Paulsen F, Pufe T, Conradi L, Varoga D, Tsokos M, Papendieck J, et al. Antimicrobial peptides are expressed and produced in healthy and inflamed human synovial membranes. J Pathol. 2002;198(3):369–77.
Scher JU, Abramson SB. The microbiome and rheumatoid arthritis. Nat Rev Rheumatol. 2011;7(10):569–78.
Diamanti AP, Manuela Rosado M, Laganà B, D’Amelio R. Microbiota and chronic inflammatory arthritis: an interwoven link. J Transl Med. 2016;14:233.
Fattahi S, Kazemipour N, Hashemi M, Sepehrimanesh M. Alpha-1 antitrypsin, retinol binding protein and keratin 10 alterations in patients with psoriasis vulgaris, a proteomic approach. Iran J Basic Med Sci. 2014;17(9):651–5.
Kazemipour N, Qazizadeh H, Sepehrimanesh M, Salimi S. Biomarkers identified from serum proteomic analysis for the differential diagnosis of systemic lupus erythematosus. Lupus. 2015;24(6):582–7.
Fischer R, Trudgian DC, Wright C, Thomas G, Bradbury LA, Brown MA, et al. Discovery of candidate serum proteomic and metabolomic biomarkers in ankylosing spondylitis. Mol Cell Proteom MCP. 2012;11(2):M111.013904-M111.
Eerola E, Mottonen T, Hannonen P, Luukkainen R, Kantola I, Vuori K, et al. Intestinal flora in early rheumatoid arthritis. Br J Rheumatol. 1994;33(11):1030–8.
Rogier R, Evans-Marin H, Manasson J, van der Kraan PM, Walgreen B, Helsen MM, et al. Alteration of the intestinal microbiome characterizes preclinical inflammatory arthritis in mice and its modulation attenuates established arthritis. Sci Rep. 2017;7:15613.
Hansson C, Eriksson C, Alenius G-M. S-Calprotectin (S100A8/S100A9): a potential marker of inflammation in patients with psoriatic arthritis. J Immunol Res. 2014;2014:5.
Cretu D, Prassas I, Saraon P, Batruch I, Gandhi R, Diamandis EP, et al. Identification of psoriatic arthritis mediators in synovial fluid by quantitative mass spectrometry. Clin Proteom. 2014;11(1):1–12.
Freudenberg J, Lee HS, Han BG, Shin HD, Kang YM, Sung YK, et al. Genome-wide association study of rheumatoid arthritis in Koreans: population-specific loci as well as overlap with European susceptibility loci. Arthritis Rheum. 2011;63(4):884–93.
Bamshad M, Lin RC, Law DJ, Watkins WC, Krakowiak PA, Moore ME, et al. Mutations in human TBX3 alter limb, apocrine and genital development in ulnar-mammary syndrome. Nat Genet. 1997;16(3):311–5.
Peres J, Prince S. The T-box transcription factor, TBX3, is sufficient to promote melanoma formation and invasion. Mol Cancer. 2013;12(1):117.
Willmer T, Hare S, Peres J, Prince S. The T-box transcription factor TBX3 drives proliferation by direct repression of the p21WAF1 cyclin-dependent kinase inhibitor. Cell Div. 2016;11(1):6.
Perlman H, Bradley K, Liu H, Cole S, Shamiyeh E, Smith RC, et al. IL-6 and matrix metalloproteinase-1 are regulated by the cyclin-dependent kinase inhibitor p21 in synovial fibroblasts. J Immunol (Baltimore, Md: 1950). 2003;170(2):838–45.
Mavers M, Cuda CM, Misharin AV, Gierut AK, Agrawal H, Weber E, et al. The Cyclin dependent kinase inhibitor p21 is essential for resolution of murine inflammatory arthritis via its C-terminal domain. Arthritis Rheum. 2012;64(1):141–52.
Abrahams VM, Cambridge G, Lydyard PM, Edwards JC. Induction of tumor necrosis factor alpha production by adhered human monocytes: a key role for Fcgamma receptor type IIIa in rheumatoid arthritis. Arthritis Rheum. 2000;43(3):608–16.
Morgan AW, Barrett JH, Griffiths B, Subramanian D, Robinson JI, Keyte VH, et al. Analysis of Fcγ receptor haplotypes in rheumatoid arthritis: FCGR3A remains a major susceptibility gene at this locus, with an additional contribution from FCGR3B. Arthritis Res Ther. 2005;8(1):R5.
Cretu D, Liang K, Saraon P, Batruch I, Diamandis EP, Chandran V. Quantitative tandem mass-spectrometry of skin tissue reveals putative psoriatic arthritis biomarkers. Clin Proteom. 2015;12(1):1.
Manolis AS, Melita-Manolis H, Stefanadis C, Toutouzas P. Plasma level changes of fibrinopeptide A after uncomplicated coronary angioplasty. Clin Cardiol. 1993;16(7):548–52.
Liu X, Piela-Smith TH. Fibrin(ogen)-Induced Expression of ICAM-1 and chemokines in human synovial fibroblasts. J Immunol. 2000;165(9):5255.
Ho PP, Lee LY, Zhao X, Tomooka BH, Paniagua RT, Sharpe O, et al. Autoimmunity against fibrinogen mediates inflammatory arthritis in mice. J Immunol. 2010;184(1):379.
Hardin JA, Cronlund M, Haber E, Bloch KJ. Activation of blood clotting in patients with systemic lupus erythematosus. Am J Med. 1978;65(3):430–6.
Edwards RL, Levine JB, Green R, Duffy M, Mathews E, Brande W, et al. Activation of blood coagulation in Crohn’s disease. Increased plasma fibrinopeptide A levels and enhanced generation of monocyte tissue factor activity. Gastroenterology. 1987;92(2):329–37.
Xu S, Zhao J, Liu J, Gou W. Fibrinopeptide a induces expression of C-reactive protein via the ROS-ERK1/2/P38-NF-κB signal pathway in vascular smooth muscle cells. Cell Physiol Biochem. 2018;47(1):266–78.
Ebert MPA, Niemeyer D, Deininger SO, Wex T, Knippig C, Hoffmann J, et al. Identification and confirmation of increased fibrinopeptide a serum protein levels in gastric cancer sera by magnet bead assisted MALDI-TOF mass spectrometry. J Proteom Res. 2006;5(9):2152–8.
Ogdie A, Kay McGill N, Shin DB, Takeshita J, Jon Love T, Noe MH, et al. Risk of venous thromboembolism in patients with psoriatic arthritis, psoriasis and rheumatoid arthritis: a general population-based cohort study. Eur Heart J. 2017;39:3608–14.
van den Broek MF, van Bruggen MC, Koopman JP, Hazenberg MP, van den Berg WB. Gut flora induces and maintains resistance against streptococcal cell wall-induced arthritis in F344 rats. Clin Exp Immunol. 1992;88(2):313–7.
Ko Y-P, Flick MJ. Fibrinogen is at the interface of host defense and pathogen virulence in staphylococcus aureus infection. Semin Thromb Hemost. 2016;42(4):408–21.
Baiwir D, Mazzucchelli G, Smargiasso N, Quesada-Calvo F, De Pauw E, Malaise M, et al. Comparison of serum fractionation methods by data independent label-free proteomics. EuPA Open Proteom. 2015;9:14–22.
Whiteaker JR, Zhang H, Eng JK, Fang R, Piening BD, Feng LC, et al. Head-to-head comparison of serum fractionation techniques. J Proteom Res. 2007;6(2):828–36.
We thank Ihor Batruch for his support in mass spectrometric analysis.
This work was supported by The Krembil Foundation.
Ethics approval and consent to participate
Human research ethics board approval was received for the study from the University Health Network, Mount Sinai Hospital and the University of Calgary. Informed consent was obtained from all patients.
Consent for publication
The authors declare that they have no competing financial interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Mahendran, S.M., Keystone, E.C., Krawetz, R.J. et al. Elucidating the endogenous synovial fluid proteome and peptidome of inflammatory arthritis using label-free mass spectrometry. Clin Proteom 16, 23 (2019). https://doi.org/10.1186/s12014-019-9243-3
- Synovial fluid
- Mass spectrometry
- Inflammatory arthritis
- Rheumatoid arthritis
- Psoriatic arthritis