|
|
||||||||
Departments of 1Biomedical Engineering and 2Medicine, University of Alabama at Birmingham, Birmingham, Alabama
| ABSTRACT |
|---|
|
|
|---|
reentry; optical mapping
Periodic wave sources in VF are generally thought to arise from reentry, which occurs when a wave circles back and reactivates tissue it has already passed through. For reentry to exist in the absence of an anatomical obstacle (i.e., a hole in the tissue or an inexcitable region), the pivoting end of the wave must contain a phase singularity (PS), sometimes called a wave break, which is a convergence point of tissue at all phases of the action potential (8, 28). Phase, in this context, is an angular variable that specifies a patch of tissue's progression through the action potential. In two-dimensional mapping data, a wave front is an isoline of the phase value that corresponds to the upstroke. In three dimensions, a wave front is an isosurface, and the phase singularity is a one-dimensional filament snaking through the tissue (17, 28). The intersection of this filament with a surface such as the epicardium is the PS observed in two-dimensional mapping data. A PS has a chirality that is determined by the rotation direction of its wave front. A singular filament with a wave rotating about it is called a rotor. We refer to rotors whose central filament intersects the epicardium as epicardial rotors.
One step in substantiating the mother rotor hypothesis in large hearts is to determine if at least one stable rotor is consistently present during VF. Previous studies that mapped portions of the swine right ventricle (22, 23, 27), left ventricle (16, 21, 22, 27), or septum (10) were unable to document such a rotor, even though relatively long-lived rotors were sometimes observed (21). However, all of these studies mapped only a limited region, typically 2040 cm2, admitting the possibility that a mother rotor was present but located in an unmapped region.
The objective of the present study was to map from the entire epicardial surface of a large heart preparation to determine by direct observation if stable epicardial rotors are consistently present during VF. This was accomplished by using a newly developed optical mapping system that merges fluorescence data acquired from four high-speed video cameras into a single continuous data set (12) and newly developed algorithms that identify and track wave fronts and PSs as they move over the epicardium. In addition, we introduce the notion of a compound rotor in which the central PS can change and describe an algorithm for automatically identifying such patterns.
| METHODS |
|---|
|
|
|---|
Panoramic optical mapping. We used a newly developed optical mapping system that records electrical activity from nearly the entire ventricular epicardium (the tip of the apex is not well imaged and is generally obscured by dripping perfusate). Figure 1 shows a schematic of the system. The shape of the epicardium is acquired with the rotating geometry camera. Electrical activity is recorded by using a voltage-sensitive dye and the mapping cameras. The movies from the mapping cameras are merged into one dataset by mapping them onto the epicardial geometry. Details are as previously described (12) with the following exceptions.
|
2) We replaced the original dual-mapping cameras with four high-speed CCD cameras (iXon DV860, Andor Technology, South Windsor, CT). The cameras are positioned around the heart in 90° increments (128 x 64 pixels in each view). Final spatial resolution with this system is
1.6 mm, and temporal sampling is 750 frames/s. The cameras are synchronized by using a computer-controlled PCI timing board (PCI-CTR05, Measurement Computing). Each mapping camera is fitted with a 2.2-mm, F number/1.0, 1/4-in. format video lens (Pentax, Golden, CO) and positioned
10 cm from the heart. The cameras are mounted on a track so they can be backed out of the way of the geometry camera while it circles the heart and then be returned accurately to the mapping position.
3) Excitation light is now provided by 32 blue light-emitting diodes (LEDs, Luxeon V Star, 470 nm, Lumileds Lighting, LLC, San Jose, CA). Each LED dissipates 5 W with luminous flux of 48 lumens. Six LEDs are mounted on the front of each mapping camera, and two are mounted on each of four positionable arms.
Study protocol.
Hearts were immobilized with 2,3-butanedione monoxime (BDM, 20 mM) to eliminate contraction artifacts in the optical signals and stained with a 57 ml bolus of the voltage-sensitive dye di-4-ANEPPS (15 µM). The hearts were periodically restained during the study if signal amplitude decreased. After an
30-min interval for the preparation to stabilize (12), the following procedure was repeated three times: 1) scan the epicardial geometry; 2) optically map 4 s of a paced rhythm; 3) induce VF by applying a 9-V battery to the right ventricle; 4)
20 s after induction, optically map 4 s of VF; and 5) defibrillate with a minimum reliable-strength shock (520 J) delivered through external paddles. One induction resulted in a monomorphic tachycardia rather than VF; this episode was excluded from analysis. We therefore analyzed a total of 17 VF episodes. At the end of each experiment, all cameras were calibrated as previously described (12).
Geometric reconstruction and decimation.
By using images collected by the geometry camera at 5° increments, the epicardial models were constructed with submillimeter resolution (12). However, because analysis of VF wave fronts requires only about 1- to 2-mm spatial resolution (1), we subsequently reduced the spatial resolution of each model to
1.6 mm. In this procedure, 2,300 points are randomly scattered onto the surface of the high-resolution model. The points are then moved about the surface under the influence of an iterative mutual repulsion function to achieve approximately equal spacing (26). A new triangular mesh is then formed from this set of points using the alpha shapes suite of programs (6). Triangles in the new mesh were typically spaced by
1.6 mm (centroid to centroid). This procedure occasionally produced topological defects in the output mesh (e.g., more than 2 triangles meeting along an edge). These defects were repaired with an interactive program. The model was further corrected so that the three vertexes for each triangle in the mesh were listed counterclockwise when viewed from the outside. This ensured that the normal vectors for each triangle pointed away from the center of the model.
Texture mapping of fluorescence data. Background-subtracted fluorescence data from the high-speed cameras were texture mapped onto the decimated epicardial models by using a modification of our previously described procedure (12). In the modified method, the vertexes of each triangle in the mesh are projected onto the image plane of the camera(s) that optimally map the triangle [by using previously described optimization rules (12)]. The fluorescence assigned to the triangle is a weighted average of fluorescence from each pixel that falls completely or partially within the triangle in the image plane. Pixel weights are determined by the percentage of the pixel that falls within the projected triangle. The fluorescence for triangles that are optimally mapped by more than one camera is a weighted average of the triangle data from each camera image, where weights are determined by the distance of the triangle to the edge of the heart (12).
In one typical geometry, the resolution of the mapping camera images was 1.42 (SD 0.05) mm (mean and SD). The triangles in the mesh were spaced by 1.61 (SD 0.16) mm (centroid to centroid). Triangles received contributions from 8.6 (SD 2.6) image plane pixels. Thus texture mapping imposed some spatial smoothing on the fluorescence data. After texture mapping, the signal for each triangle was temporally smoothed using a low-pass Butterworth filter (55-Hz pass band, 65-Hz stop band, 3-dB pass band loss, 10-dB stop band attenuation) and then normalized to the range 0:100. Spatial and temporal smoothing is useful to reduce high-frequency shot noise present in the optical recordings. Without smoothing, spurious signal fluctuations can cause temporary wavebreaks to appear in activation wave fronts and inappropriately increase the number of very short-lived PSs registered by our analysis algorithms. Filtering may obscure fine spatial and temporal details of the activation patterns but should not affect the macroscopic, long-lived patterns that are the primary focus of this study.
All optical recordings were mapped onto the geometric model acquired most recently before the mapping run. The atria and regions with poor signals acquired during pacing were clipped from the model by using an interactive program. Such regions typically included small fat-obscured areas at the atrioventricular groove, and the tip of the apex, which was poorly imaged by all mapping cameras and was also obscured by dripping perfusate.
Wave front and PS analysis. The resulting VF datasets were processed with newly developed algorithms that automatically identify PSs and wave fronts and construct data structures that describe wave front-wave front, PS-wave front, and PS-PS relationships (20). Wave front-wave front interactions are codified in a directed graph (2) in which the graph edges are wave fronts and the vertexes are wave front initiation, termination, or contact events (20, 24). Contact events include fragmentations, in which one parent wave front breaks up into two or more child wave fronts, and collisions, in which two or more parent wave fronts collide and coalesce to produce a single child wave front. We call this data structure a wave front graph. Figure 2 shows a simple example.
|
|
8 mm), a graph edge is added between the two corresponding PSs. The edges are directed from continuing-to-beginning or ending-to-continuing PSs. The singularity graph for the black singularities in Fig. 3, AC, is shown in Fig. 3E. There are two nodes (i and ii) and a single edge directed from i to ii. Although we have described the algorithm with a fragmentation event, it works in the same way for collisions. For example, if the flow of time is reversed in Fig. 3, wave front 3 appears de novo on the epicardium and coalesces with wave front 2 to form wave front 1. In this situation, PS i is again a continuing PS and PS ii is an ending PS. In the singularity graph, the arrow is directed from ii to i.
A singularity graph may contain multiple subgraphs, each of which is disconnected from all other subgraphs. In graph theory, each of these subgraphs is called a component of the full graph (2). In the limiting case of a VF episode in which no events such as the one depicted in Fig. 3, AC, occurred, both singularity graphs (one for each chirality) would consist only of nodes with no connecting edges. Each node would therefore be its own component. The singularity graph in Fig. 3E is the opposite limiting case: all nodes are connected by edges, and the entire graph therefore consists of a single component spanning all PSs. We treat each singularity graph component as a compound rotor. A compound rotor contains one or more PSs, all of which have the same chirality. The number of PSs in the compound rotor is the number of nodes in the singularity graph component. Its lifetime spans the lifetimes of all of its constituent PSs.
| RESULTS |
|---|
|
|
|---|
20 s after induction. Figure 4 shows a snapshot of a panoramically mapped VF episode. Figure 4, A and B, shows fluorescence, and Fig. 4, C and D, shows the wave fronts and PSs found by our algorithms. Supplemental movies S11 and S22 animate the entire episode.
|
|
Figure 6 shows the distribution of compound rotor lifetimes exceeding 200 ms averaged across episodes. Compound rotors lasting longer than 1 s were rare. The lifetime of the longest-lived compound rotor from each episode is listed in Table 1. One episode contained two compound rotors that persisted for the entire mapped interval (episode 11; the listed time is slightly less than 4.0 s because a few temporal samples were removed from either end of the episode as part of the filtering process). Supplemental movie S33 animates this episode. None of the other episodes contained a compound rotor that persisted for the entire mapped interval.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
20 s after induction. This corresponds to Wiggers' stage II VF, which is clinically relevant because it is the time during which implantable defibrillators deliver therapy. Our major finding is that stable epicardial rotors are not consistently present during VF in this experimental model. Compound rotors. In the strictest sense, the sequence shown in Fig. 3 contains two counterclockwise rotors: one that rotates about PS i and terminates with the termination of wave front 3, and another that rotates about PS ii and is formed when wave front 1 fragments into wave fronts 2 and 3. We argue that sequences such as this should be regarded as a single rotating entity even though the central PS changes. We therefore introduced the notion of a compound rotor, which allows multiple PSs to be associated with a single rotor, and described an objective and automatic algorithm for identifying these entities. In essence, this algorithm allows small-scale fragmentation and collision events near rotor tips to be neglected when associating PSs and wave fronts with rotors. This guards against such events (which may be noise artifacts) from resulting in inappropriately abbreviated rotor lifetimes. The use of this relaxed definition of a rotor, which has the effect of registering rotor lifetimes that are longer than those for traditional single-PS rotors, is particularly important for the present study because it amplifies our major finding that individual rotor lifetimes are generally shorter than the duration of VF.
Rotor lifetimes. All 17 VF episodes contained compound rotors that lasted for hundreds of milliseconds, and each had at least one compound rotor that lasted for 1 s or more. However, a compound rotor cannot drive the VF episode (i.e., be a mother rotor) unless it persists for the entire episode. Only 1 of the 17 VF episodes we analyzed contained compound rotors that lasted for the entire mapped interval. Previous epicardial mapping studies were unable to document mother rotors in either the left (16, 21, 22) or right (22, 23) ventricle. However, these studies could not exclude epicardial mother rotors because only limited epicardial regions were sampled. In the present study, we used a newly developed optical mapping system (12) that can identify and track VF wave fronts across nearly the entire epicardium. Thus we extend the previous results and exclude stable epicardial rotors (with the possible exception of a rotor at the very tip of the apex) as a requirement for VF maintenance in this experimental model.
The possibility still exists that VF in pig hearts is driven by a stable mother rotor whose filament does not intersect the epicardium and therefore cannot be observed directly with epicardial mapping. However, stable mother rotors were not observed in studies mapping limited nonepicardial regions, including the septum (10) and sections across the right (27) and left (21, 27) ventricular free walls. Further studies are needed to confirm or exclude the presence of stable nonepicardial mother rotors.
Despite the paucity of compound rotors that persisted for the entire mapped interval, shorter-lived rotors were common. We found that at any given instant, roughly a dozen (with lifetimes of at least 200 ms) were simultaneously present on the epicardium. These data suggest that the formation of new rotors is critical for VF maintenance because once started, rotors will eventually die out by themselves. This implies that rotor formation may be a more attractive target than rotor termination in the development of new anti-VF therapies.
Limitations.
In the present study, we imaged VF activation patterns with high spatial (
1.6 mm) and temporal (750 frames/s) resolution over the great majority of the swine epicardial surface. However, patterns within the ventricular walls were not recorded. The study focused on VF relatively early after induction (Wiggers' stage II). Activation patterns are known to change as VF proceeds (9, 29), so our findings may not apply to more prolonged VF. We mapped isolated perfused swine hearts that were exposed to BDM to eliminate motion artifacts in the optical recordings. Other species were not studied. We previously found that VF patterns in this preparation are slower and more regular than those in the same hearts before excision (18). Other authors have found that BDM converts VF to a stable tachycardia in isolated ventricular tissue from dogs (19) and pigs (13). In addition, Evans et al. (7) used BDM to stabilize reentry to create a model of monomorphic ventricular tachycardia in rabbit hearts (7). Thus, because of the slowing and stabilization effects of heart isolation and BDM on VF, we expect that rotor lifetimes are substantially shorter in intact hearts than they were in the present study.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
Present address of M. W. Kay: Dept. of Physiology and Pharmacology, The George Washington University, Washington, DC.
Present address of J. D. Gladden: Dept. of Medicine, Univ. of Alabama at Birmingham, Birmingham, AL.
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 Supplemental movie S1, available with the online version of this article, shows animation of fluorescence data from one VF episode (episode 6 from Table 1) texture-mapped onto a model of the epicardial geometry. Figure 4 in the main manuscript shows one frame from this sequence. Blue tissue is resting, and red is activated. ![]()
2 Supplemental movie S2, available with the online version of this article, shows wave fronts and PSs identified in the same VF episode shown in Supplemental Movie S1. Wave fronts are cyan. PSs are yellow or red, depending on chirality. The geometry is displayed with a Hammer map projection. ![]()
3 Supplemental movie S3, available with the online version of this article, shows a VF episode with persistent rotors (episode 11 from Table 1). The rotors are located at top right of the image. The format is the same as in Supplemental Movie S2. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. H. Keldermann, K. H. W. J. ten Tusscher, M. P. Nash, R. Hren, P. Taggart, and A. V. Panfilov Effect of heterogeneous APD restitution on VF organization in a model of the human ventricles Am J Physiol Heart Circ Physiol, February 1, 2008; 294(2): H764 - H774. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Nanthakumar, J. Jalife, S. Masse, E. Downar, M. Pop, J. Asta, H. Ross, V. Rao, S. Mironov, E. Sevaptsidis, et al. Optical mapping of Langendorff-perfused human hearts: establishing a model for the study of ventricular fibrillation in humans Am J Physiol Heart Circ Physiol, July 1, 2007; 293(1): H875 - H880. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. H.W.J. Ten Tusscher, R. Hren, and A. V. Panfilov Organization of Ventricular Fibrillation in the Human Heart Circ. Res., June 22, 2007; 100(12): e87 - e101. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Masse, E. Downar, V. Chauhan, E. Sevaptsidis, and K. Nanthakumar Ventricular fibrillation in myopathic human hearts: mechanistic insights from in vivo global endocardial and epicardial mapping Am J Physiol Heart Circ Physiol, June 1, 2007; 292(6): H2589 - H2597. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Rogers, G. P. Walcott, J. D. Gladden, S. B. Melnick, and M. W. Kay Panoramic Optical Mapping Reveals Continuous Epicardial Reentry during Ventricular Fibrillation in the Isolated Swine Heart Biophys. J., February 1, 2007; 92(3): 1090 - 1095. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |