Research Papers: Gerotarget (Focus on Aging):
Identification of ageing-associated naturally occurring peptides in human urine
Metrics: PDF 2405 views | HTML 1667 views | ?
Esther Nkuipou-Kenfack1,2, Akshay Bhat1,3, Julie Klein4,5, Vera Jankowski6, William Mullen7, Antonia Vlahou8,9, Mohammed Dakna1, Thomas Koeck1, Joost P. Schanstra5,6, Petra Zürbig1, Karl L. Rudolph10, Björn Schumacher11, Andreas Pich2 and Harald Mischak1,7
1 Mosaiques Diagnostics GmbH, Hannover, Germany
2 Hannover Medical School, Core Facility Proteomics, Hannover, Germany
3 Charité-Universitätsmedizin Berlin, Med. Klinik IV, Berlin, Germany
4 Institut National de la Santé et de la Recherche Médicale (INSERM), Institut of Cardiovascular and Metabolic Disease, Toulouse, France
5 Université Toulouse III Paul-Sabatier, Toulouse, France
6 University Hospital RWTH Aachen, Institute for molecular cardiovascular research (IMCAR), Aachen, Germany
7 BHF Glasgow Cardiovascular Research Centre, University of Glasgow, Glasgow, United Kingdom
8 Biotechnology Division, Biomedical Research Foundation Academy of Athens, Athens, Greece
9 School of Biomedical and Healthcare Sciences, Plymouth University, Plymouth, UK
10 Leibniz Institute of Age Research, Fritz Lipmann Institute, Jena, Germany
11 Institute for Genome Stability in Ageing and Disease and Cologne Excellence Cluster for Cellular Stress Responses in Aging-Associated Diseases (CECAD) Research Center, University of Cologne, Cologne, Germany
Petra Zürbig, email:
Keywords: ageing, urine, peptidomics, collagen, proteases, systems biology, Gerotarget
Received: September 10, 2015 Accepted: September 16, 2015 Published: September 29, 2015
To assess normal and pathological peptidomic changes that may lead to an improved understanding of molecular mechanisms underlying ageing, urinarypeptidomes of 1227 healthy and 10333 diseased individuals between 20 and 86 years of age were investigated. The diseases thereby comprised diabetes mellitus, renal and cardiovascular diseases. Using age as a continuous variable, 116 peptides were identified that significantly (p < 0.05; |ρ|≥0.2) correlated with age in the healthy cohort. The same approach was applied to the diseased cohort. Upon comparison of the peptide patterns of the two cohorts 112 common age-correlated peptides were identified. These 112 peptides predominantly originated from collagen, uromodulin and fibrinogen. While most fibrillar and basement membrane collagen fragments showed a decreased age-related excretion, uromodulin, beta-2-microglobulin and fibrinogen fragments showed an increase. Peptide-based in silico protease analysis was performed and 32 proteases, including matrix metalloproteinases and cathepsins, were predicted to be involved in ageing. Identified peptides, predicted proteases and patient information were combined in a systems biology pathway analysis to identify molecular pathways associated with normal and/or pathological ageing. While perturbations in collagen homeostasis, trafficking of toll-like receptors and endosomal pathways were commonly identified, degradation of insulin-like growth factor-binding proteins was uniquely identified in pathological ageing.
Normal physiological ageing is a complex, multi-mechanistic systemic process that is influenced by genetic and environmental factors. It leads to a gradual decline in biological functions. Key molecular mechanisms identified in ageing include genomic instability, telomere attrition, loss of proteostasis and mitochondrial dysfunction . However, information on normal physiological ageing may be blurred by alterations associated with pathologies (acute and chronic) developing in parallel with ageing and it is still often unclear whether an observed molecular change is due to ageing, or is (partially) due to concomitant diseases. It is thus obvious that more efforts should be invested into understanding molecular pathways underlying ageing in both healthy and diseased individuals. These may lead to strategies for the management of pathological complications during ageing.
As ageing is a complex systemic process, “omics” approaches aiming at studying multiple features at once, have been applied with the aim to unravel novel underlying molecular processes. Proteomic studies confirmed that oxidative stress occurs ubiquitously during ageing while other events were shown to be more tissue-specific (reviewed in ). However, a shortcoming in most of these studies was the use of animal models . The scarcity of human studies can be largely attributed to the inability in obtaining appropriate tissue samples. Thus, a way forward could be the investigation of readily available body fluids.
In a first small scale study, we investigated the urinary proteome in a cohort of 324 healthy individuals between 2 to 73 years of age showing the feasibility to obtain high resolution molecular information from readily available body fluids such as urine . Meanwhile, we have accumulated multiple high-resolution urine peptidomics datasets that enable the investigation of ageing-associated changes in a large cohort . In the present study, we therefore investigated the unique urinary proteome profiles of 11560 individuals in an attempt to identify specific ageing-associated alterations and investigate pathological derailment of normal ageing. This showed that perturbations in collagen homeostasis, trafficking of toll-like receptors and endosomal pathways were associated to healthy ageing, while degradation of insulin-like growth factor-binding proteins was uniquely identified in pathological ageing
Age-correlation analysis in the healthy group
Among the 11560 individual urinary peptidomes, 1227 originated from individuals without disease and were thus considered healthy (age 20-86). Correlation analysis of 2223 individual sequenced peptides with age performed in the healthy peptidomes identified 116 significantly ageing-associated peptides (p ≤ 0.05) (Supplemental Table 1). These peptides predominantly included fragments of collagen, fibrinogen, and uromodulin. Collagen fragments comprised 83 (72%) out of the 116 peptides identified. Amongst collagen fragments, most peptides originated from fibrillar collagens (89%) including type I collagen (47%) and type III collagen (11%) while basement membrane type IV collagens alpha-1and -3 showed a low abundance (2%).
The majority of peptides (65%) showed a negative correlation with age (Supplemental Table 1). The two most negatively age-correlated peptides were two type I collagen alpha-1 fragments (ρ = -0.324, p < 0.0001 and ρ = -0.315, p < 0.0001, Supplemental Table 1) and 93% of the type I collagen fragments decreased during ageing. Other negatively age-correlated peptides originated from 5-AMP-activated protein kinase subunit gamma-3 (PRKAG3), AMP/ATP-binding subunit of AMP-activated protein kinase (AMPK) and blood-derived proteins (beta-2-microglobulin, fibrinogen alpha and beta chains). Contrarily, the two most positively age-correlated peptides were type IV collagen alpha-3 and type II collagen alpha-1 fragments (ρ = 0.504, p < 0.0001 and ρ = 0.451, p < 0.0001 respectively, Supplemental Table 1). Additionally an age-dependent increase in almost 50% of type III collagen and 83% in type II collagen fragments was observed. Other positively age-correlated peptides originated from clusterin, haptoglogin, cystatin-B, retinol-binding protein 4, CD99 antigen, and the kidney-specific peptide uromodulin.
Interestingly, several peptides that were negatively correlated with age became positively correlated upon methionine oxidation. This observation was consistent for two fragments of type I collagen alpha-1 as well as fragments of type IX collagen alpha-3, type XXV collagen alpha-1, sodium/potassium-transporting ATPase subunit gamma and retinol-binding protein 4 (Supplemental Table 1).
Age-correlation analysis in the diseased group
Next we studied the correlation of urinary peptides with ageing in the 10333 peptidomes of diseased individuals to determine potential discrepant and concerted correlations compared to healthy individuals. Individuals with pathological conditions were more likely to be older compared to healthy individuals (Table 1). Out of the 116 age-correlated peptides in healthy individuals, 112 were also found to correlate in diseased individuals. However, lower correlation coefficients were observed in the diseased compared to the healthy group (Supplemental Table 1). This observation was expected, given the assumed increased heterogeneity as a result of various underlying pathologies. The 4 peptides not confirmed in the diseased group comprised three collagen fragments and a fibrinogen alpha chain fragment and were not considered for further investigations.
Table 1: Patient characteristics
N (number of individuals)
38.6 ± 12.4
54.4 ± 15.3
* p-value <0.0001
To determine if the 112 peptides were able to distinguish between young and old individuals in both healthy and diseased groups, the abundance of these peptides was studied in a dichotomous analysis in subpopulations of young versus old (Figure 1). Proteome profiles of young compared to old healthy individuals presented more visual differences than the profiles of young compared to old diseased individuals. However, the 112 age-correlated peptides were still able to distinguish between young and old individuals in healthy and diseased individuals. Interestingly, some peptides showed similar mean amplitudes in the healthy and diseased groups including for instance collagen alpha-1 (XXV) chain (Figure 1, green arrows) whereas other peptides such as collagen alpha-1 (III) chain (Figure 1, red arrows) depicted different amplitude profiles. Differences in age-correlated peptides were further investigated by comparing the correlation coefficients of the 112 peptides in both groups. As a result, peptides could be arranged into two groups: disease-unaffected and disease-affected peptides. These were defined by a non-significant (disease-unaffected) and a significant (disease-affected) p-value in the comparison of correlation coefficients between healthy and diseased groups (Supplemental Table 1, column “healthy vs diseased”). For instance among the best correlated peptides, the correlation coefficients for collagen alpha-1 (II) chain in healthy (rho = 0.451, p < 0.0001) and diseased (rho = 0.439, p < 0.0001) individuals did not differ significantly (p = 6.21E-01) (Supplemental Table 1). An example for a disease-affected peptide is a fragment of collagen alpha-3(IV) chain with correlation coefficients in healthy (rho = 0.504, p < 0.0001) and diseased (rho = 0.420, p < 0.0001) individuals that differed significantly (p < 0.0001) (Supplemental Table 1). Further disease-unaffected peptides comprised fragments of type II alpha-1 and type III alpha-1 collagen (Figure 2A), while fragments of retinol-binding protein 4 and type I collagen alpha-1 were further disease-affected peptides (Figure 2B). Overall, 27 peptides widely represented by collagen fragments (89%) were disease-unaffected, while disease-affected peptides totalled with 85 and only included 66% of collagen fragments (Figure 3).
Figure 1: Urinary peptide marker pattern for the differentiation between healthy and diseased individuals. A. Healthy young between 20-29 years of age. B. Diseased young between 20-29 years of age. C. Healthy old from 60 years old of age and above. D. Diseased old from 60 years old of age and above. Only the mean intensity for each peptide was represented.
Figure 2: Correlation analysis of individual urinary peptides in healthy and diseased groups with age. A. Disease-unaffected peptides, collagen alpha-1(II) chain (ρHealthy = 0.451, p < 0.0001 and ρDiseased = 0.439, p < 0.0001) and collagen alpha-1(I) chain (ρHealthy = -0.224, p < 0.0001 and ρDiseased = -0.251, p < 0.0001). B. Disease-affected peptides, retinol-binding protein 4 (ρHealthy = 0.311, p < 0.0001 and ρDiseased = 0.149, p < 0.0001) and collagen alpha-1(I) chain (ρHealthy = -0.308, p < 0.0001 and ρDiseased = -0.045, p < 0.0001).
Figure 3: Comparison of age-correlated peptides identified in the healthy and diseased groups. A. Disease-unaffected peptides. B. Disease-affected peptides.
Pathology-specific investigation of age-correlated peptides in different subgroups
As the diseased group of 10333 individuals included heterogeneous pathologies (Table 2), the pathology-specific age-association of the 112 peptides was investigated in three different more homogenous disease subgroups selected from the 10333 diseased individuals and then compared to the healthy group. The cardiovascular diseases (CVD, n = 1681) subgroup included individuals with heart failure, coronary artery disease and acute coronary syndrome. The chronic kidney diseases (CKD, n = 2154) subgroup included individuals with several kidney disorders such as vasculitis and glomerulopathies, whereas the diabetes mellitus (DM, n = 1560) subgroup consisted of type 1 and type 2 DM individuals with no detectable kidney disease.
The comparison of correlation analyses in all disease subgroups and the healthy group based on the 112 identified peptides provided an assessment of the distribution of age-correlated peptides. Sixty-six, 100 and 54 of the 112 age-correlated peptides were significantly correlated to age in individuals with CVD, CKD and DM, respectively (Supplemental Table 2). Thereby the number of disease-unaffected age-correlated peptides in the disease subgroups was reduced in CVD (n = 10) and DM (n = 6) compared to CKD (n = 35) (Supplemental Table 2). There was no overlap among these disease-unaffected peptides. As seen in the full diseased cohort of 10333 individuals, disease-affected non-collagen peptides were also almost two times as abundant as disease-unaffected ones in the three disease subgroups (Figure S1). In regard of collagen fragments, in the CVD subgroup 60% of disease-unaffected peptides originated from type I collagen compared to 46% of the disease-affected peptides. In comparison, type III collagen fragments comprised 30% of disease-unaffected peptides but only 12.5% disease-affected peptides. Type I collagen fragments in the CKD subgroup represented 60% of diseased-unaffected peptides compared to 43% in the disease-affected peptides (Supplemental Figure 1).
Table 2: Different pathological conditions represented in the diseased group.
N (number of individuals)
Virus-triggered diseases (e.g. hepatitis, HIV)
Kidney diseases (transplanted)
Polycystic ovary syndrome
Pathologies related to the prostate
Prediction of protease activities
Based on the N- and C-terminal sequences of naturally occurring peptides, protease activity responsible for their generation can be predicted . The in silico prediction of ageing-related changes in the activity of proteases potentially involved in the generation of the 112 peptides was based on the cleavage site consensus sequences of proteases and mean peptide intensities in individual healthy study subjects (n = 1227). The analysis resulted in 674 protease/cleavage associations related to 37 unique proteases. Amongst those, 32 proteases showed a significant correlation of their predicted activity with age in the healthy group (Supplemental Table 3). This comprised positive age-correlated activities e.g. of meprin A beta subunit, kallikrein 5, and thrombin as well as negative age-correlated activities e.g. of neprilysin, cathepsin L1, and matrix metalloproteinase-14 (MMP-14). We next compared predicted protease activities targeting disease-affected peptides between the healthy group and the disease subgroups (Figure 4). However, while we did not observe any significant differences between healthy individuals and individuals with CVD or DM, differences in age-related activities of A disintegrin metalloproteinase with thrombospondin motifs 4 (Adamts4) and MMPs appeared to be present in individuals with CKD (Figure 4, arrows).
Figure 4: Comparison of age-correlated proteases between healthy individuals and disease subgroups. A. Cardiovascular diseases (CVD). B. Diabetes Mellitus (DM). C. Chronic kidney diseases (CKD). Arrows underscore the main changes in predicted protease activity between age-correlated disease-affected peptides in the healthy and the disease subgroups. ADAMTS4: A disintegrin and metalloproteinase with thrombospondin motifs 4; CTSB: cathepsin B; CTSK: cathepsin K; CTSL1: cathepsin L1; F2: thrombin; GZMB: granzyme B; MEP1A: meprin A subunit alpha; MEP1B: meprin A subunit beta; MME: neprilysin; PREP: prolyl endopeptidase.
Pathway enrichment analysis
Reactome pathway analysis for the identified disease-unaffected or disease-affected age-correlated peptides combined with the predicted proteases (in gene symbols) using ClueGO and CluePedia software resulted in several molecular pathways being significantly affected in the context of ageing (Figure 5, Supplemental Table 4). For disease-unaffected pathway analysis, the 27 peptides were combined with 19 predicted proteases that generated these peptide sequences. The same approach was performed for the 85 disease-affected peptides with the 32 corresponding predicted proteases. The network illustrates each pathway as individual nodes, while edges between pathways denote an approximation of biological interaction between the pathways based on the cross-pathway feature overlap. The analysis using disease-unaffected peptides revealed 6 molecular pathways associated with ageing including degradation of the extracellular matrix (ECM), activation of matrix metalloproteinases, collagen degradation, assembly of collagen fibrils, trafficking and processing of endosomal Toll-like receptors (TLRs) and endosomal/vacuolar pathway being enriched. However the analysis using disease-affected peptides, these six pathways were confirmed, and, in addition, degradation of IGF binding proteins was enriched in addition to the other 6 pathways (Figure 5).
Figure 5: Molecular pathways associated with ageing. The network represents each pathway as individual octagonal node, while the circled nodes represent the predicted proteases that were targeted from the identified urinary peptides denoted in purple diamond nodes. The edges (links) between pathways denote an approximation of biological interaction between the pathways based on the cross-pathway feature overlap. Legends for the diamond nodes with a suffix of “-C/N” represent the peptide’s cleavage site; i.e. “-C” for C-terminus and “-N” for the N-terminus.”
The urinary proteome profiles of a unique cohort of 11560 individuals with an age ranging from 20 to 86 years were analysed with the aim of detecting specific ageing-associated urinary peptides and thus expand the current knowledge on the protein level and investigate the proteomic transition from normal ageing to age-related pathological complications.
The most prominent finding of the study was that increased age is associated with a decrease in the urinary excretion of fragments from collagens forming the fibrillar structure of the ECM, including type I, II, III and V . This finding is consistent with a study reported by Zürbig et al. (2009) . While 49 fragments of mainly type I and type III collagen fragments were found to be significantly age-associated, only 15 of these peptides including fibrinogen chain and several collagen fragments of the Zürbig et al. study were sequenced and fulfilled the quality criteria of the current study. A decrease in type I and type III collagens was also observed in a study evaluating the effect of ageing on skin in a Caucasian female population . Although the study cohort comprised only 218 healthy women between 33 and 77 years of age, findings may well be extrapolated to a male population. A decrease in fibrillar collagens observed in this study may result from several processes including but not limited to impaired collagen synthesis and/or impaired degradation causing aberrant ECM remodelling . However, while the urinary excretion of the majority of fibrillar collagen fragments identified in the current study decreased with age, we also observed an age-associated increase in the urinary excretion of a few specific fibrillar collagen fragments with increasing age, especially of type II and III. In regards to type II collagen, these findings may indicate on one hand tissue and organ-dependent differences in homeostasis since type II collagen is mainly present in cartilage  and on the other hand the increased likelihood to develop osteoarthritis with advanced age as it has been shown that urinary levels of type II collagen fragments increased with osteoarthritis . The increased excretion of some type III collagen fragments could be attributed to homeostasis as type III collagen plays an important role in type I fibrillogenesis and cardiovascular development . Furthermore, all collagen fragments containing oxidised methionine were positively correlated with age. This is a novel finding that indicates an accumulation of oxidative modifications associated with age, which may lead to increased degradation. Of note: the corresponding unmodified peptide showed a decrease in urinary abundance with age. These findings may potentially also reflect progressive loss of control of oxidative stress during advancing ageing .
In addition to fibrillar collagens, the excretion of peptides from basement membrane collagens including type IV collagen alpha-1 (COL4A1) and alpha-3 (COL4A3) chains was also found to be altered. The excretion of a COL4A1-derived peptide was decreased and that of a peptide derived from COL4A3 was increased in advanced age. These findings are in agreement with the literature since COL4A1 was commonly found in the glomerular basement membrane of younger individuals whereas COL4A3 appears to be more common in adult individuals . Increased urinary excretion of type IV collagen has furthermore been associated with renal dysfunction in patients with type 2 diabetes mellitus . Hence, alterations of the basement membrane, readily observed in urine, are an important molecular event observed in ageing and renal impairment.
Of the non-collagenous peptides associated to ageing, fragments of uromodulin, beta-2-microglobulin and fibrinogen alpha and beta chains were most prevalent. The urinary excretion of most of these peptides showed a positive correlation with age. Fibrinogen, a glycoprotein involved in inflammation, and uromodulin, a kidney-specific protein, were shown to participate in renal fibrosis [15, 16]. The protein beta-2-microglobulin (B2M) is expressed in all nucleated cells and part of the light chain subunit of the major histocompatibility complex class I molecules . Plasma and serum elevations of B2M were found to be associated with a plethora of pathological conditions including renal diseases  and cardiovascular diseases . These peptides, that are in many cases also significantly associated with CKD , may reflect the reduction in kidney function observed in ageing .
The comparison of age correlation coefficients between healthy and diseased individuals enabled us to distinguish between disease-unaffected peptides reflecting processes of normal or healthy ageing and disease-affected peptides indicating events of pathological ageing. A correlation analysis in the diseased subgroups revealed greater similarity in age-correlated peptide excretions in urine between normal ageing and CKD compared to CVD and DM. This may indicate that urine does reflect the “status” of the kidney to a large degree . Potential similarities in molecular alterations were suggested by decreases in excreted fibrillar collagen fragments, which is often indicative of alterations in the ECM turnover in the diseased kidney eventually resulting in fibrosis . Furthermore, the protease analysis revealed a greater influence of CKD on age compared to CVD and DM. Our findings show that fibrosis developing in advanced age and CKD are similar whereas the similarity is less pronounced with fibrosis developing in CVD and DM.
A pathway enrichment analysis incorporating the 112 identified peptides and 32 predicted proteases suggested molecular pathways that are affected in normal and pathological ageing. Processes affected in normal ageing included perturbations in the collagen homeostasis, trafficking of toll-like receptors (TLRs) and endosomal pathways. As expected based on the abundance of collagen fragments, most of the molecular pathways found to be affected during ageing were involved in collagen homeostasis. Findings suggested accumulation of ECM or formation of fibrosis during ageing caused by a decrease in ECM degradation and an increase formation of collagen fibrils. These events result in a decrease of collagen fragments in the urine. Fibrosis is observed in renal ageing progressively degrading kidney function which potentially results in CKD . Fibrosis in the heart can cause ventricular stiffening and impairment of heart function leading to cardiovascular diseases . The enrichment analysis also indicated an impaired processing and trafficking of TLRs based-on the predicted negative age-correlation of the activities of cathepsin K, L1 and S. TLRs recognise molecular patterns that are broadly shared by pathogens and are essential for innate immune response by releasing cytokines and chemokines . It was reported that cleavage of TLRs by cathepsins is crucial for the activation of TLRs signalling . Therefore, attenuated cleavage of TLRs can contribute to perturbations in immunity in advanced age. Furthermore, cathepsins are endosomal proteases participating in diverse cellular processes including apoptosis, autophagy and necrosis . Hence perturbations of the endosomal pathway should be more investigated in ageing.
Besides, the pathway analysis also enabled the identification of processes affected in pathological ageing. In addition to the molecular mechanisms affected in normal ageing, pathway enrichment analysis suggested an elevated degradation of insulin-like growth factor (IGF)-binding proteins (IGFBPs). The increased degradation of IGFBPs was predicted based on the activities of plasmin, thrombin and matrix metalloproteinase-12. Interactions between IGFBPs and IGFs generally have inhibitory effects on IGF-dependent signalling pathways potentially leading to augmented oxidative stress and inhibition of cellular proliferation, cellular differentiation and apoptosis [28, 29]. Interactions observed between ECM proteins, proteases including plasmin and thrombin  and IGFBPs contribute to the regulation of the bioavailability of IGFs . Furthermore, increase in IGFBPs have been reported in patients with severe kidney failure  suggesting the involvement of the IGF pathway in pathology. IGF-1 has indeed been showed to decrease during ageing . Thus, the activation of IGFBPs may be an important molecular event in ageing and further investigations are well justified to elucidate interactions between IGFBPs, plamin, thrombin and the ECM.
In conclusion, urinary proteome analysis enabled the detection of ageing-associated peptides thereby generating considerable information about molecular pathways associated with normal ageing and pathological ageing. Perturbations in collagen homeostasis and trafficking of TLRs and endosomal pathways were generally observed in both normal and pathological ageing. However, increased degradation of the IGFBPs was additionally identified for the first time in ageing using urine samples. Besides, the comparison of urinary proteome profiles between healthy individuals and several diseased individuals revealed that protein fragments excreted in urine better depict similarities between normal ageing and CKD than CVD and DM. Findings demonstrated that with the help of appropriate technologies, urine can be used as a powerful biological fluid in ageing research.
MATERIALS AND METHODS
The study was designed and conducted fulfilling all of the requisites of the laws on the protection of individuals collaborating in medical research and was in accordance with the principles of the Declaration of Helsinki.
Patient characteristics and CE-MS analysis
Patient data were retrieved from the “Human urinary database” dedicated to naturally occurring urinary proteins and peptides [4, 33]. All datasets included in the study were from previous studies, and all data were anonymised. The approach, employing anonymised proteomics data from previous studies, was approved by the local ethics committee. Datasets from 11560 individuals between 20 and 86 years of age were extracted (Table 1). The present cohort was divided into two groups: healthy and diseased. The healthy group included 1227 individuals and the diseased group 10333 individuals predominantly suffering from diabetes, cardiovascular and renal diseases (Table 2).
Sample preparation and capillary electrophoresis coupled to mass spectrometry analysis
For proteomic analysis, a 0.7 mL aliquot of urine was thawed immediately before use and diluted with 0.7 mL of 2 M urea, 10 mM NH4OH containing 0.02% SDS. To remove higher molecular mass proteins, such as albumin and immunoglobulin G, the sample was ultra-filtered using Centrisart ultracentrifugation filter devices (20 kDa MWCO; Sartorius, Goettingen, Germany) at 3,000 rcf until 1.1 ml of filtrate was obtained. This filtrate was then applied onto a PD-10 desalting column (GE Healthcare, Uppsala, Sweden) equilibrated in 0.01% NH4OH in HPLC-grade in H2O (Roth, Germany) to decrease matrix effects by removing urea, electrolytes, salts, and to enrich polypeptides present. Finally, all samples were lyophilised, stored at 4°C, and suspended in HPLC-grade H2O shortly before capillary electrophoresis coupled to mass spectrometry (CE-MS) analyses, as described .
CE-MS analyses were performed using a P/ACE MDQ capillary electrophoresis system (Beckman Coulter, Fullerton, USA) on-line coupled to a microTOF MS (Bruker Daltonics, Bremen, Germany) as described previously [34, 35]. The ESI sprayer (Agilent Technologies, Palo Alto, CA, USA) was grounded, and the ion spray interface potential was set between -4 and -4.5 kV. Data acquisition and MS acquisition methods were automatically controlled by the CE via contact-close-relays. Spectra were accumulated every 3 s, over a range of m/z 350 to 3000. Accuracy, precision, selectivity, sensitivity, reproducibility, and stability of the CE-MS measurements were demonstrated elsewhere .
Mass spectral peaks representing identical molecules at different charge states were deconvoluted into single masses using MosaiquesVisu software. Only signals with z>1 observed in a minimum of 3 consecutive spectra with a signal-to-noise ratio of at least 4 were considered. Reference signals of 1770 urinary polypeptides were used for CE-time calibration by locally weighted regression. For normalisation of analytical and urine dilution variances, signal intensities were normalised relative to 29 ‘‘housekeeping’’ peptides . The obtained peak lists characterise each polypeptide by its molecular mass [Da], normalised CE migration time [min] and normalised signal intensity. All detected peptides were deposited, matched, and annotated in a Microsoft SQL database allowing further statistical analysis . For clustering, peptides in different samples were considered identical if mass deviation was < 50 ppm. CE migration time was controlled to be below 0.35 minutes after calibration.
For sequencing of peptides the urine samples were analysed on a Dionex Ultimate 3000 RSLC nano flow system (Dionex, Camberly, UK) coupled to an Orbitrap Velos MS instrument (Thermo Fisher Scientific) as described in . Data files were analysed using Proteome Discoverer 1.2 (activation type: HCD; min-max precursor mass: 790-6,000; precursor mass tolerance: 10 ppm; fragment mass tolerance: 0.05 Da; S/N threshold: 1) and were searched against the Uniprot human non-redundant database without enzyme specificity. No fixed modifications were selected, oxidation of methionine, lysine and proline were selected as variable modifications. The peptide data were extracted using high confidence peptides, defined by an Xcorr ≥ 1.9, a delta mass between experimental and theoretical mass ± 5 ppm, absence of cysteines in the sequence (since cysteines without reduction and alkylation form disulphide bonds), absence of oxidised proline in protein precursors other than collagens or elastin, and top one peptide rank filters. For further validation of obtained peptide identifications, the strict correlation between peptide charge at the working pH of 2 and CE-migration time was used to prevent false identifications . Only the sequenced peptides were further considered.
Correlation and statistical analyses
As peptide profiles across the samples were not normally distributed, we used the non-parametric Spearman’s rank correlation coefficient for estimating the correlation of individual peptides using age as a continuous variable. All peptides present in the whole population were included in the correlation analysis since a frequency threshold was not set. The statistical significance was assumed at p < 0.05. The p-value was adjusted by applying Benjamini-Hochberg [40, 41]. A cut-off value was set for the correlation analysis and the coefficient of ≥ 0.2 or ≤ -0.2 (|ρ|≥0.2) was considered for further analysis. The analysis was performed using proprietary software (R-based statistic software, version 2.15.3) and verified with MedCalc version 18.104.22.168 (MedCalc Software, Mariakerke, Belgium). Graphs were generated using MedCalc. To discriminate between peptides affected by a diseased and those unaffected, the Spearman’s rank correlation coefficient of a peptide was compared using MedCalc in healthy and diseased individuals.
In silico protease prediction
In order to link urinary fragments to the proteases involved, in silico protease mapping to urinary peptides was generated using Proteasix software as previously described . Briefly, for each of the peptides, associated-proteases were predicted for both N and C-terminal cleavage sites. In parallel, a list of >6000 random octapeptide sequences was mapped using the same protocol in order to determine the specificity of the prediction. Only protease/cleavage site associations with higher prediction score than associations with random sequences were kept for further analysis. For each protease, predicted activity in each patient was calculated based on the mean of associated peptide intensities. A parametric Pearson correlation analysis between a predicted protease and the age was performed based on the mean intensities of cleaved peptides with the age of each individual.
Pathway enrichment analysis
For elucidating molecular pathways being associated to ageing, the age-correlated peptides and the significant proteases predicted by Proteasix were subjected to the Cytoscape’s plug-ins ClueGo and CluePedia [42, 43]. Reactome pathway  served as the clustering criterion using a two-sided hypegeometry test followed by Bonferroni correction (significance level of 0.05) for identifying significantly affected pathways .
The authors would like to thank Dr. Jochen Metzger, Dr. Adela Ramírez-Torres, Justyna Siwy and Maria Frantzi for suggestions for data analysis.
This work was supported by the grant GA 31354 (CodeAge) FP7-PEOPLE-2011-ITN.
CONFLICTs OF INTEREST
The authors have declared the following potential conflict of interest: H. Mischak is the founder and co-owner of Mosaiques Diagnostics GmbH, who developed the CE-MS technology for clinical application. E. Nkuipou-Kenfack, A. Bhat, M. Dakna, T. Koeck, and P. Zürbig are employees of Mosaiques Diagnostics GmbH.
1. Lopez-Otin C, Blasco MA, Partridge L, Serrano M and Kroemer G. The hallmarks of aging. Cell. 2013;153:1194-1217.
2. Nkuipou-Kenfack E, Koeck T, Mischak H, Pich A, Schanstra JP, Zürbig P and Schumacher B. Proteome analysis in the assessment of ageing. Ageing Res Rev. 2014;18C:74-85.
3. Zürbig P, Decramer S, Dakna M, Jantos J, Good DM, Coon JJ, Bandin F, Mischak H, Bascands JL and Schanstra JP . The human urinary proteome reveals high similarity between kidney aging and chronic kidney disease. Proteomics. 2009;9:2108-2117.
4. Stalmach A, Albalat A, Mullen W and Mischak H. Recent advances in capillary electrophoresis coupled to mass spectrometry for clinical proteomic applications. Electrophoresis. 2013;34:1452-1464.
5. Klein J, Eales J, Zurbig P, Vlahou A, Mischak H and Stevens R. Proteasix: A tool for automated and large-scale prediction of proteases involved in naturally occurring peptide generation. Proteomics. 2013;13:1077-1082.
6. Mouw JK, Ou G and Weaver VM. Extracellular matrix assembly: a multiscale deconstruction. Nat Rev Mol Cell Biol. 2014;15:771-785.
7. Baroni ER, Biondo-Simoes ML, Auersvald A, Auersvald LA, Montemor Netto MR, Ortolan MC and Kohler JN. Influence of aging on the quality of the skin of white women: the role of collagen. Acta Cir Bras. 2012;27:736-740.
8. Bonnans C, Chou J and Werb Z. Remodelling the extracellular matrix in development and disease. Nat Rev Mol Cell Biol. 2014;15:786-801.
9. Wilusz RE, Zauscher S and Guilak F. Micromechanical mapping of early osteoarthritic changes in the pericellular matrix of human articular cartilage. Osteoarthritis Cartilage. 2013;21:1895-1903.
10. Ishijima M, Kaneko H and Kaneko K. The evolving role of biomarkers for osteoarthritis. Ther Adv Musculoskelet Dis. 2014;6:144-153.
11. Liu X, Wu H, Byrne M, Krane S and Jaenisch R. Type III collagen is crucial for collagen I fibrillogenesis and for normal cardiovascular development. Proc Natl Acad Sci U S A. 1997;94:1852-1856.
12. Stadtman ER, Van RH, Richardson A, Wehr NB and Levine RL. Methionine oxidation and aging. Biochim Biophys Acta. 2005;1703:135-140.
13. Hudson BG, Tryggvason K, Sundaramoorthy M and Neilson EG. Alport’s syndrome, Goodpasture’s syndrome, and type IV collagen. N Engl J Med. 2003;348:2543-2556.
14. Araki S, Haneda M, Koya D, Isshiki K, Kume S, Sugimoto T, Kawai H, Nishio Y, Kashiwagi A, Uzu T and Maegawa H . Association between urinary type IV collagen level and deterioration of renal function in type 2 diabetic patients without overt proteinuria. Diabetes Care. 2010;33:1805-1810.
15. Docherty NG and Godson C. Fibrinogen as a damage-associated mitogenic signal for the renal fibroblast. Kidney Int. 2011;80:1014-1016.
16. Scolari F, Izzi C and Ghiggeri GM. Uromodulin: from monogenic to multifactorial diseases. Nephrol Dial Transplant. 2015;30:1250-1256.
17. Chiou SJ and Chen CH. Decipher beta2-microglobulin: gain- or loss-of-function (a mini-review). Med Sci Monit Basic Res. 2013;19:271-273.
18. Astor BC, Shafi T, Hoogeveen RC, Matsushita K, Ballantyne CM, Inker LA and Coresh J. Novel markers of kidney function as predictors of ESRD, cardiovascular disease, and mortality in the general population. Am J Kidney Dis. 2012;59:653-662.
19. Kals J, Zagura M, Serg M, Kampus P, Zilmer K, Unt E, Lieberg J, Eha J, Peetsalu A and Zilmer M. beta2-microglobulin, a novel biomarker of peripheral arterial disease, independently predicts aortic stiffness in these patients. Scand J Clin Lab Invest. 2011;71:257-263.
20. Schanstra JP, Zurbig P, Alkhalaf A, Argiles A, Bakker SJ, Beige J, Bilo HJ, Chatzikyrkou C, Dakna M, Dawson J, Delles C, Haller H, Haubitz M, et al. Diagnosis and Prediction of CKD Progression by Assessment of Urinary Peptides. J Am Soc Nephrol. 2015;26:1999-2010.
21. Zhou XJ, Rakheja D, Yu X, Saxena R, Vaziri ND and Silva FG. The aging kidney. Kidney Int. 2008;74:710-720.
22. Pejchinovski M, Hrnjez D, Ramirez-Torres A, Bitsika V, Mermelekas G, Vlahou A, Zürbig P, Mischak H, Metzger J and Koeck T. Capillary zone electrophoresis on-line coupled to mass spectrometry: a perspective application for clinical proteomics. Proteomics Clin Appl. 2015;9:453-468.
23. Anderson S, Halter JB, Hazzard WR, Himmelfarb J, Horne FM, Kaysen GA, Kusek JW, Nayfield SG, Schmader K, Tian Y, Ashworth JR, Clayton CP, Parker RP, et al. Prediction, progression, and outcomes of chronic kidney disease in older adults. J Am Soc Nephrol. 2009;20:1199-1209.
24. Biernacka A and Frangogiannis NG. Aging and Cardiac Fibrosis. Aging Dis. 2011;2:158-173.
25. Ma J, Chadban SJ, Zhao CY, Chen X, Kwan T, Panchapakesan U, Pollock CA and Wu H. TLR4 activation promotes podocyte injury and interstitial fibrosis in diabetic nephropathy. PLoS ONE. 2014;9:e97985.
26. Garcia-Cattaneo A, Gobert FX, Muller M, Toscano F, Flores M, Lescure A, Del Nery E and Benaroch P. Cleavage of Toll-like receptor 3 by cathepsins B and H is essential for signaling. Proc Natl Acad Sci U S A. 2012;109:9053-9058.
27. Turk B and Turk V. Lysosomes as “suicide bags” in cell death: myth or reality? J Biol Chem. 2009;284:21783-21787.
28. Kim KS, Seu YB, Baek SH, Kim MJ, Kim KJ, Kim JH and Kim JR. Induction of cellular senescence by insulin-like growth factor binding protein-5 through a p53-dependent mechanism. Mol Biol Cell. 2007;18:4543-4552.
29. Gong Z, Kennedy O, Sun H, Wu Y, Williams GA, Klein L, Cardoso L, Matheny RW, Hubbard GB, Ikeno Y, Farrar RP, Schaffler MB, Adamo ML, et al. Reductions in serum IGF-1 during aging impair health span. Aging Cell. 2014;13:408-418.
30. Booth BA, Boes M, Dake BL, Knudtson KL and Bar RS. IGFBP-3 binding to endothelial cells inhibits plasmin and thrombin proteolysis. Am J Physiol Endocrinol Metab. 2002;282:E52-E58.
31. Ivarsen P, Chen JW, Tietze I, Christiansen JS, Flyvbjerg A and Frystyk J. Marked reductions in bioactive insulin-like growth factor I (IGF-I) during hemodialysis. Growth Horm IGF Res. 2010;20:156-161.
32. Higashi Y, Sukhanov S, Anwar A, Shai SY and Delafontaine P. Aging, atherosclerosis, and IGF-1. J Gerontol A Biol Sci Med Sci. 2012;67:626-639.
33. Siwy J, Mullen W, Golovko I, Franke J and Zürbig P. Human urinary peptide database for multiple disease biomarker discovery. Proteomics Clin Appl. 2011;5:367-374.
34. Mischak H and Schanstra JP. CE-MS in biomarker discovery, validation, and clinical application. Proteomics Clin Appl. 2011;5:9-23.
35. Arab S, Gramolini AO, Ping P, Kislinger T, Stanley B, van Eyk J, Ouzounian M, MacLennan DH, Emili A and Liu PP. Cardiovascular proteomics: tools to develop novel biomarkers and potential applications. J Am Coll Cardiol. 2006;48:1733-1741.
36. Jantos-Siwy J, Schiffer E, Brand K, Schumann G, Rossing K, Delles C, Mischak H and Metzger J . Quantitative Urinary Proteome Analysis for Biomarker Evaluation in Chronic Kidney Disease. J Proteome Res. 2009;8:268-281.
37. Gu YM, Thijs L, Liu YP, Zhang Z, Jacobs L, Koeck T, Zürbig P, Lichtinghagen R, Brand K, Kuznetsova T, Olivi L, Verhamme P, Delles C, et al. The urinary proteome as correlate and predictor of renal function in a population study. Nephrol Dial Transplant. 2014;29:2260-2268.
38. Kuznetsova T, Mischak H, Mullen W and Staessen JA. Urinary proteome analysis in hypertensive patients with left ventricular diastolic dysfunction. Eur Heart J. 2012;33:2342-2350.
39. Zürbig P, Renfrow M.B., Schiffer E, Novak J, Walden M, Wittke S, Just I, Pelzing M, Neusüss C, Theodorescu D, Root KE, Ross MM and Mischak H. Biomarker discovery by CE-MS enables sequence analysis via MS/MS with platform-independent separation. Electrophoresis. 2006;27:2111-2125.
40. Benjamini Y and Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc B (Methodological). 1995;57:125-133.
41. Benjamini Y and Yekutieli D. The control of the false discovery rate under dependency. Ann Stat. 2001;29:1165-1188.
42. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, Fridman WH, Pagès F, Trajanoski Z and Galon J. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25:1091-1093.
43. Bindea G, Galon J and Mlecnik B. CluePedia Cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinformatics. 2013;29:661-663.
44. Haw R, Hermjakob H, D’Eustachio P and Stein L. Reactome pathway analysis to enrich biological discovery in proteomics data sets. Proteomics. 2011;11:3598-3613.
45. Saito M, Kimoto M, Araki T, Shimada Y, Fujii R, Oofusa K, Hide M, Usui T, Yoshizato K. Proteome analysis of gelatin-bound urinary proteins from patients with bladder cancers. Eur Urol. 2005;48:865-871.
All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 License.