Fiber Organization has Little Effect on Electrical Activation Patterns during Focal Arrhythmias in the Left Atrium.

Over the past two decades there has been a steady trend towards the development of realistic models of cardiac conduction with increasing levels of detail. However, making models more realistic complicates their personalization and use in clinical practice due to limited availability of tissue and cellular scale data. One such limitation is obtaining information about myocardial fiber organization in the clinical setting. In this study, we investigated a chimeric model of the left atrium utilizing clinically derived patient-specific atrial geometry and a realistic, yet foreign for a given patient fiber organization. We discovered that even significant variability of fiber organization had a relatively small effect on the spatio-temporal activation pattern during regular pacing. For a given pacing site, the activation maps were very similar across all fiber organizations tested.

Fiber Organization Has Little Effect on Electrical Activation Patterns During Focal Arrhythmias in the Left Atrium.

Over the past two decades there has been a steady trend towards the development of realistic models of cardiac conduction with increasing levels of detail. However, making models more realistic complicates their personalization and use in clinical practice due to limited availability of tissue and cellular scale data. One such limitation is obtaining information about myocardial fiber organization in the clinical setting. In this study, we investigated a chimeric model of the left atrium utilizing clinically derived patient-specific atrial geometry and a realistic, yet foreign for a given patient fiber organization. We discovered that even significant variability of fiber organization had a relatively small effect on the spatio-temporal activation pattern during regular pacing. For a given pacing site, the activation maps were very similar across all fiber organizations tested.

Atrial fibrillation source area probability mapping using electrogram patterns of multipole catheters.

BACKGROUND: Catheter ablation therapy involving isolation of pulmonary veins (PVs) from the left atrium is performed to terminate atrial fibrillation (AF). Unfortunately, standalone PV isolation procedure has shown to be a suboptimal success with AF continuation or recurrence. One reason, especially in patients with persistent or high-burden paroxysmal AF, is known to be due to the formation of repeating-pattern AF sources with a meandering core inside the atria. However, there is a need for accurate mapping and localization of these sources during catheter ablation. METHODS: A novel AF source area probability (ASAP) mapping algorithm was developed and evaluated in 2D and 3D atrial simulated tissues with various arrhythmia scenarios and a retrospective study with three cases of clinical human AF. The ASAP mapping analyzes the electrograms collected from a multipole diagnostic catheter that is commonly used during catheter ablation procedure to intelligently sample the atria and delineate the trajectory path of a meandering repeating-pattern AF source. ASAP starts by placing the diagnostic catheter at an arbitrary location in the atria. It analyzes the recorded bipolar electrograms to build an ASAP map over the atrium anatomy and suggests an optimal location for the subsequent catheter location. ASAP then determines from the constructed ASAP map if an AF source has been delineated. If so, the catheter navigation is stopped and the algorithm provides the area of the AF source. Otherwise, the catheter is navigated to the suggested location, and the process is continued until an AF-source area is delineated. RESULTS: ASAP delineated the AF source in over 95% of the simulated human AF cases within less than eight catheter placements regardless of the initial catheter placement. The success of ASAP in the clinical AF was confirmed by the ablation outcomes and the electrogram patterns at the delineated area. CONCLUSION: Our analysis indicates the potential of the ASAP mapping to provide accurate information about the area of the meandering repeating-pattern AF sources as AF ablation targets for effective AF termination. Our algorithm could improve the success of AF catheter ablation therapy by locating and subsequently targeting patient-specific and repeating-pattern AF sources inside the atria.

Locating Atrial Fibrillation Rotor and Focal Sources Using Iterative Navigation of Multipole Diagnostic Catheters.

PURPOSE: Multi-polar diagnostic catheters are used to construct the 3D electro-anatomic mapping of the atrium during atrial fibrillation (AF) ablation procedures; however, it remains unclear how to use the electrograms recorded by these catheters to locate AF-driving sites known as focal and rotor source types. The purpose of this study is to present the first algorithm to iteratively navigate a circular multi-polar catheter to locate AF focal and rotor sources without the need to map the entire atria. METHODS: Starting from an initial location, the algorithm, which was blinded to the location and type of the AF source, iteratively advanced a Lasso catheter based on its electrogram characteristics. The algorithm stopped the catheter when it located of an AF source and identified the type. The efficiency of the algorithm is validated using a set of simulated focal and rotor-driven arrhythmias in fibrotic human 2D and 3D atrial tissue. RESULTS: Our study shows the feasibility of locating AF sources with a success rate of greater than 95.25% within average 7.56 ± 2.28 placements independently of the initial position of the catheter and the source type. CONCLUSIONS: The algorithm could play a critical role in clinical electrophysiology laboratories for mapping patient-specific ablation of AF sources located outside the pulmonary veins and improving the procedure success.

Iterative navigation of multipole diagnostic catheters to locate repeating-pattern atrial fibrillation drivers.

INTRODUCTION: Targeting repeating-pattern atrial fibrillation (AF) sources (reentry or focal drivers) can help in patient-specific ablation therapy for AF; however, the development of reliable and accurate tools for locating such sources remains a major challenge. We describe iterative catheter navigation (ICAN) algorithm to locate AF drivers using a conventional circular Lasso catheter. METHODS AND RESULTS: At each step, the algorithm analyzes 10 bipolar electrograms recoded at a given catheter location and the history of previous catheter movements to determine if the source is inside the catheter loop. If not, it calculates new coordinates and selects a new position for the catheter. The process continues until a source is located. The algorithm was evaluated in a computer model of atrial tissue with various degrees of fibrosis under a broad range of arrhythmia scenarios. The latter included slow and fast reentry, macroreentry, figure-of-eight reentry, and fibrillatory conduction. Depending on the initial distance of the catheter from the source and scenario, it took about 3 to 16 steps to localize an AF source. In 94% of cases, the identified location was within 4 mm from the source, independently of the initial position of the catheter. The algorithm worked equally well in the presence of patchy fibrosis, low-voltage areas, fragmented electrograms, and dominant-frequency gradients. CONCLUSIONS: AF repeating-pattern sources can be localized using circular catheters without the need to map the entire tissue. The proposed algorithm has the potential to become a useful tool for patient-specific ablation of AF sources located outside the pulmonary veins.

Developing an Iterative Tracking Algorithm to Guide a Catheter Towards Atrial Fibrillation Rotor Sources in Simulated Fibrotic Tissue.

Locating atrial fibrillation (AF) rotor sources can help target ablation therapy for AF. Our aim was to develop a catheter-tracking algorithm to locate AF rotor sources using a conventional 20-electrode circular catheter. We simulated rotor-driven arrhythmias in homogeneous and fibrotic human atrial tissue and evaluated the algorithm for different initial catheter positions. The algorithm guided and detected a rotor with a success rate of greater than 97.9% independently of the initial position of the catheter with an accuracy of greater than 2.3±1.4 mm.

Development of a Rotor-Mapping Algorithm to Locate Ablation Targets During Atrial Fibrillation.

Catheter ablation therapy involving isolation of pulmonary veins (PVs) remains the cornerstone procedure to treat AF. However, due to the sub-optimal success rates of PV isolation, there is a need for new ablation techniques to locate AF ablation targets known as rotors, outside of the PVs. In this paper, we developed a novel rotor-mapping algorithm that uses a conventional diagnostic catheter, Lasso, to locate a rotor source. The algorithm, called the Region of Rotor (ROR) Mapping, utilizes the characteristics of local bipolar electrograms to navigate the catheter's iterative placements while generating a map, overlaid on the atrial anatomy, that displays the potential rotor region. We evaluated the developed ROR mapping algorithm using a 2D simulation of AF on a tissue with heterogeneous conduction properties. The results demonstrated a significant success rate of 93% in accurately locating the region of the rotor with a mean distance of 1.4mm from the ground truth trajectory. The algorithm could play a critical role in mapping non-PV AF ablation targets and improving the outcome of AF ablation.

Determining the light scattering and absorption parameters from forward-directed flux measurements in cardiac tissue.

We describe a method to accurately measure the light scattering model parameters from forward-directed flux (FDF) measurements carried out with a fiber-optic probe (optrode). Improved determination of light scattering parameters will, in turn, permit better modeling and interpretation of optical mapping in the heart using voltage-sensitive dyes. Using our optrode-based system, we carried out high spatial resolution measurements of FDF in intact and homogenized cardiac tissue, as well as in intralipid-based tissue phantoms. The samples were illuminated with a broad collimated beam at 660 and 532 nm. Measurements were performed with a plunge fiber-optic probe (NA=0.22) at a spatial resolution of up to 10 µm. In the vicinity of the illuminated surface, the FDF consistently manifested a fast decaying exponent with a space constant comparable with the decay rate of ballistic photons. Using a Monte Carlo model, we obtained a simple empirical formula linking the rate of the fast exponent to the scattering coefficient, the anisotropy parameter g, and the numerical aperture of the probe. The estimates of scattering coefficient based on this formula were validated in tissue phantoms. Potential applications of optical fiber-based FDF measurements for the evaluation of optical parameters in turbid media are discussed.

A novel catheter-guidance algorithm for localization of atrial fibrillation rotor and focal sources.

Locating atrial fibrillation (AF) focal and rotor sources can help improve target ablation therapy for AF. However, it remains unclear how to use the information provided by multi-polar diagnostic catheters (MPDC) to locate AF sources. Our aim was to develop a catheter-guidance algorithm to locate AF focal and rotor sources using a conventional MPDC. We simulated a 10 cm × 10 cm atrial tissue with focal and rotor sources using the Nygren et al. ionic model. We modeled a Lasso MPDC with 20-unipole electrodes placed with a spacing of 4.5-1-4.5 mm (diameter, d=15 mm) along a circle to obtain 10-bipole electrograms. Staring from an initial location, the algorithm, which was blinded to the location and type of the AF source, iteratively advanced the MPDC by moving its center to the location of the first activated bipole (FAB). The algorithm located an AF source if a stopping condition for either source was satisfied using bipole electrogram characteristics extracted from the MPDC placement. We tested the algorithm for a single rotor and focal source for all possible initial positions on the simulated tissue and repeated it for a random placement with a maximum of 20 possible placements at every trial. The algorithm located the AF source for 100% of trials and on average required 5.99 ± 1.92 placements to an AF source. This algorithm may be used to iteratively direct an MPDC towards an AF source and allow the AF source to be localized for customized AF ablation.

Catheter simulator software tool to generate electrograms of any multi-polar diagnostic catheter from 3D atrial tissue.

Simulations are excellent tools for assessing new therapeutic strategies and are often conducted before implementing new therapy options in a clinical practice. For patients suffering from a heart arrhythmia, the main source of information comes from an intracardiac catheter. One of the common catheters is a Lasso multi-pole diagnostic catheter, which is a catheter that has 20 electrodes in a circular pattern. In this paper, we developed algorithm and simulation software that allows the users to place a multi-pole catheter on the atrial endocardial surface and record electrograms. In 3D atrial tissue, the plane of principal curvature is determined using eigenvectors of catheter vertices, from where the normals are projected and registered to the surface using 3D geodesic distance. This tool provides a platform for performing customized virtual cardiac experiments.

Characterization of Electrograms from Multipolar Diagnostic Catheters during Atrial Fibrillation.

Atrial fibrillation (AF) is the most common arrhythmia in USA with more than 2.3 million people affected annually. Catheter ablation procedure is a method for treatment of AF, which involves 3D electroanatomic mapping of the patient's left atrium (LA) by maneuvering a conventional multipolar diagnostic catheter (MPDC) along the LA endocardial surface after which pulmonary vein (PV) isolation is performed, thus eliminating the AF triggers originating from the PVs. However, it remains unclear how to effectively utilize the information provided by the MPDC to locate the AF-sustaining sites, known as sustained rotor-like activities (RotAs). In this study, we use computer modeling to investigate the variations in the characteristics of the MPDC electrograms, namely, total conduction delay (TCD) and average cycle length (CL), as the MPDC moves towards a RotA source. Subsequently, a study with a human subject was performed in order to verify the predictions of the simulation study. The conclusions from this study may be used to iteratively direct an MPDC towards RotA sources thus allowing the RotAs to be localized for customized and improved AF ablation.

Refraction of scroll-wave filaments at the boundary between two reaction-diffusion media.

We explore the shape and the dynamics of scroll-wave filaments in excitable media with an abruptly changing diffusion tensor, important for cardiac applications. We show that, similar to a beam of light, the filament refracts at the boundary separating domains with different diffusion. We derive the laws of filament refraction and test their validity in computational experiments. We discovered that at small angles to the interface, the filament can become unstable and develop oscillations. The nature of the observed instabilities, as well as overall theoretical and experimental significance of the findings, is discussed.

Evolution of action potential alternans in rabbit heart during acute regional ischemia.

This study investigates the development of the spatiotemporal pattern of action potential alternans during acute regional ischemia. Experiments were carried out in isolated Langendorff-perfused rabbit heart using a combination of optical mapping and microelectrode recordings. The alternans pattern significantly changed over time and had a biphasic character reaching maximum at 6-9 min after occlusion. Phase I (3-11 minutes of ischemia) is characterized by rapid increase in the alternans magnitude and expansion of the alternans territory. Phase I is followed by gradual decline of alternans (Phase II) in both magnitude and territory. During both phases we observed significant beat-to-beat variations of the optical action potential amplitude (OAPA) alternans. Simultaneous microelectrode recordings from subepicardial and subendocardial layers showed that OAPA alternans coincided with intramural 2 : 1 conduction blocks. Our findings are consistent with the modeling studies predicting that during acute regional ischemia alternans can be driven by 2 : 1 conduction blocks in the ischemic region.

Cardiac response to low-energy field pacing challenges the standard theory of defibrillation.

BACKGROUND: The electric response of myocardial tissue to periodic field stimuli has attracted significant attention as the basis for low-energy antifibrillation pacing, potentially more effective than traditional single high-energy shocks. In conventional models, an electric field produces a highly nonuniform response of the myocardial wall, with discrete excitations, or hot spots (HS), occurring at cathodal tissue surfaces or large coronary vessels. We test this prediction using novel 3-dimensional tomographic optical imaging. METHODS AND RESULTS: Experiments were performed in isolated coronary perfused pig ventricular wall preparations stained with near-infrared voltage-sensitive fluorescent dye DI-4-ANBDQBS. The 3-dimensional coordinates of HS were determined using alternating transillumination. To relate HS formation with myocardial structures, we used ultradeep confocal imaging (interrogation depths, >4 mm). The peak HS distribution is located deep inside the heart wall, and the depth is not significantly affected by field polarity. We did not observe the strong colocalization of HS with major coronary vessels anticipated from theory. Yet, we observed considerable lateral displacement of HS with field polarity reversal. Models that de-emphasized lateral intracellular coupling and accounted for resistive heterogeneity in the extracellular space showed similar HS distributions to the experimental observations. CONCLUSIONS: The HS distributions within the myocardial wall and the significant lateral displacements with field polarity reversal are inconsistent with standard theories of defibrillation. Extended theories based on enhanced descriptions of cellular scale electric mechanisms may be necessary. The considerable lateral displacement of HS with field polarity reversal supports the hypothesis of biphasic stimuli in low-energy antifibrillation pacing being advantageous.

Embryonic heart morphogenesis from confocal microscopy imaging and automatic segmentation.

Embryonic heart morphogenesis (EHM) is a complex and dynamic process where the heart transforms from a single tube into a four-chambered pump. This process is of great biological and clinical interest but is still poorly understood for two main reasons. On the one hand, the existing imaging modalities for investigating EHM suffered from either limited penetration depth or limited spatial resolution. On the other hand, current works typically adopted manual segmentation, which was tedious, subjective, and time consuming considering the complexity of developing heart geometry and the large size of images. In this paper, we propose to utilize confocal microscopy imaging with tissue optical immersion clearing technique to image the heart at different stages of development for EHM study. The imaging method is able to produce high spatial resolution images and achieve large penetration depth at the same time. Furthermore, we propose a novel convex active contour model for automatic image segmentation. The model has the ability to deal with intensity fall-off in depth which is characterized by confocal microscopy images. We acquired the images of embryonic quail hearts from day 6 to day 14 of incubation for EHM study. The experimental results were promising and provided us with an insight view of early heart growth pattern and also paved the road for data-driven heart growth modeling.

Optical mapping at increased illumination intensities.

Voltage-sensitive fluorescent dyes have become a major tool in cardiac and neuro-electrophysiology. Achieving high signal-to-noise ratios requires increased illumination intensities, which may cause photobleaching and phototoxicity. The optimal range of illumination intensities varies for different dyes and must be evaluated individually. We evaluate two dyes: di-4-ANBDQBS (excitation 660 nm) and di-4-ANEPPS (excitation 532 nm) in the guinea pig heart. The light intensity varies from 0.1 to 5 mW/mm2, with the upper limit at 5 to 10 times above values reported in the literature. The duration of illumination was 60 s, which in guinea pigs corresponds to 300 beats at a normal heart rate. Within the identified duration and intensity range, neither dye shows significant photobleaching or detectable phototoxic effects. However, light absorption at higher intensities causes noticeable tissue heating, which affects the electrophysiological parameters. The most pronounced effect is a shortening of the action potential duration, which, in the case of 532-nm excitation, can reach ∼30%. At 660-nm excitation, the effect is ∼10%. These findings may have important implications for the design of optical mapping protocols in biomedical applications.

Anchoring of drifting spiral and scroll waves to impermeable inclusions in excitable media.

Anchoring of spiral and scroll waves in excitable media has attracted considerable interest in the context of cardiac arrhythmias. Here, by bombarding inclusions with drifting spiral and scroll waves, we explore the forces exerted by inclusions onto an approaching spiral and derive the equations of motion governing spiral dynamics in the vicinity of inclusion. We demonstrate that these forces nonmonotonically depend on distance and can lead to complex behavior: (a) anchoring to small but circumnavigating larger inclusions; (b) chirality-dependent anchoring.

Extracting surface activation time from the optically recorded action potential in three-dimensional myocardium.

Optical mapping has become an indispensible tool for studying cardiac electrical activity. However, due to the three-dimensional nature of the optical signal, the optical upstroke is significantly longer than the electrical upstroke. This raises the issue of how to accurately determine the activation time on the epicardial surface. The purpose of this study was to establish a link between the optical upstroke and exact surface activation time using computer simulations, with subsequent validation by a combination of microelectrode recordings and optical mapping experiments. To simulate wave propagation and associated optical signals, we used a hybrid electro-optical model. We found that the time of the surface electrical activation (t(E)) within the accuracy of our simulations coincided with the maximal slope of the optical upstroke (t(F)*) for a broad range of optical attenuation lengths. This was not the case when the activation time was determined at 50% amplitude (t(F50)) of the optical upstroke. The validation experiments were conducted in isolated Langendorff-perfused rat hearts and coronary-perfused pig left ventricles stained with either di-4-ANEPPS or the near-infrared dye di-4-ANBDQBS. We found that t(F)* was a more accurate measure of t(E) than was t(F50) in all experimental settings tested (P = 0.0002). Using t(F)* instead of t(F50) produced the most significant improvement in measurements of the conduction anisotropy and the transmural conduction time in pig ventricles.

Imaging electrical excitation inside the myocardial wall.

Cardiac arrhythmias are often triggered by ectopic membrane depolarization originating deep inside the myocardial wall. Here we propose a new method utilizing a novel near-infrared voltage-sensitive fluorescent dye DI-4-ANBDQBS to determine the three-dimensional (3D) coordinates of the sources of such depolarization. We tested the method in live preparations of pig left and right ventricular myocardium (thickness 8-18 mm) and phantoms imitating the optical properties of myocardial tissue. The method utilizes an alternating transillumination approach that involves comparing pairs of simultaneously recorded broad-field epifluorescence and transillumination images produced at two alternating directions of illumination. Recordings were taken simultaneously by two CCD cameras facing the endocardial and epicardial surfaces of the heart at a frame rate up to 3 KHz. In live preparations, we were able to localize the origin of the depolarization wave with a precision of ±1.3mm in the transmural direction and 3 mm in the image plane. The accuracy of detection was independent of the depth of the source inside ventricular wall.

Probing field-induced tissue polarization using transillumination fluorescent imaging.

Despite major successes of biophysical theories in predicting the effects of electrical shocks within the heart, recent optical mapping studies have revealed two major discrepancies between theory and experiment: 1), the presence of negative bulk polarization recorded during strong shocks; and 2), the unexpectedly small surface polarization under shock electrodes. There is little consensus as to whether these differences result from deficiencies of experimental techniques, artifacts of tissue damage, or deficiencies of existing theories. Here, we take advantage of recently developed near-infrared voltage-sensitive dyes and transillumination optical imaging to perform, for the first time that we know of, noninvasive probing of field effects deep inside the intact ventricular wall. This technique removes some of the limitations encountered in previous experimental studies. We explicitly demonstrate that deep inside intact myocardial tissue preparations, strong electrical shocks do produce considerable negative bulk polarization previously inferred from surface recordings. We also demonstrate that near-threshold diastolic field stimulation produces activation of deep myocardial layers 2-6 mm away from the cathodal surface, contrary to theory. Using bidomain simulations we explore factors that may improve the agreement between theory and experiment. We show that the inclusion of negative asymmetric current can qualitatively explain negative bulk polarization in a discontinuous bidomain model.