Donate Help Contact The AHA Sign In Home
American Heart Association
Circulation Research
Search: search_blue_button Advanced Search
Circulation Research. 2002;91:339-345
Published online before print August 1, 2002, doi: 10.1161/01.RES.0000031801.84308.F4
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
91/4/339    most recent
01.RES.0000031801.84308.F4v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Choi, B.-R.
Right arrow Articles by Salama, G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Choi, B.-R.
Right arrow Articles by Salama, G.
Related Collections
Right arrow Electrophysiology
Right arrow Arrythmias-basic studies
Right arrow Arrhythmias, clinical electrophysiology, drugs
(Circulation Research. 2002;91:339.)
© 2002 American Heart Association, Inc.


Integrative Physiology

Life Span of Ventricular Fibrillation Frequencies

Bum-Rak Choi, Wonchul Nho, Tong Liu, Guy Salama

From the Department of Cell Biology and Physiology (B.-R.C., T.L., G.S.) and Department of Electrical Engineering (W.N.), University of Pittsburgh, Pittsburgh, Pa.

Correspondence to Guy Salama, PhD, Department of Cell Biology and Physiology, University of Pittsburgh, School of Medicine, S314 Biomedical Science Tower, 3500 Terrace St, Pittsburgh, PA 15261. E-mail gsalama{at}pitt.edu


*    Abstract
up arrowTop
*Abstract
down arrowIntroduction
down arrowMaterials and Methods
down arrowResults
down arrowDiscussion
down arrowReferences
 
The nature and organization of electrical activity during ventricular fibrillation (VF) are important and controversial subjects dominated by 2 competing theories: the wavebreak and the dominant mother rotor hypothesis. To investigate spatiotemporal characteristics of ventricular fibrillation (VF), transmembrane potentials (Vm) were recorded from multiple sites of perfused rabbit hearts using a voltage-sensitive dye and a photodiode array or a CCD camera, and the time-frequency characteristics of Vm were analyzed by short-time fast Fourier transform (FFT) or generalized time-frequency representation with a cone-shaped kernel. The analysis was applied to all pixels to track VF frequencies in time and space. VF consisted of blobs, which are groups of contiguous pixels with a common frequency and an ill-defined shape. At any time t, several VF frequency blobs coexisted in the field of view, and the number of coexisting blobs was on average 5.9±2.1 (n=8 hearts) as they appeared and disappeared discontinuously with time and were not fixed in space. The life span of frequency blobs from birth to either annihilation or breakup to another frequency had a half-life of 0.39±0.13 second (n=4 hearts). The Ca2+ channel blocker nifedipine increased the stability of VF frequencies and reduced the number of frequency blobs progressing to a single frequency. In conclusion, VF consists of dynamically changing frequency blobs, which have a short life span and can be modified by pharmacological interventions, suggesting that VF is maintained by dynamically changing multiple wavelets.


Key Words: ventricular fibrillation • ventricular tachycardia • time-frequency analysis • optical mapping • L-type Ca2+ channel


*    Introduction
up arrowTop
up arrowAbstract
*Introduction
down arrowMaterials and Methods
down arrowResults
down arrowDiscussion
down arrowReferences
 
Ventricular fibrillation (VF) has been linked to the development of vortex-like reentry or spiral waves that have been studied in computer simulations,1,2 tissue slices, and perfused hearts.3,4 In this context, nonstationary (meandering or drifting) spiral waves and/or their turbulence could account for the complex morphology of electrocardiograms (ECGs) in VF. Experimental studies mapped electrical activation using multiple electrodes or optical probes of membrane potential to investigate the mechanisms underlying ECG signals seen in VF. VF has been traditionally investigated by analyzing activation maps, but the complex waveforms recorded in VF and the algorithms used to derive activation maps have made it difficult to interpret the underlying mechanisms. Another approach is to analyze voltage oscillations in the frequency domain using fast Fourier transforms (FFT) and to represent reentrant circuits as sources of periodic wave formation.57 However, FFT spectra must be interpreted with caution because they measure the averaged contribution of each frequency component with no information on their time of occurrence. For instance, simple FFT spectra with a single dominant peak can be generated from continuously appearing and disappearing frequency sources; conversely, complex FFT spectra can be obtained from a single source with an increasing and/or decreasing frequency.

To overcome the limitations of frequency analysis, this study applied time-frequency analysis where membrane potential oscillations from different regions of the heart are transferred to time and frequency domains to visualize the creation and annihilation of each frequency component in VF in time and space.


*    Materials and Methods
up arrowTop
up arrowAbstract
up arrowIntroduction
*Materials and Methods
down arrowResults
down arrowDiscussion
down arrowReferences
 
Heart Preparations
New Zealand White rabbits (female, 1.6 to 2.5 kg, n=8; Myrtle’s Rabbitry, Thompson Station, Tenn) were injected with pentobarbital (35 mg/kg, IV) plus heparin (200 U/kg), and the heart was excised and retrogradely perfused through the aorta with (in mmol/L) NaCl 130, NaHCO3 24, MgCl2 1.0, KCl 4.0, NaH2PO4 1.2, dextrose 20, and CaCl2 1.5, at pH 7.4, and gassed with 95% O2 and 5% CO2. In 5 of the 8 rabbit hearts, the L-type Ca2+ blocker nifedipine (2 µmol/L) was perfused during VF, and the electrical activity was monitored continuously to track changes in the structure of VF caused by nifedipine. This investigation conformed to the current Guide for the Care and Use of Laboratory Animals published by the National Institutes of Health (NIH publication No. 85-23, revised 1996).

Optical Apparatus
Optical action potentials (APs) were recorded from Langendorff-perfused rabbit hearts stained with di4-ANEPPS. Hearts were mounted in a chamber to image the anterior surface with either a photodiode array (16x16 matrix, 1000 frames per second; Hamamatsu, Bridgewater, NJ) or a 14-bit digital CCD camera (72x78 pixels, 1000 frames per second; FastOne, PixelVision, Tigard, Oreg). The optical apparatus, chamber, computer interface, and software were previously described.7 Each diode viewed a 1.1x1.1-mm2 area, and each pixel on the CCD camera viewed a 0.055x0.055-mm2 area. The integration time of CCD camera was set to 0.2 ms to reduce possible motion blurring caused by rapid wave propagation. VF was induced by burst stimulation for 2 seconds; voltage recordings were taken continuously for 5 minutes with the photodiode array after the onset of VF (n=8). The recordings of the CCD camera were taken for 8-second intervals from different hearts (n=4). In the subgroup of experiments with nifedipine (n=5), electrical activity was continuously mapped for 5 minutes after the addition of the drug.

Time-Frequency Analysis
To analyze Vm signals during VF in time and frequency domains, a spectrogram (short-time Fourier transform) was generated for each pixel by calculating the FFT spectrum for a brief gaussian window of 1.5 to 2.0 seconds, then shifting the window stepwise in time ({Delta}t=1 ms) and remeasuring the FFT spectrum at successive t intervals. The frequency resolution is 0.5 to 0.66 Hz and is inversely related to the window size. Transmembrane potential oscillations during VF events are thus transformed to time and frequency domains to visualize the evolution of VF frequencies.

The spectral peaks in time-frequency representations can be enhanced using Cohen’s bilinear class,8 convolution of the Wigner distribution with a smoothing kernel. A generalized time-frequency representation C(t, f;{Phi}) of the signal x(t) with the kernel {Phi}(t, {tau}) is as follows8:

equation


where t is time, {tau} is the delay in the autocorrelation function,

equation


and the asterisk (*) denotes complex conjugation. This study applied a cone-shaped kernel, 9 as follows:

equation


The cone-shaped kernel has a constraint on a cone-shaped nonzero region in t and {tau} domains to maintain good resolution for fast-changing spectral peaks and to preserve the onset of time at which the signals arise.9 An efficient algorithm to calculate discrete data was previously reported.9 The cone-shaped smoothing kernel resulted in considerably sharper frequency resolution compared with the standard spectrogram, as previously applied to ECG recordings.10

The time-frequency distribution is represented as a contour map where the horizontal axis represents time and the vertical axis represents frequency. The corresponding optical traces are shown on top of the frequency distribution and overall power spectra of optical traces are shown on the left. The contour lines are drawn every 12.5% of maximum, and the higher the intensity the darker the color. The time-frequency representation with the cone-shaped kernel was filtered further in the time domain to reduce small high-frequency noise using a low-pass Butterworth filter (cutoff <60 Hz).

Segmentation of VF Frequency Blobs
Time-frequency analysis was extended to all 252 channels (256-4 channels at each corner), and the distribution of VF frequencies was analyzed in 4 dimensions, space (x, y), frequency (Hz), and time (t). Because of the massive computational and memory requirements, time-frequency analysis for all pixels was subdivided into 2-minute intervals. The distribution of VF frequencies P(x,y, f) at time t is visualized as an isosurface plot in space (x, y) and frequency (z) axis where P(x, y,f) is above an arbitrary threshold (see Figure 1A). The threshold ({approx}12.5%) was chosen from the histogram of time-frequency distribution so that only the background noise can be eliminated without losing small-frequency peaks.



View larger version (19K):
[in this window]
[in a new window]
 
Figure 1. A, Frequency distribution. The frequency distribution on the epicardium can be visualized as blobs of 3D isosurface plot in space (x-, y-axis) and frequency (f-axis). The connected region is given a unique identifier using a blob-coloring algorithm (see Materials and Methods for detail). B, Tracking VF frequencies in space and time. VF frequencies are traced as an object in 4 dimensions (2 spatial, frequency, and time), and VF frequencies near 16 Hz are shown in space (x-, y-axis) and time (horizontal axis). The life span was calculated from the initial time that a blob appears to the time when it either disappears or splits.

The isosurface plots show that frequency peaks are not uniformly distributed but that certain frequencies appear in small contiguous regions of ill-defined shape. We define this contiguous object in space and frequency domains as a VF frequency blob, because a common approach to identify such objects uses a blob-coloring algorithm. Individual VF frequency blobs at time t are given a unique identification number (region index) as follows. The function describing the distribution of VF frequencies, P(x,y,f), is set to 0 when the pixel value P(x,y,f) is below threshold. For pixels above threshold, each pixel is assigned the same unique identification number as its neighboring pixels when the neighboring pixels are also above threshold (nonzero). VF frequency blobs are then formed by the pixels that are contiguous in space and frequency, with the same identification number (see Figure 1A).

Life Span of VF Frequency Blobs
A similar algorithm can be applied to calculate the life span of VF frequency blobs but now in 4 dimensions to include x, y, frequency, and time t. Blobs connected in x, y, f, and t dimensions are given the same identification number, and each blob was tested to determine whether it subdivided into several blobs in space and frequency domains. Therefore, the life span of VF frequency blobs can be measured until they either split or disappear (see Figure 1B). The life span of VF frequency blobs was calculated from 4 rabbit hearts (2-minute interval of VF) such that >500 VF frequency blobs were identified, tracked in time and space to calculate their life span. A histogram of life span was generated and fitted to the first-order exponential decay curve (y=Ae-t/{tau}) by nonlinear least-squares fitting using Microcal Origin 6.0 (Microcal Software, Inc, Northampton, Mass), and half-life was calculated from the exponential decay constant.


*    Results
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMaterials and Methods
*Results
down arrowDiscussion
down arrowReferences
 
FFT Spectra From Photodiode Array and CCD Camera
Electrical activity measured during VF was found to have complex multiple component FFT spectra when recorded optically with a photodiode array.7,11 The interpretation of multiple-component FFT spectra was challenged because of the lower spatial resolution of the photodiode array compared with the CCD sensors.12 Figure 2A illustrates the field of view of the 16x16 photodiode array and the image (72x78) recorded by the CCD camera. Figure 2B compares the optical recording from a diode from an array (top) to a pixel recording from the CCD camera (bottom) and their FFT spectra (Figure 2C). Voltage transients and FFT spectra from a 0.8x0.8-mm2 area of the heart (top trace) were equally complex as those recorded at a higher spatial resolution from a 50x50-µm2 area (bottom trace) from a single CCD pixel.



View larger version (35K):
[in this window]
[in a new window]
 
Figure 2. Spectral analysis of optical traces in VF. A, Schematic of the heart and the field of view of the photodiode array (top) and an image of the heart taken with the CCD camera (bottom). B, Optical traces from a photodiode on the array (top) and a CCD pixel (bottom) with their associated FFT spectra in panel C. Integration time was set at 0.2 ms in the CCD camera to reduce blurring effects of fast propagating waves.

Time-Frequency Representation of VF Signals
To determine the time when individual FFT peaks occurred, voltage oscillations were represented in time and frequency domains by generating a standard spectrogram (see Materials and Methods). Transmembrane potential oscillations during VF are thus transformed to time and frequency domains to visualize the evolution of VF frequencies. Figure 3A illustrates a spectrogram from a single photodiode where the frequency peaks tend to appear and disappear rapidly in time. Note that the bands of VF frequencies are broad and lack sharp boundaries. The poor resolution is due to an inherent trade-off between time and frequency in the spectrogram algorithm that is best suited for slowly varying frequencies.



View larger version (40K):
[in this window]
[in a new window]
 
Figure 3. Time-frequency analysis of VF. A, Spectrogram. The spectrogram was calculated by sliding a 1.5-second gaussian window over normalized optical traces. Top, Optical trace. Left, Overall FFT spectra. Contour map, Spectrogram with isolines drawn every 12.5% of maximum. Spectrogram plots frequency (ordinate) versus time (abscissa) and is shown for 14 seconds of VF; the darker the color, the higher the energy density at that frequency. B, Time-frequency distribution using a cone-shaped smoothing kernel. C, Time-frequency distribution from one pixel of CCD data using a cone-shaped kernel.

In Figure 3B, the frequency peaks in the spectrogram were enhanced by applying the general time-frequency distribution analysis with a cone-shaped kernel9 (see Materials and Methods for detail). The cone-shaped kernel (Figure 3B) was applied to the same signal to produce a contour map with similar general features but considerably sharper frequency resolution compared with the standard spectrogram (Figure 3A), as previously described. 10 Enhanced time-frequency representation shows separate discontinuous bands, suggesting that VF frequency peaks appeared and disappeared abruptly.

A contour map of time-frequency representation in Figure 3C is obtained from a CCD signal after applying the cone-shaped kernel and shows the same general features as from a photodiode, indicating that the multicomponents in the time-frequency domain were not due to lower spatial resolution of the photodiode array.

Spatial and Temporal Characteristics of VF Frequencies
In many intervals of VF (see Figure 3), VF frequencies were not stable but appeared and disappeared, forming discrete bands. To study the dynamics of VF frequencies along the surface of the heart, time-frequency analysis was extended to all the pixels and reconstructed in space, time, and frequency domains. Figure 4A shows snapshots of VF frequency distribution as isosurface plots (see Materials and Methods) by plotting frequencies (z-axis) against position on the heart (x-, y-axis). The results show that at any given time t, several frequencies coexisted, occupying different regions. It should be noted that successive isosurface plots (Figure 4A) are significantly different from each other.



View larger version (35K):
[in this window]
[in a new window]
 
Figure 4. Spatial distribution of VF frequencies. A, Time-frequency analysis was extended to all photodiodes, and the spatial distribution of VF frequencies at time t was reconstructed as isosurface plots. Note that their spatial locations are not fixed but change with time. B, Number of VF frequency blobs versus time. VF frequency blobs at any given time (gray) and averaged number of VF frequency blobs (red) are plotted as a function of time for a 2-minute interval of VF. C, Life span of VF frequency blobs. The majority of frequency blobs lasted <1 second (half life=0.39±0.13 second).

VF frequency distribution consisted of separate blobs occupying contiguous regions of epicardium, and these were assigned a unique identifier using a blob-coloring algorithm, as described in Materials and Methods. We further investigated how many VF frequency blobs coexist during VF and how long the frequency blobs exist. A plot of the number of frequency blobs versus time (Figure 4B) shows the dynamic nature of VF frequencies and that 2 to 12 blobs (average 5.9±2.1, n=8) coexist at any given time in the field of view. Each blob was traced in time to investigate its life span, as described in Materials and Methods. The life span of individual VF frequencies from their time of birth to either their death or subdivision to several VF frequency blobs was measured and the number of VF blobs with a life span >= time t was displayed as a histogram (Figure 4C). The majority of VF frequencies lasted <1 second, with the longest time of 8 seconds in this episode of VF. The life span plot was curve-fitted to a single-exponential decay function with a time constant {tau}=0.57±0.19 second (n=4, half-life=0.39 second). VF frequencies were spread over an area of 32±34 mm2 during their life time, where the larger the area occupied by a frequency blob, the longer the lifetime of the frequency blob with correlation coefficient of 0.59 (P<0.0001).

Effect of L-Type Ca2+ Channel Blocker on VF Frequencies
L-type Ca2+ channel blockers such as verapamil have been shown to change VF to ventricular tachycardia (VT).13,14 Time-frequency analysis was applied to better understand and visualize how a Ca2+ channel blocker influences the behavior of VF dynamics in time and space. Nifedipine (2 µmol/L) was added during VF, and transmembrane potentials were recorded continuously for 5 minutes for time-frequency analysis. Figure 5A shows snapshots of VF frequency isosurface plots recorded during perfusion of the heart with nifedipine. The complex spatial distribution of VF frequencies was consistently converted to a large stable frequency blob by nifedipine (n=5 of 5 hearts). In this experiment, nifedipine converted VF to VT, which appeared at steady state, as a stable doughnut. The activation map of corresponding isosurface plots is shown in Figure 5B. Time-frequency plots from a single pixel in Figure 5C also showed that the complex dynamics of VF were converted to a single frequency on the addition of nifedipine. Figure 5D shows that the number of frequency blobs was gradually reduced by nifedipine, showing that the shift from VF to VT was due to a reduction of the number of VF frequencies that coexisted in the field of view.



View larger version (46K):
[in this window]
[in a new window]
 
Figure 5. Dynamics of VF during perfusion with nifedipine. A, Snapshots of isosurface plot during VF and its conversion to VT by adding nifedipine (2 µmol/L) to the perfusate. Complex multicomponent frequency distribution gradually converted to a single dominant frequency. B, Activation map of corresponding doughnut-shaped frequency distribution. Isochronal lines are 2 ms apart. C, Time-frequency analysis from one pixel. D, Number of VF frequency blobs versus time during nifedipine perfusion.


*    Discussion
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMaterials and Methods
up arrowResults
*Discussion
down arrowReferences
 
Technological advances now provide high-resolution maps of electrical activity that have been used to investigate the structure of VF using frequency analysis to correlate a particular frequency with a source of reentry.5,6 However, FFT analysis does not provide information regarding time-varying frequencies, which is essential to understand the life and death of VF frequencies in time and space. To overcome this limitation, we analyzed signals in time, frequency, and space at high spatial and temporal resolution, making it possible to track VF dynamics in great detail. Major findings are that VF is composed of several VF frequency blobs that have ill-defined shapes, are unstable in space and time, and are created, subdivided, and destroyed with a brief life span. The dynamics of VF frequency blobs can be modified by a blocker of voltage-gated Ca2+ channel, suggesting that VF is maintained by the continuous creation and annihilation of periodic sources.

Reentrant wave propagation has been thought to be the major mechanism underlying complex transmembrane oscillations during VF, but the maintenance of VF is still under debate. As early as 1914, Garrey 15,16 suggested that atrial fibrillation (AF) is caused by reentry and that the multiple waves maintain AF. However, Lewis17 proposed that AF is caused by a single, rapidly rotating central wave, and its rapid activation causes an irregular impulse propagation across partially refractory tissue (see review18). Nearly 80 years later, the hypotheses first proposed by Garrey and Lewis are still controversial and have been expanded to multiple wavelets and mother rotor hypotheses (see reviews19,20). According to the multiple wavelets hypothesis, continuous wavebreaks create meandering fragmented wavelets that cause complex meandering, annihilation, and creation of new wavelets.

However, the mother rotor hypothesis suggests that a high-frequency stable rotor generates fragmented conduction or fibrillatory conduction. The mother rotor hypothesis was based on the observation that electrical activity during VF recorded optically with a CCD camera exhibited dominant frequency peaks that were relatively stable and distributed in discrete regions of the heart.5,6 In smaller rabbit and guinea pig hearts, a single peak frequency was found to cover the entire left ventricle, supporting the mother rotor hypothesis. In contrast, Wiggers21 reported a considerably more complex frequency distribution in hearts under VF with marked instabilities and heterogeneities of frequencies across the epicardium. Several groups2224 have reproduced these results and demonstrated that frequencies are spatially heterogeneous and unstable in time, which is opposite to the premise of the mother rotor hypothesis. Furthermore, FFT spectra were more complex, showing multiple peaks instead of an FFT with a single peak and a narrow bandwidth. It is important to emphasize that the detection of multiple frequency peaks requires an analysis with high spectral resolution, which in turn requires longer time intervals of FFT analysis.24

A major limitation of FFT analysis is that it fails to provide temporal information of when a specific FFT peak (or VF source) occurs in time. One possibility is that multiple frequency peaks arise from multiple frequency sources. Alternatively, multiple frequency peaks can arise from a single rotor if the amplitude of the signals oscillates or if frequency signals change gradually, as can occur by the Doppler effect from the drifting of a single rotor. Therefore, VF signals should be investigated in time and frequency domains, and the analysis should be extended to the multiple locations to obtain a complete picture of VF dynamics.

Generalized time-frequency distribution based on Cohen’s bilinear class has been applied to study time-varying biological signals including AF and VF.23,25,26 For example, using grids of extracellular electrodes, Moghe et al23 showed that VF frequencies are unstable, varying with time and location, in different regions of the fibrillating heart. Patwardhan et al 26 performed a time-frequency analysis on ECG signals and tested the cross correlations between different leads to determine whether the time-frequency distribution in different regions of the heart can be explained by a time shift. They found that after applying a time shift to one ECG lead, they could not improve the cross correlation with respect to another lead, such that the different VF frequencies from different leads could not be attributed to the same frequency source. These findings were not consistent with the existence of a single large, drifting rotor. However, these studies were limited to the number of electrodes that could be simultaneously recorded and by the difficulties of interpreting the surface electrograms during VF.

We applied time-frequency analysis to the surface of hearts with high spatial and temporal resolution. Optical recordings of voltage changes, whose shape is equivalent to intracellular microelectrode measurements, were recorded from the anterior surface of the heart. Signals were analyzed in time, frequency, and space at high spatial and temporal resolution, making it possible to track VF frequency peaks in great detail. It should be noted that signals were recorded for a long period of time (5 minutes). Data analysis was performed in 2-minute segments of VF recordings, owing to a limitation of computing power. Nevertheless, the period of analysis is markedly longer than in previous studies, in which only a few periods of the dominant frequency were analyzed.

Our time-frequency analysis found discrete bands, suggesting that VF frequencies are not fixed but change abruptly in time. It is possible that gradual changes in frequency may appear as step changes in cases with poor spectral resolutions. However, in most cases, the step changes between discrete frequency bands exceeded the frequency resolution (0.5 Hz), and the amplitude of the signals between frequency bands is negligible, such that the discrete bands are most likely due to abrupt changes in VF frequencies. The brief life span and the spatially discontinuous appearance and disappearance of VF frequency blobs suggest that VF is maintained by the seemingly chaotic creation and annihilation of periodic sources.

The general consensus is that the source of a VF frequency is a reentrant circuit or a rotor; however, the relationship between frequency and its underlying rotor has not been demonstrated rigorously. Computer simulations and mapping studies2732 found complex activation patterns with wavefronts of varying amplitude and activation intervals. However, complete reentrant circuits are rarely delineated from activation maps. Previous analysis of activation patterns in an open-chest canine model of VF showed reentrant circuits with a short life span (<0.5 second per wavefront).33,34 This study used an array of bipolar electrodes and is in excellent agreement with our life-span measurements based on time-frequency analysis and optical measurements of electrical activity on an isolated rabbit heart. It is reasonable, therefore, to speculate that VF frequency blobs identified by frequency distribution analysis represent reentrant circuits with virtually identical life spans. However, it is possible that single VF frequency blobs are composed of several rotors with the same frequency. Alternatively, frequency blobs can be due to a complex meandering reentrant circuit. Previous studies5,7,23,35 indicated that the action potential duration (APD) and refractory period at each site influences the VF frequencies and activation intervals recorded from the same site. Hence, conditions that reduce the dispersion of APDs are expected to depress the number of VF frequency blobs because either a single dominant rotor is formed or multiple rotors in different location share the same frequency. Further studies may reveal the correlation between wave propagation patterns and VF frequencies using a combination of wave pattern and time-frequency analysis.

Theoretical studies predicted that wave breakup can be caused by spontaneous wavelength oscillations of reentry circuits.36,37 For instance, APDs exhibit a dynamic adaptation to rate changes and diastolic intervals (period from the recovery of the last AP to the onset of the next AP), a phenomenon called restitution kinetics. Experimental studies showed a marked correlation between the complexity of voltage oscillations in VF and the slope of the restitution curve.20,34,38,39 When the slope of the restitution curve is >1, small perturbations in the APD become amplified with successive reentrant cycles, and when the APD or wavelength is too short, conduction fails, resulting in wavebreak. The slope of the restitution curve has therefore been suggested as a major predictor of wave breakup. Of interest would be to look for a correlation between restitution and life span of VF frequencies.

In conclusion, the analysis provides new insights regarding the birth and evolution of periodic sources during VF, which are tracked in space and time. We found that VF frequencies are not anchored to certain regions on the epicardium and are not stable in time, lasting <1 second before being annihilated, suggesting that VF is maintained by multiple wavelets as a result of continuous wave breakups. The visualization of VF dynamics in space, time, and frequency domains in various pathological conditions (heart failure, ischemia, and long-QT-related arrhythmias) can provide new insights of the nature of these arrhythmias and new strategies for their interruption.

Limitations of the Study
A small change in VF frequency, which is due to changes in several frequency sources, may appear as a single VF frequency blob in cases of poor frequency resolution, causing an overestimation of life span. Although time-frequency representation with a cone-shaped kernel provided a marked improvement of peak resolution, further improvements could be achieved with more sophisticated algorithms. In the future, increasing computer power may make it possible to calculate optimized time-frequency representations for VF in relatively brief computer times.

Optical recordings were limited to the anterior region of the epicardium. Therefore, it is possible that stable frequency sources in time and space can still exist transmurally. Further studies may reveal whether a stable frequency source can exist in different regions.


*    Acknowledgments
 
This work was supported by grant awards from the NIH: R01 HL57929 and HL59614 to G. Salama and a postdoctoral fellowship from the Western Pennsylvania Affiliate of the American Heart Association to Drs B.-R. Choi and T. Liu.

Received June 4, 2002; revision received July 17, 2002; accepted July 22, 2002.


*    References
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMaterials and Methods
up arrowResults
up arrowDiscussion
*References
 
1. Winfree AT. Electrical instability in cardiac muscle: phase singularities and rotors. J Theor Biol. 1989; 138: 353–405.[CrossRef][Medline] [Order article via Infotrieve]

2. Krinsky VI. Mathematical models of cardiac arrhythmias (spiral waves). Pharmacol Ther [B]. 1978; 3: 539–555.[Medline] [Order article via Infotrieve]

3. Pertsov AM, Davidenko JM, Salomonsz R, Baxter WT, Jalife J. Spiral waves of excitation underlie reentrant activity in isolated cardiac muscle. Circ Res. 1993; 72: 631–650.[Abstract/Free Full Text]

4. Davidenko JM, Pertsov AM, Salomonsz R, Baxter WT, Jalife J. Stationary and drifting spiral waves of excitation in isolated cardiac muscle. Nature. 1991; 355: 349–351.[CrossRef]

5. Zaitsev AV, Berenfeld O, Mironov SF, Jalife J, Pertsov AM. Distribution of excitation frequencies on the epicardial and endocardial surfaces of fibrillating ventricular wall of the sheep heart. Circ Res. 2000; 86: 408–417.[Abstract/Free Full Text]

6. Chen J, Mandapati R, Berenfeld O, Skanes AC, Jalife J. High-frequency periodic sources underlie ventricular fibrillation in the isolated rabbit heart. Circ Res. 2000; 86: 86–93.[Abstract/Free Full Text]

7. Choi BR, Liu T, Salama G. The distribution of refractory periods influences the dynamics of ventricular fibrillation. Circ Res. 2001; 88: e49–e58.[Medline] [Order article via Infotrieve]

8. Cohen L. Time-frequency distributions: a review. Proceedings of the IEEE: Institute of Electrical and Electronics Engineers. 1989; 77: 941–981.

9. Zhao Y, Atlas L, Marks R. The use of cone-shaped kernels for generalized time-frequency representations of nonstationary signals. IEEE Trans Acoust. 1990; 38: 1084–1091.[CrossRef]

10. Millet-Roig J, Rieta-Ibanez JJ, Vilanova E, Mocholi A, Chorro FJ. Time-frequency analysis of a single ECG: to discriminate between ventricular tachycardia and ventricular fibrillation. Comput Cardiol. 1999; 26: 711–714.

11. Choi BR, Liu T, Salama G. Ventricular fibrillation: mother rotor or multiple wavelets? Circ Res. 2001; 89: e30.Letter.[Medline] [Order article via Infotrieve]

12. Berenfeld O, Pertsov AM, Jalife J. What is the organization of waves in ventricular fibrillation? Circ Res. 2001; 89: e22.Letter.[Medline] [Order article via Infotrieve]

13. Watanabe Y, Uchida H. Verapamil-induced sustained ventricular tachycardia in isolated, perfused rabbit hearts. Jpn Circ J. 1987; 51: 188–195.[Medline] [Order article via Infotrieve]

14. Merillat JC, Lakatta EG, Hano O, Guarnieri T. Role of calcium and the calcium channel in the initiation and maintenance of ventricular fibrillation. Circ Res. 1990; 67: 1115–1123.[Abstract/Free Full Text]

15. Garrey WE. Auricular fibrillation. Physiol Rev. 1924; 4: 215–250.[Free Full Text]

16. Garrey WE. The nature of fibrillary contraction of the heart: its relation to tissue mass and form. Am J Physiol. 1914; 33: 215–250.

17. Lewis T. Observations upon flutter and fibrillation, part IV: impure flutter; theory of circus movement. Heart. 1920; 7: 293–331.

18. Wit AL, Cranefield PF. Reentrant excitation as a cause of cardiac arrhythmias. Am J Physiol. 1978; 235: H1–H17.[Medline] [Order article via Infotrieve]

19. Samie FH, Jalife J. Mechanisms underlying ventricular tachycardia and its transition to ventricular fibrillation in the structurally normal heart. Cardiovasc Res. 2001; 50: 242–250.[Abstract/Free Full Text]

20. Weiss JN, Chen P-S, Qu Z, Karagueuzian HS, Lin S-F, Garfinkel A. Electrical restitution and cardiac fibrillation. J Cardiovasc Electrophysiol. 2002; 13: 292–295.[CrossRef][Medline] [Order article via Infotrieve]

21. Wiggers CJ. The mechanisms and nature of ventricular fibrillation. Am Heart J. 1940; 20: 399–412.[CrossRef]

22. Choi B-R, Salama G. Spatial distribution of activation intervals and conduction velocities (CV) during ventricular fibrillation (VF). Pacing Clin Electrophysiol. 2000; 23: 609.Abstract.

23. Moghe SA, Qu F, Leonelli FM, Patwardhan AR. Time-frequency representation of epicardial electrograms during ventricular fibrillation. Biomed Sci Instrum. 2000; 36: 45–50.[Medline] [Order article via Infotrieve]

24. Valderrábano M, Yang J, Omichi C, Kil J, Lamp ST, Qu Z, Lin SF, Karagueuzian HS, Garfinkel A, Chen PS, Weiss JN. Frequency analysis of ventricular fibrillation in swine ventricles. Circ Res. 2002; 90: 213–222.[Abstract/Free Full Text]

25. Stridh M, Sornmo L, Meurling CJ, Olsson SB. Characterization of atrial fibrillation using the surface ECG: time-dependent spectral properties. IEEE Trans Biomed Eng. 2001; 48: 19–27.[CrossRef][Medline] [Order article via Infotrieve]

26. Patwardhan A, Moghe S, Wang K, Leonelli F. Frequency modulation within electrocardiograms during ventricular fibrillation. Am J Physiol Heart Circ Physiol. 2000; 279: H825–H835.[Abstract/Free Full Text]

27. Chialvo DR, Gilmour RF, Jalife J. Low dimensional chaos in cardiac tissue. Nature. 1990; 343: 653–657.[CrossRef][Medline] [Order article via Infotrieve]

28. Garfinkel A, Chen PS, Walter DO, Karagueuzian HS, Kogan B, Evans SJ, Karpoukhin M, Hwang C, Uchida T, Gotoh M, Nwasokwa O, Sager P, Weiss JN. Quasiperiodicity and chaos in cardiac fibrillation. J Clin Invest. 1997; 99: 305–314.[Medline] [Order article via Infotrieve]

29. Huang J, Rogers JM, Kenknight BH, Rollins DL, Smith WM, Ideker RE. Evolution of the organization of epicardial activation patterns during ventricular fibrillation. J Cardiovasc Electrophysiol. 1998; 9: 1291–1304.[Medline] [Order article via Infotrieve]

30. Moe G, Rheinboldt W, Abildskov J. A computer model of atrial fibrillation. Am Heart J. 1964; 67: 200–220.[CrossRef][Medline] [Order article via Infotrieve]

31. Qu Z, Weiss JN, Garfinkel A. Cardiac electrical restitution properties and stability of reentrant spiral waves: a simulation study. Am J Physiol. 1999; 276: H269–H283.[Medline] [Order article via Infotrieve]

32. Rogers JM, Huang J, Smith WM, Ideker RE. Incidence, evolution, and spatial distribution of functional reentry during ventricular fibrillation in pigs. Circ Res. 1999; 84: 945–954.[Abstract/Free Full Text]

33. Cha YM, Birgersdotter-Green U, Wolf PL, Peters BB, Chen PS. The mechanism of termination of reentrant activity in ventricular fibrillation. Circ Res. 1994; 74: 495–506.[Abstract/Free Full Text]

34. Garfinkel A, Kim YH, Voroshilovsky O, Qu Z, Kil JR, Lee MH, Karagueuzian HS, Weiss JN, Chen PS. Preventing ventricular fibrillation by flattening cardiac restitution. Proc Natl Acad Sci U S A. 2000; 97: 6061–6066.[Abstract/Free Full Text]

35. Wang L, Li CY, Yong AC, Kilpatrick D. Fast Fourier transform analysis of ventricular fibrillation intervals to predict ventricular refractoriness and its spatial dispersion. Pacing Clin Electrophysiol. 1998; 21: 2588–2595.[CrossRef][Medline] [Order article via Infotrieve]

36. Karma A. Spiral breakup in model equations of action potential propagation in cardiac tissue. Phys Rev Lett. 1993; 71: 1103–1106.[CrossRef][Medline] [Order article via Infotrieve]

37. Karma A. Electrical alternans and spiral wave breakup in cardiac tissue. Chaos. 1994; 4: 461–472.[CrossRef][Medline] [Order article via Infotrieve]

38. Riccio ML, Koller ML, Gilmour RF Jr. Electrical restitution and spatiotemporal organization during ventricular fibrillation. Circ Res. 1999; 84: 955–963.[Abstract/Free Full Text]

39. Gilmour RF, Chialvo DR. Electrical restitution, critical mass, and the riddle of fibrillation. J Cardiovasc Electrophysiol. 1999; 10: 1087–1089.[Medline] [Order article via Infotrieve]




This article has been cited by other articles:


Home page
J. Physiol.Home page
B. London, L. C. Baker, P. Petkova-Kirova, J. M. Nerbonne, B.-R. Choi, and G. Salama
Dispersion of repolarization and refractoriness are determinants of arrhythmia phenotype in transgenic mice with long QT
J. Physiol., January 1, 2007; 578(1): 115 - 129.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
A. Baher, Z. Qu, A. Hayatdavoudi, S. T. Lamp, M.-J. Yang, F. Xie, S. Turner, A. Garfinkel, and J. N. Weiss
Short-term cardiac memory and mother rotor fibrillation
Am J Physiol Heart Circ Physiol, January 1, 2007; 292(1): H180 - H189.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
M. W. Kay, G. P. Walcott, J. D. Gladden, S. B. Melnick, and J. M. Rogers
Lifetimes of epicardial rotors in panoramic optical maps of fibrillating swine ventricles
Am J Physiol Heart Circ Physiol, October 1, 2006; 291(4): H1935 - H1941.
[Abstract] [Full Text] [PDF]


Home page
CirculationHome page
R. E. Ideker and J. M. Rogers
Human Ventricular Fibrillation: Wandering Wavelets, Mother Rotors, or Both?
Circulation, August 8, 2006; 114(6): 530 - 532.
[Full Text] [PDF]


Home page
CirculationHome page
M. P. Nash, A. Mourad, R. H. Clayton, P. M. Sutton, C. P. Bradley, M. Hayward, D. J. Paterson, and P. Taggart
Evidence for Multiple Mechanisms in Human Ventricular Fibrillation
Circulation, August 8, 2006; 114(6): 536 - 542.
[Abstract] [Full Text] [PDF]


Home page
CirculationHome page
T. H. Everett IV, E. E. Wilson, S. Foreman, and J. E. Olgin
Mechanisms of Ventricular Fibrillation in Canine Models of Congestive Heart Failure and Ischemia Assessed by In Vivo Noncontact Mapping
Circulation, September 13, 2005; 112(11): 1532 - 1541.
[Abstract] [Full Text] [PDF]


Home page
CirculationHome page
S. P. Thomas, A. Thiagalingam, E. Wallace, P. Kovoor, and D. L. Ross
Organization of Myocardial Activation During Ventricular Fibrillation After Myocardial Infarction: Evidence for Sustained High-Frequency Sources
Circulation, July 12, 2005; 112(2): 157 - 163.
[Abstract] [Full Text] [PDF]


Home page
Circ. Res.Home page
I. R. Efimov, V. P. Nikolski, and G. Salama
Optical Imaging of the Heart
Circ. Res., July 9, 2004; 95(1): 21 - 33.
[Abstract] [Full Text] [PDF]


Home page
CirculationHome page
P.-S. Chen, T.-J. Wu, C.-T. Ting, H. S. Karagueuzian, A. Garfinkel, S.-F. Lin, and J. N. Weiss
A Tale of Two Fibrillations
Circulation, November 11, 2003; 108(19): 2298 - 2303.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
91/4/339    most recent
01.RES.0000031801.84308.F4v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Choi, B.-R.
Right arrow Articles by Salama, G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Choi, B.-R.
Right arrow Articles by Salama, G.
Related Collections
Right arrow Electrophysiology
Right arrow Arrythmias-basic studies
Right arrow Arrhythmias, clinical electrophysiology, drugs