Theoretical considerations for mapping activation in human cardiac fibrillation

The physics of heart rhythm disorders

The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.

Control of the chirality of spiral waves and recreation of spatial excitation patterns through optogenetics

Spiral waves lead to dangerous arrhythmias in the cardiac system. In 2015 Burton et al. demonstrated the reversal of the spiral wave chirality through the rotating spiral-shaped illumination on the optogenetically modified cardiac monolayers. We show that this process entails the recreation of a spiral wave. We show how this methodology can be used to control and create the desired spatial excitation pattern. We found that the control is sensitive to the area of illuminated region but independent of the phase difference of the existing spiral wave and the applied spiral-shaped light. We also discovered that our methodology can temporarily resynchronize a turbulent system. The results offer numerical evidence for the control of spatial pattern in biological excitable systems with optogenetics.

Finding type and location of the source of cardiac arrhythmias from the averaged flow velocity field using the determinant-trace method

Life threatening cardiac arrhythmias result from abnormal propagation of nonlinear electrical excitation waves in the heart. Finding the locations of the sources of these waves remains a challenging problem. This is mainly due to the low spatial resolution of electrode recordings of these waves. Also, these recordings are subjected to noise. In this paper, we develop a different approach: the AFV-DT method based on an averaged flow velocity (AFV) technique adopted from the analysis of optical flows and the determinant-trace (DT) method used for vector field analysis of dynamical systems. This method can find the location and determine all important types of sources found in excitable media such as focal activity, spiral waves, and waves rotating around obstacles. We test this method on in silico data of various wave excitation patterns obtained using the Luo-Rudy model for cardiac tissue. We show that the method works well for data with low spatial resolutions (up to 8×8) and is stable against noise. Finally, we apply it to two clinical cases and show that it can correctly identify the arrhythmia type and location. We discuss further steps on the development and improvement of this approach.

Mapping Technologies for Catheter Ablation of Atrial Fibrillation Beyond Pulmonary Vein Isolation

Catheter ablation remains the most effective and relatively minimally invasive therapy for rhythm control in patients with AF. Ablation has consistently shown a reduction of arrhythmia-related symptoms and significant improvement in patients’ quality of life compared with medical treatment. The ablation strategy relies on a well-established anatomical approach of effective pulmonary vein isolation. Additional anatomical targets have been reported with the aim of increasing procedure success in complex substrates. However, larger ablated areas with uncertainty of targeting relevant regions for AF initiation or maintenance are not exempt from the potential risk of complications and pro-arrhythmia. Recent developments in mapping tools and computational methods for advanced signal processing during AF have reported novel strategies to identify atrial regions associated with AF maintenance. These novel tools – although mainly limited to research series – represent a significant step forward towards the understanding of complex patterns of propagation during AF and the potential achievement of patient-tailored AF ablation strategies for the near future.

Three dimensional reconstruction to visualize atrial fibrillation activation patterns on curved atrial geometry

BACKGROUND

The rotational activation created by spiral waves may be a mechanism for atrial fibrillation (AF), yet it is unclear how activation patterns obtained from endocardial baskets are influenced by the 3D geometric curvature of the atrium or 'unfolding' into 2D maps. We develop algorithms that can visualize spiral waves and their tip locations on curved atrial geometries. We use these algorithms to quantify differences in AF maps and spiral tip locations between 3D basket reconstructions, projection onto 3D anatomical shells and unfolded 2D surfaces.

METHODS

We tested our algorithms in N = 20 patients in whom AF was recorded from 64-pole baskets (Abbott, CA). Phase maps were generated by non-proprietary software to identify the tips of spiral waves, indicated by phase singularities. The number and density of spiral tips were compared in patient-specific 3D shells constructed from the basket, as well as 3D maps from clinical electroanatomic mapping systems and 2D maps.

RESULTS

Patients (59.4±12.7 yrs, 60% M) showed 1.7±0.8 phase singularities/patient, in whom ablation terminated AF in 11/20 patients (55%). There was no difference in the location of phase singularities, between 3D curved surfaces and 2D unfolded surfaces, with a median correlation coefficient between phase singularity density maps of 0.985 (0.978-0.990). No significant impact was noted by phase singularities location in more curved regions or relative to the basket location (p>0.1).

CONCLUSIONS

AF maps and phase singularities mapped by endocardial baskets are qualitatively and quantitatively similar whether calculated by 3D phase maps on patient-specific curved atrial geometries or in 2D. Phase maps on patient-specific geometries may be easier to interpret relative to critical structures for ablation planning.

Time resolution for wavefront and phase singularity tracking using activation maps in cardiac propagation models

The dynamics of cardiac fibrillation can be described by the number, the trajectory, the stability, and the lifespan of phase singularities (PSs). Accurate PS tracking is straightforward in simple uniform tissues but becomes more challenging as fibrosis, structural heterogeneity, and strong anisotropy are combined. In this paper, we derive a mathematical formulation for PS tracking in two-dimensional reaction-diffusion models. The method simultaneously tracks wavefronts and PS based on activation maps at full spatiotemporal resolution. PS tracking is formulated as a linear assignment problem solved by the Hungarian algorithm. The cost matrix incorporates information about distances between PS, chirality, and wavefronts. A graph of PS trajectories is generated to represent the creations and annihilations of PS pairs. Structure-preserving graph transformations are applied to provide a simplified description at longer observation time scales. The approach is validated in 180 simulations of fibrillation in four different types of substrates featuring, respectively, wavebreaks, ionic heterogeneities, fibrosis, and breakthrough patterns. The time step of PS tracking is studied in the range from 0.1 to 10 ms. The results show the benefits of improving time resolution from 1 to 0.1 ms. The tracking error rate decreases by an order of magnitude because the occurrence of simultaneous events becomes less likely. As observed on PS survival curves, the graph-based analysis facilitates the identification of macroscopically stable rotors despite wavefront fragmentation by fibrosis.

Wavefront Field Mapping Reveals a Physiologic Network Between Drivers Where Ablation Terminates Atrial Fibrillation

BACKGROUND

Localized drivers are proposed mechanisms for persistent atrial fibrillation (AF) from optical mapping of human atria and clinical studies of AF, yet are controversial because drivers fluctuate and ablating them may not terminate AF. We used wavefront field mapping to test the hypothesis that AF drivers, if concurrent, may interact to produce fluctuating areas of control to explain their appearance/disappearance and acute impact of ablation.

METHODS

We recruited 54 patients from an international registry in whom persistent AF terminated by targeted ablation. Unipolar AF electrograms were analyzed from 64-pole baskets to reconstruct activation times, map propagation vectors each 20 ms, and create nonproprietary phase maps.

RESULTS

Each patient (63.6±8.5 years, 29.6% women) showed 4.0±2.1 spatially anchored rotational or focal sites in AF in 3 patterns. First, a single (type I; n=7) or, second, paired chiral-antichiral (type II; n=5) rotational drivers controlled most of the atrial area. Ablation of 1 to 2 large drivers terminated all cases of types I or II AF. Third, interaction of 3 to 5 drivers (type III; n=42) with changing areas of control. Targeted ablation at driver centers terminated AF and required more ablation in types III versus I (P=0.02 in left atrium).

CONCLUSIONS

Wavefront field mapping of persistent AF reveals a pathophysiologic network of a small number of spatially anchored rotational and focal sites, which interact, fluctuate, and control varying areas. Future work should define whether AF drivers that control larger atrial areas are attractive targets for ablation.

Understanding the Beat-to-Beat Variations of P-Waves Morphologies in AF Patients During Sinus Rhythm: A Scoping Review of the Atrial Simulation Studies

The remarkable advances in high-performance computing and the resulting increase of the computational power have the potential to leverage computational cardiology toward improving our understanding of the pathophysiological mechanisms of arrhythmias, such as Atrial Fibrillation (AF). In AF, a complex interaction between various triggers and the atrial substrate is considered to be the leading cause of AF initiation and perpetuation. In electrocardiography (ECG), P-wave is supposed to reflect atrial depolarization. It has been found that even during sinus rhythm (SR), multiple P-wave morphologies are present in AF patients with a history of AF, suggesting a higher dispersion of the conduction route in this population. In this scoping review, we focused on the mechanisms which modify the electrical substrate of the atria in AF patients, while investigating the existence of computational models that simulate the propagation of the electrical signal through different routes. The adopted review methodology is based on a structured analytical framework which includes the extraction of the keywords based on an initial limited bibliographic search, the extensive literature search and finally the identification of relevant articles based on the reference list of the studies. The leading mechanisms identified were classified according to their scale, spanning from mechanisms in the cell, tissue or organ level, and the produced outputs. The computational modeling approaches for each of the factors that influence the initiation and the perpetuation of AF are presented here to provide a clear overview of the existing literature. Several levels of categorization were adopted while the studies which aim to translate their findings to ECG phenotyping are highlighted. The results denote the availability of multiple models, which are appropriate under specific conditions. However, the consideration of complex scenarios taking into account multiple spatiotemporal scales, personalization of electrophysiological and anatomical models and the reproducibility in terms of ECG phenotyping has only partially been tackled so far.

Robust approach for rotor mapping in cardiac tissue

The motion of and interaction between phase singularities that lie at the centers of spiral waves capture many qualitative and, in some cases, quantitative features of complex dynamics in excitable systems. Being able to accurately reconstruct their position is thus quite important, even if the data are noisy and sparse, as in electrophysiology studies of cardiac arrhythmias, for instance. A recently proposed global topological approach [Marcotte and Grigoriev, Chaos 27, 093936 (2017)] promises to meaningfully improve the quality of the reconstruction compared with traditional, local approaches. Indeed, we found that this approach is capable of handling noise levels exceeding the range of the signal with minimal loss of accuracy. Moreover, it also works successfully with data sampled on sparse grids with spacing comparable to the mean separation between the phase singularities for complex patterns featuring multiple interacting spiral waves.

Far-field contributions in multi-electrodes atrial recordings blur distinction between anatomical and functional reentries and may cause imaginary phase singularities - A computational study

BACKGROUND

Atrial fibrillation (AF) is the most common cardiac arrhythmia and the most important cause of embolic stroke, requiring new technologies for its better understanding and therapies. Recent approaches to map the electrical activity during AF with multi-electrode systems aim at localizing patient-specific ablation targets of reentrant patterns. However, there is a critical need to determine the accuracy of those mapping systems. We performed computer simulations as a numerical approach of systematically evaluating the influence of far-field sources on the electrical recordings and detection of rotors.

METHODS

We constructed 2 computer models of atrial tissue: (i) a 2D sheet model with varying non-active cells area in its center, and (ii) a whole realistic 3D atrial model. Phase maps were built based on the Hilbert transform of the unipolar electrograms recorded by virtual 2D and 3D multi-electrode systems and rotors were tracked through phase singularities detections.

RESULTS

Analysis of electrograms recorded away from the 2D atrial model shows that the larger the distance between an electrode and the tissue model, the stronger the far-field sources contribution to the electrogram is. Importantly, even if an electrode is positioned in contact with the tissue, the electrogram contains significant contributions from distal sources that blur the distinction between anatomical and functional reentries. Moreover, when mapping the 3D atrial model, remote activity generated false phase singularities at locations without local reentrant excitation patterns.

CONCLUSIONS

Far-field contributions to electrograms during AF reduce the accuracy of detecting and interpreting reentrant activity.

Characterizing Electrogram Signal Fidelity and the Effects of Signal Contamination on Mapping Human Persistent Atrial Fibrillation

Objective: Determining accurate intracardiac maps of atrial fibrillation (AF) in humans can be difficult, owing primarily to various sources of contamination in electrogram signals. The goal of this study is to develop a measure for signal fidelity and to develop methods to quantify robustness of observed rotational activity in phase maps subject to signal contamination.

Methods: We identified rotational activity in phase maps of human persistent AF using the Hilbert transform of sinusoidally recomposed signals, where localized ablation at rotational sites terminated fibrillation. A novel measure of signal fidelity was developed to quantify signal quality. Contamination is then introduced to the underlying electrograms by removing signals at random, adding noise to computations of cycle length, and adding realistic far-field signals. Mean tip number N and tip density δ, defined as the proportion of time a region contains a tip, at the termination site are computed to compare the effects of contamination.

Results: Domains of low signal fidelity correspond to the location of rotational cores. Removing signals and altering cycle length accounted for minor changes in tip density, while targeted removal of low fidelity electrograms can result in a significant increase in tip density and stability. Far-field contamination was found to obscure rotation at the termination site.

Conclusion: Rotational activity in clinical AF can produce domains of low fidelity electrogram recordings at rotational cores. Observed rotational patterns in phase maps appear most sensitive to far-field activation. These results may inform novel methods to map AF in humans which can be tested directly in patients at electrophysiological study and ablation.

Not all rotors, effective ablation targets for nonparoxysmal atrial fibrillation, are included in areas suggested by conventional indirect indicators of atrial fibrillation drivers: ExTRa Mapping project

Background

Effects of nonparoxysmal atrial fibrillation (non-PAF) ablation targeting complex fractionated atrial electrogram (CFAE) areas and/or low voltage areas (LVAs) are still controversial.

Methods and Results

A recently developed online real-time phase mapping system (ExTRa Mapping) was used to conduct LVA mapping and simultaneous ExTRa and CFAE mapping in 28 non-PAF patients after pulmonary vein isolation (PVI). Nonpassively activated areas, in the form of meandering rotors and/or multiple wavelets assumed to contain non-PAF drivers, partly overlapped with CFAE/LVAs but not always coincided with them.

Conclusion

Real-time rotor imaging, rather than conventional indirect indicators only, might be very useful for detecting non-PAF drivers.

Inflection point analysis: A machine learning approach for extraction of IEGM active intervals during atrial fibrillation

OBJECTIVE

In this paper, we propose a novel algorithm to extract the active intervals of intracardiac electrograms during atrial fibrillation.

METHODS

First, we show that the characteristics of the signal waveform at its inflection points are prominent features that are implicitly used by human annotators for distinguishing between active and inactive intervals of IEGMs. Then, we show that the natural logarithm of features corresponding to active and inactive intervals exhibits a mixture of two Gaussian distributions in three dimensional feature space. An Expectation Maximization algorithm for Gaussian mixtures is then applied for automatic clustering of the features into two categories.

RESULTS

The absolute error in onset and offset estimation of active intervals is 6.1ms and 10.7ms, respectively, guaranteeing a high resolution. The true positive rate for the proposed method is also 98.1%, proving the high reliability.

CONCLUSION

The proposed method can extract the active intervals of IEGMs during AF with a high accuracy and resolution close to manually annotated results.

SIGNIFICANCE

In contrast with some of the conventional methods, no windowing technique is required in our approach resulting in significantly higher resolution in estimating the onset and offset of active intervals. Furthermore, since the signal characteristics at inflection points are analyzed instead of signal samples, the computational time is significantly low, ensuring the real-time application of our algorithm. The proposed method is also robust to noise and baseline variations thanks to the Laplacian of Gaussian filter employed for extraction of inflection points.

Factors affecting basket catheter detection of real and phantom rotors in the atria: A computational study

Anatomically based procedures to ablate atrial fibrillation (AF) are often successful in terminating paroxysmal AF. However, the ability to terminate persistent AF remains disappointing. New mechanistic approaches use multiple-electrode basket catheter mapping to localize and target AF drivers in the form of rotors but significant concerns remain about their accuracy. We aimed to evaluate how electrode-endocardium distance, far-field sources and inter-electrode distance affect the accuracy of localizing rotors. Sustained rotor activation of the atria was simulated numerically and mapped using a virtual basket catheter with varying electrode densities placed at different positions within the atrial cavity. Unipolar electrograms were calculated on the entire endocardial surface and at each of the electrodes. Rotors were tracked on the interpolated basket phase maps and compared with the respective atrial voltage and endocardial phase maps, which served as references. Rotor detection by the basket maps varied between 35-94% of the simulation time, depending on the basket's position and the electrode-to-endocardial wall distance. However, two different types of phantom rotors appeared also on the basket maps. The first type was due to the far-field sources and the second type was due to interpolation between the electrodes; increasing electrode density decreased the incidence of the second but not the first type of phantom rotors. In the simulations study, basket catheter-based phase mapping detected rotors even when the basket was not in full contact with the endocardial wall, but always generated a number of phantom rotors in the presence of only a single real rotor, which would be the desired ablation target. Phantom rotors may mislead and contribute to failure in AF ablation procedures.

A statistical model of false negative and false positive detection of phase singularities

The complexity of cardiac fibrillation dynamics can be assessed by analyzing the distribution of phase singularities (PSs) observed using mapping systems. Interelectrode distance, however, limits the accuracy of PS detection. To investigate in a theoretical framework the PS false negative and false positive rates in relation to the characteristics of the mapping system and fibrillation dynamics, we propose a statistical model of phase maps with controllable number and locations of PSs. In this model, phase maps are generated from randomly distributed PSs with physiologically-plausible directions of rotation. Noise and distortion of the phase are added. PSs are detected using topological charge contour integrals on regular grids of varying resolutions. Over 100 × 10⁶ realizations of the random field process are used to estimate average false negative and false positive rates using a Monte-Carlo approach. The false detection rates are shown to depend on the average distance between neighboring PSs expressed in units of interelectrode distance, following approximately a power law with exponents in the range of 1.14 to 2 for false negatives and around 2.8 for false positives. In the presence of noise or distortion of phase, false detection rates at high resolution tend to a non-zero noise-dependent lower bound. This model provides an easy-to-implement tool for benchmarking PS detection algorithms over a broad range of configurations with multiple PSs.

Demonstration of cardiac rotor and source mapping techniques in embryonic chick monolayers

Excitable media, such as the heart, display propagating waves with different geometries including target patterns and rotors (spiral waves). Collision of two waves leads to annihilation of both. We present algorithms for data processing and analysis to identify the core of rotors. In this work, we show that as the spatial sampling resolution decreases it becomes increasingly difficult to identify rotors-there are instances of false negatives and false positives. These observations are relevant to current controversies concerning the role of rotors in the initiation, maintenance, and treatment of cardiac arrhythmias, especially atrial fibrillation. Currently some practitioners target the core of rotors for ablation, but the effectiveness of this procedure has been questioned. In view of the difficulties inherent in the identification of rotors, we conclude that methods to identify rotors need to first be validated prior to assessing the efficacy of ablation.

Spatial Resolution Requirements for Accurate Identification of Drivers of Atrial Fibrillation

Supplemental Digital Content is available in the text.

Background—

Recent studies have demonstrated conflicting mechanisms underlying atrial fibrillation (AF), with the spatial resolution of data often cited as a potential reason for the disagreement. The purpose of this study was to investigate whether the variation in spatial resolution of mapping may lead to misinterpretation of the underlying mechanism in persistent AF.

Methods and Results—

Simulations of rotors and focal sources were performed to estimate the minimum number of recording points required to correctly identify the underlying AF mechanism. The effects of different data types (action potentials and unipolar or bipolar electrograms) and rotor stability on resolution requirements were investigated. We also determined the ability of clinically used endocardial catheters to identify AF mechanisms using clinically recorded and simulated data. The spatial resolution required for correct identification of rotors and focal sources is a linear function of spatial wavelength (the distance between wavefronts) of the arrhythmia. Rotor localization errors are larger for electrogram data than for action potential data. Stationary rotors are more reliably identified compared with meandering trajectories, for any given spatial resolution. All clinical high-resolution multipolar catheters are of sufficient resolution to accurately detect and track rotors when placed over the rotor core although the low-resolution basket catheter is prone to false detections and may incorrectly identify rotors that are not present.

Conclusions—

The spatial resolution of AF data can significantly affect the interpretation of the underlying AF mechanism. Therefore, the interpretation of human AF data must be taken in the context of the spatial resolution of the recordings.

Spatiotemporal Progression of Early Human Ventricular Fibrillation

Objectives

The objective of this study was to evaluate the spatio-temporal organization and progression of human ventricular fibrillation (VF) in the left (LV) and right (RV) ventricles.

Background

Studies suggest that localized sources contribute to VF maintenance, but the evolution of VF episodes has not been quantified.

Methods

Synchrony between electrograms recorded from 25 patients with induced VF is computed and used to define the Asynchronous Index (ASI), indicating regions which are out-of-step with surrounding tissue. Computer simulations show that ASI can identify the location of VF-maintaining sources, where larger values of ASI_max correlate with more stable sources.

Results

Automated synchrony analysis shows elevated values of ASI in a majority of self-terminating episodes (LV: 8/9, RV: 7/8) and sustained episodes (LV: 11/11, RV: 12/12). The locations of ASI_max in sustained episodes co-localize with rotor cores when rotational activity is simultaneously present in phase maps (LV: 8/8, RV: 5/7, p<.05). The distribution of ASI_max differentiates self-terminating from sustained episodes (mean ASI_max = 0.60±0.14 and 0.70±0.16, respectively; p=0.01). Across sustained episodes the LV exhibits an increase in ASI_max with time.

Conclusions

Quantitative analysis identifies localized asynchronous regions that correlate with sources in VF, with sustained episodes evolving to exhibit more stable activation in the LV. This successive increase in stability indicates a stabilizing agent may be responsible for perpetuating fibrillation in a "migrate-and-capture" mechanism in the LV.

Basket-Type Catheters: Diagnostic Pitfalls Caused by Deformation and Limited Coverage

Whole-chamber mapping using a 64-pole basket catheter (BC) has become a featured approach for the analysis of excitation patterns during atrial fibrillation. A flexible catheter design avoids perforation but may lead to spline bunching and influence coverage. We aim to quantify the catheter deformation and endocardial coverage in clinical situations and study the effect of catheter size and electrode arrangement using an in silico basket model. Atrial coverage and spline separation were evaluated quantitatively in an ensemble of clinical measurements. A computational model of the BC was implemented including an algorithm to adapt its shape to the atrial anatomy. Two clinically relevant mapping positions in each atrium were assessed in both clinical and simulated data. The simulation environment allowed varying both BC size and electrode arrangement. Results showed that interspline distances of more than 20 mm are common, leading to a coverage of less than 50% of the left atrial (LA) surface. In an ideal in silico scenario with variable catheter designs, a maximum coverage of 65% could be reached. As spline bunching and insufficient coverage can hardly be avoided, this has to be taken into account for interpretation of excitation patterns and development of new panoramic mapping techniques.

Determining conduction patterns on a sparse electrode grid: Implications for the analysis of clinical arrhythmias

We present a general method of utilizing bioelectric recordings from a spatially sparse electrode grid to compute a dynamic vector field describing the underlying propagation of electrical activity. This vector field, termed the wave-front flow field, permits quantitative analysis of the magnitude of rotational activity (vorticity) and focal activity (divergence) at each spatial point. We apply this method to signals recorded during arrhythmias in human atria and ventricles using a multipolar contact catheter and show that the flow fields correlate with corresponding activation maps. Further, regions of elevated vorticity and divergence correspond to sites identified as clinically significant rotors and focal sources where therapeutic intervention can be effective. These flow fields can provide quantitative insights into the dynamics of normal and abnormal conduction in humans and could potentially be used to enhance therapies for cardiac arrhythmias.

Anti-arrhythmic strategies for atrial fibrillation: The role of computational modeling in discovery, development, and optimization

Atrial fibrillation (AF), the most common cardiac arrhythmia, is associated with increased risk of cerebrovascular stroke, and with several other pathologies, including heart failure. Current therapies for AF are targeted at reducing risk of stroke (anticoagulation) and tachycardia-induced cardiomyopathy (rate or rhythm control). Rate control, typically achieved by atrioventricular nodal blocking drugs, is often insufficient to alleviate symptoms. Rhythm control approaches include antiarrhythmic drugs, electrical cardioversion, and ablation strategies. Here, we offer several examples of how computational modeling can provide a quantitative framework for integrating multiscale data to: (a) gain insight into multiscale mechanisms of AF; (b) identify and test pharmacological and electrical therapy and interventions; and (c) support clinical decisions. We review how modeling approaches have evolved and contributed to the research pipeline and preclinical development and discuss future directions and challenges in the field.

Maximum likelihood cardiac conduction velocity estimation from sequential mapping in the presence of activation time noise with unknown variances

The cardiac conduction velocity (CV) can be estimated by analysing the activation times (ATs) and the locations of the electrodes that are used for the intracardiac electrogram (IEGM) recording. Here, we study the problem of the CV estimation in sequential mapping without using any independent electrogram as a time alignment reference. We assume that the IEGMs are sequentially recorded from several sites, where at each site, at least two of the catheter's electrodes are in contact with the cardiac tissue. We consider the planar wavefront with stable CV that propagates within the recording sites throughout our data collection period. Assuming the zero-mean Gaussian AT estimation error, we derive the maximum likelihood estimations of the CV and AT at a desired location on the cardiac shell. The CV is estimated when the variance of the AT estimation error and the time delay between the sequential recordings are unknown variables. Our simulation results show that the proposed method can precisely estimate the CV of the planar wavefront.

Cardiac conduction velocity estimation from sequential mapping assuming known Gaussian distribution for activation time estimation error

In this paper, we study the problem of the cardiac conduction velocity (CCV) estimation for the sequential intracardiac mapping. We assume that the intracardiac electrograms of several cardiac sites are sequentially recorded, their activation times (ATs) are extracted, and the corresponding wavefronts are specified. The locations of the mapping catheter's electrodes and the ATs of the wavefronts are used here for the CCV estimation. We assume that the extracted ATs include some estimation errors, which we model with zero-mean white Gaussian noise values with known variances. Assuming stable planar wavefront propagation, we derive the maximum likelihood CCV estimator, when the synchronization times between various recording sites are unknown. We analytically evaluate the performance of the CCV estimator and provide its mean square estimation error. Our simulation results confirm the accuracy of the proposed method and the error analysis of the proposed CCV estimator.

Organized Sources Are Spatially Conserved in Recurrent Compared to Pre-Ablation Atrial Fibrillation: Further Evidence for Non-Random Electrical Substrates

INTRODUCTION

Recurrent atrial fibrillation (AF) after ablation is associated with reconnection of initially isolated pulmonary vein (PV) trigger sites. Substrates are often targeted in addition to PVI, but it is unclear how substrates progress over time. We studied if substrates in recurrent AF are conserved or have developed de novo from pre-ablation AF.

METHODS AND RESULTS

Of 137 patients undergoing Focal Impulse and Rotor Mapping (FIRM) at their index procedure for AF, 29 consecutive patients (60 ± 8 years, 79% persistent) recurred and were also mapped at repeat procedure (21 ± 20 months later) using carefully placed 64-pole baskets and RhythmView(TM) (Topera, Menlo Park, CA, USA) to identify AF sources and disorganized zones. Compared to index AF, recurrent AF had a longer cycle length (177 ± 21 vs. 167 ± 19 milliseconds, P = 0.01). All patients (100%) had 1 or more conserved AF rotors between procedures with surrounding disorganization. The number of sources was similar for recurrent AF post-PVI versus index AF (3.2 ± 1.4 vs. 3.1 ± 1.0, P = 0.79), but was lower for recurrent AF after FIRM+PVI versus index AF (4.4 ± 1.4 vs. 2.9 ± 1.7, P = 0.03). Overall, 81% (61/75) of AF sources lay in conserved regions, while 19% (14/75) were detected de novo.

CONCLUSION

Electrical propagation patterns for recurrent AF after unsuccessful ablation are similar in individual patients to their index AF. These data support temporospatial stability of AF substrates over 1-2 years. Trials should determine the relative benefit of adding substrate mapping and ablation to PVI for recurrent AF.

Ventricular fibrillation: triggers, mechanisms and therapies

Ventricular fibrillation (VF) is a common, life-threatening arrhythmia responsible for significant morbidity and mortality. Due to challenges in safely mapping VF, a comprehensive understanding of its mechanisms remains elusive. Recent findings have provided new insights into mechanisms that sustain early VF. Notably, the central role of electrical rotors and catheter-based ablation of VF rotor substrate have been recently reported. In this article, we will review data regarding four stages of VF: initiation, transition, maintenance and evolution. We will discuss the particular mechanisms for each stage and therapies targeting these mechanisms. We also examine inherited arrhythmia syndromes, including the mechanisms and therapies specific to each. We hope that the overview of VF outlined in this work will assist other investigators in designing future therapies to interrupt this life-threatening arrhythmia.

The effect of electrode density on the interpretation of atrial activation patterns in epicardial mapping of human persistent atrial fibrillation

BACKGROUND

Mechanisms sustaining human persistent atrial fibrillation (AF) remain debated, with significant differences between high-density epicardial and global endocardial mapping studies. A key difference is the density of recording electrodes.

OBJECTIVE

We aimed to determine the differences in the prevalence of different atrial activation patterns, and specifically in the prevalence of rotational activations, with varying densities of bipolar electrodes.

METHODS

Epicardial mapping was performed in 10 patients undergoing cardiac surgery, with bipolar electrograms recorded using a triangular plaque (6.75 cm(2) area; 117 bipoles; 2.5-mm inter-bipole spacing) applied to the left atrial posterior wall or right atrial free wall. Dynamic wavefront mapping based on the timing of atrial electrograms was applied to 2 discrete 10-second AF segments. The spacing between bipolar electrode locations was increased from 2.5 × 3.5 mm in the horizontal and oblique directions to 5.0 × 3.5, 5.0 × 7.1, and 7.5 × 10.6 mm, with wavefront mapping repeated at each density.

RESULTS

As density reduced, there was a significant change in relative proportions of the various activation patterns (F=3.69; P < .001). Simple broad wavefront activations became more prevalent (20% ± 8% to 54% ± 8%; P < .05) and complex patterns became less prevalent (48% ± 8% to 9% ± 8%; P < .05) with reducing density. The prevalence of rotational activity declined with bipole density, from median 5.0% (range 0.9%-12.1%) to 0% (range 0%-1.5%) (P = .03). The largest change occurred between inter-bipole spacings of 5.0 × 3.5 and 5.0 × 7.1 mm.

CONCLUSION

Apparent activation patterns in persistent AF vary significantly with electrode density. Low density underestimates the prevalence of complex and rotational patterns. The largest difference occurs between an inter-bipole spacing of 5.0 × 3.5 and a spacing of 5.0 × 7.1 mm. This may have important implications for mapping technology design.

Phase synchrony reveals organization in human atrial fibrillation

It remains unclear if human atrial fibrillation (AF) is spatially nonhierarchical or exhibits a hierarchy of organization sustained by sources. We utilize activation times obtained at discrete locations during AF to compute the phase synchrony between tissue regions, to examine underlying spatial dynamics throughout both atria. We construct a binary synchronization network and show that this network can accurately define regions of coherence in coarse-grained in silico data. Specifically, domains controlled by spiral waves exhibit regions of high phase synchrony. We then apply this analysis to clinical data from patients experiencing cardiac arrhythmias using multielectrode catheters to simultaneously record from a majority of both atria. We show that pharmaceutical intervention with ibutilide organizes activation by increasing the size of the synchronized domain in AF and quantify the increase in temporal organization when arrhythmia changes from fibrillation to tachycardia. Finally, in recordings from 24 patients in AF we show that the level of synchrony is spatially broad with some patients showing large spatially contiguous regions of synchronization, while in others synchrony is localized to small pockets. Using computer simulations, we show that this distribution is inconsistent with distributions obtained from simulations that mimic multiwavelet reentry but is consistent with mechanisms in which one or more spatially conserved spiral waves is surrounded by tissue in which activation is disorganized.

Mechanisms for the Termination of Atrial Fibrillation by Localized Ablation: Computational and Clinical Studies

BACKGROUND

Human atrial fibrillation (AF) can terminate after ablating localized regions, which supports the existence of localized rotors (spiral waves) or focal drivers. However, it is unclear why ablation near a spiral wave tip would terminate AF and not anchor reentry. We addressed this question by analyzing competing mechanisms for AF termination in numeric simulations, referenced to clinical observations.

METHODS AND RESULTS

Spiral wave reentry was simulated in monodomain 2-dimensional myocyte sheets using clinically realistic rate-dependent values for repolarization and conduction. Heterogeneous models were created by introduction of parameterized variations in tissue excitability. Ablation lesions were applied as nonconducting circular regions. Models confirmed that localized ablation may anchor spiral wave reentry, producing organized tachycardias. Several mechanisms referenced to clinical observations explained termination of AF to sinus rhythm. First, lesions may create an excitable gap vulnerable to invasion by fibrillatory waves. Second, ablation of rotors in regions of low-excitability (from remodeling) produced re-entry in more excitable tissue allowing collision of wavefront and back. Conversely, ablation of rotors in high-excitability regions migrated spiral waves to less excitable tissue, where they detached to collide with nonconducting boundaries. Third, ablation may connect rotors to nonconducting anatomic orifices. Fourth, reentry through slow-conducting channels may terminate if ablation closes these channels.

CONCLUSIONS

Limited ablation can terminate AF by several mechanisms. These data shed light on how clinical AF may be sustained in patients' atria, emphasizing heterogeneities in tissue excitability, slow-conducting channels, and obstacles that are increasingly detectable in patients and should be the focus of future translational studies.

Practical Considerations of Mapping Persistent Atrial Fibrillation With Whole-Chamber Basket Catheters

OBJECTIVES

This study sought to evaluate basket catheter deployment, catheter-tissue contact, and time-space stability of unipolar atrial electrograms (aEGMs) recorded in persistent atrial fibrillation (AF) patients.

BACKGROUND

Panoramic mapping of human AF using multiple-electrode basket catheters may identify AF sources. Although clinical results using this technique are provocative, questions remain about its effectiveness.

METHODS

Data were collected from patients (N = 25) undergoing catheter ablation for AF during the multicenter STARLIGHT (Signal Transfer of Atrial Fibrillation Data to Guide Human Treatment) trial (NCT01765075). Left and right aEGM signals were recorded using basket catheters during baseline AF, following ablation and during sinus rhythm. Data were analyzed for basket deployment, peak-to-peak voltage, and electrogram stability and organization. Electrogram stability and organization were evaluated via time-frequency analysis (TFA).

RESULTS

Basket catheters displayed equatorial bunching when deployed in atria. Interspline spacing ranged from 1.7 to 64.0 mm in the right atrial and from 1.5 to 85.08 mm in the left atrial basket. Approximately one-third of mapping electrodes failed to demonstrate a median peak-to-peak voltage >2× the low-voltage threshold. Time-space stability and organization was observed in 13 of 22 (59.09%) right atrial and 10 of 22 (45.45%) left atrial baskets.

CONCLUSIONS

Despite poor deployment and a large number of low-voltage electrodes, stability and organization was observed in about one-half of the mapped patients. Although this study suggests that basket catheters have limitations for patient-specific AF mapping, concordant activation occurs in some persistent AF patients, which may be amenable to high-density mapping techniques.

Modifying Ventricular Fibrillation by Targeted Rotor Substrate Ablation: Proof-of-Concept from Experimental Studies to Clinical VF

INTRODUCTION

Recent work has suggested a role for organized sources in sustaining ventricular fibrillation (VF). We assessed whether ablation of rotor substrate could modulate VF inducibility in canines, and used this proof-of-concept as a foundation to suppress antiarrhythmic drug-refractory clinical VF in a patient with structural heart disease.

METHODS AND RESULTS

In 9 dogs, we introduced 64-electrode basket catheters into one or both ventricles, used rapid pacing at a recorded induction threshold to initiate VF, and then defibrillated after 18±8 seconds. Endocardial rotor sites were identified from basket recordings using phase mapping, and ablation was performed at nonrotor (sham) locations (7 ± 2 minutes) and then at rotor sites (8 ± 2 minutes, P = 0.10 vs. sham); the induction threshold was remeasured after each. Sham ablation did not alter canine VF induction threshold (preablation 150 ± 16 milliseconds, postablation 144 ± 16 milliseconds, P = 0.54). However, rotor site ablation rendered VF noninducible in 6/9 animals (P = 0.041), and increased VF induction threshold in the remaining 3. Clinical proof-of-concept was performed in a patient with repetitive ICD shocks due to VF refractory to antiarrhythmic drugs. Following biventricular basket insertion, VF was induced and then defibrillated. Mapping identified 4 rotors localized at borderzone tissue, and rotor site ablation (6.3 ± 1.5 minutes/site) rendered VF noninducible. The VF burden fell from 7 ICD shocks in 8 months preablation to zero ICD therapies at 1 year, without antiarrhythmic medications.

CONCLUSIONS

Targeted rotor substrate ablation suppressed VF in an experimental model and a patient with refractory VF. Further studies are warranted on the efficacy of VF source modulation.

Computationally efficient method for localizing the spiral rotor source using synthetic intracardiac electrograms during atrial fibrillation

Atrial fibrillation (AF), the most common sustained cardiac arrhythmia, is an extremely costly public health problem. Catheter-based ablation is a common minimally invasive procedure to treat AF. Contemporary mapping methods are highly dependent on the accuracy of anatomic localization of rotor sources within the atria. In this paper, using simulated atrial intracardiac electrograms (IEGMs) during AF, we propose a computationally efficient method for localizing the tip of the electrical rotor with an Archimedean/arithmetic spiral wavefront. The proposed method deploys the locations of electrodes of a catheter and their IEGMs activation times to estimate the unknown parameters of the spiral wavefront including its tip location. The proposed method is able to localize the spiral as soon as the wave hits three electrodes of the catheter. Our simulation results show that the method can efficiently localize the spiral wavefront that rotates either clockwise or counterclockwise.

Ablation of atrial fibrillation

Ablation is increasingly used to treat AF, since recent trials of pharmacological therapy for AF have been disappointing. Ablation has been shown to improve maintenance of sinus rhythm compared to pharmacological therapy in many multicenter trials, although success rates remain suboptimal. This review will discuss several trends in the field of catheter ablation, including studies to advance our understanding of AF mechanisms in different patient populations, innovations in detecting and classifying AF, use of this information to improve strategies for ablation, technical innovations that have improved the ease and safety of ablation, and novel approaches to surgical therapy and imaging. These trends are likely to further improve results from AF ablation in coming years as it becomes an increasingly important therapeutic option for many patients.

Role of Rotors in the Ablative Therapy of Persistent Atrial Fibrillation

Atrial fibrillation (AF) ablation is increasingly used to maintain sinus rhythm yet its results are sub-optimal, especially in patients with persistent AF or prior unsuccessful procedures. Attempts at improvement have often targeted substrates that sustain AF after it is triggered, yet those mechanisms are debated. Many studies now challenge the concept that AF is driven by self-sustaining disordered wavelets, showing instead that localised drivers (rotors) may drive disorder via a process known as fibrillatory conduction. Novel mapping using wide-area recordings, physiological filtering and phase analysis demonstrates rotors in human AF. Contact mapping with focal impulse and rotor modulation (FIRM) shows that localised ablation at sources can improve procedural success in many populations on long-term follow up and some newer approaches to rotor mapping are qualitatively similar. This review critically evaluates the data on rotor mapping and ablation, which advances our conceptual understanding of AF and holds the promise of substantially improving ablative outcomes in patients with persistent AF.

Mechanisms Underlying AF: Triggers, Rotors, Other?

OPINION STATEMENT

There is ongoing debate regarding the precise mechanisms underlying atrial fibrillation (AF). An improved understanding of these mechanisms is urgently needed to improve interventional strategies to suppress and eliminate AF, since the success of current strategies is suboptimal. At present, guidelines for AF ablation focus on pulmonary vein (PV) isolation for the prevention of arrhythmia. Additional targets are presently unclear, and include additional linear ablation and electrogram-guided substrate modification, without clear mechanistic relevance. PV and non-PV triggers are likely central in the first few seconds of AF initiation. Rapid activation from such triggers interacts with transitional mechanisms including conduction velocity slowing, action potential duration (APD) alternans, and steep APD restitution to cause conduction block and initiate functional reentry. However, complete suppression of potential triggers has proven elusive, and the intra-procedural mapping and targeting of transitional mechanisms has not been reported. A growing body of research implicates electrical rotors and focal sources as central mechanisms for the maintenance of AF. In several recent series, they were observed in nearly all patients with sustained arrhythmia. Ablation of rotor and focal source sites, prior to pulmonary vein isolation, substantially modulated atrial fibrillation in a high proportion of patients, and improved ablation outcomes versus pulmonary vein isolation alone. These results have subsequently been confirmed in multicenter series, and the improved outcomes have been found to persist to a mean follow-up of 3 years. Recently, rotors have been observed by multiple groups using diverse technologies. These findings represent a paradigm shift in AF, focusing on sustaining mechanisms, as is currently done with other arrhythmias such as atrioventricular node reentrant tachycardia. Studies are currently underway to assess the optimal strategy for the application of rotor-based ablation in AF management, including clinical trials on the relative efficacy of rotor-only ablation versus PVI-only ablation, which will inform future practice guidelines.

Tracking rotors with minimal electrodes: modulation index-based strategy

BACKGROUND

High-frequency periodic sources during cardiac fibrillation can be detected by phase mapping techniques. To enable practical therapeutic options for modulating periodic sources (existing techniques require high density multielectrode arrays and real time simultaneous mapping capability), a method to identify electrogram morphologies colocalizing to rotors that can be implemented on few electrograms needs to be devised.

METHOD AND RESULTS

Multichannel ventricular fibrillation electrogram data from 7 isolated human hearts using Langendorff setup and intraoperative clinical data from 2 human hearts were included in the analysis. The spatial locations of rotors were identified using phase maps constructed from 112 electrograms. Electrograms were analyzed for repeating patterns and discriminating signal morphologies around the locations of rotors and nonrotors were identified and quantified. Features were extracted from the unipolar electrogram patterns, which corroborated well with the spatial location of rotors. The results suggest that using the proposed modulation index feature, and as low as 1 sample point in the vicinity of the rotors, an accuracy as high as 86% (P<0.001) was obtained in separating rotor locations versus nonrotor locations. The analysis of bipolar electrogram signatures in the vicinity of the rotor locations suggest that 62.5% of the rotors occur at locations where the bipolar electrogram demonstrates continuous activities during ventricular fibrillation.

CONCLUSIONS

Unipolar electrogram extracted modulation index-based detection of rotors is feasible with few electrodes and has greater detection rate than bipolar approach. This strategy may be suitable for nonarray-based single mapping catheter enabled detection of rotors.

The role of rotors in atrial fibrillation

Despite significant advances in our understanding of atrial fibrillation (AF) mechanisms in the last 15 years, ablation outcomes remain suboptimal. A potential reason is that many ablation techniques focus on anatomic, rather than patient-specific functional targets for ablation. Panoramic contact mapping, incorporating phase analysis, repolarization and conduction dynamics, and oscillations in AF rate, overcomes many prior difficulties with mapping AF. This approach provides evidence that the mechanisms sustaining human AF are deterministic, largely due to stable electrical rotors and focal sources in either atrium. Ablation of such sources (Focal Impulse and Rotor Modulation: FIRM ablation) has been shown to improve ablation outcome compared with conventional ablation alone; independent laboratories directly targeting stable rotors have shown similar results. Clinical trials examining the role of stand-alone FIRM ablation are in progress. Looking forward, translating insights from patient-specific mapping to evidence-based guidelines and clinical practice is the next challenge in improving patient outcomes in AF management.