Heterogeneity of proteome dynamics between connective tissue phases of adult tendon

  1. Howard Choi
  2. Deborah Simpson
  3. Ding Wang
  4. Mark Prescott
  5. Andrew A Pitsillides
  6. Jayesh Dudhia
  7. Peter D Clegg
  8. Peipei Ping
  9. Chavaunne T Thorpe  Is a corresponding author
  1. Department of Physiology and Medicine, David Geffen School of Medicine, UCLA, United States
  2. Centre for Proteome Research, Biosciences Building, Institute of Integrative Biology, University of Liverpool, United Kingdom
  3. Department of Comparative Biomedical Sciences, Royal Veterinary College, United Kingdom
  4. Department of Clinical Sciences and Services, Royal Veterinary College, United Kingdom
  5. Department of Musculoskeletal Biology, Institute of Ageing and Chronic Disease, University of Liverpool, United Kingdom
8 figures, 1 table and 2 additional files

Figures

Metabolic labelling of rats using deuterium.

(A) Schematic showing 2H2O labelling of rats and details of sample collection, and interfascicular matrix isolation by laser capture microdissection (scale bar = 100 µm). (B) 2H enrichment in serum occurred rapidly, reaching a plateau of 5.6% by day 4 and remained constant throughout the study. Data are shown as mean ± SD and raw data provided in Figure 1—source data 1. (C) Example extracted ion chromatograms from tendon samples demonstrating increased abundance of higher mass isotopomer peaks over the course of the study (M1 – M5). (D) Example curve showing the relative abundance of the unlabelled monoisotopic peak (M0) for the cartilage oligomeric matrix protein (COMP) peptide QMEQTYWQANPFR, calculated by ProTurn software and plotted as a function of time.

Figure 1—source data 1

Enrichment of serum water as measured by GC-MS.

https://cdn.elifesciences.org/articles/55262/elife-55262-fig1-data1-v1.xlsx
Protein-protein interaction map of proteins identified in whole tendon.

Unconnected nodes were removed to provide clarity of the interactome. The total cluster was built with STRING (Szklarczyk et al., 2019) allowing for experimentally verified and predicted protein-protein interactions at high confidence levels (0.700). Line thickness indicates the strength of supporting data. Source data are provided in Figure 2—source data 1.

Figure 2—source data 1

Proteins identified in whole tendon samples used to create protein-protein interaction maps.

https://cdn.elifesciences.org/articles/55262/elife-55262-fig2-data1-v1.xlsx
Protein-protein interaction map of proteins identified in the interfascicular matrix.

Unconnected nodes were removed to provide clarity of the interactome. The total cluster was built with STRING (Szklarczyk et al., 2019) allowing for experimentally verified and predicted protein-protein interactions at high confidence levels (0.700). Line thickness indicates the strength of supporting data. Source data are provided in Figure 3—source data 1.

Figure 3—source data 1

Proteins identified in the interfascicular matrix used to create protein-protein interaction maps.

https://cdn.elifesciences.org/articles/55262/elife-55262-fig3-data1-v1.xlsx
Protein-protein interaction map of proteins identified in the fascicular matrix.

Unconnected nodes were removed to provide clarity of the interactome. The total cluster was built with STRING (Szklarczyk et al., 2019) allowing for experimentally verified and predicted protein-protein interactions at high confidence levels (0.700). Line thickness indicates the strength of supporting data. Source data are provided in Figure 4—source data 1.

Figure 4—source data 1

Proteins identified in the fascicular matrix used to create protein-protein interaction maps.

https://cdn.elifesciences.org/articles/55262/elife-55262-fig4-data1-v1.xlsx
Calculation of protein turnover rate from mass isotopomer distribution over time.

Relative abundance of M0 in selected decorin (A) and collagen type 1, alpha one chain (B) peptides, and resulting k values calculated by non-steady state curve fitting using ProTurn software. Fractional synthesis rates (FSR), calculated from peptide k values, demonstrate more rapid turnover of decorin (C) compared to collagen type 1, alpha one chain (D). Source data, generated by ProTurn are available in Source data 1.

Peptide turnover rate constants and corresponding protein half-lives in whole tendon.

(A) The peptide rate constants (k) for individual proteins are plotted in descending order on a logarithmic scale, with the median value represented by a red line. (B) The median k values for each protein were used to calculate protein half-life, assuming a first order reaction. (C) Protein turnover rates plotted against MatrisomeDB and PANTHER categories, with the median value represented by a red line. Due to space constraints, gene, rather than protein, names are displayed in parts A and B. Source data, generated by ProTurn are available in Source data 1.

Peptide turnover rate constants and corresponding protein half-lives in tendon phases.

(A) Rate constants (k) for collagen type I peptides identified in both tendon phases were significantly greater in the IFM compared to the FM (n = 39). (B) Median peptide decay curves for Col1a1 and Col1a2 in the FM (▲; solid line) and IFM (●; dashed line), showing goodness of fit (r2) and standard error (S.E.). (C) Turnover rate constants (k) and corresponding half-lives for proteins identified in each phase. *** indicates p<0.0001. Source data, generated by ProTurn are available in Figure 7—source data 1 for the FM and Figure 7—source data 2 for the IFM.

Figure 7—source data 1

ProTurn output for FM.

hl tab contains protein half-life information organized by peptide sequence and hl-data tab contains isotopomer relative abundance at each time point.

https://cdn.elifesciences.org/articles/55262/elife-55262-fig7-data1-v1.xlsx
Figure 7—source data 2

ProTurn output for IFM.

hl tab contains protein half-life information organized by peptide sequence and hl-data tab contains isotopomer relative abundance at each time point.

https://cdn.elifesciences.org/articles/55262/elife-55262-fig7-data2-v1.xlsx
Manual calculation of turnover rate of selected proteins.

(A) Relative abundance of M0 as a function of time, and resultant non-linear curve fitting for peptides corresponding to decorin identified in the FM (▲; solid line) and IFM (●; dashed line). (B) Resultant k values for decorin peptides were significantly greater in the IFM than in the FM (p=0.042). (C) Manually calculated turnover rate constants (k) and corresponding half-lives for tendon proteins of interest. Source data are available in Figure 8—source data 1.

Figure 8—source data 1

GraphPad Prism output showing the manually calculated relative abundance of M0 at different time points for decorin, fibromodulin, biglycan and Col3a1 peptides, and resultant K values.

https://cdn.elifesciences.org/articles/55262/elife-55262-fig8-data1-v1.xlsx

Tables

Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional
information
Strain, strain background (Rattus Norvegicus, Female)WistarCharles RiverRRID:RGD_13508588Female
Commercial assay or kitPierce protein assayThermoFisher22660
Chemical compound, drug[2H]2OCK isotopes LtDDLM-2259
Chemical compound, drugAcetonitrile with 0.1% formic acid LCMS gradeFisher Scientific10723857
Chemical compound, drugAmmonium bicarbonateSigma09830
Chemical compound, drugChondroitinase ABCSigmaC3667
Chemical compound, drugDithiothreitolMelford laboratoriesMB1015
Chemical compound, drugFormic acid (0.1% v/v) LCMS gradeFisher Scientific10188164
Chemical compound, drugGuanidine hydrocholorideSigmaG3272
Chemical compound, drugIodoacetamideSigma AldrichI1149
Chemical compound, drugRapiGest SFWaters186001861
Chemical compound, drugTrifluoroacetic acid OptimaFisher Scientific10723857
Chemical compound, drugTrifluoroacetic acid (0.1% v/v) LCMS gradeFisher Scientific10516625
Chemical compound, drugTrypsin Gold MS gradePromegaV5280
Software, algorithmPeaks Studio v8.5Bioinformatics Solutionswww.bioinfor.com/peaks-studio
Software, algorithmSTRING v11.0PMID:30476243string-db.org
Software, algorithmMatrisomeDBPMID:21937732matrisomeproject.mit.edu
Software, algorithmPANTHERPMID:23868073www.pantherdb.org
Software, algorithmProLuCIDPMID:26171723fields.scripps.edu/yates/wp/?page_id = 821
Software, algorithmProteoWizard v3PMID:28188540proteowizard.sourceforge.net/index.html
Software, algorithmProTurnPMID:24614109proturn.heartproteome.org
Software, algorithmPrism v8.2GraphPadwww.graphpad.com
OtherFilter units, Vivacon 500 10000 MWCOSartoriusVN01H02

Additional files

Source data 1

Source data for Figures 5 and 6.

ProTurn output for whole tendon. hl tab contains protein half-life information organized by peptide sequence and hl-data tab contains isotopomer relative abundance at each time point

https://cdn.elifesciences.org/articles/55262/elife-55262-data1-v1.xlsx
Transparent reporting form
https://cdn.elifesciences.org/articles/55262/elife-55262-transrepform-v1.docx

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Howard Choi
  2. Deborah Simpson
  3. Ding Wang
  4. Mark Prescott
  5. Andrew A Pitsillides
  6. Jayesh Dudhia
  7. Peter D Clegg
  8. Peipei Ping
  9. Chavaunne T Thorpe
(2020)
Heterogeneity of proteome dynamics between connective tissue phases of adult tendon
eLife 9:e55262.
https://doi.org/10.7554/eLife.55262