The human heart pumps oxygenated blood to the organs and extremities in order to maintain normal physiologic function, while simultaneously pumping deoxygenated blood to the lungs for reoxygenation. Coordinated contraction of individual cardiac myocytes provides the mechanical force necessary to produce sufficient pressure and ensure that distant organs and extremities remain oxygenated. Before cardiac myocytes may contract, they must undergo excitation in order to begin the sequence of events which results in an intracellular calcium (Cai) rise, which in turn precipitates actin-myosin binding and ultimately results in contraction. The electrical signature of this series of events is reflected in the cardiac action potential (AP), a segment of a transmembrane voltage (Vm) recording which indicates electrical excitation (depolarization) and relaxation (repolarization) of the myocardium.
The duration, amplitude, upstroke velocity (dVm/dt), and overall morphology of the cardiac AP are important markers of the electrical status of the heart. Studies of the cardiac AP have provided important insights into the mechanisms which drive the transition from a normal, healthy heartbeat toward a deadly cardiac arrhythmia.
Early recordings of the cardiac AP were obtained using microelectrodes (Coraboeuf & Weidmann, 1949a; Coraboeuf & Weidmann, 1949b; Draper & Weidmann, 1951; Sano et al., 1959; Sano et al., 1960; Weidmann, 1951). Although this method was highly effective in tracking temporal changes in the Vm of individual cells, the method could not be easily applied to the problem of tracking excitation over a region of tissue. Extracellular electrode mapping offered a partial solution to this problem and was sufficient to determine activation times in regions of tissue, but with this method the details of repolarization were lost and had to be estimated using indirect indicators. Further, this method required that the electrodes be in direct contact with the tissue. This made defibrillation studies difficult, since large amplitude defibrillation shocks typically obscure the details of activation during electrical recordings. Monophasic action potential (MAP) recordings were capable of elucidating the details of repolarization without damaging tissue, and have even been recorded in the beating human heart using a cardiac catheter (Shabetai et al., 1968). However they too were restricted by having little or no spatial resolution and could not be placed in close contact with each other. As with extracellular electrodes, MAP recordings also require that the electrodes be placed in contact with the tissue.
With the emergence of Vm-sensitive dyes in the 70’s, it became possible to interrogate cardiac tissue optically (Salama, 1976), and soon afterward optical methods were developed to interrogate multiple spots simultaneously in a small (~cm2) region of tissue. Since then the field of cardiac optical mapping (COM) has greatly expanded in scope, from relatively simple early recordings using one or relatively few spots (Morad & Dillon, 1981; Salama, 1976) to highly complex optical systems. These include high spatiotemporal resolutions systems (Choi et al., 2007), panoramic systems (Kay et al., 2004; Rogers et al., 2007), and systems which are capable of interrogating electrophysiological activity beneath the surface (Byars et al., 2003). In addition, several labs have used photodiode-based optical mapping systems to map Vm and Cai simultaneously, on both the whole heart (Choi & Salama, 2000; Lakireddy et al., 2006; Laurita & Singal, 2001; Pruvot et al., 2004) and in monolayer cell cultures of cardiac myocytes (Fast, 2005; Fast & Ideker, 2000; Lan et al., 2007).
Cardiac optical mapping systems have greatly increased our understanding in nearly all areas of cardiac electrophysiology, from basic studies of conduction patterns (Cabo et al., 1994; Knisley & Hill, 1995) and effects of fiber geometry (Knisley & Baynham, 1997; Knisley et al., 1994; Knisley et al., 1999; Neunlist & Tung, 1995) to more clinical studies of defibrillation (Al-Khadra et al., 2000; Fast et al., 2002; Federov et al., 2008; Tung & Cysk, 2007) and ablation therapy (Himel et al., 2007; Perez et al., 2006). Although COM has not yet led to a widely accepted method of three-dimensional cardiac tissue interrogation, there have been significant advances in this area as well. Investigators have successfully used optical surface recordings to determine wavefront orientation beneath the surface (Hyatt et al., 2005; Zemlin et al., 2008), and also to interrogate deeper layers of tissue using transillumination methods (Baxter et al., 2001) and deeper-penetrating, near-infrared fluorescing dyes (Matiukas et al., 2006; Matiukas et al., 2007; Salama et al., 2005).
Photodiode sensors were used in some of the earliest optical recordings of cardiac APs (Morad & Salama, 1979; Salama, 1976), and continue to be used today (Cheng, 2006; Sakai, 2008). Photodiodes function by transferring incoming photonic energy to bound electrons in a semi-conductive material in a transistor configuration. These energized electrons may then cross from one side of the transistor to the other, resulting in a voltage difference between the two sides.
If a wire is connected from one side of the photodiode to the other while the photodiode is receiving photonic energy, current will flow in a linear fashion with respect to the input intensity of the collected light (Scherz, 2007). This makes photodiodes an excellent choice as a detector in COM systems, and this fact has been reflected by their widespread use over the past 30 years. Examples of optical APs and activation maps recorded with a photodiode array-based system are shown in figure 1.
Although other technologies such as CMOS and CCD cameras have recently gained popularity due to their higher spatial resolution, photodiode systems remain in use due to their ruggedness, high signal-to-noise ratios, excellent temporal resolution, versatility, and low cost. Recently, for example, photodiodes and photodiode arrays (PDAs) have been used in the construction of optrodes, a novel technique used to record optical signals from deeper intramural regions within the ventricular wall (Byars et al., 2003; Caldwell et al., 2005; Hooks et al., 2001; Kong et al., 2007).
2. Basic principles of cardiac optical mapping
Epi-illumination occurs when the fluorescence emission detector is placed on the same side of the tissue as the excitation source, whereas with trans-illumination the detector and excitation source are placed on opposite sides. For monolayer mapping systems, both epi and trans-illumination are possible since cardiac monolayers are typically only a few tens of micrometers thick. For whole-heart mapping systems which map excitation on the surface of the intact heart preparations, epi-illumination is the preferred method since very little fluorescence is transmitted through the relatively thick myocardial wall.
The tissue being mapped must be illuminated using an excitation source, which excites at least one parameter-sensitive dye in order to elicit a fluorescent signal. Changes in a targeted physiological parameter cause changes in the properties of the dye (e.g., a conformational change in the dye molecules). This results in a change in the emission spectrum of the dye, which is then recorded by a detector (e.g., a PDA), digitized, and stored on a PC for post-experimental analysis.
Changes in fluorescence due to changes in the physiological parameter are often measured as a fraction of the baseline fluorescence. This is an important parameter in optical mapping, and is known as fractional fluorescence (ΔF/F). Fractional fluorescence is useful because it is a way to measure the effectiveness of a particular dye in transducing a physiological change into recordable fluorescent signal. Fractional fluorescence also indicates the general effectiveness of the system, and higher ΔF/F values are typically accompanied by higher signal-to-noise ratios. Transmembrane voltage is the most commonly studied physiological parameter in optical mapping, but intracellular calcium transients (CaiT) have also been studied extensively.
There are several variations of the COM system, however there are basic components that are common to all systems. These basic components include an excitation source, detector, and electronic components used for digitization, filtration, and multiplexing. A schematic for a typical whole-heart mapping system is shown in figure 2.
2.1. Excitation source
The excitation source may be either focused (i.e. laser light) or broadfield illumination (using halogen, tungsten, or more recently, high-power LED sources). Laser and some LED light sources have sufficiently narrow bands so as not to interfere with the fluorescence emission, however broadfield sources should be pre-processed using optical filters in order to decrease the width of their wavelength spectrum before illuminating the target (e.g., the heart). In general, brighter excitation sources lead to higher ΔF/F values, however the intensity of the excitation source cannot be increased without regard for
There are several types of detectors that are currently in use for COM, however the focus of this review is upon those detectors which are photodiode-based. Other detector types will be discussed for the purpose of comparison.
Current alternatives to PDAs include photomultiplier (PMT) systems, charge-coupled device (CCD) cameras, and complimentary metal-oxide semiconductor (CMOS) cameras. PMT systems have an extremely high gain (up to 108 increase in intensity), and can even be used in a process call “photon counting” whereby the release of individual photons may be recorded. PMT systems are capable of extremely high sensitivity, and typically have a quicker response time than photodiode systems. However for normal optical mapping applications, the sensitivity and response time of photodiodes is more than sufficient. Typically, PMT systems employ sequential rather than simultaneous recordings of adjacent spots. However, due to their rapid response time they may be coupled with a laser scanner to produce acquisition rates equivalent to those seen in photodiode-based systems. For example, a laser scanner which scans at 256 kHz with a 256-spot grid is essentially equivalent to a PDA-based 256-spot system which scans at a rate of 1 kHz. For this example, there will be a time difference of ~1 ms between the 1st and 256th spots scanned. This difference can often be neglected during normal propagation in cardiac tissue. The basis for the photomultiplier tube is the
Cameras using CCD technology typically have impressive spatial resolution, however these systems have not been able to achieve the same temporal resolution as PDA-based systems. The CardioCCD-SMQ (RedShirt Imaging, Decatur, GA) currently offers a spatial resolution of 80x80 while achieving a temporal resolution of 2000 frames per second (fps). The CCD camera may achieve higher temporal resolution and/or higher signal-to-noise ratios by a process known as “binning,” however this process decreases the spatial resolution of the camera, the major advantage of CCD technology. When using the CCD camera for cardiac mapping experiments, light saturation of the detector may be an issue and thus gain must be carefully controlled. Typically this is less of a problem with PDA systems. The basis of operation of the CCD detector is a shift registration process. Charge which is accumulated due to photons falling upon light-sensitive regions within the detector is transferred along the 2-dimensional detector array until the charge within a given potential well reaches the final readout electrode. The transfer of charge is accomplished by changing the voltage in an adjacent pixel, causing electrical charge to flow in the desired direction.
Although more costly than either CCD cameras or PDAs, the CMOS camera boasts exciting technology that has recently become capable of delivering extremely high spatiotemporal resolution. RedShirt Imaging lists the CardioCMOS-128f as being capable of recording 128x128 spots at an acquisition rate of 10,000 fps (RedShirt Imaging, LLC, Decatur, GA). CMOS technology typically allows for a very large well depth and a large dynamic range, although these cameras still lack the DC coupling ability of the PDA system. The architecture of the CMOS detector is what sets it apart from the CCD detector. The CMOS detector uses specialized manufacturing techniques to create micro-arrays of photodetectors, each with their own dedicated amplifiers and independent circuitry. Thus CMOS detectors are capable of performing signal processing functions on a pixel-by-pixel basis. The ability of the CMOS detector pixel to record and process signals on a pixel-by-pixel basis is in contrast to the CCD detector, which must transfer signals from individual pixels to be processed by downstream electrical circuitry. This fundamental difference in circuitry architecture is reflected in the higher acquisition rates typically observed in CMOS detectors.
The PDA, as well as the individual photodiode, remains a cost-effective and rugged solution to a wide variety of problems within the field of COM. Photodiode arrays boast a wide spectral response, high dynamic range, high temporal resolution, and the largest well depth of the COM detectors. The PDA is typically a rugged device and can operate in high-light conditions typical in most laboratory experiments, while still delivering a high signal-to-noise ratio. Although low-light conditions can be achieved in COM applications, this typically places limitations on ΔF/F. Another interesting and unique feature of the PDA as a COM detector is that they may be AC coupled. This removal of the DC component of the signal allows the entire dynamic range of the detector to be used on the signal itself, rather than on both the signal and the background baseline fluorescence. This is particularly helpful when imaging with dyes that have a large background fluorescence signal (e.g., di-4-ANEPPS). Due to their robust nature and versatility in a wide variety of applications, photodiode-based systems remain the “workhorse” detector in COM. The PDA operates on principles based on the individual photodiode. Currents from individual photodiodes in the PDA are converted into voltages, which are read out of the detector in parallel and processed by a PC.
The typical noise levels, sensitivity, speed, spatial resolution, and ability to respond to high light levels without saturation (i.e. well depth) are summarized in the following table:
|Noise Levels||Sensitivity||Speed||Spatial Resolution||Well Depth|
For the above chart, more asterisks means better performance, which in the case of COM translates into lower noise, greater sensitivity, greater speed, higher spatial resolution, and a larger well depth. A few caveats apply to the above chart. Although the PMT has extraordinary speed and sensitivity, it is limited by the fact that there is typically no spatial resolution, and must be used in a sequential scanning method. Also, the authors have tried to take into account the
2.3. Filtration, digitization, and multiplexing
Optical signals are subject to several types of noise which must be removed in order to accurately study the details of the cardiac AP. We will briefly review the types of noise most relevant to COM systems. Various types of white noise are ubiquitous throughout all types of electronics, and are typically of frequencies well above those of cardiac signals. Thus low pass filters are typically used to help remove white noise. Sixty-cyle is another type of noise that is often encountered when collecting optical signals. This noise may contaminate signals by way of electromagnetic waves from nearby power outlets, or may be introduced if equipment used in optical mapping experiments is powered using AC power (i.e., if the equipment is not isolated). It may be alleviated by the use of a Faraday cage and/or the use of a band-stop (i.e., notch) filter centered at 60-Hz. Mechanical vibrations may also affect optical signals, and can range from fluctuations in air current to vibrations due to foot traffic. Sources of mechanical vibrations are highly varied in nature, and must be dealt with on a case-by-case basis. A research-grade optical table with active isolation should be sufficient to suppress most sources of mechanical noise. Optical recordings may also contain drifts in basline voltage due to several sources. These include photobleaching, dye washout, and dye internalization into the inner leaflet of the cell membrane. One way to reduce the impact of these noise sources is to employ the technique of ratiometry (Knisley et al., 2000), which is discussed in the
Many optical devices, including photodiodes and photomultiplier tubes, record analog signals that must be digitized before being stored on a PC. Digitization equipment must have sufficient speed and throughput in order to follow the high spatiotemporal resolution required for optical mapping applications. A review article by Entcheva et al. summarizes the state of the art in this sub-field of COM (Entcheva & Bien, 2006).
When data is recorded from a two-dimensional grid of sites simultaneously, the most intuitive storage method is a two-dimensional matrix of values. However, prior to storage on a PC, this data must be routed from the digitization equipment to the PC. This requires arranging the data in a sequential fashion, a process known as multiplexing. For a 16x16 element PDA, data for a single millisecond might be arranged from sites 1 to 256 and then be sequentially sent to the PC for storage. Following this, the data could then be demultiplexed and arranged as a 2D matrix, a more logical form for creating activation maps and other graphics to assist with visualization of the data.
2.4. Pre/post-conditioning of optical data
In contrast to electrical signals, optical signals are highly sensitive to heart motion. Various methods have been developed in order to reduce “motion artifacts” which are often present in optical signals. These motion artifacts are thought to be the result of a change in the location of the mapped region on the heart surface, where a fluorescence gradient typically exists (Himel & Knisley, 2006). Methods to reduce motion artifact include physically restraining the heart (Efimov et al., 1996; Girouard et al., 1996), the use of electromechanical uncouplers (Federov et al., 2007; Jalife et al., 1998; Li & Nattel, 2007; Wu et al., 1998), and the technique of ratiometry (Hooks et al., 2001; Knisley et al., 2000; Kong et al., 2003). Physically restraining the heart reduces motion artifacts simply by limiting the extent to which the heart can move during contraction, thus limiting the amount by which the mapped region moves with respect to its original position on the heart. Electromechanical uncouplers work by a variety of methods, but most have an effect upon the actin-myosin cytoskeleton which is responsible for contraction. Electromechanical uncouplers should be used with care, as some studies have shown that these agents can affect various parameters of the cardiac AP and may also affect the dynamics of ventricular fibrillation (VF) (Baker et al., 2004; Biermann et al., 1998; Hayashi et al., 2003; Lee et al., 2001).
Ratiometry is a signal processing method that requires the use of two recorded optical signals, each of a different band of the wavelength spectrum. We will refer to the longer wavelength signal as “red” and the shorter wavelength signal as “green”. When the Vm-sensitive dye di-4-ANEPPS is excited, the peak of the emission spectrum of the dye shifts toward shorter wavelengths (green). Thus the green signal would show an increase in fluorescence intensity while the red signal would show a corresponding decrease in intensity. Using this method, an upright cardiac AP would be recorded in the green signal while an inverted (or “upside-down”) AP would be recorded in the red signal. The important thing to consider is that the emission signal corresponding to Vm is emitted at a relatively narrow frequency band. Contrastingly, emission due to motion is not heavily wavelength-dependent, and will cause the change in fluorescence signals in the same direction regardless of the wavelength band of the collected signal. Since the motion signals are common to both collected wavelengths, we may reduce motion artifacts by simply taking the ratio of the green signal to the red on a point-by-point basis. This will cause a significant reduction in the motion artifact, and will help us to isolate the electrical signal. This technique could be achieved in the laboratory by using dual PDAs and separating fluorescence emission into two wavelength bands, one above and one below the peak emission wavelength of the dye of interest. In addition to motion artifact removal, ratiometric signals have also been used to study motion artifacts optically. This may be achieved by subtracting the electrical signal from the signal containing both electrical and motion components (Himel et al., 2006).
Since fibrillation (both ventricular and atrial) is a topic of great clinical and theoretical interest, considerable effort has been expended in order to analyze data recorded from the fibrillating heart. This data is challenging to analyze and interpret, since recordings of fibrillation often have a chaotic appearance when viewed with time as the horizontal axis. Thus a variety of alternate methods have been used to gain insight into the nature of fibrillation, including dominant frequency analysis (Caldwell et al., 2007; Choi et al., 2003; Choi et al., 2006; Joel & Hsia, 2005; Moreno et al., 2005; Wu et al., 2004; Wu et al., 2006; Zaitsev et al., 2003) and mutual information (Omichi et al., 2004, Wu et al., 2005). More recently the use of a metric known as spatiotemporal entropy has been used to analyze oscillatory dynamics in cardiac and neural systems (Bub et al., 2005; Jung et al., 2000; Himel et al., 2009).
Dominant frequency analysis involves use of the Fourier transform to examine the frequency content of fibrillation recordings. Some groups have used dominant frequency analysis to support the theory of a “mother rotor” (i.e., a high-frequency region of the heart that drives fibrillatory activity) (Chen et al., 2000; Jalife et al., 1998; Zaitsev et al., 2000), however this issue remains controversial (Berenfeld et al., 2001). Others believe that fibrillation is maintained by the constant creation and annihilation of wavelets which occur due to functional and anatomical heterogeneities (Choi et al., 2002; Lee et al., 1996; Moe, 1962; Rogers et al., 1999; Valderrábano et al., 2002). More recently, some groups have concluded that fibrillatory activity may be driven by both mechanisms, depending on specific conditions in the heart (Chen et al., 2003; Liu et al., 2004; Nash et al., 2006; Wu et al., 2002).
Mutual information techniques have been used to examine the relationship between Vm and Cai. Algorithms assign a numerical value for individual signals in order to quantify the degree of similarity between Vm and Cai during fibrillation in order to give insight into the mechanisms of arrhythmogenesis and the maintenance of fibrillation; however, mutual information must be calculated for individual signal pairs and by itself does not indicate
Spatiotemporal entropy has been used to quantify the degree of uncertainty in both time and space by considering them as lumped parameters, and analyzing activations in the context of space-time cubes (i.e., stacked two-dimensional optical maps with time as the third dimension). Spatiotemporal entropy has been used to analyze neural simulations as well as oscillatory dynamics in cultured cell monolayers (Bub et al., 2005; Jung et al., 2000). Spatiotemporal entropy analysis is appealing when analyzing optical mapping data, since one of the main strengths of COM data lies in its spatial resolution.
3. Studies in cardiac electrophysiology research using photodiode arrays
This section will showcase three recent optical mapping studies from this lab which examine cardiac arrhythmia mechanisms in the context of global ischemia (Himel et al., 2009; Lakireddy et al., 2005; Lakireddy et al., 2006). These studies examine the dynamic relationship between Vm and Cai over the course of ischemia/reperfusion injury. These studies used a photodiode-based system which simultaneously records Vm and Cai with two separate 16x16 PDAs (figure 2 shows 16x16 element photoactive region of the PDA, see figure 3 for a schematic of the simultaneous dual-measurement system). This system was designed by B.R. Choi and G. Salama, and uses two Hamamatsu C4675-103 detectors. Please see the excellent review by Salama et al. for more details regarding this system (Salama et al., 2009). The whole-heart guinea pig (GP) Langendorff model was used, where the aorta was cannulated and hung vertically within a surrounding bath of physiological solution referred to as Tyrode’s solution. Heated (~37° C), oxygenated Tyrode’s solution was pumped retrogradely through the aorta in order to nourish and oxygenate the heart via the coronary arteries. These studies simulated the condition of global ischemia (i.e., oxygen deprivation) and reperfusion by temporarily interrupting perfusion into the coronary arteries and then restarting the perfusion. The guinea pig animal model is similar to human physiology in terms of its AP morphology and intracellular calcium handling.
Study 1 (Himel et al., 2009): Ventricular tachycardia is frequently observed in the clinical environment. VT may either be brief (<30 seconds) or may continue indefinitely, and patient survival is critically dependent on the duration of the VT episode. Spontaneous termination of the VT episode within a few seconds leaves the patient unharmed. If the VT episode continues however, the patient is typically dependent upon a defibrillatory shock to terminate the episode, without which the patient may not survive. Although extremely important clinically, the mechanism which determines the duration of the VT episode is not well understood. It is thought that uncoupling between Vm and Cai coupling may play a role in determining whether an episode of VT will terminate spontaneously or continue indefinitely. This study sought to use the concept of spatiotemporal entropy (E) as a metric to determine the degree of uncoupling between Vm and Cai during the early phases of a VT episode, and to test the predictive power of E for VT duration. During normal sinus rhythm, propagation of excitation waves in the heart is uniform and wavebreak is typically not observed. During episodes of VT, wavebreak is common, and fragmentation of excitation waves in the heart is frequently observed. Moreover, differences between Vm and Cai signals are often observed. These differences in Vm/Cai coupling may be quantified using the absolute value of the difference in E, symbolized by Ed.
This study examined several groups of VT episodes which were divided according to whether or not they terminated spontaneously. Self-terminating episodes of VT were further classified as short (<5 seconds) or long (>5 seconds). Ed was determined for the first 500 ms of all VT episodes. Ed values for non self-terminating episodes of VT were significantly greater than self-terminating VT episodes. Further, Ed values for long self-terminating episodes of VT were significantly greater than those for short self-terminating episodes.
Wavebreak and differences in Vm and Cai spatial fluorescence maps were consistently identified during periods of high E and high Ed. The results are illustrated by figures 4 and 5, which show typical examples of low and high Ed during VT episodes.
High Ed correlated with a greater duration of a VT episode. This may be related to destabilization of propagation and uncoupling between Vm and Cai activation wavefronts.
Study 2 (Lakireddy et al., 2005): This study examines spatial dispersion of repolarization in the context of global ischemia, and also the role spatial dispersion plays in the development of electrical alternans. Electrical alternans is a term used to describe beat-to-beat alterations in AP morphology. For example, a one-to-one APD alternans occurs when a normal AP is followed by a short-duration AP, which is then followed by normal AP, short-duration AP, and so on. Electrical alternans are considered to be a strong marker of electrical instability, and often precede malignant arrhythmias such as ventricular tachycardia (VT) and VF (Hohnsloser et al., 1998; Gold et al., 2000; Ikeda et al., 2006; Pastore et al., 1999; Pham et al., 2003; Rashba et al., 2004; Rosenbaum et al., 1994).
In this study by Lakireddy et al., ischemia-induced changes in APD and intracellular calcium transient duration (CaiT-D) were determined, and their relationship with electrical alternans was investigated. Recordings show that ischemia resulted in a significant decrease in APD, but resulted in a significant increase in CaiT-D. In addition, changes in APD were spatially heterogeneous while changes in CaiT-D were relatively homogeneous (see figure 6). Sites with less shortening of APD displayed alternans in both CaiT-D and APD, while sites with more shortening of APD displayed CaiT-D alternans but little or no APD alternans, leading to a condition of significant spatial dispersion of the APD. The condition of increased spatial dispersion due to ischemia is thought to account for the vulnerability of the heart to alternans.
Study 3 (Lakireddy et al., 2006): In a second study by Lakireddy et al., the association between arrhythmogenesis and spontaneous calcium oscillations (S-CaOs) was examined in the intact heart. It is known that ischemia/reperfusion leads to elevated Cai and an alteration in Cai kinetics (Bers, 2002; Marban et al., 1990; Steenbergen et al., 1987). This alteration in normal Cai kinetics can lead to S-CaOs. Under such conditions, the normal master/slave relationship between Vm/CaiT signals is reversed (i.e., calcium signals precede and drive Vm).
The goal of this study was to investigate the correlation between S-CaOs and arrhthmogenesis using an experimental GP model with 15 minutes of no flow ischemia followed by 15 minutes of reperfusion. Changes in Cai and Vm in a limited zone of the epicardial surface of the GP heart were simultaneously recorded and carefully examined. The study provided evidence of a linkage between S-CaOs and arrhythmogenesis in the setting of ischemia/reperfusion (I/R). In the intact heart during I/R, spontaneous premature beats (PBs) occurred and were ubiquitous. Some PBs initiated a VT or VF (see figure 7), while others remained confined to their site of origin and did not result in an arrhythmia (see figure 8). Two important observations had to be made in order to link an arrhythmia to S-CaOs in the experimental model. First, the beginning of S-CaOs preceded the onset of the simultaneously recorded membrane depolarization by 2-15 ms at a very restricted site in the optical field. In recordings obtained further away from the focal site of origin, the relative amplitude of the S-CaOs gradually decreased and the start of membrane depolarization preceded the onset of S-CaOs. Second, the presence of some degree of conduction block, which by definition is the failure of S-CaOs to trigger a fully propagated response, was essential for the localization of the focal site of origin. Thus S-CaOs may remain concealed (and hence benign) or may manifest as PBs, VT or VF.
4. Recent advances in cardiac optical mapping
We will now turn to briefly discuss a few of the important recent advances in COM. Optrodes are bundles of microscopic fiberoptic cables which are inserted into cardiac tissue in order to interrogate intramural activation patterns. They are similar to plunge needle electrodes in their usage; however optrodes are capable of measuring complete APs, including the repolarization phase, whereas plunge electrodes measure extracellular potentials only. It is thought that in the future optrodes will play an important role in more carefully examining transmural dispersion of repolarization, an important factor in arrhythmogenesis in a variety of cardiac diseases (Antzelevitch, 2007; El-Sherif et al., 1996; Milberg et al., 2005; Shimizu et al., 1997). Optrodes may also be important in the study of the dynamics of arrhythmic circuits, since they are often present deeper in the myocardial wall (Allison et al., 2007; Li et al., 2008; Valderrábano et al., 2001). However, like plunge electrodes optrodes must also be inserted into the tissue and therefore cause damage which may by itself alter activation patterns. Thus the effects of this insertion must be carefully
Several groups have recently begun to use multiple cameras to simultaneously interrogate opposing sides of the ventricular wall (Evertson et al., 2008; Kay & Rogers, 2006; Kay et al., 2004; Kay et al., 2006; Rogers et al., 2007). In addition, some of these groups use additional cameras to recreate the geometry of the heart in order to properly orient optical maps from several cameras on the epicardial surface (Kay et al., 2004; Evertson et al., 2008). Most Panoramic optical mapping systems are based on CCD technology, however systems have also been built using multiple PDAs (Qu et al., 2007). Panoramic optical mapping does not address the problem of lost depth information, but does provide a significant improvement over traditional optical mapping which only maps a limited region on the epicardial surface.
The use of monolayer cell cultures in COM also represents an important advance, allowing for highly controlled studies of basic conduction as well as studies to elucidate fundamental arrhythmic mechanisms (Bub et al., 1998; Entcheva et al., 2000; Fast et al., 2000; Iravanian et al., 2003; Tung & Cysyk, 2007). An appealing aspect of the cardiac monolayer is that it allows us to study conduction in cardiac tissue without the complexity associated with the three-dimensional whole-heart Langendorff model. Since the cardiac monolayer is essentially two-dimensional (only tens of micrometers thick while being tens of millimeters in diameter), the entire monolayer may be mapped; therefore data interpretation is not complicated by the absence of missing depth information. And although the monolayer is technically three-dimensional, typical optical mapping systems interrogate at sufficient depths so that no information is lost beneath the surface (Ding et al., 2001). Despite being similar to whole-heart mapping in many respects, the actual practice of monolayer mapping carries with it significant challenges, and is in many respects more difficult than whole-heart mapping (Entcheva & Bien, 2006).
Photodiodes have played an essential role in the development of the field of COM. They were used in the earliest COM systems and continue to have widespread use today, both in typical applications as well as more modern designs such as optrodes and panoramic systems. Applications for photodiodes within COM continue to emerge, and will likely remain a vital part of this important and ever-expanding branch of cardiac electrophysiology research.
6. List of abbreviations
AP – action potential
AP-A – anthopleurin-A
APD – action potential duration
Cai – intracellular calcium
CaiT – intracellular calcium transient
CaiT-D – intracellular calcium transient duration
CCD - charge-coupled device
CL – cycle length
CMOS - complimentary metal-oxide semiconductor
COM – cardiac optical mapping
GP – guinea pig
I/R – ischemia/reperfusion
LQTS – long QT syndrome
LQT3 – long QT syndrome 3
PB – premature beat
PDA - photodiode array
PMT - photomultiplier tube
TdP – Torsades de Pointes
VF – ventricular fibrillation
Vm – transmembrane voltage
VT – ventricular tachycardia
Allison J. S. Qin H. Dosdall D. J. Huang J. Newton J. C. JD Allred Smith. W. M. Ideker R. E. 2007The transmural activation sequence in porcine and canine left ventricle is markedly different during long-duration ventricular fibrillation. J Cardiovasc Electrophysiol; 18 1306 1312
Al-Khadra A. Nikolski V. Efimov I. R. 2000The role of electroporation in defibrillation. Circ Res; 87 797 804
Antzelevitch C. 2007 Role of spatial dispersion of repolarization in inherited and acquired sudden cardiac death syndromes. Am J Physiol Heart Circ Physiol;293:H2024-H2038 2024
Baker L. C. Wolk R. Choi B. R. Watkins S. Plan P. Shah A. Salama G. 2004Effects of mechanical uncouplers, diacetyl monoxime, and cytochalasin-D on the electrophysiology of perfused mouse hearts. J Physiol Heart Circ Physiol;287:H 1771H1779.
Baxter W. T. Mirinov S. F. Zaitsev A. V. Jalife J. Pertsov A. M. 2001Visualizing excitation waves inside cardiac muscle using transillumination. Biophys J; 80 516 530
Berenfeld O. Pertsov A. M. Jalife J. 2001Letter to the editor: What is the organization of waves in ventricular fibrillation. Circ Res;89:e22.
Bers DM. 2002Calcium and cardiac rhythms: physiological and pathophysiological. Circ Res; 90 14 17
Biermann M. Rubart M. Moreno A. Wu J. Josiah-Durant A. Zipes D. P. 1998Differential effects of cytochalasin D and 2,3 butanedione monoxime on isometric twitch force and transmembrane action potential in isolated ventricular muscle: implications for optical measurements of cardiac repolarization. J Cardiovasc Electrophysiol; 9 1348 1357
Bub G. Glass L. Publicover N. G. Shrier A. 1998Bursting calcium rotors in cultured cardiac myocyte monolayers. Proc Natl Acad Sci USA; 95 10283 10287
Bub G. Shrier A. Glass L. 2005Global organization of dynamics in oscillatory heterogeneous excitable media. Physical Review Lett;94:028105.
Byars JL, Smith WM, Ideker RE, Fast VG. 2003Development of an optrode for intramural multisite optical recordings of Vm in the heart. J Cardiovasc Electrophysiol; 14 1196 1202
Cabo C. Pertsov A. M. Baxter W. T. Davidenko J. M. Gray R. A. Jalife J. 1994Wave-front curvature as a cause of slow conduction and block in isolated cardiac muscle. Circ Res; 75 1014 1028
Caldwell J. Burton F. L. Smith G. L. Cobbe S. T. 2007Heterogeneity of ventricular fibrillation dominant frequency during global ischemia in isolated rabbit hearts. J Cardiovasc Electrophysiol; 18 854 861
Caldwell B. J. Legrice I. J. Hooks D. A. Tai D. Pullan A. J. Smaill B. H. 2005Intramural measurement of transmembrane potential in the isolated pig heart: Validation of a novel technique. J Cardiovasc Electrophysiol; 16 1001 1010
Chen J. Mandapati R. Berenfeld O. Skanes A. C. Jalife J. 2000High-frequency periodic sources underlie ventricular fibrillation in the isolated rabbit heart. Circ Res; 86 86 93
Chen PS, Wu TJ, Ting CT, Karagueuzian HS, Garfinkel A, Lin SF, Weiss JN. 2003 A tale of two fibrillations. Circulation;108:2298-2303.
Cheng Y. 2006Optical mapping of shock-induced arrhythmogenesis in the rabbit heart with healed myocardial infarction: fluorescent imaging with a photodiode array. Methods Mol Med; 129 149 61
Choi B. R. Burton F. Salama G. 2002Cytosolic Ca2+ triggers early afterdepolarizations and torsade de pointes in rabbit hearts with type 2 long QT syndrome. J Physiol; 543 615 631
Choi B. R. Hatton W. J. Hume J. R. Liu T. Salama G. 2006Low osmolarity transforms ventricular fibrillation from complex to highly organized, with a dominant high-frequency source. Heart Rhythm; 3 1210 1220
Choi B. R. Jang W. Salama G. 2007Spatially discordant voltage alternans cause wavebreaks in ventricular fibrillation. Heart Rhythm; 4 1057 1068
Choi B. R. Liu T. Lavasani M. Salama G. 2003Fiber orientation and cell-cell coupling influence ventricular fibrillation dynamics. J Cardiovasc Electrophysiol; 14 851 860
Choi B. R. Salama G. 2000Simultaneous maps of optical action potentials and calcium transients in guineapig hearts: mechanisms underlying concordant alternans. J Physiol; 529 171 188
Coraboeuf E. Weidmann S. 1949Potentiels d’action du muscle obtenus à l’aide de microélectrodes intracellulaires. Présence d’une inversion de potential. CR Soc Biol (Paris); 143 1360 1360
Coraboeuf E. Weidmann S. 1949Potentiels de repos et potentials d’action du muscle cardiaque, mesurés á l’aide d’électrodes intracellulaires. CR Soc Biol (Paris); 143 1329 1331
Draper M. H. Weidmann S. 1951Cardiac resting and action potentials recorded with an intracellular electrode. J Physiol; 115 74 94
Ding L. Splinter R. Knisley S. B. 2001Quantifying spatial localization of optical mapping using Monte Carlo simulations. IEEE Trans Biomed Eng; 48 1098 1107
Efimov I. R. Ermentrout B. Huang D. T. Salama G. 1996Activation and repolarization patterns are governed by different structural characteristics of ventricular myocardium: Experimental study with voltage-sensitive dyes and numerical simulations. J Cardiovasc Electrophysiol; 7 512 530
El -Sherif N. 2007The challenge of cardiac tridimensional mapping. Heart Rhythm; 4 1437 1440
El -Sherif N. Caref E. B. Yin H. Restivo M. 1996The electrophysiological mechanism of ventricular arrhythmias in the long QT syndrome. Circ Res; 79 474 492
Entcheva E. Bien H. 2006Macroscopic optical mapping of excitation in cardiac cell networks with ultra-high spatiotemporal resolution. Progress in Biophysics and Molecular Biology; 92 232 257
Entcheva E. Lu S. N. Troppman R. H. Sharma V. Tung L. 2000Contact fluorescence imaging of reentry in monolayers of cultured neonatal rat ventricular myocytes. J Cardiovasc Electrophysiol; 11 665 676
Evertson D. W. Holcomb M. R. Eames M. D. C. MA Bray Sidorov. V. Y. Xu J. Wingard H. Dobrovolny H. M. Woods M. C. Gauthier D. J. Wikswo J. P. 2008High-resolution high-speed panoramic cardiac imaging system. IEEE Trans Biomed Eng; 55 1241 1243
Fast VG. 2005Simultaneous optical imaging of membrane potential and intracellular calcium. J Electrocardiol; 38 107 112
Fast VG, Ideker RE. 2000Simultaneous optical mapping of transmembrane potential and intracellular calcium in myocyte cultures. J Cardiovasc Electrophysiol; 11 547 556
Fast VG, Sharifov OF, Cheek ER, Newton JC, Ideker RE. 2002Intramural virtual electrodes during defibrillation shocks in left ventricular wall assessed by optical mapping of membrane potential. Circulation; 106 1007 1014
Fedorov V. V. Kostecki G. Hemphill M. Efimov I. R. 2008Atria are more susceptible to electroporation than ventricles: Implications for atrial stunning, shock-induced arrhythmia and defibrillation failure. Heart Rhythm; 5 593 604
Fedorov V. V. Lozinsky I. T. Sosunov E. A. Anyukhovsky E. P. Rosen M. R. Balke W. Efimov I. R. 2007Application of blebbistatin as an excitation- contraction uncoupler for electrophysiologic study of rat and rabbit hearts. Heart Rhythm; 4 619 626
Girouard SD, Pastore JM, Laurita KR, Gregory KW, Rosenbaum DS. 1996Optical mapping in a new guinea pig model of ventricular tachycardia reveals mechanisms for multiple wavelengths in a single reentrant circuit. Circulation; 93 603 613
Gold MR, Bloomfield DM, MD, Anderson KP, El-Sherif N, Wilber DJ, Groh WJ, Estes M, Kaufman ES, Greenberg ML, Rosenbaum DS. 2000 A comparison of T-wave alternans, signal averaged electrocardiography and programmed ventricular stimulation for arrhythmia risk stratification. JACC;36:2247-2253.
Hayashi H. Miyauchi Y. Chou C. C. Karagueuzian H. S. Chen P. S. Lin S. F. 2003Effects of cytochalasin D on electrical restitution and the dynamics of ventricular fibrillation in isolated rabbit heart. J Cardiovasc Electrophysiol; 14 1077 1084
Himel I. V. H. D. Bub G. Yue Y. El -Sherif N. 2009Early voltage/calcium uncoupling predestinates the duration of ventricular tachyarrhythmias during ischemia/reperfusion. Heart Rhythm; 6 1359 1365
Himel IV, HD, Dumas JH, Kiser AC, Knisley SB. 2007Translesion stimulus-excitation delay indicates quality of linear lesions produced by radiofrequency ablation in rabbit hearts. Physiol. Meas.; 28 611 623
Himel I. V. H. D. Knisley S. B. 2006Imaging of cardiac movement using ratiometric and nonratiometric optical mapping: effects of ischemia and 2, 3-butaneodione monoxime. IEEE Trans Med Imaging; 25 122 127
Hohnloser S. H. Klingenheben T. Li Y. G. Zabel M. Peetermans J. Cohen R. J. 1998T wave alternans as a predictor of recurrent ventricular tachyarrhythmias in ICD recipients: Prospective comparison with conventional risk markers. J Cardiovasc Electrophysiol; 9 1258 1268
Hooks DA, LeGrice IJ, Harvey JD, Smaill BH. 2001Intramural multisite recording of transmembrane potential in the heart. Biophys J; 81 2671 2680
Hyatt CJ, Mirinov SF, Vetter FJ, Zemlin CW, Pertsov AM. 2005Optical Action Potential Upstroke Morphology Reveals Near-Surface Transmural Propagation Direction. Circ Res; 97 277 284
Ikeda T, Yoshino H, Sugi K, Tanno K, Shimizu H, Watanabe J, Kasamaki Y, Yoshida A, Kato T. 2006 Predictive value of microvolt T-wave alternans for sudden cardiac death in patients with preserved cardiac function after acute myocardial infarction. JACC;48:2268-2274.
Iravanian S. Nabutovsky Y. Kong C. R. Saha S. Bursac N. Tung L. 2003Functional reentry in cultured monolayers of neonatal rat cardiac cells. Am J Physiol Heart Circ Physiol;285:H 449H456.
Jalife J. Berenfeld O. Skanes A. Mandapati R. 1998Mechanisms of atrial fibrillation: mother rotors or multiple daughter wavelets, or both? J Cardiovasc Electrophysiol;9:S 2 12
Joel SE, Hsia PW. 2005Discovery of gradient pattern in dominant frequency maps during fibrillation: implication of rotor drift and epicardial conduction velocity changes. J Electrocardiol; 38 159 165
Jung P, Wang J, Wackerbauer R, Showalter K. 2000 Coherent structure analysis of spatiotemporal chaos. Physical Review E; 61 2095 2098
Kay MW, Amison PM, Rogers JM. 2004Three-dimensional surface reconstruction and panoramic optical mapping of large hearts. IEEE Trans Biomed Eng; 51 1219 1229
Kay MW, Rogers JM. 2006 Epicardial rotors in panoramic optical maps of fibrillating swine ventricles. Conf Proc IEEE Med Biol Soc; 1 2268 2271
Kay MW, Walcott GP, Gladden JD, Melnick SB, Rogers JM. 2006 Lifetimes of epicardial rotors in panoramic optical maps of fibrillating swine ventricles. Am J Physiol Heart Circ Physiol; 291 1935 1941
Knisley S. B. Baynham T. 1997Line stimulation parallel to myofibers enhances regional uniformity of transmembrane voltage changes in rabbit hearts. Circ Res; 81 229 241
Knisley SB, Hill BC. 1995Effects of bipolar point and line stimulation in anisotropic rabbit epicardium: assessment of the critical radius of curvature for longitudinal block. IEEE Trans Biomed Eng; 42 957 966
Knisley SB, Hill BC, Ideker RE. 1994Virtual electrode effects in myocardial fibers. Biophys J; 66 719 728
Knisley S. B. Justice R. K. Kong W. Johnson P. L. 2000Ratiometry of transmembrane voltage-sensitive fluorescent dye emission in hearts. Am J Physiol Heart Circ Physiol; 279: H 1421H1433.
Knisley S. B. Trayanova N. Aguel F. 1999Roles of electric field and fiber structure in cardiac electric stimulation. Biophys J; 77 1404 1417
Kong W. Fakhari N. Sharifov O. F. Ideker R. E. Smith W. M. Fast V. G. 2007Optical measurements of intramural action potentials in isolated porcine hearts using optrodes. Heart Rhythm; 4 1430 1436
Kong W. Walcott G. P. Smith W. M. Johnson P. L. Knisley S. B. 2003Emission ratiometry for simultaneous calcium and action potential measurements with coloaded dyes in rabbit hearts: reduction of motion and drift. J Cardiovasc Electrophysiol; 14 76 82
Lakireddy V. Baweja P. Syed A. Bub G. Boutjdir M. El -Sherif N. 2005Contrasting effects of ischemia on the kinetics of membrane voltage and intracellular calcium transient underlie electrical alternans. Am J Physiol Heart Circ Physiol; 288 400 407
Lakireddy V. Bub G. Baweja P. Syed A. Boutjdir M. El -Sherif N. 2006The kinetics of spontaneous calcium oscillations and arrhythmogenesis in the in vivo heart during ischemia/reperfusion. Heart Rhythm; 3 58 66
Lan DZ, Pollard AE, Knisley SB. 2007Optical mapping of Vm and Cai 2+ in a model of arrhythmias induced by local catecholamine application in patterned cell cultures. Eur J Physiol; 453 871 877
Laurita KR, Singal A. 2001 Mapping action potentials and calcium transients simultaneously from the intact heart. Am J Physiol Heart Circ Physiol;280: H2053-H2060 2053
Lee J. J. Kamjoo K. Hough D. Hwang C. Fan W. Fishbein M. C. Bonometti C. Ikeda T. Karagueuzian H. S. Chen P. S. 1996Reentrant wave fronts in wiggers’ stage II ventricular fibrillation. Circ Res; 78 660 675
Lee M. H. Lin S. F. Ohara T. Omichi C. Okuyama Y. Chudin E. Garfinkel A. Weiss J. N. Karagueuzian H. S. Chen P. S. 2001Effects of diacetyl monoxime and cytochalasin D on ventricular fibrillation in swine right ventricles. Am J Physiol Heart Circ Physiol;280: H 2689H2696.
Li D. Nattel S. 2007Pharmacological elimination of motion artifacts during optical imaging of cardiac tissues: Is blebbistatin the answer? Heart Rhythm; 4 627 628
Li L. Jin Q. Huang J. Cheng K. A. Ideker R. E. 2008Intramural foci during long duration fibrillation in the pig ventricle. Circ Res; 102 1256 1264
Liu Y. B. Pak H. N. Lamp S. T. Okuyama Y. Hayashi H. Wu T. J. Weiss J. N. Chen P. S. Lin S. F. 2004Coexistence of two types of ventricular fibrillation during acute regional ischemia in rabbit ventricle. J Cardiovasc Electrophysiol; 15 1433 1440
Marban E. Kitakaze M. Koretsune Y. Yue D. T. Chacko V. P. MM Pike 1990Quantification of [Ca2+]i in perfused hearts: critical evaluation of the 5F-BAPTA and nuclear magnetic resonance method as applied to the study of ischemia and reperfusion. Circ Res; 66 1255 1267
Matiukas A. Mitrea B. G. Pertsov A. M. Wuskell J. P. Wei M. Watras J. Millard A. C. Loew L. M. 2006New near-infrared optical probes of cardiac electrical activity. Am J Physiol Heart Circ Physiol;290: H 2633H2643.
Matiukas A. Mitrea B. G. Qin M. Pertsov A. M. Shvedko A. G. MD Warren Zaitsev. A. V. Wuskell J. P. Wei M. Watras J. Loew L. M. 2007Near-infrared voltage-sensitive fluorescent dyes optimized for optical mapping in blood-perfused myocardium. Heart Rhythm; 4 1441 1451
Milberg P. Reinsch P. Wasmer K. Mönnig G. Stypmann P. Osada N. Breithardt G. Haverkamp W. Eckardt L. 2005Transmural dispersion of repolarization as a key factor of arrhythmogenicity in a novel intact heart model of LQT3. Circ Res; 65 397 404
Moe GK. 1962On the multiple wavelet hypothesis of atrial fibrillation. Arch Int Pharmacodyn Ther; 140 183 188
Morad M. Dillon S. 1981A new laser scanning system for measuring action potential propagation in the heart. Science; 214 453 456
Morad M. Salama G. 1979Optical probes of membrane potential in heart muscle. J Physiol; 292 267 295
Moreno J. Zaitsev A. V. Warren M. Berenfeld O. Kalifa J. Lucca E. Mironov S. Guha P. Jalife J. 2005Effect of remodelling, stretch and ischaemia on ventricular fibrillation frequency and dynamics in a heart failure model. Cardiovasc Res; 65 158 166
Nash M. P. Mourad A. Clayton R. H. Sutton P. M. CP Bradley Hayward. M. Paterson D. J. Taggart P. 2006Evidence for multiple mechanisms in human ventricular fibrillation. Circulation; 114 536 542
Neunlist M, Tung L. 1995 Spatial distribution of cardiac transmembrane potentials around an extracellular electrode: Dependence on fiber orientation. Biophys J; 68 2310 2322
Omichi C, Lamp ST, Lin SF, Yang J, Baher A, Zhou S, Attin M, Lee MH, Karagueuzian HS, Kogan B, Qu Z, Garfinkel A, Chen PS, Weiss JN. 2004 Intracellular Ca dynamics in ventricular fibrillation. Am J Physiol Heart Circ Physiol;286: H1836-H1844
Pastore JM, Girouard SD, Laurita KR, Akar FG, Rosenbaum DS. 1999Mechanism linking T-wave alternans to the genesis of cardiac fibrillation. Circulation; 99 1385 1394
Perez FJ, Wood MA, Schubert CM. 2006Effects of gap geometry on conduction through discontinuous radiofrequency lesions. Circulation; 113 1723 1729
Pham Q. Quan K. J. DS Rosenbaum 20. T-wave alternans. marker mechanism. methodology for. predicting sudden. cardiac death. 2003 Twave alternans: marker, mechanism, and methodology for predicting sudden cardiac death. J Electrocardiol; 36 75 81
Pruvot EJ, Katra RP, Rosenbaum DS, Laurita KR. 2004Role of calcium cycling versus restitution in the mechanism of repolarization alternans. Circ Res; 94 1083 1090
Qu F. Ripplinger C. M. Nikolski V. P. Grimm C. Efimov I. R. 2007Three-dimensional panoramic imaging of cardiac arrhythmias in rabbit heart. J Biomed Opt.;12:044019.
Rashba E. J. Osman A. F. Macmurdy K. MM Kirk Sarang. S. E. Peters R. W. Shorofsky S. R. Gold M. R. 2004Enhanced detection of arrhythmia vulnerability using T wave alternans, left ventricular ejection fraction, and programmed ventricular stimulation: A prospective study in subjects with chronic ischemic heart disease. J Cardiovasc Electrophysiol; 15 170 176
Rogers J. M. Huang J. Smith W. M. Ideker R. E. 1999Incidence, evolution, and spatial distribution of functional reentry during ventricular fibrillation in pigs. Circ Res; 84 945 954
Rogers JM, Walcott GP, Gladden JD, Melnick SB, Kay MW. 2007Panoramic optical mapping reveals continuous epicardial reentry during ventricular fibrillation in the isolated swine heart. Biophys J; 92 1090 1095
DS Rosenbaum Jackson. L. E. Smith J. M. Garan H. Ruskin J. N. Cohen R. J. 1994Electrical alternans and vulnerability to ventricular arrhythmias. N Engl J Med; 330 235 241
Sakai T. 2008Optical mapping analysis of the spatiotemporal pattern of experimental tachyarrhythmia in improved isolated rat atrium preparation. J Physiol Sci; 58 87 97
Salama G. 1976Merocyanine 540 as an optical probe of transmembrane electrical activity in the heart. Science; 191 485 487
Salama G. Choi B. R. Azour G. Lavasani M. Tumbev V. BM Salzberg Patrick. MJ Ernst L. A. AS Waggoner 2005Properties of new, long-wavelength, voltage-sensitive dyes in the heart. J. Membrane Biol; 208 125 140
Salama G. Hwang S. M. 2009Simultaneous optical mapping of intracellular free calcium and action potentials from langendorff perfused hearts. Curr Protoc Cytom; Chapter 12: Unit 12.17.
Sano T. Ohtsuka E. Shimamoto T. 1960Unidirectional” atrioventricular conduction studied by microelectrodes. Circ Res; 8 600 608
Sano T. Takayama N. Shimamoto T. 1959Directional difference of conduction velocity in the cardiac ventricular syncytium studied by microelectrodes. Circ Res; 7 262 267
Scherz P. 2007Practical Electronics for inventors. Second edition. New York, NY: McGraw-Hill;.
Shabetai R. Surawicz B. Hammill W. 1968Monophasic Action Potentials in Man. Circulation; 38 341 352
Shimizu W, Antzelevitch C. 1997 Sodium channel block with mexiletine is effective in reducing dispersion of repolarization and preventing torsade de pointes in LQT2 and LQT3 models of the long-QT syndrome. Circulation;96:2038-2047.
Steenbergen C. Murphy E. Levy L. London R. E. 1987Elevation in cytosolic free calcium concentration early in myocardial ischemia in perfused rat heart. Circ Res; 60 700 707
Tung L. Cysyk J. 2007Imaging fibrillation: defibrillation in a dish. J Electrocardiol;40:S 62S65.
Tung L. Zhang Y. 2006Optical imaging of arrhythmias in tissue culture. J Electrocardiol;39:S 2S6.
Valderrábano M. Lee M. H. Ohara T. Lai A. C. Fishbein M. C. Lin S. F. Karagueuzian H. S. Chen P. S. 2001Dynamics of intramural and transmural reentry during ventricular fibrillation in isolated swine ventricles. Circ Res; 88 839 848
Valderrábano M. Yang J. Omichi C. Kil J. Lamp S. T. Qu Z. Lin S. F. Karagueuzian H. S. Garfinkel A. Chen P. S. Weiss J. N. 2002Frequency analysis of ventricular fibrillation in swine ventricles. Circ Res; 90 213 222
Weidmann S. 1951Effect of current flow on the membrane potential of cardiac muscle. J Physiol; 115 227 236
Wu J. Biermann M. Rubart M. Zipes D. P. 1998Cytochalasin D as excitation-contraction uncoupler for optically mapping action potentials in wedges of ventricular myocardium. J Cardiovasc Electrophysiol; 9 1336 1347
Wu S. Weiss J. N. Chou C. C. Attin M. Hayashi H. Lin S. F. 2005Dissociation of membrane potential and intracellular calcium during ventricular fibrillation. J Cardiovasc Electrophysiol; 16 186 192
Wu TJ, Lin SF, Baher A, Qu Z, Garfinkel A, Weiss JN, Ting CT, Chen PS. 2004 Mother rotors and the mechanisms of D600-induced type 2 ventricular fibrillation. Circulation;110:2110-2118.
Wu TJ, Lin SF, Hsieh YC, Ting CT, Chen PS. 2006Ventricular fibrillation during no-flow global ischemia in isolated rabbit hearts. J Cardiovasc Electrophysiol; 17 1112 1120
Zaitsev A. V. Guha P. K. Sarmast F. Kolli A. Berenfeld O. Pertsov A. M. de Groot J. R. Coronel R. Jalife J. 2003Wavebreak formation during ventricular fibrillation in the isolated, regionally ischemic pig heart. Circ Res; 92 546 553
Wu TJ, Lin SF, Weiss JN, Ting CT and Chen PS. 2002 Two types of ventricular fibrillation in isolated rabbit hearts: Importance of excitability and action potential duration restitution. Circulation;106:1859-1866.
Zaitsev A. V. Berenfeld O. Mirinov S. F. Jalife J. Pertsov A. M. 2000Distribution of excitation frequencies on the epicardial and endocardial surfaces of fibrillating ventricular wall of the sheep heart. Circ Res; 86 408 417
Zemlin C. W. Bernus O. Matiukas A. Hyatt C. J. Pertsov A. M. 2008Extracting intramural wavefront orientation from optical upstroke shapes in whole hearts. Biophys J; 95 942 950