Research | Open | Published:
A proteomic analysis of chondrogenic, osteogenic and tenogenic constructs from ageing mesenchymal stem cells
Stem Cell Research & Therapyvolume 7, Article number: 133 (2016)
Mesenchymal stem cells (MSCs) have prospective applications in regenerative medicine and tissue engineering but to what extent phenotype and differentiation capacity alter with ageing is uncertain. Consequently, any loss in functionality with age would have profound consequences for the maintenance of tissue viability and the quality of tissues. Proteomics enables the set of proteins responsible for a particular cell phenotype to be identified, as well as enabling insights into mechanisms responsible for age-related alterations in musculoskeletal tissues. Few proteomic studies have been undertaken regarding age-related effects on tissue engineered into cartilage and bone, and none for tendon. This study provides a proteome inventory for chondrogenic, osteogenic and tenogenic constructs synthesised from human MSCs, and elucidates proteomic alterations as a consequence of donor age.
Human bone-marrow derived MSCs from young (n = 4, 21.8 years ± 2.4SD) and old (n = 4, 65.5 years ± 8.3SD) donors were used to make chondrogenic, osteogenic and tenogenic tissue-engineered constructs. We utilised an analytical method relying on extracted peptide intensities as a label-free approach for peptide quantitation by liquid chromatography–mass spectrometry. Results were validated using western blotting.
We identified proteins that were differentially expressed with ageing; 128 proteins in chondrogenic constructs, 207 in tenogenic constructs and four in osteogenic constructs. Differentially regulated proteins were subjected to bioinformatic analysis to ascertain their molecular functions and the signalling pathways. For all construct types, age-affected proteins were involved in altered cell survival and death, and antioxidant and cytoskeletal changes. Energy and protein metabolism were the principle pathways affected in tenogenic constructs, whereas lipid metabolism was strongly affected in chondrogenic constructs and mitochondrial dysfunction in osteogenic constructs.
Our results imply that further work on MSC-based therapeutics for the older population needs to focus on oxidative stress protection. The differentially regulated proteome characterised by this study can potentially guide translational research specifically aimed at effective clinical interventions.
Mesenchymal stem cells (MSCs) are a heterogeneous population of multipotent cells with the ability to differentiate into cells of mesenchymal origin to accommodate tissue homeostasis and repair. These cells have prospective applications in clinical tissue regenerative strategies, due to their potential to differentiate into musculoskeletal cell lineages, including cartilage, bone and tendon [1, 2], high ex-vivo expansion capacity and relative ease of culture following isolation . The principles of tissue engineering involve a complex interplay of factors and to what extent MSC phenotype and differentiation capacity alter with ageing is unknown. Consequently, any loss in functionality with age would have profound consequences for the maintenance of tissue viability and the quality of tissues. The capacity of MSCs to differentiate into various types of tissue seems to change with age [4, 5]. For example, MSCs in older age are impaired with respect to osteogenic differentiation, preferring the adipogenic pathway of differentiation . However, many current research findings are conflicting. For instance, relating to chondrogenic differentiation of MSCs, one study identified a reduction in glycosaminoglycan with age  whereas another found no change . Furthermore, in osteogenic differentiation one group identified an increase in alkaline phosphatase with age  whilst another demonstrated a decrease . It is thought these discrepancies could be due to the heterogeneous population which is propagated within and amongst donor populations. One recent study demonstrated that commercially isolated MSCs had a shorter doubling time than freshly isolated cultures. This may be due to commercially isolated cells having a more homogeneous population, eliminating uncertainty that heterogeneous freshly isolated populations enact on MSC characterisation .
The regenerative capacity of MSCs is believed to be due to their high proliferation and differentiation capacities, paracrine action (such as the secretion of plasminogen activators and matrix metalloproteases; important for extracellular matrix (ECM) remodelling), and immunological privilege. However, therapeutic effectiveness of MSCs for most clinical applications remains limited, possibly due to the attenuation of their regenerative potential in aged patients with chronic diseases. Nevertheless, MSCs have already been used in clinical trials of cell therapy for cartilage repair and osteoarthritis (reviewed in ), bone fracture healing  and in a limited number of tendon therapies [13, 14].
Because of the interest in MSCs for regenerative medicine, and with an ever ageing population it is important to understand how the sets of proteins in MSC-derived musculoskeletal ‘tissue’ change with age, and what the potential effects are on clinical outcome. In a recent review, the basic protein inventory of a characteristic MSC has been assembled which encompasses six functional groups of proteins . Further work has identified the protein changes in osteogenic , myogenic , chondrogenic  and adipogenic differentiation  from bone-marrow-derived MSCs (BM-MSCs).
Proteomics enables the set of proteins responsible for a particular cell phenotype to be identified, as well as enabling insights into mechanisms responsible for age-related alterations in musculoskeletal tissues . Label-free quantification proteomics methods are based on the direct measurement of the mass spectrum (mass to charge ratio of ions in the gas phase) acquired signal. When constituent peptides are produced following protein digestion and are converted into ions, the most abundant proteins will produce the most ions and thus the greatest signal intensities . This type of workflow coupled with a high mass resolution mass spectrometer and bioinformatics enables hundreds of proteins to be identified and quantified in a single sample with high specificity and sensitivity .
Little information is available concerning the cellular and ECM protein changes that occur when BM-MSCs from different biologically aged donors are differentiated into chondrogenic or osteogenic constructs and no data are available for tenogenic constructs. This is important because tissue engineering aims to develop biomimetic tissues that recapitulate biological, structural and functional characteristics of native tissue. Age-related changes thus have potential implications for the tissue-engineering strategies used for enhancing musculoskeletal repair. This study aims to evaluate and compare the proteome of chondrogenic, osteogenic and tenogenic constructs derived from young and old human BM-MSCs in order to determine similar and distinct changes in the construct proteome with ageing.
All chemicals are supplied by Sigma unless stated otherwise. Human MSCs from young and old donors were purchased from Stem Cell Technologies (Grenoble, France) and Promocell (Heidelberg, Germany).
Cell culture and differentiation
Human bone-marrow-derived commercially available MSCs from young (n = 4, 21.8 years ± 2.4SD) and old (n = 4, 65.5 years ± 8.3SD) donors were grown to passage 4, and each donor was differentiated into chondrogenic, osteogenic  and tenogenic  constructs as described previously and used in all subsequent experiments . Engineered tendon constructs were formed in six-well plates coated with ~1.5 ml SYLGARD (WPI, Hertfordshire, UK) and pinned with staples. For each construct, 0.4 ml of 0.6 × 106 cells/ml was suspended in 80 μl of 20 mg/ml fibrinogen to which 8 μl of 200 U/ml thrombin was added. Then 160 μl of mixture was placed in each well to form an even covering. Each cell-embedded fibrin gel was cultured in 2 ml DMEM supplemented with 100 units/ml penicillin/streptomycin, 10 % FBS, 500 ng/ml amphotericin, 2 mM l-glutamine, 200 μM l-ascorbic acid 2-phosphate, non-essential amino acid at 10 μl/ml concentration and aprotinin at 10 μl/ml. Every other day the edges of the constructs were scored using a fine pipette tip and the media were changed. All constructs were fully contracted between the staples and harvested at 28 days post seeding . All tissue culture was undertaken in 5 % oxygen and media were changed every third day. MSCs from each donor were analysed separately in the three lineages and proteomic analysis was carried out on all preparations for each of eight donors.
Validation of differentiation
Differentiation was assessed for chrondrogenesis and osteogenesis by comparison with control MSCs treated identically except with maintenance media (complete Dulbecco’s Modified Eagle’s Media; Gibco) using histology and quantitative real-time PCR (qRT-PCR). Calcium depositions, a marker of bone formation and differentiation, were determined in osteogenic constructs using Alizarin red staining as described previously . Chondrogenic pellets were paraffin embedded and 4 μm sections taken and further stained with Alcian Blue/Nuclear Fast Red. Tendon constructs were fixed in 4 % paraformaldehyde, longitudinally embedded in paraffin and 4 μm sectioned on polylysine slides. Staining was undertaken with Masson’s Trichrome (collagen) .
Transmission electron microscopy (TEM) of tendon constructs was performed by fixation in 2.5 % glutaraldehyde in 0.1 M sodium cacodylate buffer for 8 hours, followed by buffer washing procedure and second fixation and contrast stain with 0.1 % osmium tetroxide for 90 minutes. Samples were stained with 8 % uranyl acetate in 0.69 % maleic acid for 90 minutes, dehydrated in ascending ethanol concentrations and embedded in epoxy resin. Ultrathin cross-sections (60–90 nm) were cut with a Reichert-Jung Ultracut on an ultramicrotome using a diamond knife. Cut cross-sections were then mounted on 200 mesh copper grids and stained with ‘Reynold’s Lead citrate’ stain for 4 minutes. Images were viewed in Philips EM208S Transmission Electron Microscope at 80 k.
Total RNA was prepared from constructs using 0.5 ml Tri Reagent (Ambion, Warrington, UK) per construct with homogenisation using a fine needle and syringe. The guanidinium thiocyanate–phenol–chloroform extraction technique was used as described previously . M-MLV reverse transcriptase and random hexamer oligonucleotides were used to synthesise cDNA from RNA (Promega, Southampton, UK) in a 25 ml reaction. Aliquots of 1 ml were amplified by PCR in 20 ml reaction volumes on an ABI 7300 Sequence Detector using a SYBR Green PCR mastermix (Applied Biosystems, Warrington, UK). cDNA was used for lineage-specific gene expression markers using qRT-PCR relative to GAPDH . Steady-state transcript abundance of potential endogenous control genes was measured in the RNAseq data (unpublished data). GAPDH was selected as the most stable endogenous control gene. Primers used are presented in Additional file 1: Table S1.
Protein extraction and sample preparation
Proteins were extracted from the constructs using either guanidine hydrochloride as described previously  (chondrogenic and osteogenic) or 0.1 % Rapigest™  (tenogenic), following optimisation of protein extraction methods for each construct type (data not shown). Protein extracts were normalised following protein assay using the Bradford assay with Coomassie Plus™ protein assay reagent (Thermo Scientific, Rockford, IL, USA) read at 660 nm. In-solution trypsin digestion was undertaken on all samples as described previously . Samples were desalted using C18 tips (Merck Millipore, Watford, UK) .
Construct soluble extracts of MSCs and constructs were analysed by one-dimensional sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) to assess gross quantitative/qualitative differences in protein profiles . Additionally, Rapigest™ extracts of mesenchymal stems cells prior to differentiation were also evaluated. Then 30 μg was loaded according to equal volumes after ethanol precipitation and resolubilisation in SDS loading buffer (Invitrogen) and stained with Coomassie.
Mass spectrometry and label-free quantification
Liquid chromatography tandem mass spectrometry (LC-MS/MS) was performed using a NanoAcquity™ ultraperformance LC (Waters, Manchester, UK) online to an LTQ-Orbitrap Velos mass spectrometer (Thermo-Fisher Scientific, Hemel Hempstead) as described previously . The proteomics data were deposited to the ProteomeXchange Consortium  via the PRIDE partner repository with the dataset identifier PXD001952.
For label-free quantification the Thermo raw files of the acquired spectra from in-solution tryptic digests were analysed by the ProgenesisQI™ software (Version 1; Waters, Manchester, UK) . Briefly, the top five spectra for each feature were exported from ProgenesisQI™ and utilised for peptide identification in PEAKS® 7 PTM (Bioinformatics Solutions Inc., Ontario, Canada) using the reviewed Uniprot human database. Search parameters used were: 10 ppm peptide mass tolerance and 0.6 Da fragment mass tolerance; one missed cleavage allowed; fixed modification, carbamidomethylation; and variable modifications, methionine, proline, lysine oxidation. Proteins were identified with a false discovery rate (FDR) of 1 % and a minimum of two peptides per protein. The resulting peptide-spectrum matches were imported into ProgenesisQI™ for label-free relative quantification. Differentially expressed (DE) proteins were defined with FDR p < 0.05 and ±2-fold regulation.
For neopeptide determination of tenogenic constructs, mass spectrometry data from the in-solution tryptic digests were analysed. Neopeptides were identified by searches against the reviewed Uniprot human database using Mascot with search parameters as already stated except: enzyme, semi-tryptic. Neopeptides were only included if they passed a conservatively applied threshold with e > 0.001. Patterns of fragmentation were determined for collagens, proteoglycans and glycoproteins.
Functional analysis of proteomics data
To determine gene ontology, functional analyses, networks, canonical pathways and protein–protein interactions of age-related DE proteins in each construct type we performed the analyses using the Protein ANalysis THrough Evolutionary Relationships (PANTHER) Classification System , the functional analysis and clustering tool from the Database for Annotation, Visualisation, and Integrated Discovery (DAVID bioinformatics resources 6.7)  and Ingenuity Pathway Analysis (IPA) . Proteins were further characterised and classified using MatrisomeDB .
Western blotting validation of mass spectrometry data
To confirm the mass spectrometry data, the relative abundance of cartilage oligomeric matrix protein (COMP) and biglycan in tenogenic and superoxide dismutase 1 (SOD1) in the chrondrogenic protein extracts was determined using automated western blotting (Wes Simple Western Analysis; ProteinSimple, San Jose, CA, USA) as described previously . All reagents were supplied by ProteinSimple. Simple western analysis was performed according to the user manual. Protocols were optimised for antibody and protein loading. In brief, protein extracts were mixed with a master mix to give a final concentration of 0.004–0.2 mg/ml total protein, 1× sample buffer, 1× fluorescent molecular weight markers and 40 mM dithiothreitol. Samples were heated at 95 °C for 5 min. Samples, blocking solution, primary antibodies, horseradish peroxidase-conjugated secondary antibodies, chemiluminescent substrate, and separation and stacking matrices were loaded into designated wells in a microplate. After plate loading, fully automated electrophoresis and immunodetection took place with the capillary system. Proteins were separated by molecular weight at 375 V for 25 minutes, and primary and secondary antibodies incubated for 30 minutes. All antibodies were diluted in antibody diluent to the required concentrations (Additional file 2: Table S2). Chemiluminescence was captured by a charge-coupled device camera, and the digital image was analysed using ProteinSimple Compass software. The relative amount of each protein, relative to total protein content, was calculated based on the peak area. A 29 kDa system control antibody was spiked into each sample to provide within-capillary normalisation. Statistical significance was determined using a Mann–Whitney test.
Statistical analysis was undertaken with Mann–Whitney U tests for qRT-PCR and neopeptide analysis using GraphPad Prism version 6.0 (GraphPad Software, San Diego, CA, USA).
Characterisation of tissue constructs
To confirm chondrogenic induction of MSCs, mature markers of chondrocytes were assessed; Alcian Blue staining for glycosaminoglycans and aggrecan, COL2A1 and SOX9 gene expression. In line with previous reports  we identified an increase in Alcian Blue staining and aggrecan, COL2A1 and SOX9 expression  (Fig. 1a–d), demonstrating chondrogenic differentiation of MSCs. Osteogenic differentiation was evaluated with Alizarin Red and RUNX2 gene expression. There was a significant increase in staining with Alizarin Red both visually and using quantitative analysis (Fig. 1e, f) and increased RUNX2 expression (Fig. 1g), demonstrating osteogenic differentiation of MSCs. Tenogenic differentiation was evaluated histologically following Masson’s Trichrome staining, indicating areas of organised and disorganised collagen fibril formation within the constructs. This was confirmed with TEM and with gene expression of COL1A1 (Fig. 1 h, i). There was no qualitative difference in the collagen organisation of tendon constructs derived from young and old MSCs. Further genes previously identified as markers of tendon expression  were significantly increased following tenogenic differentiation but were not age related; serpin peptidase inhibitor F (SERPINF1) (Fig. 1j) and thrombospondin 4 (THBS4) (Fig. 1k).
SDS-PAGE comparative analysis of protein extracts
One-dimensional SDS-PAGE of the soluble protein extracts demonstrated differences in the intensity of the staining and the number of bands between construct type and MSCs but not with age. Osteogenic profiles were most similar to undifferentiated MSCs, showing more cellular profiles. The chondrogenic profiles were the least complex, as demonstrated by few stained bands (Fig. 2).
A total of 2226 proteins for chondrogenic constructs, 2233 proteins for osteogenic constructs and 615 proteins for tenogenic constructs were identified within each construct type. Additional file 3: Table S3 presents detailed information on the identification of proteins in each. Each dataset was subjected to gene ontology using PANTHER which identified few differences in terms of gene ontological terms between each construct type (data not shown). Proteins in the ECM fraction were classified and further investigated using MatrisomeDB  (Additional file 3: Table S3). The overlap of proteins shared between construct types is shown Additional file 4: Figure S1. The percentage of matrisomal proteins compared with all proteins was 4 % for chondrogenic and osteogenic constructs and 11 % for tenogenic constructs.
Label-free relative quantification
To compare relative protein levels between young and old chondrogenic, osteogenic and tenogenic constructs, samples were processed for LC-MS/MS, and label-free quantitative analysis undertaken using ProgenesisQI™. Each young and old sample from each lineage was analysed separately and then comparisons were made between all young and old donors for each lineage. Principal component analysis of all of the proteins identified revealed that the proteins clustered according to the age of the donor with a principal component for chondrogenic constructs of 71 %, osteogenic constructs of 66 % and tenogenic constructs of 83 %. The number of DE proteins with at least a 2-fold change, p < 0.05 and q < 0.05 was 128 in chondrogenic constructs (28 higher in old, 100 lower), 207 in tenogenic constructs (201 higher in old, six lower) and four in osteogenic constructs (three higher in old – gamma-adducin (ADDG), uveal autoantigen with coiled-coil domains and ankyrin repeats (HOYNH8), and heterogeneous nuclear ribonucleoprotein (ROA0) – one lower – aldehyde dehydrogenase X(ALB1)) (Additional file 5: Table S4).
Gene ontology of DE proteins
Using DAVID analysis we identified ontology terms which were shared and distinct between construct types. For gene ontology analysis of the DE proteins for chondrogenic and tenogenic constructs we used the parameters 2-fold change, p < 0.05 and q < 0.05. For osteogenic analysis, because only four proteins were DE at 2-fold change with p < 0.05 and q < 0.05, we investigated gene ontology terms and undertook IPA analysis for the 63 DE proteins with 2-fold change and p < 0.05 (Additional file 6: Table S5). A significant term common to all constructs was actin cytoskeleton organisation. For chondrogenic constructs DAVID terms were also principally involved in tissue morphogenesis and cell adhesion, and for tenogenic constructs these were principally glycolysis and protein metabolism. For osteogenic constructs the single significant term ‘ribonucleoproteins’ was identified (Additional file 7: Table S6). A similar number of DE matrisomal proteins was identified in chondrogenic and tenogenic constructs but none in osteogenic constructs (Additional file 3: Table S3). Four matrisomal proteins were DE in both chondrogenic and tenogenic constructs with age; COL4A2, matrix metalloproteinase 14 (MMP14), matrix remodelling associated 5 (MXRA5) and thrombospondin 1(THBS1). In osteogenic constructs, four matrisomal proteins were DE: WNT1 inducible signalling pathway protein 2 (WISP2), collagen type XV alpha 1 (COL15A1), proline/arginine-rich end leucine-rich repeat protein (PRELP) and serpin peptidase inhibitor (SERPINE2) (Additional file 3: Table S3). Proteins DE for chondrogenic and tenogenic constructs (2-fold change, p < 0.05 and q < 0.05) and osteogenic constructs (2-fold change and p < 0.05) were further classified for gene ontology with PANTHER, demonstrating both similarities and differences in the type of proteins affected with age (Fig. 3).
Differential expressed genes and network analysis
For chondrogenic and tenogenic constructs the sets of DE proteins associated with ageing were analysed with IPA. The top canonical pathways, networks and diseases and biological functions with ageing are presented in Tables 1, 2 and 3 respectively. From the DE proteins associated with ageing input into IPA, a number of potential upstream regulators were identified in chondrogenic, tenogenic and osteogenic constructs (Additional file 8: Table S7). The activation z score is used to infer likely activation states of upstream regulators based on comparison with a model that assigns random regulation directions. Only upstream regulators with significant activation z scores were investigated to increase confidence in the data. Interesting upstream regulatory analysis revealed that targets of SMAD-2, SMAD-3, SMAD-4 and transforming growth factor beta (TGFβ) are regulated in chondrogenic constructs. Here TGFβ was predicted to be inhibited with ageing, with effects on tissue development and chondrogenic differentiation. In contrast, TGFβ was identified in tenogenic constructs as an activating upstream regulator in ageing. Targets of hypoxia-induced factor-1 alpha (HIF1α) were activated or inhibited in the same way as TGFβ in chondrogenic and tenogenic constructs. Targets of epidermal growth factor receptor family (ERBB4) were differentially regulated in osteogenic constructs (Additional file 8: Table S7). Interestingly, age-related changes in signalling pathways for chondrogenic constructs were principally in lipid metabolism (Fig. 4a), for tenogenic constructs were in glucose metabolism (Fig. 4b) and for osteogenic constructs were in mitochondrial dysfunction (Fig. 4c). Furthermore, age-related alterations in functions in tenogenic constructs included protein metabolism (Fig. 4d). Table 4 summarises the functional categories, canonical pathways and networks identified from the changes in protein expression in each construct with ageing using IPA.
Validation of mass spectrometry results by automated western blotting
Abundance of cartilage oligomeric matrix protein (COMP) and biglycan (tenogenic) and superoxide dismutase 1 (SOD1) was validated by western blotting (Fig. 5). Apparent molecular weights were larger than expected for COMP and biglycan. This is probably due to samples not being deglycosylated prior to western blotting, but this may also be due to evident differences in molecular weight (particularly for glycosylated proteins) using the ProteinSimple system for the majority of antibodies due to a difference in the separation matrix . In agreement with the mass spectrometry data, for tenogenic constructs COMP and biglycan were higher in abundance (p = 0.05) in old constructs. SOD1 was higher in abundance (p = 0.05) in the old chondrogenic constructs, also supporting the mass spectrometry results.
Identification of ECM fragmentation patterns in tendon constructs
Neopeptides were identified in all samples. Significantly more neopeptides were identified in the tendon constructs derived from old MSCs than from young MSCs (14 ± 2.5 versus 2.5 ± 1; p = 0.049). Together with significant pathways of 'activation of metabolism' and 'protein metabolism', this observation indicates increased turnover in tendon constructs derived from older MSC donors. A summary of the number of neopeptides identified in each condition is presented in Table 5. Very few neopeptides were identified for non-collagenous proteins whilst a large number of collagen proteins demonstrated neopeptides.
MSCs are an appealing source for cell-based treatment of musculoskeletal diseases and injury. Ageing is associated with various altered cellular phenotypes. Furthermore, the regeneration potential of MSCs is reduced with increasing age and is correlated with changes in cellular functions [5, 24]. This study of chondrogenic, osteogenic and tenogenic constructs derived from young and old MSCs provides a comprehensive proteomic analysis of tissue constructs whilst concurrently enabling improved understanding of the age-related functional and biological variations, which may affect their applications to regenerative medicine. In addition, our approach enables common and tissue-specific pathways of musculoskeletal ageing in an in-vitro system to be identified.
Characterisation of the tissue constructs was undertaken using standard methods following chondrogenic and osteogenic differentiation. For tenogenic constructs we used histological staining with Masson’s Trichrome and TEM to ascertain the presence and organisation of a collagenous matrix together with collagen type I, THBS4  and SERPINEF1 gene expression. The latter two have recently been identified as being the most DE genes in tendon differentiation but with low expression in chondrogenic differentiation . There was heterogeneity in the response to differentiation of these markers which was not age related, similar to that in aggrecan expression in chondrogenic constructs. The conditions of differentiation may impact on results. Because differentiation involves culture with different factors, which change according to the method used, it is conceivable that there might an age-related change in response to alteration of these factors which could affect results. Further work is required in this area.
The proteomic profiles of the constructs demonstrated using 1D gels showed that there was no gross difference in the profiles within construct type with ageing. When these results were compared with the number of identified proteins there were a higher number of protein identifications within the osteogenic (2233 proteins) and chondrogenic (2226 proteins) constructs, indicating a more complex proteome in these tissues. This was similar to the number of proteins identified in MSCs (2347 proteins) . A Rapigest™-based protein extraction workflow was utilised for tenogenic (615 proteins) constructs whereas a guanidine-based extraction was used for chondrogenic and osteogenic constructs, because Rapigest™ provided superior results compared with guanidine for tenogenic constructs in terms of proteins identified. The difference in number of protein identifications was probably due to the altered cell to ECM ratios within constructs, evident from histological sections. One limitation of the study is that it cannot be ruled out that some of the proteome pathway differences evident could be due to the culture format or duration.
We used label-free quantification to identify age-related DE proteins within each construct type. For our initial analysis we filtered data using p and q values (false-discovery adjusted p values for multiple testing) in order to reduce the number of false positives . This produced 128 DE proteins in chondrogenic constructs, 207 DE proteins in tenogenic constructs but only four proteins for osteogenic constructs. This indicates that at this level of filtering the age of the MSC donor has little effect on the proteome of osteogenic 2D constructs, whilst the protein composition of tendon constructs is most affected by MSC donor age. This disparity with the other two construct types could be due to their differentiation in 3D. We used standard 3D construct differentiation techniques for chondrogenic constructs due to problems with dedifferentiation into monolayers . The pellet characteristics closely mimic cartilage . In tenogenic differentiation, uniaxial tension in 3D is a requirement for differentiation . For osteogenic constructs we used the standard 2D system because it is the best described method and we could foresee problems with the mass spectrometry compatibility of materials included in many of the 3D osteogenic systems . However, 2D techniques inadequately produce the in-vivo environment for stem cells established by extrinsic and intrinsic cell signalling affecting biological function and differentiation capacity over time . Therefore, the effect of MSC donor age on the 3D osteogenic construct proteome should be studied in future.
For the musculoskeletal constructs it is essential to demonstrate the cellular phenotype and tissue composition, especially the ECM molecules that play a structural role and that contribute to the resulting mechanical properties. Therefore we identified the compositional and age-related DE matrisomal proteins in constructs. A number of matrisomal proteins were shared between all constructs such as COMP, TIMP1, decorin and biglycan, whilst some were shared between some types such as TIMP3 between chondrogenic and tenogenic constructs. This demonstrates that, similar to native tissues, the constructs have contrasting ECM profiles [20, 47]. Furthermore, when DE matrisomal proteins were investigated some proteins again shared age-related changes (COL4A2, MXRA5, THBS1 and MMP14) in chondrogenic and tenogenic constructs. Others such as plasminogen in chondrogenic constructs, and lumican in tenogenic constructs, were construct distinct. Interestingly, in agreement with results from 1D gels, the tenogenic constructs contained the most ECM matrisomal proteins as a percentage of all proteins identified within the construct type. Our findings revealed that the age of the donor MSCs had distinct or similar effects on construct ECM, depending on the differentiation lineage. The consequence of these altered matrices on the mechanical competence of the constructs requires further work.
Gene ontology revealed that metabolic processes were overrepresented in tenogenic constructs. This was also evident using IPA, which identified an increase in glucose and protein metabolism (both identified as activated in ageing), the latter related specifically to protein expression, proteolysis, catabolism and anabolism. Protein metabolism was demonstrated, for example, by an increase in abundance in old tendon constructs of TIMP-1, TIMP-3, MMP-2 and MMP-14. Furthermore, this was validated by neopeptide analysis, an indicator of protein turnover [20, 30]. Neopeptides represent ECM fragments produced by tissue remodelling. We have previously shown that neopeptide expression is altered in normal tendon ageing and disease . The changes could be attributed to altered metabolism within the cells present in ageing tendon constructs. However, the increased abundance of proteases could also be due to release of intracellular proteases because of cell death. The DNA content of old tendon constructs was reduced despite all constructs being seeded at an equivalent rate at the start of the experiments. This could be due to cell death in older constructs or reduced proliferation capacity. However, because there was a concomitant increase in DE ECM proteins, these protein metabolism changes seem to be due to a dysregulation of protein metabolism in ageing tendon constructs. These findings may also help understand how the tendon undergoes physiological remodelling that is evident in ageing.
Gene ontology also identified that DE age-related proteins in chondrogenic constructs were higher for the cellular component ECM and extracellular region proteins compared with the other construct types, indicating that donor age affects matrix proteins of chondrogenic constructs the most. These age-related changes in the ECM could have important implications for the quality of engineered tissue.
There were protein changes in lipid metabolism-related proteins in chondrogenic constructs. In cartilage, lipids are a source of energy and are incorporated into structural components and signalling molecules. Chondrocytes express several proteins for fatty acid metabolism and cholesterol biosynthesis and these molecules are increased during chondrogenesis . Chondrocyte lipid peroxidation has been suggested to have a role in cartilage ageing . We identified age-related changes in proteins involved with LXR activation and cholesterol biosynthesis. Given previous findings, it would seem that MSC donor age has an impact on lipid metabolism which could affect their chondrogenic potential further given that our culture conditions were hypoxic.
A further interesting feature derived from pathway analysis was the demonstration of an age-related inflammatory response in chondrogenic constructs, a significant feature in chondrogenic constructs. Ageing MSCs are known to undergo inflammageing and tissue-engineered cartilage may be more predisposed compared with other tissue types, similar to native ageing cartilage . Because MSC-derived chondrocytes are a potential treatment for chondral lesions , the use of allogeneic MSCs from younger donors may be beneficial in treating older patients. These findings demonstrate that our use of young and old donor-derived MSCs to produce musculoskeletal constructs could be a useful model to study musculoskeletal ageing.
A principal age-related feature of osteogenic constructs was mitochondrial dysfunction, when reactive oxygen species (ROS)-mediated oxygen stress overpowers the antioxidant defence system. Oxidative damage affects replication and transcription of mitochondrial DNA, leading to a decline in mitochondrial function and enhanced ROS production with further damage to mitochondrial DNA. In all constructs we demonstrated age-related protein changes involved in cell death and survival, and the alterations in oxidative stress probably contribute to this. Our results imply that MSC-derived tissue engineering from older donors must focus on oxidative stress protection.
Pathway analysis revealed that actin cytoskeleton changes were common to all ageing constructs. Others have identified an age-affected alteration in cytoskeletal organisation in rat MSCs . In all ageing constructs, similar to MSCs  we hypothesise that there is a decline in responsiveness to mechanical and biological signals due to a less dynamic cytoskeleton. The age-related network of proteins involved in cell migration and movement were also affected in all construct types. In total these findings are consistent with altered cytoskeletal dynamics affecting cell movement through coupling to actin organisation and turnover .
A number of upstream regulators of DE proteins were identified in each construct type. Our analysis revealed a number of significant transcriptional regulators potentially responsible for the protein changes. These data provide a starting point for future studies in MSC-derived tissue engineering of musculoskeletal constructs from older patients. One interesting finding was the contrasting roles of TGFβ and HIF1α in the DE proteins from chondrogenic and tenogenic constructs. TGFβ was significantly predicted to affect protein changes relating to tissue development in chondrogenic and tenogenic constructs and relating to differentiation in chondrogenic constructs, but in opposite directions; inhibited in chondrogenic but activated in tenogenic. Whilst there is only indirect evidence for a role of TGFβ in ageing, it is an important growth factor in development and differentiation. In our study, TGFβ signalling was predicted to be inhibited in chondrogenic constructs and activated in tenogenic constructs. This could be due to differing responses to oxidative stress previously identified here or distinctive requirements for culture conditions dependent on construct type. Furthermore, expression of the sets of age and tissue-specific transcriptional regulators may explain these findings, leading to similar molecular scenarios for some pathways (antioxidant, cell survival and cytoskeleton) but contrasting in others (e.g. protein and energy metabolism in tenogenic constructs and lipid metabolism in chondrogenic constructs).
Finally, the functions of the DE proteins identified which relate to cell death and survival, and antioxidant and cytoskeletal changes, are associated and important for chondrogenic, osteogenic and tenogenic differentiation. Chondrogenesis is characterised by changes in cell shape  and actin organisation is essential . In tenogenesis, cytoskeletal organisation is also paramount . Furthermore, osteogenesis is tightly regulated by ROS (reviewed in ). Age-related proteomic changes will thus affect the ability and quality of tissue-engineered constructs.
The study of musculoskeletal ageing in bone, cartilage and tendon is generally undertaken in isolation and it is often difficult to attain aged matched tissue samples in humans. We propose our approach as a model for musculoskeletal ageing that could be probed further to identify factors that may aid in recapitulation of a younger tissue phenotype. This is because as musculoskeletal tissues age they become more prone to age-related musculoskeletal disease such as osteoarthritis, tendinopathy and osteoporosis. We have identified some shared age-related characteristics (inflammageing, oxidative stress, cytoskeletal) which raise the prospects that common therapeutic targets could be developed to prevent these diseases. Understanding what drives these changes in diverse tissues could lead to the development of new therapeutic methods, which are advantageous to the musculoskeletal system in general.
Mass spectrometry-based proteomics provides an efficient method to monitor the complete profile of cellular and ECM molecules from tissue-engineered constructs. The differentially regulated proteome characterised by this study can potentially guide translational research specifically aimed at effective clinical interventions. Our data have provided valuable clues for our better understanding of the underlying mechanisms that are responsible for age-related changes in tissue-engineered constructs, thus assisting in the application of MSCs in cell-based therapy for cartilage, bone and tendon regeneration. These results also have significant implications for therapeutic cell source decisions (autologous or allogeneic), revealing the necessity of approaches to improve functionality of ageing MSCs.
Bone-marrow-derived mesenchymal stem cell
Cartilage oligomeric matrix protein
Database for Annotation Visualisation and Integrated Discovery
False discovery rate adjusted
hypoxia-induced factor-1 alpha
Ingenuity pathway analysis
Liquid chromatography tandem mass spectrometry
Mesenchymal stem cell
Matrix remodelling associated
Protein ANalysis THrough Evolutionary
Proline/arginine-rich end leucine-rich repeat protein
Qualitative real-time PCR
Runt-related transcription factor 2
Sodium dodecyl sulphate polyacrylamide gel electrophoresis
Serpin peptidase inhibitor
Transmission electron microscopy
Transforming growth factor beta
WNT inducible signalling pathway protein 2
Pittenger MF, et al. Multilineage potential of adult human mesenchymal stem cells. Science. 1999;284(5411):143–7.
Minguell JJ, Erices A, Conget P. Mesenchymal stem cells. Exp Biol Med (Maywood). 2001;226(6):507–20.
Kassem M, Abdallah BM. Human bone-marrow-derived mesenchymal stem cells: biological characteristics and potential role in therapy of degenerative diseases. Cell Tissue Res. 2008;331(1):157–63.
Sethe S, Scutt A, Stolzing A. Aging of mesenchymal stem cells. Ageing Res Rev. 2006;5(1):91–116.
Stolzing A, et al. Age-related changes in human bone marrow-derived mesenchymal stem cells: consequences for cell therapies. Mech Ageing Dev. 2008;129(3):163–73.
Laschober GT, et al. Age-specific changes of mesenchymal stem cells are paralleled by upregulation of CD106 expression as a response to an inflammatory environment. Rejuvenation Res. 2011;14(2):119–31.
Zaim M, et al. Donor age and long-term culture affect differentiation and proliferation of human bone marrow mesenchymal stem cells. Ann Hematol. 2012;91(8):1175–86.
Leskela HV, et al. Osteoblast recruitment from stem cells does not decrease by age at late adulthood. Biochem Biophys Res Commun. 2003;311(4):1008–13.
Alt EU, et al. Aging alters tissue resident mesenchymal stem cell properties. Stem Cell Res. 2012;8(2):215–25.
Okolicsanyi RK, et al. Human mesenchymal stem cells retain multilineage differentiation capacity including neural marker expression after extended in vitro expansion. PLoS One. 2015;10(9):e0137255.
Orth P, et al. Current perspectives in stem cell research for knee cartilage repair. Stem Cells Cloning. 2014;7:1–17.
Gomez-Barrena E, et al. Bone fracture healing: cell therapy in delayed unions and nonunions. Bone. 2015;70:93–101.
Mazzocca AD, et al. Bone marrow-derived mesenchymal stem cells obtained during arthroscopic rotator cuff repair surgery show potential for tendon cell differentiation after treatment with insulin. Arthroscopy. 2011;27(11):1459–71.
Docheva D, et al. Biologics for tendon repair. Adv Drug Deliv Rev. 2015;84:222–39.
Park HW, Shin JS, Kim CW. Proteome of mesenchymal stem cells. Proteomics. 2007;7(16):2881–94.
Sun HJ, et al. A proteomic analysis during serial subculture and osteogenic differentiation of human mesenchymal stem cell. J Orthop Res. 2006;24(11):2059–71.
Kurpinski K, et al. Proteomic profiling of mesenchymal stem cell responses to mechanical strain and TGF-beta1. Cell Mol Bioeng. 2009;2(4):606–14.
Ji YH, et al. Quantitative proteomics analysis of chondrogenic differentiation of C3H10T1/2 mesenchymal stem cells by iTRAQ labeling coupled with on-line two-dimensional LC/MS/MS. Mol Cell Proteomics. 2010;9(3):550–64.
Ju AHJ, Ko SP, Hyung SP. Membrane proteomic analysis of human mesenchymal stromal cells during adipogenesis. Proteomics. 2007;7(22):4181–91.
Peffers MJ, et al. Proteomic analysis reveals age-related changes in tendon matrix composition, with age- and injury-specific matrix fragmentation. J Biol Chem. 2014;289(37):25867–78.
Peffers MJ, Beynon RJ, Clegg PD. Absolute quantification of selected proteins in the human osteoarthritic secretome. Int J Mol Sci. 2013;14(10):20658–81.
Fabre B, et al. Label-free quantitative proteomics reveals the dynamics of proteasome complexes composition and stoichiometry in a wide range of human cell lines. J Proteome Res. 2014;13(6):3027–37.
Kapacee Z, et al. Synthesis of embryonic tendon-like tissue by human marrow stromal/mesenchymal stem cells requires a three-dimensional environment and transforming growth factor beta3. Matrix Biol. 2010;29(8):668–77.
Peffers MJ, et al. Age-related changes in mesenchymal stem cells identified using a multi-omics approach. Eur Cell Mater. 2016;31:136–59.
Paul H, Reginato AJ, Schumacher HR. Alizarin red S staining as a screening test to detect calcium compounds in synovial fluid. Arthritis Rheum. 1983;26(2):191–200.
Bancroft JD, Gamble M. Theory and practice of histological techniques. Philadelphia: Elsevier Health Sciences; 2008.
Chomczynski P, Sacchi N. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987;1987(162):156–9.
Peffers MJ, et al. Transcriptome analysis of ageing in uninjured human Achilles tendon. Arthritis Res Ther. 2015;17:33.
Little D, et al. Proteomic differences between male and female anterior cruciate ligament and patellar tendon. PLoS One. 2014;9(5):e96526.
Peffers MJ, Thornton DJ, Clegg PD. Characterization of neopeptides in equine articular cartilage degradation. J Orthop Res. 2015;34:106–120.
Vizcaino JA, et al. ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat Biotechnol. 2014;32(3):223–6.
Mi H, Muruganujan A, Thomas PD. PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Res. 2013;41(Database issue):D377–86.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
URL3. Ingenuity systems. Ingenuity pathway analysis. http://www.ingenuity.com/
Hynes RO, Naba A. Overview of the matrisome—an inventory of extracellular matrix constituents and functions. Cold Spring Harb Perspect Biol. 2012;4(1):a004903.
Chen JQ, Wakefield LM, Goldstein DJ. Capillary nano-immunoassays: advancing quantitative proteomics analysis, biomarker assessment, and molecular diagnostics. J Transl Med. 2015;13:182.
Oldershaw RA, et al. Notch signaling through Jagged-1 is necessary to initiate chondrogenesis in human bone marrow stromal cells but must be switched off to complete chondrogenesis. Stem Cells. 2008;26(3):666–74.
Mueller AM. A systems biology approach to musculoskeletal tissue engineering: transcriptomic and proteomic analysis of cartilage and tendon cells, in musculoskeletal biology. Liverpool: University of Liverpool; 2015. p. 634.
Peffers MJ. Personal communication. San Jose, CA, USA: ProteinSimple; 2016.
Jelinsky SA, et al. Tendon-selective genes identified from rat and human musculoskeletal tissues. J Orthop Res. 2010;28(3):289–97.
Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003;100(16):9440–5.
von der Mark K, et al. Relationship between cell shape and type of collagen synthesised as chondrocytes lose their cartilage phenotype in culture. Nature. 1977;267(5611):531–2.
Bosnakovski D, et al. Chondrogenic differentiation of bovine bone marrow mesenchymal stem cells in pellet cultural system. Exp Hematol. 2004;32(5):502–9.
Kuo CK, Tuan RS. Mechanoactive tenogenic differentiation of human mesenchymal stem cells. Tissue Eng Part A. 2008;14(10):1615–27.
Cordonnier T, et al. 3D environment on human mesenchymal stem cells differentiation for bone tissue engineering. J Mater Sci Mater Med. 2010;21(3):981–7.
Park E, Patel AN. Changes in the expression pattern of mesenchymal and pluripotent markers in human adipose-derived stem cells. Cell Biol Int. 2010;34(10):979–84.
Peffers M. Protein and transcriptomic sigantures of cartilage ageing and disease, in Musculoskeletal Biology. Liverpool: University of Liverpool; 2013. p. 338.
Bernstein P, et al. Expression pattern differences between osteoarthritic chondrocytes and mesenchymal stem cells during chondrogenic differentiation. Osteoarthr Cartil. 2010;18(12):1596–607.
Tiku ML, Shah R, Allison GT. Evidence linking chondrocyte lipid peroxidation to cartilage matrix protein degradation. Possible role in cartilage aging and the pathogenesis of osteoarthritis. J Biol Chem. 2000;275(26):20069–76.
Greene MA, Loeser RF. Aging-related inflammation in osteoarthritis. Osteoarthr Cartil. 2015;23(11):1966–71.
Kasper G, et al. Insights into mesenchymal stem cell aging: involvement of antioxidant defense and actin cytoskeleton. Stem Cells. 2009;27(6):1288–97.
Webb DJ, Parsons JT, Horwitz AF. Adhesion assembly, disassembly and turnover in migrating cells—over and over and over again. Nat Cell Biol. 2002;4(4):E97–100.
Jin EJ, et al. Akt signaling regulates actin organization via modulation of MMP-2 activity during chondrogenesis of chick wing limb bud mesenchymal cells. J Cell Biochem. 2007;102(1):252–61.
Xu B, et al. RhoA/ROCK, cytoskeletal dynamics, and focal adhesion kinase are required for mechanical stretch-induced tenogenic differentiation of human mesenchymal stem cells. J Cell Physiol. 2012;227(6):2722–9.
Atashi F, Modarressi A, Pepper MS. The role of reactive oxygen species in mesenchymal stem cell adipogenic and osteogenic differentiation: a review. Stem Cells Dev. 2015;24(10):1150–63.
The authors wish to thank the Wellcome Trust and the Medical Research Council (MRC) and Arthritis Research UK as part of the MRC-Arthritis Research UK Centre for Integrated research into Musculoskeletal Ageing (CIMA) for funding.
MJP was supported by a Wellcome Trust Veterinary Integrated Research Fellowship (WT088557MA) and supported by the Medical Research Council (MRC) and Arthritis Research UK as part of the MRC-Arthritis Research UK Centre for Integrated research into Musculoskeletal Ageing (CIMA).
Availability of supporting data
All proteomics data were deposited at the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD001952.
MJP was responsible for conception and design, collection and assembly of data, data analysis and interpretation, manuscript writing and final approval of the manuscript. JC was responsible for collection and assembly of data, manuscript writing and final approval of the manuscript. JL was responsible for provision of study material, conception and design, manuscript writing and final approval of the manuscript. CP was responsible for conception and design, manuscript writing and final approval of the manuscript. PDC was responsible for conception and design, manuscript writing and final approval of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Ethical approval was not required because all cells were purchased.
Presenting the primer sequences. (XLSX 11 kb)
Presenting details of the antibodies used for validation studies. (XLSX 10 kb)
Presenting the proteins identified in PEAKS PTM using the human reviewed database. (XLSX 425 kb)
Showing a Venn diagram of the ECM proteins identified in constructs using MatrisomeDB. (DOCX 100 kb)
Presenting the classification of proteins in chondrogenic, osteogenic and tenogenic constructs using MatrisomeDB. (XLSX 23 kb)
Presenting results of the DE proteins identified in constructs made from young and old MSCs using label-free quantification with ProgenesisQI™. (XLSX 52 kb)
Complete list of significantly expressed proteins and DAVID analysis of the expression patterns. (XLSX 297 kb)
Presenting upstream log ratio molecule predicted activation p values of target molecules in the dataset. (XLSX 25 kb)