- Research article
- Open Access
Two independent proteomic approaches provide a comprehensive analysis of the synovial fluid proteome response to Autologous Chondrocyte Implantation
Arthritis Research & Therapyvolume 20, Article number: 87 (2018)
Autologous chondrocyte implantation (ACI) has a failure rate of approximately 20%, but it is yet to be fully understood why. Biomarkers are needed that can pre-operatively predict in which patients it is likely to fail, so that alternative or individualised therapies can be offered. We previously used label-free quantitation (LF) with a dynamic range compression proteomic approach to assess the synovial fluid (SF) of ACI responders and non-responders. However, we were able to identify only a few differentially abundant proteins at baseline. In the present study, we built upon these previous findings by assessing higher-abundance proteins within this SF, providing a more global proteomic analysis on the basis of which more of the biology underlying ACI success or failure can be understood.
Isobaric tagging for relative and absolute quantitation (iTRAQ) proteomic analysis was used to assess SF from ACI responders (mean Lysholm improvement of 33; n = 14) and non-responders (mean Lysholm decrease of 14; n = 13) at the two stages of surgery (cartilage harvest and chondrocyte implantation). Differentially abundant proteins in iTRAQ and combined iTRAQ and LF datasets were investigated using pathway and network analyses.
iTRAQ proteomic analysis confirmed our previous finding that there is a marked proteomic shift in response to cartilage harvest (70 and 54 proteins demonstrating ≥ 2.0-fold change and p < 0.05 between stages I and II in responders and non-responders, respectively). Further, it highlighted 28 proteins that were differentially abundant between responders and non-responders to ACI, which were not found in the LF study, 16 of which were altered at baseline. The differential expression of two proteins (complement C1s subcomponent and matrix metalloproteinase 3) was confirmed biochemically. Combination of the iTRAQ and LF proteomic datasets generated in-depth SF proteome information that was used to generate interactome networks representing ACI success or failure. Functional pathways that are dysregulated in ACI non-responders were identified, including acute-phase response signalling.
Several candidate biomarkers for baseline prediction of ACI outcome were identified. A holistic overview of the SF proteome in responders and non-responders to ACI has been profiled, providing a better understanding of the biological pathways underlying clinical outcome, particularly the differential response to cartilage harvest in non-responders.
Identification of putative biomarkers that can be used to predict patient outcome prior to treatment for cartilage injury has been highlighted as a key initiative for the prevention of osteoarthritis (OA) by the Osteoarthritis Research Society International . Further, in the United Kingdom, the National Health Service has increased the need to identify accurate prognostic biomarkers for application of the recent National Institute for Health and Care Excellence (NICE) recommendation for use of the cell therapy called autologous chondrocyte implantation (ACI) .
We recently published the first study , to our knowledge, in which a proteomic approach has been used with the aim of identifying candidate biomarkers to predict the success of ACI, a cellular therapy for the treatment of traumatic cartilage injury [4, 5]. This therapy is composed of a two-stage procedure: During the initial surgery (stage I), healthy cartilage is harvested from a minor load-bearing region of the joint, then chondrocytes are isolated and culture is expanded for 3–4 weeks prior to a second surgery (stage II), in which the chondrocytes are implanted into the cartilage defect [5, 6]. Approximately 500 patients have been treated with ACI in our centre, and despite an 81% success rate , we have yet to fully understand why some individuals do not respond well. We have identified a biomarker, aggrecanase-1, that, when its activity is undetectable pre-operatively, can be used together with known demographic and injury-associated risk factors to help predict ACI success [8, 9]. However, we have yet to identify a biomarker (or panel of biomarkers) that can be used to accurately predict ACI failure. The identification of such a biomarker(s) for ACI and other cartilage repair strategies would allow for the better stratification of patients prior to joint surgery and may provide candidates for therapies to improve ACI success.
Proteomic analyses remain one of the most widely used methods to identify novel biomarker candidates and have previously been used to identify biomarkers of OA progression (as summarised by Hsueh et al. in 2014 ). The synovial fluid (SF) provides an attractive biological fluid for biomarker identification because it bathes the injured joint and therefore contains proteins that might reflect the whole joint environment. Proteomic profiling of the SF, however, is technically difficult owing to the broad dynamic range of proteins present within it [7, 8]. Several unbiased global proteomic studies aimed at the identification of biomarkers within the SF have been completed. Nevertheless, the number of protein ‘hits’ has been somewhat limited, because researchers either have tended to profile SF with no pre-treatment to account for the wide range of proteins [11,12,13,14,15,16] or have depleted high-abundance proteins [17,18,19,20,21,22], meaning that the altered quantities of these proteins cannot be considered.
Isobaric tags for absolute and relative quantitation (iTRAQ) is reported to be the most accurate labelling method for quantifying comparative abundance of proteins . When compared with label-free quantitation (LF) proteomics, iTRAQ quantitation has traditionally been considered a more accurate technique ; however, as mass spectrometers are improved, these techniques are becoming more comparable, and LF is becoming increasingly popular . Unlike LF proteomics, iTRAQ uses isobaric tags to label the primary amines at the peptide level prior to pooling the samples to enable simultaneous identification and quantitation of the proteins. Fourplex and eightplex labels are available, enabling quantitation of up to eight conditions in a single analysis, thus minimising the number of mass spectrometry runs that can be cost-effective and time-efficient. However, when compared with LF, in which any number of samples can be analysed and compared, iTRAQ labelling limits the number of samples that can be compared, meaning biological replicate samples are often pooled together into relevant biological conditions. iTRAQ proteomics is a commonly used tool for the identification of biomarkers in a plethora of diseases. This proteomic approach has been used to profile the SF proteome [20, 26], successfully identifying differentially abundant protein biomarker candidates for several diseases/conditions.
Our previous study highlighted the potential of using protein equalisation to study low-abundance proteins in human SF, but this identified few differentially abundant proteins in baseline SF, when comparing individuals who did or did not do well following cartilage repair therapy . The aim of the present study therefore was to increase the number of protein biomarker candidates that could be identified for the pre-operative prediction of clinical outcome following ACI and to allow for the assessment of high-abundance proteins that may also strengthen the understanding of the biological processes underlying treatment success.
SF collection and storage
SF was collected as described previously [3, 8, 27] from the knee joints of patients who provided informed consent and following local research ethics committee approval. Immediately prior to both ACI surgeries, stage I (cartilage harvest) and stage II (chondrocyte implantation), 20 ml of saline was injected into the joint and 20 rounds of leg flexion and extension were carried out to allow aspiration of as much SF as possible [3, 27]. SF was then centrifuged at 6000 × g for 15 minutes at 4 °C and split into aliquots for long-term storage in liquid nitrogen. The dilution factor of the SF samples was calculated by comparing urea content in SF with matched blood plasma using a QuantiChrom™ Urea Assay Kit (BioAssay Systems, Hayward, CA, USA) according to the manufacturer’s instructions and as described previously [3, 8, 28], and SF samples with a dilution factor > 10 were excluded from the study.
Clinical responders to ACI were defined as individuals who demonstrated a Lysholm score increase of ≥ 10 points at 12 months post-treatment compared with their baseline score, as has been used previously [29,30,31]. The Lysholm score is a validated  patient-self assessment score encompassing knee pain and joint function that ranges from 0 to 100, with 100 representing ‘perfect’ knee function [32, 33]. Thirteen patients were considered as non-responders to ACI, demonstrating a mean decrease in Lysholm score of 14 points (range − 4 to − 46), and 14 SF donors were considered responders with a mean improvement of 33 points (range 17–54).
Sample preparation and analysis using iTRAQ proteomics (iTRAQ nanoLC-MS/MS)
Total protein was quantified using a Pierce™ 660 nm Protein Assay (Thermo Fisher Scientific, Hemel Hempstead, UK) , and a total of 200 μg of SF protein was pooled equally from the donors in each of the following experimental groups: stage I responders (n = 8), stage I non-responders (n = 7), stage II responders (n = 12), and stage II non-responders (n = 12) (Table 1). The pooled samples were then precipitated in six volumes of ice-cold acetone overnight at − 20 °C. The precipitates were pelleted by centrifugation at 13,000 × g for 10 minutes at 4 °C before being re-suspended in 200 μl of triethylammonium bicarbonate buffer. Eighty-five micrograms of protein for each experimental sample were then subjected to reduction, alkylation (as instructed in the iTRAQ labelling kit; Applied Biosystems, Bleiswijk, The Netherlands). Sequencing Grade Modified Trypsin (10 μg/85 μg of protein; Promega, Madison, WI, USA) was then added to the samples for overnight digestion at 37 °C. Tryptic digests were labelled with the iTRAQ tags according to the manufacturer’s instructions before being pooled into one microcentrifuge tube prior to being dried in a vacuum centrifuge: 114 tag- stage II responders, 115 tag- stage II non-responders, 116 tag- stage I responders and 117 tag- stage I non-responders.
iTRAQ-labelled peptides were resuspended in 0.6 ml of loading buffer (10 mM monopotassium phosphate [KH2PO4], 20% acetonitrile [MeCN], pH 3.0), followed by sonication. The pH was adjusted to 3.0 with 0.5 M orthophosphoric acid (H3PO4). The peptides were separated by strong cation exchange chromatography (SCX) as described previously . A total of 14 SCX fractions were analysed by nano-electrospray ionisation-LC-MS/MS using a TripleTOF 5600 tandem mass spectrometer (AB Sciex, Framingham, MA, USA) as described previously .
The raw mass spectrometry data file was subsequently analysed using ProteinPilot 4.5 software with the Paragon™ and ProGroup™ algorithms (AB Sciex) against the human sequences in the UniProtKB/Swiss-Prot database (downloaded in December 2012). Searches were performed using the pre-set iTRAQ settings in ProteinPilot. Trypsin was selected as the cleavage enzyme and methyl methanethiosulphonate for the modification of cysteines with a ‘thorough ID’ search effort. ProteinPilot’s bias correction assumes that most proteins do not change in expression. Finally, detected proteins were reported with a protein threshold [unused ProtScore (confidence)] > 0.05 and used in the quantitative analysis if they were identified with two or more unique peptides with 95% confidence or above. p Values and false discovery rates for the iTRAQ ratios were calculated using the ProteinPilot software. Proteins with iTRAQ ratios with p values ≤ 0.05 and with differential abundance of greater than or equal to ± 2.0-fold change (FC) were used in further analysis.
Verification of iTRAQ nanoLC-MS/MS results using enzyme linked immunosorbent assay
Two proteins of biological relevance were measured by enzyme-linked immunosorbent assay (ELISA) in the non-pooled samples to verify the MS findings. Firstly, complement C1S subcomponent (C1s) was selected because this protein demonstrated differential abundance between responders and non-responders to ACI within the baseline SF (prior to stage I surgery) and therefore could have potential as a biomarker of outcome prediction. C1s was assessed using a human ELISA (CUSABIO, Houston, TX, USA). Samples were first assayed using a 1:100 dilution in assay sample diluent, and for those samples that were undetectable in the assay, the assay was repeated using undiluted samples. Secondly, matrix metalloproteinase 3 (MMP3) was selected to investigate the differential response to stage I surgery (i.e., the proteomic shift between stages I and II) in non-responders to ACI. MMP3 was assessed using a human Quantikine® ELISA (R&D Systems, Abingdon, UK). Samples were diluted 1:100 in assay kit diluent prior to assessment. Both ELISAs were carried out according to the manufacturer’s instructions, and protein concentrations were normalised to the sample dilution factor. Statistical analysis was performed using Prism version 6.0 software (GraphPad Software, La Jolla, CA, USA). Student’s t tests were used to assess differential abundance.
Assessment of protein overlap identified using the two proteomic approaches
To assess whether the use of two independent proteomic approaches allows for a greater number of significant protein changes to be identified, the datasets from this study (iTRAQ nanoLC-MS/MS [nLC-MS/MS]) and our previously published study assessing the same patient samples (LF LC-MS/MS ) were compared with one another. Venn diagrams were plotted using VENNY 2.1.0 software  to assess the overlap of differentially abundant proteins that were identified via the two approaches.
Pathway and network analysis of proteomic datasets
The datasets generated using both proteomic approaches were combined. Specifically, proteins that were differentially expressed (≥ 1.2 FC; p ≤ 0.05) in each biological comparison (e.g., stage I responders versus non-responders) in either proteomic approach were merged into a single dataset. A modest FC cutoff was used to ensure that the greatest number of differentially abundant proteins could be included in the pathway and network analyses, as has been done previously [3, 18]. The iTRAQ nLC-MS/MS dataset independently, as well as when merged with the LF dataset, was analysed using pathway enrichment analysis (Ingenuity Pathway Analysis; Qiagen Bioinformatics, Redwood City, CA, US) to identify and visualise affected canonical pathways. Pathways with a significance level of p ≤ 0.005 were considered statistically significant (Fisher’s exact test).
The merged LF and iTRAQ nLC-MS/MS datasets of proteomic response to cartilage harvest (e.g., differential abundance between stages I and II) in responders and non-responders were assessed using interactome network analysis, which is an unbiased mathematical method of visualising and interpreting complex interactions between large numbers of molecules . Interactome networks are made up of nodes (the individual objects being studied, such as proteins) and edges (the connections between the objects, such as known protein-protein interactions) . By studying groups of proteins that are highly interconnected, known as modules, key functions within an interactome network can be highlighted . Conducting interactome network analysis alongside pathway enrichment analysis allows for greater confidence in the selection of candidate pathways or molecules for further study, because these represent two independent methods of mapping the data: known protein-protein interactions and text mining, respectively. The interactions between the differentially abundant proteins were assessed using the PINA4MS (Protein Interaction Network Analysis For Multiple Sets) app  in Cytoscape version 3.0 to generate network models based on protein-protein interactions. These models were based either on only those proteins identified in the proteomic analyses (non-inferred nodes) or on proteins identified in the proteomic analyses alongside their inferred interactions (inferred nodes) . The ModuLand (version 2.8.3) algorithm  was applied to the interactome networks in Cytoscape version 3.0 to identify highly connected clusters of proteins (modules) that demarcate the hierarchical structure of the interactome network. The biological function of each module was assessed by analysing the proteins identified within each module using the pathway analysis tool in Reactome software [43, 44]. The significance of the pathway functions identified in Reactome was determined by Fisher’s exact test, and p ≤ 0.05 was considered statistically significant.
The proteomic data derived from this study have been deposited in the PRoteomics IDEntifications (PRIDE) ProteomeXchange and can be accessed using the identifier [PXD008321].
Identification of proteins to predict ACI outcome prior to stage I or stage II
iTRAQ nLC-MS/MS highlighted 16 proteins (greater than or equal to ± 2.0 FC; p ≤ 0.05) that were differentially abundant between responders and non-responders to ACI at baseline (immediately prior to stage I) (Table 2). Prior to stage II of the ACI procedure, 12 proteins displayed differential abundance between responders and non-responders (Table 3).
At both stages of treatment, SF analysed using iTRAQ nLC-MS/MS identified a greater number of differentially abundant proteins between individuals who did or did not respond well to ACI compared with SF that had undergone protein normalisation using the ProteoMiner™ protein enrichment kit (Bio-Rad Laboratories, Hercules, CA, USA) and LF LC-MS/MS analysis . Further, the two proteomic techniques identified no common differentially abundant proteins. The two proteins selected and assessed by ELISA (C1s and MMP3) could verify the iTRAQ nLC-MS/MS (Fig. 1).
Differential abundance of proteins at stage II compared with stage I of ACI
Proteomic profiling of the SF using iTRAQ nLC-MS/MS highlighted a considerable effect of the cartilage harvest procedure (stage I) in both responders and non-responders, with 70 and 54 proteins being differentially abundant between stages I and II, respectively, thus strengthening the similar findings derived from the analysis of these samples using LF LC-MS/MS . Interestingly, the iTRAQ nLC-MS/MS and LF LC-MS/MS identified no common protein differences between stage I and stage II in the clinical responders (70 differentially abundant proteins identified by iTRAQ nLC-MS/MS and 14 identified by LF LC-MS/MS) (Table 4). This lack of overlap between the two proteomic techniques is highlighted in Fig. 2. There were, however, six proteins (gelsolin, vitamin K-dependent protein S, C4b-binding protein alpha chain, fibrinogen alpha chain, fibrinogen beta chain and fibrinogen gamma chain) that were identified by both proteomic techniques in the non-responders, all of which showed commonality in the direction of protein shift across the MS platforms, with iTRAQ nLC-MS/MS consistently resulting in greater differences in abundance than those identified from the LF LC-MS/MS data. A total of 54 protein abundance changes between stages I and II in non-responders were identified using iTRAQ nLC-MS/MS, and 55 protein differences were identified by LF LC-MS/MS (Table 5 and Fig. 2).
iTRAQ nLC-MS/MS confirmed a significant response to cartilage harvest procedure (stage I) in nonresponders to ACI
Pathway analysis of the iTRAQ nLC-MS/MS-identified proteins, using the pathway enrichment tools in Ingenuity Pathway Analysis, suggested that the proteins which were differentially abundant at stage II compared with stage I in non-responders are likely to impact numerous canonical pathways, many of which were confirmatory of the previously published functional pathways identified from the LF nLC-MS/MS-derived proteins . These functional pathways included acute-phase response signalling (p = 2.93 × 10− 1), the complement system (p = 2.11 × 10− 1) and liver X receptor/retinoic X receptor signalling (p = 1.95 × 10− 1). Moreover, many more functional pathways were affected as a result of the proteins that were differentially abundant in response to stage II compared with stage I in non-responders compared with responders (Additional file 1: Tables S1 and S2), reiterating that the SF proteomic response to cartilage harvest is more distinct in non-responders to ACI.
Similar pathways were identified from the differentially abundant proteins identified in iTRAQ nLC-MS/MS and LF LC-MS/MS analyses
Both iTRAQ nLC-MS/MS and LF LC-MS/MS analyses resulted in acute-phase response signalling being highlighted as one of the most significantly affected pathways in response to cartilage harvest in non-responders to ACI; therefore, this pathway was further assessed. Figure 3 highlights that analysis of the SF proteome using the two independent proteomic techniques resulted in a greater number of differentially abundant downstream proteins being identified. In addition, many complementary proteins have been identified when comparing these datasets, with the vast majority of proteins that are predicted to be increased in the plasma (the standard bodily fluid referred to in Ingenuity Pathway Analysis) during the acute-phase response being more abundant in the SF at stage II than at stage I and vice versa.
Because the results of the two proteomic approaches seem to be complementary to one another, the two datasets were combined to generate a more comprehensive profile of the SF proteome. Ingenuity Pathway Analysis again identified many functional pathways similar to those identified via the independent LF LC-MS/MS and iTRAQ nLC-MS/MS datasets. The most significant canonical pathways associated with the non-responder response to cartilage harvest (stage II versus stage I) were acute-phase response signalling (p = 1.10 × 10− 9), intrinsic prothrombin activation pathway (p = 3.43X10− 7) and the complement system (p = 1.22 × 10− 6). Further, analysis of upstream regulators to these dysregulated proteins included those identified using the LF LC-MS/MS analysis data alone, such as transforming growth factor-β1 (p = 2.05 × 10− 13), dihydrotestosterone (p = 4.48 × 10− 11) and peroxisome proliferator-activated receptor-α (p = 1.09 × 10− 9) .
The combined datasets were then used to generate unbiased interactome networks that represent the differentially abundant proteins (non-inferred networks), their likely interacting proteins (inferred networks) and how these proteins interact with one another, resulting in models of systemic protein response to cartilage harvest in either the responders or non-responders to ACI. Based on proteins that were differentially abundant between stages I and II of ACI in non-responders, an interactome network consisting of 115 nodes (proteins) and 40 edges (protein-protein interactions) was generated. Further, an inferred network consisting of 2893 proteins and 35,576 protein-protein interactions was generated on the basis of the addition of proteins that are likely to interact with the differentially abundant proteins (PINA4MS interactome database). Proteins that were differentially abundant in response to cartilage harvest in responders to ACI were used to generate interactome networks (non-inferred, 83 nodes and 118 edges; inferred, 2084 nodes and 54,007 edges). The ModuLand algorithm was applied to each of these networks to identify modules within the network that can be hierarchically ranked to identify groups of proteins that are the most fundamental in the functioning of the network. Figure 4 highlights the top ten modules from each of the networks generated. These modules again highlight the disparity between the ACI responder and non-responder response to cartilage harvest, with only modules centred on the proto-oncogene tyrosine-protein kinase (Src) protein being identified in the inferred networks of both non-responder and responder groups. Interestingly, assessment of the functional pathways related to the ModuLand identified modules in the non-responder networks again highlighted regulation of the complement cascade (p = 1.68 × 10− 8 by Fisher’s exact test), thus providing confidence in its importance based on identification via two independent bioinformatics approaches.
On the basis of its recent technology appraisal of ACI, NICE has recommended this treatment for a specific subset of patients with cartilage injury in the knee . The identification of novel biomarkers that can strengthen current patient demographic risk factors in predicting clinical outcomes , as well as development of a greater understanding of the underlying biology associated with success and failure, will be beneficial, particularly because this treatment option is likely to be implemented on a wider scale in the near future. The present study builds upon our previously published work [3, 8], highlighting a number of novel protein candidates that have potential as biomarkers to predict ACI outcome. Moreover, comprehensive proteomic profiling of SF has further highlighted proteomic differences between responders and non-responders to ACI.
In the majority of studies in which the SF proteome has been profiled, either high-abundance proteins [11,12,13,14,15,16] or low-abundance proteins [17,18,19,20,21,22] have been assessed via depletion or non-depletion of abundant proteins prior to proteomic analysis. Our study highlights that the use of both a proteomic dynamic range compression technique (ProteoMiner™)  in tandem with analysis of non-depleted SF samples can provide a more holistic overview of proteomic changes, because both iTRAQ nLC-MS/MS and LF LC-MS/MS highlighted large numbers of differentially abundant proteins between stages I and II of ACI, with little crossover between techniques. This type of all-inclusive approach to unbiased whole-proteome analysis of biological fluids may therefore be more successful in the identification of candidate biomarkers for treatments/disease states beyond those we investigated.
A limitation of our previous study  was that very few proteins were identified as differentially abundant between responders and non-responders at baseline. In order for biomarkers aimed at predicting ACI success to be most clinically useful, patients who are likely to fail or respond to this procedure need to be identified prior to any surgical intervention. Interestingly, analysis of non-dynamic range compressed proteins with iTRAQ nLC-MS/MS analysis was able to detect a greater number of differentially abundant proteins between responders and non-responders prior to stage I surgery. The protein with most altered abundance in responders compared with non-responders at stage I was C1s. This higher abundance in responders was confirmed in individual patient samples using a biochemical assay. C1s is a major constituent of the trimeric complement C1 protein, which triggers the classical complement pathway. Once activated, the classical complement pathway promotes inflammation to enable the removal of damaged cells and/or microbes. Moreover, C1s has been shown to cleave insulin-like growth factor 1 (IGF-1)  and insulin-like growth factor binding protein 5 (IGFBP-5) . Both IGF-1 and IGFBP-5 are chondroprotective when in their intact state [45, 47], and inhibition of C1s activity within the canine SF reduced cleavage of IGFBP-5 and IGF-1, resulting in reduced cartilage damage following anterior cruciate ligament rupture . These studies indicate that high C1s activity levels are likely detrimental to cartilage repair. Further, the complement cascade is known to be important in the pathogenesis of OA, with patients with OA demonstrating increased gene expression of complement agonists compared to inhibitors . OA-related pathogenesis, such as the release of cartilage extracellular matrix molecules and the production of inflammatory mediators, induces complement activation . The increased pre-operative levels that we identified in individuals who responded well to ACI perhaps indicate that ACI has potential to be successful in individuals who may have developed an early OA phenotype.
Analysis of the iTRAQ nLC-MS/MS and LF LC-MS/MS datasets, both independently and when combined, highlighted that there is a marked proteomic shift in response to cartilage harvest (i.e., between stages I and II of ACI). This analysis resulted in a plethora of candidate biomarkers that may have the potential to be informative regarding whether an individual is likely to respond well to ACI prior to chondrocytes being implanted during stage II. The proteoglycan, collagens II-, IX- and X-degrading enzyme, MMP3  has been biochemically validated as one of these candidate proteins that is significantly increased at stage II compared with stage I only in non-responders to ACI. Use of these biomarkers could have the potential to prevent the burden of a second surgery in a patient for whom this therapy is likely to be unsuccessful and could indicate that a greater period of time should be left from when the cartilage harvest procedure takes place to when the cells are implanted or that a tailored cartilage implantation procedure would be more efficacious.
To investigate the significant proteome shift that exists in response to cartilage harvest, pathway analyses were performed to better distinguish the underlying biological mechanisms that dictate whether an individual will respond to ACI. The acute-phase response was the pathway predicted to be most significantly differentially regulated in response to cartilage harvest in non-responders to ACI. In-depth assessment of individual protein changes within this pathway again highlighted the benefit of using independent proteomic techniques to profile the SF, because a large number of proteins were differentially abundant between stages I and II, only three of which were identified using both techniques. The acute-phase response is the body’s first systemic response to immunological stress, trauma and surgery . At the site of injury/trauma, pro-inflammatory cytokines are normally released, activating inflammatory cells and ultimately resulting in inflammatory mediators and cytokines being released into the extracellular fluid compartment to be circulated in the blood . Interestingly, previous bioinformatics analyses of the proteome of patients with late OA compared with healthy control subjects highlighted a dysregulated acute-phase response in the end-stage OA cohort . The exacerbated activation of the acute-phase response in non-responders following initial surgery could indicate that these patients have a greater immune response to surgery and that they have a lesser ability to dampen the acute-phase following surgery or that they have already developed an advanced OA phenotype, deeming a therapy to repair cartilage injury unsuitable.
Finally, the datasets of combined iTRAQ nLC-MS/MS and LF LC-MS/MS identified proteins were used to generate interactome models that represent the systemic proteomic response to cartilage harvest which exists within the SF of both ACI responders and non-responders, from which biological functional pathways could be further studied. Biological functional pathways that were identified using this approach, as well as using Ingenuity Pathway Analysis can most confidently be taken forward as candidates for further study because they have been identified by independent bioinformatic methods. Furthermore, given the complexity of the knee joint environment, it is likely that the responder/non-responder phenotype is the result of many subtle protein changes which together contribute to overall dysfunction of a biological network, rather than being the result of an individual biological molecule or pathway per se. Therefore, the interactome models generated in this study provide an important opportunity to consider how these proteins interact with one another and result in such phenotypes, and they also provide a platform for further studies to investigate how potential modifications to the ACI procedure (e.g., using co-incidental anti-inflammatory drugs in non-responders at stage II) may alter these biological networks. Thus, these models may provide a potential in silico tool for predicting ACI outcome, as is commonly used in drug development strategies .
This study highlights the advantage of using two independent proteomic techniques to profile a holistic overview of the SF proteome, ideal for unbiased identification of biomarker candidates. iTRAQ nLC-MS/MS analysis of SF samples from individuals who have either responded well or very poorly to ACI has highlighted proteins that, with further validation, have the potential to predict clinical outcome prior to treatment. We have confirmed that there is a marked SF proteome shift following cartilage injury, which is exacerbated in non-responders. Network and pathway analyses have demonstrated the complexity of the biological response underlying this proteome shift in non-responders, with several biological pathways identified that may act as targets for therapeutic intervention.
Autologous chondrocyte implantation
Body mass index
Complement C1S subcomponent
Enzyme-linked immunosorbent assay
Insulin-like growth factor 1
Insulin-like growth factor-binding protein 5
Isobaric tagging for relative and absolute quantitation
Lateral femoral condyle
Lateral tibial plateau
Medial femoral condyle
Matrix metalloproteinase 3
National Institute for Health and Care Excellence
Protein Interaction Network Analysis For Multiple Sets
Strong cation exchange chromatography
Proto-oncogene tyrosine-protein kinase protein
Kraus VB, Blanco FJ, Englund M, Henrotin Y, Lohmander LS, Losina E, et al. OARSI Clinical Trials Recommendations: soluble biomarker assessments in clinical trials in osteoarthritis. Osteoarthritis Cartilage. 2015;23(5):686–97.
National Institute for Health and Care Excellence (NICE). Autologous chondrocyte implantation for treating symptomatic articular cartilage defects of the knee. Technology appraisal guidance [TA477]. 4 Oct 2017. https://www.nice.org.uk/guidance/ta477/resources/autologous-chondrocyte-implantation-for-treating-symptomatic-articular-cartilage-defects-of-the-knee-pdf-82604971061701. Accessed 25 Nov 2017.
Hulme CH, Wilson EL, Peffers MJ, Roberts S, Simpson DM, Richardson JB, et al. Autologous chondrocyte implantation-derived synovial fluids display distinct responder and non-responder profiles. Arthritis Res Ther. 2017;19:150.
Gillogly SD, Voight M, Blackburn T. Treatment of articular cartilage defects of the knee with autologous chondrocyte implantation. J Orthop Sport Phys Ther. 1998;28(4):241–51.
Richardson JB, Caterson B, Evans EH, Ashton BA, Roberts S. Repair of human articular cartilage after implantation of autologous chondrocytes. J Bone Joint Surg Br. 1999;81(6):1064–8.
Wright KT, Mennan C, Fox H, Richardson JB, Banerjee R, Roberts S. Characterization of the cells in repair tissue following autologous chondrocyte implantation in mankind: a novel report of two cases. Regen Med. 2013;8:699–709.
Bhosale AM, Kuiper JH, Johnson WE, Harrison PE, Richardson JB. Midterm to long-term longitudinal outcome of autologous chondrocyte implantation in the knee joint: a multilevel analysis. Am J Sports Med. 2009;37(Suppl 1):131S–8S.
Wright KT, Kuiper JH, Richardson JB, Gallacher P, Roberts S. The absence of detectable ADAMTS-4 (aggrecanase-1) activity in synovial fluid is a predictive indicator of autologous chondrocyte implantation success. Am J Sports Med. 2017;45(8):1806–14.
Dugard MN, Kuiper JH, Parker J, Roberts S, Robinson E, Harrison P, et al. Development of a tool to predict outcome of autologous chondrocyte implantation. Cartilage. 2017;8:119–30.
Hsueh MF, Önnerfjord P, Kraus VB. Biomarkers and proteomic analysis of osteoarthritis. Matrix Biol. 2014;39:56–66.
Chiaradia E, Pepe M, Tartaglia M, Scoppetta F, D’Ambrosio C, Renzone G, et al. Gambling on putative biomarkers of osteoarthritis and osteochondrosis by equine synovial fluid proteomics. J Proteomics. 2012;75:4478–93.
Pan X, Huang L, Chen J, Dai Y, Chen X. Analysis of synovial fluid in knee joint of osteoarthritis: 5 proteome patterns of joint inflammation based on matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Int Orthop. 2012;36(1):57–64.
Sohn DH, Sokolove J, Sharpe O, Erhart JC, Chandra PE, Lahey LJ, et al. Plasma proteins present in osteoarthritic synovial fluid can stimulate production via Toll-like receptor 4. Arthritis Res Ther. 2012;14(1):R7.
Liao W, Li Z, Wang H, Wang J, Fu Y, Bai X. Proteomic analysis of synovial fluid: insight into the pathogenesis of knee osteoarthritis. Int Orthop. 2013;37(6):1045–53.
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(1):119–26.
Liao W, Li Z, Zhang H, Li J, Wang K, Yang Y. Proteomic analysis of synovial fluid as an analytical tool to detect candidate biomarkers for knee osteoarthritis. Int J Clin Exp Pathol. 2015;8(9):9975–89.
Mateos J, Lourido L, Fernández-Puente P, Calamia V, Fernández-López C, Oreiro N, Ruiz-Romero C, Blanco FJ. Differential protein profiling of synovial fluid from rheumatoid arthritis and osteoarthritis patients using LC–MALDI TOF/TOF. J Proteomics. 2012;75(10):2869–78.
Ritter SY, Subbaiah R, Bebek G, Crish J, Scanzello CR, Krastins B, et al. Proteomic analysis of synovial fluid from the osteoarthritic knee: comparison with transcriptome analyses of joint tissues. Arthritis Rheum. 2013;65(4):981–92.
Balakrishnan L, Nirujogi RS, Ahmad S, Bhattacharjee M, Manda SS, Renuse S, et al. Proteomic analysis of human osteoarthritis synovial fluid. Clin Proteomics. 2014;11(1):6.
Balakrishnan L, Bhattacharjee M, Ahmad S, Nirujogi RS, Renuse S, Subbannayya Y, et al. Differential proteomic analysis of synovial fluid from rheumatoid arthritis and osteoarthritis patients. Clin Proteomics. 2014;11(1):1.
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(10):4377–87.
Bhattacharjee M, Balakrishnan L, Renuse S, Advani J, Goel R, Sathe G, et al. Synovial fluid proteome in rheumatoid arthritis. Clin Proteomics. 2016;13:12.
Wu WW, Wang G, Baek SJ, Shen RF. Comparative study of three proteomic quantitative methods, DIGE, cICAT, and iTRAQ, using 2D gel- or LC-MALDI TOF/TOF. J Proteome Res. 2006;5(3):651–8.
Bantscheff M, Schirle M, Sweetman G, Rick J, Kuster B. Quantitative mass spectrometry in proteomics: a critical review. Anal Bioanal Chem. 2007;389(4):1017–31.
Lilley KS, Beynon RJ, Eyers CE, Hubbard SJ. Focus of quantitative proteomics. Proteomics. 2015;15(18):3101–3.
Herr MM, Fries KM, Upton LG, Edsberg LE. Potential biomarkers of temporomandibular joint disorders. J Oral Maxillofac Surg. 2011;69(1):41–7.
Roberts S, Evans H, Wright K, van Niekerk L, Caterson B, Richardson JB, et al. ADAMTS-4 activity in synovial fluid as a biomarker of inflammation and effusion. Osteoarthritis Cartilage. 2015;23(9):1622–6.
Kraus VB, Stabler TV, Kong SY, Varju G, McDaniel G. Measurement of synovial fluid volume using urea. Osteoarthritis Cartilage. 2007;15(10):1217–20.
Ehrich E, Davies G, Watson D, Bolohnese J, Seidenberg B, Bellamy N. Minimal perceptible clinical improvement with the Western Ontario and McMaster Universities Osteoarthritis Index questionnaire and global assessments in patients with osteoarthritis. J Rheumatol. 2000;27(11):2635–41.
Roos E, Lohmander L. The Knee Injury and Osteoarthritis Outcome Score (KOOS): from joint injury to osteoarthritis. Health Qual Life Outcomes. 2003;1:64.
Saris D, Vanlauwe J, Victor J, Almqvist K, Verdonk R, Bellemans J, et al. Treatment of symptomatic cartilage defects in the knee: characterized chondrocyte implantation results in better clinical outcome at 36 months in a randomized trial compared to microfracture. Am J Sports Med. 2009;37(Suppl 1):10S–9S.
Smith HJ, Richardson JB, Tennant A. Modification and validation of the Lysholm Knee Scale to assess articular cartilage damage. Osteoarthritis Cartilage. 2009;17:53–8.
Lysholm J, Gillquist J. Evaluation of knee ligament surgery results with special emphasis on use of a scoring scale. Am J Sports Med. 1982;10(3):150–4.
Stoscheck CM. Protein assay sensitive at nanogram levels. Anal Biochem. 1987;160(2):301–5.
Fuller HR, Mandefro B, Shirran SL, Gross AR, Kaus A, Botting CH, et al. Spinal muscular atrophy patient iPSC-derived motor neurons have reduced expression of proteins important in neuronal development. Front Cell Neurosci. 2016;9:506.
Fuller H, Slade R, Jovanov-Milošević N, Babić M, Sedmak G, Šimić G, et al. Stathmin is enriched in the developing corticospinal tract. Mol Cell Neurosci. 2015;69:12–21.
Oliveros JC. VENNY: an interactive tool for comparing lists with Venn diagrams. 2007. http://bioinfogp.cnb.csic.es/tools/venny/. Accessed 28 Sept 2017.
Albert R, Barabasi AL. Statistical mechanics of complex networks. Rev Mod Phys. 2002;74(1):47–97.
Vidal M, Cusick ME, Barabási AL. Interactome networks and human disease. Cell. 2011;144(6):986–98.
Cowley MJ, Pinese M, Kassahn KS, Waddell N, Pearson JV, Grimmond SM, et al. PINA v2.0: mining interactome modules. Nucleic Acids Res. 2012;40(Database issue):D862–5.
Barabási AL, Gulbahce N, Loscalzo J. Network medicine: a network-based approach to human disease. Nat Rev Genet. 2011;12(1):56–68.
Nepusz T, Yu H, Paccanaro A. Detecting overlapping protein complexes in protein-protein interaction networks. Nat Methods. 2012;9:471–2.
Vastrik I, D’Eustachio P, Schmidt E, Gopinath G, Croft D, de Bono B, et al. Reactome: a curated knowledge base of biologic pathways and processes. Genome Biol. 2008;8:R39.
Croft D, Mundo AF, Haw R, Milacic M, Weiser J, Wu G, et al. The Reactome pathway knowledgebase. Nucleic Acids Res. 2014;42(Database issue):D472–7.
Clemmons DR, Busby WH, Garmong A, Schultz DR, Howell DS, Altman RD, et al. Inhibition of insulin-like growth factor binding protein 5 proteolysis in articular cartilage and joint fluid results in enhanced concentrations of insulin-like growth factor 1 and is associated with improved osteoarthritis. Arthritis Rheum. 2002;46(3):694–703.
Busby WH, Yocum SA, Rowland M, Kellner D, Lazerwith S, Sverdrup F, et al. Complement 1s is the serine protease that cleaves IGFBP-5 in human osteoarthritic joint fluid. Osteoarthritis Cartilage. 2009;17(4):547–55.
van Osch GJVM, van den Berg WB, Hunziker EB, Hauselmann HJ. Differential effects of IGF-1 and TGFβ-2 on the assembly of proteoglycans in pericellular and territorial matrix by cultured bovine articular chondrocytes. Osteoarthritis Cartilage. 1998;6:187–95.
Wang Q, Rozelle AL, Lepus CM, Scanzello CR, Song JJ, Larsen DM, et al. Identification of a central role for complement in osteoarthritis. Nat Med. 2012;17(12):1674–9.
Murphy G. Matrix metalloproteinases and their inhibitors. Acta Orthop Scand Suppl. 1995;266:55–60.
Gruys E, Toussaint MJM, Niewold TA, Koopmans SJ. Acute phase reaction and acute phase proteins. J Zheihang Univ Sci B. 2005;6(11):1045–56.
Guney E, Menchem J, Vidal M, Barabasi AL. Network-based in silico drug efficacy screening. Nat Commun. 2016;7:10331.
We thank Professor Rob Beynon for his advice on study design and help with analysis of samples using label-free quantitation proteomic analysis.
We thank Arthritis Research UK for supporting this work via grants 19429, 20815 and 21122. The sponsors had no involvement in the study design, data collection and interpretation, or preparation of the manuscript. MJP is supported through a Wellcome Trust Clinical Intermediate Fellowship. This work was supported by Wellcome Trust grant 094476/Z/10/Z, which funded the purchase of the TripleTOF 5600 mass spectrometer at the Biomedical Sciences Research Complex Mass Spectrometry and Proteomics Facility, University of St. Andrews (Fife, UK).
Availability of data and materials
Proteomic data derived from this study have been deposited in the PRoteomics IDEntifications (PRIDE) ProteomeXchange and can be accessed using the identifier [PXD008321].
Ethics approval and consent to participate
SF samples from patients undergoing ACI were collected under three independent ethical approvals: Investigating the potential for cells and molecules isolated from orthopaedic patients for modelling and understanding pathogenic conditions and developing diagnostic markers and therapies for musculoskeletal disorders and spinal cord injury (National Research Ethics Service [NRES] committee North West – Liverpool East 11/NW/0875), Autologous cell therapy for Osteoarthritis: An evaluation of the safety and efficacy of autologous transplantation of articular chondrocytes and/or bone marrow derived stromal cells to repair chondral/osteochondral lesions of the knee (NRES committee West Midlands – Coventry and Warwickshire 11/WM/0175), and Arthritis and cartilage repair study (Shropshire and Staffordshire-Shropshire local research ethics committee 06/Q6201/9). All patients gave valid informed consent prior to their samples’ being collected.
Consent for publication
Not applicable for this study.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Canonical pathways altered in the synovial fluid of clinical nonresponders at stage I compared with stage II of ACI, identified using Ingenuity Pathway Analysis based on proteins that were identified using iTRAQ proteomics (≥ 1.2-fold change). Significance was assessed using a right-sided Fisher’s exact test; therefore, the most significant canonical pathways represent those that are the least likely to have been identified because of molecules being in the canonical pathway by random chance. The z-score represents canonical pathways that are likely activated or inhibited (based on the pattern of differentially abundant proteins); NaN means no prediction could be made based on the number of differentially abundant proteins in the pathway. Table S2. Canonical pathways altered in the synovial fluid of clinical responders at stage I compared with stage II of ACI, identified using Ingenuity Pathway Analysis, based on proteins which were identified using iTRAQ proteomics (≥ 1.2-fold change). Significance was assessed using a right-sided Fisher’s exact test; therefore, the most significant canonical pathways represent those that are the least likely to have been identified because of molecules being in the canonical pathway by random chance. The z-score represents canonical pathways that are likely activated or inhibited (based on the pattern of differentially abundant proteins); NaN means no prediction could be made based on the number of differentially abundant proteins in the pathway. (XLSX 19 kb)
About this article
- Autologous chondrocyte implantation (ACI)
- iTRAQ proteomics
- Label-free quantitation proteomics
- Synovial fluid
- Cartilage repair
- Complement C1S subcomponent
- Matrix metalloproteinase 3