| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Circulation. 2007;116:888-893.)
© 2007 American Heart Association, Inc.
Coronary Heart Disease |
From the Department of Cardiovascular Medicine (M.S.L.), the Department of Thoracic and Cardiovascular Surgery (D.M., H.I., E.H.B.), and the Department of Quantitative Health Sciences (H.I., E.H.B.), The Cleveland Clinic Foundation, and the Department of Epidemiology and Biostatistics (M.S.L.), Case Western Reserve University School of Medicine, Cleveland, Ohio.
Correspondence to Dr Michael S. Lauer, Division of Prevention and Population Science, National Heart, Lung, and Blood Institute, 6701 Rockledge Dr, Room 10122, Bethesda, Md 20892. E-mail lauer{at}nhlbi.nih.gov
Received April 12, 2006; accepted June 8, 2007.
| Abstract |
|---|
|
|
|---|
Methods and Results— For 6 years we followed 8166 patients who underwent primary isolated coronary artery bypass grafting between 1990 and 2003, all of whom had routine preoperative ECGs. With use of specialized digital software, quantitative measures were recorded on ventricular rate, P duration, PR interval, QRS duration, QT interval, QRS axis, Sokolow-Lyon and Cornell voltages, and ST-segment depression and slope. There were 1516 deaths. After adjustment for age, gender, clinical characteristics, left ventricular ejection fraction, and other confounders, death was independently predicted by ventricular rate (adjusted hazard ratio [AHR] for 90 versus 60 beats per minute, 1.34; 95% confidence interval [CI], 1.21 to 1.50; P<.0001), PR interval (AHR for 200 versus 150 ms, 1.05; 95% CI, 1.00 to 1.10; P<.0001), QRS duration (AHR for 120 versus 80 ms, 1.24; 95% CI, 1.07 to 1.44; P<.0001), Sokolow-Lyon voltage (AHR for 3.5 versus 1.5 mV, 1.18; 95% CI, 1.05 to 1.31; P<.0001), and ST-segment slope (AHR for –0.1 versus 0 mV, 1.16; 95% CI, 1.02 to 1.31; P<.0001). We derived a quantitative ECG score and demonstrated that, with the exception of age, it was the most powerful predictor of long-term death.
Conclusions— Quantitative ECG measures of left ventricular rate, mass, and repolarization are predictive of mortality among patients who underwent isolated coronary artery bypass grafting. These findings suggest that quantitative electrocardiography may be valuable for risk stratification in patients with severe coronary artery disease.
Key Words: coronary disease electrocardiography prognosis surgery
| Introduction |
|---|
|
|
|---|
Clinical Perspective p 893
The value of quantitative electrocardiography in other routine clinical settings is less clear. As in most institutions, our hospital routinely obtains a preoperative ECG in patients who must undergo coronary artery bypass grafting (CABG). Given the prognostic associations of quantitative ECG measures of outcome in population-based and heart failure cohorts, we asked whether quantitative ECG measures obtained as part of routine clinical care may be predictive of long-term outcome in patients who undergo CABG, one of the most common cardiovascular procedures.
| Methods |
|---|
|
|
|---|
The cohort for the present study included 8166 patients who underwent primary isolated CABG between 1990 and 2003. All patients had an ECG obtained at the Cleveland Clinic as part of their routine preoperative evaluation within 2 weeks of surgery. Patients with atrial fibrillation, prior pacemaker implantation, Wolff-Parkinson-White syndrome, or pathological Q waves on their initial ECG were excluded.
Permission was obtained from the Institutional Review Board to analyze registry data for outcomes research. The requirement for written informed consent was waived.
Digital Electrocardiography
During the study period, all ECGs obtained at the Cleveland Clinic were recorded with a Marquette MUSE System. Digital ECG files were retrieved and analyzed with use of General Electrics Magellan Software System (GE Healthcare, Menomonee Falls, Wisc). The software provides detailed data on the duration and amplitudes of all segments of the P wave, QRS complex, ST segment, and T wave in all 12 leads, with amplitudes recorded to the nearest 100th of a millivolt and times recorded to the nearest millisecond.
We focused on global QRS duration as well as measures of left ventricular mass and ST segment changes. Sokolow-Lyon left voltage was calculated by addition of the amplitude of the S wave in lead V1 to the amplitude of the maximum R wave in lead V5 or V6.10 Cornell voltage was calculated by addition of the amplitude of the R wave in lead AVL to the amplitude of the S wave in lead V3.10,11 ST segment deviation at the J point, at the midpoint of the ST segment, and at the end of ST segment were extracted. ST slope was calculated as the difference between ST segment deviation at the end of the ST segment and at the J point. Presence or absence of left or right bundle-branch block was determined according to standard criteria.
End Points
The primary end point was time-related all-cause mortality relative to surgery, an end point which we have previously argued is a wholly objective, clinically relevant, and unbiased measure.12,13 Mortality was assessed by use of the Social Security Death Index. We have previously shown that among patients systematically followed in the Cardiovascular Information Registry, the Social Security Death Index has a sensitivity of 97%.14 High specificity of the Social Security Death Index has also been documented.15 A secondary end point was 30-day mortality.
Statistical Analyses
Continuous variables are described as median (interquartile range [IQR]) whereas categorical variables are described as number (percent). Cumulative mortality plots were constructed by use of the Kaplan-Meier method,16 with differences according to quartiles of different ECG measures tested by the log-rank
2 statistic.
Assessment of the association between quantitative ECG measures and time to death was performed with Cox proportional hazards regression.17 The proportional hazard assumption was tested by scaled Schoenfeld residuals as well as inspection of hazard ratio plots over time.18 Linear assumptions for continuous variables were relaxed by restricted cubic splines.18 Patients were grouped according to quintiles of predicted 10-year survival probabilities to assess model calibration and to compare predicted Kaplan-Meier survival to actual Kaplan-Meier survival; this process was repeated over 100 bootstrap resamplings. An additional set of survival analyses was performed with wholly parametric, time-decomposed, multiphase hazard regression as described previously19; these analyses yielded results that were not materially different from the Cox analyses. Analyses of the secondary end point, 30-day mortality, were performed with logistic regression with backward variable selection.
An ECG score was generated on the basis of the estimated restricted cubic spline function for the Cox model. Ideally the validity of this score would be determined by its application to an external data set. However, as we are not aware of any external data that contain the kind of detailed clinical and quantitative ECG information we have, we opted for an alternative approach, namely out-of-bagging, a process similar in concept to the leave-one-out bootstrap.20 We obtained bootstrap samples and used each sample to compute a prediction model that incorporated the ECG score and potential confounders; all demographic, clinical, and angiographic variables listed in Table 1 were included. Each bootstrap sample left out, on average,
37% of the data, which we refer to as the out-of-bag data. The prediction model was applied to the out-of-bag data to calculate an error rate, namely 1 – the c index.21 The c index is the generalized form of the c statistic for censored data21; despite the presence of heavy censoring in this index, which specifically incorporates censoring information, is arguably reasonable and reliable.22 Use of the c index as a measure of prediction error has been previously reported.23
|
To determine the relative importance of each variable, we recalculated the prediction error after randomly permutation of that variable in the out-of-bag data; a variable with a high degree of importance would be expected to yield a greater change in prediction error. The process was repeated 100 times to estimate relative importance values for the ECG score and other key covariates.
All analyses were performed with SAS version 9.1.3 (SAS Institute Inc., Cary, NC) or S-Plus 2000 Professional Software (Insightful, Inc, Seattle, Wash) with Harrells Design and Hmisc Libraries.18 All probability values were 2-sided and only considered significant if P<0.05.
The authors had full access to and take full responsibility for the integrity of the data. All authors have read and agree to the manuscript as written.
| Results |
|---|
|
|
|---|
ECG Findings and Mortality
During a median follow-up of 6 years, 1516 patients (19%) died. Unadjusted associations of ECG measures and mortality are shown in Figures 1 through 5![]()
![]()
![]()
. Correlates of increased mortality risk included longer QRS duration (Figure 1), higher Sokolow-Lyon voltage (Figure 2), higher Cornell voltage (Figure 3), greater ST segment depression (Figure 4), and greater negativity of ST segment slope (Figure 5).
|
|
|
|
|
Multivariable Analyses
After accounting for all the potential covariates listed in Table 1, higher heart rate, longer PR interval, longer QRS duration, greater Sokolow-Lyon voltage, and greater negativity of ST segment slope were independent ECG predictors of death (P<0.0001 for all) (Table 2). We found no statistically significant interactions between ECG measures or between ECG measures and a number of preselected clinical covariables.
|
ECG Score and Prognostic Validation
An adjusted ECG score was generated with the confounder-adjusted parameter coefficients for each of the independently predictive ECG variables (see Appendix for equation). A strong gradient was apparent whereby increased ECG score was strongly predictive of death (Figure 6), even after accounting for all the covariates listed in Table 1 (Figure 7). By out-of-bagging validation, ECG score was the second strongest predictor of death, exceeded only by age (Figure 8).
|
|
|
Thirty-Day Mortality
During the first 30 days after surgery, 111 deaths (1.4%) occurred. By both univariable and multivariable analyses, the only ECG predictor of 30-day mortality was the amount of ST-segment depression in lead V5 (by quartiles of most ST depression to least ST depression, mortality of 2.0%, 1.4%, 1.2%, and 0.9%; P for trend, 0.002). In a multivariable logistic regression model that included the variables listed in Table 1, the only independent predictors of death were intraaortic balloon pump use (adjusted odds ratio, 13.88; 95% confidence interval, 8.11 to 23.76; P<0.0001), no thoracic artery graft (adjusted odds ratio, 2.75; 95% confidence interval, 1.26 to 6.00; P=0.011), and amount of ST segment depression in lead V5 (adjusted odds ratio for 0.2 mV compared with none, 1.07; 95% confidence interval, 1.04 to 1.10; P<0.0001). The association of ST segment depression with 30-day mortality was markedly nonlinear, with a threshold effect noted at 0 mV; ST segment elevation was not predictive.
| Discussion |
|---|
|
|
|---|
The present results parallel recently reported findings among patients with chronic hypertension.3,7 As discussed elsewhere, ECG left ventricular hypertrophy is reflected by both QRS voltage and QRS duration.7 By taking advantage of digital quantitative data, we were able to show that ECG left ventricular hypertrophy is not a dichotomous phenomenon. Both by univariable and multivariable analyses, Sokolow-Lyon voltage predicted death according to a J-shape with higher risks at highest and lowest values (Figure 2), whereas Cornell voltage followed a more typical continuous pattern (Figure 3); the reasons for this are unclear and will require further investigation. Our findings that link negative ST segment depression, a marker of abnormal ventricular repolarization, with mortality also parallel previous reports based on population-based cohorts.1,6 We generated a composite confounder-adjusted ECG risk score and demonstrated its strong prognostic validity by use of the out-of-bagging technique; of note, except for age the ECG risk score was the strongest predictor of death (Figure 8). Taken together, the present study provides further evidence to support routine use of quantitative ECG measures for clinical risk stratification.
Some important limitations need to be acknowledged. Our population was a relatively healthy one for patients who underwent cardiac surgery and excluded patients with important valvar disease, atrial fibrillation, and prior Q wave myocardial infarction. Because of software limitations, we could not measure time-voltage areas, which have been shown to be strong correlates of left ventricular hypertrophy.25 We considered only baseline preoperative ECGs and did not systematically obtain routine follow-up measurements, as was done, for example, in the Losartan Intervention for Endpoint Reduction (LIFE) trial.3 Another arguable clinical limitation is that, unlike QRS duration, quantitative measures of Sokolow-Lyon voltage and ST segment depression and slope are not routinely reported on standard ECGs. However, current ECG software has the capability to report them. Although QRS duration and voltage and ST segment depression are predictive of long-term mortality, it is not known whether specific pharmacological intervention, such as the use of angiotensin receptor blockers, can modify outcome specifically among CABG patients.26 Recently, some patients with prolonged QRS intervals may have been treated with cardioverter-defibrillator and/or cardiac resynchronization therapy27; the beneficial effects of these treatments, if any, would have led to an underestimation of the importance of prolonged QRS duration as a predictor of long-term outcome in patients who undergo isolated CABG.
Despite these limitations, we found that routine preoperative quantitative ECG measures of heart rate, conduction, left ventricular mass, and repolarization are independently predictive of long-term mortality among patients who undergo isolated CABG. Future research will be needed to confirm these findings and to determine how best to use ECG measures for risk stratification and to modify long-term postoperative management.
| Appendix |
|---|
|
|
|---|
ECG score = 10 x (0.0020805218 x ST segment slope + 0.021994013 x ventricular rate – 1.1811993e–05 x MAX (ventricular rate – 53,0)3 + 1.9686654e–05 x MAX (ventricular rate – 67,0)3 – 7.8746617e–06 x MAX (ventricular rate – 88,0)3 – 0.0020579794 x PR interval + 8.2754258e–07 x MAX (PR interval – 136,0)3 – 1.4068224e–06 x MAX (PR interval – 164,0)3 + 5.7927981e–07 x MAX (PR interval – 204,0)3 + 0.0063538696 x QRSD duration – 0.00012443304 x Sokolow-Lyon voltage + 5.1974407e–11 x MAX (Sokolow-Lyon voltage – 1129,0)3 – 9.4513461e–11 x MAX (Sokolow-Lyon voltage – 2060,0)3 + 4.2539054e–11 x MAX (Sokolow-Lyon voltage – 3197.5,0)3)
| Acknowledgments |
|---|
The present work was supported by National Institutes of Health grants R01 HL-66004–2, R01 HL-072771–02, and P50 HL-77107–1 (to Drs Lauer and Blackstone).
Disclosures
None.
| References |
|---|
|
|
|---|
2. Levy D, Salomon M, DAgostino RB, Belanger AJ, Kannel WB. Prognostic implications of baseline electrocardiographic features and their serial changes in subjects with left ventricular hypertrophy. Circulation. 1994; 90: 1786–1793.
3. Okin PM, Devereux RB, Jern S, Kjeldsen SE, Julius S, Nieminen MS, Snapinn S, Harris KE, Aurup P, Edelman JM, Wedel H, Lindholm LH, Dahlof B. Regression of electrocardiographic left ventricular hypertrophy during antihypertensive treatment and the prediction of major cardiovascular events. JAMA. 2004; 292: 2343–2349.
4. Iuliano S, Fisher SG, Karasik PE, Fletcher RD, Singh SN. QRS duration and mortality in patients with congestive heart failure. Am Heart J. 2002; 143: 1085–1091.[CrossRef][Medline] [Order article via Infotrieve]
5. Kannel WB, Gordon T, Offutt D. Left ventricular hypertrophy by electrocardiogram. Prevalence, incidence, and mortality in the Framingham study. Ann Intern Med. 1969; 71: 89–105.
6. Daviglus ML, Liao Y, Greenland P, Dyer AR, Liu K, Xie X, Huang CF, Prineas RJ, Stamler J. Association of nonspecific minor ST-T abnormalities with cardiovascular mortality: the Chicago Western Electric Study. Jama. 1999; 281: 530–536.
7. Oikarinen L, Nieminen MS, Viitasalo M, Toivonen L, Jern S, Dahlof B, Devereux RB, Okin PM. QRS duration and QT interval predict mortality in hypertensive patients with left ventricular hypertrophy: the Losartan Intervention for Endpoint Reduction in Hypertension Study. Hypertension. 2004; 43: 1029–1034.
8. Shenkman HJ, Pampati V, Khandelwal AK, McKinnon J, Nori D, Kaatz S, Sandberg KR, McCullough PA. Congestive heart failure and QRS duration: establishing prognosis study. Chest. 2002; 122: 528–534.[CrossRef][Medline] [Order article via Infotrieve]
9. Lauer MS, Lytle B, Pashkow F, Snader CE, Marwick TH. Prediction of death and myocardial infarction by screening with exercise-thallium testing after coronary-artery-bypass grafting. Lancet. 1998; 351: 615–622.[CrossRef][Medline] [Order article via Infotrieve]
10. Sundstrom J, Lind L, Arnlov J, Zethelius B, Andren B, Lithell HO. Echocardiographic and electrocardiographic diagnoses of left ventricular hypertrophy predict mortality independently of each other in a population of elderly men. Circulation. 2001; 103: 2346–2351.
11. Casale PN, Devereux RB, Alonso DR, Campo E, Kligfield P. Improved sex-specific criteria of left ventricular hypertrophy for clinical and computer interpretation of electrocardiograms: validation with autopsy findings. Circulation. 1987; 75: 565–572.
12. Lauer MS, Blackstone EH, Young JB, Topol EJ. Cause of death in clinical research: time for a reassessment? J Am Coll Cardiol. 1999; 34: 618–620.
13. Lauer MS, Topol EJ. Clinical trials: multiple treatments, multiple end points, and multiple lessons. JAMA. 2003; 289: 2575–2577.
14. Nishime EO, Cole CR, Blackstone EH, Pashkow FJ, Lauer MS. Heart rate recovery and treadmill exercise score as predictors of mortality in patients referred for exercise ECG. JAMA. 2000; 284: 1392–1398.
15. Newman TB, Brown AN. Use of commercial record linkage software and vital statistics to identify patient deaths. J Am Med Inform Assoc. 1997; 4: 233–237.
16. Kaplan E, Meier P. Nonparametric estimation from incomplete observations. J Am Stat Assoc. 1958; 53: 457–481.[CrossRef]
17. Cox D. Regression models and life tables. J R Stat Soc. 1972; 34: 187–220.
18. Harrell FE Jr. Regression Modeling Strategies: With Applications to Linear Models, Logistic Regression, and Survival Analysis. New York: Springer-Verlag; 2001.
19. Blackstone EH, Naftel DC, Turner MEJ. The decomposition of time-varying hazard into phases, each incorporating a separate stream of concomitant information. J Am Stat Assoc. 1986; 81: 615–624.[CrossRef]
20. Efron B, Tibshirani R. Improvements on cross-validation: the 0.632+bootstrap method. J Am Stat Assoc. 1997; 92: 548–560.[CrossRef]
21. Harrell FE, Jr., Califf RM, Pryor DB, Lee KL, Rosati RA. Evaluating the yield of medical tests. JAMA. 1982; 247: 2543–2546.
22. May M, Royston P, Egger M, Justice AC, Sterne JA. Development and validation of a prognostic model for survival time data: application to prognosis of HIV positive patients treated with antiretroviral therapy. Stat Med. 2004; 23: 2375–2398.[CrossRef][Medline] [Order article via Infotrieve]
23. Kattan MW, Hess KR, Beck JR. Experiments to determine whether recursive partitioning (CART) or an artificial neural network overcomes theoretical limitations of Cox proportional hazards regression. Comput Biomed Res. 1998; 31: 363–373.[CrossRef][Medline] [Order article via Infotrieve]
24. Holmvang L, Clemmensen P, Lindahl B, Lagerqvist B, Venge P, Wagner G, Wallentin L, Grande P. Quantitative analysis of the admission electrocardiogram identifies patients with unstable coronary artery disease who benefit the most from early invasive treatment. J Am Coll Cardiol. 2003; 41: 905–915.
25. Okin PM, Roman MJ, Devereux RB, Pickering TG, Borer JS, Kligfield P. Time-voltage QRS area of the 12-lead electrocardiogram: detection of left ventricular hypertrophy. Hypertension. 1998; 31: 937–942.
26. Okin PM, Devereux RB, Jern S, Kjeldsen SE, Julius S, Nieminen MS, Snapinn S, Harris KE, Aurup P, Edelman JM, Dahlof B. Regression of electrocardiographic left ventricular hypertrophy by losartan versus atenolol: the Losartan Intervention for Endpoint reduction in Hypertension (LIFE) Study. Circulation. 2003; 108: 684–690.
27. Bristow MR, Saxon LA, Boehmer J, Krueger S, Kass DA, De Marco T, Carson P, DiCarlo L, DeMets D, White BG, DeVries DW, Feldman AM. Cardiac-resynchronization therapy with or without an implantable defibrillator in advanced chronic heart failure. N Engl J Med. 2004; 350: 2140–2150.
![]() |
E. Hsich, E. Z. Gorodeski, R. C. Starling, E. H. Blackstone, H. Ishwaran, and M. S. Lauer Importance of Treadmill Exercise Time as an Initial Prognostic Screening Tool in Patients With Systolic Left Ventricular Dysfunction Circulation, June 30, 2009; 119(25): 3189 - 3197. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. W. Magnani, M. A. Williamson, P. T. Ellinor, K. M. Monahan, and E. J. Benjamin P Wave Indices: Current Status and Future Directions in Epidemiology, Clinical, and Research Applications Circ Arrhythmia Electrophysiol, February 1, 2009; 2(1): 72 - 79. [Full Text] [PDF] |
||||
![]() |
B. B. Pavri, M. B. Hillis, H. Subacius, G. E. Brumberg, A. Schaechter, J. H. Levine, A. Kadish, and for the Defibrillators in Nonischemic Cardiomyopat Prognostic Value and Temporal Behavior of the Planar QRS-T Angle in Patients With Nonischemic Cardiomyopathy Circulation, June 24, 2008; 117(25): 3181 - 3186. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2007 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |