Urine proteomes of healthy aging humans reveal extracellular matrix (ECM) alterations and immune system dysfunction.

Aging is a complex physiological process that poses considerable conundrums to rapidly aging societies. For example, the risk of dying from cardiovascular diseases and/or cancer steadily declines for people after their 60s, and other causes of death predominate for seniors older than 80 years of age. Thus, physiological aging presents numerous unanswered questions, particularly with regard to changing metabolic patterns. Urine proteomics analysis is becoming a non-invasive and reproducible diagnostic method. We investigated the urine proteomes in healthy elderly people to determine which metabolic processes were weakened or strengthened in aging humans. Urine samples from 37 healthy volunteers aged 19-90 years (19 men, 18 women) were analyzed for protein expression by liquid chromatography-tandem mass spectrometry. This generated a list of 19 proteins that were differentially expressed in different age groups (young, intermediate, and old age). In particular, the oldest group showed protein changes reflective of altered extracellular matrix turnover and declining immune function, in which changes corresponded to reported changes in cardiovascular tissue remodeling and immune disorders in the elderly. Thus, urinary proteome changes in the elderly appear to reflect the physiological processes of aging and are particularly clearly represented in the circulatory and immune systems. Detailed identification of "protein trails" creates a more global picture of metabolic changes that occur in the elderly.

Introduction

Aging is a complex physiological process that poses many problems and conundrums for rapidly aging societies. For example, the risk of dying from cardiovascular diseases and/or cancer steadily declines for people after their 60s, and in seniors older than 80 years of age, other causes of death are predominant (de los Campos et al. 2012). Thus, diagnoses and research in which heart disease and cancer are the focus are ineffectual for half of the patients older than 80 years, and physiological aging raises many unanswered questions, particularly with regard to changing metabolic patterns.

Urine proteomics has become a non-invasive, reproducible, and easy-to-frequent repetition diagnostic method. In this study, we show that there are urine proteome changes with age in healthy humans and that this may be useful for understanding certain pathogenic processes in aging. This investigation was based on the following considerations.

Fluid biopsy

Urine proteome analysis is a rapidly developing analytical tool that is aimed at establishing new, non-invasive “fluid biopsy” diagnostic methods for renal and urinary tract diseases (Tyan et al. 2006; Varghese et al. 2007). Although this methodology is still in its infancy, knowledge on urine proteomes in health and disease is growing rapidly due to advancements in instrumentation, data analysis techniques, and the relative ease of acquiring clinically relevant samples. Thus, urinary proteome analysis is the subject of a rapidly growing number of studies, and this topic is often reviewed (Pisitkun et al. 2006).

Numerous potential markers for different pathologies, not only for the kidney, but also for the urinary tract, vascular system, and other organs, have been found (van Eyk and Dunn 2008). For example, the urine proteomes of small cohorts of patients studied by capillary electrophoresis coupled with mass spectrometry have established potential markers for ADPKD (Pavik et al. 2012; Kistler et al. 2009). The first urine proteomic study was done in 1979 using 2-D PAGE electrophoresis (Anderson et al. 1979). Later, numerous hybrid techniques were used for different diseases (including uronephrological cases), such as MALDI-TOF, SELDI-TOF, CE–MS/MS, and LC–ESI-TOF-MS/M (M'Koma et al. 2007; Mosley et al. 2006; Tyan et al. 2006; Julian et al. 2007). Depending on the methodology used, the number of detected proteins can vary greatly from 200 to 2,500 (Pisitkun et al. 2006; Weissinger et al. 2004; Barratt and Topham 2007). Thus, urine proteomics may provide a more “global view” of underlying physiological processes that are not necessarily disease specific.

Consequences of aging in general

Aging influences many biological processes (vB Hjelmborg et al. 2006). Changes in these processes are directed by several guiding physiologic principles: diminished flexibility (Lipsitz and Goldberger 1992), as manifested for example by less variable heart rate responses; changes in circadian patterns (Veldhuis 1997); and a loss of physiologic reserves required to cope with challenges to homeostasis (Cowdry 1942). While the progression of certain biological changes is linear over time, such as graying of hair, reduced skin elasticity, and vascular aging, other biological variables, such as lung capacity and sensory and cognitive functions, do not always correlate with chronological age (Bulpitt et al. 1994). Aging changes progress in a non-uniform manner, and organs age at different rates and are influenced by many factors, including genetics, lifestyle, environmental factors, and a general loss of functional reserves (vB Hjelmborg et al. 2006).

In the elderly, these progressive changes result in increased susceptibility to numerous diseases (Martin and Sheaff 2007). For example, with age, the renal mass tends to become smaller and function declines. The older kidney is more prone to nephrotoxicity related to medications (Michel and Kelly 1998) or application of intravenous contrast dye (Rich and Crecelius 1990). However, much of the research on age-related changes in humans has focused on particular disease processes, in which processes are often quite complex (e.g., cardiovascular disease, immunological disorders) and are affected by numerous factors wholly unrelated to human aging. Thus, more investigations are needed on “normal aging” in humans to determine which physiological processes are most affected by this natural chronological phenomenon.

Currently, more and more studies are using urinary proteome analysis applied to clinical research. The results have often been distant from what was expected, probably due to the presence of proteins/peptides whose levels change significantly with aging. More basic research on the effects of factors unrelated to disease, such as aging, on urine proteomes may aid in selecting control groups for clinical trials. More importantly, a more global view of urine proteome changes with human aging may highlight previously unexpected physiological changes. Thus, we investigated if urine proteomes changed significantly in healthy elderly people and if these would highlight metabolic processes that weakened or strengthened in aging humans.

Material and methods

Urine samples

Urine samples were collected from 37 healthy volunteers aged 19–90 years (19 men and 18 women). The exclusion criteria included the following: current infection of the urinary tract, macroscopic hematuria, diabetes mellitus, malignancy of the urinary tract or generalized malignancy of another system, post-organ transplantation status, and current pregnancy.

Study protocol was approved by our local ethics committee. Informed consent was obtained from all participants. This study was performed in accordance with the Declaration of Helsinki Principles.

Urine collection

Samples were collected from the 37 healthy donors using a uniform protocol. The second or third morning midstream urine was collected from all participants at 1 and 3 h after previous micturition. Sterile urine containers were used to collect samples. Sample pH was stabilized at 7.2 by adding one tenth of the volume of 1 M HEPES (pH 7.2) immediately after collection. Further sample preparation steps were carried out within 1 h after collection during which time the sample was kept at room temperature. Samples were vortexed for 2 min, centrifuged (3000×g, room temperature) for 10 min to clear debris, filtered through a 0.4-μm filter (Rotilabo-Spritzenfilter, P819.1, Roth), and separated into 1-ml aliquots. Sample aliquots were stored at −80 °C for later use.

Sample filtration

Membrane filters with a 10 kDa cutoff (Amicon Ultra-0.5, UFC501096, Millipore) were washed twice with MilliQ water prior to use. Urine was centrifuged through the membrane at 14,000×g for 15 min. Next, 500 μl of MQ was added to the retentate and centrifugation was repeated. To recover a concentrated and desalted sample, the filter was placed upside down in a clean microtube and centrifuged for 2 min at 1,000×g. Protein concentration was determined by the Bradford method. Samples aliquots were stored in siliconized microtubes (Sigma-Aldrich) at −80 °C for later use.

Urine proteome preparation

A total of 4 μg of protein from each sample was diluted in 100 mM NH4HCO3 to a final volume of 10 μl. DTT (100 mM) was added to a sample for a final concentration of 10 mM and incubated at 56 °C for 40 min. To block reduced cysteines, 0.5 M iodoacetamide at a final concentration of 50 mM was added, and the sample was incubated at room temperature for 30 min in the dark. Trypsin (Promega) was added at a 1:20 v/v ratio and incubated at 37 °C overnight. Finally, trifluoroacetic acid was added to digested protein samples to reduce the pH to 2 and inactivate trypsin. Resulting peptides were analyzed by liquid chromatography–tandem mass spectrometry (LC–MS/MS) for peptide identification and LC–MS for relative quantitation. Data were acquired in separate experiments as described below.

LC–MS settings

LC–MS proteomic analyses of urine samples used a LTQ-Orbitrap Velos mass spectrometer (Thermo Scientific) coupled with a nanoACQUITY (Waters Corporation) LC system. Spectrometer parameters were as follows: polarity mode, positive; and capillary voltage, 1.5 kV. A sample was first applied to the nanoACQUITY UPLC Trapping Column (Waters) using water containing 0.1 % formic acid as the mobile phase. Next, the peptide mixture was transferred to the nanoACQUITY UPLC BEH C18 Column (Waters, 75-μm inner diameter; 250-mm long) using an acetonitrile gradient (5–35 % acetonitrile over 160 min) in the presence of 0.1 % formic acid with a flow rate of 250 nl/min and eluted directly to the ion source of the mass spectrometer. Each LC run was preceded by a blank run to ensure that there was no carry-over of material from the previous analysis.

Qualitative analyses (i.e., peptide and protein identification) used pooled urine samples in data-dependent MS-to-MS/MS acquisition mode. Collision-induced dissociation was used, and up to five MS/MS processes were allowed for each MS scan. To increase the numbers of peptide identifications, three LC–MS/MS runs were performed per pooled sample, with each covering one of three m/z value ranges: 300–600, 500–900, and 800–2,000. This approach substantially improved the coverage for protein identification.

Quantitative analyses for individuals used separate survey scan LC–MS runs with an m/z measurement range of 300–2,000 using the same acetonitrile gradient as for qualitative LC–MS/MS runs. The data-dependent MS-to-MS/MS switch was disabled, and the spectrometer resolution was set to 15,000.

Qualitative MS data processing and database search

The acquired MS/MS raw data files were processed to produce peak lists with Mascot Distiller software (version 2.2.1, Matrix Science). The resulting ion lists were searched using the Mascot search engine (version 2.2.03, Matrix Science) against a database comprising all human protein entries from the SwissProt database (20,273 sequences) and their reversed versions. The search parameters were as follows: enzyme specificity, semitrypsin; fixed modification, carbamidomethylation (C); variable modifications, oxidation (M) and carbamidomethylation (K); and protein mass, unrestricted. The peptide and fragment ion mass tolerances (6 ppm and 0.4 Da, respectively) were established separately for individual LC–MS/MS runs using a procedure involving two database searches separated by a mass measurement error recalibration step, as previously described in Mikula et al. (2010).

Mascot search result processing used MScan, a proprietary software tool implemented in the Java programming language (http://proteom.ibb.waw.pl/). The statistical significance of peptide identifications was assessed using a joined target/decoy database search approach and a procedure that provided q-value estimates for each peptide spectrum match (PSMs) in the data set (Mikula et al. 2010; Käll et al. 2008). Only PSMs with q-values of ≤0.01 and proteins represented by at least two peptides were regarded as confidently identified. Proteins identified by a subset of peptides from another protein were excluded from analysis, and proteins matching the same set of peptides were clustered into single groups.

Quantitative MS data processing

Peptide identifications from all LC–MS/MS runs were merged into a common list, which was next overlaid onto 2-D maps generated from the LC–MS profile data of individual samples. A more detailed description of the feature extraction procedure was previously reported (Bakun et al. 2009). Briefly, the list of identified peptides was used to tag the corresponding peptide-related ion spectra based on m/z differences, deviations from the predicted elution times, and the match between the theoretical and observed isotopic envelopes. The maximum deviation accepted in m/z and retention time were established separately for each of the processed LC–MS spectra, in order to take into account possible variations in mass measurement accuracy and chromatographic separation between runs. First, an initial search with wide tolerances (±25 ppm and ±10 min, respectively, for m/z and RT) and restrictive parameters of isotopic envelope fits (at least four peaks detected, coefficient of determination ≥ 0.95) was performed. Next, non-linear mass and time calibration function were calculated using locally weighted regression smoother (LOESS), and the search was repeated with narrowed tolerances and relaxed fit requirements. The mean values of the estimated tolerances in this study were equal to ±17 ppm and ±6 min. Finally, the relative abundances of peptide ions were determined as the heights of 2-D fits to the most prominent peaks of the tagged isotopic envelopes.

After the quantitative feature extraction step, a filter was applied to the data set to remove peptide ions with missing values in >20 % of the samples. To avoid removing strictly differential peptides (i.e., detected or missing only in one of the studied groups of samples), a χ 2 statistic was used to assess significant inequalities in the distribution of missing values in the groups of samples. Next, a modified k-nearest neighbors procedure was used to impute the remaining missing values. The applied modifications included the use of variable-size neighborhoods and correlation-based metrics rather than Euclidean distances.

To minimize the effects of non-biological sources of variation, log-transformed peptide abundances were normalized by fitting a robust LOESS between individual samples and a median pseudo-sample. The parameters for the fit were established using a set of features that exhibited low variance in the non-normalized data and then applied to the whole data set. Finally, the normalized peptide-level data were rolled up to relative protein abundances. The procedure used involved rescaling abundances of peptides that originated from the same protein to a common level, followed by computing their median value.

Statistical analysis of quantitative MS results

To select differentially expressed proteins, we used a non-parametric ANOVA-based resampling statistical test. The resulting p values were then corrected for multiple-hypothesis testing using a two-step Benjamini–Hochberg procedure that controlled for the false discovery rate (Benjamini and Hochberg 1995). Relative protein abundances with adjusted p values ≤ 0.05 and fold-change (FC) values ≥ 1.5 were considered as significantly changed in at least one of the studied groups.

For graphical summaries and to evaluate the relationships among the studied samples, we used principal component analysis and Ward's linkage hierarchical clustering. All statistical analyses used proprietary software running in the MATLAB (MathWorks) environment (MStat, available at http://proteom.ibb.waw.pl/).

Results

Second morning urine samples were collected from 37 healthy volunteers aged 19–90 in the Clinical Hospital of Medical University of Warsaw. The protein fractions obtained after filtration, reduction, alkylation, and trypsin cleavage (described in “Materials and methods”) were appropriately aliquoted to avoid multiple freeze–thaw cycles and stored at −80 °C until use.

Protein identification

LC–MS/MS measurements of pooled urine samples were conducted for protein identification, resulting in the acquisition of 18,336 MS/MS spectra. A search against the SwissProt database using the Mascot engine confidently identified a set of 2,714 peptides with an estimated false discovery rate of 0.01 (Supplementary Table S1). In total, 578 proteins were identified, of which 257 were represented by at least two peptides (Supplementary Table S2).

Quantitative analysis

To obtain quantitative information, 0.5 μg protein of each of the 37 samples was subjected to LC–MS analysis. The samples were assigned to three separate groups, according to the age distribution of our study subjects (Fig. 1): age range 19–26 (average 21, 13 samples), age range 45–54 (average 49, 12 samples), and age range 72–90 (average 79, 12 samples). Demographic characteristics and renal function results for the three age groups (further referred as, respectively, AVG_21, AVG_49, and AVG_79) are summarized in Table 1.

Fig. 1

Three sample groups established according to the age distribution of the 37 healthy volunteers participating in the study: AVG_21 (range, 19–26 years; mean, 21), AVG_49 (range, 45–54 years; mean, 49), and AVG_79 (range, 72–90 years; mean, 79)

Table 1

Demographic characteristics and renal function of three study groups

	AVG_21	AVG_49	AVG_79
n	13	12	12
M/F	7/6	6/6	6/6
Age (years)	21 ± 2	49 ± 3	79 ± 6
GFR (ml/min/1.73 m2)	112 ± 9	100 ± 39	83 ± 13
Serum creatinine (μmol/l)	68 ± 9	69 ± 19	70 ± 15
CRP	2 ± 1	2 ± 1	2 ± 2

Label-free quantitative analysis was based on unique peptides (i.e., peptides unambiguously assigned to proteins) with confident MS/MS identifications. Proteins represented by at least two such peptides were quantified after overlying the qualitative information obtained from LC–MS/MS runs onto the LC–MS spectra of individual samples (see “Material and methods”). The resulting relative peptide and protein abundances are supplied in supplementary materials (Supplementary Tables S3 and S4).

Further statistical analyses were conducted using a resampling ANOVA-based test. Only proteins with corresponding p values ≤ 0.05 (after the adjustment for multiple-hypothesis testing) and FC values ≥ 1.5 were deemed significant. A simultaneous comparison of all studied groups, followed by three pair-wise comparisons of individual groups, resulted in the selection of 24 age-related proteins (Table 2). The largest set, consisting of 22 proteins, differentiated the two extreme age groups: AVG_21 and AVG_79. Of them, 8 showed increased abundance in the AVG_79 group, while the levels of the remaining 14 proteins were lower in this group. The comparison of the AVG_49 group and the AVG_79 group revealed nine differential proteins. Seven of these proteins were in common with those selected in the previous comparison and exhibited a consistent directionality of changes with respect to the AVG_79 group in both comparisons (Table 2 and Fig. 2). Under the applied threshold of significance, no differentially expressed proteins were identified in the last pair-wise comparison, involving AVG_21 and AVG_49 groups.

Table 2

Differentially expressed proteins in urine proteomes based on statistical analysis of three age groups

SwissProt ACC	Protein name	AVG_49 vs. AVG_21		AVG_79 vs. AVG_49		AVG_79 vs. AVG_21
		FC	Adjusted p value	FC	Adjusted p value	FC	Adjusted p value
P01877	Ig alpha-2 chain C region	–	–	2.13	0.0478	3.36	0.0033
P02763	Alpha-1-acid glycoprotein 1	–	–	3.6	0.0041	2.98	0.0096
P24855	Deoxyribonuclease-1	–	–	0.55	0.0209	0.5	0.0096
P26992	Ciliary neurotrophic factor receptor alpha	–	–	0.61	0.0078	0.56	0.0096
P25311	Zinc-alpha-2-glycoprotein	–	–	3.47	0.0041	2.4	0.0150
P19652	Alpha-1-acid glycoprotein 2	–	–	2.58	0.0041	1.68	0.0239
Q6UXB8	Peptidase inhibitor 16	–	–	2.73	0.0041	1.94	0.0396
P04217	Alpha-1B-glycoprotein	–	–	3.75	0.0213	–	–
P01833	Polymeric immunoglobulin receptor	–	–	0.52	0.0494	–	–
O75882	Attractin	–	–	–	–	0.6	0.0096
Q6GTX8	Leukocyte-associated immunoglobulin-like receptor 1	–	–	–	–	0.41	0.0004
P09486	SPARC	–	–	–	–	0.62	0.0117
P10451	Osteopontin	–	–	–	–	0.51	0.0239
Q9BY67	Cell adhesion molecule 1	–	–	–	–	0.55	0.0265
P05060	Secretogranin-1	–	–	–	–	0.38	0.0239
P00747	Plasminogen	–	–	–	–	0.51	0.0239
Q12805	EGF-containing fibulin-like extracellular matrix protein 1	–	–	–	–	1.51	0.0239
Q6UXG3	CMRF35-like molecule 9	–	–	–	–	0.62	0.0239
P01860	Ig gamma-3 chain C region	–	–	–	–	0.55	0.0370
O95967	EGF-containing fibulin-like extracellular matrix protein 2	–	–	–	–	0.59	0.0301
P01133	Pro-epidermal growth factor	–	–	–	–	0.49	0.0361
P00450	Ceruloplasmin	–	–	–	–	1.55	0.0480
Q9NQ38	Serine protease inhibitor Kazal-type 5	–	–	–	–	0.62	0.0480
P02760	AMBP protein	–	–	–	–	1.52	0.0396

Fig. 2

Boxplots of seven urine proteins commonly identified as differentially expressed in two pair-wise comparisons: AVG_79 vs. AVG_21 and AVG_79 vs. AVG_49

The results of the statistical tests indicate that the observed differential expression can be mostly attributed to large changes in protein levels between the two extreme groups. This statement is also supported by the results of hierarchical clustering (Fig. 3) and principal component analysis (Supplementary Figure 1). In both cases, only the samples belonging to AVG_21 and AVG_79 groups are well separated, while AVG_49 samples tend to cluster with those of younger people.

Fig. 3

Hierarchical clustering of urine samples from groups: AVG_21 (red, 13 samples), AVG_49 (blue, 12 samples), and AVG_79 sample sets (green, 12 samples). Relative abundances of 24 proteins showing statistical significance (adjusted p value ≤ 0.05; FC ≥ 1.5) in at least one pair-wise comparison between groups (Table 2) were used to generate the cluster tree

Due to limited access to informed consent samples, no confirmatory sample set could be extracted, and the necessary confirmation step assuming the use of an independent set of samples could not be carried out within the frame of the present study. Lack of such confirmatory data set is a shortcoming of the present study.

Discussion

Our analysis of the urine proteomes of 37 healthy volunteers of both sexes divided into three discontinuous age groups (Table 1) revealed 24 proteins that differentiated between the age groups of 21 years (range, 19–26 years), 49 years (range, 45–54 years), and 79 years (range, 72–90 years) (Table 2). Of them, 22 proteins were differentially expressed in the urine proteome of the youngest group compared with that of the oldest group, and 9 proteins were differentially expressed in the intermediate age group compared with the oldest group.

The profiles of the proteins that differed in the elderly included those involved in tissue remodeling processes, increased likelihood of coagulation and fibrinolysis, low-grade inflammation, and immune and metabolic dysregulation. These changes are in keeping with the observed increases in incidence with age of cardiovascular diseases, autoimmune diseases, and metabolic disorders. They are also in line with current understanding of the pathophysiology of tissue remodeling, prothrombotic status, and the occurrence of low-grade inflammation in the elderly (Meyer et al. 1998) and age-associated dysregulation of immune system function (Yao et al. 2011) and autoimmune disease.

Our major findings come under the following broad categories.

Tissue remodeling

We found increased expression of peptidase inhibitor 16 (PI16) in the urine proteome in our group of clinically healthy elderly people who were free of chronic diseases, current infections, inflammatory disorders, and malignancies as compared to the youngest and intermediate age groups. These elderly people also had reduced expressions of EGF-containing fibulin-like extracellular matrix protein 2 and secreted protein acidic and rich in cysteine (SPARC). These changes are consistent with current understanding of the mechanisms of aging tissues, particularly in the context of cardiovascular tissue aging (Chen and Frangogiannis 2010). Physiologically, about one third of the cells in the myocardium comprise contractile tissue. The remaining tissue comprises connective tissue components, including fibroblasts and extracellular matrix proteins, such as collagen types I and III, elastin, and fibronectin (Lindsey et al. 2005).

It is thought that the disproportionate increase in the volume of non-contractile cells (predominantly fibroblasts) and an increased accumulation of collagen that is modified both in the intercellular and perivascular spaces may be responsible for ventricular dysfunction and disturbances in coronary hemodynamics (Weber 1989). These are characterized by an increase in extracellular matrix proteins in the components of connective tissue. Under physiological conditions, the metabolism of extracellular matrix proteins and collagen contents is stable. However, with arterial hypertension and the aging process, an imbalance occurs between collagen synthesis and degradation. With high blood pressure, there is an increase in both the synthesis and degradation of extracellular matrix components, with an emphasis on synthesis. In contrast, during the course of aging, these two activities are reduced, particularly with regard to protein degradation, which also leads to a reduced susceptibility of cardiovascular structures to injury and myocardial fibrosis.

These processes depend on the metabolic levels of cytokines and growth factors, including PDGF and VEGF. These are important for proper wound healing and regulate tissue regeneration. SPARC is a PDGF and VEGF antagonist that reduces their mitogenic activities. A SPARC deficiency may be the cause of the excessive accumulation of extracellular matrix proteins, inadequate regulation of healing along with subsequent fibrosis, and reduced tissue susceptibility (Jørgensen et al. 2009). Also, reduced amounts of EGF-containing fibulin-like extracellular matrix protein 2, which is necessary for the proper incorporation of elastin into extracellular matrix, may be a factor in old age for weakening blood vessels and skin connective tissues (cutis laxa) and result in greater tissue elasticity, which includes more collagen converted to elastin (Dasouki et al. 2007).

These processes may also be dependent on protease inhibitors. One of these is protease inhibitor 16 (PI16); the expression of which increased in our oldest group. In animal studies, increased PI 16 expression was associated with cardiac hypertrophy. PI16 is rapidly secreted following stimulation by serum in primary cardiomyocytes in vitro. In vivo, transgenic mice that overexpress PI16 in cardiomyocytes show normal heart function but have smaller hearts with hypotrophic cardiomyocytes. PI16 expression is enhanced in the hypertrophic and failing heart and suppresses cardiomyocyte growth (Frost and Engelhardt 2007). PI16 is a putative serine protease inhibitor that may also be used as a marker following prostatectomy for prostate cancer (Reeves et al. 2006).

Prothrombotic status

Our finding of reduced plasminogen expression in the urine of elderly people is in line with the reported decreases in the concentrations of plasminogen and plasminogen activator inhibitor in the blood of elderly people (Tait et al. 1992). This phenomenon can be explained by both a decrease in plasminogen activator inhibitor concentrations and an increased prothrombotic tendency in old age. Plasminogen is a major protein of the fibrinolytic system and has a high molecular weight. Its relatively high concentrations in serum and its long half-life would not support its role as regulating factor of fibrinolysis (Kołodziejczyk and Wachowicz 2009). These characteristics may qualify plasminogen as an aging marker as opposed to a plasminogen activator inhibitor, whose concentrations and half-life are low.

Low-grade inflammation

Proteins that could distinguish older people also have roles in inflammatory processes. Features of low-grade inflammation are found as a trait of aging in lung tissue (Meyer et al. 1998). The urine proteome of older people significantly differed from that of younger people for proteins that included alpha 1 acidic glycoprotein 1 and alpha 1 acidic glycoprotein 2. These belong to acute phase proteins, and their expression was increased. Serine protease inhibitor Kazal-type 5 is probably an important anti-inflammatory agent that protects mucous membranes and defends against the invasion of microorganisms. It also maintains the structure and function of the skin by regulating the activity of proteases that are responsible for the actions of defense and desquamation-involved proteases (Jayakumar et al. 2005; Deraison et al. 2007). Its reduced levels may be indicative of subclinical inflammation in old age, as described for lung tissue (Meyer et al. 1998). This is in line with lower SPARC expression accompanied by a pro-fibrotic cytokine profile (Jørgensen et al. 2009).

Dysregulation of the immune system

Changes in the urinary proteome of the elderly may also reflect immune system disturbances. We found increased expression of IgA-1 chain C region in the proteome of older people compared to younger people and reduced IgG-3 chain C region. IgA is the major immunoglobulin class that is secreted into body fluids. IgA has a dual role of combating local infections and preventing foreign antigens' access to the general immune system (Corthésy 2013). IgG antibodies are the primary humoral defense mechanism in the context of acquired immunity (Senatorski and Hryniewiecka 2011). IgG binds antigens. This could be of importance, for example, in the direct neutralization of bacterial toxins or preventing virus penetration into cells. IgG3 activates the complement cascade, which supports the humoral arm of innate immunity (Quintana et al. 2008).

Another protein with reduced expression in the elderly belonged to the immunoglobulin superfamily leukocyte-associated immunoglobulin-like receptor 1, also known as CD305 or LAIR1. This is found on peripheral mononuclear cells, such as T cells, NK cells, and B cells. As an inhibitor, it can regulate immune responses and prevent against the lysis of cells identified as self (Meyaard 2008).

Autoimmune diseases

Abnormal cell activation underlies autoimmune diseases (Tsukumo and Yasutomo 2004). One reason for systemic lupus is the presence of antinuclear autoantibodies. Their production may result from overexposure to nucleosomal antigens (dsDNA, native DNA). One hypothesis involves the reduced clearance of nucleosomal antigens and their accumulation in old age concomitant with a reduced activity of deoxyribonuclease-1 in the blood, which is consistent with our research that demonstrated reduced expression of DNase-1 in the urine proteome of the elderly (Lachmann 2003). This also correlates with an increased incidence of systemic lupus erythematosus after 55 years of age.

Proteins with multidirectional effects

Another protein with reduced expression in the elderly is osteopontin (OPN), also known as bone sialoprotein I (Wang and Denhardt 2008). OPN plays an important role in bone remodeling as it anchors osteoclasts to the bone mineral matrix, which promotes the development of osteoclasts' corrugated surfaces that generate a low pH environment and locally initiate bone resorption. OPN in urine prevents the formation of urinary tract deposits. OPN also has immunological functions. By binding to receptors on lymphocytes, α4β1, α9β1, and α9β4, it affects the adhesion, migration, and survival of cells. A diverse range of these receptors' expression is found on macrophages, neutrophils, dendritic cells, and T and B cells.

As an immunomodulating agent, OPN has chemotactic properties for enhancing the recruitment of inflammatory cells. In animal experiments, inflammation is observed with prolonged persistence of macrophage clusters. OPN acts as an adhesion protein, which promotes cell adhesion and healing processes. In addition, OPN mediates cell activation and cytokine production and is a regulator of cellular apoptosis. It has effects similar to IL-10, a Th2 cytokine that exacerbates Th1 cell responses, and indirectly leads to increased production of antibodies by B lymphocytes. Thus, OPN plays a role in the pathogenesis of autoimmune diseases (rheumatoid arthritis, multiple sclerosis, AIH), cancer (lung, breast, colorectal cancers), and inflammatory diseases (glomerulonephritides, interstitial renal diseases, plaque formation) (Tachibana et al. 2012). Reduced OPN expression in older people could intensify allergy and atopic phenomena, as well as mitigate the course of autoimmune diseases, which is consistent with clinical observations (Jamroz-Wiśniewska et al. 2007).

Attractin (dipeptidyl peptidase IV (DP-IV)/CD26-like molecule) is a serine protease that has four copies of an EGF-like repeat motif and a C-type lectin domain and is involved in multiple intercellular communications and apoptosis, and regulates the effects of other cytokines and signaling molecules, including RANTES (Oravecz et al. 1997). It also inhibits TGF-β, thus contributing to extracellular matrix turnover, fibrosis, wound healing and tissue repair (Waller and Nicholson 2001), cell growth and proliferation, and anti-inflammatory properties. TGF-β also plays an essential role in the pathogenesis of fibrotic diseases of the kidney, liver, and lung, as well as atherosclerosis, artery remodeling in hypertension and diabetes mellitus, rheumatoid arthritis, and cancer (Harris et al. 2007; Zegarska et al. 2006; Jain et al. 2000; Reinhold et al. 1997). Decreased DP-IV expression in the urinary proteome may explain the lack of TGF-β inhibition and indicate a role for DP-IV and TGF-β in the pathogenesis of aging and may explain its reduced ability to regulate the immune system as well as increased susceptibility to type 2 diabetes seen in the elderly (Laudes et al. 2010).

Secretogranin-1 (chromogranin B, secretoneurin) is a neuroendocrine protein involved in the formation of secretory granules in the pituitary, adrenal medulla, and hypothalamus (Natori et al. 1998). It also has chemotactic properties for monocytes, eosinophils (strength comparable to IL8), and endothelial cells. It is a known regulator of endothelial cell proliferation and is a potent angiogenic factor. What its reduced expression in the urine proteome in the elderly means is not known.

Conclusions

Characteristic changes in healthy elderly peoples' urine proteomes appear to reflect certain physiological processes associated with aging. We confirmed results from our previous experimental and clinical studies regarding the influence of the aging process on the activity of serum proteolytic enzymes and that the activities of individual proteolytic enzymes in the serum change with age (Paczek et al. 2008, 2009). Changes in the circulatory system and the immune system were particularly clearly represented. Detailed identification of these “protein trails” provides for creating a more “global picture” of the metabolic changes that occur in the elderly.

Our research may help to identify active and structural proteins whose levels change significantly during the natural process of aging. These results may help identify clinically important pathogenic pathways and build a multiscale model of human aging. This, in turn, could enable the development of intervention methods with the intent of slowing down the aging process.

Below is the link to the electronic supplementary material.

(XLS 2.63 mb)