Skip to main content

A proteomic evaluation of urinary changes associated with cardiopulmonary bypass

Abstract

Background

The urinary proteome of patients undergoing cardiopulmonary bypass (CPB) may provide important insights into systemic and renal changes associated with the procedure. Such information may ultimately provide a basis to differentiate changes or properties associated with the development of acute kidney injury. While mass spectrometry (MS) analysis offers the potential for in-depth compositional analysis it is often limited in coverage and relative quantitation capacity. The aim of this study was to develop a process flow for the preparation and comparison of the intraoperative urinary proteome.

Methods

Urines were collected from patients at the start of CPB and 1-h into CPB. Pooled samples (n = 5) from each time point were processed using a modified Filter Assisted Sample Preparation protocol. The resulting peptides were analyzed by 2D-LC–MS/MS and by 1D-LC–MS/MS SWATH (Sequential Window acquisition of All Theoretical fragment ion spectra).

Results

The 2D-LC–MS/MS analysis identified 1324 proteins in the two pools, of which 744 were quantifiable. The SWATH approach provided quantitation for 730 proteins, 552 of which overlapped with the common population from the 2D-IDA results. Intensity correlation filtering between the two methods gave 475 proteins for biological interpretation. Proteins displaying greater than threefold changes (>log2 1.59) at 1-hour CPB relative to the initiation of CPB (26 down-regulated and 22 up-regulated) were selected for further analysis. Up-regulated proteins were enriched in GO terms related to humoral immune response, predominantly innate immunity (C4b, lactotransferrin, protein S100-A8, cathelicidin, myeloperoxidase) and extracellular matrix reorganization (e.g. MMP-9).

Conclusions

This study describes a scheme for processing urine from patients undergoing CPB for mass spectrometry-based analysis. The introduction of SWATH into the workflow offers a sample and instrument sparing approach to obtaining consistent in-depth sample analysis. The design of the methodology is such that it can be readily applied to large numbers of clinical samples with the potential for automation. The results also suggest that activation of the innate immune responses occur during cardiac bypass surgery.

Background

Urine proteomic analysis is a powerful tool that may provide insights into protein changes that are predictive of disease development and to define underlying pathophysiologic processes. Urine provides a non-invasive, rich source of proteins that are potentially informative of systemic and renal processes. Indeed while there are unique aspects to the urine proteome, it has been shown to have approximately 61 % overlap with both the kidney and plasma proteomes [1]. These considerations have been the rationale for a number of proteomic studies of urinary proteins [25]. However, urine proteomic analysis is complicated by protein concentration variability, the presence of interfering salts, metabolites and other substances [6]. These considerations often require significant sample processing which can impact the reproducibility of quantitative analysis of urine [7].

Different urine proteomics approaches have been applied to clinical scenarios in the intensive care unit [5, 8], renal transplants [9, 10], and pre/post-operative cardiac surgery [2]. However, proteomic evaluation of intraoperative urines is particularly challenging because of the presence of unique inhibitory materials that interfere with chromatography. Only targeted SELDI TOF–MS analysis has been successfully applied to intra-operative cardiac surgery urine samples [3, 11, 12], however this approach offers extremely limited proteome coverage and does not provide specific protein identification.

Shotgun proteomics on complex biological samples using in-depth 2D-LC–MS/MS analysis provides maximum proteome coverage and protein identification. However, the stochastic nature of DDA/IDA peptide selection for analysis often results in variable identification coverage between analyses [1, 13, 14]. Sequential Window acquisition of All Theoretical fragment ion spectra Mass-Spectrometry (SWATH-MS) has been proposed as an alternative approach [15]. SWATH-MS offers accurate and reproducible protein quantification from complex samples using data-independent acquisition. A significant advantage of SWATH analysis is that it reports the relative abundance of selected proteins in different samples. SWATH-MS also provides a permanent digital file representing the continuous MS-measurable proteome of a sample that can subsequently be re-analyzed as new approaches and insights become available [16].

This study describes a workflow for the processing and analysis of urines from adult patients undergoing cardiac bypass surgery. A comparison of DDA/IDA and SWATH-based analysis is provided suggesting that the latter may represent a useful adjunct for the comparison of clinical samples.

Methods

Patients and samples

The study protocol was approved by the Health Research Ethics Board, University of Manitoba and all patients provided informed consent. Patients with stable renal function (n = 5) were selected for analysis from a larger prospective observational cohort of adult cardiac surgery patients using the following criteria: baseline estimated glomerular filtration rate (eGFR) ≥60 mL/min and post-operative serum creatinine rise <10 % from baseline. Serial urine samples were collected at the initiation of CPB and after 1-h on CPB. Urines were spun at 870 g for 6 min, and the supernatants stored at −80 °C for further analysis. The baseline clinical characteristics of the patient cohort are described in Table 1.

Table 1 Baseline patient characteristics

Sample preparation for mass-spectrometry analysis

Total urine protein was quantified using a total Protein Kit, Micro Pyrogallol Red Method (Sigma, USA). A pool of urines from 5 patients (40 µg protein/patient) was prepared for each time point. The pooled samples were processed using a modified FASP protocol [17]. Pooled urines were treated for 5 min with 0.1 M DTT at 100 °C. The cooled samples were brought to 8 M urea, 100 mM Tris pH 8.0 and diluted with an equal volume of the same solution. The samples were transferred to the indicated Amicon MWCO ultra-15 centrifugal filter (Millipore) (i.e. 3, 10 or 30 kDa cut-off) and centrifuged at 4000×g for 15 min. The retentate was then diluted to 15 ml with 8 M urea, 100 mM Tris pH 8.0 and centrifuged for 30–50 min at 4000×g to remove excess DTT. This was repeated once after which 2 ml of 50 mM iodoacetamide was added to the sample and incubated in the dark for 15 min. The sample was centrifuged for 15 min at 4000×g and 2 ml of 20 mM DTT was added to react with residual iodoacetamide. The sample was centrifuged for 15 min at 4000×g, after which the retentate was brought to 500 μl with 50 mM ammonium bicarbonate. Trypsin was added at 1:50 ratio based on the starting protein amount and the sample was incubated overnight at 28 °C with shaking. The digest was brought to 500 mM NaCl and gently mixed. The samples were centrifuged for 15 min at 4000×g. The MWCO filter was washed with 200 μl of water by centrifugation for 15 min at 4000×g. The peptides were eluted from the membrane by centrifugation (3000×g, 15 min) with 500 μl of 15 % acetonitrile, 0.1 % trifluoroactic acid in water. The eluted peptides were dried using a speedvac to remove the organic solvents used for elution from the MWCO membranes. The sample was dissolved in 500 μl of 0.5 % trifluoroacetic acid and desalted with a C18-SD extraction disc cartridge (Sigma, USA). The eluted peptides were dried and reconstituted in 40 μl of 0.1 % formic acid in water. An estimation of peptide quantity was made using a Nano-drop spectrophotometer (Thermo Scientific Nanodrop 2000) at 205 nm absorbance where 1.0 absorbance unit corresponded to a peptide concentration of ~1 μg/μl.

2D-HPLC–MS/MS data acquisition and peptide identification

The processed trypsin digested peptides were separated using two-dimensional liquid chromatographic method [18]. C-18 cleaned peptides obtained after modified FASP procedure were gradient fractionated on a C18 X-Terra column, (1 × 100 mm, 3 µm, 100 Å, Waters, Milford, MA). Both eluents A (water) and B (90 % acetonitrile) contained 20 mM ammonium formate buffer pH 10.0. A total of 30 fractions were collected using a gradient of 1–44 % of solvent B in 30 min, at a flow rate of 300 µl/min. The resulting fractions were concatenated to eight fractions, dried and dissolved in 50 µl of eluent A (0.1 % Formic Acid). One microgram from each concatenated fraction was injected into a splitless nano-flow Tempo LC system (Eksigent, Dublin, CA) via a PepMap100 trap column (0.3 × 5 mm, 5 µm, 100 Å, Dionex, Sunnyvale, CA) and a 100 μm × 150 mm analytical column packed with 3 μm Luna C18(2) (Phenomenex, Torrance, CA). Both eluents A (2 % acetonitrile in water) and B (98 % acetonitrile) contained 0.1 % formic acid as ion-pairing modifier. A 0.44 % acetonitrile per minute linear gradient (0–35 % B in 80 min, 500 nl/min) was used for peptide elution, followed by 5 min wash with 80 % B.

Samples were analyzed on a 5600 TripleTOF mass spectrometer (SCIEX, Canada) in standard MS/MS data dependent acquisition (DDA) mode also referred as information dependent acquisition (IDA). Each DDA/IDA cycle included a 250 ms MS scan (400–1600 m/z) and up to 20 MS/MS (100 ms each, 100–1600 m/z) for ions state of +2 to +5 and an intensity of minimum 300 cps. Selected ions and their isotopes were dynamically excluded from further fragmentation for 12 s.

For each pooled sample the resulting eight instrument-level WIFF files were converted to MGF format using the conversion tool bundled with Analyst. These were then sequentially concatenated into a single MGF file for each sample pool. Peptide identification and quantitation was done using X!tandem (cyclone 2012.10.01.1) [19]. Peptides were identified against a human database of the June 2015 release of SwissProt with the following search settings: constant modification C +57.021 (cysteine protection); variable modifications: M, W +15.995 (oxidation) or +31.989 (double oxidation); S, T, Y +79.966 (phosphorylation), N, Q +0.984 (deamidation); parent mass error: ±20 PPM, fragment mass error: 0.1 Da. Peptide identification Log10 expectation values were computed using a survival function approach [20, 21].

SWATH-MS acquisition

SWATH-MS, also referred to as data independent acquisition (DIA), runs were conducted using a 5600 TripleTOF mass-spectrometer, which operated in a looped product ion mode. This mode helps instrument to allow quadrupole resolution of 25 amu per mass selection. Precursor selection windows had an overlap of 1 Da with their adjacent selection window to ensure complete isotope coverage between SWATH blocks. Each cycle of SWATH analysis consisted of a 100 ms MS scan and a 100 ms MS/MS scan in 25 m/z blocks in 400–1250 m/z range, and a total of 34 SWATH blocks were collected for each scan. An accumulation time of 100 ms was used for each fragment ion scan and for the survey scans acquired at the beginning of each cycle, resulting in a total cycle time of 3.5 s. The collision energy for each window was determined on the basis of the collision energy for a 2+ ion centered in the respective window (equation: 0.0625 × m/z −3.5) with a collision energy spread of 15 eV. The mass spectrometer was always operated in high sensitivity mode. Peptides were separated using the same LC system configuration as for the second dimension DDA/IDA acquisition, using a 0.30 % acetonitrile per minute linear gradient (0–33 % B in 110 min, 500 nl/min), followed by a 5 min wash with 80 % B.

SWATH library and analysis

For SWATH library construction all non-modified peptides across a large dataset 2D DDA/IDA runs were concatenated into a non-redundant collection using methods described earlier [22]. We opted to include all peptides from a broader study of eight urine pool samples, regardless of their expectation values or number of peptides per protein detected in each run, to yield maximum SWATH depth. The peptide entries from each DDA/IDA run were retention-aligned against a system average value using the retention times extracted for the custom peptides spiked into each sample.

To enhance our potential SWATH extraction coverage, we developed a method to supplement this library with transitions from the comprehensive human SWATH library [23]. These “10K-library” transitions represented peptide fragment ion patterns detected in SWATH runs, which differ somewhat from their DDA/IDA counterparts due to differences in instrument fragmentation settings across the acquisition methods. This library also provided transition versions for multiple potential charge states of a given peptide, a feature that varies between 1D and 2D runs as a product of finite instrument acquisition time per MS/MS cycle.

This combined library was retention aligned against average retention values acquired from prior DDA/IDA runs under identical chromatographic settings. This library drove the SWATH analysis component of Peakview (1.2.0.4; AB Sciex 2012). The SWATH parameters were a 6 min extraction window (±3 min), a mass accuracy 50 PPM, a search depth of 100 peptides per protein, 6 transitions per peptide, and a 1 % false discovery rate. Extraction was performed on all eight SWATH runs of our larger study simultaneously to permit optimal “gap filling”, yielding a modest transition detection rate of ~11 % (21,205 of 193,112 targeted transitions) spanning 6050 peptides.

The utility of the custom-peptide [24] based retention alignments is illustrated in (Additional file 1: Table S1). The differences between the library target retention times and their actual extracted retention times in both of our SWATH runs are well within our SWATH retention extraction time window with no apparent systematic skewing of the retention differences, indicating good control over the chromatography across the SWATH runs and between the methods.

Results and discussion

Optimization of urine sample preparation for mass spectrometry

The aim of the present study was to develop methods for the reliable comparison of the urinary proteome during CPB. Preliminary studies indicated that urine samples collected during cardiac surgery contained substance(s) that interfered with processing and chromatographic separations of the peptides derived from these samples. However, samples collected post operatively from the same patient group did not appear to contain the chromatographic-interfering substances. Although unproven it was felt that the most likely cause of these interfering effects were due to lipophilic substances possibly propofol, an anaesthetic used during the surgery.

After examining several approaches we established that a modified version of the filter-aided sample preparation (FASP) procedure [17] provided a suitable approach. This procedure involves sample processing with subsequent concentration of the sample using a molecular weight cut off filter. The concentrated sample is then processed in situ and the resulting peptides are eluted for analysis. This approach allows for sample concentration and clean-up in a single reaction vessel, limiting protein loss and minimizing technical variability. The SDS and urea used in this procedure also offered a means of removing hydrophobic interfering substances which were present in the samples collected from patients during surgery. Preliminary experiments determined that samples prepared with either 3 or 10 kDa cut off filters proposed in the original FASP method [17] still contained compounds that altered the chromatographic properties of the recovered peptides (Fig. 1a). The retention times of peptides isolated under these conditions were reduced and there was considerable peak spreading which resulted in sub-optimal separations. Increasing the pore size of the filters to 30 kDa molecular weight cut off membrane eliminated this effect suggesting that the interfering materials were removed under these conditions (Fig. 1a). However, there was concern that the increased pore size could also result in the selective loss of lower molecular weight species in the urine samples. A comparative analysis of the 30 kDa membrane processed samples with those produced with 3 or 10 kDa membranes indicated that the 30 kDa sample provided a greater number of proteins in the 1–30 kDa range and overall more protein identifications than the other samples (Fig. 1b). While these results could not exclude the possibility of some protein loss using the 30 kDa molecular weight cut off membranes, it was apparent that the recovery and range of proteins was better with these membranes than with either the 3 or 10 kDa membranes. This increased recovery was demonstrated by the higher numbers of protein identifications and the numbers of proteins with molecular weights less than 30 kDa in samples processed using the 30 kDa MWCO filters (502 and 162, respectively) as compared with the 3 kDa (290 and 97) or 10 kDa (394 and 123) MWCO processed samples. The recent results of Berger et al. suggest that perhaps even larger pore PVDF based membranes may be used for sample preparation because of the efficiency of protein capture [25].

Fig. 1
figure 1

a The chromatographic properties of CPB urine samples processed using 3, 10 or 30 kDa MWCO filters in the modified-FASP method. The total ion chromatogram of each sample is given. Note the peak spreading and altered retention times in the 3 and 10 kDa MWCO membrane processed samples. b The predicted molecular weight and pI distribution of proteins identified from the same urine sample FASP processed with 3, 10 or 30 kDa molecular weight cut-off membranes. The values beside each panel indicate the number of proteins with molecular weights less than 30 kDa and the total number of proteins identified in sample

Differential protein analysis

Digests of the samples were analyzed by 2D-LC–MS and peptides were quantified by summing the intensities of all fragments in their associated CID spectra. Assignment of peptides into their source proteins was performed by X!tandem’s algorithm, and we chose to use its primary protein assignments only. Protein level intensities were the sum of the intensities of member peptides, expressed a log2 scale for simple differential analysis.

High confidence protein identifications were based on detection of at least two distinct peptides each with log10 expectation scores of −1.5 or lower, roughly equivalent to a protein-level log10 expectation score of −3.0 or lower. This yielded the identification of 1154 proteins at the start of CPB and 915 proteins at 1-h into CPB surgery. The total number of unique proteins observed was 1325, with 744 proteins identified in both samples (Additional file 2: Table S2). However, over a third of the population of the proteins identified (581) were seen exclusively at either the start of CPB or 1-h into CPB, rendering them ambiguous for quantitative analysis. This result highlights the challenge of comparative analysis of complex samples by mass spectrometry which arises as a result of the stochastic nature of instrument dependent data acquisition. Given the significant amount of time and effort committed to performing 2D DDA/IDA analysis and the variable complementarity of coverage between samples this was felt to be a suboptimal approach because such outcomes make it unclear if “missing” proteins in a sample were exclusively expressed in one set of the samples or if the peptides were simply not selected in one of the analyses. This motivated us to examine the feasibility of using a hyper reaction monitoring data independent approach [15].

As described in the Methods and Materials, we generated a large SWATH library by combining an existing human protein transition library [23] (henceforth called the “10K-library”) and the results of our larger-scale multiple 2D-LC–MS analysis of urines from cardiac bypass patients. Our approach was to replace our observed DDA/IDA transitions with their 10 K-library counterparts for the peptides common to both library modes (and across their multiple charge state versions), and retain transition patterns for peptides observed in our samples’ library but outside the 10K-library. The 10K-library’s unitless iRT retention scale [26] was linearly mapped into retention time using a simple regression, exhibiting an excellent correlation of R2 > 0.99. The combined library contained 18,628 peptides spanning 3620 proteins; of these 3429 peptides spanning 934 proteins were not part of the original 10K-library but were detected in our 2D analysis of urine samples (Additional file 3: Table S3).

SWATH analysis

We applied the same two-peptide rule for SWATH quantitation as we did for the 2D DDA/IDA analysis, giving 4755 peptides spanning 730 proteins. This result suggested a remaining 1295 proteins were potentially quantifiable in SWATH, but by only a single peptide each. However, as quantitation based on a single peptide is less reliable these peptides were excluded from subsequent analysis. The complete peptide extraction report is available for additional exploration (see Availability of data and materials).

A total of 730 SWATH-quantifiable proteins were identified in the two sample pools, 552 of these also had quantitation values in both of the 2D LC runs. In addition there were 131 proteins that could be quantified using SWATH but not using the DDA/IDA derived data. These “gap filling” readouts indicated that SWATH presented the potential to differentially access approximately 22 % of the proteins that the 2D IDA dataset saw as incomplete (Additional file 3: Table S3).

Combined DDA/IDA and SWATH quantitation approaches

The combined DDA/IDA and SWATH quantitation data were used to complement the quantitative analysis of each approach. The correlation between log2 SWATH and DDA/IDA protein expression values within a time point was moderate (R2 = 0.794 for start of CPB across 651 common proteins; R2 = 0.736 for 1-h into CPB across 584 common proteins), suggesting a general agreement between the readout technologies (Additional file 4: Table S4).

Applying a population filtering to these “acquisition technology-replicates” improved our quantitation reliability with a modest loss of total quantified proteins. By expressing the difference between the SWATH and DDA/IDA readouts (at the same time point) as normalized Z-scores, then discarding proteins with |Z| > 1.65 in either population, the resulting 475 common proteins (86 % of the 552 protein initial list) improved the cross-technology difference correlation to R2 = 0.549, while the correlation of difference values for the 77 proteins rejected by this filtering step is low (R2 = 0.01) (Fig. 2).

Fig. 2
figure 2

The use of protein level intensity correlations between 2D DDA/IDA and SWATH runs to improve quantitation reliability. Comparison of log2 protein intensities a START of CPB b 1-hour CPB. c The expression of differences between SWATH and DDA/IDA readouts (at the same time point) as normalized Z-scores identifies proteins that show discordance between the two methods. This provides a basis for discarding the proteins with |Z| > 1.65 in either population (as shown by the red dots), the resulting 475 common proteins (86 % of the 552 protein initial list; the blue dots). This improved the cross-technology difference correlation to R2 = 0.549. d The correlation of difference values for the 77 proteins rejected by this filtering step, shown is very low (R2 = 0.01)

These results were encouraging considering that SWATH runs consumed approximately one-sixth of the instrument acquisition time, with one-third amount less sample preparation time as their 2D-DDA/IDA counterparts and fewer steps to potentially compromise quantitation. Going forward, the SWATH ion library will be applied to broader AKI studies for the analysis of individual patient urine samples, where the subsequent “large-N” readout sets are amendable to more sophisticated statistical methods.

The supportive filtered population of 475 proteins formed the basis of our biological interpretation. In the absence of sampling statistics, we chose a fairly severe cut-off of threefold (log2 1.59) protein expression difference required in both DDA/IDA and SWATH readouts. This gave 22 up-regulating proteins and 26 down-regulating proteins 1-h into CPB versus the start of CPB (Table 2).

Table 2 Correlation filtered 2D DDA/IDA and SWATH protein differences of magnitude greater than 1.59 (threefold up regulating at CPB-1-hour) from both methods marked in (A) and less than −1.59 (threefold down regulating at CPB-1-hour) marked in (B)

Comparative analysis of the urinary proteome in patients undergoing CPB

The proteomic changes observed in these cardiac surgery patients were felt to reflect physiological responses to the trauma of the surgery and the associated bypass process, including ischemia reperfusion injury.

The majority of proteins that displayed altered expression after 1-h of CPB were localised to exosomes (Table 3). These findings are consistent with the observation that exosomes are derived from multivesicular bodies which fuse with renal tubular epithelial cell membranes before being secreted into the urine [27]. Urinary exosomes are typically isolated by differential centrifugation prior to proteomic analysis, and identification of exosomal proteins here suggests that these methods may identify low abundance proteins; although proteome overlap between whole urine and exosomes has been reported [28]. These results could suggest that the repertoire of exosomal proteins changed during CPB but the process of exosome generation continued.

Table 3 GO cellular components enriched in differentially regulated proteins

The up-regulated proteins were enriched for humoral immune responses (Table 4), which included predominantly innate features as outlined by the list of proteins involved (e.g. complement, lactotransferrin, myeloperoxisae and cathelicidin anti-microbial peptide). It is noteworthy that Chitinase-3-like protein 1 (CHI3L1, also called YKL-40) was also increased at 1-h post initiation of CPB. CHI3L1 is a member of the chilectin subfamily of the 18 glycosyl hydrolase family. This protein has retained the ability to bind heparin, collagen and chitin but it is enzymatically inactive [29]. It has been shown to have a number of immunostimulatory properties for the adaptive immune system as well as pro-inflammatory activities [29]. Increased levels of CHI3L1 have been associated with loss of renal function in diabetes, and renal or cardiac transplantation patients [3033]. Our present observations suggest that CHI3L1 levels also increase in patients undergoing cardiac bypass surgery potentially as a result of intra operative subclinical ischemia–reperfusion injury.

Table 4 GO processes enriched in upregulated proteins

These up-regulated proteins also contained a number of extracellular components, although these proteins appeared to derive more from vesicle or granule contents. ECM reorganization was a process characterized by the up-regulated proteins, including elevation of urinary MMP9. These findings are consistent with reports that demonstrate elevated urinary MMP7 and renal MMP2 expression in ischemia reperfusion injury of renal allografts [4, 34]. While there was one diabetic patient in the cohort, it is unlikely that the acute intraoperative elevation of urinary MMP9 observed over 1-h of CPB was related to undocumented diabetic nephropathy [35].

Significantly, the up-regulated proteins also reflected general activation/regulation of innate immunity (C4b, lactotransferrin, protein S100-A8, myeloperoxidase) which is consistent with rapid, generic responses to cardiac surgery-induced tissue injury. Interestingly, there were also changes suggesting that the adaptive immune response might be down-regulated following 1-h on CPB, a process not previously described. Taken together, these findings corroborate pre-clinical models of the immune responses to renal ischemia reperfusion injury, but further expand on them to demonstrate that repair and re-modeling processes may be initiated sooner than previously thought [36].

It should be emphasised that the criteria used for the selection of proteins with significant changes between the two time points was based on the requirement that the protein was identified in both SWATH and DDA/IDA analyses of the same sample, and the changes had to be of comparable magnitudes. These represent very stringent criteria to provide confidence regarding the presence and magnitude of the change using two independent analytical methods. This was not meant to imply that other quantitative changes predicted by observations made in only one of the approaches (DDA/IDA or SWATH) may not also be valid. However, the confidence of such changes is lower and would require validation using alternate methodologies.

Conclusions

Our results provide a method for analyzing the urinary proteomic changes in adult patients undergoing cardiac surgery. Unlike the original FASP method the current procedure does not use SDS. The method requires limited sample processing to remove interfering substances and generates mass spectrometry ready samples. The design is such that it can be readily applied to large numbers of clinical samples with the potential for automation. The introduction of SWATH into the workflow offers several benefits to the analysis. The most significant aspect is the ability to consistently report on the same panel of >700 proteins in each analysis. This allows for more confident comparisons between samples, while keeping both the sample quantity requirements and the mass spectrometer instrument time relatively low. Although SWATH does not currently offer the analytical depth provided by 2D-LC–MS/MS analyses, it may be possible to extend the level of quantitation through longer chromatographic run times or via pre-fractionation steps prior to SWATH-MS. There have been several comparisons to suggest that SWATH provides comparable quantitative analysis to other methods [22]. The approach outlined here is particularly relevant to the analysis of sequential urines from patients in various clinical conditions, including intraoperative samples.

Abbreviations

CPB:

cardiopulmonary bypass

MS:

mass-spectrometry

FASP:

filter assisted sample preparation

SWATH:

Sequential Window acquisition of All Theoretical fragment ion spectra

eGFR:

estimated glomerular filtration rate

MWCO:

molecular weight cut-off ultra-centrifugal filters

HPLC:

high-performance liquid chromatography

2D-LC–MS/MS:

two dimensional liquid chromatography mass-spectrometry

References

  1. Farrah T, Deutsch EW, Omenn GS, Sun Z, Watts JD, Yamamoto T, et al. State of the human proteome in 2013 as viewed through PeptideAtlas: comparing the kidney, urine, and plasma proteomes for the biology- and disease-driven Human Proteome Project. J Proteome Res. 2014;13(1):60–75.

    Article  CAS  PubMed  Google Scholar 

  2. Aregger F, Pilop C, Uehlinger DE, Brunisholz R, Carrel TP, Frey FJ, et al. Urinary proteomics before and after extracorporeal circulation in patients with and without acute kidney injury. J Thorac Cardiovasc Surg. 2010;139(3):692–700.

    Article  CAS  PubMed  Google Scholar 

  3. Ho J, Lucy M, Krokhin O, Hayglass K, Pascoe E, Darroch G, et al. Mass spectrometry-based proteomic analysis of urine in acute kidney injury following cardiopulmonary bypass: a nested case-control study. Am J Kidney Dis Off J Natl Kidney Found. 2009;53(4):584–95.

    Article  CAS  Google Scholar 

  4. Ho J, Rush DN, Krokhin O, Antonovici M, Gao A, Bestland J, et al. Elevated urinary matrix metalloproteinase-7 detects underlying renal allograft inflammation and injury. Transplantation. 2016;100(3):648–54.

    Article  CAS  PubMed  Google Scholar 

  5. Maddens B, Ghesquiere B, Vanholder R, Demon D, Vanmassenhove J, Gevaert K, et al. Chitinase-like proteins are candidate biomarkers for sepsis-induced acute kidney injury. Mol Cell Proteomics MCP. 2012;11(6):M111 013094.

    Article  PubMed  Google Scholar 

  6. Thongboonkerd V. Practical points in urinary proteomics. J Proteome Res. 2007;6(10):3881–90.

    Article  CAS  PubMed  Google Scholar 

  7. Loftheim H, Nguyen TD, Malerod H, Lundanes E, Asberg A, Reubsaet L. 2-D hydrophilic interaction liquid chromatography-RP separation in urinary proteomics–minimizing variability through improved downstream workflow compatibility. J Sep Sci. 2010;33(6–7):864–72.

    Article  CAS  PubMed  Google Scholar 

  8. Zhou H, Pisitkun T, Aponte A, Yuen PS, Hoffert JD, Yasuda H, et al. Exosomal Fetuin-A identified by proteomics: a novel urinary biomarker for detecting acute kidney injury. Kidney Int. 2006;70(10):1847–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Sigdel TK, Kaushal A, Gritsenko M, Norbeck AD, Qian WJ, Xiao W, et al. Shotgun proteomics identifies proteins specific for acute renal transplant rejection. Proteomics Clin Appl. 2010;4(1):32–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Sigdel TK, Salomonis N, Nicora CD, Ryu S, He J, Dinh V, et al. The identification of novel potential injury mechanisms and candidate biomarkers in renal allograft rejection by quantitative proteomics. Mol Cell Proteomics MCP. 2014;13(2):621–31.

    Article  CAS  PubMed  Google Scholar 

  11. Devarajan P, Krawczeski CD, Nguyen MT, Kathman T, Wang Z, Parikh CR. Proteomic identification of early biomarkers of acute kidney injury after cardiac surgery in children. Am J Kidney Dis Off J Natl Kidney Found. 2010;56(4):632–42.

    Article  CAS  Google Scholar 

  12. Nguyen MT, Dent CL, Ross GF, Harris N, Manning PB, Mitsnefes MM, et al. Urinary aprotinin as a predictor of acute kidney injury after cardiac surgery in children receiving aprotinin therapy. Pediatr Nephrol. 2008;23(8):1317–26.

    Article  PubMed  Google Scholar 

  13. Adachi J, Kumar C, Zhang Y, Olsen JV, Mann M. The human urinary proteome contains more than 1500 proteins, including a large proportion of membrane proteins. Genome Biol. 2006;7(9):R80.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Nagaraj N, Mann M. Quantitative analysis of the intra- and inter-individual variability of the normal urinary proteome. J Proteome Res. 2011;10(2):637–45.

    Article  CAS  PubMed  Google Scholar 

  15. Gillet LC, Navarro P, Tate S, Rost H, Selevsek N, Reiter L, et al. Targeted data extraction of the MS/MS spectra generated by data-independent acquisition: a new concept for consistent and accurate proteome analysis. Mol Cell Proteomics MCP. 2012;11(6):O111 016717.

    Article  PubMed  Google Scholar 

  16. Guo T, Kouvonen P, Koh CC, Gillet LC, Wolski WE, Rost HL, et al. Rapid mass spectrometric conversion of tissue biopsy samples into permanent quantitative digital proteome maps. Nat Med. 2015;21(4):407–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wisniewski JR, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods. 2009;6(5):359–62.

    Article  CAS  PubMed  Google Scholar 

  18. Dwivedi RC, Spicer V, Harder M, Antonovici M, Ens W, Standing KG, et al. Practical implementation of 2D HPLC scheme with accurate peptide retention prediction in both dimensions for high-throughput bottom-up proteomics. Anal Chem. 2008;80(18):7036–42.

    Article  CAS  PubMed  Google Scholar 

  19. Craig R, Cortens JP, Beavis RC. The use of proteotypic peptide libraries for protein identification. Rapid Commun Mass Spectrom RCM. 2005;19(13):1844–50.

    Article  CAS  PubMed  Google Scholar 

  20. Fenyo D, Beavis RC. A method for assessing the statistical significance of mass spectrometry-based protein identifications using general scoring schemes. Anal Chem. 2003;75(4):768–74.

    Article  PubMed  Google Scholar 

  21. Fenyo D, Eriksson J, Beavis R. Mass spectrometric protein identification using the global proteome machine. Methods Mol Biol. 2010;673:189–202.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. McQueen P, Spicer V, Schellenberg J, Krokhin O, Sparling R, Levin D, et al. Whole cell, label free protein quantitation with data independent acquisition: quantitation at the MS2 level. Proteomics. 2015;15(1):16–24.

    Article  CAS  PubMed  Google Scholar 

  23. Rosenberger G, Koh CC, Guo T, Rost HL, Kouvonen P, Collins BC, et al. A repository of assays to quantify 10,000 human proteins by SWATH-MS. Sci Data. 2014;1:140031.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Krokhin OV, Spicer V. Peptide retention standards and hydrophobicity indexes in reversed-phase high-performance liquid chromatography of peptides. Anal Chem. 2009;81(22):9522–30.

    Article  CAS  PubMed  Google Scholar 

  25. Berger ST, Ahmed S, Muntel J, Cuevas Polo N, Bachur R, Kentsis A, et al. MStern blotting-high throughput polyvinylidene fluoride (PVDF) membrane-based proteomic sample preparation for 96-well plates. Mol Cell Proteomics MCP. 2015;14(10):2814–23.

    Article  CAS  PubMed  Google Scholar 

  26. Escher C, Reiter L, MacLean B, Ossola R, Herzog F, Chilton J, et al. Using iRT, a normalized retention time for more targeted measurement of peptides. Proteomics. 2012;12(8):1111–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Pisitkun T, Shen RF, Knepper MA. Identification and proteomic profiling of exosomes in human urine. Proc Natl Acad Sci USA. 2004;101(36):13368–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Salih M, Demmers JA, Bezstarosti K, Leonhard WN, Losekoot M, van Kooten C, et al. Proteomics of urinary vesicles links plakins and complement to polycystic kidney disease. J Am Soc Nephrol. 2016. [Epub ahead of print].

  29. Di Rosa M, Distefano G, Zorena K, Malaguarnera L. Chitinases and immunity: ancestral molecules with new functions. Immunobiology. 2016;221(3):399–411.

    Article  PubMed  Google Scholar 

  30. De Loor J, Decruyenaere J, Demeyere K, Nuytinck L, Hoste EA, Meyer E. Urinary chitinase 3-like protein 1 for early diagnosis of acute kidney injury: a prospective cohort study in adult critically ill patients. Crit Care. 2016;20(1):38.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Malyszko J, Koc-Zorawska E, Malyszko J. YKL-40, a marker of cardiovascular disease and endothelial dysfunction, in kidney transplant recipients. Transplant Proc. 2014;46(8):2651–3.

    Article  CAS  PubMed  Google Scholar 

  32. Przybylowski P, Janik L, Wasilewski G, Nowak E, Kozlik P, Malyszko J. YKL-40, a novel marker of cardiovascular complications, is related to kidney function in heart transplant recipients. Transplant Proc. 2014;46(8):2860–3.

    Article  CAS  PubMed  Google Scholar 

  33. Zurawska-Plaksej E, Lugowska A, Hetmanczyk K, Knapik-Kordecka M, Adamiec R, Piwowar A. Proteins from the 18 glycosyl hydrolase family are associated with kidney dysfunction in patients with diabetes type 2. Biomarkers. 2015;20(1):52–7.

    Article  CAS  PubMed  Google Scholar 

  34. Wanga S, Ceron CS, Delgado C, Joshi SK, Spaulding K, Walker JP, et al. Two distinct isoforms of matrix metalloproteinase-2 are associated with human delayed kidney graft function. PLoS One. 2015;10(9):e0136276.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Li SY, Huang PH, Yang AH, Tarng DC, Yang WC, Lin CC, et al. Matrix metalloproteinase-9 deficiency attenuates diabetic nephropathy by modulation of podocyte functions and dedifferentiation. Kidney Int. 2014;86(2):358–69.

    Article  CAS  PubMed  Google Scholar 

  36. Rabb H, Griffin MD, McKay DB, Swaminathan S, Pickkers P, Rosner MH, et al. Inflammation in AKI: current understanding, key questions, and knowledge gaps. J Am Soc Nephrol JASN. 2016;27(2):371–9.

    Article  PubMed  Google Scholar 

  37. Jensen LJ, Kuhn M, Stark M, Chaffron S, Creevey C, Muller J, et al. STRING 8—a global view on proteins and their functional interactions in 630 organisms. Nucleic Acids Res. 2009;37(Database issue):D412–6.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

RD conducted the experiments, data analysis and helped in manuscript writing; MN helped conduct the experiments; NC collected the patient samples and completed the patient database; VS performed the bioinformatics and data processing; OK provided technical input on the mass-spectrometry and chromatography platform; CR and RA designed and implemented the prospective observational cohort of cardiopulmonary bypass patients; JH and JW were involved in the overall study design, data analysis, manuscript writing and obtaining funding. All authors read and approved the final manuscript.

Acknowledgements

We would like to thank Peyman Ezzati for his role in method development and for his contributions to the sample analysis and mass spectrometry. The authors would like to gratefully acknowledge funding for this study from the Kidney Foundation of Manitoba, the Canadian Institutes of Health Research and Research Manitoba.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The supplementary data for this work is stored at the University of California, San Diego’s MassIVE archive (massive.ucsd.edu) under index MSV000079693. The 2D DDA/IDA run’s MGF files CPBSTART.mgf.gz and CPB1HR.mgf.gz are concatenated conversions from the raw instrument format, for the start of CPB and 1-hour CPB time points, respectively. The SWATH peptide-level extraction report is stored under CPBSWATH-Peptides.txt, while CPB-system.txt contains an overall log2 protein expression matrix for both the 2D DDA/IDA and SWATH data.

Ethics approval and consent to participate

This study protocol was approved by the University of Manitoba health research ethics board (HS15221) and all patients provided informed consent.

Funding

This study was funded by support from the Kidney Foundation of Canada (JH, JAW), Canadian Institutes of Health Research (JH, JAW) and Research Manitoba (JH). JH is supported by a CIHR New Investigator Salary Award. The funding body did not have any role in the study design, collection, analysis, data interpretation or writing of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to John A. Wilkins.

Additional files

12014_2016_9118_MOESM1_ESM.xlsx

Additional file 1: Table S1. The use of custom retention alignment peptides to map peptide retention times from 2D DDA/IDA runs into the SWATH extraction library, and their corresponding SWATH extracted retention times.

Additional file 2: Table S2. The log(2) protein values for DDA/IDA at time-points START and CPB-1-hour.

12014_2016_9118_MOESM3_ESM.xlsx

Additional file 3: Table S3. (A) An overview of the peptide and protein identifications and resulting differential quantitation for the 2D-LC DDA/IDA and SWATH run datasets. (B) The impact of the application of a DDA/IDA and SWATH intensity correlation filtering to the differential analysis.

Additional file 4: Table S4. Correlation filtered 2D DDA/IDA and SWATH protein difference values.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Dwivedi, R.C., Navarrete, M., Choi, N. et al. A proteomic evaluation of urinary changes associated with cardiopulmonary bypass. Clin Proteom 13, 17 (2016). https://doi.org/10.1186/s12014-016-9118-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12014-016-9118-9

Keywords