Unstable QT Interval Dynamics Precedes Ventricular Tachycardia Onset in Patients With Acute Myocardial InfarctionClinical Perspective
A Novel Approach to Detect Instability in QT Interval Dynamics From Clinical ECG
Background—Instability in ventricular repolarization in the presence of premature activations (PA) plays an important role in arrhythmogenesis. However, such instability cannot be detected clinically. This study developed a methodology for detecting QT interval (QTI) dynamics instability from the ECG and explored the contribution of PA and QTI instability to ventricular tachycardia (VT) onset.
Methods and Results—To examine the contribution of PAs and QTI instability to VT onset, ECGs of 24 patients with acute myocardial infarction, 12 of whom had sustained VT (VT) and 12 nonsustained VT (NSVT), were used. From each patient ECG, 2 10-minute-long ECG recordings were extracted, 1 right before VT onset (onset epoch) and 1 at least 1 hour before it (control epoch). To ascertain how PA affects QTI dynamics stability, pseudo-ECGs were calculated from an MRI-based human ventricular model. Clinical and pseudo-ECGs were subdivided into 1-minute recordings (minECGs). QTI dynamics stability of each minECG was assessed with a novel approach. Frequency of PAs (fPA) and the number of minECGs with unstable QTI dynamics (Nus) were determined for each patient. In the VT group, fPA and Nus of the onset epoch were larger than in control. Positive regression relationships between fPA and Nus were identified in both groups. The simulations showed that both fPA and the PA degree of prematurity contribute to QTI dynamics instability.
Conclusions—Increased PA frequency and QTI dynamics instability precede VT onset in patients with acute myocardial infarction, as determined by novel methodology for detecting instability in QTI dynamics from clinical ECGs.
- ventricular tachycardia
- QT interval
- acute myocardial infarction
- premature activation
- repolarization instability
Instability in the dynamics of ventricular repolarization plays an important role in the mechanisms of arrhythmia, especially when a premature activation (PA) is presented. At the cellular level, ventricular repolarization is measured by the action potential duration (APD), the latter often expressed as a function of the preceding diastolic interval (DI).1 This relationship is known as APD restitution. Over the past decade, much emphasis has been placed on the maximum slope of the restitution curve as a major factor in both the onset of arrhythmias and the dynamic destabilization of reentrant waves that underlie arrhythmias. A large (>1) maximum APD restitution slope indicates instability in APD dynamics1–3; it is observed in the short DI range, which is typically associated with PAs. The shorter the DI, the larger the maximum restitution slope. Research has demonstrated that unstable APD dynamics results in spatial gradients in APD,4,5 leads to ventricular tachycardia (VT) after PA,6–8 and causes the transition from VT to ventricular fibrillation.2,3,5
Clinical Perspective on p 866
To test the value of the restitution slope as an arrhythmia risk indicator in the clinic, several studies have measured this slope in patients using MAP catheters.9,10 These studies reported increased APD restitution slopes in diseased human hearts, as identified by invasive constant-rate pacing protocols. The use of such pacing protocols eliminated short-term memory, which is the dependence of APD on activation history. Activation history is typically represented by several preceding APDs and DIs.11,12 Because short-term memory can enhance or diminish instability in APD dynamics,13,14 APD restitution slope cannot be used as an accurate measure of APD instability in the clinical setting.15
The QT interval (QTI) of an ECG is the global manifestation of ventricular APD. Based on the concept of APD restitution, larger-than-1 slope of QTI restitution (dependence of QTI on the preceding TQ interval, TQI), which indicates instability in QTI dynamics, has been tested as an index of arrhythmia risk.16,17 However, similar to APD restitution, QTI restitution slope has to be measured under invasive constant-rate pacing protocols, and similarly eliminates the effect of short-term memory on QTI dynamics. Currently, there is no robust methodology that can detect instability in QTI dynamics noninvasively from the clinical ECG and without canceling the important contribution of short-term memory. Furthermore, it remains unknown whether instability in QTI dynamics, if detected, precedes the onset of arrhythmia. Finally, the contributions of PA to instability in QTI dynamics and the onset of arrhythmias are also unknown.
In the present study, we propose a novel methodology for detecting instability in QTI dynamics. The algorithm defines QTI dynamics as a function of both its preceding R-R interval (RRI) and activation history and assesses the stability of this function. Because the effect of short-term memory is included in the algorithm, the presented approach has an advantage over restitution slope in measuring repolarization instability. The methodology uses the clinical ECG and does not require a specific pacing protocol; it is thus noninvasive. In this study, we test the hypothesis that instability in QTI dynamics, as detected by this methodology, precedes the onset of sustained VT, but not that of nonsustained VT (NSVT). Furthermore, we test the additional hypothesis that instability in QTI dynamics depends on both the PA degree of prematurity (DOP) and the frequency of PAs (fPA). The results of this study reveal the important contributions of PA and QTI instability to VT onset. Confirming these hypotheses would indicate that the present methodology for detecting instability in QTI dynamics could be used to develop a novel index that predicts the risk of arrhythmias.
We tested the first hypothesis using clinical ECG recordings collected from patients with acute myocardial infarction (AMI). However, an ECG recording could have different fPA, and each PA could be of different DOP, which makes directly testing the second hypothesis difficult. Instead, we used pseudo-ECG recordings generated by a novel biophysically detailed, MRI-based electrophysiological model of the human ventricles, so that fPA and the DOP of PAs could be controlled.
The details of the clinical data collection, the computer simulations with the MRI-based human ventricular model, and the algorithm for detecting instability in QTI dynamics are presented in the following sections.
The ECG recordings used in this study are part of a clinical ECG database at the Johns Hopkins Hospital.18 We used 24 AMI patients from this database; 12 patients had sustained and 12 nonsustained monomorphic VT events (VT and NSVT groups, respectively). The clinical demographics of both groups, such as age, sex, diagnosis, β-blocker usage, and antiarrhythmic drug therapy are presented in Table 1.
From the multilead recordings of each patient, the recording with the best signal-to-noise ratio was chosen for analysis. VT events were identified from these ECG recordings by the cardiologist. Two 10-minute-long ECG recordings, extracted immediately before a chosen VT onset and at least 1 hour before that onset, were assembled into an onset epoch and a control epoch, respectively, with 1 recording per patient in each epoch. Each control ECG had to be at least 1 hour after any prior cardiac events. Each 10-minute ECG trace was then divided into 10 1-minute ECG recordings (minECGs).
For each minECG, Q, R, and T waves were annotated and QTIs, TQIs, and RRIs were extracted (Figure 1). The QRS complex was first identified, followed by the T-wave and the isoelectric line, using the approach of Laguna et al19 The end of the T-wave was found as the intersection of the isoelectric line with the tangent to the T-wave at the point of maximum slope. Annotations were visually checked to ensure accuracy. The number of PAs in each minECG was counted. In this study, PA could be a premature atrial contraction (PAC) or a premature ventricular contraction (PVC). Following the approach of Huikuri et al,20 a PA was counted each time when RRI of a beat was shortened by at least 100 ms with respect to that of the preceding beat.
Computer Simulations With the MRI-Based, Biophysically Detailed Model of Human Ventricular Electrophysiology
We used a novel, anatomically realistic model of electrophysiology in the human ventricles (Figure 2A), reconstructed in our laboratory from MRI and diffusion-tensor (DT) MRI data.21,22 The model was used to examine how fPA and the DOP of a PA affected stability of QTI dynamics. In this study, the DOP of a beat was defined as the shortening of the cycle length (CL) with respect to the one of the preceding beat.
Description of our pipeline for heart model construction was published previously.22,23 Mathematical description of the electrophysiological behavior of cardiac tissue was based on the monodomain representation, ensuring speed of execution of the simulations. The membrane dynamics of the human ventricular cell was represented by the ten Tusscher et al ionic model.24 To obtain an ECG consistent with the clinical signal, we included transmural heterogeneities in cellular properties24 in the model ventricles. The spatial distribution of endocardial, epicardial, and M-cells across the ventricular walls (Figure 2B) was implemented as in a previous study,25 based on data by Drouin et al.26
A 1-minute-long pacing train was delivered epicardially at the apex of the MRI-based human ventricular model (Figure 2A). The pacing intervals in this train were the same as the RRIs of one of the PA-free minECGs in the VT group onset epoch. To construct pacing trains with different fPA and DOP, the timings of 4, 8, or 12 randomly selected beats of this PA-free pacing train were shortened to simulate PAs. All PAs in a pacing train had the same DOP, either 150 ms or 200 ms (Figure 2C). The number of beats between the PAs in the pacing train was chosen such that it represented the most frequent pattern of distribution of PAs and sinus beats among the minECGs in the onset epoch of the VT group.
After each 1-minute pacing protocol, the MRI-based human ventricular model was used to generate 1-minute-long pseudo-ECGs (pseudominECGs), as done previously.25 We refer to the simulated ECGs as pseudo-ECGs since the model ventricles were not surrounded by torso25 (Figure 2A), as well as to distinguish them from the clinical minECGs. The model ECG leads were configured to simulate lead II of the standard 12-lead ECG (Figure 2A). The pseudominECGs were annotated in the same way as the minECGs, except that QTI was defined as the interval between Rbegin and Tend (Figure 2D), because pseudo-ECGs do not have a Q-wave.25
Detection of Instability in QTI Dynamics
For each minECG and pseudominECG, the dependence of a QTI on the preceding QTIs and RRIs was expressed as Equation 1 below, which constitutes a linear autoregressive model with exogenous input (ARX). This model has been used previously to describe APD dynamics.12 Equation 1 In Eq1, n is the beat number in the minECG or pseudominECG; ai and bi (i=1,.,M) are the weights (constants) with which each preceding QTI and RRI, respectively, contribute to QTIn. M is the number of beats representing the extent of the activation history. The autoregressive term is (dependence of a given QTI on preceding QTIs), whereas the term represents the exogenous input. In Equation 1, we chose to use RRIs and not TQIs (as done in restitution studies16,17), because TQI is affected by the preceding QTI, and thus it is not an independent exogenous input.
The parameters ai and bi (i=1,.,M) of each ARX model were evaluated with a Matlab function. For each ARX model, the value of M was determined by increasing it from 1, in steps of 1, and examining, at each step, whether an accurate prediction of QTI dynamics in the minECG was achieved. The value of M at which the prediction reached a desired accuracy was denoted as Mmax. The chosen accuracy was that the mean square error between the predicted and the measured QTI was smaller than 5 ms2. We also used residual analysis (95% confidence interval) to validate the ARX modeling at each M value.
We used a well-established method27 to assess the stability of each ARX (ie, the stability of Equation 1) and expressed it as a stability index, Pm. The definition and calculation of Pm can be found in the online-only Data Supplement. Based on the value of Pm, the ARX model was found stable (if Pm <1) or unstable (if Pm ≥1).27 When the ARX model was unstable, the degree of instability was manifested in the value of Pm; the larger the value of Pm, the larger the degree of instability in QTI dynamics.
To compare the performance of our methodology in predicting VT onset with that of other measures of QTI dynamics, we also calculated the QTI restitution slope for each minECG and the QT variability index (QTVI)28 over the entire ECG recording (10 minutes) for each patient.
In the present report, normally distributed data are presented as mean±SD and nonnormally distributed data are presented as median and interquartile range. For comparisons between 2 sets of data, we used the t test for normal continuous data, Wilcoxon rank-sum test for nonnormal continuous data, and Fisher exact test for categorical data. The significance level of all these tests was 0.05. For comparisons between more than 2 sets of data, we used 1-way ANOVA for normal data and Kruskal-Wallis test for nonnormally distributed data. The significance level was 0.05. If significant difference was identified in the preliminary test, post hoc multiple comparisons test (Tukey test) was used to perform pairwise comparisons between the datasets.
Using the above algorithm, each minECG (or pseudominECG) was tagged as either stable or unstable based on its Pm value. To test the first hypothesis, the number of unstable minECGs (Nus) was counted for each ECG recording and was compared between the 4 epochs, that is, the onset and control epochs of the VT and NSVT groups. Other variables, such as fPA, Mmax, and QTVI, were compared in the same way. The fPA in each ECG recording was obtained as the number of PAs divided by the duration of the recording (10 minutes). The relationship between Nus and fPA within each epoch was determined with linear regression. We also examined whether QTI instability evolved before VT onset. For each minute within the onset epochs of the VT and NSVT groups, the mean value of Pm was calculated and analyzed with linear regression and generalized estimating equation model to examine the evolution of Nus in time before VT onset.
The second hypothesis was tested by comparing the stability of pseudominECGs with different fPA and different DOPs.
Using the clinical recordings from the AMI patients, an ARX model (Equation 1) was constructed for each minECG. The value of each QTI in the minECG was predicted with the ARX model and the accuracy of the prediction assessed. As illustrated in Figure 3A, an accurate prediction of QTI dynamics of an example minECG was achieved for Mmax=24. Figure 3B shows the discrepancy between the predicted and the measured QTI dynamics when M=3. The dependence of the mean square error in the prediction on the value of M is presented in Figure 3C, demonstrating that as M increases toward Mmax, the accuracy of the prediction increases.
The value of Mmax was determined from each minECG and is presented in Table 2. One-way ANOVA identified significant differences (P<0.01) in Mmax between the 4 epochs. Post hoc multiple comparisons test further revealed that within both VT and NSVT groups, Mmax was significantly larger during the onset epoch as compared with the control epoch. These findings indicate that in both groups, before VT onset, QTI dynamics has a longer short-term memory than in control.
Analysis of the Stability of QTI Dynamics
Next, the stability in QTI dynamics of each minECG (ie, each ARX) was assessed. We found that instability in QTI dynamics, when present in a minECG, can be observed at different M values, the minimum of which was denoted as Mmin. For the minECG used in Figure 3, Mmin was 3. Although at that value of M the prediction of QTI dynamics was not accurate (Figure 3B), an incomplete activation history (ie, M<Mmax) was nonetheless sufficient to detect instability present in the QTI dynamics of this minECG. The values of Mmin of the 4 epochs are presented in Table 2, and were compared using 1-way ANOVA. The results did not show any significant difference (P=0.07).
The number of unstable minECGs in each ECG recording, Nus, was counted for all 4 epochs. The medians and interquartile ranges of Nus are presented in Table 2. The data points and box plots of Nus of the VT group are presented in Figure 4A and 4B. The Kruskal-Wallis test showed that the 4 epochs are significantly different from each other (P=0.03) in the value of Nus. Multiple comparisons test (Table 3) further identified, within the VT group, a significantly larger value of Nus in the onset epoch compared with the control epoch. This difference was not observed in the NSVT group. This finding indicates that an increase in Nus took place only before the onset of sustained VT but not NSVT.
Figure 5 presents plots of the mean value of Pm at each minute of the onset epochs of the VT and NSVT groups. Linear regression analysis of these data did not identify any trend in QTI instability; R2 was 0.38 (VT) and 0.17 (NSVT). The generalized estimating equation model of Pm also showed that Pm is not dependent on time (P=0.39) or group (P=0.27). These findings indicate that QTI instability did not evolve within the 10 minutes before VT onset.
fPA and Its Relation to QTI Instability and VT Onset
The medians and interquartile ranges of fPA of the 4 epochs are presented in Table 2. The data points and box plots of fPA of the VT group are presented in Figure 4C and 4D. One-way ANOVA identified significant differences (P<0.01) in the values of fPA between the 4 epochs. Post hoc multiple comparisons showed that within the VT group, fPA of the onset epoch was significantly larger than that of the control (Table 4). This difference was not observed in the NSVT group (Table 4). These findings indicate that fPA increased before VT onset in the VT group. Among all minECGs that presented PAs, 78% had PVCs and 35% had PACs. However, in 37.5% of the minECGs that had PACs, PVCs were also present. Similarly, in 16.8% of the minECGs that had PVCs, PACs were also present. Thus, the contributions of PACs and PVCs to QTI instability could not be evaluated separately.
Finally, we found that QTI dynamics in all PA-free minECGs was stable. In all 4 epochs, a positive regression relationship was identified between Nus and fPA (Figure 6), suggesting a cause-and-effect relationship between the 2 variables.
QTI Restitution, QTVI, and CL
For most of the minECGs QTI restitution curves could not be constructed due to the large scatter in QTI (TQI) data. An example QTI restitution curve was constructed from a minECG from the onset epoch of the VT group (Figure 7). PAs were present in this minECG, and our methodology assessed its QTI dynamics as unstable. Figure 7 presents a monoexponential function fitted through the data. Although the maximum slope of the fitted restitution curve was >1 (1.04), the small R2 value (0.38) indicates that the quality of curve fit was very low.
QTVI and CL values (presented in Table 2) were compared among the 4 epochs with 1-way ANOVA. The results showed significant differences in the value of CL (P<0.01) but not in QTVI (P=0.22). Post hoc multiple comparisons demonstrated that CL of the onset epoch was significantly larger than that of the control in the VT group but not in the NSVT group.
QTI Instability in PseudominECGs
The simulation results (Figure 8) demonstrate that at a given DOP, increasing fPA increased QTI instability. For a given fPA, increasing DOP also increased QTI instability, indicating that QTI stability depends on both DOP and fPA.
This study presents a novel methodology for detecting instability in QTI dynamics from the clinical ECG recording. The methodology does not require a specific pacing protocol; it is thus noninvasive and does not exclude the contribution of short-term memory to QTI dynamics. Using this methodology, we demonstrated that QTI dynamics becomes unstable before the onset of sustained VT, and found that this instability is related to the DOP and frequency of PAs. The results also show that fPA increased before the onset of VT but not NSVT. It was also found that in all 4 epochs Nus was positively correlated with fPa, which suggests that QTI instability is unmasked by PAs. These findings shed light on the mechanisms of arrhythmogenesis, and indicate a possibility for the development of an index that predicts arrhythmia risk from the clinical ECG.
Contribution of PAs and Instability in QTI Dynamics to VT Onset
Theoretical studies have shown that instability in APD (resulting from steep restitution), manifested after PA occurrence (which brings the operational point of the system to the portion of the restitution curve with slope >1), leads to heterogeneous distribution of DIs in the heart, and consequentially, to the induction of conduction block, reentry, and VT.6–8 However, demonstrating in the clinical setting that repolarization instability precedes VT onset has proven difficult. Heart rhythm before arrhythmia onset is typically variable, rendering not applicable all restitution-based methodologies for assessing repolarization instability. Applying our novel methodology for assessing repolarization instability, which does not eliminate short-term memory, we were able to demonstrate, from clinical ECG recordings of AMI patients, that unstable QTI dynamics precedes VT onset.
Advantage of the Present Methodology Over Other Approaches
With the inclusion of multiple beats comprising activation history in the ARX model, the present algorithm allows for the detection of instability in QTI dynamics directly from the clinical ECG, without the need for invasive pacing. Accounting for short-term memory provides a great advantage of the present methodology over the traditional restitution approach.
Another QTI-based index, QTVI, has been a valuable tool in arrhythmia risk stratification. QTVI characterizes, in a statistical manner, the relationship between QTI and RRI dynamics over the entire ECG recording and provides an overall estimation of QTI variability normalized by the magnitude of RRI variability.28 Elevated QTVI, which has been reported in diseased hearts, indicates repolarization dynamics that is out of proportion to RRI dynamics.28,29 Our approach evaluates the contribution to QTI dynamics not only of RRIs but also of prior QTIs; it is well understood that the latter contributes to arrhythmia initiation.13,14 In addition, PAs (ie, ectopic beats, which could be either PACs or PVCs) were included in the present study since they could uncover arrhythmogenic unstable repolarization dynamics in the heart; the present methodology captures this instability. These differences between QTVI and the present approach help to explain why our approach successfully predicted the onset of VT while QTVI did not.18
Using the present methodology, we were able to advance the knowledge of the relationship between ectopic beats and arrhythmia mechanisms. Research using heart rate turbulence (HRT) analysis has demonstrated that a single ectopic beat modulates heart rate dynamics through autonomic inputs, with the dynamics of RRI after the ectopic beat being depressed in the diseased heart.30 However, HRT algorithms have not been able to unravel the contribution of ectopic beats and the subsequent altered heart rate dynamics to the onset of arrhythmia. The results of our study revealed that the RRI dynamics after PA underlies the unstable QTI dynamics preceding VT onset.
In the present study, we propose a novel methodology for detecting instability in QTI dynamics, which uses the clinical ECG and does not require a specific pacing protocol. The methodology could be applied in both clinical monitoring of events preceding arrhythmias and arrhythmia risk stratification. While frequent PVCs and short coupling intervals are known risk factors for VT in the AMI patients,31 this study provides another measure of arrhythmia risk in this particular clinical setting, namely unstable QT dynamics. Because all patients, the ECG recordings of whom were used in this study, had AMI and VT (either sustained or nonsustained), the positive regression relationship between fPA and NUS reported here may not apply to the healthy heart. If this is proven to be the case, a positive regression relationship between fPA and NUS in the ECG would indicate arrhythmia risk, and could be used as a risk stratification index. Prospective studies need to be carried out to test this hypothesis.
Both PVCs and PACs were included as PAs in this study. The QRS width of PVC is typically different from that of a sinus beat; sinus beats and PVCs could also be associated with different T-wave morphologies. This could affect QTI annotation. Therefore, analysis of QTI dynamics in ECG recordings with large number of PVCs might be less accurate than that of recordings without PVCs. An additional limitation in this study is that short-term memory was represented by a series of preceding QTIs and RRIs. This representation ignores the fact that the same QTI may be associated with different T-wave shapes.
A limitation of the simulations with the MRI-based electrophysiological model of the human ventricles is that epiventricular pacing, rather than sinus rhythm, was used. Inclusion of the Purkinje system would have complicated dramatically the model and thus was not implemented in this study.
Finally, certain clinical demographics, such as left ventricular ejection fraction, duration of AMI, the use of thrombolytics, and baseline ECG features could not be accessed from patient records and thus could not be reported here because the Johns Hopkins database was constructed under an institutional review board protocol that prohibited retention of any patient-specific identifier.
Sources of Funding
This study was supported by National Institutes of Health grant R01HL082729 and a research grant from Medtronic to Dr Trayanova.
Dr Trayanova is a cofounder of CardioSolv LLC. CardioSolv LLC was not involved in this research.
The online-only Data Supplement is available at http://circep.ahajournals.org/lookup/suppl/doi:10.1161/CIRCEP.110.961763/-/DC1.
- Received December 23, 2010.
- Accepted July 11, 2011.
- © 2011 American Heart Association, Inc.
- Nolasco JB,
- Dahlen RW
- Garfinkel A,
- Kim Y,
- Voroshilovsky O,
- Qu Z,
- Kil J,
- Lee M,
- Karagueuzian H,
- Weiss J,
- Chen P
- Qu Z,
- Garfinkel A,
- Chen PS,
- Weiss JN
- Fox JJ,
- Riccio ML,
- Drury P,
- Werthman A,
- Gilmour RF
- Koller ML,
- Maier SKG,
- Gelzer AR,
- Bauer WR,
- Meesmann M,
- Robert F.,
- Gilmour J
- Huang J,
- Zhou X,
- Smith WM,
- Ideker RE
- Fossa AA,
- Wisialowski T,
- Crimin K,
- Wolfgang E,
- Couderc J-P,
- Hinterseer M,
- Kaab S,
- Zareba W,
- Badilini F,
- Sarapa N
- Huikuri HV,
- Valkama JO,
- Airaksinen KE,
- Seppanen T,
- Kessler KM,
- Takkunen JT,
- Myerburg RJ
- ten Tusscher KH,
- Noble D,
- Noble PJ,
- Panfilov AV
- Lathi BP
- Berger RD,
- Kasper EK,
- Baughman KL,
- Marban E,
- Calkins H,
- Tomaselli GF
- Desanctic R,
- Block P,
- Hutter A
Unstable repolarization dynamics is known to affect the onset of arrhythmia. However, it is very difficult to measure the instability in repolarization dynamics in the clinic. In this report, we present a novel methodology that measures the instability in QT interval (QTI) dynamics from the clinical ECG which includes the contributions of restitution and short-term memory effect on repolarization. This new methodology successfully captured the instability in QTI dynamics before the onset of sustained ventricular tachycardia and identified a positive regression relationship between the frequency of premature activation (fPA) and QTI instability. The reported findings suggest a new way to predict ventricular tachycardia onset by monitoring the development of QTI instability. These findings also indicate that an arrhythmia risk stratification index can be constructed based on the measurement of the relationship between QTI instability and fPA. Compared with QTI restitution slope, which has been previously used as an index of QTI instability, the present methodology included the contributions of both restitution and short-term memory effect and thus is a more comprehensive approach in the measurement of QTI instability. This methodology is noninvasive because it does not need constant pacing protocol, and thus is more feasible in the clinic.