|Home | About | Journals | Submit | Contact Us | Français|
As part of a clinical proteomics programme focused on diabetes and its complications, it was our goal to investigate the proteome of plasma in order to find improved candidate biomarkers to predict diabetic nephropathy.
Proteins derived from plasma from a cross-sectional cohort of 123 type 1 diabetic patients previously diagnosed as normoalbuminuric, microalbuminuric or macroalbuminuric were enriched with hexapeptide library beads and subsequently pooled within three groups. Proteins from the three groups were compared by online liquid chromatography and tandem mass spectrometry in three identical repetitions using isobaric mass tags (iTRAQ). The results were further analysed with ingenuity pathway analysis. Levels of apolipoprotein A1, A2, B, C3, E and J were analysed and validated by a multiplex immunoassay in 20 type 1 diabetic patients with macroalbuminuria and 10 with normoalbuminuria.
A total of 112 proteins were identified in at least two out of three replicates. The global protein ratios were further evaluated by ingenuity pathway analysis, resulting in the recognition of apolipoprotein A2, B, C3, D and E as key nodes in the top-rated network. The multiplex immunoassay confirmed the overall protein expression patterns observed by the iTRAQ analysis.
The candidate biomarkers discovered in this cross-sectional cohort may turn out to be progression biomarkers and might have several clinical applications in the treatment and monitoring of diabetic nephropathy; however, they need to be confirmed in a longitudinal cohort.
The online version of this article (doi:10.1007/s12014-010-9053-0) contains supplementary material, which is available to authorized users.
Approximately 8% of type 1 diabetic patients develop end-stage renal disease because of diabetic nephropathy 30 years after diagnosis . As a consequence, diabetic nephropathy is the most common cause of renal failure in the developed world [2, 3]. Diabetic nephropathy advances through a number of recognisable steps from subclinical disease to the first measurable stage of microalbuminuria (MIC), defined as persistent urinary albumin excretion rates (UAER) of 30–300 mg/24 h, to macroalbuminuria/diabetic nephropathy (DMN) with UAER > 300 mg/24 h. DMN is characterised by declining renal function and ultimately end-stage renal disease. Although positive effects on the development and progression of diabetic nephropathy through strict control of blood glucose , blood pressure  and, in particular, blockade of the renin–angiotensin system (6, 7) have been reported, it still has not been enough to prevent the high incidence of end-stage kidney damage caused by diabetes. The urinary albumin excretion rate is used to identify patients at risk, but only a fraction of all patients with microalbuminuria progress to DMN. New biomarkers that can help in identifying normo- or microalbuminuric patients at risk for progression to DMN are needed as early intervention could further curb progression to DMN and end-stage renal disease.
Access to plasma samples from a cross-sectional study at the Steno Diabetes Center of type 1 diabetic patients with various degrees of diabetic nephropathy formed the basis for this study. The same cohort was previously investigated by fractionation coupled with surface-enhanced laser desorption/ionisation time-of-flight MS (SELDI-TOF-MS), which resulted in the discovery of differential levels of apolipoprotein C1 and apolipoprotein A1 along with transthyretin and cystatin C .
In this study, the cohort was further analysed using iTRAQ labelling (Applied Biosystems Inc.) and LC-MS/MS. In order to identify low abundant proteins in the plasma, we subjected the samples to hexapeptide library bead-based fractionation before proteolysis, iTRAQ labelling and LC-MS/MS. The hexapeptide library beads have been reported to have the capacity to reduce the most abundant proteins in plasma whilst simultaneously concentrating the most dilute species . The library contains linear hexapeptides with millions of different ligands. The diversity and number of peptides should result in specific baits for a large number of proteins present in a complex proteome [10, 11]. Because the library has the same quantity of each ligand, hypothetically, the maximum amount of each protein bound is the same under overloading conditions [12, 13]. The iTRAQ labelling simplified the analysis and included protein identification information not present in previous studies. Evaluation of the iTRAQ results by ingenuity pathway analysis identified apolipoproteins A2, B, C3, D and E as central nodes in a network that was highly rated. For that reason, the apolipoprotein profile of type 1 diabetic patients were investigated further in a cohort of patient samples that were chosen based on the rate of loss of kidney function over a time period of 4 years. We compared the ten patients with the fastest rate of decline with those with the slowest rate of decline. Ten samples came from patients with DMN and a high loss of renal filtration based on the measurement of GFR (annual difference in ΔGFR > 5.2 ml min−1/1.73 m2, termed RD), and ten samples came from DMN patients with a relatively slow decline of renal filtration, ΔGFR < 3.3 (termed SlD). In addition, ten samples from normoalbuminuric patients were analysed. The samples were analysed by the cardiovascular disease (CVD) immunoplex assay from Merck (Whitehouse Station, NJ, USA) detecting apolipoprotein A1, A2, B, C3, E, H and J.
By exploring the plasma proteome using a different approach, we aimed to expand our list of candidate biomarkers for diabetic nephropathy. Biomarkers identifying MIC or DMN in this cross-sectional study can be evaluated as potential predictors of the development of diabetic nephropathy in longitudinal studies of diabetic patients.
The cohort was collected as a part of the GENEX Study in 2004 [14, 15], and a total of 123 samples from Caucasian type 1 diabetic patients examined and diagnosed at Steno Diabetes Center had the adequate volume for the MS analysis. The patients were divided into three groups, normoalbuminuric (N, n=42), MIC (n=40) and DMN (n=41), according to the UACR in at least two out of three consecutive urines collected in the outpatient clinic. The time between DMN diagnosis and sample collection for the GENEX cohort was variable. DMN patients were selected for the presence of DMN irrespective of its duration. All patients in the DMN group, except for three, were being treated to reduce their blood pressure, and consequently, their urinary albumin (U-albumin) at the time of the blood sample collection was not necessarily representative of the stage of the disease. At the time of plasma sampling, 26 out of 41 DMN patients still had high U-albumin levels, whilst 13 had intermediate levels and two had levels low enough to place them in the N group. Of the 15 patients whose levels were lower than at the point of diagnosis, five were put on a short break from their antihypertensive medication and three of them returned to DMN levels before continuing their medication. The two other patients had significant level increases of DMN values and support the original diagnosis. In a follow-up period 4 years after sampling, there were only three patients (two from MIC and one from N) who had progressed to diabetic nephropathy. The low progression rate may be explained by the stable disease state of this cohort at baseline without diabetic nephropathy despite long-standing diabetes.
The three diabetic groups were matched at group level by gender, age and duration of diabetes. Investigations and blood sampling were performed in the morning after an overnight fast and has been thoroughly described previously . Comparisons of clinical data between the three groups were made with Prism 5 software (GraphPad Software, CA, USA). Values of p were two-sided, and a p value <0.05 was considered statistically significant. The Kruskal–Wallis (KW) test was used for ordinal data (gender). The distributions of continuous data were tested for normality by the D’Agostino–Pearson omnibus normality test. If distributions were normal, data were compared with one-way ANOVA, and if distributions were not normal, data were compared with a KW test.
The study was approved by the local ethics committee and all patients gave their informed consent.
Patient plasma was fractioned with ProteoMiner beads as described previously . Plasma proteins within each group were pooled and 100 μg protein was aliquoted three times for triplicate studies. Each replicate was precipitated with six times the volume of cold acetone separately, yielding a total of nine individual samples. The samples were reduced and blocked as suggested by the manufacturers. Proteins were digested overnight with Trypsin (2%, w/w), and the resulting peptides were subsequently labelled using iTRAQ 4-plex. The three samples from the DMN group were labelled with iTRAQ 114, the MIC group samples iTRAQ 115 and the N group samples iTRAQ 116, and dried in a vacuum centrifuge. The samples were reconstituted in 500 μL SCX loading buffer (Applied Biosystems) and mixed 1:1:1 to give a final of three mixed samples (each containing one iTRAQ 114-labelled DMN sample, one 115 iTRAQ-labelled MIC sample and so forth). The pH of the samples was adjusted to 2.5–3 with HCl and fractioned on a SCX column with 50, 100, 150, 200, 250, 300, 350, 400, 500 mM and 1 M KCl. The pH values of the fractions were adjusted to 2–3 and cleaned on C18 Stagetips (Proxeon, Odense, Denmark) using 50 μL 100% ACN as wetting solution, 50 μL 0.1% TFA for equilibrating the tips and washing the sample after loading, and, finally, elution of the peptides with 50 μL 0.1% TFA and 50% ACN. Afterwards, the samples were completely dried.
Peptides from the SCX fractions were dissolved in 0.3 μL 100% formic acid and diluted to 5 μL in 0.05% TFA. The peptides were loaded onto an EASY-nano LC system (Proxeon). Peptides were loaded directly onto a 20-cm, 100-μm inner diameter and 360-μm outer diameter, ReproSil-Pur C18 AQ 3 μm (Dr. Maisch, Ammerbuch-Entringen, Germany) reversed phase capillary column. The peptides were eluted using a gradient from 100% phase A (0.1% formic acid) to 40% phase B (0.1% formic acid, 80% acetonitrile) over 100 min at 200 nL/min directly into a LTQ-Orbitrap XL mass spectrometer (Thermo Scientific, San Jose, CA, USA). The LTQ-Orbitrap XL was operated in a data-independent mode automatically switching between MS, MS/MS and higher-energy collision dissociation (HCD) using a threshold of 30,000 for ion selection. For each MS scan, the two most abundant precursor ions were selected for fragmentation using MS/MS (normalised collision energy, 35; activation time, 30 ms) to facilitate peptide identification and then subsequently HCD (normalised collision energy, 55) to obtain reliable iTRAQ quantitative data . Raw data were viewed in Xcalibur v2.1 (Thermo Scientific).
The data were processed using Proteome Discoverer, version 220.127.116.11. The files were subsequently submitted to an in-house MASCOT server (version 2.2.05, Matrix Science Ltd., London, UK) for database searching through the Proteome Discoverer programme. The data were searched against the human sequence database (NCBInr 20100602 (11186807 sequences; 3815639892 residues); taxonomy: Homo sapiens (human, 232199 sequences)). The search was performed choosing trypsin as specific enzyme. A maximum of one missed cleavage was allowed. Oxidation (M), iTRAQ 4plex (K) and iTRAQ 4plex (N-term) were chosen as variable modifications. The data obtained on the LTQ-orbitrap XL were searched with a peptide mass tolerance of 10 ppm and a fragment mass tolerance of 0.8 Da. Protein ratios were normalised using the overall median ratio for all the peptides in the sample for each separate ratio in every individual experiment. The ratio for a given protein was calculated by taking the average of all the peptides ratios that identified the protein. The final list of protein ratios was an average of the protein ratios of the three experiments and consisted only of proteins discovered in at least two of the three experiments.
The iTRAQ results were evaluated using ingenuity pathway analysis (IPA; Ingenuity Systems, Mountain View, CA; www.ingenuity.com). This software analyses protein expression data in the context of known biological response and regulatory networks as well as other higher-order response pathways. IPA constructs hypothetical protein interaction clusters on the foundation of a regularly updated Ingenuity Pathways Knowledge Base, which is a very large database that consists of millions of individual relationships between proteins gathered from the biologic literature . The data set that contained the differentially expressed proteins identified in the iTRAQ experiment was converted by IPA to ‘fold change’ and uploaded into IPA. No expression value cutoff was selected, and we were interested in looking at both up- and downregulated proteins. Hypothetical networks were generated from these proteins and other proteins from the database that were needed to fill out a protein cluster. Network generation was optimised for the inclusion of as many proteins from the input expression profile as possible and aimed for highly connected networks.
IPA computes a score for each possible network according to the fit of that network to the input proteins. The score is calculated as the negative base-10 logarithm of the p value that indicates the probability of the input proteins in a given network being found together as a result of random chance. Consequently, scores of 2 or higher have at least a 99% confidence of not being generated by random chance alone .
Concentrations of Apo A1, A2, B, C3, E, H and J were simultaneously evaluated using a multiplex bead-based sandwich immunoassay kits (Novagen WideScreen ™ Human CVD Panel 1, Merck Chemicals Ltd. Nottingham, UK). Assays were performed following the manufacturer’s instructions. Briefly, seven distinct sets of fluorescently dyed beads loaded with capture monoclonal antibodies specific for each apolipoprotein to be tested were used. Plasma samples diluted 1:2,500 with dilution buffer (30 μL/well) or standards (30 μL/well) along with 10 μL of biotinylated antigens in blocking buffer were incubated with 10 μl of pre-mixed bead sets into the wells of a pre-wet 96-well microtitre plate. After incubation and washing, 20 μL of fluorescent detection antibody mixture (Streptavidin PE) was added for 30 min and the samples were washed and resuspended in assay buffer. The formation of different sandwich immunocomplexes on distinct bead sets was measured and quantified using the Bio-Plex Protein Array System (Bio-Rad Laboratories). A 50 μL volume was sampled from each well and the fluorescent signal of a minimum of 100 beads per region (chemokine/apolipoprotein) was evaluated and recorded. Values presenting a coefficient of variation beyond 15% were discarded before the final data analysis.
The clinical data were evaluated by the Prism software. Bio-plex data were analysed using the Bio-Plex Manager software version 3.0 (Bio-Rad Laboratories). Standard levels between 80% and 120% of the expected values were considered to be accurate and were used. Values were expressed as nanograms per millilitre and presented as mean (SD).
The clinical data of the GENEX cohort are summarised in Table 1. Patient groups in the cross-sectional cohort were matched with respect to gender, duration of diabetes and body mass index (BMI), but differed slightly by age (p=0.01). The DMN group had significantly lower estimated glomerular filtration rates (eGFR) compared to the other groups (p<0.0001), whilst equivalent eGFR values were observed in the N and MIC groups. There were no significant differences in systolic blood pressure or diastolic blood pressure for comparisons of all groups or for levels of serum cholesterol between the groups. Haemoglobin A1c (HbA1c) differed significantly for an all group comparison (p=0.02); this was caused by a significant difference between the N and DMN groups.
Protein identifications were initially accepted when based on peptide identifications with Mascot scores of ≥25 and a false discovery rate of <1%. This resulted in the identification of 3,365, 3,667 and 3,495 peptides, respectively, from the iTRAQ triplicate experiments. In total, 130, 133 and 122 proteins were identified in experiments 1, 2 and 3, respectively. The relative quantitation was measured as two pairwise ratios: the MIC group against the DMN group (iTRAQ115/iTRAQ114) and the N group against the DMN group (iTRAQ116/iTRAQ114). A total of 112 proteins were identified in at least two out of three iTRAQ runs. The molecular identity and abundance ratios of the 112 proteins are listed in Electronic Supplementary Materials. Examples of MS/MS spectra of peptides from apo E and A1 are presented in Fig. 1.
Figures 2 and and3a3a shows one of the largest protein clusters that were generated by the pathway analysis of the protein ratios between N/DMN. This network consists of a cluster of 35 proteins, of which 15 were not included in our list. The network was suggested to be involved in cardiovascular disease, genetic disorders and metabolic diseases, and had a score of 36. Apo A2, B, C3, D and E were all included in this protein network, and the results inspired us to look further into these apolipoproteins in plasma.
For further analysis of the aforementioned apolipoproteins and validation of the iTRAQ results, 30 plasma samples were selected for quantitative measurements of apolipoproteins by a multiplex immunoassay. The patient samples were chosen based on the rate of loss of kidney function over a time period of 4 years. Ten samples came from patients with DMN and a high loss of renal filtration based on the measurement of GFR (annual difference in ΔGFR > 5.2, termed RD), and ten samples came from DMN patients with a relatively slow decline of renal filtration, ΔGFR < 3.3 (termed SlD). In addition, ten samples from normoalbuminuric patients were analysed.
The clinical data for the 30 patients chosen for the multiplex immunoassay is listed in Table 2.
There was no significant difference in GFR measured in 2004 between the RD and SlD groups. The distribution of gender, systolic and diastolic blood pressure, BMI, HbA1c, cholesterol and triglyceride values were not significant different across all three groups. The RD and SlD groups varied in the decline of renal function (ΔGFR) over a subsequent 4-year period (p=0.0002). At follow-up, the three groups varied slightly in age and duration of diabetes. This was caused by one young patient (aged 27) in the RD group and one patient with only 2 years of diabetes also in the RD group.
Plasma concentrations of apo A1, A2, B, C3, E, H and J were detected in all the patients. Table 3 shows the concentrations of the measured apolipoproteins in the N, SlD and RD groups. The concentrations of the apolipoproteins did not differ between the three groups. Apo A1 concentrations correlated with systolic blood pressure (Spearman ρ coefficient = −0.45, p=0.014), Apo A2 with triglycerides (Spearman ρ coefficient = 0.41, p=0.0245), Apo B with age (Spearman ρ coefficient = 0.3642, p=0.0478) and Apo C3 with cholesterol (Spearman ρ coefficient = 0.4361, p=0.0160). There was no correlation between ΔGFR and concentrations of apolipoproteins in the SlD and RD groups, nor any between the GFR in 2004 and apolipoproteins in all three groups.
In order to compare the concentrations of the apolipoproteins from the multiplex immunoassay with the ratios of the apolipoproteins from the iTRAQ experiment, the average of the specified protein from the N group was divided by the average of the two macroalbuminuric (N/SlD + RD) groups and the SlD (N/SlD) and RD (N/RD) groups separately (Fig. 4). The ratios of the apolipoprotein in the multiplex immunoassay had, overall, the same pattern as first revealed by the iTRAQ analysis, except for apo C3 and apo J where the reversed relationship was seen between iTRAQ and the multiplex immunoassay.
Identification of biomarkers in the blood of patients with high risk of developing diabetic nephropathy may have several clinical applications, including more aggressive control of glucose or arterial blood pressure, earlier onset of renoprotective treatment and a more tailored management of the disease.
Although several studies of urine, plasma and serum from T2D patients have been conducted [18–21], the plasma of T1D patients has not been fully examined. In this report, we present data to show that plasma obtained from three different stages of diabetic nephropathy in T1D patients contains candidate biomarkers for diabetic nephropathy.
Plasma samples from a cross-sectional study conducted at Steno Diabetes Center in 2004 formed the foundation of the investigations presented here. Previously, the plasma samples were fractioned with ProteoMiner beads and analysed on SELDI-TOF-MS . The SELDI-TOF-MS approach relied on a side-by-side comparison of spectra generated from each patient and revealed a number of protein peaks significantly different between patients divided into groups according to albuminuria. The disadvantage of this methodology was the time-consuming job of identifying the various protein peaks in the spectra by immunoprecipitation, sequencing and Western blots. The combination of isobaric labelling using iTRAQ and LC-MS/MS offered a way to overcome this situation. Peptides are identified in the MS/MS mode, whilst quantitation is obtained from the iTRAQ reporter ion using higher-energy collision dissociation (HCD) and collision-induced dissociation (CID). A total of 112 proteins were identified in at least two out of three repetitions of the experiments. The protein number is relatively low compared to other published studies, such as those reported for a head and neck cancer biomarker discovery study and protein expression profiling of plasma from patients with placental diseases where 811 and 166 proteins were identified, respectively [22, 23]. An explanation could be the initial treatment with ProteoMiner beads; previously, Pernemalm et al.  evaluated the ProteoMiner technique on crude plasma and identified only 150 unique proteins by direct LC-MS/MS and GeLC-MS/MS with a start plasma volume comparable to ours. Additionally, a recent study demonstrated that iTRAQ isobaric tags significantly increased the average ion charge state of peptides, resulting in a considerable reduction in the number of peptides identified from a given sample when using CID and HCD peptide fragmentation .
Some of the more interesting proteins identified in this study are the lipoproteins where the majority are present in higher levels in MIC and N compared to DMN. Evaluation of the iTRAQ data by ingenuity pathway analysis identified apo A1, A2, B, C3, E and J (Clusterin) to contribute in the top-rated network (Fig. 2). Diabetes has been shown to be associated with the amounts of different classes of apolipoproteins, and accordingly, the apolipoproteins might contribute to the development of diabetic nephropathy and cardiovascular disease . Apo B, Apo C1 and C3 are according to Hirano et al.  increased in patients with diabetic nephropathy compared to healthy individuals. We observed the same tendency in our results for apo C3 where higher levels of this protein were found in the groups with macroalbumninuria. Several reports have investigated different isoforms of apo E caused by three different alleles of the gene coding for apo E [28–30], and it has been suggested that the apo E4 isoform protects against diabetic nephropathy whilst the apo E2 isoform promotes the disease [31, 32]. We did not see any differences in the overall level of apo E in the different patient groups, suggesting that diagnostic potential of the protein potentially should be found in the different isoforms in contrast to the generally protein level. Apo A1 levels were lower in patients with macroalbuminuria in our study, although not statistically significant. Tolonen et al.  reported that low apo A1 levels predicted the progression of patients with macroalbuminuria to end-stage renal disease, and it might be reasonable to speculate that low levels of this protein in some way is associated to the development and progression of diabetic nephropathy. Whether the abnormal levels of the apolipoproteins are causing the development of diabetic nephropathy and therefore could be used as potential predictive markers or are just a consequence of the disease is difficult to say based on these cross-sectional studies presented here. The issue would ideally be addressed through further investigation in a longitudinal study.
To verify our results obtained with the iTRAQ investigation and ingenuity pathway analysis, we examined apo A1, A2, C3, E, H and J on an orthogonal platform in a series of patients with macroalbuminuria that had either a fast decline in kidney function (GFR, the RD group) and patients with macroalbuminuria and a slower decline in GFR (the SlD group) over 4 years. The ratios obtained in the iTRAQ analysis and the multiplex immuno-analysis appear quite similar. Differences could be explained by the initial treatment with ProteoMiner beads of the samples for the iTRAQ experiment, whereas the multiplex immunoassay was conducted on untreated plasma. Another explanation could be the difference in the number of plasma samples that constitutes the pool, i.e. the fact that the iTRAQ analysis is conducted in pools of plasma samples from a larger group of patients and therefore individuals with very high or low levels of protein would be diluted. In the multiplex immunoassay experiment, one patient had very high levels of apo A1, A2, C3 and H, and this resulted in a general higher average of these proteins. Another explanation could be potentially different isoforms of the apolipoproteins, with differences in which isoforms can be detected by the method employed. A recent paper by Borges et al.  reviews the microheterogenity of the human proteome. Evidently, polymorphisms, transcript variants and posttranslational modifications give rise to structurally similar but functionally distinct proteins whose detection is platform-dependent.
Originally, we hypothesised that apolipoprotein levels would be different between patients with macroalbuminuria that had a fast decline in kidney function (GFR, the RD group) and patients with macroalbuminuria and a lower fall in GFR (the SlD group). This was not observed in this study, possibly because of the relatively low number of samples in each group.
The candidate biomarkers discovered in this cross-sectional cohort may turn out to be progression biomarkers, but they need to be confirmed in a longitudinal cohort. The patients in the cohort all represent a stable disease state due to their long duration, and the proteins discovered in this study are indeed indicative of their current state, but temporal information about the changes needs to be investigated. The nature of the current cohort with a spread in U-albumin in the DMN and MIC group brings robustness and persistence to the actually obtained candidate biomarkers and reflects everyday clinical setting.
In this study, we investigated the plasma protein profiles from T1D patients with various degrees of albuminuria using hexapeptide library fractionation, iTRAQ labelling and LC-MS/MS, with the result of increasing our list of candidate biomarkers for diabetic nephropathy. Future research is aimed at analysing larger groups of samples. After establishing a model based on proteomic patterns for the separation of the three groups, the model will be validated in longitudinal sample sets to determine the diagnostic and predictive value of the obtained protein profiles. It is our goal to discover biomarkers able to predict, detect and monitor the progression of DMN and the effect of reno-protective intervention in diabetic patients.
Below is the link to the electronic supplementary material.
(DOC 182 kb)
This work was supported by funding from the Juvenile Diabetes Research Foundation, Paul and Erna Sehested Hansen Foundation, Steno Diabetes Center and Novo Nordisk A/S. We thank Peter Iversen for technical assistance with the Luminex immunoassays.
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.