Skip to main content

Serum protein profiles predict coronary artery disease in symptomatic patients referred for coronary angiography

Abstract

Background

More than a million diagnostic cardiac catheterizations are performed annually in the US for evaluation of coronary artery anatomy and the presence of atherosclerosis. Nearly half of these patients have no significant coronary lesions or do not require mechanical or surgical revascularization. Consequently, the ability to rule out clinically significant coronary artery disease (CAD) using low cost, low risk tests of serum biomarkers in even a small percentage of patients with normal coronary arteries could be highly beneficial.

Methods

Serum from 359 symptomatic subjects referred for catheterization was interrogated for proteins involved in atherogenesis, atherosclerosis, and plaque vulnerability. Coronary angiography classified 150 patients without flow-limiting CAD who did not require percutaneous intervention (PCI) while 209 required coronary revascularization (stents, angioplasty, or coronary artery bypass graft surgery). Continuous variables were compared across the two patient groups for each analyte including calculation of false discovery rate (FDR ≤ 1%) and Q value (P value for statistical significance adjusted to ≤ 0.01).

Results

Significant differences were detected in circulating proteins from patients requiring revascularization including increased apolipoprotein B100 (APO-B100), C-reactive protein (CRP), fibrinogen, vascular cell adhesion molecule 1 (VCAM-1), myeloperoxidase (MPO), resistin, osteopontin, interleukin (IL)-1β, IL-6, IL-10 and N-terminal fragment protein precursor brain natriuretic peptide (NT-pBNP) and decreased apolipoprotein A1 (APO-A1). Biomarker classification signatures comprising up to 5 analytes were identified using a tunable scoring function trained against 239 samples and validated with 120 additional samples. A total of 14 overlapping signatures classified patients without significant coronary disease (38% to 59% specificity) while maintaining 95% sensitivity for patients requiring revascularization. Osteopontin (14 times) and resistin (10 times) were most frequently represented among these diagnostic signatures. The most efficacious protein signature in validation studies comprised osteopontin (OPN), resistin, matrix metalloproteinase 7 (MMP7) and interferon γ (IFNγ) as a four-marker panel while the addition of either CRP or adiponectin (ACRP-30) yielded comparable results in five protein signatures.

Conclusions

Proteins in the serum of CAD patients predominantly reflected (1) a positive acute phase, inflammatory response and (2) alterations in lipid metabolism, transport, peroxidation and accumulation. There were surprisingly few indicators of growth factor activation or extracellular matrix remodeling in the serum of CAD patients except for elevated OPN. These data suggest that many symptomatic patients without significant CAD could be identified by a targeted multiplex serum protein test without cardiac catheterization thereby eliminating exposure to ionizing radiation and decreasing the economic burden of angiographic testing for these patients.

Peer Review reports

Background

Coronary heart disease is the most prevalent chronic disease and the leading cause of death in the US, with more than half a million newly diagnosed coronary artery disease (CAD) patients annually [1, 2]. Cardiac catheterization and coronary angiography are often necessary for definitive evaluation of coronary artery anatomy, the presence of coronary atherosclerosis and to determine the need for interventional therapy. Despite the high prevalence of CAD, approximately half of patients undergoing invasive cardiac catheterization either have no significant coronary lesions or do not require any mechanical or surgical form of revascularization [3–5]. Thus, the procedure could be eliminated in many cases if alternative, non-invasive tools were available to assess the presence or absence of significant CAD and confirm the need for angiography.

Clinical assessment of CAD represents a significant medical and economic challenge comprising more than a million coronary angiograms annually in the US alone with demographics of aging and obesity forecasting growing demand [2–5]. The risk and expense of cardiac catheterization (ionizing radiation, contrast media, morbidity) and the large number of patients with normal coronary arteries or 'non-significant' CAD undergoing invasive angiography warrant development of alternative tests for CAD without cardiac catheterization [5]. While progress has been made using non-invasive computed tomography (CT) particularly for its negative predictive value, CT incorporates significant exposure to ionizing radiation with considerably lower resolution than catheterization-based angiography [6].

Efforts to identify circulating biomarkers for CAD have shown promise by interrogating transcriptional profiles of patient blood cells and plasma for unique mRNA and microRNA signatures [7, 8]. Since extracellular RNA undergoes rapid degradation, prospective mRNA signatures were derived predominantly from nucleated blood cells while the miRNAs identified in plasma were likely protected in circulating vesicles or bound to protective protein complexes [9]. Consequently, the utility of RNA as an indicator of CAD is constrained by its selective cell source in the bloodstream, the friability of the ribonucleotide targets and the arduous process of RNA capture, purification, amplification and analysis. In contrast, circulating proteins are more stable than RNA in blood and serum with several individual markers identified previously as prospective biomarkers for the presence of atherosclerosis, myocardial infarction, heart failure, or markers of pathways involved in these cardiac conditions such as inflammation, thrombosis, plaque stability, and oxidative stress, for example, troponin C, pro-brain natriuretic peptide (BNP) and C-reactive protein (CRP) [10, 11]. However, circulating biomarkers have proven to be of limited value in clinical tests to diagnose coronary artery disease antecedent to a cardiac event, primarily because most studies have focused on single or at most a few markers to make the diagnosis [12]. The difficulty in identifying predictive factors of CAD in blood or serum is compounded by the multifactorial etiology of coronary artery disease which makes early diagnosis by a single, endpoint marker unlikely prior to activation of a common ischemic pathway or until significant myocardial compromise has occurred.

The hypothesis underlying the current study was that coronary artery disease status can be assessed via individual and/or combinatorial protein changes in serum that assess multiple pathways of atherosclerosis as a low-risk, non-invasive approach for screening of symptomatic patients, that is, patients referred for cardiac catheterization. The study targeted patients who were referred for a clinically-indicated cardiac catheterization from either the emergency room or outpatient cardiac clinic in a major metropolitan center who presented with symptoms of heart disease. All patients had blood drawn prior to coronary angiography and revascularization. By analyzing a compendium of 41 circulating protein targets associated with atherogenesis, inflammation, thrombosis, and plaque vulnerability, we discovered 12 diverse proteins expressed across a broad dynamic range that were significantly different concomitant with the need for these patients to undergo therapeutic intervention including stent placement, angioplasty, or coronary artery bypass graft surgery (CABG). We also tested multiplex biomarker signatures for the potential to discriminate patients lacking significant coronary artery disease from patients with CAD requiring corrective interventional therapy. In particular, the ability to rapidly and decisively rule out clinically significant coronary artery disease using a potentially low cost, low risk blood test in even a small percentage of patients with normal coronary arteries could be highly beneficial.

Methods

Patient group

Samples comprised serum from among 359 subjects referred for a clinically-indicated cardiac catheterization for symptoms of CAD. The study was performed according to the Department of Health and Human Services Code of Federal Regulations (45 CFR 46) for the protection of human research subjects including ethical considerations consistent with the Office for Human Research Protections. Blood was collected following a genetic banking protocol (#990835) approved by the University of Pittsburgh Institutional Review Board (IRB). Only patients who signed the approved IRB voluntary informed consent document for this study were included (11 January 2000 to 21 July 2004). Venous blood (5 ml) was drawn into a red top vacutainer serum tube (Becton Dickinson #366430, Franklin Lakes, NJ, USA) and placed upright 30 to 60 minutes until clot formation. The tubes were centrifuged in a swinging bucket rotor (1,300 g × 20 min) and the serum was pipetted into 1.5 ml cryovials for storage at -80°C. All 359 patients underwent diagnostic coronary angiography and 209 required interventional therapy comprising stent placement, angioplasty or CABG. The other patients had normal or clinically insignificant coronary artery disease, that is, angiography revealed absence of any vessel obstruction or non-critical, < 50% obstruction in any epicardial vessel. Thus, serum samples comprised one patient group based on clinical symptoms but two outcome groups based on therapeutic intervention.

Proteomics analysis

Samples underwent a first thaw on ice to apportion into 200 μl aliquots for -80°C storage until analysis. An exploratory study of 56 samples was performed using fluorokine multianalyte profiling (xMAP) of 33 proteins (Luminex 100; Luminex, Austin, TX, USA) to determine serum dilution factors and rule out targets lacking statistical discrimination. The assay used polystyrene microspheres incorporating differing ratios of two fluorophores yielding different spectrally addressed bead sets each conjugated with a biotinylated protein-specific capture antibody. Assays were processed in duplicate in a randomized, blinded manner regarding patient outcomes, including generation of a standard curve using recombinant target proteins. Each 96 well microplate was incubated overnight at 4°C on a microtiter shaker. Wells were washed with buffer (3 ×) and a secondary antibody added to each well for incubation (2 h, room temperature) followed by streptavidin-phycoerythrin (0.5 h, room temperature, agitation). The wells were then washed (2 ×), assay buffer was added, and samples were analyzed using the Bio-Plex suspension array system and Bio-Plex Manager software 4.0 (Bio-Rad Laboratories, Hercules, CA, USA). Absolute quantities were determined by comparison to the five-point standard curve for each analyte.

The Searchlight Protein Array System (Aushon Biosystems, Inc, Billerica, MA, USA) was used to interrogate patient serum samples at 2 different stages (stage 1: 239 samples, 24 analytes; stage 2: 120 samples, 10 analytes). First, 239 samples were evaluated for 24 analytes over concentration ranges defined by the preliminary study of 56 samples. The assay comprised a multiplex sandwich ELISA of monoclonal capture antibodies spotted in planar arrays in 96-well microtiter plates. After serum incubation and washing, a second biotinylated monoclonal antibody to a different site from the capture epitope was introduced and streptavidin-horseradish peroxidase (HRP) subsequently bound to the biotin site. Luminol enhancer/peroxidase solution was added and the HRP catalyzed oxidation of luminol to 3-aminophthalate resulting in light emission at 428 nm. A chemiluminescent image was acquired and processed using a four-parameter curve fit program (SearchLight Array Analyst Software) to compare the experimental samples to the recombinant calibration curve run in parallel wells to derive absolute concentrations adjusted for dilution and quality values.

The largest SearchLight panel simultaneously evaluated seven analytes diluted 1:1 (volume/volume) (dilution factor (df) = 2 ×) in assay buffer (RPMI1640 without phenol red + 10% heat inactivated FBS) including interferon γ (IFNγ), interleukin 1β (IL-1β), IL-6, IL-10, matrix metalloproteinase protein 1 (MMP1), thrombomodulin (TM) and tumor necrosis factor α (TNFα). Leptin, platelet endothelial cell adhesion molecule 1 (PECAM-1), endothelial leukocyte adhesion molecule 1 (E-selectin), monocyte chemoattractant protein 1 (MCP-1), MMP7 and vascular cell adhesion molecule 1 (VCAM-1) were assayed together at 25 × dilution factor. Tissue inhibitor of metalloproteinase 1 (TIMP-1), fibrinogen, resistin, leukocyte selectin (L-selectin) and myeloperoxidase (MPO) (df = 1,000 ×) were analyzed in a five-analyte panel. Adiponectin (ACRP-30) and C-reactive protein (CRP) were assayed together at a dilution factor of 5,000 ×. Apolipoprotein A1 (APO-A1, df = 50,000 ×), apolipoprotein B100 (APO-B100, df = 10,000 ×), osteopontin (OPN, df = 10 ×) and N-terminal fragment protein precursor brain natriuretic peptide (NT-pBNP, df = 2 ×) were interrogated independently.

A second stage study of 120 additional serum samples was repeated twice to validate the previous findings from the 239 sample set and to test assay reproducibility across different reagent and planar array lots. The sample prep, quality control (QC), methodological protocols for recombinant protein calibration profiles, serial dilutions and serum assays were performed as before but using fewer panels and smaller analyte configurations. These included MPO, fibrinogen and resistin (df = 1,000 ×) in a three-multiplex configuration, ACRP-30 and APO-B100 together (df = 10,000 ×), MMP7 and VCAM-1 together (df = 25 ×) and osteopontin (df = 10 ×), IFN-γ (df = 2 ×) and IL-1β (df = 2 ×) separately.

Statistical analysis

Patients were operationally defined as 'symptomatic' by referral for a clinically indicated catheterization. Based on coronary angiography outcome, the serum samples were classified from patients with 'normal' coronary arteries, that is, no clinically significant coronary artery disease (n = 150) or patients with coronary disease requiring therapy including stent placement, angioplasty or CABG (n = 209). The hypothesis undergoing statistical testing was that serum proteins were significantly different between the two patient outcome classifications. Statistical analysis was initially performed on 239 samples in stage 1 evaluated for 24 analytes comprising 101 serum samples from patients with clinically normal coronary arteries and 138 samples from patients requiring percutaneous intervention (PCI). These samples were also used to develop and train a scoring function algorithm. A second stage validation study (n = 120) interrogating 10 analytes was subsequently performed to validate the algorithm. Results from all of these studies were combined for statistical comparison. Continuous variables were compared (Partek Genomics Suite, St. Louis, MO, USA) using the unpaired Student's t test across the two patient groups for each analyte including calculation of a false discovery rate (FDR) and Q value as the minimum positive FDR for rejecting a statistic [13]. Significant differential expression of proteins was defined at an FDR of ≤ 1% with Q = 0.01 and statistical significance for P values was adjusted to ≤ 0.01. Categorical variables were compared using the Pearson's χ2 test.

Algorithm development and validation of selected markers

Data for all 24 markers interrogated in stage 1 were evaluated as randomized, multimarker signatures to classify patients with CAD requiring treatment versus patients without clinically significant CAD. A scoring function (SF) algorithm was generated for all protein combinations as disease 'signatures' including 24 'artificial' markers derived by randomly scrambling the data (see Additional File 1). The SF for each signature was a linear combination of natural logarithms of marker concentrations generated by iterative computation. Monte Carlo optimization determined coefficients that provided highest diagnostic accuracy, that is, specificity (SP: identification of negatives for significant CAD) for detecting patients with normal coronary arteries while maintaining 95% sensitivity (SN: identification of CAD requiring interventional therapy) for patients with coronary artery disease. We ranked > 2 million combinations of 2 to 5 marker signatures comprising the 24 actual and 24 'artificial' markers for ability to classify patients since combinations of 6 or more proteins with high classification strength commonly contained an artificial marker (see Additional file 2 Table S1). For each signature of 2 to 5 markers, the top 50 panels with highest SP for normal (while correctly detecting at least 95% of the CAD patients) underwent cross validation testing where 80% of participants were randomly selected as a training set to build the optimal SF and the remaining 20% of participants were thenclassified using this SF. The crossvalidation procedure was repeated 500 times and average SP and SN were used to identify best performing signatures.

Independent verification of the scoring function algorithm was performed in 2 repeated studies of 120 serum specimens from an additional cohort of symptomatic patients with clinical characteristics matching the previous 239 patients. Concentration values for these samples were entered into the algorithm in a macro subroutine program using the offset, coefficients and cut-offs to detect CAD based on patient outcome. The results of the 120-sample validation study were compared to the diagnostic classification of each patient after coronary catheterization and follow-up therapy to determine the sensitivity and specificity of prospective signatures.

Results

Diagnostic coronary angiography revealed that 209 of the patients in this study exhibited significant coronary artery disease requiring therapeutic intervention while 150 patients did not exhibit clinically significant coronary artery disease despite symptoms or other findings that led to referral for cardiac catheterization. These two distinct outcome groups were otherwise identical upon admission regarding clinical symptoms and physical characteristics including gender, diabetic status, smoking history, body surface area, basal metabolic rates, cholesterol, LDL and creatinine values (see Table 1). Among continuous variables, there were small albeit significant differences in age, HDL levels and ejection fraction between groups; but the differences were of minimal diagnostic value and all patients required coronary angiography. Regarding categorical variables, there were no significant differences in gender or diabetes between the two groups; however, the number of patients with hypertension was significantly higher in the CAD group.

Table 1 Clinical characteristics of the patient groups

All serum samples were collected, processed, stored and analyzed in an identical manner to limit the effect of preanalytic variability including differential protein degradation among specimens. Significant differences were detected in 12 serum proteins (Q value = 0.01; P < 0.01) between patients diagnosed as having CAD requiring intervention and those with non-significant CAD after diagnostic coronary angiography. The differences detected in the stage 1 study (n = 239) were reinforced by the additional samples from the validation study (n = 120) (see Table 2). APO-A1 and APO-B100 were among the highest expressed proteins overall averaging approximately 300 μg/ml of serum (Figure 1). APO-A1 fell significantly in patients with significant CAD versus non-significant CAD while APO-B100 was significantly increased. Within the same concentration range, fibrinogen was present at levels typically exceeding 1 μg/ml with values fivefold higher in patients with significant CAD (Figure 1). At serum concentrations from 10 ng to 1 μg/ml serum, five proteins were significantly higher in CAD patients. Specifically, CRP, VCAM-1, MPO, resistin and osteopontin were 1.2 to 3.1-fold higher than patients with non-significant CAD (Figure 2). Four analytes, IL-6, IL-1β, IL-10 and NT-pBNP were significantly higher in the CAD group among analytes detected in a range from 1 pg/ml to 1 ng/ml (Figure 3). There were no significant statistical correlations between any of these 12 analytes and age, ejection fraction or hypertension status which were significantly but incrementally different between the patient outcome groups. No other analytes among the 41 interrogated targets were significantly altered between the 2 groups of patient samples using either bead-based or planar platforms.

Table 2 Multiplex proteomics analysis of the coronary artery disease (CAD) and normal (NOR) patient groups
Figure 1
figure 1

Significant differences in apolipoprotein A1 (APO-A1), apolioprotein B100 (APO-B100) and fibrinogen in serum from normal and coronary artery disease (CAD) patients. Solid bars are values expressed as average plus 1 SD for APO-A1, APO-B100, and fibrinogen obtained from patients without clinically significant coronary artery disease (normal, n = 150) based on coronary angiographic evaluation. Open bars are results obtained from patients with CAD requiring interventional therapy (n = 209). Values are expressed in μg/ml on a logarithmic ordinate scale and each was significantly different (*) between groups (see Table 2 for individual P values).

Figure 2
figure 2

Significant differences in vascular cell adhesion molecule, myeloperoxidase, C-reactive protein, resistin and osteopontin in serum from normal and coronary artery disease (CAD) patients. Normal and CAD data are displayed according to Figure 1 but expressed in ng/ml on a logarithmic ordinate scale. All comparisons represent significant statistical differences delineated in Table 2 (*) for vascular cell adhesion molecule (VCAM-1), myeloperoxidase (MPO), C-reactive protein (CRP), resistin and osteopontin (OPN).

Figure 3
figure 3

Significant differences in interleukin (IL)-6, IL-1β, IL-10 and N-terminal fragment pro-brain natriuretic peptide (NT-pBNP) in serum from normal and coronary artery disease (CAD) patients. Normal and CAD data are displayed according to Figure 1 but expressed in pg/ml on a logarithmic ordinate scale. All comparisons represent significant statistical differences (*) reported in Table 2 for IL-6, IL-1β, IL-10 and NT-pBNP.

We identified 14 multiplex signatures of 2 to 5 proteins each with the highest acuity to detect patients without significant CAD (22.6% to 58.4% SP) while detecting 95% of the significant CAD group (95% SN) in the stage 1 study (see Additional file 2 Table S1). A total of 11 distinct proteins were shared among the 14 signatures with osteopontin (14 of 14), and resistin (10 of 14) most frequently represented. There was a trend for protein signatures with increased numbers of analytes to detect more normal patients at a fixed sensitivity for CAD patients (95%) (two proteins = 39.3% ± 0.3% vs five proteins = 50.0% ± 0.01% of normal patients). However, a performance plateau was reached at five biomarkers based on crossvalidated classifier performance and the frequency of appearance of 'artificial' markers in test signatures exceeding five proteins. Receiver operating characteristics analysis indicated that these signatures were effective in discerning patients without significant CAD. The area under the curve (AUC) for the top signatures ranged from a low of 0.839 ± 0.028 (mean ± SD) for a two-protein signature (OPN, resistin) to a maximum AUC of 0.845 using three biomarkers (OPN, resistin, APO-B100) (Figure 4). These ROC curves were compared to those generated by the Bayesian compound covariate predictor algorithm for the same data set. The area under the curve using the scoring function algorithm exceeded that obtained by the Bayesian predictor in every case. A clinical validation test of 120 additional serum samples (49 normal, 71 patients requiring intervention) was performed to test the performance of the scoring function algorithm. In two separate studies, the best performing multiplex signatures contained five proteins (OPN, resistin, MMP7, IFNγ with either CRP or ACRP-30) and were able to correctly classify 88% and 92% of patients requiring percutaneous intervention while delineating 33% and 36% of the patients with normal coronary arteriograms.

Figure 4
figure 4

Receiver operating characteristics (ROC) for 2 to 5 protein panels for identification of normal patients with 95% specificity for detection of coronary artery disease (CAD) patients. The ROC curves are derived from 4 separate panels optimized to detect 101 normal patients (true positives in this figure) at highest specificity while maintaining a sensitivity of 95% for patients with CAD (138 samples). The ROC curves are obtained by iteratively testing each biomarker panel for classification of a randomly excluded portion (20%) of the dataset. The areas under the curve (AUC) were comparable as indicated in the curves for two proteins (osteopontin (OPN) and resistin: AUC = 0.839), three proteins (OPN, resistin, apolioprotein B100 (APO-B100): AUC = 0.845), four proteins (OPN, resistin, matrix metalloproteinase 7 (MMP7) and interferon γ (IFNγ): AUC = 0.839) and five proteins (OPN, resistin, MMP7, IFNγ and C-reactive protein (CRP): AUC = 0.827). The predicted specificity for detection of normal patients at 95% sensitivity for CAD patients was two proteins = 50%, three proteins = 52%, four proteins = 63% and five proteins = 64%.

Discussion

Proteins were selected for evaluation in this study based on their roles in mechanisms underlying atherogenesis, atherosclerosis and plaque instability including vascular inflammation, thrombosis, aberrant lipid regulation, metabolism hormones, and vascular smooth muscle and extracellular matrix (ECM) remodeling [14]. The 41 preliminary targets we interrogated were restricted by availability of monoclonal antibody pairs optimized for use in the commercial assay platforms. IL-1β, IL-6, IL-10 and VCAM-1, were significantly elevated in patients with CAD in the present study consistent with an injury-induced, inflammatory response [15, 16]. Elevated IL-1β and IL-6 have been associated previously with acute phase protein induction and may explain the concomitant significant increases in fibrinogen and CRP concentration we detected. CRP has long been proposed as a surrogate marker for inflammatory mediators in predicting coronary events while NT-pBNP has been used as an indicator of left ventricular dysfunction in CAD patient cohorts comparable to this study [11, 17, 18]. Both analytes were significantly elevated in the present study among patients requiring therapeutic intervention and CRP was among the best single molecule classifiers delineating 19% of normal samples while detecting 95% of the patients with significant CAD.

Significant reciprocal changes were detected in APO-A1 and APO-B100 in CAD patients consistent with reports defining aberrant lipid transport and accumulation as contributory to atherosclerosis [19]. Mutations in the APO-B100 gene cause autosomal dominant, hereditary familial hypercholesterolemia and premature coronary artery disease due to defective ligand binding [19, 20]. At the same time, elevated APO-A1 is associated with a cardioprotective effect and enhancement of APO-A1 expression has been proposed as a therapeutic strategy to inhibit atheroma formation [19, 21]. The increased APO-B100 and decreased APO-A1 levels in our patients requiring PCI versus normal controls support these previous findings. Myeloperoxidase was also significantly increased in CAD patients associated with its role as a catalyst for lipid peroxidation at inflammation sites and as a marker of plaque instability [22, 23]. Resistin levels were elevated in the PCI patients indicative of 1) metabolic shifts in lipid utilization and adipogenesis and/or 2) an inflammatory response with resistin secreted from macrophages concomitant with the release of proinflammatory cytokines [24].

Many targets traditionally associated with vascular smooth muscle and ECM remodeling were not significantly altered among these patient groups including matrix metalloproteinases 1, 2, 3, 7, 9 and tissue inhibitors of metalloproteinases 1, 2, 3 and 4. Only osteopontin, which acts as a negative regulator of calcification in bone remodeling, was elevated within this category with the rejoinder that OPN also may act as a chemokine in the cell mediated type 1 immune response associated with inflammatory cell accumulation rather than as a substrate for cell adhesion [25]. Thus, the proteins that demarcated our patient outcome groups were predominantly associated with processes of inflammation and lipid regulation rather than cellular aggregation and ECM remodeling. However, we recognize that the domain of proteins susceptible to interrogation in this study was limited to analytes for which high affinity antibody pairs precisely characterized to two different epitopes were available. The involvement of additional proteins and pathways associated with CAD will likely be reinforced and/or revealed as the inventory of immunoassays becomes more comprehensive.

Our data indicate multiplex proteomics analyses using monoclonal antibodies provide relevant information regarding circulating serum analyte concentrations when assayed at a dilution that allows direct comparison to parallel recombinant calibration standards. Advantages include small serum volumes (< 100 μl) collected by standard clinical protocols, rapid turnaround times (minutes to hours), high sensitivity (pg) and a broad dynamic range (8 logs). Disadvantages include high assay cost, limited target availability and poor concurrence of concentration measurements across dilutions and commercial platforms associated with variations in antibodies, buffers, diluents and capture structures. In the present study, 15 targets were tested at identical serum dilutions using bead-based (Luminex) and planar (Aushon) technologies in 56 identical samples, albeit with different aliquots and in serial studies. A total of 12 assays concurred in detection of statistically significant differences between the 2 patient outcome groups. These results suggest that multiplex immunochemical assays of serum may provide information of diagnostic relevance but that protocols and reagents must be optimized and standardized prior to routine clinical application.

The results of this study were somewhat surprising both for discovery of unique proteins as discriminants of CAD and for the absence of statistically significant differences in many targets with established roles in atherosclerosis. For example, osteopontin has been only indirectly associated with atherosclerosis yet exhibited the greatest statistical difference between patient groups (P = 1.75 × 10-12). Osteopontin was first identified as a sialoprotein from mineralized bone matrix and only recently was associated with calcification of plaques in cardiac valves and vessels [25–27]. Similarly, resistin has been linked only indirectly to CAD through a role in metabolic homeostasis and insulin sensitivity [28]. On the other hand, multiple growth factors (VEGF, leptin, ghrelin), lipoproteins (APO-A2, E, serum amyloid A: SAA), cell adhesion molecules (thrombospondin, PECAM-1, ICAM-1, selectins E, L, P) and MMP and TIMP targets associated with ECM remodeling exhibited no statistically significant differences. There are several potential explanations for the latter findings: (1) a rigorous statistical standard was utilized to avoid multiple testing errors and while MMP1, MMP7, ACRP-30, and leptin were borderline for statistical significance (P = 0.015, 0.045, 0.027, 0.027 respectively) they failed to reach the Q = 0.01 level established for significance with adjusted P values ≤ 0.01 in this study; (2) serum may not be an effective transducer of deleterious protein changes participating in structural rearrangements within the coronary vascular anatomy and the extracellular matrix; and (3) the patients comprised a diverse range of coronary obstruction and plaque vulnerability since they were selected for symptoms upon emergent presentation requiring diagnostic coronary angiography without the occurrence of a clinically obvious myocardial infarct or an 'event'. A subset of patients selected for advanced disease might reveal additional protein changes but stray from the intended focus of this study.

A scoring function algorithm was developed, tested and validated for the ability to classify patients symptomatic for heart disease consistent with the outcome of coronary angiography studies and need for interventional therapy. We minimized selection bias by testing a hypothesis driven biomarker panel and avoided overfitting by performing cross validation and follow-up testing utilizing additional serum samples from the cohort. The algorithm was designed to be 'tuned' to increase the sensitivity for capture of patients who required coronary revascularization at the expense of detecting fewer patients who did not require coronary revascularization. All serum signatures with highest classification strength from the training trial (239 samples) included osteopontin and signatures containing 4 or 5 proteins performed best during both training and validation phases. The most efficacious protein signature in validation studies comprised OPN, resistin, MMP7 and IFNγ as a four-marker panel while the addition of either CRP or ACRP-30 yielded comparable results in five protein signatures.

Further validation of the diagnostic accuracy of this approach will require extensive testing in greater numbers of patients at multiple locations as well as a prognostic cohort. It is possible that inclusion of clinical variables and risk factors in the biomarker algorithm or using the algorithm as part of a clinical scoring system will enhance both the fidelity and the efficacy of this approach for diagnostic purposes [29, 30]. In that context, we calculated 10-year Framingham Coronary Heart Disease (CHD) Risk Scores for patients where all clinical variables (gender, age, total cholesterol, HDL, systolic blood pressure, smoking and diabetes status, use of antihypertensive medication) were acquired prior to coronary angiography [31]. This represented 91 patients who subsequently required therapeutic revascularization (CAD: CHD Score = 14.9 ± 8.5) versus 63 patients who were determined to be free of significant coronary artery disease (no CAD: CHD score = 10.2 ± 6.7). The Framingham CHD scores were statistically different between groups (P < 0.001, unpaired Student's t test) but they classified only 16% of the subjects without significant CAD (10 of 63) at a 95% sensitivity for patients with CAD. In contrast, our algorithm incorporating serum values for OPN, RES, CRP, MMP7 and IFNγ identified 63% of the subjects without significant CAD (40 of 63) at 95% sensitivity for patients with CAD. Thus, our multiplex serum protein classifier correctly identified four times as many patients as the Framingham index. The strength of adding clinical variables to our scoring function remains to be determined, but the ability to exempt significant numbers of patients with normal coronary arteries or non-significant CAD from cardiac catheterization with a blood test represents a major economic and health benefit given the growing epidemic of CAD in the US and abroad.

Conclusions

The results of the present study indicate that a serum, multiplex biomarker assay may provide a clinically useful tool in combination with other standardized clinical tests to facilitate the decision-making process for performing cardiac catheterization in symptomatic patients. The development of highly sensitive monoclonal antibodies to additional pertinent targets along with the formulation of novel predictive algorithms will likely improve the efficacy of this approach. The long-term potential benefits include reduced patient exposure to ionizing radiation and minimization of the rapidly escalating healthcare costs associated with the use of invasive angiography to rule out coronary artery disease.

Abbreviations

ACRP-30:

adiponectin

APO:

apolipoprotein

AUC:

area under the curve

CABG:

coronary artery bypass graft surgery

CAD:

coronary artery disease

CHD:

coronary heart disease

CRP:

C-reactive protein

CT:

computed tomography

ECM:

extracellular matrix

E-selectin:

endothelial leukocyte adhesion molecule

ELISA:

enzyme linked immunosorbant assay

FDR:

false discovery rate

HDL:

high-density lipoprotein

HRP:

horseradish peroxidase

IFN:

interferon

IL:

interleukin

IRB:

institutional review board

LDL:

low-density lipoprotein

L-selectin:

leukocyte selectin

MCP:

monocyte chemoattractant protein

MMP:

matrix metalloproteinase protein

PCI:

percutaneous intervention

QC:

quality control

ROC:

receiver operator characteristics

SF:

scoring function

SN:

sensitivity

SP:

specificity

MPO:

myeloperoxidase

NT-pBNP:

N-terminal fragment protein precursor brain natriuretic peptide

OPN:

osteopontin

PECAM:

platelet endothelial cell adhesion molecule

RES:

resistin

SAA:

serum amyloid A

TIMP:

tissue inhibitor of metalloproteinase

TM:

thrombomodulin

TNF:

tumor necrosis factor

VCAM:

vascular cell adhesion molecule.

References

  1. Rodger VL, Go AS, Lloyd-Jones D, Adams RJ, Benjamin EJ, Berry JD, Borden WB, Bravata DM, Dai S, Ford ES, Fox CS, Fullerton HJ, Gillespie C, Hailpern SM, Heit JA, Howard VJ, Kissela BM, Kittner SJ, Lackland DT, Lichtman JH, Lisabeth LD, Makuc DM, Marcus GM, Marelli A, Matchar DB, Moy CS, Mozaffarian D, Mussolino ME, Nichol G, Paynter NP, et al: Heart disease and stroke statistics--2012 update: a report from the American Heart Association. Circulation. 2012, 125: e2-e220.

    Article  Google Scholar 

  2. Xu JQ, Kochanek KD, Murphy SL, Tejada-Vera B: Deaths: final data for 2007 National vital statistics reports. 2010, Hyattsville, MD: National Center for Health Statistics

    Google Scholar 

  3. Scanlon PJ, Faxon DP, Audet AM, Carabello B, Dehmer GJ, Eagle KA, Legako RD, Leon DF, Murray JA, Nissen SE, Pepine CJ, Watson RM, Ritchie JL, Gibbons RJ, Cheitlin MD, Gardner TJ, Garson A, Russell RO, Ryan TJ, Smith SC: ACC/AHA guidelines for coronary angiography. A report of the American College of Cardiology/American Heart Association Task Force on practice guidelines (Committee on Coronary Angiography). Developed in collaboration with the Society for Cardiac Angiography and Interventions. J Am Coll Cardiol. 1999, 33: 1756-1824. 10.1016/S0735-1097(99)00126-6.

    Article  CAS  PubMed  Google Scholar 

  4. Kolansky DM: Acute coronary syndromes: morbidity, mortality, and pharmacoeconomic burden. Am J Manag Care. 2009, 15 (Suppl): S36-41.

    PubMed  Google Scholar 

  5. Patel MR, Peterson ED, Dai D, Brennan JM, Redberg RF, Anderson HV, Brindis RG, Douglas PS: Low diagnostic yield of elective coronary angiography. N Engl J Med. 2010, 362: 886-895. 10.1056/NEJMoa0907272.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Hacioglu Y, Gupta M, Budoff MJ: Noninvasive anatomical coronary artyery imaging versus myocardial perfusion imaging: which confers superior diagnostic and prognostic information?. J Comput Assist Tomogr. 2010, 34: 637-644. 10.1097/RCT.0b013e3181e10525.

    Article  PubMed  Google Scholar 

  7. Rosenberg S, Elashoff MR, Beineke P, Daniels SE, Wingrove JA, Tingley WG, Sager PT, Sehnert AJ, Yau M, Kraus WE, Newby LK, Schwartz RS, Voros S, Ellis SG, Tahirkheli N, Waksman R, McPherson J, Lansky A, Winn ME, Schork NJ, Topol EJ, PREDICT (Personalized Risk Evaluation and Diagnosis in the Coronary Tree) Investigators: Multicenter validation of the diagnostic accuracy of a blood-based gene expression test for assessing obstructive coronary artery disease in nondiabetic patients. Ann Intern Med. 2010, 153: 425-434.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Fichtlscherer S, De Rosa S, Fox H, Schwietz T, Fischer A, Liebetrau C, Weber M, Hamm CW, Röxe T, Müller-Ardogan M, Bonauer A, Zeiher AM, Dimmeler S: Circulating microRNAs in patients with coronary artery disease. Circ Res. 2010, 107: 677-684. 10.1161/CIRCRESAHA.109.215566.

    Article  CAS  PubMed  Google Scholar 

  9. Wang K, Zhang S, Weber J, Baxter D, Galas DJ: Export of microRNAs and microRNA-protective protein by mammalian cells. Nucl Acids Res. 2010, 38: 7248-7259. 10.1093/nar/gkq601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Katus HA, Remppis A, Neumann FJ, Scheffold T, Diederich KW, Vinar G, Noe A, Matern G, Kuebler W: Diagnostic efficiency of troponin T measurements in acute myocardial infarction. Circulation. 1991, 83: 902-912. 10.1161/01.CIR.83.3.902.

    Article  CAS  PubMed  Google Scholar 

  11. James SK, Lindahl B, Siegbahn A, Stridsberg M, Venge P, Armstrong P, Barnathan ES, Califf R, Topol EJ, Simoons ML, Wallentin L: N-terminal pro-brain natriuretic peptide and other risk markers for the separate prediction of mortality and subsequent myocardial infarction in patients with unstable coronary artery disease: a Global Utilization of Strategies To Open occluded arteries (GUSTO)-IV substudy. Circulation. 2003, 108: 275-281. 10.1161/01.CIR.0000079170.10579.DC.

    Article  CAS  PubMed  Google Scholar 

  12. Moriates C, Maisel A: The utility of biomarkers in sorting out the complex patient. Am J Med. 2010, 123: 393-399. 10.1016/j.amjmed.2009.07.034.

    Article  CAS  PubMed  Google Scholar 

  13. Storey JD: The positive false discovery rate: A bayesian interpretation and the q-value. Ann Stat. 2003, 31: 2013-2035. 10.1214/aos/1074290335.

    Article  Google Scholar 

  14. Ross R: The pathogenesis of atherosclerosis: a perspective for the 1990s. Nat Rev. 1993, 362: 801-809.

    CAS  Google Scholar 

  15. Erren M, Reinecke H, Junker R, Fobker M, Schulte H, Schurek JO, Kropf J, Kerber S, Breithardt G, Assmann G, Cullen P: Systemic inflammatory parameters in patients with atherosclerosis of the coronary and peripheral arteries. Arterioscler Thromb Vasc Biol. 1999, 19: 2355-2363. 10.1161/01.ATV.19.10.2355.

    Article  CAS  PubMed  Google Scholar 

  16. Tedgui A, Mallat Z: Cytokines in atherosclerosis: pathogenic and regulatory pathways. Physiol Rev. 2006, 86: 515-581. 10.1152/physrev.00024.2005.

    Article  CAS  PubMed  Google Scholar 

  17. Jernberg T, Stridsberg M, Venge P, Lindahl B: N-terminal pro brain natriuretic peptide on admission for early risk stratification of patients with chest pain and no ST-segment elevation. J Am Coll Cardiol. 2002, 40: 437-445. 10.1016/S0735-1097(02)01986-1.

    Article  CAS  PubMed  Google Scholar 

  18. Di Angelantonio E, Chowdhury R, Sarwar N, Ry KK, Gobin R, Saleheen D, Thompson A, Gudnason V, Sattar N, Danesh J: B-Type natriuretic peptides and cardiovascular risk. Systematic review and meta-analysis of 40 prospective studies. Circulation. 2009, 120: 2177-2187. 10.1161/CIRCULATIONAHA.109.884866.

    Article  CAS  PubMed  Google Scholar 

  19. Walldius G, Jungner I, Aastveit AH, Holme I, Curt D, Furberg CD, Sniderman AD: The apoB/apoA-I ratio is better than the cholesterol ratios to estimate the balance between plasma proatherogenic and antiatherogenic lipoproteins and to predict coronary risk. Clin Chem Lab Med. 2004, 42: 1355-1363.

    Article  CAS  PubMed  Google Scholar 

  20. Pullinger CR, Hennessy LK, Chatterton JE, Liu W, Love JA, Mendel CM, Frost PH, Malloy MJ, Schumaker VN, Kane JP: Familial ligand-defective apolipoprotein B. Identification of a new mutation that decreases LDL receptor binding affinity. J Clin Invest. 1995, 95: 1225-1234. 10.1172/JCI117772.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Haas MJ, Mooradian AD: Therapeutic interventions to enhance apolipoprotein A-1 mediated cardioprotection. Drugs. 2010, 70: 805-821. 10.2165/11535410-000000000-00000.

    Article  CAS  PubMed  Google Scholar 

  22. Zhang R, Brennan ML, Shen Z, MacPherson JC, Schmitt D, Molenda CE, Hazen SL: Myeloperoxidase functions as a major enzymatic catalyst for initiation of lipid peroxidation at sites of inflammation. J Biol Chem. 2002, 277: 46116-46122. 10.1074/jbc.M209124200.

    Article  CAS  PubMed  Google Scholar 

  23. Naruko T, Furukawa A, Yunoki K, Komatsu R, Nakagawa M, Matsumura Y, Shirai N, Sugioka K, Takagi M, Hozumi T, Itoh A, Haze K, Yoshiyama M, Becker AE, Ueda M: Increased expression and plasma levels of myeloperoxidase are closely related to the presence of angiographically-detected complex lesion morphology in unstable angina. Heart. 2010, 96: 1716-1722. 10.1136/hrt.2009.187609.

    Article  CAS  PubMed  Google Scholar 

  24. Stofkova A: Resistin and visfatin: regulators of insulin sensitivity, inflammation and immunity. Endocr Regul. 2010, 44: 25-36. 10.4149/endo_2010_01_25.

    Article  CAS  PubMed  Google Scholar 

  25. Ashkar S, Weber GF, Panoutsakopoulou V, Sanchirico ME, Jansson M, Zawaideh S, Rittling SR, Denhardt DT, Glimcher MJ, Cantor H: Eta-1 (osteopontin): an early component of type-1 (cell-mediated) immunity. Science. 2000, 287: 860-864. 10.1126/science.287.5454.860.

    Article  CAS  PubMed  Google Scholar 

  26. Franzén A, Heinegård D: Isolation and characterization of two sialoproteins present only in bone calcified matrix. Biochem J. 1985, 232: 715-724.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Abdel-Azeez HA, Al-Zaky M: Plasma osteopontin as a predictor of coronary artery disease: association with echocardiographic characteristics of atherosclerosis. J Clin Lab Anal. 2010, 24: 201-206. 10.1002/jcla.20378.

    Article  CAS  PubMed  Google Scholar 

  28. Ye F, Than A, Zhao Y, Goh KH, Chen P: Vesicular storage, vesicle trafficking, and secretion of leptin and resistin: the similarities, differences, and interplays. J Endocrinol. 2010, 206: 27-36. 10.1677/JOE-10-0090.

    Article  CAS  PubMed  Google Scholar 

  29. Feng Z, Prentice R, Srivastava S: Research issues and strategies for genomic and proteomic biomarker discovery and validation: a statistical perspective. Pharmacogenom Rev. 2004, 5: 709-719. 10.1517/14622416.5.6.709.

    Article  CAS  Google Scholar 

  30. Bossuyt PM, Reitsma JB, Bruns DE, Gatsonis CA, Glasziou PP, Irwig LM, Moher D, Rennie D, deVet HCW, Lijmer JG: The STARD statement for reporting studies of diagnostic accuracy: explanation and elaboration. Ann Int Med. 2003, 138: 40-44.

    Article  PubMed  Google Scholar 

  31. Wilson PW, D'Agostino RB, Levy D, Belanger AM, Silbershatz H, Kannel WB: Prediction of coronary heart disease using risk factor categories. Circulation. 1998, 97: 1837-1847. 10.1161/01.CIR.97.18.1837.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

This study was supported by an NIH/NCI Cancer Center Support Grant (WALaF: co-investigator: 3P30 CA047904-20S1) and funding from Prevencio LLC (454 N 34th Street, Seattle, WA, USA) to the University of Pittsburgh and RD for the proteomics assays. The funding agencies had no role in the design and performance of the study, the interpretation of the data and the writing of the manuscript. The investigators thank Dr Robert D Guthrie for his thoughts, advice and critical input during the performance of this study, interpretation of the data and preparation of the manuscript. Many of the ideas underlying this study are derived from creative insights he provided more than a decade ago.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to William A LaFramboise.

Additional information

Competing interests

The University of Pittsburgh has applied for a patent related to the content of this manuscript including DMMcN, OCM, RVM, SRM, AL and WALaF as co-inventors. WALaF and RVM are founders and stockholders in Prevencio LLC, which provided partial support for performance of the proteomics assays. Neither WALaF nor RVM were directly involved in the selection of samples for analysis, performance of the primary proteomics assays or generation of the raw data reported in this study. RD, LAK, PP, JMK-B, CMS, MAL-W and URC have no competing interests associated with this study.

Authors' contributions

RD, WALaF, OCM, RVM, SRM and DMMcN conceived the study, developed the experimental design and analyte target list, analyzed the data and wrote the manuscript. OCM, SRM and DMMcN were responsible for acquisition and selection of all serum samples utilized in this study. LAK, JMK-B, MAL-W, CSM and PP performed all laboratory assays under the direction of RD, compiled the database and performed primary statistical analysis. JMK-B, AL, WALaF and URC performed secondary statistical analysis. AL and WALaF developed and tested the scoring function algorithm. All authors have read and approved the final version of this manuscript.

Electronic supplementary material

12916_2012_650_MOESM1_ESM.DOCX

Additional file 1: Algorithm development. The analytical form of the scoring function is provided. The testing that was performed to derive and optimize the algorithm is described including the use of artificial markers and cross validation testing. (DOCX 25 KB)

12916_2012_650_MOESM2_ESM.DOCX

Additional file 2: Predictive panel selection. The process of testing and final derivation of the panels comprising two, three, four or five markers that best predicted the classification outcome of the patients based on coronary angiography is provided. (DOCX 24 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

LaFramboise, W.A., Dhir, R., Kelly, L.A. et al. Serum protein profiles predict coronary artery disease in symptomatic patients referred for coronary angiography. BMC Med 10, 157 (2012). https://doi.org/10.1186/1741-7015-10-157

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1741-7015-10-157

Keywords