- Research article
- Open Access
Matrix assisted laser desorption ionization mass spectrometry imaging identifies markers of ageing and osteoarthritic cartilage
Arthritis Research & Therapy volume 16, Article number: R110 (2014)
Cartilage protein distribution and the changes that occur in cartilage ageing and disease are essential in understanding the process of cartilage ageing and age related diseases such as osteoarthritis. The aim of this study was to investigate the peptide profiles in ageing and osteoarthritic (OA) cartilage sections using matrix assisted laser desorption ionization mass spectrometry imaging (MALDI-MSI).
The distribution of proteins in young, old and OA equine cartilage was compared following tryptic digestion of cartilage slices and MALDI-MSI undertaken with a MALDI SYNAPT™ HDMS system. Protein identification was undertaken using database searches following multivariate analysis. Peptide intensity differences between young, ageing and OA cartilage were imaged with Biomap software. Analysis of aggrecanase specific cleavage patterns of a crude cartilage proteoglycan extract were used to validate some of the differences in peptide intensity identified. Immunohistochemistry studies validated the differences in protein abundance.
Young, old and OA equine cartilage was discriminated based on their peptide signature using discriminant analysis. Proteins including aggrecan core protein, fibromodulin, and cartilage oligomeric matrix protein were identified and localised. Fibronectin peptides displayed a stronger intensity in OA cartilage. Age-specific protein markers for collectin-43 and cartilage oligomeric matrix protein were identified. In addition potential fibromodulin and biglycan peptides targeted for degradation in OA were detected.
MALDI-MSI provided a novel platform to study cartilage ageing and disease enabling age and disease specific peptides in cartilage to be elucidated and spatially resolved.
Osteoarthritis (OA) is an age related joint disease characterized by alterations in the chondrocytes and loss of cartilage extracellular matrix (ECM) . Age is the most common risk factor for its initiation and progression, with symptomatic OA affecting 10 to 20% of people more than 50 years old . ‘Wear and tear’ injuries due to mechanical loading over many years, genetic polymorphisms and other factors contribute to the development of the disease.
Although much work has been undertaken to investigate the pathogenesis of OA the molecular mechanisms involved are not fully understood, with few validated markers for disease diagnosis and progression being available. Treatment of OA is limited to pain relief with few disease-modifying therapeutics available to patients. Hence, methodologies which permit the detailed study of aging and OA cartilage organization could improve the knowledge of how aging increases the risk of OA. Previously the comparison of young, old and OA tissue has been challenging as the provision of normal old and old OA is difficult and few markers have been identified that differentiate between age or disease-related change in cartilage except urinary collagen type II .
Measurement tools enabling the capture of the dynamic and complex interplay between proteins, lipids, DNA and other molecules are necessary. This interplay involves the appearance, interaction and disappearance of many species on varying time scales. Capturing this spatial and temporal information will help in the quest to understand the pathogenesis of OA and ultimately provide disease modifying treatments. Mass spectrometry (MS) is an analytical tool that enables the accurate measurement of the mass-to-charge ratio of ions. Matrix assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) of tissue samples is a powerful technique that allows for spatially resolved, comprehensive and specific characterization of hundreds of unknown molecular species (proteins, peptides, lipids or metabolites) in situ in a single molecular imaging experiment . The label-free nature of the technique combined with high spatial resolving power better than 50 micrometer makes it an ideal tool to study the molecular aspects of aging in cartilage tissue . Specific protein patterns revealed by MSI have been demonstrated as prognostic  and diagnostic indicators . MSI techniques have only recently been implemented in cartilage research. One study employed a time of flight secondary ion mass spectrometry (TOF-SIMS) workflow in order to acquire molecular-specific spatial distribution of lipids in normal and OA cartilage . A further study used MALDI-MSI to identify and study OA specific peptides and proteins in cartilage . In this and other works, a direct correlation between the detected peptides and the distribution and identity of the original precursor protein was achieved .
The horse provides an excellent model for the study of cartilage in man as it suffers clinical joint diseases similar to man  and, as such, has been used as a model of naturally occurring OA  due to the extensive knowledge of its pathogenesis and clinical experience of the disease . Indeed, the incidence of equine metacarpophalangeal OA in young racehorses  in training is similar to the incidence of post traumatic OA in man . Additionally, the articular cartilage thickness is also comparable between the species .
The aim of this study was to establish peptide profiles in young, aging and OA horse cartilage of full thickness cartilage from surface to bone, with a high spatial distribution to determine changing molecular events distinct between aging and disease. We hypothesize that many alterations in specific peptide abundance are due to proteinase driven degradation.
Sample collection, preparation and processing
Ethical approval was not required for this study as samples were obtained from an abattoir. Full thickness equine cartilage slices were removed from the mid condyle region of metacarpophalangeal joints collected from an abattoir and snap frozen in liquid nitrogen within hours. All samples were scored macroscopically using Kawcak scoring for pathological grading of the distal condyles of metacarpophalageal III . For age-related studies, samples were taken from skeletally mature young (four years old, n = 3) and old (greater than fifteen years old, n = 3) horses. For OA studies, old normal samples were taken from macroscopically and microscopically normal condylar areas adjacent to diseased tissue (n = 3) (greater than 15 years old). For young horses, one year is equivalent to an age of approximately 3.5 years of a human . Hence, horses more than 15 years old used in this study equate to humans of older than 52 years. For OA studies, skeletally mature donors >15 years were chosen with mild macroscopic OA changes.
Cartilage was sectioned perpendicular to the joint surface at 12 μm thicknesses on a cryostat Microm HM 525 (Microm International, Walldorf, Germany) and thaw mounted on glass slides. All samples were dried in a vacuum desiccator for 10 minutes prior to further processing.
Tissue digestion, matrix deposition and MALDI-MSI
Young and old cartilage was studied in duplicate while old and OA cartilage was examined in triplicate and then compared by MALDI-MSI. The sections were washed in 70% ethanol and chloroform for 30 seconds each. Tissue digestion and matrix deposition and MALDI-MSI were undertaken as previously described .
Multivariate analysis and data interpretation following MALDI-MSI
Data pre-treatment, principal component analysis (PCA) and discriminant analysis (DA) were performed using AMOLF in-house build MATLAB (The MathWorks, Natick, MA, USA) software tools . The first 20 principal components (PC) were used as input for discriminant analysis (DA). PCA was performed on the raw spectral data and was used as a data compression and noise filtering step before application of DA. In short, PCA is an unsupervised statistical method that aims at pooling together a maximum amount of variance in a minimum number of independent variables. These new variables (principal components) are a construct of the original variables (mass channels) weighed by their correspondence with the new variable. This weighing factor is referred to as the PC loading. DA is a supervised multivariate method that calculates the combination of variables (in this case PCs) which gives maximum separation between the pre-defined groups. DA maximizes the between group - within group variance.
Biomap 22.214.171.124 software (Novartis Pharma AG, Basel, Switzerland) was used to generate ion images. Normalization of the intensity of all m/z channels was performed using the intensity of the m/z 190 matrix peak. P values for statistical differences found in MALDI-MSI experiments were calculated with one-way analysis of variance (ANOVA) with a Bonferroni post hoc test using SPSS (SPSS Inc. Chicago IL, USA) following normality testing. Differences were considered to be statistically significant at P ≤0.05. The data were expressed as mean intensity ± standard error.
Tissue digestion and matrix deposition for profiling experiments
Protein identification was undertaken with profiling experiments directly from each donor tissue by applying 10 μL of trypsin diluted in 25 mM ammonium bicarbonate 0.05 μg/μL directly and incubating overnight at 37°C. α-Cyano-4-hydroxycinnamic acid (CHCA) matrix was applied as previously described . Data dependent analysis (DDA) of tryptic peptides was performed with the MALDI SYNAPT™ HDMS system. Every MS survey scan was followed by collisional fragmentation of the most intense ions with subsequent collection of tandem mass spectroscopy (MS/MS) spectra. Direct MS/MS fragmentations were performed directly from the tissue on the peptides that differentiate young, old or OA tissue. These target peptides were found following DA. In addition, fragmentation was also undertaken on a peptide list generated from bibliographical analysis of previous MALDI-MSI studies in cartilage . The resulting data files were submitted to an in-house Mascot server (Matrix Sciences, London, UK) and searched against the Unihorse and Swissprot databases. Search parameters used were: peptide mass tolerances 50 ppm, fragment mass tolerance of 0.5 Da, 1+ ions, missed cleavages; 3, and instrument type MALDI-Q-TOF. Modifications used were variable oxidation of methionine.
Histological staining and immunohistochemistry studies
Following MSI the tissue sections were washed in 70% ethanol for five minutes to remove matrix. The slides were immersed in Harris hematoxylin solution (Sigma-Aldrich, Dorset, UK) for one minute. After washing with water for 15 minutes, they were rinsed in 95% ethanol and counterstained in an eosin solution for 30 seconds. Digital images were acquired with the Mirax system (Carl Zeiss, Sliedrecht, The Netherlands) after dehydrating steps. A score based on the modified Mankin’s scoring system for semi-quantitative histological assessment of all equine cartilage samples was undertaken .
Macroscopic grading was used to allocate samples into normal or OA groups, and scores of greater than 0 were assigned OA. Only macroscopically normal cartilage was used for age related studies. All samples for age related studies had a modified Mankin’s grade of 0 (data not shown). The macroscopic grading and Mankin’s score for each sample used in the old normal versus OA studies are shown in Additional file 1: Table S1.
To perform immunohistochemistry studies, paraffin-embedded tissues were cut into 4 μm sections by a microtome. Protease antigen retrieval was undertaken with a bacterial protease solution. Sections were incubated with polyclonal rabbit anti-human fibronectin (Dako, Heverlee, Belgium) at 1:2400 dilution, overnight at 4°C. Goat anti-rabbit biotinylated secondary antibody at 1:100 was incubated at room temperature for 30 minutes. Immunostaining was detected with avidine-biotin complex (Vector Laboratories, Peterborough, UK) and chromogenic substrate diaminobenzidine (Sigma-Aldrich, Dorest, UK). Counterstaining was undertaken with Papanicolau stain (Merck Millipore, Watford, UK). Negative control experiments were carried out with omission of the primary antibody and substitution with non-immune rabbit immunoglobulin G (IgG) (Abcam, Cambridge, UK). No staining was observed in negative control experiments. Digital images were visualized and acquired with a Nikon DS-L2 stand-alone control unit using a Nikon eclipse 80i microscope. Analysis was undertaken with ImageJ. ImageJ 1.42 image software and band densitometry was analyzed by arbitrarily drawing a rectangular box around the selected protein bands of interest. The same rectangular box was used to measure all bands of interest and band intensities were quantified by conversion to profile plot histograms. By specifically selecting the peak of interest and closing off the peak with the line tool, background-subtracted density of each peak of interest is quantified. All densitometry data are normalized to the band intensity of the loading control used.
Proteoglycan-enriched fraction isolation from cartilage
In order to validate our hypothesis that the reduction in specific peptides identified by MALDI-MSI was due to degradation of those specific peptides we undertook proteinase digestion of proteoglycan-enriched cartilage. Equine articular cartilage from the grossly normal metacarpophalangeal joint of three horses, mean age 10 ± 1 years, obtained from an abattoir was pooled and pulverized with a dismembranator (Miko, S-Braun, Bethlehem, PA, USA). Proteins were extracted with cartilage extraction buffer containing 4 M guanidinium chloride, 50 mM sodium acetate and proteinase inhibitors (complete protease inhibitors, ethylenediaminetetraacetic acid (EDTA)-free, Roche, Lewes, UK), pH 6.0 using end-over-end mixing for 20 hours at 4°C. After extraction the soluble fraction (proteoglycan-enriched) was removed following centrifugation for 15 minutes at 13,000 g at 4°C. Following dialysis in a 14,000-kD cut-off membrane (Spectrapor, Breda, NL) for 24 hours at 4°C against 0.1 M sodium acetate, pH 6.0 in the presence of proteinase inhibitors, the extract was clarified using centrifugation for 15 minutes at 13,000 g at 4°C. A crude proteoglycan-enriched extract was isolated by associative cesium chloride (CsCl) density gradient centrifugation , also in the presence of the proteinase inhibitors. The supernatant was fractionated in an associative CsCl density gradient (starting density 1.5 g/ml) for 60 hours at 100,000 g in an ultracentrifuge (Beckman 50Ti, Galway, Ireland). The tube was fractionated into quarters, A1-A4 , and the A1-A2 fraction dialyzed first against 0.1 M sodium acetate for 48 hours at 4°C and then against ultrapure water for 36 hours at 4°C, both in the presence of proteinase inhibitors. The samples were then lyophilized.
Proteinase digestion of the proteoglycan-enriched fraction in vitro
A reconstituted A1-A2 aliquot was digested in proteinase digestion buffer (50 mM Tris HCl, 100 mM NaCl, 10 mM CaCl2, pH 7.5) with 0.014 nmol truncated human recombinant ADAMTS-4 (a distintegrin and metalloproteinase with thrombospondin motifs 4) (Calbiochem, La Jolla, CA, USA) for seven hours at 37°C. A control aliquot was incubated under the same conditions in the presence of recombinant protein formulation buffer. The enzymatic digestion reactions (1 ml) were stopped by addition of EDTA.
In-solution tryptic digestion
The crude proteoglycan extracts following from the proteinase digestion were easily trypsin digested as previously described . Samples were desalted and purified using C18 resin in the form of a ZipTip® (Merck Millipore, Rockland, MA, USA).
LC-MS/MS analysis of proteinase digested cartilage proteoglycan-enriched fraction
Liquid chromatography (LC)-MS/MS analysis was performed using NanoAcquity™ Ultraperformance LC (Waters, Manchester, UK) on line to an LTQ-Orbitrap Velos (Thermo-Fisher Scientific, Hemel Hempstead, UK) as described previously .
Data analysis proteinase digested cartilage proteoglycan-enriched fraction
In order to identify potential cleavage sites of ADAMTS-4 identification raw spectra were converted to mascot generated files (mgf) using Proteome Discoverer software (Thermo, Hemel Hempstead, UK). The resulting mgf files were searched against the Unihorse database using an in-house Mascot  server (Matrix Science, London, UK). Search parameters used were: enzyme; none, peptide mass tolerances 10 ppm, fragment mass tolerance of 0.6 Da, 1+, 2+ and 3+ ions, missed cleavages; 1, and instrument type ESI-TRAP. Modifications included were: fixed; carbamidomethyl cysteine and variable; oxidation of methionine.
Patterns of fragmentation were determined for biglycan and fibromodulin. The probability that a match was correct (P <0.05) was determined using the Mascot derived ion score where P was the probability that the observed match was a random event.
To validate the results further, raw data files were loaded into PEAKS® Studio 6.0 (Bioinformatics Solutions Inc., Waterloo, Canada) and de novo sequencing and protein identification performed. PEAKS® software employs multiple analytical algorithms and is able to identify every peptide in the data, and can validate database searches using de novo sequencing results. Estimate false discovery rate (FDR) function was used in order to create a ‘decoy fusion’ database based on the Unihorse database. Data prerefinement was undertaken by choosing peak centroiding, charge deconvolution, and deisotope options. The quality value was set greater than 0.65. The search parameters for the PEAKS® software were identical to those for the Mascot search. Results generated using PEAKS® Studio were manually curated against the MASCOT search engine results. Only peptides with a 10lgP score of >20 were considered as significant .
MS/MS profiling experiments on cartilage slices were undertaken in order to determine proteins that could be identified directly from the tissue (Table 1). Tryptic peptides of cartilage oligomeric matrix protein (COMP) and fibromodulin were identified in all groups. Following PCA and DA between young, old and OA groups, young and old, and old and OA, the peaks with the highest absolute loadings in the DF1 spectra were targeted directly from tissue slices for MS fragmentation and database searching in order to identify the protein it pertained to. In addition, a list of masses was also selected that had been previously identified in human cartilage MALDI-MSI studies (Table 1) .
MALDI-MSI reveals the peptide signature in young, old and OA cartilage
Young, old and OA cartilage
After imaging experiments, data from young, old and OA cartilage groups were analyzed together. The spectra of each group revealed different profiles. A combined spectrum of representative digested equine young, old and OA samples is demonstrated in Additional file 2: Figure S1. PCA analysis and a scatter plot of DF1-DF2 revealed differences between each sample type (Additional file 3: Figure S2A and B). In addition, a comparison between biological repeats of the same subgroup was undertaken to assess the reproducibility of the score distribution of independent samples (data not shown).
Following DA a number of peaks were detected as specific for young, old or OA groups (Figure 1A). A peptide specific to young samples with m/z 2414.9 was subsequently identified as collectin-43 using fragmentation targeting of this mass (Figure 1B).
In addition, comparisons were also made using PCA and DA between young and old, and old and OA groups separately.
Young and old cartilage
Differential peptide distributions between young and old cartilage were assessed by MALDI-MSI followed by DA. The resulting DF1 classified the data into young and old (Figure 1C). The DF1 spectra (Figure 1D) identified that the peaks of the positive part of the DF1 were specific from young cartilage samples, such as collectin-43, and the peaks of the negative part were more abundant in cartilage from old donors.
The peak distribution intensity differences in young, old and OA cartilage of the identified peptides in Table 1 were analyzed by Biomap software. While there was no difference in the intensity of the majority of the peptides between groups, including collagen type II, a significantly different intensity of peptides for COMP, m/z 2256.1 and biglycan 2027.2 was observed (Figure 2).
Old and OA cartilage
The data were assigned to two groups (OA and old) and DA was performed (Figure 3). The positive part of DF1 showed the masses specific to OA, the negative part the masses specific to the old group.
Of these, peptide m/z 1366.5 was found to be imaged with Biomap at greater intensities in all OA samples (mean intensity 1.01 ± 0.03), P <0.001, young or old versus OA). Mean intensity for young was 0.05 ± 0.003 and old 0.06 ± 0.04 (Additional file 4: Figure S3). In addition, we checked other peptides as found in other studies on human samples . Distinctive intensity differences were observed for the previously identified fibromodulin peptide ELHLDHNQISR, m/z 1361.7. This peptide is homologous to the horse and was significantly reduced in OA equine cartilage (Figure 4A and B). Other peptides, such as m/z 1349.7, identified previously as a fibronectin peptide and homologous to the horse sequence were visualized. The highest intensity differences were between young (mean intensity 0.11 ± 0.009 and OA (mean intensity 0.36 ± 0.12) samples (P = 0.018, young versus OA). However, the intensity between old (0.12 ± 0.007) and OA was also significantly different (Figure 4C and D). Interestingly, a peptide m/z 1401.7 identified as fibronectin and previously nominated as a human OA marker  was visualized in young, old and OA samples (Figure 4E and F). The highest intensity differences for this peptide were between young (mean intensity 0.09 ± 0.009) and OA (mean intensity 0.17 ± 0.11) samples (P = 0.02, young versus OA). The higher intensity of fibronectin peaks in OA cartilage identified by MSI was orthogonally validated using fibronectin immunohistochemical staining (Figure 5). There was an increase in fibronectin staining in OA cartilage compared to old cartilage (mean pixel intensity 29.5 ± 1.9 old versus 41.5 ± 1.9 OA, n = 3, P <0.05). Primary and secondary antibody controls indicated that the primary antibody was specific to the antigen and the label was specific to the primary antibody (data not shown).
Finally, we measured the intensity of the distribution of peptides in different areas of the cartilage in young, old and OA cartilage in order to assess if the distributions were homogenous. However, we did not find differences in peptide distribution.
Validation of fibromodulin and biglycan fragmentation patterns
In order to validate our MALDI-MSI data and test our hypothesis that alterations in specific peptide abundance identified by Biomap were due to proteinase driven degradation, we identified ADAMTS-4 derived fragmentation patterns in fibromodulin and biglycan. A number of ADAMTS-4 derived neopeptides were identified for fibromodulin and biglycan which may account for the reduction in peptide intensity identified by Biomap (Table 2).
MALDI-MSI technology has been previously used in biomarker discovery  and to localize proteins, peptides and lipids in normal and OA human cartilage [8, 9]. A recent study used MALDI-MSI to investigate synovial tissue from OA and rheumatoid arthritis patients resulting in the identification of differentially expressed proteins between the two diseases . While these studies were undertaken on cartilage from normal and late OA adult human knees, nothing has been described until now in relation to age-related changes in cartilage and MALDI-MSI. In addition, this is the first time potential biomarkers of age and OA have been identified in cartilage using this technology.
Following MALDI-MSI and DA the peptide profile of young, old and OA horse cartilages could be distinguished. There was also a marked difference between the peptides identified in old versus OA as well as young versus old groups. Interestingly, when DA was applied to all three groups together (young, old, OA) there was a contribution of the old samples to young and OA indicating that peptides within old samples were also present in young and OA samples. This is not surprising as the majority of the proteins in young, old and OA cartilage will be the same ECM proteins . Furthermore, previous studies have shown that tissue classified as OA still contains areas or cells with a healthy signature [8, 9]. When a peak profile comparison was made between old and OA a distinction between OA related peptides was made which suggests that some age-related changes in peptides occur which do not contribute to disease related changes.
The major findings in this study were the identification of potential degradation sites in OA cartilage in biglycan and fibromodulin together with a number of potential markers of both age-related and disease-associated changes in cartilage. In addition, potential degradation sites within specific COMP peptides were identified. COMP is a major protein with a role in cartilage structural integrity through its collagen I, II, IX and chondrocyte binding capacity, interaction with other ECM proteins including matrilin-3  and role in fibrillogenesis . It has been proposed as a biomarker for arthritis . Interestingly, the peptides identified were from the collagen binding C-terminal of COMP which others have demonstrated is important in intra- and extracellular processes . Indeed, mutations in the genes encoding COMP and matrilin-3 result in multiple epiphyseal dysplasias [31, 32]. Measurement of intact COMP and fragments thereof in synovial fluid or serum correlates to cartilage destruction in rheumatoid arthritis (RA) and OA patient studies . Interestingly, studies in rats have found that the plasma levels of COMP are age dependent . We observed a decrease in both young and OA compared to old cartilage in the distribution of the COMP peptide FYEGELVADSNVVLDTTMR. This peptide is located in the C-terminal end of COMP which binds collagen I, IX and II, and regulates fibril formation . Our previous studies following cytokine stimulation of mature equine cartilage explants identified a neopeptide which indicated degradation of this peptide at Asn712–713Thr (data not shown). This suggests that in OA there is degradation within this tryptic peptide resulting in its reduced expression, demonstrated when OA cartilage was imaged. In young cartilage it may represent reduced synthesis or cartilage remodelling. The peptide represents a possible marker of age, but not disease, related changes in cartilage ECM. Future studies will investigate this finding further with the use of specifically designed monoclonal antibodies. It should be noted that MS intensities do not always reflect analyte concentration. It is possible that alterations in the nature of the tissue in aging and disease could influence ionization efficiency and, hence, cause ion suppression effects . However, as our previous studies identified a degradation product within this peptide, this is less likely.
In addition, a marker of young cartilage, peak 2414.9 identified as collectin-43, was demonstrated. Collectin-43 is a C-type serum lectin with collagenous regions and a member of the collectin family of soluble proteins that are effector molecules in innate immunity . Biglycan and decorin have been identified as binding collectin-43 and may have an important role in the resolution of C1q-mediated inflammatory processes in cartilage. Biglycan and decorin may down-regulate proinflammatory effects mediated by the collectins .
The study also demonstrated a reduction in some peptides for biglycan and fibromodulin, in both young and old compared to OA cartilage. Biglycan and fibromodulin are members of the small leucine rich repeat proteoglycan family with important collagen binding properties [38, 39]. Structural changes related to aging are evident in biglycan, where there appears to be a cleavage in the amino terminal domain resulting in a ‘no-glycan’ biglycan as the terminal peptide containing the glycosaminoglycan chain separates from the protein core . In order to study these findings further a cartilage-derived proteoglycan-enriched fraction was treated with the recombinant aggrecanase ADAMTS-4 and subjected to trypsin digestion and LC-MS/MS. In OA, the degradation of cartilage is characterized by a loss of ECM caused by secreted proteases, principally matrix metalloproteases and aggrecanases (reviewed ). We hypothesized that a reduction in the intensity of the fibromodulin tryptic peptide ELHLDHNQISR m/z 1361.7 and the biglycan mid-region tryptic peptide NHLVEIPPNLPSSLVELR m/z 2027.2 in OA samples was due to proteinase driven fragmentation within the peptide. The reduction in intensity of these peptides either represents a reduction in the synthesis, degradation from disease or differences in ionization. However, as another peptide identified in biglycan IQAIELEDLLR m/z 1312.7 (also mid-region) does not reveal reduced intensity in OA it is more likely that it is lost following degradation. Thus in OA, it is likely that there is cleavage within these tryptic peptides resulting in a loss of this mass and so a reduction in peak intensity.
Fibromodulin peptides were identified and their distribution imaged. The peptides m/z 1955.2 and 1361.7 correspond to the fibromodulin tryptic peptides identified in our human OA MALDI-MSI study . Discrepancies in the abundance of the peptides may be due to the stage of OA in the two different studies. In the human study, the OA samples represented late stage OA while here we used lower grade early OA samples. One of a number of peptides, identified in MALDI-MSI studies of human cartilage , and with the same mass in the horse investigated was the fibromodulin peptide ELHLDHNQISR m/z 1361.7. This was significantly reduced in OA samples despite the distribution of other fibromodulin peptides identified being unchanged, for example with m/z 1557.0 or 1955.2 identified in profiling experiments between categories by either DA or following peak intensity analysis and statistical testing (data not shown). This, together with the identification of an ADAMTS-4 driven cleavage site between Asn186–187Gln would again indicate degradation of this peptide as opposed to a reduction in synthesis of fibromodulin. ADAMTS-4 is a pertinent enzyme in the pathogenesis of OA and although fibromodulin has been previously identified as a substrate for ADAMTS-4 , this was at the Tyr 44–45Ala bond . Degraded fragments of the core fibromodulin protein have been observed in OA cartilage  and with age . Removal of this portion of fibromodulin would result in weaker interactions of collagen fibers with surrounding structures. This illustrates the potential usefulness of MALDI-MSI in identifying and spatially resolving novel cleavage sites with pathological relevance. Indeed, using two different methodological approaches, this study indicates that cleavage of the fibromodulin peptide ELHLDHNQISR and the biglycan peptide NHLVEIPPNLPSSLVELR is disease, and not age, related. Insights such as this may aid in the understanding of the age-related, but not age-distinct, disease OA.
A number of tentative OA markers were detected. The intensity of the differences in abundance of peptides with m/z 1366.5 and 1349.7 in young, old and OA samples found their abundance increased in OA as detected by MALDI-MSI imaging. The m/z 1349.7 had been identified from human studies  as being derived from fibronectin and its sequence homolog confirmed between the horse and man. Furthermore, we identified the fibronectin peptide m/z 1401.7 as being more abundant in OA cartilage compared to young but not old cartilage. In a human OA MALDI-MSI study this peptide was elevated in OA . Here, the fibronectin peptide m/z 1401.7 appears to be affected by age-related changes. Fibronectin, an ECM glycoprotein, and fibronectin fragments have been associated with enhanced levels of catabolic cytokines and up regulation of MMPs involved in both normal homeostasis and arthritic diseases. Fibronectin fragments and fibronectin-aggrecan complexes have previously been suggested as biomarkers of OA [46, 47]. Therefore, these peptides may provide promising biomarkers of OA as they are not affected by age-related changes. Fibronectin may provide a key species for potential diagnostic and drug targets. The increase in some fibronectin peptides identified in OA cartilage is most likely to be due to an increase in synthesis, according to other authors .
The protein melanoma inhibitory activity protein 3 (MIA3) was also identified for the first time in MALDI-MSI studies of cartilage. This is a secreted protein expressed by chondrocytes with a fundamental role in maintaining the chondrocyte phenotype . A similar distribution of a peptide with the same mass as MIA3, m/z 1186.74 in young, old and OA tissue was evident (data not shown) which was not surprising as even mature chondrocytes secrete MIA3 .
This study undertook analysis on equine cartilage from the metacarpophalangeal joint which is a high motion joint similar to the knee in man. However, cartilage from different joints and different species may represent a different set of matrix changes; hence, it was not surprising that there were some key differences between this study and a previous MALDI-MSI study in end stage OA human cartilage . Although previous studies enabled differential peptide resolution between superficial and deep layers of cartilage to be identified, this was not possible using equine cartilage from the metacarpal joint. This was due to the thickness of the equine cartilage from the metacarpophalangeal III bone which measures 400 to 450 μm, considerably thinner than the human cartilage used. The cartilage from the stifle joint of the horse varies in thickness from 1,760 to 2,215 μm  and may provide an alternative source of cartilage for studies attempting to resolve changes in superficial and deep layers of equine cartilage. Secondly, there was no evidence for heterogeneity of peptide distribution that had been previously identified in OA cartilage. This was probably because the former study used severe OA cartilage while here only mild histological changes were evident.
The ex vivo imaging of aged and diseased cartilage provided ‘label-free’ and stain-free information about its biomolecular composition. We identified a number of potential degradation sites in some key ECM proteins as well as markers of cartilage aging. With the resolving power of MALDI-MSI certain to increase in the future and improved methods to identify peptides in situ, further MALDI-MSI studies using a greater age range may help in understanding why aging cartilage is more prone to OA.
a distintegrin and metalloproteinase with thrombospondin motifs
cartilage oligomeric matrix protein
data dependant analysis
matrix assisted desorption ionization
melanoma inhibitory activity protein
tandem mass spectrometry
mass spectrometry imaging
principal component analysis
time of flight secondary ion mass spectrometry
- 95% CI:
95% confidence interval.
Goldring MB: The role of the chondrocyte in osteoarthritis. Arthritis Rheum. 2000, 43: 1916-1926.
Lawrence RC, Felson DT, Helmick CG, Arnold LM, Choi H, Deyo RA, Gabriel S, Hirsch R, Hochberg MC, Hunder GG, Jordan JM, Katz JN, Kremers HM, Wolfe F, National Arthritis Data Workgroup: Estimates of the prevalence of arthritis and other rheumatic conditions in the United States. Part II. Arthritis Rheum. 2008, 58: 26-35.
Mouritzen U, Christgau S, Lehmann HJ, Tanko LB, Christiansen C: Cartilage turnover assessed with a newly developed assay measuring collagen type II degradation products: influence of age, sex, menopause, hormone replacement therapy, and body mass index. Ann Rheum Dis. 2003, 62: 332-336.
Seeley EH, Caprioli RM: Molecular imaging of proteins in tissues by mass spectrometry. Proc Natl Acad Sci U S A. 2008, 105: 18126-18131.
Amstalden van Hove ER, Smith DF, Heeren RM: A concise review of mass spectrometry imaging. J Chromatogr A. 2010, 1217: 3946-3954.
Balluff B, Rauser S, Meding S, Elsner M, Schone C, Feuchtinger A, Schuhmacher C, Novotny A, Jutting U, Maccarrone G, Sarioglu H, Ueffing M, Braselmann H, Zitzelsberger H, Schmid RM, Höfler H, Ebert MP, Walch A: MALDI imaging identifies prognostic seven-protein signature of novel tissue markers in intestinal-type gastric cancer. Am J Pathol. 2011, 179: 2720-2729.
Gustafsson JO, Oehler MK, Ruszkiewicz A, McColl SR, Hoffmann P: MALDI Imaging Mass Spectrometry (MALDI-IMS)-application of spatial proteomics for ovarian cancer classification and diagnosis. Int J Mol Sci. 2011, 12: 773-794.
Cillero-Pastor B, Eijkel G, Kiss A, Blanco Garcia FJ, Heeren RM: A time of flight secondary ion mass spectrometry based molecular distribution distinguishes healthy and osteoarthritic human cartilage. Anal Chem. 2012, 84: 8909-8916.
Cillero-Pastor B, Eijkel GB, Kiss A, Blanco FJ, Heeren RM: Matrix-assisted laser desorption ionization-imaging mass spectrometry: a new methodology to study human osteoarthritic cartilage. Arthritis Rheum. 2013, 65: 710-720.
Cillero-Pastor B, Heeren RM: Matrix assisted laser desorption ionisation mass spectrometry imaging for peptide and protein analyses: a critical review of on-tiisue digestion. J Proteom Res. 2013, 13: 325-335.
Innes JF, Clegg P: Comparative rheumatology: what can be learnt from naturally occurring musculoskeletal disorders in domestic animals?. Rheumatology (Oxford). 2010, 49: 1030-1039.
McIlwraith CW, Frisbie DD, Kawcak CE: The horse as a model of naturally occurring osteoarthritis. Bone Joint Res. 2012, 1: 297-309.
Goodrich LR, Nixon AJ: Medical treatment of osteoarthritis in the horse - a review. Vet J. 2006, 171: 51-69.
Neundorf RH, Lowerison MB, Cruz AM, Thomason JJ, McEwen BJ, Hurtig MB: Determination of the prevalence and severity of metacarpophalangeal joint osteoarthritis in Thoroughbred racehorses via quantitative macroscopic evaluation. Am J Vet Res. 2010, 71: 1284-1293.
Dirschl DR, Marsh JL, Buckwalter JA, Gelberman R, Olson SA, Brown TD, Llinias A: Articular fractures. J Am Acad Orthop Surg. 2004, 12: 416-423.
Brommer H, Laasanen MS, Brama PA, van Weeren PR, Helminen HJ, Jurvelin JS: Functional consequences of cartilage degeneration in the equine metacarpophalangeal joint: quantitative assessment of cartilage stiffness. Equine Vet J. 2005, 37: 462-467.
Kawcak CE, Frisbie DD, Werpy NM, Park RD, McIlwraith CW: Effects of exercise vs experimental osteoarthritis on imaging outcomes. Osteoarthritis Cartilage. 2008, 16: 1519-1525.
Human-horse age chart. http://equineresources.com/pfizer-animal-health/259-pfizer-horse- human-age-chart ,
McIlwraith CW, Frisbie DD, Kawcak CE, Fuller CJ, Hurtig M, Cruz A: The OARSI histopathology initiative - recommendations for histological assessments of osteoarthritis in the horse. Osteoarthritis Cartilage. 2010, 18: S93-S105.
Lark MW, Bayne EK, Flanagan J, Harper CF, Hoerrner LA, Hutchinson NI, Singer II, Donatelli SA, Weidner JR, Williams HR, Mumford RA, Lohmander LS: Aggrecan degradation in human cartilage. Evidence for both matrix metalloproteinase and aggrecanase activity in normal, osteoarthritic, and rheumatoid joints. J Clin Invest. 1997, 100: 93-106.
Heinegard D: Polydispersity of cartilage proteoglycans. Structural variations with size and buoyant density of the molecules. J Biol Chem. 1977, 252: 1980-1989.
Peffers MJ, Beynon RJ, Clegg PD: Absolute quantification of selected proteins in the human osteoarthritic secretome. Int J Mol Sci. 2013, 14: 20658-20681.
Yan J, Liu L, Li X, Wang F, Zhu T, Yuan P, Zhang Q: A potential use of collagen-hyaluronan-chondroitin sulfate tri-copolymer scaffold for cartilage tissue engineering. Zhongguo Xiu Fu Chong Jian Wai Ke Za Zhi. 2006, 20: 130-133. In Chinese
Reed CC, Iozzo RV: The role of decorin in collagen fibrillogenesis and skin homeostasis. Glycoconj J. 2002, 19: 249-255.
Minerva L, Ceulemans A, Baggerman G, Arckens L: MALDI MS imaging as a tool for biomarker discovery: methodological challenges in a clinical setting. Proteomics Clin Appl. 2012, 6: 581-595.
Kriegsmann M, Seeley EH, Schwarting A, Kriegsmann J, Otto M, Thabe H, Dierkes B, Biehl C, Sack U, Wellmann A, Kahaly GJ, Schwamborn K, Caprioli RM: MALDI MS imaging as a powerful tool for investigating synovial tissue. Scand J Rheumatol. 2012, 41: 305-309.
Peffers MP: Proteomic and Transcriptomic Signatures in Cartilage Ageing and Disease. 2013, Liverpool: The University of Liverpool
Mann HH, Ozbek S, Engel J, Paulsson M, Wagener R: Interactions between the cartilage oligomeric matrix protein and matrilins. Implications for matrix assembly and the pathogenesis of chondrodysplasias. J Biol Chem. 2004, 279: 25294-25298.
Halasz K, Kassner A, Morgelin M, Heinegard D: COMP acts as a catalyst in collagen fibrillogenesis. J Biol Chem. 2007, 282: 31166-31173.
Tseng S, Reddi AH, Di Cesare PE: Cartilage oligomeric matrix protein (COMP): a biomarker of arthritis. Biomark Insights. 2009, 4: 33-44.
Briggs MD, Chapman KL: Pseudoachondroplasia and multiple epiphyseal dysplasia: mutation review, molecular interactions, and genotype to phenotype correlations. Hum Mutat. 2002, 19: 465-478.
Chapman KL, Mortier GR, Chapman K, Loughlin J, Grant ME, Briggs MD: Mutations in the region encoding the von Willebrand factor A domain of matrilin-3 are associated with multiple epiphyseal dysplasia. Nat Genet. 2001, 28: 393-396.
Saxne T, Heinegard D: Cartilage oligomeric matrix protein: a novel marker of cartilage turnover detectable in synovial fluid and blood. Br J Rheumatol. 1992, 31: 583-591.
Vingsbo-Lundberg C, Saxne T, Olsson H, Holmdahl R: Increased serum levels of cartilage oligomeric matrix protein in chronic erosive arthritis in rats. Arthritis Rheum. 1998, 41: 544-550.
Rohner TC, Staab D, Stoeckli M: MALDI mass spectrometric imaging of biological tissue sections. Mech Ageing Dev. 2005, 126: 177-185.
Holmskov U, Malhotra R, Sim RB, Jensenius JC: Collectins: collagenous C-type lectins of the innate immune defense system. Immunol Today. 1994, 15: 67-74.
Groeneveld TW, Oroszlan M, Owens RT, Faber-Krol MC, Bakker AC, Arlaud GJ, McQuillan DJ, Kishore U, Daha MR, Roos A: Interactions of the extracellular matrix proteoglycans decorin and biglycan with C1q and collectins. J Immunol. 2005, 175: 4715-4723.
Wiberg C, Klatt AR, Wagener R, Paulsson M, Bateman JF, Heinegard D, Morgelin M: Complexes of matrilin-1 and biglycan or decorin connect collagen VI microfibrils to both collagen II and aggrecan. J Biol Chem. 2003, 278: 37698-37704.
Hedbom E, Heinegard D: Interaction of a 59-kDa connective tissue matrix protein with collagen I and collagen II. J Biol Chem. 1989, 264: 6898-6905.
Roughley PJ, Nguyen Q, Mort JS, Hughes CE, Caterson B: Proteolytic degradation in human articular cartilage: its relationship to stromelysin. Agents Actions Suppl. 1993, 39: 149-159.
Bondeson J, Wainwright S, Hughes C, Caterson B: The regulation of the ADAMTS4 and ADAMTS5 aggrecanases in osteoarthritis: a review. Clin Exp Rheumatol. 2008, 26: 139-145.
Gendron C, Kashiwagi M, Lim NH, Enghild JJ, Thogersen IB, Hughes C, Caterson B, Nagase H: Proteolytic activities of human ADAMTS-5: comparative studies with ADAMTS-4. J Biol Chem. 2007, 282: 18294-18306.
Fushimi K, Troeberg L, Nakamura H, Lim NH, Nagase H: Functional differences of the catalytic and non-catalytic domains in human ADAMTS-4 and ADAMTS-5 in aggrecanolytic activity. J Biol Chem. 2008, 283: 6706-6716.
Roughley PJ, White RJ, Cs-Szabo G, Mort JS: Changes with age in the structure of fibromodulin in human articular cartilage. Osteoarthritis Cartilage. 1996, 4: 153-161.
Cs-Szabo G, Roughley PJ, Plaas AH, Glant TT: Large and small proteoglycans of osteoarthritic and rheumatoid articular cartilage. Arthritis Rheum. 1995, 38: 660-668.
Zack MD, Arner EC, Anglin CP, Alston JT, Malfait AM, Tortorella MD: Identification of fibronectin neoepitopes present in human osteoarthritic cartilage. Arthritis Rheum. 2006, 54: 2912-2922.
Scuderi GJ, Golish SR, Cook FF, Cuellar JM, Bowser RP, Hanna LS: Identification of a novel fibronectin-aggrecan complex in the synovial fluid of knees with painful meniscal injury. J Bone Joint Surg Am. 2011, 93: 336-340.
Bosserhoff AK, Kondo S, Moser M, Dietz UH, Copeland NG, Gilbert DJ, Jenkins NA, Buettner R, Sandell LJ: Mouse CD-RAP/MIA gene: structure, chromosomal localization, and expression in cartilage and chondrosarcoma. Dev Dyn. 1997, 208: 516-525.
Lougheed JC, Holton JM, Alber T, Bazan JF, Handel TM: Structure of melanoma inhibitory activity protein, a member of a recently identified family of secreted proteins. Proc Natl Acad Sci U S A. 2001, 98: 5515-5520.
Cillero-Pastor B, Eijkel B, Kiss A, Blanco RJ, Heeren RMA: Matrix assisted laser desorption ionization imaging mass spectrometry: A new technique to study human osteoarthritic cartilage. ArthritisRhem. 2013, 65: 710-720.
Frisbie DD, Cross MW, McIlwraith CW: A comparative study of articular cartilage thickness in the stifle of animal species used in human pre-clinical studies compared to articular cartilage thickness in the human knee. Vet Comp Orthop Traumatol. 2006, 19: 142-146.
Mandy Peffers was supported by a Wellcome Veterinary Integrated Research Fellowship WT088557MA. This work was funded by an EU Farm Animal Proteomics COST initiative (COST Action FA1002). In addition we acknowledge Foundation of the Fundamental Research on Matter- Netherlands Organisation for Scientific Research (FOM-NWO) for their assistance. This publication was supported by the Dutch national program COMMIT.
The authors declare they have no competing interests.
MJP, BCP, PDC and RMAH conceived and designed the study. MJP and BCP undertook the acquisition and analysis of data. GBE contributed to the analysis and interpretation of the data. All authors contributed to the drafting of the article. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Microscopic (A) and macroscopic (B) evaluation of normal and OA cartilage samples. Tables containing gross and microscopic scores for the normal old and OA cartilage. (PDF 199 KB)
Additional file 2: Combined spectrum of representative digested equine young, old and OA samples reveal different profiles. Representative spectrum from each group. (PDF 129 KB)
Additional file 3: A) PCA scatter plot B. Discriminant function plot. PCA scatter plot of the first 20 principal components. B) DF1 and DF2 scores. (PDF 277 KB)
Additional file 4: m/z 1366.5. Biomap image and histogram for a m/z 1366.5 peptide in representative samples of young, old and OA samples. (PDF 162 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Peffers, M.J., Cillero-Pastor, B., Eijkel, G.B. et al. Matrix assisted laser desorption ionization mass spectrometry imaging identifies markers of ageing and osteoarthritic cartilage. Arthritis Res Ther 16, R110 (2014). https://doi.org/10.1186/ar4560
- Cartilage Oligomeric Matrix Protein
- Cartilage Oligomeric Matrix Protein
- Laser Desorption Ionization Mass Spectrometry
- Multiple Epiphyseal Dysplasia
- Matrix Assisted Laser Desorption Ionization Mass Spectrometry Imaging