- Research article
- Open Access
The association of spatial T wave axis deviation with incident coronary events. The ARIC cohort
BMC Cardiovascular Disorders volume 5, Article number: 2 (2005)
Although current evidence suggests that the spatial T wave axis captures important information about ventricular repolarization abnormalities, there are only a few and discordant epidemiologic studies addressing the ability of the spatial T wave axis to predict coronary heart disease (CHD) occurrence.
This prospective study analyzed data from 12,256 middle-aged African American and white men and women, from the Atherosclerosis Risk in Communities Study (ARIC). Following a standardized protocol, resting standard 12-lead, 10-second electrocardiograms were digitized and analyzed with the Marquette GE program. The median follow-up time was 12.1 years; incident coronary heart disease comprised fatal and non-fatal CHD events.
The incidence rate of CHD was 4.26, 4.18, 4.28 and 5.62 per 1000 person-years respectively, across the spatial T wave axis quartiles. Among women for every 10 degrees increase in the spatial T wave axis deviation, there was an estimated increase in the risk of CHD of 1.16 (95% CI 1.04–1.28). After adjustment for age, height, weight, smoking, hypertension, diabetes, QRS axis and minor T wave abnormalities, this hazard rate ratio for women fell to 1.03 (0.92–1.14). The corresponding crude and adjusted hazard ratios for men were 1.05 (95% CI 0.96–1.15) and 0.95 (0.86–1.04) respectively.
In conclusion, this prospective, population-based, bi-ethnic study of men and women free of coronary heart disease at baseline shows that spatial T wave axis deviation is not associated with incident coronary events during long-term follow up. It is doubtful that spatial T wave axis deviation would add benefit in the prediction of CHD events above and beyond the current traditional risk factors.
Ventricular repolarization abnormalities play an important role in the determination of arrhythmia and sudden cardiac death [1, 2] and may reflect subclinical myocardial ischemia changes. The process of repolarization at rest is routinely quantified from standard 12-lead electrocardiogram (ECG), either as time-domain indexes such as the QT interval and its derivations, or as abnormalities of the ST segment or of the T wave. While these indexes have all been reported to be associated to some degree with incident coronary heart disease events, they have some limitations. The widely used QT interval reflects only the temporal aspect of the repolarization, and as currently defined and measured from the 12-lead ECG, QT dispersion has major conceptual and technical limitations [3, 4].
Theoretical and experimental studies suggest that ventricular repolarization occurs in a nonlinear and inhomogeneous fashion [5–10]. As a consequence, spatial measures of repolarization that take into account T-wave complexity using the T-wave vector (axis) should be more accurate and useful surface ECG markers of repolarization abnormalities than simple scalar intervals from the ECG, such as the QT interval or QT dispersion [11–16].
Clinical studies have shown that the T wave axis reflects changes associated with autonomic adaptive or maladaptive influences , systemic hypertension [18, 19], coronary occlusion , and microalbuminuria in individuals without diabetes mellitus . Moreover, it has recently been shown that spatial T wave axis deviation has good measurement properties and is repeatable , a finding that supports its use in clinical and epidemiological research.
Two prospective population-based reports on T wave axis deviation measured from standard 12-lead ECGs in older populations (the Rotterdam study , and Cardiovascular Health Study (CHS)  suggest that it is an indicator of increased risk of coronary heart disease and total mortality, independent of other cardiovascular risk factors. In the cohort of high-risk, middle-aged men from the Multiple Risk Factor Intervention Trial (MRFIT), baseline spatial T wave axis deviation was not significantly associated with incident coronary events, although the change over time in the spatial T wave axis deviation was reported to be associated with incident events on long-term follow-up . These studies suggest that spatial T wave axis deviations capture changes in the ventricular repolarization process that are of potential clinical and epidemiologic importance.
The aim of this study was to assess whether a single base-line measurement of spatial T wave axis deviation is associated with prospectively ascertained coronary events in a population-based, bi-ethnic cohort of middle-aged men and women, to determine if such an association is independent of other risk factors, and to compare its association with that of other indexes of altered repolarization, such as a prolonged QT interval and minor T wave abnormalities.
Between 1987 and 1989, the Atherosclerosis Risk in Communities (ARIC) study examined population-based samples of residents aged 45 to 64 years from 4 communities in North Carolina, Mississippi, Minnesota, and Maryland. From the initial sample of ARIC baseline participants (n = 15,792) we excluded those not having an ECG (n = 205), and because of small numbers, self-identified ethnicity other than white or African American (n = 48). Further, we excluded participants likely to present secondary repolarization abnormalities such as myocardial infarction (self-reported or ECG-diagnosed, n = 790), use of digitalis or anti-arrhytmic drugs (n = 115), or pathologies likely to alter the measurement or interpretation of the spatial T wave axis as identified by the Minnesota code (MC), such as major Q or QS waves (MC 1.1, 1.2), ST depression or elevation (MC 4.1 to 4.4 and MC 9.2), negative T waves (MC 5.1 or 5.2), WPW pattern (MC 6.6), ventricular conduction defects (MC 7.1, 7.2, 7.4, and atrial fibrillation or flutter (MC 8.3), (n = 2093). Participants who had ECG evidence or history of myocardial infarction, coronary bypass, or angioplasty) (n = 162) at baseline were excluded from analysis. Participants whose angina status was positive or unknown by the Rose questionnaire were not excluded because the questionnaire's validity, especially in women, has been questioned . The final study sample was 12,256 participants.
ECG measurements, processing, and definition of the spatial T wave axis
During a fifteen-minute supine rest, trained and certified technicians positioned disposable Ag/AgCl electrodes and recorded a supine, resting 12-lead, 10-second electrocardiogram using a MAC Personal Cardiographer™ (Marquette Electronics, Inc., Jupiter, FL). Detailed procedures used for electrode placement and skin preparations are described in the Operation Manual . ECGs were digitized at 250 Hz and sent daily via modem to the Epidemiological Cardiology Research (EPICARE) Center, North Carolina. The EPICARE Center, blinded to participant identity processed the ECGs using the most recent version of the Marquette GE program, version 12SL.
Quality assurance procedures were designed specifically for ECG acquisition and processing . These programs included: standardized electrode positions were assured by marking of the skin using a standardized flexible ruler according to a detailed procedure manual. Electrical noise, overall and beat-to-beat electrocardiogram drift were monitored and quality-scored using a five level grading system to assist technicians in identifying unacceptable electrocardiograms and maintaining quality. Technical performance criteria were established and monitored by the EPICARE Center. All technicians were specifically trained and certified. All ECGs were read according to the Minnesota code  without knowledge of clinical or demographic data. All resting 12 lead ECGs coded as having key abnormalities by Minnesota codes including any 1-code, any 1-2 code or 2-2, 5-1 or 2 or any 9-2, 6-4, 7-1-1 or 7-2-1 codes and a 10% random sample of all other ECGs read by computer also were visually coded at the Minnesota ECG Reading Center. Adjudication of discrepancies was performed at the Minnesota Center. The computer-assigned codes are used as study data except where adjudication resulted in a code different from the original computer codes.
Spatial T-wave axis was calculated as previously described  from integrated T-wave amplitudes of the XYZ leads. Briefly, the inverse transformation by Dower et al.  was used to derive Frank's XYZ leads, with the polarity of the Z lead inverted in order to generate QRS and T wave patterns with a more familiar waveform in that lead. Thus, the positive direction of the Z axis is in the anterior direction. Spatial T-wave axis was calculated from the scalar product between the T vector and a unit vector in a normal reference direction (x = 1/√3, y = 1/√3, and z = -1/√3, where x, y and z are the unit vector components in the X, Y and Z directions). T axis expresses spatial T vector deviation from the approximate normal direction of the T vector 45° anteriorly in the XZ plane and at 45° elevation from the Y axis.
The repeatability of the spatial T wave axis deviation was determined using identical study procedures in a group of apparently healthy volunteers with demographic characteristics similar to the ARIC participants. The repeatability of the spatial T wave axis deviation was high with an intra-class correlation coefficient of 0.87 . The Cornell voltage index was determined as the sum of R wave amplitude in leads aVL and the S wave amplitude in lead V3 . The QT interval from the digital 12-lead ECG was determined by the NOVACODE program . An overall QT interval was calculated from the common QRS onset and T offset for all 12 leads together. To attempt correction for the heart rate dependence, we used several approaches: Bazett's and Fridericia's formulae as well as the QT Index. For comparability with other studies we report the QT interval "corrected "by the Bazett's formula , in spite of its well-known limitations.
The baseline examination was carried out in 1987–1989 and consisted of a home interview of all potential cohort members including items on cardiovascular risk factors, socioeconomic factors and family medical history. After obtaining informed consent, the clinic examination consisted of medical history interview, blood pressure and anthropometric measurements, venipuncture for blood samples and a 12-lead standard electrocardiogram. Anthropometrics were measured with participants wearing scrub suits and no shoes. Height, measured to the nearest centimeter and weight, measured to the nearest pound, were used to calculate body mass index (BMI). Waist-to-hip ratio was defined as the waist girth at the umbilicus (centimeters) divided by the maximum girth of the hips (centimeters). Smoking status was defined as "current smoker" if the person answered "yes" to both of the following questions: "Have you ever smoked cigarettes? and " Do you now smoke cigarettes?". Cholesterol and triglycerides were measured enzymatically. LDL-cholesterol was calculated. Hypertension was defined as having systolic blood pressure values equal or higher than 140 mmHg, or diastolic blood pressure values equal or higher than 90 mmHg or use of blood pressure lowering medication use in the past two weeks .
Ascertainment and classification of incident CHD cases
Deaths and hospitalization events were ascertained by annual follow-up calls to the cohort members, review of vital records, and community surveillance of hospitalized and fatal events. CHD death was defined as lacking a probable non-CHD cause, and occurring in the context of a recent myocardial infarction, chest pain within 72 hours of death, or a history of CHD. Events were classified independently by two members of the Mortality and Morbidity Classification Committee and discrepancies were adjudicated by a third member. Descriptions of event ascertainment and classification have been published [33, 34]. For the present study, we included CHD events occurring between the ARIC baseline examination and December 31, 2000. The median follow-up time was 12.1 years (maximum of 14.1). We defined CHD incidence as (1) a definite or probable MI, (2) a silent MI between examinations ascertained by ECG, or (3) a definite CHD death.
We treated spatial T wave axis deviation as a continuous variable. For the purpose of comparison, however we categorized the spatial T wave axis deviation using thresholds employed in previously published reports: normal (<30 degrees), borderline (30 to 45 degrees), and abnormal deviation (>=45 degrees). Cox proportional hazards regression  was used to model time-to-event in the presence of censoring. Time-at-risk (time-to-event or time-to-censoring) was calculated from the date of the baseline examination to the earliest of the following: date of coronary event, date of death, date of last follow-up contact, or December 31, 2000. Hazard ratios with 95%confidence intervals were derived for each variable. When treating spatial T wave axis as a categorical variable, dummy variables were created. Interactions between spatial T wave axis and gender, QRS axis, minor T wave abnormalities and QT interval, were tested as cross-product terms in the multivariable model. Significance of interaction was assessed based on the likelihood ratio test (p < 0.10). Statistical analyses were conducted for each gender to allow for well-known differences in ventricular repolarization by gender, and because of a statistically significant interaction between spatial T axis and gender on CHD events. SAS version 8.1 (SAS Institute, Cary, North Carolina) was used in all analyses.
The study population consisted of 12,256 participants, of which 7,143 (58.3%) were women and 3,050 (24.9%) were African American. The mean (SD) age of the study population at baseline was 53.8 (5.7) years. The cohort was followed up for a maximum of 14.1 years [mean (SD) 11.6 (2.3) years]. During this time 653 incident coronary heart disease events occurred. Of the total number of events, 250 coronary heart disease events occurred among women (cumulative incidence 3.5%) and 403 (cumulative incidence 7.9%) among men. The average annual incidence rate was 4.6 per 1000 person-years (95% CI, 4.2–4.9 per 1000 person-years).
Correlates of spatial T wave axis deviation
The distribution of the spatial T wave axis deviation at baseline was similar for men and women, with a mean (SD) of 23.9 (10.7) degrees in men and 23.2 (11.6) degrees in women. Several traditional cardiovascular risk factors showed graded associations across the spatial T wave axis deviation quartiles, with similar patterns among both men and women (Table 1). Participants of African American ethnicity, patients with hypertension and diabetes, and those with higher Cornell voltage presented larger values of spatial T wave axis deviation; lower values for the spatial T wave axis were observed among participants with larger values of the QRS axis, T wave axis in the frontal plane and among smokers. We considered the frontal plane T wave axis for ensuring comparability with other published studies. Other covariates analyzed were anthropometric measures (weight, BMI, waist-to-hip ratio), systolic and diastolic blood pressure, glucose levels, total-, and LDL-cholesterol, QT interval, the prevalence of minor T wave abnormalities. All showed statistically significant, positive associations with spatial T wave axis quartile. Age was not associated with spatial T wave axis deviation.
Of particular relevance is the association of the spatial T wave axis with hypertension and with the heart rate. The relationship between the spatial T wave axis and hypertension status (according to the JNC-VII classification) is presented in Table 2 by gender and adjusted for age, height and weight. The more pronounced the hypertension status the larger the mean values for the spatial T wave axis. Among the 3860 individuals defined as hypertensive, the mean (SD) values for the spatial T wave axis adjusted for age, height and weight were higher among those with uncontrolled HTN compared to those having the blood pressure levels below the treatment goal: 26.73 (0.27) degrees vs. 25.90 (0.27) (p = 0.03).
We found a small but significant inverse association between spatial T wave axis deviation and heart rate. From the first to the fourth T wave axis quartile the mean (SE) values for the heart rate were: 70.24 (0.30), 69.05 (0.23), 69.20 (0.23), 67.70 (0.23) beats per minute for women and 67.26 (0.28), 65.74 (0.28), 65.56 (0.28), 66.67 (0.28) beats per minute for men.
Differences between CHD cases and non-cases
Compared to those who remained free of disease, men and women who subsequently developed CHD tended to have higher mean age, blood pressure, glucose and atherogenic lipids at baseline. They also were more likely to be smokers, have hypertension, diabetes, a higher Cornell voltage, minor T wave abnormalities and a lower mean values for the QRS axis in frontal plane. The mean QT interval and the mean T wave axis in frontal plane were not significantly different between cases and non-cases. In contrast to men, women who developed CHD had higher mean values for body mass index, waist to hip ratio, and heart rate at baseline. Spatial T wave axis deviation was only slightly increased among cases compared to non-cases (25.2 degrees versus 23.2 degrees among women, and 24.4 degrees versus 23.9 degrees among men). Mean values for the QT interval, heart rate, or T wave axis in frontal plane were not statistically different.
Incident CHD events and relative risks of CHD
Overall, the cumulative incidence of CHD was almost uniform across the distribution of the spatial T wave axis deviation, with an increase in the last quartile. The incidence rate was 4.26, 4.18, 4.28 and 5.62 per 1000 person-years across quartiles of spatial T wave axis deviation respectively.
As Table 3 illustrates, the incidence rate ratio was 1.62 (95%CI of 1.15 to 2.27) for women in the upper quartile of the spatial T wave axis deviation compared with women in the bottom quartile. Adjustment for age, height and weight attenuated this ratio to 1.44 (95%CI of 1.02 to 2.04). Among men, the incidence rate ratio comparing the highest with the lower quartile was 1.09 (0.83–1.42), which was unchanged (1.08, 0.83–1.42) after adjustment for age, height and weight.
Cox regression-modeling of the unadjusted association between CHD and spatial T wave axis deviation among women showed that for every 10 degrees increase in the spatial T wave axis deviation, there was a 1.16 (1.04 and 1.28)-fold increase in the risk of CHD (Table 4). When treating spatial T wave axis deviation categorically, borderline and abnormal deviation was associated with a 1.59 (1.20–2.10) and 1.69 (1.04–2.75)-fold higher risk of developing CHD compared to the normal category. When adjusting for age, height, weight, smoking, hypertension, diabetes mellitus, QRS axis and T wave minor abnormalities, this association for women was no longer statistically significant (hazard rate ratio of 1.03 (0.92–1.14). Neither heart rate or the QT interval (regardless of the formula used for heart rate correction), did qualify as confounders, as they were not associated with incident events in this study population.
Among men, there was no significant association between spatial T wave axis deviation and CHD events. This lack of association persisted after adjusting for covariates. T wave axis in the frontal plane was not associated with coronary events either; the crude associations were nil and persisted as such in multivariate analysis in both women and men.
When using the same T wave axis deviation categories as in the CHS and Rotterdam studies, borderline and abnormal deviation was associated with a 1.11 (0.89–1.40) and 1.27 (0.79–2.04)-fold higher risk for developing CHD compared to the normal category. After adjusting for age, height, and weight and as well for smoking, hypertension, diabetes, QRS axis and T wave minor abnormalities, the association remained statistically non-significant. The prevalence of an abnormal spatial T wave axis deviation (larger than 45 degrees) in the ARIC population was extremely low (4.5%) so these risk estimates need to be interpreted with caution.
Minor T wave axis abnormalities were only slightly associated with the risk of developing a coronary heart disease event, in both men and women. The hazard rate ratios for CHD in relation to minor T wave abnormalities were 1.96 (1.37–2.81) for women and 1.31 (0.92–1.88) for men. After adjustment for age, height, weight, smoking, hypertension, diabetes and QRS axis, these estimates became 1.25 (0.86–1.80) for women and 1.02 (0.71–1.47) for men.
The measured (uncorrected for heart rate) QT interval was not associated with coronary events, in either gender group. Incidence rates for coronary events in the first and the last gender-specific QT interval quartiles were 8.04 and 8.00 per 1000 person-years among men and 3.8 and 4.5 per 1000 person-years among women. Adjustment for age and heart rate (either by using the heart rate corrected QT interval by Bazett's formula, QT Index or by including heart rate in the Cox proportional models) did not change this lack of association.
This prospective, population-based, bi-ethnic study of men and women free of coronary heart disease at baseline shows that spatial T wave axis deviation is not associated with incident coronary events during long-term follow-up. The QT interval was not associated, and minor T wave abnormalities were only weakly associated, with the incidence of coronary events after adjustment for various covariates.
Given the lack of association between spatial T wave axis and coronary events in this study, several issues must be considered. Potentially, incomplete ascertainment or misclassification of CHD events could bias the results and obscure an underlying association. The excellent quality control strategies employed by the ARIC study, the highly standardized protocol for data collection, independent case ascertainment through the cohort follow-up, and the epidemiologic surveillance system in place at each of the four ARIC study communities make the possibility of misclassification of coronary heart disease very unlikely.
Inadequate statistical power to detect a meaningful association is not likely, due to the large sample size, long follow-up and relatively large number of events. The minimum effect size that the present study could have missed was computed by dichotomizing spatial T wave axis deviation, then contrasting "borderline" and "abnormal " deviation with the "normal" category (nQuery Advisor, version 4.0 ). Given the observed distribution of T wave axis deviation in the present study, a log-rank test of the survival curves at 12 years of follow up would have approximately 79 to 99% power to detect differences between the survival curves corresponding to hazard ratios between 1.11 and 1.22, at a 0.05 two-sided significance level.
Our results do not confirm the findings of the Rotterdam study  and the CHS  with regard to the spatial T wave axis. The risk of CHD associated with abnormal deviation (more than 45 degrees) versus normal T wave axis deviation (less than 30 degrees) was 2.9 (95% CI 2.0–4.3)) in the Rotterdam study and 1.58 (95%CI 1.25–1.99) in the CHS study. In contrast, findings from the present study are similar to those from the Multiple Risk Factor Intervention Trial (MRFIT) , which did not find a significant association between spatial T wave axis deviation at baseline and incident CHD events, although there was a significant relation between change in spatial T axis and long-term mortality from CHD.
The discrepancy between our results and those reported in previous articles has several potential explanations, involving true variation between study populations (differences between study populations in terms of exposure distribution, demographic characteristics, absolute cardiovascular risk level or effect modification by genetic or environmental factors) and different patterns of confounding by other cardiovascular risk factors.
Compared to participants in the present study, those in studies that found a significant association between T wave axis deviation and CHD were older (Rotterdam study and CHS), and more likely to have a history of myocardial infarction, angina pectoris (Rotterdam study) or silent (ECG-diagnosed) myocardial infarction and hypertension at baseline (CHS study). These study populations also had more ECG evidence of ischemia and infarction as reflected by the presence of Q-waves and major ST-T abnormalities (Rotterdam study) or T wave inversion (CHS study) but this was not accounted for these confounders in the analysis. By contrast, the MRFIT study and the present ARIC study excluded these participants, preserving only those with minor T wave abnormalities, and found only weak or non-significant associations with CHD events. Studies reporting stronger associations also had shorter follow-up time. To explore the potential impact of these dissimilarities, we performed several post hoc analyses under different scenarios.
i) without excluding any participants
In comparing individuals without CHD at baseline with prevalent CHD cases at baseline, the spatial T wave axis deviation was significantly different: 26.66 (16.81) and 44.88 (27.29) degrees respectively (p value < .0001). The entire distribution of the spatial T wave axis deviation among those with CHD at baseline was shifted to the right of the distribution of participants without CHD at baseline. This difference persisted statistically significant even after adjustment for age and several standard cardiovascular risk factors (anthropometry, smoking, hypertension, diabetes, heart rate, QRS axis and QT interval) with adjusted mean values of 26.79 and 41.47 degrees respectively. However, among those without CHD at baseline in our study, there were various conditions or markers of clinical interest, some of established predictive value, such as ST segment depression, negative T waves or atrial fibrillation. In the clinical setting, the need for improved prediction and discriminatory ability of a new potentially useful marker is particularly useful in the absence of such markers. After excluding Minnesota codes representing secondary repolarization abnormalities, this difference was markedly attenuated, with mean spatial T wave axis deviation of 23.84 (12.36) for those free of CHD and 24.08 (15.34) degrees for those with CHD. This suggests that spatial T wave axis has no ability to discriminate between CHD cases and CHD-free individuals beyond and above secondary repolarization abnormalities. This might be one of possible explanations for the discrepant results between different cohort studies exploring T wave axis, as different studies have used different exclusion criteria.
As an alternative to the exclusion of Minnesota codes reflecting secondary repolarization abnormalities, we performed a separate analysis this time including QRS duration as a covariate. The relationship of the spatial T wave axis deviation with incident CHD events followed the same pattern as in the restricted dataset, i.e. a slight increase in the incidence rates in the fourth quartile of the T wave axis, for each gender. The crude hazard rate ratios (95%CI) for incident CHD events for a 10 degree increase in the spatial T wave axis were 1.18 (1.13–1.23) for women and 1.14 (1.09–1.19) for men. After adjustment for the same covariates as in the restricted data set and for the QRS duration, the hazard rate ratios (95%CI) were 1.05 (0.99–1.16) and 1.04 (0.98–1.1).
ii) restricted to fatal events only
Ventricular repolarization abnormalities are closely related to the propensity for arrhythmia. Although none of the published population-based studies on T wave axis deviation attempted to ascertain arrhythmic events, it is plausible that the arrhythmic component of the CHD events was the driving force of the observed associations in some published studies. This suggestion is supported by the fact that in both, Rotterdam and CHS study, the association between T wave axis and coronary events was stronger for fatal than non-fatal events. Similarly, the spatial QRS-T angle was recently reported to be associated with fatal but not with non-fatal cardiac events . To explore this issue we restricted our analysis to fatal coronary events. During an average follow-up time of 11.6 years, 143 fatal events were recorded in our study population, 62 among women and 81 among men, with an overall cumulative incidence of 1.17%. The relationship of fatal events with the spatial T wave axis deviation followed the same pattern as the combined outcome, i.e. a slight increase in the incidence rates for CHD only in the fourth quartile of the spatial T wave axis deviation, for each gender. As expected, in the multivariate models (Table 5), the crude hazard rate ratios for fatal events were slightly larger than for the combined outcome. After adjustment for a limited number of risk factors (due to the small number of events), these modest associations were no longer statistically significant. However, these patterns must be interpreted with caution. While in the Rotterdam and CHS studies, the incidence rate for fatal events was 4.2 (95% CI 3.3–5.1) and 5.8 (95% CI 5.0–6.7) per 1000 person year respectively, this rate was much lower in our study: 1.01 per 1000 person-years (95% CI, 0.84–1.17).
iii) similarly, we repeated the analysis without excluding those with negative T waves corresponding to MC 5.2, restricted to the age group greater than 55 years, restricted to participants with both hypertension and diabetes at baseline, restricted to CHD events occurring after a shorter follow-up (4–6 years). While the incidence rates were larger in absolute value in each T wave axis quartile, in relative terms no association with CHD events was detected.
Experimental models to link the spatial T wave axis deviation with coronary evens are still sparse, but it can be speculated that an abnormal spatial T wave axis may possibly reflect disturbances in the repolarization process caused by subclinical myocardial disease, with or without an increased propensity for arrhythmic events, and thus an increased risk for fatal cardiac events. This is supported by the fact that subclinical disease is likely more prevalent in the older population of the CHS and Rotterdam study (which found an association between T wave axis and incident coronary events) than in the younger populations of the ARIC and MRFIT studies (both reporting a lack of association). Some indicators of subclinical disease such as the left ventricular hypertrophy (LVH) are clearly more prevalent in the Rotterdam than in the ARIC study. The presence of detected or undetected subclinical disease (LVH, patchy myocardial fibrosis) may lead to ventricular electrical instability and a higher proportion of cardiac death related to primary arrhythmic events in a population. Thus, measured or unmeasured, subclinical disease may play a confounding role, and suggests one more explanation for the discrepancy of results between different studies.
The ARIC, CHS and MRFIT study protocols all used the same definition of the spatial T wave axis deviation, while the Rotterdam study used the frontal plane T wave axis. There are no data available to explain the discrepancy of results solely on the basis of different measurement of the T wave axis.
The comparison of spatial T wave axis deviation with its temporal counterpart, the QT interval raises an intriguing question. Incident coronary events were strongly associated with T wave axis in the cohort of older adults of the Rotterdam study, which examined a population of older adults in whom a prolonged QT interval was previously reported to be predictive. In the current work, neither spatial T wave axis deviation nor the QT interval (as measured or heart-rate corrected by different formulae) was associated with incident coronary events.
In our study, a statistical interaction between QT interval and spatial T wave axis deviation was not detected within the range of the QT interval in this population (data not shown). It is important to mention that in the Rotterdam study the mean QT interval (measured with identical methodology as in our study) was considerably larger than in the present study. From electrophysiological point of view, an abnormal ventricular repolarization can alter the QT interval, the T wave axis or both. It is therefore possible that a given magnitude of an altered spatial T wave axis in the presence of a QT interval within "the normal" ranges, may not be sufficient to elicit any causal association with incident coronary events.
Our study found an inverse association between spatial T wave axis deviation and heart rate. Heart rate influences have been reported on the frontal plane T wave angle  and on the spatial ST-T vector . The physiologic explanation of the dependency of an angular measure on heart rate is not clear. It is plausible that this phenomenon is similar to some degree with that observed for another index of altered repolarization, namely T wave alternans, characterized by the change in amplitude, morphology and axis of the spatial T wave and associated with an increased risk of sudden cardiac death. Even if the mechanism of heart rate dependency of the spatial T wave axis remains elusive, it raises questions about the need for "rate correction" of this measurement.
Our study confirms the results of other studies [18, 19] in finding a strong association between spatial T wave axis and hypertension status. The explanation is likely related to the presence of an increased left ventricular wall in the context of the left ventricular hypertrophy. It is also possible that increased spatial Twave axis deviation is related to other processes present in the hypertensive myocardium, such as microfocal areas of fibrous tissue and/or increased alteration of ionic channels. These findings suggest that spatial T wave axis deviation may serve as an auxiliary early marker of repolarization abnormalities in hypertensive individuals. Due to the low prevalence of left ventricular hypertrophy in our study population, we were not able to further explore this issue.
Strengths and limitations
The present study has several strengths. The highly standardized data collection procedure and ECG protocol increase the internal validity of the findings. The ECG recording was performed by trained technicians using a standardized protocol for lead placement. The events ascertainment in the ARIC study is ensured by a highly standardized quality assurance protocol. The inclusion of apparently healthy participants in a population-based study, of African American and white men and women, the large sample size and the geographic distribution of the study population enhance generalizability.
This study also has some limitations. The narrow age range precluded an extensive investigation of the role of age. The study also was limited by the lack of accurate measurements for characterizing the elliptical shape of the thoracic cavity. Constitutional and anthropometric variables such as the chest transverse diameter or chest shape have been associated with a shift of the T wave axis  or changes in average T wave potential amplitudes. A slender body with a limited spatial sampling of the thorax directly overlaying the heart can cause lower amplitudes in women than in men . While our study attempted to account for this fact by adjusting for height and weight, a more accurate measurement of thoracic cavity, shape or body surface area would probably result in a better ascertainment of spatial T axis. Nonetheless, previous studies reporting positive findings did not account for the shape of thoracic cavity either, and it is unlikely that such measurement would have influenced the findings of the present study.
Ventricular replarization is a complex process, which cannot be fully captured by only a limited number of indexes such as QT interval and spatial T wave axis, exploring either the temporal aspects or the general direction of this process. Other indexes or more complex descriptors of the T wave loop have been described to be accurate markers of an altered repolarization . The spatial QRS-T angle has recently regained attention and several studies have attested its usefulness in post infarction  and hypertensive patients . The spatial QRS-T angle was also the strongest predictor of fatal cardiac events, even after controlling for the frontal plane T axis, in the large population-based cohort of men and women aged 55 years or older of the Rotterdam study . Unfortunately, this index is not currently available in our database.
In conclusion, this prospective, population-based, bi-ethnic study of men and women free of coronary heart disease at baseline shows that spatial T wave axis deviation is not associated with incident coronary events during long-term follow up. It is doubtful that T wave axis deviation would add benefit in the prediction of CHD events above and beyond the current traditional risk factors. The QT interval was not associated with incident coronary events, and minor T wave abnormalities were only weakly associated among women.
Zareba W, Maison-Blanche P, Locati E: Noninvasive Electrocardiology in Clinical Practice. 2001, Armonk, NY, Futura Publishing Company
Surawicz B: The pathogenesis and clinical significance of primary T-wave abnormalities. Advances in Electrocardiography. Edited by: Schlant RC. 1972, Grune & Stratton, New York
Rautaharju PM: QT and dispersion of ventricular repolarization: the greatest fallacy in electrocardiography in the 1990s. Circulation. 1999, 99: 2477-8.
Rautaharju PM: Why did QT dispersion die?. Card Electrophysiol Rev. 2002, 6: 295-301. 10.1023/A:1016397529393.
Merri M, Benhorin J, Alberti M, Locati E, Moss AJ: Electrocardiographic quantitation of ventricular repolarization. Circulation. 1989, 80: 1301-8.
Antzelevitch C, Sicouri S, Litovsky SH, Lukas A, Krishnan SC, Di Diego JM, Gintant GA, Liu DW: Heterogeneity within the ventricular wall. Electrophysiology and pharmacology of epicardial, endocardial, and M cells. Circulation Research. 1991, 69: 1427-49.
Franz MR, Bargheer K, Rafflenbeul W, Haverich A, Lichtlen PR: Monophasic action potential mapping in human subjects with normal electrocardiograms: direct evidence for the genesis of the T wave. Circulation. 1987, 75: 379-86.
Lux RL, Urie PM, Burgess MJ, Abildskov JA: Variability of the body surface distributions of QRS, ST-T and QRST deflection areas with varied activation sequence in dogs. Cardiovascular Research. 1980, 14: 607-12.
Malik M, Batchvarov VN: Measurement, interpretation and clinical potential of QT dispersion. Journal of the American College of Cardiology. 2000, 36: 1749-66. 10.1016/S0735-1097(00)00962-1.
Harumi K, Burgess MJ, Abildskov JA: A theoretic model of the T wave. Circulation. 1966, 34: 657-68.
Zabel M, Portnoy S, Franz MR: Electrocardiographic indexes of dispersion of ventricular repolarization: an isolated heart validation study. Journal of the American College of Cardiology. 1995, 25: 746-52. 10.1016/0735-1097(94)00446-W.
Zabel M, Acar B, Klingenheben T, Franz MR, Hohnloser SH, Malik M: Analysis of 12-lead T-wave morphology for risk stratification after myocardial infarction. Circulation. 2000, 102: 1252-7.
Badilini F, Fayn J, Maison-Blanche P, Leenhardt A, Forlini MC, Denjoy I, Coumel P, Rubel P: Quantitative aspects of ventricular repolarization: relationship between three-dimensional T wave loop morphology and scalar dispersion. Ann Noninv Electrocardiol. 1997, 2: 146-157.
Priori SG, Mortara DW, Napolitano C, Diehl L, Paganini V, Cantu F, Cantu G, Schwartz PJ: Evaluation of the spatial aspects of T-wave complexity in the long-QT syndrome. Circulation. 1997, 96: 3006-12.
Malik M, Acar B, Gang Y, et al: QT dispersion does not represent electrocardiographic interlead heterogeneity of ventricular repolarization. J Cardiovasc Electrophysiol. 2000, 11: 835-843.
De Ambroggi L, Aime E, Ceriotti C, Rovida M, Negroni S: Mapping of ventricular repolarization potentials in patients with arrhythmogenic right ventricular dysplasia: principal component analysis of the ST-T waves. Circulation. 1997, 96: 4314-8.
Batchvarov V, Dilaveris P, Farbom P, Ghuran A, Acar B, Hnatkova K, Camm AJ, Malik M: New descriptors of homogeneity of the propagation of ventricular repolarization. Pacing Clin Electrophysiol. 2000, 23: 1968-72.
Dilaveris P, Gialafos E, Poloniecki J, Hnatkova K, Richter D, Andrikopoulos G, Lazaki E, Gialafos J, Malik M: Changes of the T-wave amplitude and angle: an early marker of altered ventricular repolarization in hypertension. Clinical Cardiology. 2000, 23: 600-6.
Dilaveris P, Gialafos E, Pantazis A, Synetos A, Triposkiadis F, Gialafos J: The spatial QRS-T angle as a marker of ventricular repolarisation in hypertension. Journal of Human Hypertension. 2001, 15: 63-70. 10.1038/sj.jhh.1001129.
Nowinski K, Jensen S, Lundahl G, Bergfeldt L: Changes in ventricular repolarization during percutaneous transluminal coronary angioplasty in humans assessed by QT interval, QT dispersion and T vector loop morphology. Journal of Internal Medicine. 2000, 248: 126-36. 10.1046/j.1365-2796.2000.00709.x.
Diercks GF, Hillege HL, van Boven AJ, Kors JA, Janssen WM, Grobbee DE, Crijns HJ, van Gilst WH: Relation between albumin in the urine and electrocardiographic markers of myocardial ischemia in patients without diabetes mellitus. American Journal of Cardiology. 2001, 88: 771-4. 10.1016/S0002-9149(01)01849-5.
Vaidean G, Schroeder EB, Whitsel EA, Prineas RJ, Chambless LE, Perhac JS, Heiss G, Rautaharju PM: Short-term repeatability of electrocardiographic spatial T wave axis and QT interval. Journal of Electrocardiology.
Kors JA, de Bruyne MC, Hoes AW, van Herpen G, Hofman A, van Bemmel JH, Grobbee DE: T axis as an indicator of risk of cardiac events in elderly people. Lancet. 1998, 352: 601-5. 10.1016/S0140-6736(97)10190-8.
Rautaharju PM, Nelson JC, Kronmal RA, Zhang ZM, Robbins J, Gottdiener JS, Furberg CD, Manolio T, Fried L: Usefulness of T-axis deviation as an independent risk indicator for incident cardiac events in older men and women free from coronary heart disease (the Cardiovascular Health Study). American Journal of Cardiology. 2001, 88: 118-23. 10.1016/S0002-9149(01)01604-6.
Prineas RJ, Grandits G, Rautaharju PM, Cohen JD, Zhang ZM, Crow RS, MRFIT Research Group: Long-term prognostic significance of isolated minor electrocardiographic T-wave abnormalities in middle-aged men free of clinical cardiovascular disease (The Multiple Risk Factor Intervention Trial [MRFIT]). American Journal of Cardiology. 2002, 90: 1391-5. 10.1016/S0002-9149(02)02881-3.
Garber CE, Carleton RA, Heller GV: Comparison of "Rose Questionnaire Angina" to exercise thallium scintigraphy: different findings in males and females. Journal of Clinical Epidemiology. 1992, 45 (7): 715-20. 10.1016/0895-4356(92)90048-R.
NHLBI, National Heart, Lung, and Blood Institute: Atherosclerosis Risk in Communities (ARIC) Study. Operations manual, no.5. Electrocardiography. version 1.0. 1987, Chapel Hill: ARIC Coordinating Center, School of Public Health, University of North Carolina
Prineas RJ, Crow RS, Blackburn H: The Minnesota Code Manual of Electrocardiographic Findings: Standards and Procedures for Measurement and Classification. 1982, Littleton, Mass: John Wright-PSG Inc
Edenbrandt L, Pahlm O: Vectorcardiogram synthetized from a 12-lead ECG: superiority of the inverse Dower transform matrix. J Electrocardiol. 1988, 21: 361-367. 10.1016/0022-0736(88)90113-6.
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-72.
Rautaharju PM, MacInnis PJ, Warren JW, et al: Methodology of ECG interpretation in the Dalhousie program: NOVACODE ECG classification procedures for clinical trials and population health surveys. Methods Inf Med. 1990, 29: 362-374.
Bazett HC: An analysis of the time-relations of electrocardiograms. Heart. 1920, 7: 353-370.
White AD, Folsom AR, Chambless LE, Sharret AR, Yang K, Conwill D, Higgins M, Williams OD, Tyroler HA: Community surveillance of coronary heart disease in the Atherosclerosis Risk in Communities (ARIC) Study: methods and initial two years' experience. Journal of Clinical Epidemiology. 1996, 49: 223-33. 10.1016/0895-4356(95)00041-0.
The ARIC investigators: The Atherosclerosis Risk in Communities (ARIC) Study: design and objectives. The ARIC investigators. American Journal of Epidemiology. 1989, 129: 687-702.
Cox D: Regression models and life-tables. J Royal Stat Soc. 1972, 34: 187-220.
Elashoff JD: nQuerry Advisor Version 4.0 User's Guide. 2000, Los Angeles, CA
Kardys I, Kors JA, van der Meer IM, Hofman A, van der Kuip DA, Witteman JC: Spatial QRS-T angle predicts cardiac death in a general population. Eur Heart J. 2003, 24 (14): 1357-64. 10.1016/S0195-668X(03)00203-3.
Lehmann MH, Yang H: Sexual dimorphism in the electrocardiographic dynamics of human ventricular repolarization: characterization in true time domain. Circulation. 2001, 104 (1): 32-8.
Ishikawa K: Correlation coefficients for electrocardiographic and constitutional variables. Am Heart J. 1976, 92: 152-61.
van Oosterom A, Hoekema R, Uijen GJ: Geometrical factors affecting the interindividual variability of the ECG and the VCG. Electrocardiol. 2000, 33 (Suppl): 219-27.
Rubulis A, Jensen J, Lundahl G, Tapanainen J, Wecke L, Bergfeldt L: T vector and loop characteristics in coronary artery disease and during acute ischemia. Heart Rhythm. 2004, 1: 317-325. 10.1016/j.hrthm.2004.03.076.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2261/5/2/prepub
The authors thank the staff and participants in the ARIC study for their important contributions.
The author(s) declare that they have no competing interests.
GV performed the overall study design, statistical analysis and manuscript preparation.
PR designed the data collection protocol, participated in the interpretation of the ECG data, provided comments on the manuscript.
RP participated in the quality control of the ECG data, participated in the interpretation of the findings and with comments on the manuscript.
EAW participated in the interpretation of the findings and with comments on the manuscript, provided assistance in the manuscript preparation.
LEC provided oversight of the statistical analysis, verified the appropriateness of statistical models, provided assistance in the manuscript preparation.
ARF participated in the interpretation of the findings and with comments on the manuscript, provided assistance in the manuscript preparation.
WDR collaborated in the study design, participated with comments on the manuscript.
ZMZ managed the automated coding of ECGs and the quality assurance.
RC participated with comments on the manuscript.
GH collaborated in the design of the study, provided comments on the manuscript.
About this article
Cite this article
Vaidean, G.D., Rautaharju, P.M., Prineas, R.J. et al. The association of spatial T wave axis deviation with incident coronary events. The ARIC cohort. BMC Cardiovasc Disord 5, 2 (2005). https://doi.org/10.1186/1471-2261-5-2
- Coronary Heart Disease
- Coronary Heart Disease Event
- Rotterdam Study
- Cardiovascular Health Study
- Wave Axis