Electrophysiological Rotor Ablation in In-Silico Modeling of Atrial Fibrillation: Comparisons with Dominant Frequency, Shannon Entropy, and Phase Singularity

Overcoming Uncertainties in Electrogram-Based Atrial Fibrillation Mapping: A Review

In clinical rhythmology, intracardiac bipolar electrograms (EGMs) play a critical role in investigating the triggers and substrates inducing and perpetuating atrial fibrillation (AF). However, the interpretation of bipolar EGMs is ambiguous due to several aspects of electrodes, mapping algorithms and wave propagation dynamics, so it requires several variables to describe the effects of these uncertainties on EGM analysis. In this narrative review, we critically evaluate the potential impact of such uncertainties on the design of cardiac mapping tools on AF-related substrate characterization. Literature suggest uncertainties are due to several variables, including the wave propagation vector, the wave's incidence angle, inter-electrode spacing, electrode size and shape, and tissue contact. The preprocessing of the EGM signals and mapping density will impact the electro-anatomical representation and the features extracted from the local electrical activities. The superposition of multiple waves further complicates EGM interpretation. The inclusion of these uncertainties is a nontrivial problem but their consideration will yield a better interpretation of the intra-atrial dynamics in local activation patterns. From a translational perspective, this review provides a concise but complete overview of the critical variables for developing more precise cardiac mapping tools.

Restitution Slope Affects the Outcome of Dominant Frequency Ablation in Persistent Atrial Fibrillation: CUVIA-AF2 Post-Hoc Analysis Based on Computational Modeling Study

Introduction

Although the dominant frequency (DF) localizes the reentrant drivers and the maximal slope of the action potential duration (APD) restitution curve (Smax) reflects the tendency of the wave-break, their interaction has never been studied. We hypothesized that DF ablation has different effects on atrial fibrillation (AF) depending on Smax.

Methods

We studied the DF and Smax in 25 realistic human persistent AF model samples (68% male, 60 ± 10 years old). Virtual AF was induced by ramp pacing measuring Smax, followed by spatiotemporal DF evaluation for 34 s. We assessed the DF ablation effect depending on Smax in both computational modeling and a previous clinical trial, CUVIA-AF (170 patients with persistent AF, 70.6% male, 60 ± 11 years old).

Results

Mean DF had an inverse relationship with Smax regardless of AF acquisition timing (p < 0.001). Virtual DF ablations increased the defragmentation rate compared to pulmonary vein isolation (PVI) alone (p = 0.015), especially at Smax <1 (61.5 vs. 7.7%, p = 0.011). In post-DF ablation defragmentation episodes, DF was significantly higher (p = 0.002), and Smax was lower (p = 0.003) than in episodes without defragmentation. In the post-hoc analysis of CUVIA-AF2, we replicated the inverse relationship between Smax and DF (r = −0.47, p < 0.001), and we observed better rhythm outcomes of clinical DF ablations in addition to a PVI than of empirical PVI at Smax <1 [hazard ratio 0.45, 95% CI (0.22–0.89), p = 0.022; log-rank p = 0.021] but not at ≥ 1 (log-rank p = 0.177).

Conclusion

We found an inverse relationship between DF and Smax and the outcome of DF ablation after PVI was superior at the condition with Smax <1 in both in-silico and clinical trials.

Rotor Localization and Phase Mapping of Cardiac Excitation Waves Using Deep Neural Networks

The analysis of electrical impulse phenomena in cardiac muscle tissue is important for the diagnosis of heart rhythm disorders and other cardiac pathophysiology. Cardiac mapping techniques acquire local temporal measurements and combine them to visualize the spread of electrophysiological wave phenomena across the heart surface. However, low spatial resolution, sparse measurement locations, noise and other artifacts make it challenging to accurately visualize spatio-temporal activity. For instance, electro-anatomical catheter mapping is severely limited by the sparsity of the measurements, and optical mapping is prone to noise and motion artifacts. In the past, several approaches have been proposed to create more reliable maps from noisy or sparse mapping data. Here, we demonstrate that deep learning can be used to compute phase maps and detect phase singularities in optical mapping videos of ventricular fibrillation, as well as in very noisy, low-resolution and extremely sparse simulated data of reentrant wave chaos mimicking catheter mapping data. The self-supervised deep learning approach is fundamentally different from classical phase mapping techniques. Rather than encoding a phase signal from time-series data, a deep neural network instead learns to directly associate phase maps and the positions of phase singularities with short spatio-temporal sequences of electrical data. We tested several neural network architectures, based on a convolutional neural network (CNN) with an encoding and decoding structure, to predict phase maps or rotor core positions either directly or indirectly via the prediction of phase maps and a subsequent classical calculation of phase singularities. Predictions can be performed across different data, with models being trained on one species and then successfully applied to another, or being trained solely on simulated data and then applied to experimental data. Neural networks provide a promising alternative to conventional phase mapping and rotor core localization methods. Future uses may include the analysis of optical mapping studies in basic cardiovascular research, as well as the mapping of atrial fibrillation in the clinical setting.

Atrial fibrillation driver identification through regional mutual information networks: a modeling perspective

Purpose

Effective identification of electrical drivers within remodeled tissue is a key for improving ablation treatment for atrial fibrillation. We have developed a mutual information, graph-based approach to identify and propose fault tolerance metric of local efficiency as a distinguishing feature of rotational activation and remodeled atrial tissue.

Methods

Voltage data were extracted from atrial tissue simulations (2D Karma, 3D physiological, and the Multiscale Cardiac Simulation Framework (MSCSF)) using multi-spline open and parallel regional mapping catheter geometries. Graphs were generated based on varied mutual information thresholds between electrode pairs and the local efficiency for each graph was calculated.

Results

High-resolution mapping catheter geometries can distinguish between rotational and irregular activation patterns using the derivative of local efficiency as a function of increasing mutual information threshold. The derivative is decreased for rotational activation patterns comparing to irregular activations in both a simplified 2D model (0.0017 ± 1 × 10⁻⁴ vs. 0.0032 ± 1 × 10⁻⁴, p < 0.01) and a more realistic 3D model (0.00092 ± 5 × 10⁻⁵ vs. 0.0014 ± 4 × 10⁻⁵, p < 0.01). Average local efficiency derivative can also distinguish between degrees of remodeling. Simulations using the MSCSF model, with 10 vs. 90% remodeling, display distinct derivatives in the grid design parallel spline catheter configuration (0.0015 ± 5 × 10⁻⁵ vs. 0.0019 ± 6 × 10⁻⁵, p < 0.01) and the flower shaped open spline configuration (0.0011 ± 5 × 10⁻⁵ vs. 0.0016 ± 4 × 10⁻⁵, p < 0.01).

Conclusion

A decreased derivative of local efficiency characterizes rotational activation and varies with atrial remodeling. This suggests a distinct communication pattern in cardiac rotational activation detectable via high-resolution regional mapping and could enable identification of electrical drivers for targeted ablation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10840-021-01101-z.

Evaluating spatial disparities of rotor sites and high dominant frequency regions during catheter ablation for PersAF patients targeting high dominant frequency sites using non-contacting mapping

Purpose: Several studies have emphasised the significance of high dominant frequency (HDF) and rotors in the perpetuation of AF. However, the co-localisation relationship between both attributes is not completely understood yet. In this study, we aim to evaluate the spatial distributions of HDF regions and rotor sites within the left atrium (LA) pre and post HDF-guided ablation in PersAF. Methods: This study involved 10 PersAF patients undergoing catheter ablation targeting HDF regions in the LA. 2048-channels of atrial electrograms (AEG) were collected pre- and post-ablation using a non-contact array (EnSite, Abbott). The dominant frequency (DF, 4-10 Hz) areas with DF within 0.25 Hz of the maximum out of the 2048 points were defined as "high" DF (HDF). Rotors were defined as PSs that last more than 100 ms and at a similar location through subsequent phase frames over time. Results: The results indicated an extremely poor spatial correlation between the HDF regions and sites of the rotors in pre-versus post-ablation cases for the non-terminated (pre: CORR; 0.05 ± 0.17. vs. post: CORR; -0.030 ± 0.19, and with terminated patients (pre: CORR; -0.016 ± 0.03. post: CORR; -0.022 ± 0.04). Rotors associated with AF terminations had a long-lasting life-span post-ablation (non-terminated vs. terminated 120.7 ± 6.5 ms vs. 139.9 ± 39.8 ms), high core velocity (1.35 ± 1.3 mm/ms vs. 1.32 ± 0.9 mm/ms), and were less meandering (3.4 ± 3.04 mm vs. 1.5 ± 1.2 mm). Although the results suggest a poor spatial overlapping between rotors' sites and sites of AFCL changes in terminated and non-terminated patients, a higher correlation was determined in terminated patients (spatial overlapping percentage pre: 25 ± 4.2% vs. 17 ± 3.8% vs. post: 8 ± 4.2% vs. 3.7 ± 1.7% p < 0.05, respectively). Conclusion: Using non-contact AEG, it was noted that the correlation is poor between the spatial distribution of HDF regions and sites of rotors. Rotors were longer-lasting, faster and more stationary in patients with AF termination post-ablation. Rotors sites demonstrated poor spatial overlapping with sites of AFCL changes that lead to AF termination.

Clinical Outcomes of Computational Virtual Mapping-Guided Catheter Ablation in Patients With Persistent Atrial Fibrillation: A Multicenter Prospective Randomized Clinical Trial

Background: Clinical recurrence after atrial fibrillation catheter ablation (AFCA) still remains high in patients with persistent AF (PeAF). We investigated whether an extra-pulmonary vein (PV) ablation targeting the dominant frequency (DF) extracted from electroanatomical map–integrated AF computational modeling improves the AFCA rhythm outcome in patients with PeAF.

Methods: In this open-label, randomized, multi-center, controlled trial, 170 patients with PeAF were randomized at a 1:1 ratio to the computational modeling-guided virtual DF (V-DF) ablation and empirical PV isolation (E-PVI) groups. We generated a virtual dominant frequency (DF) map based on the atrial substrate map obtained during the clinical AF ablation procedure using computational modeling. This simulation was possible within the time of the PVI procedure. V-DF group underwent extra-PV V-DF ablation in addition to PVI, but DF information was not notified to the operators from the core lab in the E-PVI group.

Results: After a mean follow-up period of 16.3 ± 5.3 months, the clinical recurrence rate was significantly lower in the V-DF than with E-PVI group (P = 0.018, log-rank). Recurrences appearing as atrial tachycardias (P = 0.145) and the cardioversion rates (P = 0.362) did not significantly differ between the groups. At the final follow-up, sinus rhythm was maintained without any AADs in 74.7% in the V-DF group and 48.2% in the E-PVI group (P < 0.001). No significant difference was found in the major complication rates (P = 0.489) or total procedure time (P = 0.513) between the groups. The V-DF ablation was independently associated with a reduced AF recurrence after AFCA [hazard ratio: 0.51 (95% confidence interval: 0.30–0.88); P = 0.016].

Conclusions: The computational modeling-guided V-DF ablation improved the rhythm outcome of AFCA in patients with PeAF.

Clinical Trial Registration: Clinical Research Information Service, CRIS identifier: KCT0003613.

Applications of multimodality imaging for left atrial catheter ablation

Atrial arrhythmias, including atrial fibrillation and atrial flutter, may be treated through catheter ablation. The process of atrial arrhythmia catheter ablation, which includes patient selection, pre-procedural planning, intra-procedural guidance, and post-procedural assessment, is typically characterized by the use of several imaging modalities to sequentially inform key clinical decisions. Increasingly, advanced imaging modalities are processed via specialized image analysis techniques and combined with intra-procedural electrical measurements to inform treatment approaches. Here, we review the use of multimodality imaging for left atrial ablation procedures. The article first outlines how imaging modalities are routinely used in the peri-ablation period. We then describe how advanced imaging techniques may inform patient selection for ablation and ablation targets themselves. Ongoing research directions for improving catheter ablation outcomes by using imaging combined with advanced analyses for personalization of ablation targets are discussed, together with approaches for their integration in the standard clinical environment. Finally, we describe future research areas with the potential to improve catheter ablation outcomes.

Computational modeling of atrial fibrillation

With the aging society, the prevalence of atrial fibrillation (AF) continues to increase. Nevertheless, there are still limitations in antiarrhythmic drugs (AAD) or catheter interventions for AF. If it is possible to predict the outcome of AF management according to various AADs or ablation lesion sets through computational modeling, it will be of great clinical help. AF computational modeling has been utilized for in-silico arrhythmia research and enabled high-density entire chamber mapping, reproducible condition control, virtual intervention, not possible clinically or experimentally, in-depth mechanistic research. With the recent development of computer science and technology, more sophisticated and faster computational modeling has become available for clinical application. In particular, it can be applied to determine the extra-PV target of persistent AF catheter ablation or to select the AAD with the best effect. AF computational modeling combined with artificial intelligence is expected to contribute to precision medicine for more diverse uses in the future. Therefore, in this review, we will deal with the history, development, and various applications of computation modeling.

Current progress of computational modeling for guiding clinical atrial fibrillation ablation

Atrial fibrillation (AF) is one of the most common arrhythmias, associated with high morbidity, mortality, and healthcare costs, and it places a significant burden on both individuals and society. Anti-arrhythmic drugs are the most commonly used strategy for treating AF. However, drug therapy faces challenges because of its limited efficacy and potential side effects. Catheter ablation is widely used as an alternative treatment for AF. Nevertheless, because the mechanism of AF is not fully understood, the recurrence rate after ablation remains high. In addition, the outcomes of ablation can vary significantly between medical institutions and patients, especially for persistent AF. Therefore, the issue of which ablation strategy is optimal is still far from settled. Computational modeling has the advantages of repeatable operation, low cost, freedom from risk, and complete control, and is a useful tool for not only predicting the results of different ablation strategies on the same model but also finding optimal personalized ablation targets for clinical reference and even guidance. This review summarizes three-dimensional computational modeling simulations of catheter ablation for AF, from the early-stage attempts such as Maze III or circumferential pulmonary vein isolation to the latest advances based on personalized substrate-guided ablation. Finally, we summarize current developments and challenges and provide our perspectives and suggestions for future directions.

Spatial Changes in the Atrial Fibrillation Wave-Dynamics After Using Antiarrhythmic Drugs: A Computational Modeling Study

Background: We previously reported that a computational modeling-guided antiarrhythmic drug (AAD) test was feasible for evaluating multiple AADs in patients with atrial fibrillation (AF). We explored the anti-AF mechanisms of AADs and spatial change in the AF wave-dynamics by a realistic computational model. Methods: We used realistic computational modeling of 25 AF patients (68% male, 59.8 ± 9.8 years old, 32.0% paroxysmal AF) reflecting the anatomy, histology, and electrophysiology of the left atrium (LA) to characterize the effects of five AADs (amiodarone, sotalol, dronedarone, flecainide, and propafenone). We evaluated the spatial change in the AF wave-dynamics by measuring the mean dominant frequency (DF) and its coefficient of variation [dominant frequency-coefficient of variation (DF-COV)] in 10 segments of the LA. The mean DF and DF-COV were compared according to the pulmonary vein (PV) vs. extra-PV, maximal slope of the restitution curves (Smax), and defragmentation of AF. Results: The mean DF decreased after the administration of AADs in the dose dependent manner (p < 0.001). Under AADs, the DF was significantly lower (p < 0.001) and COV-DF higher (p = 0.003) in the PV than extra-PV region. The mean DF was significantly lower at a high Smax (≥1.4) than a lower Smax condition under AADs. During the episodes of AF defragmentation, the mean DF was lower (p < 0.001), but the COV-DF was higher (p < 0.001) than that in those without defragmentation. Conclusions: The DF reduction with AADs is predominant in the PVs and during a high Smax condition and causes AF termination or defragmentation during a lower DF and spatially unstable (higher DF-COV) condition.

Prediction of Cardiac Mechanical Performance From Electrical Features During Ventricular Tachyarrhythmia Simulation Using Machine Learning Algorithms

In ventricular tachyarrhythmia, electrical instability features including action potential duration, dominant frequency, phase singularity, and filaments are associated with mechanical contractility. However, there are insufficient studies on estimated mechanical contractility based on electrical features during ventricular tachyarrhythmia using a stochastic model. In this study, we predicted cardiac mechanical performance from features of electrical instability during ventricular tachyarrhythmia simulation using machine learning algorithms, including support vector regression (SVR) and artificial neural network (ANN) models. We performed an electromechanical tachyarrhythmia simulation and extracted 12 electrical instability features and two mechanical properties, including stroke volume and the amplitude of myocardial tension (ampTens). We compared predictive performance according to kernel types of the SVR model and the number of hidden layers of the ANN model. In the SVR model, the prediction accuracies of stroke volume and ampTens were the highest when using the polynomial kernel and linear kernel, respectively. The predictive performance of the ANN model was better than that of the SVR model. The prediction accuracies were the highest when the ANN model consisted of three hidden layers. Accordingly, we propose the ANN model with three hidden layers as an optimal model for predicting cardiac mechanical contractility in ventricular tachyarrhythmia. The results of this study are expected to be used to indirectly estimate the hemodynamic response from the electrical cardiac map measured by the optical mapping system during cardiac surgery, as well as cardiac contractility under normal sinus rhythm conditions.

Canine atrial fibrillation: Pathophysiology, epidemiology and classification

Atrial fibrillation (AF) is the most common non-physiological arrhythmia in dogs and humans. Its high prevalence in both species and the impact it has on survival time and quality of life of affected patients, makes it a very relevant topic of medical research. Significant developments in understanding the mechanisms underlying this arrhythmia in humans has occurred over the last decades and some of this knowledge is being applied to veterinary medicine, despite the many differences between species. This article reviews the current understanding of the pathophysiology of AF. The epidemiology and classification of AF in dogs will also be discussed.

In situ procedure for high-efficiency computational modeling of atrial fibrillation reflecting personal anatomy, fiber orientation, fibrosis, and electrophysiology

We previously reported the feasibility and efficacy of a simulation-guided clinical catheter ablation of atrial fibrillation (AF) in an in-silico AF model. We developed a highly efficient realistic AF model reflecting the patient endocardial voltage and local conduction and tested its clinical feasibility. We acquired > 500 endocardial bipolar electrograms during right atrial pacing at the beginning of the AF ablation procedures. Based on the clinical bipolar electrograms, we generated simulated voltage maps by applying fibrosis and local activation maps adjusted for the fiber orientation. The software's accuracy (CUVIA2.5) was retrospectively tested in 17 patients and feasibility prospectively in 10 during clinical AF ablation. Results: We found excellent correlations between the clinical and simulated voltage maps (R = 0.933, p < 0.001) and clinical and virtual local conduction (R = 0.958, p < 0.001). The proportion of virtual local fibrosis was 15.4, 22.2, and 36.9% in the paroxysmal AF, persistent AF, and post-pulmonary vein isolation (PVI) states, respectively. The reconstructed virtual bipolar electrogram exhibited a relatively good similarities of morphology to the local clinical bipolar electrogram (R = 0.60 ± 0.08, p < 0.001). Feasibility testing revealed an in situ procedural computing time from the clinical data acquisition to wave-dynamics analyses of 48.2 ± 4.9 min. All virtual analyses were successfully achieved during clinical PVI procedures. We developed a highly efficient, realistic, in situ procedural simulation model reflective of individual anatomy, fiber orientation, fibrosis, and electrophysiology that can be applied during AF ablation.

Multiple factors influence the morphology of the bipolar electrogram: An in silico modeling study

Although bipolar electrograms (Bi-egms) are commonly used for catheter mapping and ablation of cardiac arrhythmias, the accuracy and reproducibility of Bi-egms have not been evaluated. We aimed to clarify the influence of the catheter orientation (CO), catheter contact angle (CA), local conduction velocity (CV), scar size, and catheter type on the Bi-egm morphology using an in silico 3-dimensional realistic model of atrial fibrillation. We constructed a 3-dimensional, realistic, in silico left atrial model with activation wave propagation including bipolar catheter models. Bi-egms were obtained by computing the extracellular potentials from the distal and proximal electrodes. The amplitude and width were measured on virtual Bi-egms obtained under different conditions created by changing the CO according to the wave direction, catheter-atrial wall CA, local CV, size of the non-conductive area, and catheter type. Bipolar voltages were also compared between virtual and clinically acquired Bi-egms. Bi-egm amplitudes were lower for a perpendicular than parallel CO relative to the wave direction (p<0.001), lower for a 90° than 0° CA (p<0.001), and lower for a CV of 0.13m/s than 0.48m/s (p<0.001). Larger sized non-conductive areas were associated with a decreased bipolar amplitude (p<0.001) and increased bipolar width (p<0.001). Among three commercially available catheters (Orion, Pentaray, and Thermocool), those with more narrowly spaced and smaller electrodes produced higher voltages on the virtual Bi-egms (p<0.001). Multiple factors including the CO, CA, CV, and catheter design significantly influence the Bi-egm morphology. Universal voltage cut-off values may not be appropriate for bipolar voltage-guided substrate mapping.

Cardiac arrhythmias are rhythm disorders of the heart leading to abnormal heart function. For the diagnosis and treatment of the arrhythmias, clinicians insert catheters into the heart and examine the electrical signal propagation in the heart. Among different type of catheters, bipolar catheters have two electrodes at the tip of the catheter with the signal being the difference between the two electrodes, which provides sharper signal than unipolar catheter. However, bipolar electrogram is dependent on many factors including catheter design and orientation, and consequently, knowledge of the determinants of the bipolar electrogram is needed for proper interpretation of the signal. In this study, we examined the effects of many factors on bipolar electrogram using computer simulation. Computer simulation is very useful in this type of study because, in clinical settings, it is not feasible to control each factor precisely. We quantitatively demonstrated the effects of catheter design and orientation, and cardiac wave propagation speed on bipolar electrogram.

Entropy Mapping Approach for Functional Reentry Detection in Atrial Fibrillation: An In-Silico Study

Catheter ablation of critical electrical propagation sites is a promising tool for reducing the recurrence of atrial fibrillation (AF). The spatial identification of the arrhythmogenic mechanisms sustaining AF requires the evaluation of electrograms (EGMs) recorded over the atrial surface. This work aims to characterize functional reentries using measures of entropy to track and detect a reentry core. To this end, different AF episodes are simulated using a 2D model of atrial tissue. Modified Courtemanche human action potential and Fenton–Karma models are implemented. Action potential propagation is modeled by a fractional diffusion equation, and virtual unipolar EGM are calculated. Episodes with stable and meandering rotors, figure-of-eight reentry, and disorganized propagation with multiple reentries are generated. Shannon entropy (ShEn), approximate entropy (ApEn), and sample entropy (SampEn) are computed from the virtual EGM, and entropy maps are built. Phase singularity maps are implemented as references. The results show that ApEn and SampEn maps are able to detect and track the reentry core of rotors and figure-of-eight reentry, while the ShEn results are not satisfactory. Moreover, ApEn and SampEn consistently highlight a reentry core by high entropy values for all of the studied cases, while the ability of ShEn to characterize the reentry core depends on the propagation dynamics. Such features make the ApEn and SampEn maps attractive tools for the study of AF reentries that persist for a period of time that is similar to the length of the observation window, and reentries could be interpreted as AF-sustaining mechanisms. Further research is needed to determine and fully understand the relation of these entropy measures with fibrillation mechanisms other than reentries.

Clinical Usefulness of Computational Modeling-Guided Persistent Atrial Fibrillation Ablation: Updated Outcome of Multicenter Randomized Study

OBJECTIVE

Catheter ablation of persistent atrial fibrillation (AF) is still challenging, no optimal extra-pulmonary vein lesion set is known. We previously reported the clinical feasibility of computational modeling-guided AF catheter ablation.

METHODS

We randomly assigned 118 patients with persistent AF (77.8% men, age 60.8 ± 9.9 years) to the computational modeling-guided ablation group (53 patients) and the empirical ablation group (55 patients) based on the operators' experience. For virtual ablation, four virtual linear and one electrogram-guided lesion sets were tested on patient heart computed tomogram-based models, and the lesion set with the fastest termination time was reported to the operator in the modeling-guided ablation group. The primary outcome was freedom from atrial tachyarrhythmias lasting longer than 30 s after a single procedure.

RESULTS

During 31.5 ± 9.4 months, virtual ablation procedures were available in 95.2% of the patients (108/118). Clinical recurrence rate was significantly lower after a modeling-guided ablation than after an empirical ablation (20.8 vs. 40.0%, log-rank p = 0.042). Modeling-guided ablation was independently associated with a better long-term rhythm outcome of persistent AF ablation (HR = 0.29 [0.12-0.69], p = 0.005). The rhythm outcome of the modeling-guided ablation showed better trends in males, non-obese patients with a less remodeled atrium (left atrial dimension < 50 mm), ejection fraction ≥ 50%, and those without hypertension or diabetes (p < 0.01). There were no significant differences between the groups for the total procedure time (p = 0.403), ablation time (p = 0.510), and major complication rate (p = 0.900).

CONCLUSION

Among patients with persistent AF, the computational modeling-guided ablation was superior to the empirical catheter ablation regarding the rhythm outcome.

CLINICAL TRIAL REGISTRATION

This study was registered with the ClinicalTrials.gov, number NCT02171364.

Analytical approaches for myocardial fibrillation signals

Atrial and ventricular fibrillation are complex arrhythmias, and their underlying mechanisms remain widely debated and incompletely understood. This is partly because the electrical signals recorded during myocardial fibrillation are themselves complex and difficult to interpret with simple analytical tools. There are currently a number of analytical approaches to handle fibrillation data. Some of these techniques focus on mapping putative drivers of myocardial fibrillation, such as dominant frequency, organizational index, Shannon entropy and phase mapping. Other techniques focus on mapping the underlying myocardial substrate sustaining fibrillation, such as voltage mapping and complex fractionated electrogram mapping. In this review, we discuss these techniques, their application and their limitations, with reference to our experimental and clinical data. We also describe novel tools including a new algorithm to map microreentrant circuits sustaining fibrillation.

A Computational Framework to Benchmark Basket Catheter Guided Ablation in Atrial Fibrillation

Catheter ablation is a curative therapeutic approach for atrial fibrillation (AF). Ablation of rotational sources based on basket catheter measurements has been proposed as a promising approach in patients with persistent AF to complement pulmonary vein isolation. However, clinically reported success rates are equivocal calling for a mechanistic investigation under controlled conditions. We present a computational framework to benchmark ablation strategies considering the whole cycle from excitation propagation to electrogram acquisition and processing to virtual therapy. Fibrillation was induced in a patient-specific 3D volumetric model of the left atrium, which was homogeneously remodeled to sustain reentry. The resulting extracellular potential field was sampled using models of grid catheters as well as realistically deformed basket catheters considering the specific atrial anatomy. The virtual electrograms were processed to compute phase singularity density maps to target rotor tips with up to three circular ablations. Stable rotors were successfully induced in different regions of the homogeneously remodeled atrium showing that rotors are not constrained to unique anatomical structures or locations. Density maps of rotor tip trajectories correctly identified and located the rotors (deviation < 10 mm) based on catheter recordings only for sufficient resolution (inter-electrode distance ≤3 mm) and proximity to the wall (≤10 mm). Targeting rotor sites with ablation did not stop reentries in the homogeneously remodeled atria independent from lesion size (1-7 mm radius), from linearly connecting lesions with anatomical obstacles, and from the number of rotors targeted sequentially (≤3). Our results show that phase maps derived from intracardiac electrograms can be a powerful tool to map atrial activation patterns, yet they can also be misleading due to inaccurate localization of the rotor tip depending on electrode resolution and distance to the wall. This should be considered to avoid ablating regions that are in fact free of rotor sources of AF. In our experience, ablation of rotor sites was not successful to stop fibrillation. Our comprehensive simulation framework provides the means to holistically benchmark ablation strategies in silico under consideration of all steps involved in electrogram-based therapy and, in future, could be used to study more heterogeneously remodeled disease states as well.

Targeting the Substrate in Ablation of Persistent Atrial Fibrillation: Recent Lessons and Future Directions

While isolation of the pulmonary veins is firmly established as effective treatment for the majority of paroxysmal atrial fibrillation (AF) patients, there is recognition that patients with persistent AF have substrate for perpetuation of arrhythmia existing outside of the pulmonary veins. Various computational approaches have been used to identify targets for effective ablation of persistent AF. This paper aims to discuss the clinical aspects of computational approaches that aim to identify critical sites for treatment. Various analyses of electrogram characteristics have been performed with this aim. Leading techniques for electrogram analysis are Complex Fractionated Atrial Electrograms (CFAE) and Dominant Frequency (DF). These techniques have been the subject of clinical trials of which the results are discussed. Evaluation of the activation patterns of atria in AF has been another avenue of research. Focal Impulse and Rotor Modulation (FIRM) mapping and forms of Body Surface Mapping aim to characterize multiple atrial wavelets, macro-reentry and focal sources which have been proposed as basic mechanisms perpetuating AF. Both invasive and non-invasive activation mapping techniques are reviewed. The presence of atrial fibrosis causes non-uniform anisotropic impulse propagation. Therefore, identification of fibrosis by imaging techniques is an avenue of potential research. The leading contender for imaging-based techniques is Cardiac Magnetic Resonance (CMR). As this technology advances, improvements in resolution and scar identification have positioned CMR as the mode of choice for analysis of atrial structure. AF has been demonstrated to be associated with obesity, inactivity and diseases of modern life. An opportunity exists for detailed computational analysis of the impact of risk factor modification on atrial substrate. This ranges from microstructural investigation through to examination at a population level via registries and public health interventions. Computational analysis of atrial substrate has moved from basic science toward clinical application. Future directions and potential limitations of such analyses are examined in this review.

Information Theory and Atrial Fibrillation (AF): A Review

Atrial Fibrillation (AF) is the most common cardiac rhythm disorder seen in hospitals and in general practice, accounting for up to a third of arrhythmia related hospitalizations. Unfortunately, AF treatment is in practice complicated by the lack of understanding of the fundamental mechanisms underlying the arrhythmia, which makes detection of effective ablation targets particularly difficult. Various approaches to AF mapping have been explored in the hopes of better pinpointing these effective targets, such as Dominant Frequency (DF) analysis, complex fractionated electrograms (CFAE) and unipolar reconstruction (FIRM), but many of these methods have produced conflicting results or require further investigation. Exploration of AF using information theoretic-based approaches may have the potential to provide new insights into the complex system dynamics of AF, whilst also providing the benefit of being less reliant on empirically derived definitions in comparison to alternate mapping approaches. This work provides an overview of information theory and reviews its applications in AF analysis, with particular focus on AF mapping. The works discussed in this review demonstrate how understanding AF from a signal property perspective can provide new insights into the arrhythmic phenomena, which may have valuable clinical implications for AF mapping and ablation in the future.

Effectiveness of atrial fibrillation rotor ablation is dependent on conduction velocity: An in-silico 3-dimensional modeling study

BACKGROUND

We previously reported that stable rotors are observed in in-silico human atrial fibrillation (AF) models, and are well represented by a dominant frequency (DF). In the current study, we hypothesized that the outcome of DF ablation is affected by conduction velocity (CV) conditions and examined this hypothesis using in-silico 3D-AF modeling.

METHODS

We integrated 3D CT images of left atrium obtained from 10 patients with persistent AF (80% male, 61.8±13.5 years old) into in-silico AF model. We compared AF maintenance durations (max 300s), spatiotemporal stabilities of DF, phase singularity (PS) number, life-span of PS, and AF termination or defragmentation rates after virtual DF ablation with 5 different CV conditions (0.2, 0.3, 0.4, 0.5, and 0.6m/s).

RESULTS

1. AF maintenance duration (p<0.001), spatiotemporal mean variance of DF (p<0.001), and the number of PS (p = 0.023) showed CV dependent bimodal patterns (highest at CV0.4m/s and lowest at CV0.6m/s) consistently. 2. After 10% highest DF ablation, AF defragmentation rates were the lowest at CV0.4m/s (37.8%), but highest at CV0.5 and 0.6m/s (all 100%, p<0.001). 3. In the episodes with AF termination or defragmentation followed by 10% highest DF ablation, baseline AF maintenance duration was shorter (p<0.001), spatiotemporal mean variance of DF was lower (p = 0.014), and the number of PS was lower (p = 0.004) than those with failed AF defragmentation after DF ablation.

CONCLUSION

Virtual ablation of DF, which may indicate AF driver, was more likely to terminate or defragment AF with spatiotemporally stable DF, but not likely to do so in long-lasting and sustained AF conditions, depending on CV.

Impact of number of co-existing rotors and inter-electrode distance on accuracy of rotor localization

BACKGROUND

Conflicting evidence exists on the efficacy of focal impulse and rotor modulation on atrial fibrillation ablation. A potential explanation is inaccurate rotor localization from multiple rotors coexistence and a relatively large (9-11mm) inter-electrode distance (IED) of the multi-electrode basket catheter.

METHODS AND RESULTS

We studied a numerical model of cardiac action potential to reproduce one through seven rotors in a two-dimensional lattice. We estimated rotor location using phase singularity, Shannon entropy and dominant frequency. We then spatially downsampled the time series to create IEDs of 2-30mm. The error of rotor localization was measured with reference to the dynamics of phase singularity at the original spatial resolution (IED=1mm). IED has a significant impact on the error using all the methods. When only one rotor is present, the error increases exponentially as a function of IED. At the clinical IED of 10mm, the error is 3.8mm (phase singularity), 3.7mm (dominant frequency), and 11.8mm (Shannon entropy). When there are more than one rotors, the error of rotor localization increases 10-fold. The error based on the phase singularity method at the clinical IED of 10mm ranges from 30.0mm (two rotors) to 96.1mm (five rotors).

CONCLUSIONS

The magnitude of error of rotor localization using a clinically available basket catheter, in the presence of multiple rotors might be high enough to impact the accuracy of targeting during AF ablation. Improvement of catheter design and development of high-density mapping catheters may improve clinical outcomes of FIRM-guided AF ablation.

A New Efficient Method for Detecting Phase Singularity in Cardiac Fibrillation

Background

The point of phase singularity (PS) is considered to represent a spiral wave core or a rotor in cardiac fibrillation. Computational efficiency is important for detection of PS in clinical electrophysiology. We developed a novel algorithm for highly efficient and robust detection of PS.

Methods

In contrast to the conventional method, which calculates PS based on the line integral of the phase around a PS point equal to ±2π (the Iyer-Gray method), the proposed algorithm (the location-centric method) looks for the phase discontinuity point at which PS actually occurs. We tested the efficiency and robustness of these two methods in a two-dimensional mathematical model of atrial fibrillation (AF), with and without remodeling of ionic currents.

Results

1. There was a significant association, in terms of the Hausdorff distance (3.30 ± 0.0 mm), between the PS points measured using the Iyer-Gray and location-centric methods, with almost identical PS trajectories generated by the two methods. 2. For the condition of electrical remodeling of AF (0.3 × I_CaL), the PS points calculated by the two methods were satisfactorily co-localized (with the Hausdorff distance of 1.64 ± 0.09 mm). 3. The proposed location-centric method was substantially more efficient than the Iyer-Gray method, with a 28.6-fold and 28.2-fold shorter run times for the control and remodeling scenarios, respectively.

Conclusion

We propose a new location-centric method for calculating PS, which is robust and more efficient compared with the conventionally used method.

The Spatiotemporal Stability of Dominant Frequency Sites in In-Silico Modeling of 3-Dimensional Left Atrial Mapping of Atrial Fibrillation

BACKGROUND

We previously reported that stable rotors were observed in in-silico human atrial fibrillation (AF) models, and were well represented by dominant frequency (DF). We explored the spatiotemporal stability of DF sites in 3D-AF models imported from patient CT images of the left atrium (LA).

METHODS

We integrated 3-D CT images of the LA obtained from ten patients with persistent AF (male 80%, 61.8 ± 13.5 years old) into an in-silico AF model. After induction, we obtained 6 seconds of AF simulation data for DF analyses in 30 second intervals (T1-T9). The LA was divided into ten sections. Spatiotemporal changes and variations in the temporal consistency of DF were evaluated at each section of the LA. The high DF area was defined as the area with the highest 10% DF.

RESULTS

1. There was no spatial consistency in the high DF distribution at each LA section during T1-T9 except in one patient (p = 0.027). 2. Coefficients of variation for the high DF area were highly different among the ten LA sections (p < 0.001), and they were significantly higher in the four pulmonary vein (PV) areas, the LA appendage, and the peri-mitral area than in the other LA sections (p < 0.001). 3. When we conducted virtual ablation of 10%, 15%, and 20% of the highest DF areas (n = 270 cases), AF was changed to atrial tachycardia (AT) or terminated at a rate of 40%, 57%, and 76%, respectively.

CONCLUSIONS

Spatiotemporal consistency of the DF area was observed in 10% of AF patients, and high DF areas were temporally variable. Virtual ablation of DF is moderately effective in AF termination and AF changing into AT.