| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Circulation. 2000;102:2152.)
© 2000 American Heart Association, Inc.
Basic Science Reports |
From the Cardiac Bioelectricity Research and Training Center (CBRTC) and the Department of Biomedical Engineering, Case Western Reserve University, Cleveland, Ohio (J.E.B., Y.R.), and the Cardiovascular Research and Training Institute, University of Utah, Salt Lake City (B.T.).
Correspondence to Yoram Rudy, Director, Cardiac Bioelectricity Center, Wickenden 509, Case Western Reserve University, Cleveland, OH 44106-7207. E-mail yxr@po.cwru.edu yxr{at}po.cwru.edu
| Abstract |
|---|
|
|
|---|
Methods and ResultsEpicardial potentials were recorded during VT with a 490 electrode sock during an open chest procedure in 2 dogs with 4-day-old myocardial infarctions. Body surface potentials were generated from these epicardial potentials in a human torso model. Realistic geometry errors and measurement noise were added to the torso data, which were then used to noninvasively reconstruct epicardial isochrones, electrograms, and potentials with excellent accuracy. ECGI reconstructed the reentry pathway and its key components, including (1) the central common pathway, (2) the VT exit site, (3) lines of block, and (4) regions of slow and fast conduction. This allowed for detailed characterization of the reentrant circuit morphology.
ConclusionsECGI can noninvasively image arrhythmic activation on the epicardium during VT to identify and localize key components of the arrhythmogenic pathway that can be effective targets for antiarrhythmic intervention.
Key Words: electrocardiography myocardial infarction reentry tachycardia
| Introduction |
|---|
|
|
|---|
Clearly, we need a true imaging modality that can noninvasively provide detailed electrophysiological information on the surface of the heart. Potential applications include screening patients at risk for sudden cardiac death, noninvasive evaluation of electrophysiological changes caused by drug therapy, and guiding nonpharmacological interventions, such as ablation and pacing, in terms of optimal positioning relative to the arrhythmogenic substrate. Furthermore, advances in genetic and molecular biology have helped determine the genetic basis and cellular mechanisms of various cardiac arrhythmias (ie, long-QT syndrome6 7 8 9 10 11 ). A true imaging modality for cardiac electrical function will be needed for the diagnosis and treatment (eg, targeted drug delivery or genetic/molecular modification) of such hereditary rhythm disorders.
The most common mechanism of clinical cardiac arrhythmias is reentry,12 13 an abnormal propagation of the excitation wavefront that forms closed loops or spiral waves14 that capture the heart and disturb the normal pattern of synchronized excitation necessary for efficient blood pumping. Reentry develops in pathological substrates, such as those associated with ischemia and infarction12 or electrophysiological remodeling,15 16 17 18 and it involves slow conduction and conduction block. In this report, we demonstrate for the first time the feasibility of the noninvasive imaging of reentrant cardiac arrhythmias.
Electrocardiographic imaging (ECGI) is a noninvasive method for reconstructing detailed electrophysiological information on the surface of the heart using body surface measurements.19 Using a recently developed and validated ECGI methodology,20 21 22 we reconstructed the reentry pathway and activation pattern of the heart during ventricular tachycardia (VT). We showed that ECGI can noninvasively reconstruct the sequence of activation on the epicardium, including the double-loop reentry pathway, regions of slow conduction, and epicardial potentials and electrograms during VT. Importantly, lines of conduction block and the common pathway of the 2 reentrant waves (a strategic site for ablation) are accurately reconstructed and located by this approach.
| Methods |
|---|
|
|
|---|
The VT was induced through programmed stimulation with a basic cycle length drive train and 1, 2, or 3 premature extrastimuli. Epicardial potentials were recorded with a temporal resolution of 1 ms using a sock containing 490 electrodes. The potentials were placed on the surface of a digitized heart, which was positioned in the correct anatomical location within the human torso model, as described in a previous study.22 The recorded epicardial potentials were used to generate electrocardiographic potentials on the human torso surface; these potentials were contaminated by measurement noise and electrode-position errors. Fifty microvolt peak-to-peak gaussian noise (0.5% of the signal) was added to the electrocardiographic torso potentials, and 2 mm geometric error with gaussian distribution was added to the torso surface points to represent realistic electrocardiographic recordings. A more detailed description of this model and discussion of its advantages and limitations was published previously.22
Computational Methods
Epicardial potentials were noninvasively reconstructed from the
electrocardiographic torso potentials over the entire surface of the
heart and throughout the cardiac cycle using previously developed and
validated ECGI mathematical methodology.21 25 Briefly,
ECGI noninvasively computes potentials on the heart surface by solving
Laplaces equation within the torso volume, using torso surface
potentials and the geometric relationship between the epicardial and
torso surfaces as inputs. The first step in this procedure is the
discretization of the relationship between the epicardial potentials
and the torso surface potentials. The discretization is obtained
through the use of the boundary element method, which results in the
linear matrix relationship of
VT=AVE, where
VT and VE represent
vectors of torso and epicardial surface potentials, respectively. The
matrix A encodes the geometric relationship between the epicardial and
torso surfaces as discretized through the boundary element method.
The second step is the inversion of the linear matrix relation to
obtain epicardial potentials from the measured torso surface
potentials. Tikhonov regularization26 is used to
stabilize the solution because this problem is ill-posed (ie, unstable
in the presence of noise). The Composite Residual and Smoothing
Operator method27 is used to determine the degree of
regularization required for an accurate estimate of the epicardial
potentials. Two consecutive reentry cycles of monomorphic VT are
presented in detail as representative data.
Reconstructions were performed for 4 additional reentrant cycles in
both dogs. The accuracy of reconstruction was similar for all cases,
and statistical measures are provided in the
Table
. In addition to spatial epicardial
potential maps, epicardial electrograms (potential over time at a
single site) were reconstructed at 490 evenly distributed epicardial
sites and used to construct isochronal maps describing the sequence
of ventricular activation during the
arrhythmia.
|
The general quality of reconstructed potential maps and electrograms
was evaluated with respect to measured data using absolute error,
relative error, and correlation coefficient, as defined
previously22 (definitions are provided in the
Table
). The average absolute error, relative error, and
correlation coefficient presented for potential maps were
averaged over all reconstructed time frames, and the corresponding
values used for electrograms were averaged over all 490 epicardial
points. Activation times were determined as the time of maximum
negative derivative in the epicardial electrograms.21
Information from neighboring electrograms was used to edit activation
times in a small number (<5%) of electrograms with multiple large
negative derivatives. Lines of block were drawn to separate sites with
activation time differences >30 ms.
| Results |
|---|
|
|
|---|
|
The noninvasively reconstructed isochrones for the first cycle show
the same reentrant pattern of activation, with the CCP and lines of
block reconstructed in approximately the same locations as their
directly measured counterparts. The crowding of isochrones around
the tips of the lines of block indicates a slowing of conduction as the
wavefronts turn around the pivot points on exiting from the CCP. The
sparse isochrones across the right ventricle (120 to 165 ms)
reflect faster conduction. These variations of conduction velocity are
captured in the noninvasive isochrone map. The measured
isochrones from the second cycle of reentry (Figure 1
, bottom) show the same reentrant pattern, with lines of block forming in
approximately the same locations as during the first cycle. Again, the
ECGI-reconstructed isochrones capture the morphology of the
reentrant circuit, allowing for a detailed noninvasive identification
of the arrhythmia pathway.
Electrograms were recorded from 12 locations around the heart
surface during 2 cycles of the VT; they are shown in Figure 2
with their noninvasively reconstructed
counterparts. Sites A through C are in the region of the reentrant
circuit CCP. The noninvasive electrograms are similar to the measured
electrograms, with low-level reconstruction noise and suppression of
peak-to-peak amplitude at site A. Sites D through F are located near
the exit region of the CCP. The timing of the maximum negative
derivatives of the electrograms (time of local activation) shows that
activation spreads from site D to sites E and F at about the same time,
indicating the formation of the superior and inferior arms,
as shown in Figure 1
. The reconstructed electrograms at these
sites closely match the measured electrograms, resulting in the same
estimated sequence of activation.
|
Electrograms from sites G-I lie over the infarct-free anterior right ventricle and have sharp RS (positive-negative) deflections, which are characteristic of the passage of a well-organized activation wavefront through normal myocardium. The reconstructed electrograms at sites G and I closely resemble the measured electrograms. At site H, the reconstructed electrogram has a reduced amplitude. However, the timing of the sharp negative deflection at this site still matches the timing of the measured deflection, providing an accurate activation time estimation. Electrograms at sites J through L, from the posterior surface of the heart, show very sharp negative deflections, indicating the passage of a well-organized wavefront under the recording site. ECGI noninvasively reconstructs the morphologies of these electrograms with good accuracy.
Note that, on average, the maximum negative derivatives of the
reconstructed electrograms are 30% smaller in magnitude (less sharp)
relative to the measured ones. However, the timing of their occurrence
is similar in both, resulting in accurate determination of local
activation times and isochrones (Figure 1
). The Table
summarizes the average errors of the 490 reconstructed electrograms for
the 6 reentrant cycles.
Potential distributions on the epicardial surface at various time
points throughout the VT provide yet another method for mapping the
reentrant circuit. Epicardial potential maps are shown in Figure 3
, together with the corresponding torso
surface electrocardiographic potentials. The location of the activation
wavefront on the epicardial surface can be estimated from the presence
of sharp potential changes (sharp spatial gradients).28 In
some cases, the wavefront also marks the boundary between negative and
positive potentials (propagation is toward positive potentials). At 5
ms (Figure 3A
), the wavefront is entering the CCP while the
previous cycle terminates over the basal right ventricle. This is
reflected in the positive regions, both at the opening of the CCP and
over the basal right ventricle, which are separated by a large negative
region over the anterior right ventricle. This potential pattern is
noninvasively reconstructed from the torso potential map using
ECGI.
|
At 50 ms (Figure 3B
), the wavefront has exited the CCP and
propagated inferiorly, as indicated by the large magnitude
positive region and steep gradients over the inferior left
ventricle. This progression of the wavefront is accurately
reconstructed in the noninvasive epicardial map. At 105 ms (Figure 3C
), the reentrant wavefront has formed 2 arms that envelope the
CCP. ECGI reconstructs this complex potential pattern, correctly
showing the 2 reentry arms of the circuit. Note that the torso
potentials at this time only show a single anterior negative region,
which does not reflect the complexity of the underlying epicardial
pattern. Despite this lack of detail, ECGI can reconstruct the
complexity of the epicardial potentials from the smooth and low
resolution torso pattern. Potential map reconstruction errors for each
reconstructed reentrant cycle are presented in the
Table
.
| Discussion |
|---|
|
|
|---|
The infarct substrates produced in this study through LAD
occlusion23 24 differed from the ethanol-induced infarct
used previously to evaluate ECGI in the context of infarct
substrates.22 The abnormal substrate produced by LAD
occlusion resulted in reentrant arrhythmias with epicardial
participation in the area over the infarct (epicardial border zone) in
both dogs studied. For all 6 reconstructed reentrant cycles, an
epicardial reentrant circuit similar to those presented in
Figures 1 through 3![]()
![]()
was recorded. The
thickness of the spared tissue in the epicardial border zones varied
from 0 to 2 mm. ECGI reconstructions accurately reproduced the
general circuit morphology, lines of conduction block, and regions of
slow and fast conduction in all cases. Although clinical reentry often
occurs in the 3D myocardium, the subepicardium plays an
important role in the maintenance of reentry in up to 33% of
patients undergoing epicardial photocoagulation or ablation
therapy.29 30 31 32 Failure to locate and ablate subepicardial
target sites has been suggested as a reason for the failure of catheter
ablation to treat patients with VT.33 34 Thus, the
noninvasive reconstruction of epicardial
electrophysiological information
demonstrated in this study would provide valuable clinical information
that is presently not available.
It should be noted that ECGI need not be limited to epicardial reentry, because previous studies have shown that information on intramural activity can be obtained from the epicardial potential patterns and their evolution in time.21 35 36 Furthermore, in the case of transmural infarction, ECGI could be used to characterize the electrophysiology of the infarct substrate22 and identify sites of epicardial breakthrough during reentrant activation. Interpretation of intramural arrhythmogenic activity could be further enhanced by direct catheter mapping37 38 or noncontact catheter reconstruction39 40 41 of electrophysiological information on the endocardial surface simultaneously with noninvasive epicardial ECGI. The combination of epicardial and endocardial electrophysiological information, with knowledge of the intramural anatomic organization of the myocardium, would provide an unprecedented ability to localize arrhythmogenic activity within the myocardial depth using only noninvasive or minimally invasive procedures.
Noninvasive localization of critical regions such as the CCP is a prerequisite for guiding noninvasive therapeutic intervention. Once a target is identified and localized, focused energy could be applied for noninvasive ablation.42 ECGI could also be used to immediately and noninvasively view reentrant circuit changes caused by drug administration. Finally, ECGIs ability to detect abnormal electrophysiological substrates22 can provide a noninvasive procedure for identifying patients at a high risk of life-threatening arrhythmias. After screening, prophylactic measures (eg, implantable defibrillators, ablation, drug therapy, or genetic/molecular modification) can be taken before sudden cardiac death occurs.
This study builds on the existing foundation of work20 21 22 43 toward the clinical implementation of ECGI through careful validation using controlled and verifiable combinations of experimental and modeling approaches. In this and previous studies,20 21 22 the exact geometry of the epicardial and torso surfaces was known, and the volume between the 2 surfaces was assumed homogeneous (eg, no lungs). A previous detailed geometric study43 suggested that successful clinical implementation of ECGI could be achieved through the fusion of existing anatomic imaging modalities (eg, CT or MRI), with ECGI for patient-specific acquisition of the torso and epicardial geometry. The assumption of a homogenous torso used in this study is valid because it has been previously demonstrated to have minimal effect on epicardial potential patterns, electrogram waveforms, or isochrones.28 44 Of course, implementing ECGI in the clinical setting will require additional validation in patients with known activation sequences; this can be done through direct comparison with epicardial maps recorded during open heart surgery.
| Acknowledgments |
|---|
| Footnotes |
|---|
Received April 10, 2000; revision received May 11, 2000; accepted June 8, 2000.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. Berger, G. Fischer, B. Pfeifer, R. Modre, F. Hanser, T. Trieb, F. X. Roithinger, M. Stuehlinger, O. Pachinger, B. Tilg, et al. Single-Beat Noninvasive Imaging of Cardiac Electrophysiology of Ventricular Pre-Excitation J. Am. Coll. Cardiol., November 21, 2006; 48(10): 2045 - 2052. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Ramanathan, P. Jia, R. Ghanem, K. Ryu, and Y. Rudy Activation and repolarization of the normal human heart under complete physiological conditions PNAS, April 18, 2006; 103(16): 6309 - 6314. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Zhang, I. Ramachandra, Z. Liu, B. Muneer, S. M. Pogwizd, and B. He Noninvasive three-dimensional electrocardiographic imaging of ventricular activation sequence Am J Physiol Heart Circ Physiol, December 1, 2005; 289(6): H2724 - H2732. [Abstract] [Full Text] [PDF] |
||||
![]() |
Members of the Sicilian Gambit New Approaches to Antiarrhythmic Therapy, Part I: Emerging Therapeutic Applications of the Cell Biology of Cardiac Arrhythmias Circulation, December 4, 2001; 104(23): 2865 - 2873. [Abstract] [Full Text] [PDF] |
||||
![]() |
Members of the Sicilian Gambit New approaches to antiarrhythmic therapy; emerging therapeutic applications of the cell biology of cardiac arrhythmias Eur. Heart J., December 1, 2001; 22(23): 2148 - 2163. [Abstract] [PDF] |
||||
![]() |
J. E. Burnes, B. Taccardi, P. R. Ershler, and Y. Rudy Noninvasive electrocardiogram imaging of substrate and intramural ventricular tachycardia in infarcted hearts J. Am. Coll. Cardiol., December 1, 2001; 38(7): 2071 - 2078. [Abstract] [Full Text] [PDF] |
||||
![]() |
Members of the Sicilian Gambit New approaches to antiarrhythmic therapy: emerging therapeutic applications of the cell biology of cardiac arrhythmias Cardiovasc Res, December 1, 2001; 52(3): 345 - 360. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. E. Burnes, R. N. Ghanem, A. L. Waldo, and Y. Rudy Imaging Dispersion of Myocardial Repolarization, I: Comparison of Body-Surface and Epicardial Measures Circulation, September 11, 2001; 104(11): 1299 - 1305. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. N. Ghanem, J. E. Burnes, A. L. Waldo, and Y. Rudy Imaging Dispersion of Myocardial Repolarization, II: Noninvasive Reconstruction of Epicardial Measures Circulation, September 11, 2001; 104(11): 1306 - 1312. [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. |