|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular & Cellular Proteomics 7:290-298, 2008.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
We hypothesized that proteomics analysis of urine should yield a panel of biomarker peptides useful as additional tools for the diagnosis and monitoring of CAD. Furthermore we aimed to obtain sequences of biomarkers of the CAD signature panel to gain insight into pathogenetic mechanisms and facilitate comparison with currently used biochemical markers. Capillary electrophoresis on-line coupled to electrospray ionization-time-of-flight mass spectrometry (CE-ESI-TOF MS) seems ideally suited for this purpose because of its non-invasive nature, high resolution, and amenability for future adaptation to clinical laboratory analysis (5).
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Vascular stiffness was assessed by two methods. First, pulse contour analysis of the diastolic pressure decay was used to estimate large (C1) and small artery compliance (C2; HDI/Pulse Wave CR2000, HDI Inc., Eagan, MN) based on a three-element Windkessel model (10). Second, the augmentation index of the central aorta was derived from the radial pulse waveform using a generalized transfer function (SphygmoCor pulse wave analysis system, AtCor Medical, West Ryde, New South Wales, Australia) (11). The augmentation index was calculated from the ratio of the pulse pressure at the second systolic peak to that at the first systolic peak.
The study was approved by the West Glasgow Ethics Committee, and all subjects gave written informed consent. This study was designed according to current guidelines for studies on clinical proteomics (12) and the minimum information about proteomics experiments (MIAPE) (13).
To exclude the effect of medication, 17 paired urine samples from age- and sex-matched patients with hypertension and type 2 diabetes, but without albuminuria, before and 12 weeks after commencing treatment with the angiotensin-converting enzyme inhibitor ramipril (5–10 mg once daily) were evaluated. To rule out center specific bias, samples from 233 new appointees at the University of Hannover who were free of self-reported illness were also analyzed. Detailed characteristics of all patients and controls are shown in Table I with additional data on 233 healthy university recruits and 18 ramipril patients shown in Table II.
|
|
Urine Sample Preparation and CE-MS Analysis—
After collection, all the spot urine samples were frozen at –80 °C until analysis. For proteomics analysis samples were prepared as described previously (15, 16). CE-MS analysis was performed as described previously (15–17) using a P/ACE MDQ capillary electrophoresis system (Beckman Coulter, Fullerton, CA) on-line coupled to a Micro-TOF MS instrument (Bruker Daltonics, Bremen, Germany). The performance of the sample preparation procedure as well as the analytical performance of the instrumental setup was evaluated (5, 18). The average recovery of the sample preparation procedure is
85% with a detection limit of
1 fmol. The monoisotopic mass signals could be resolved for z
6. The mass accuracy of the CE-TOF-MS method was determined to be <25 ppm for monoisotopic resolution and <100 ppm for unresolved peaks (z > 6). The precision of the analytical method was determined by assessing (a) the reproducibility achieved for repeated measurement of the same aliquot and (b) the reproducibility achieved for repeated preparation and measurement of the same urine sample (5, 18). The 200 most abundant polypeptides were detected with a rate of 98%. The performance of the analytical system over time was assessed with consecutive measurements of the same aliquot over a period of 24 h. No significant loss of polypeptides was observed implying the stability of the CE-MS setup, the postpreparative stability of the urine samples at 4 °C, and their resistance to e.g. oxidizing processes or precipitation (5, 18).
Data Processing and Cluster Analysis—
Data processing and cluster analysis were performed as described previously (15, 19). Only signals observed in a minimum of three consecutive spectra with a minimum signal-to-noise ratio of 4 were considered. Mass spectral ion peaks representing identical molecules at different charge states were deconvoluted into single masses using either the distance between resolved isotope peaks of the ion or according to conjugated signals for unresolved isotope peaks (MosaiquesVisu software (16, 18)). In addition, migration time and ion signal intensity (amplitude) were normalized using internal polypeptide standards (15). The resulting peak list characterizes each polypeptide by its molecular mass (kDa), normalized migration time (min), and normalized signal intensity. All detected polypeptides were deposited, matched, and annotated in a Microsoft SQL (structured query language) database, allowing further analysis and comparison of multiple samples (patient groups). Polypeptides within different samples were considered identical if the mass deviation was less than 100 ppm and the migration time deviation was less than 3%. CE-MS data of all individual samples can be accessed in the supplemental table.
Definition of Biomarkers and Sample Classification—
For biomarker panel definition, we used polypeptides that were found in more than 75% of the urine samples in at least one of the different groups of the training set (e.g. CAD or healthy controls). Polypeptides fulfilling this criterion were further evaluated using receiver operating characteristic (ROC) statistics (20). The amplitude distribution of the CE-MS data of polypeptides present in the samples was used as the ROC variable, and the affiliation to a diagnostic group (i.e. CAD or healthy control) was used as the classification variable. The obtained area under the ROC curve (AUC) value of the analysis of a given polypeptide was interpreted as a measure of its discriminatory potential. An initial list of potential marker candidates was further refined using the Mann-Whitney test with p
0.05 as the significance level. Model establishment and sample classification were performed by using a linear classifier algorithm according to F =
i ci log Ai with F as classification factor, ci as classification coefficient, and Ai as amplitude of the CE-MS signal of the marker i. The algorithm generates a classification model based on polypeptides that are best suited to discriminate between two defined sample groups. Models consist of fewer biomarkers than samples to avoid overfitting of models. The probability to have CAD at a given classification factor F taking into account the related probability for a negative diagnosis was calculated according to Equation 1.
![]() |
Pattern Composition—
For the first phase of the study a training set was established. The training set consisted of 50 urine samples from randomly selected subjects, 30 CAD patients, and 20 control subjects, respectively. The first step of biomarker selection led to a set of 187 potential CAD-specific polypeptides.
In a second step, these preselected polypeptides were compared with 233 urine samples from healthy volunteers from different centers to eliminate polypeptides that may show center-specific bias. To exclude the effect of medication on constituting markers, an additional control group of patients before (n = 15) and after (n = 17) 12-week treatment with the angiotensin-converting enzyme ramipril was used to refine the selected polypeptides. Polypeptides that showed up/down-regulation of CE-MS signal intensity in direct comparison of both groups and in addition a uniform behavior in pairwise comparison in the majority of patients were considered as medication artifacts and eliminated.
Subsequently the established pattern of 15 polypeptides was evaluated in a blinded assessment of 59 urine samples: 47 samples from patients with CAD and 12 samples from healthy controls. All samples were examined using the CAD panel. Seventy-six urine samples from follow-up examination were also evaluated using the CAD panel.
Sequencing of Polypeptides—
Peptide sequencing was performed using an LTQ-OrbitrapTM hybrid mass spectrometer (Thermo Fisher Scientific, Bremen, Germany) equipped with a Dionex Ultimate 3000 nanoflow system and a nanoelectrospray ion source. Peptide separation took place on a 5-µm C18 nanocolumn (NanoSeparations, Nieuwkoop, Netherlands) in a precolumn setup using a flow rate of 5 µl/min followed by a flow of 250 nl/min and a linear gradient (60 min) from 2 to 50% MeCN in H2O (0.1% formic acid). The mass spectrometer was operated in data-dependent mode to automatically switch between MS and MS/MS acquisition. Survey full-scan MS spectra (from m/z 300 to 2000) were acquired in the Orbitrap with resolution R = 60,000 at m/z 400 (target value of 500,000 charges in the linear ion trap). The most intense ions were sequentially isolated for fragmentation in the linear ion trap using collisionally induced dissociation and the detection took place either in the linear ion trap (parallel mode; target value 10,000) or in the Orbitrap (target value of 500,000). Orbitrap MS/MS were acquired with resolution R = 15,000 at m/z 400. General mass spectrometric conditions were: electrospray voltage, 1.6 kV; no sheath and auxiliary gas flow; ion transfer tube temperature, 225 °C; collision gas pressure, 1.3 millitorrs; normalized collision energy, 32% for MS/MS. The ion selection threshold was 500 counts for MS/MS.
Further analysis was performed using instruments with electron transfer dissociation (ETD) capability (21–23). Upon arrival, samples were resuspended (50 µl of 100 mM acetic acid) and bomb-loaded onto a 360 x 75-µm microcapillary precolumn that was connected to a 360 x 50-µm analytical column with a
1-µm tip pulled with a laser puller (both columns were packed in-house with
5–8 cm of C18 resin). Peptides were separated by nano-reversed phase HPLC (Agilent 1100; flow split by a tee to
100 nl/min) and introduced into either an ETD-enabled LTQ quadrupole linear ion trap (Thermo Fisher Scientific, San Jose, CA) or LTQ-Orbitrap (Thermo Fisher Scientific, Bremen, Germany) mass spectrometer via nano-ESI. Samples analyzed using the LTQ were run in a data-dependent manner with the five most abundant species subjected to both ETD and collision-induced dissociation fragmentation (in separate alternating scans). Ion trap instrumental parameters were used as described recently (24). LTQ-Orbitrap analyses were performed according to parameters presented previously (25). Targeted analyses were also performed with the LTQ-Orbitrap where target m/z values of interest, which were observed in CE-MS analyses but were not characterized previously, were inspected. Spectral data were converted into .dta files (DTAs) using Bioworks Browser and MakeDTA (a gift from Don Hunt) for collision-induced dissociation and ETD, respectively. Collision-induced dissociation and ETD DTAs were subsequently made into separate batches of
2000 files each. Alternatively LC-MALDI-TOF and LC-Q-TOF peptide sequencing was performed as described in detail elsewhere (26).
Identification of Peptide Sequences—
MS/MS data were submitted to MASCOT (Matrix Science) for a search against human entries in the Mass Spectrometry Protein Sequence Database (MSDB). Accepted parent ion mass deviation was 50 ppm; accepted fragment ion mass deviation was 500 ppm. All search results with a MASCOT peptide score better than 20, depending on the ion coverage as related to the main spectra features, were accepted. Data files were also searched against the National Center for Biotechnology Information (NCBI) human non-redundant database using the Open Mass Spectrometry Search Algorithm (OMSSA) with an E-value cutoff of 0.01. The number of basic and neutral polar amino acids of obtained peptide sequences was utilized to correlate peptide sequencing data to CE-MS data as described earlier (26).
Statistical Analyses—
All of the statistical analyses for patient characteristics and clinical data were performed using the Minitab (Minitab for Windows 12.21, Minitab Inc., State College, PA) software package. The Kolmogorov-Smirnov test was used to test for normal distribution of the data. Data are expressed as mean ± S.D. if normally distributed or median (interquartile range) if their distribution was not normal. Categorical variables are presented as frequency counts, and intergroup comparisons were analyzed by
2 test for smoking, gender, and medication with statins. For continuous variables, differences between the groups were evaluated using unpaired Student's t test or Mann-Whitney U test for variables that were normally distributed and those that were not normally distributed, respectively.
Estimates of sensitivity and specificity were calculated based on tabulating the number of correctly classified samples. Confidence intervals (95% CI) were calculated in MedCalc (MedCalc for Windows 8.1.1.0, Medcalc Software, Mariakerke, Belgium). The ROC plot was obtained by plotting all sensitivity values (true positive fraction) on the y axis against their equivalent (1 – specificity) values (false positive fraction) for all available thresholds on the x axis (MedCalc Software). The AUC was evaluated because it provides a single measure of overall accuracy that is not dependent upon a particular threshold (20).
| RESULTS |
|---|
|
|
|---|
At base-line evaluation, midstream urine specimens were available from 86 of 88 patients; nine of 86 urine samples did not fulfill quality control criteria (15). Of the remaining 77 samples, 30 were used to establish a training set, and 47 were evaluated in a blinded assessment (Table II). The 77 patients who entered analysis did not differ in clinical characteristics from the 88 originally included in the study. In addition to healthy controls from the same population (Glasgow, UK), a group of healthy controls from another population (Hannover, Germany) and patients before and after treatment with ramipril (Nürnberg, Germany) were used to rule out center-specific bias and effect of medication on constituting markers, respectively (Table II). At follow-up evaluation only five of the 82 urine samples failed quality control criteria.
Samples were analyzed by CE-MS, and biomarker peptide patterns distinguishing CAD from healthy controls could be established (Fig. 1 and Table III). The discriminatory ability of the classification factor F to distinguish CAD and healthy controls in the training set was assessed using ROC analysis and showed an AUC of 0.96 (Fig. 2A). A classification threshold, FCAD = 13.0, classified the training set with sensitivity of 97% (95% CI, 83–99) and specificity of 80% (95% CI, 56–94). Furthermore a test set including patients after coronary artery bypass graft surgery and patients after acute coronary syndrome was subject to a blinded study. Of 59 urine samples 48 scored positive as "CAD," and 11 scored negative as healthy controls using the threshold of FCAD = 13.0. After unblinding, 46 of 47 CAD samples and 10 of 12 healthy control samples were predicted correctly (Fig. 2C) with sensitivity of 98% (95% CI, 89–99) and specificity of 83% (95% CI, 52–97). A ROC analysis showed an AUC of 0.94 (Fig. 2B).
|
|
|
To test whether this assessment can also be used to evaluate the effects of therapeutic interventions, we analyzed patients who had undertaken different levels of physical exercise. Self-reported and physiotherapist-validated activity levels determined improvement of the classification factor between base-line and follow-up assessment. At follow-up assessment the classification factor was almost unchanged in inactive and low grade active patients (
F = –0.20), whereas the very active patient group (
F = –1.90) showed a significant improvement (p = 0.02).
To determine the identity of biomarkers in the panel (Table III) used to distinguish CAD from healthy controls we performed LC-MS/MS to obtain sequence information. Examples of sequences deduced from high resolution fragmentation spectra are shown in Fig. 3, and five sequences of biomarkers in the panel could be identified: collagen
-1(I) chains and collagen
-1(III) fragments (Table IV). In all cases, the identified collagen type I or III fragments were up-regulated in CAD samples compared with controls. From the initial marker list 38 polypeptides (supplemental table) could be identified. The majority of these sequences were collagen fragments.
|
|
| DISCUSSION |
|---|
|
|
|---|
Our polypeptide pattern distinguished between the presence and absence of disease. Furthermore we were also able to demonstrate a dynamic behavior of the polypeptide pattern in response to exercise. Inactive patients had no change in pattern over time, whereas very active patients showed significant changes toward a "healthier" biomarker pattern. These results illustrate an important difference between proteomics versus genomics analysis. Genomics analysis identifies predisposing risk factors, whereas proteomics can identify the point in time when predisposition is developing into disease. This is because the proteome is inherently dynamic and hence can better reflect changes (27). This also pertains to measuring the effects of therapeutic interventions and will assist in personalized medicine strategies. Furthermore this observation also provides a link between the kidney and CAD on a molecular basis and may help to explain why chronic kidney disease is one of the best predictors for CAD (28).
For a reliable CAD-specific polypeptide panel it is mandatory that the constituting markers are independent of medication effects because the majority of patients require multiple drugs including antihypertensives. As an example, we used an additional control group of patients before and after treatment with an angiotensin-converting enzyme inhibitor, ramipril, to refine selected polypeptides. This comparison offered a more precise and sensitive monitoring of medication effect that is independent of individual fluctuations in urinary peptide patterns. The obtained polypeptide pattern was validated using a test set including patients with CAD in addition to healthy controls.
We were also able to identify five of the polypeptides constituting the CAD-specific panel: all of them were collagen type I or III fragments. These collagens are predominant proteins in the arterial walls. They also appear together in the thickened intima of atherosclerotic lesions (29). Collagen is synthesized containing C-terminal and N-terminal propeptide sequences (PICP and PINP for type I procollagen and PIIICP and PIIINP for type III procollagen). These propeptides are used as biomarkers of the rate of collagen synthesis (PICP), collagen degradation (ICTP), and collagen turnover (PIIINP) (30–34). In the normal artery, both synthesis and degradation of extracellular matrix proteins are remarkably slow. Atherosclerosis leads to increased synthesis of many matrix components, including collagen types I and III, elastin, and several proteoglycans (35).
All sequenced collagen fragments were up-regulated in CAD samples compared with controls. In line with these sequence data suggesting elevated collagen degradation levels, increased circulating levels of collagenases, such as MMP-9, have been reported in patients with stable angiographic coronary atherosclerosis (36, 37) or intermittent claudication (38). In patients with stable CAD, circulating MMP-9 levels independently predict rapid lumen diameter reduction (39) and fatal cardiovascular events (40). Elevated MMP-9 activity has been found in unstable plaques, suggesting a crucial role in plaque rupture (41, 42).
The majority of the identified polypeptides constituting the initial marker list were also collagen fragments (supplemental table). These findings suggest that CAD-specific information is redundantly available in urine samples in the form of different detectable collagen fragments. In addition to the collagen fragments, a fragment of membrane-associated progesterone receptor component 1 was identified (supplemental table). Progesterone receptors are reported to be associated with thoracic ascending aorta, internal carotid artery, coronary artery, and left atrial appendage (43).
Several groups have reported on the application of proteomics techniques to analyze tissue or plaque specimens to study cardiovascular disease or arteriosclerosis (44–48). Although these results provide new insights into disease-related pathways, they do not allow non-invasive detection of coronary artery disease. Only a few of these studies focused on the analysis of body fluids derived through minimally invasive means. Furthermore either these studies were based on pooled blood specimens (49, 50), making individual sample classification impossible, or they included only small patient cohorts without blinded studies for validation (51). All of these limitations are avoided in our study. The proteomics analysis of urine allowed for the reproducible and standardized analysis of a non-invasively obtained body fluid for highly accurate detection of CAD, which was subsequently validated in a blinded study.
In summary, in this patient population a CAD-specific urinary proteome panel is an accurate and non-invasive predictor for CAD. These and other similar biomarkers have the potential to be used for early diagnosis and thus more efficient prophylaxis as well as monitoring of therapeutic interventions and as novel drug targets.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Published, MCP Papers in Press, October 19, 2007, DOI 10.1074/mcp.M700394-MCP200
1 The abbreviations used are: CAD, coronary artery disease; CE, capillary electrophoresis; ETD, electron transfer dissociation; ROC, receiver operating characteristic; AUC, area under the ROC curve; CI, confidence interval; DTA, .dta file; LDL, low density lipoprotein; HDL, high density lipoprotein. ![]()
* This work was supported in part by the British Heart Foundation Chair and Programme Grant BHF PG/02/128, Wellcome Trust Cardiovascular Functional Genomics Initiative 066780/2/012, and European Union InGenious HyperCare Grant LSHM-CT-2006-037093 (to A. F. D. and H. M.); by the Joint Infrastructure Fund funding for the Sir Henry Wellcome Functional Genomics Facility; and by the Biotechnology and Biological Sciences Research Council and Engineering and Physical Sciences Research Council Radical Solutions for Researching the Proteome (RASOR) grant. H. Mischak is founder and co-owner of Mosaiques Diagnostics, which developed the CE-MS technology and the MosaiquesVisu software. E. Schiffer and P. Zürbig are employees of Mosaiques Diagnostics. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. ![]()
S The on-line version of this article (available at http://www.mcponline.org) contains supplemental material. ![]()
b Both authors contributed equally to this work. ![]()
c Supported by Swiss National Science Foundation Grant PBBSB-105860 and the Lichtenstein-Stiftung, Basel, Switzerland. Present address: Medical Outpatient Dept., University Hospital, CH-8091 Zurich, Switzerland. ![]()
f Supported by a National Institutes of Health predoctoral fellowship (Biotechnology Training Program Grant NIH 5T32GM08349). ![]()
l Supported by a personal fellowship from the Deutsche Forschungsgemeinschaft (Grant DE 826/1-1). ![]()
m To whom correspondence should be addressed: BHF Glasgow Cardiovascular Research Centre, 126 University Place, University of Glasgow, Glasgow G12 8TA, UK. Tel.: 44-141-330-5420; Fax: 44-141-330-6997; E-mail: ad7e{at}clinmed.gla.ac.uk
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. Bramham, H.D. Mistry, L. Poston, L.C. Chappell, and A.J. Thompson The non-invasive biopsy--will urinary proteomics make the renal tissue biopsy redundant? QJM, June 24, 2009; (2009) hcp071v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. I. Balog, P. J. Hensbergen, R. Derks, J. J. Verweij, G. J. van Dam, B. J. Vennervald, A. M. Deelder, and O. A. Mayboroda Novel Automated Biomarker Discovery Work Flow for Urinary Peptidomics Clin. Chem., January 1, 2009; 55(1): 117 - 125. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Rossing, H. Mischak, M. Dakna, P. Zurbig, J. Novak, B. A. Julian, D. M. Good, J. J. Coon, L. Tarnow, P. Rossing, et al. Urinary Proteomics in Diabetes and CKD J. Am. Soc. Nephrol., July 1, 2008; 19(7): 1283 - 1290. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Journal of Biological Chemistry |
| Journal of Lipid Research | ASBMB Today |