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.

A Data-Driven Preprocessing Framework for Atrial Fibrillation Intracardiac Electrocardiogram Analysis

Atrial Fibrillation (AF) is the most common cardiac arrhythmia. Signal-processing approaches are widely used for the analysis of intracardiac electrograms (iEGMs), which are collected during catheter ablation from patients with AF. In order to identify possible targets for ablation therapy, dominant frequency (DF) is widely used and incorporated in electroanatomical mapping systems. Recently, a more robust measure, multiscale frequency (MSF), for iEGM data analysis was adopted and validated. However, before completing any iEGM analysis, a suitable bandpass (BP) filter must be applied to remove noise. Currently, no clear guidelines for BP filter characteristics exist. The lower bound of the BP filter is usually set to 3-5 Hz, while the upper bound (BP¯th) of the BP filter varies from 15 Hz to 50 Hz according to many researchers. This large range of BP¯th subsequently affects the efficiency of further analysis. In this paper, we aimed to develop a data-driven preprocessing framework for iEGM analysis, and validate it based on DF and MSF techniques. To achieve this goal, we optimized the BP¯th using a data-driven approach (DBSCAN clustering) and demonstrated the effects of different BP¯th on subsequent DF and MSF analysis of clinically recorded iEGMs from patients with AF. Our results demonstrated that our preprocessing framework with BP¯th = 15 Hz has the best performance in terms of the highest Dunn index. We further demonstrated that the removal of noisy and contact-loss leads is necessary for performing correct data iEGMs data analysis.

Atrial Flutter Mechanism Detection Using Directed Network Mapping

Atrial flutter (AFL) is a common atrial arrhythmia typically characterized by electrical activity propagating around specific anatomical regions. It is usually treated with catheter ablation. However, the identification of rotational activities is not straightforward, and requires an intense effort during the first phase of the electrophysiological (EP) study, i.e., the mapping phase, in which an anatomical 3D model is built and electrograms (EGMs) are recorded. In this study, we modeled the electrical propagation pattern of AFL (measured during mapping) using network theory (NT), a well-known field of research from the computer science domain. The main advantage of NT is the large number of available algorithms that can efficiently analyze the network. Using directed network mapping, we employed a cycle-finding algorithm to detect all cycles in the network, resembling the main propagation pattern of AFL. The method was tested on two subjects in sinus rhythm, six in an experimental model of in-silico simulations, and 10 subjects diagnosed with AFL who underwent a catheter ablation. The algorithm correctly detected the electrical propagation of both sinus rhythm cases and in-silico simulations. Regarding the AFL cases, arrhythmia mechanisms were either totally or partially identified in most of the cases (8 out of 10), i.e., cycles around the mitral valve, tricuspid valve and figure-of-eight reentries. The other two cases presented a poor mapping quality or a major complexity related to previous ablations, large areas of fibrotic tissue, etc. Directed network mapping represents an innovative tool that showed promising results in identifying AFL mechanisms in an automatic fashion. Further investigations are needed to assess the reliability of the method in different clinical scenarios.

The openCARP simulation environment for cardiac electrophysiology

BACKGROUND AND OBJECTIVE

Cardiac electrophysiology is a medical specialty with a long and rich tradition of computational modeling. Nevertheless, no community standard for cardiac electrophysiology simulation software has evolved yet. Here, we present the openCARP simulation environment as one solution that could foster the needs of large parts of this community.

METHODS AND RESULTS

openCARP and the Python-based carputils framework allow developing and sharing simulation pipelines which automate in silico experiments including all modeling and simulation steps to increase reproducibility and productivity. The continuously expanding openCARP user community is supported by tailored infrastructure. Documentation and training material facilitate access to this complementary research tool for new users. After a brief historic review, this paper summarizes requirements for a high-usability electrophysiology simulator and describes how openCARP fulfills them. We introduce the openCARP modeling workflow in a multi-scale example of atrial fibrillation simulations on single cell, tissue, organ and body level and finally outline future development potential.

CONCLUSION

As an open simulator, openCARP can advance the computational cardiac electrophysiology field by making state-of-the-art simulations accessible. In combination with the carputils framework, it offers a tailored software solution for the scientific community and contributes towards increasing use, transparency, standardization and reproducibility of in silico experiments.

Cycle length statistics during human atrial fibrillation reveal refractory properties of the underlying substrate: a combined in silico and clinical test of concept study

AIMS

The treatment of atrial fibrillation beyond pulmonary vein isolation has remained an unsolved challenge. Targeting regions identified by different substrate mapping approaches for ablation resulted in ambiguous outcomes. With the effective refractory period being a fundamental prerequisite for the maintenance of fibrillatory conduction, this study aims at estimating the effective refractory period with clinically available measurements.

METHODS AND RESULTS

A set of 240 simulations in a spherical model of the left atrium with varying model initialization, combination of cellular refractory properties, and size of a region of lowered effective refractory period was implemented to analyse the capabilities and limitations of cycle length mapping. The minimum observed cycle length and the 25% quantile were compared to the underlying effective refractory period. The density of phase singularities was used as a measure for the complexity of the excitation pattern. Finally, we employed the method in a clinical test of concept including five patients. Areas of lowered effective refractory period could be distinguished from their surroundings in simulated scenarios with successfully induced multi-wavelet re-entry. Larger areas and higher gradients in effective refractory period as well as complex activation patterns favour the method. The 25% quantile of cycle lengths in patients with persistent atrial fibrillation was found to range from 85 to 190 ms.

CONCLUSION

Cycle length mapping is capable of highlighting regions of pathologic refractory properties. In combination with complementary substrate mapping approaches, the method fosters confidence to enhance the treatment of atrial fibrillation beyond pulmonary vein isolation particularly in patients with complex activation patterns.

Using Machine Learning to Characterize Atrial Fibrotic Substrate From Intracardiac Signals With a Hybrid in silico and in vivo Dataset

In patients with atrial fibrillation, intracardiac electrogram signal amplitude is known to decrease with increased structural tissue remodeling, referred to as fibrosis. In addition to the isolation of the pulmonary veins, fibrotic sites are considered a suitable target for catheter ablation. However, it remains an open challenge to find fibrotic areas and to differentiate their density and transmurality. This study aims to identify the volume fraction and transmurality of fibrosis in the atrial substrate. Simulated cardiac electrograms, combined with a generalized model of clinical noise, reproduce clinically measured signals. Our hybrid dataset approach combines in silico and clinical electrograms to train a decision tree classifier to characterize the fibrotic atrial substrate. This approach captures different in vivo dynamics of the electrical propagation reflected on healthy electrogram morphology and synergistically combines it with synthetic fibrotic electrograms from in silico experiments. The machine learning algorithm was tested on five patients and compared against clinical voltage maps as a proof of concept, distinguishing non-fibrotic from fibrotic tissue and characterizing the patient's fibrotic tissue in terms of density and transmurality. The proposed approach can be used to overcome a single voltage cut-off value to identify fibrotic tissue and guide ablation targeting fibrotic areas.

Comparison of Unipolar and Bipolar Voltage Mapping for Localization of Left Atrial Arrhythmogenic Substrate in Patients With Atrial Fibrillation

Background: Presence of left atrial low voltage substrate in bipolar voltage mapping is associated with increased arrhythmia recurrences following pulmonary vein isolation for atrial fibrillation (AF). Besides local myocardial fibrosis, bipolar voltage amplitudes may be influenced by inter-electrode spacing and bipole-to-wavefront-angle. It is unclear to what extent these impact low voltage areas (LVA) in the clinical setting. Alternatively, unipolar electrogram voltage is not affected by these factors but requires advanced filtering. Objectives: To assess the relationship between bipolar and unipolar voltage mapping in sinus rhythm (SR) and AF and identify if the electrogram recording mode affects the quantification and localization of LVA. Methods: Patients (n = 28, 66±7 years, 46% male, 82% persistent AF, 32% redo-procedures) underwent high-density (>1,200 sites, 20 ± 10 sites/cm², using a 20-pole 2-6-2 mm-spaced Lasso) voltage mapping in SR and AF. Bipolar LVA were defined using four different thresholds described in literature: <0.5 and <1 mV in SR, <0.35 and <0.5 mV in AF. The optimal unipolar voltage threshold resulting in the highest agreement in both unipolar and bipolar mapping modes was determined. The impact of the inter-electrode distance (2 vs. 6 mm) on the correlation was assessed. Regional analysis was performed using an 11-segment left atrial model. Results: Patients had relevant bipolar LVA (23 ± 23 cm² at <0.5 mV in SR and 42 ± 26 cm² at <0.5 mV in AF). 90 ± 5% (in SR) and 85 ± 5% (AF) of mapped sites were concordantly classified as high or low voltage in both mapping modes. Discordant mapping sites located to the border zone of LVA. Bipolar voltage mapping using 2 vs. 6 mm inter-electrode distances increased the portion of matched mapping points by 4%. The unipolar thresholds (y) which resulted in a high spatial concordance can be calculated from the bipolar threshold (x) using following linear equations: y = 1.06x + 0.26mV (r = 0.994) for SR and y = 1.22x + 0.12mV (r = 0.998) for AF. Conclusion: Bipolar and unipolar voltage maps are highly correlated, in SR and AF. While bipole orientation and inter-electrode spacing are theoretical confounders, their impact is unlikely to be of clinical importance for localization of LVA, when mapping is performed at high density with a 20-polar Lasso catheter.