| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Circulation. 2000;102:1252.)
© 2000 American Heart Association, Inc.
Clinical Investigation and Reports |
From the Cardiology Division (M.Z.), Klinikum Benjamin Franklin, Free University, Berlin, Germany; Cardiological Sciences (B.A., M.M.), St Georges Hospital Medical School, London, England; Cardiology Division (T.K., S.H.H.), J.W. Goethe University, Frankfurt, Germany; and Cardiology Division (M.R.F.), VA Medical Center and Georgetown University, Washington, DC.
Correspondence to Marek Malik, PhD, MD, Cardiological Sciences, St Georges Hospital Medical School, London SW17 0RE, UK. E-mail m.malik{at}sghms.ac.uk
| Abstract |
|---|
|
|
|---|
Methods and ResultsIn 280 consecutive post-MI patients, a 12-lead ECG was recorded before discharge, optically scanned, and digitized. For the present study, 5 T-wave morphology descriptors were automatically calculated after singular value decomposition of the ECG signal. The total cosine R-to-T (TCRT [describes the global angle between repolarization and depolarization wavefront]) and the T-wave loop dispersion were univariately associated (P=0.0002 and P<0.002, respectively, U test) with 27 prospectively defined clinical events in 261 patients (mean follow-up 32±10 months). Kaplan-Meier event probability curves for strata above and below the median confirmed the strong risk discrimination by TCRT and T-wave loop dispersion (P<0.003 and P<0.001, respectively, log-rank test). On Cox regression analysis, with the entering of age, left ventricular ejection fraction, heart rate, QRS width, reperfusion therapy, ß-adrenergicblocker treatment, and standard deviation of R-R intervals on 24-hour Holter monitoring, TCRT (P<0.03) yielded independent predictive value, whereas T-wave loop dispersion was of borderline independence (P=0.064). Heart rate (P<0.02), left ventricular ejection fraction (P<0.02), and reperfusion therapy (P<0.02) also remained in the final model.
ConclusionsComputerized T-wave morphology analysis of the 12-lead resting ECG permits independent assessment of post-MI risk and an improved risk stratification when combined with other risk markers.
Key Words: myocardial infarction death, sudden risk factors waves electrocardiography
| Introduction |
|---|
|
|
|---|
We hypothesized that composite temporospatial measures of ECG repolarization dispersion would improve risk stratification for postmyocardial infarction (MI) ventricular arrhythmias and sudden cardiac death in comparison with stratifiers such as QTD, left ventricular ejection fraction (LVEF), or thrombolytic therapy. In particular, we applied recently developed novel quantitative repolarization analyses that are accessible from a single ECG beat.11 The total cosine R-to-T (TCRT) reflects the spatial angle between depolarization and repolarization, akin to the venerable concept of the ventricular gradient.12 T-wave loop dispersion extends this concept, reflecting variability of the T-wave vector loop. The normalized T-wave loop area measures heterogeneity of principal components of the T wave within its loop, whereas T-wave morphology dispersion expresses morphological heterogeneity within the 12-lead ECG. Conceptually similar but distinct examinations of repolarization complexity were recently successful in stratifying arrhythmic risk in the long QT syndrome13 and arrhythmogenic right ventricular dysplasia.14
We therefore set out to determine the usefulness of our analyses in temporospatial dispersion for the stratification of arrhythmic risk and death in a series of prospectively studied post-MI patients.15 This analysis resulted in a surprisingly powerful accuracy of the prediction of cardiac mortality from the 12-lead resting ECG.
| Methods |
|---|
|
|
|---|
T-Wave Morphology Descriptors
Analysis of the digital ECG recordings was
performed in a fully automatic manner with a custom-developed software
implemented on a personal computer.11 The study was
conducted in a strictly blinded manner: digital ECGs were sent from
Germany to London for processing without any clinical data. The
analytical results were returned to the German center as an unlabeled
numerical spreadsheet. Only after completion of the statistical
analyses were the German center workers given the description
of the analytical procedures.
The analysis program performs a singular value decomposition of the ECG signal into a minimum dimensional space. From singular value components, principal component analysis (PCA) was performed as recently described.13 14 Complexity ratio (CR) was the ratio of the singular value of the second most significant component to the square root of the sum of the squares of all 8 singular values.
Based on the decomposition, several descriptors were calculated of
spatial and temporal variations of T-wave morphology and repolarization
wavefront direction11 (Figure 1
). (1) The so-called T-wave loop
dispersion measures the variation of the ECG vector (ie, the variation
of the interlead relations among domain of interlead relations spanned
by the ECG vector). This variable is unitless; its maximum value is
100. (2) The so-called normalized T-wave loop area describes the shape
and irregularity of the T-wave loop by expressing its area as a
fraction of the rectangle that encompasses the loop. The variable
is unitless. (3) The so-called TCRT measures the vector deviation
between the depolarization and repolarization waves by calculating
cosine values between the 3-dimensional R- and T-wave loop vectors
within the optimized decomposition space. Negative values correspond to
large differences in the orientation of the 2 loops. The variable
is unitless. (4) The so-called T-wave morphology dispersion expresses
the dissimilarities between the T-wave shapes in individual leads,
based on the differences between reconstruction vectors of individual
ECG leads created from the 3-dimensional T-wave loop. It is calculated
as the average of angles between all possible pairs of reconstruction
vectors. A small value indicates that reconstruction vectors are
close to each other, indicating similar T-wave morphology between
leads.
|
The fully automatic processing ensures 100% reproducibility of all variables for any given ECG. For serial ECGs in the same subject, a high reproducibility, of up to 99.7%, was previously reported.11
Follow-Up
Follow-up data and information on clinical end points were
unchanged from the initial publication of the same post-MI patient
group.15 The prospectively defined primary end point
combined all-cause mortality, sustained VT, and resuscitated VF. The
secondary end points consisted of arrhythmic events: sudden cardiac
death, documented sustained VT, and resuscitated VF.
Statistical Analysis
Statistical evaluation was processed independent of the T-wave
analyses. For comparison, the results for conventional QTD
variables (QRS width, QTD, JT dispersion, T peak-to-end interval,
area under the T wave), for other risk stratifiers (SDNN from Holter
recordings, LVEF), and for clinical variables (age,
reperfusion therapy, ß-blocker treatment, heart rate from the study
ECG) were taken from the previous study of this
cohort.15
Data were analyzed with SPSS Version 7.0 for Windows for Cox
regression analyses and JMP-3.1 software (SAS Institute) for
all other statistics. Continuous values are reported as mean±SD.
Comparisons between patients with and without events during follow-up
were performed by the nonparametric U test.
Pearsons correlations between ECG and clinical variables were
used. The relation of ECG variables to categorical clinical
variables was tested by a
2 test.
Kaplan-Meier event probability curves were computed, with patient
groups stratified according to the median value of the respective
variable. The cumulative probability of events of 2 patient groups
was compared by the log-rank test. The independent correlation of
multiple variables with the timing of events during follow-up as
the dependent variable was determined by Cox regression
analysis. Because of the strong correlation of left
bundle-branch block (LBBB) with TCRT, all statistical analyses
were repeated with 9 patients with LBBB excluded. Statistical
significance was considered for P<0.05.
| Results |
|---|
|
|
|---|
|
T-Wave Morphology Analysis: Correlations With Conventional
QTD Variables and Clinical Data
Pearsons correlation coefficients were determined for pairs of
T-wave morphology variables. All r values were <0.25
and thus of no clinical relevance. Specifically, TCRT and T-wave loop
dispersion were unrelated (r=0.11, P=NS). For
pairs with conventional QTD variables and QRS width, an
intermediate inverse relationship with the average area under the T
wave was found for CR (r=-0.40, P<0.001).
Patients with LBBB (n=9) had TCRT of -0.77±0.25, whereas patients
without LBBB (n=252) had an average TCRT of 0.11±0.12
(P=0.005). T-wave loop dispersion exhibited a borderline
difference related to LBBB (33±4 with LBBB versus 35±5 without LBBB,
P=0.044), whereas all other T-wave variables were not
influenced by LBBB. All other correlations were not practically
relevant, including several significant r values between
0.17 and 0.25. Right bundle-branch block was not associated with the
T-wave morphology variables. Moreover, none of the continuous
clinical variables, including heart rate, age, SDNN from Holter,
and LVEF, were related to the T-wave morphology variables. This was
also true for other clinical variables, such as sex, use of
reperfusion therapy, ß-blocker treatment, or infarct location.
T-Wave Morphology Analysis: Univariate
Prognostic Information
Of 280 ECGs, 19 were excluded due to insufficient data quality or
missing leads. Tables 2
and 3
summarize the values for T-wave
morphology variables for 27 patients with events compared with 234
event-free patients. Tables 2
and 3
also show the same
comparison between 17 patients with and 244 patients without arrhythmic
events. CR exhibited paradoxically lower rather than theoretically
expected higher values in the event group, with a borderline
significance (P<0.07). The comparison for patients with and
without arrhythmic events was not significant. TCRT was lower in the
patient group with events (P<0.0002) and in patients with
arrhythmic events (P<0.004). Similarly, T-wave loop
dispersion was lower in patients with primary end points
(P<0.002) and arrhythmic events (P<0.003). No
difference was found for T-wave morphology dispersion and normalized
T-wave loop area.
|
|
T-Wave Morphology Analysis: Kaplan-Meier Event
Probabilities
Patients with a TCRT below the median had decreased survival
(P<0.003, Figure 2A
) and a
higher incidence of arrhythmic events (P<0.007). The median
value of T-wave loop dispersion also discriminated patients at risk of
overall end points (P<0.001, Figure 2B
) and
arrhythmic events (P<0.004). CR, T-wave morphology
dispersion, and normalized T-wave loop area did not result in different
Kaplan-Meier event probabilities. In an analysis of the
patients without LBBB (n=252), the results were similar with TCRT
(P=0.022 and P<0.04, respectively, for primary
and arrhythmic end points) and T-wave loop dispersion
(P=0.003 and P=0.014, respectively).
|
Multivariate Analysis of Risk
Factors
Clinical variables (age, reperfusion therapy, LVEF,
ß-blocker treatment, QRS width, heart rate), Holter
parameters (SDNN), and T-wave morphology descriptors (TCRT,
T-wave loop dispersion) that were univariately predictive
of follow-up end points were entered as independent variables into
a Cox regression model with stepwise backward removal (Tables 4
and 5
).
At least 1 T-wave morphology variable remained in the equation at
the last regression step, namely, TCRT for the prediction of primary
end points and T-wave loop dispersion for the prediction of arrhythmic
events. By entering only univariately predictive T-wave
morphology variables into the Cox regression model (TCRT, T-wave
loop dispersion, normalized T-wave loop area) and by considering only
patients without LBBB (n=252), we showed that T-wave loop dispersion
(P=0.009 and P=0.003, respectively) and TCRT
(P=0.006 and P=0.07, respectively) were
independently predictive of primary end points and arrhythmic
events.
|
|
| Discussion |
|---|
|
|
|---|
1 T-wave morphology descriptor remained in the
model, adding independent information to other risk stratifiers such as
heart rate, LVEF, and the administration of reperfusion therapy.
Comparison With PCA and Conventional QTD
In a comparison with PCA of the T wave,13 14 2 of the
proposed T-wave morphology descriptors were superior. At best, PCA
yielded borderline significance in the univariate
comparison for the prediction of overall events. In contrast, 2 T-wave
morphology variables were strong risk predictors for overall or
arrhythmic events, confirmed by the calculation of Kaplan-Meier curves.
This strong prognostic value contrasts with the disappointing results
with conventional QTD variables13 from the same ECGs.
A major aspect of the new variables is the technical improvement
over QTD in both reproducibility and measurement objectivity.
Importantly, the proposed algorithms do not rely on an accurate T-wave
offset measurement, which is the major source of inaccuracy in
automatic QTD measurements.16 17 18 In most studies,
however, QTD was analyzed manually, which is even less reliable
due to subjective variations in waveform assessment between centers.
The intrasubject reproducibility is also improved.11
Pathophysiological Concept of Novel T-Wave
Morphology Descriptors
The prognostic content of abnormal repolarization had been
suspected because the early experiments demonstrated the role of
ventricular repolarization pathologies in the
arrhythmogenesis.1 19 20 21 Initially, this concept was
followed with the use of body surface potential
mapping.21 22 23 24 QTD from 12-lead ECGs was proposed as a
more practical surrogate.25 After the initial
enthusiastic reports,5 6 25 26 27 meticulous methodological
studies16 28 and reviews28 29 30 cast doubt on
the true value of QTD. Finally, the first truly prospective study in a
large post-MI population established that none of the conventional QTD
variables carried post-MI prognosis.15 Although this
finding is contradicted by recent epidemiological
studies,7 8 a more sophisticated concept of repolarization
analysis for the purpose of risk stratification is clearly
required.31 It was therefore attempted to explore more
accurate repolarization qualities from 12-lead surface ECGs, such as
PCA13 14 and T-wave loop morphology,32 or to
explore dynamic QT changes on Holter monitoring.4 33 34
None of these newer repolarization markers have been tested in a
post-MI population. Although PCA had been shown to be predictive in
long QT syndrome13 and in patients with arrhythmogenic
right ventricular dysplasia,14 it was not
found to be very useful in the present study. Beyond the approach
of PCA, Acar et al11 developed a set of novel T-wave
morphology descriptors to quantify various abnormal temporospatial
repolarization indices. With an analysis of these novel T-wave
descriptors in a strictly blinded manner, the predictive usefulness of
TCRT and T-wave loop dispersion was established in this study.
TCRT reflects a 3-dimensional comparison between repolarization and depolarization wavefronts and therefore is akin to the ventricular gradient introduced by Wilson et al35 in 1934 and later expanded with QRST area distributions20 21 22 23 24 from body surface mapping. Importantly, with TCRT, the amplitude of depolarization and repolarization processes is not considered, but rather only the difference in their direction is considered and thus represents a modification of the earlier concepts. The variable was found to be associated with LBBB but not with a prolonged QRS duration alone. With or without consideration of LBBB, it proved to be a powerful risk predictor in the present study.
When values for the present post-MI population are compared with the reported normal values,11 increased T-wave morphology dispersion and lowered values for TCRT, T-wave loop dispersion, and normalized T-wave area are pathological as such and differentiate event from nonevent patients. Although a low TCRT relates to a large deviation between the QRS and T-wave loops, a large T-wave morphology dispersion is associated with increased dissimilarities among the T waves in ECG leads, and a low value of normalized T-wave area as well as T-wave loop dispersion may be explained as arrhythmogenic with a pathologically compressed and narrowed T-wave loop. Furthermore, irregularity of the loop and deviation from an ellipsoid would also result in a decreased normalized T-wave area. Beyond these considerations and due to the novelty of the algorithms used in the present study, they lack a more detailed experimental basis of the electrophysiological mechanisms involved. Although basic research is under way, other approaches to noninvasive risk stratification of arrhythmia substrates, such as PCA, T-wave alternans, or QT variability, likewise fell short of a precise pathophysiological explanation during their clinical evaluation.
Importantly, only very weak relationships between the various novel descriptors and conventional variables of QTD were found, demonstrating that the proposed variables assess as yet undetected qualities of repolarization and do not reproduce or refine the more conventional measurements. Specifically, we could not confirm the result of the study by Kors et al,36 which correlated T-wave loop characteristics with QTD in a database of 1220 ECGs.
Comparison With Other Risk Markers
In testing the independent contribution of T-wave morphology
descriptors by means of Cox regression, at least 1 T-wave morphology
marker provided independent risk stratification. LVEF, heart rate,
reperfusion therapy, and, finally, TCRT or T-wave loop dispersion; all
added independently to the risk prediction of overall events.
Study Limitations
Although a large post-MI patient population was studied, the
overall number of general end points as well as arrhythmic events was
low, which reflects the modern treatment strategies. In particular, the
statistics with regard to arrhythmic end points should be viewed with
this limitation the number of events in more strictly defined
subcategories is not sufficient for further detailed comparisons.
The present study was the first attempt to apply detailed T-wave morphology descriptors in a set of prospectively collected ECG data. Consequently, we have used the same algorithmic settings as in the technical report by Acar et al.11 It is likely that for the purposes of post-MI risk stratification, the original algorithms that quantify T-wave morphology abnormalities should be further refined. It is also possible that the values of the T-wave morphology variables underlie dynamic changes after MI, particularly during the first year. However, such an influence is not expected to yield a false-positive relationship with outcome but rather would weaken the predictive value of the ECG descriptors.
In terms of other risk variables, we restricted our present investigation to the set of risk factors used in the previous report on this population.15 Other variables, such as recently described factors of heart rate turbulence,37 were not considered. It is unlikely, however, that heart rate turbulence would be pathophysiologically linked to T-wave morphological abnormalities.
Implications
The present results are important for future risk
stratification studies such as those of prophylactic
post-MI treatment. The risk stratifiers proposed represent 2
easily accessible measurements available from a single surface ECG
recording. Technically, the algorithms can be incorporated into
commercially available ECG recorders with digital capabilities.
Additional studies are warranted that investigate the usefulness of
T-wave morphology descriptors in other populations at risk of sudden
cardiac death.
Received January 28, 2000; revision received April 12, 2000; accepted April 14, 2000.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C.-Y. Lin, L.-Y. Lin, and P.-C. Chen Analysis of T-wave morphology from the 12-lead electrocardiogram for prediction of long-term prognosis in patients initiating haemodialysis Nephrol. Dial. Transplant., September 1, 2007; 22(9): 2645 - 2652. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. D. Critchley, P. Taggart, P. M. Sutton, D. R. Holdright, V. Batchvarov, K. Hnatkova, M. Malik, and R. J. Dolan Mental stress and sudden cardiac death: asymmetric midbrain activity as a linking mechanism Brain, January 1, 2005; 128(1): 75 - 85. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. M. Okin QT interval prolongation and prognosis: further validation of the quantitative approach to electrocardiography J. Am. Coll. Cardiol., February 18, 2004; 43(4): 572 - 575. [Full Text] [PDF] |
||||
![]() |
P. M. Okin, R. B. Devereux, E. T. Lee, J. M. Galloway, and B. V. Howard Electrocardiographic Repolarization Complexity and Abnormality Predict All-Cause and Cardiovascular Mortality in Diabetes: The Strong Heart Study Diabetes, February 1, 2004; 53(2): 434 - 440. [Abstract] [Full Text] |
||||
![]() |
P. Smetana, V. N. Batchvarov, K. Hnatkova, A. J. Camm, and M. Malik Ventricular gradient and nondipolar repolarization components increase at higher heart rate Am J Physiol Heart Circ Physiol, January 1, 2004; 286(1): H131 - H136. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. V. Huikuri, T. H. Makikallio, M. J. P. Raatikainen, J. Perkiomaki, A. Castellanos, and R. J. Myerburg Prediction of Sudden Cardiac Death: Appraisal of the Studies and Methods Assessing the Risk of Sudden Arrhythmic Death Circulation, July 8, 2003; 108(1): 110 - 115. [Full Text] [PDF] |
||||
![]() |
I. Kardys, J. A Kors, I. M van der Meer, A. Hofman, D. A.M van der Kuip, and J. C.M. Witteman Spatial QRS-T angle predicts cardiac death in a general population Eur. Heart J., July 2, 2003; 24(14): 1357 - 1364. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Geelen, I. A. Brouwer, P. L. Zock, J. A. Kors, C. A. Swenne, M. B. Katan, and E. G. Schouten (n-3) Fatty Acids Do Not Affect Electrocardiographic Characteristics of Healthy Men and Women J. Nutr., October 1, 2002; 132(10): 3051 - 3054. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Malik and V. N. Batchvarov Measurement, interpretation and clinical potential of QT dispersion J. Am. Coll. Cardiol., November 15, 2000; 36(6): 1749 - 1766. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Smetana, V. N. Batchvarov, K. Hnatkova, A. J. Camm, and M. Malik Sex differences in repolarization homogeneity and its circadian pattern Am J Physiol Heart Circ Physiol, May 1, 2002; 282(5): H1889 - H1897. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. M. Okin, R. B. Devereux, R. R. Fabsitz, E. T. Lee, J. M. Galloway, and B. V. Howard Principal Component Analysis of the T Wave and Prediction of Cardiovascular Mortality in American Indians: The Strong Heart Study Circulation, February 12, 2002; 105(6): 714 - 719. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Zabel, M. Malik, K. Hnatkova, V. Papademetriou, A. Pittaras, R. D. Fletcher, and M. R. Franz Analysis of T-Wave Morphology From the 12-Lead Electrocardiogram for Prediction of Long-Term Prognosis in Male US Veterans Circulation, March 5, 2002; 105(9): 1066 - 1070. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2000 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |