Frequency domain algorithm for quantifying atrial fibrillation organization to increase defibrillation efficacy

Systematic in-silico evaluation of fibrosis effects on re-entrant wave dynamics in atrial tissue

Despite the key role of fibrosis in atrial fibrillation (AF), the effects of different spatial distributions and textures of fibrosis on wave propagation mechanisms in AF are not fully understood. To clarify these aspects, we performed a systematic computational study to assess fibrosis effects on the characteristics and stability of re-entrant waves in electrically-remodelled atrial tissues. A stochastic algorithm, which generated fibrotic distributions with controlled overall amount, average size, and orientation of fibrosis elements, was implemented on a monolayer spheric atrial model. 245 simulations were run at changing fibrosis parameters. The emerging propagation patterns were quantified in terms of rate, regularity, and coupling by frequency-domain analysis of correspondent synthetic bipolar electrograms. At the increase of fibrosis amount, the rate of reentrant waves significantly decreased and higher levels of regularity and coupling were observed (p < 0.0001). Higher spatial variability and pattern stochasticity over repetitions was observed for larger amount of fibrosis, especially in the presence of patchy and compact fibrosis. Overall, propagation slowing and organization led to higher stability of re-entrant waves. These results strengthen the evidence that the amount and spatial distribution of fibrosis concur in dictating re-entry dynamics in remodeled tissue and represent key factors in AF maintenance.

Novel Entropy-Based Metrics for Long-Term Atrial Fibrillation Recurrence Prediction Following Surgical Ablation: Insights from Preoperative Electrocardiographic Analysis

Atrial fibrillation (AF) is a prevalent cardiac arrhythmia often treated concomitantly with other cardiac interventions through the Cox-Maze procedure. This highly invasive intervention is still linked to a long-term recurrence rate of approximately 35% in permanent AF patients. The aim of this study is to preoperatively predict long-term AF recurrence post-surgery through the analysis of atrial activity (AA) organization from non-invasive electrocardiographic (ECG) recordings. A dataset comprising ECGs from 53 patients with permanent AF who had undergone Cox-Maze concomitant surgery was analyzed. The AA was extracted from the lead V1 of these recordings and then characterized using novel predictors, such as the mean and standard deviation of the relative wavelet energy (RWEm and RWEs) across different scales, and an entropy-based metric that computes the stationary wavelet entropy variability (SWEnV). The individual predictors exhibited limited predictive capabilities to anticipate the outcome of the procedure, with the SWEnV yielding a classification accuracy (Acc) of 68.07%. However, the assessment of the RWEs for the seventh scale (RWEs7), which encompassed frequencies associated with the AA, stood out as the most promising individual predictor, with sensitivity (Se) and specificity (Sp) values of 80.83% and 67.09%, respectively, and an Acc of almost 75%. Diverse multivariate decision tree-based models were constructed for prediction, giving priority to simplicity in the interpretation of the forecasting methodology. In fact, the combination of the SWEnV and RWEs7 consistently outperformed the individual predictors and excelled in predicting post-surgery outcomes one year after the Cox-Maze procedure, with Se, Sp, and Acc values of approximately 80%, thus surpassing the results of previous studies based on anatomical predictors associated with atrial function or clinical data. These findings emphasize the crucial role of preoperative patient-specific ECG signal analysis in tailoring post-surgical care, enhancing clinical decision making, and improving long-term clinical outcomes.

Vector Field Heterogeneity for the Assessment of Locally Disorganised Cardiac Electrical Propagation Wavefronts From High-Density Multielectrodes

High-density multielectrode catheters are becoming increasingly popular in cardiac electrophysiology for advanced characterisation of the cardiac tissue, due to their potential to identify impaired sites. These are often characterised by abnormal electrical conduction, which may cause locally disorganised propagation wavefronts. To quantify it, a novel heterogeneity parameter based on vector field analysis is proposed, utilising finite differences to measure direction changes between adjacent cliques. The proposed Vector Field Heterogeneity metric has been evaluated on a set of simulations with controlled levels of organisation in vector maps, and a variety of grid sizes. Furthermore, it has been tested on animal experimental models of isolated Langendorff-perfused rabbit hearts. The proposed parameter exhibited superior capturing ability of heterogeneous propagation wavefronts compared to the classical Spatial Inhomogeneity Index, and simulations proved that the metric effectively captures gradual increments in disorganisation in propagation patterns. Notably, it yielded robust and consistent outcomes for [Formula: see text] grid sizes, underscoring its suitability for the latest generation of orientation-independent cardiac catheters.

A 3D-CNN with temporal-attention block to predict the recurrence of atrial fibrillation based on body-surface potential mapping signals

Catheter ablation has become an important treatment for atrial fibrillation (AF), but its recurrence rate is still high. The aim of this study was to predict AF recurrence using a three-dimensional (3D) network model based on body-surface potential mapping signals (BSPMs). BSPMs were recorded with a 128-lead vest in 14 persistent AF patients before undergoing catheter ablation (Maze-IV). The torso geometry was acquired and meshed by point cloud technology, and the BSPM was interpolated into the torso geometry by the inverse distance weighted (IDW) method to generate the isopotential map. Experiments show that the isopotential map of BSPMs can reflect the propagation of the electrical wavefronts. The 3D isopotential sequence map was established by combining the spatial–temporal information of the isopotential map; a 3D convolutional neural network (3D-CNN) model with temporal attention was established to predict AF recurrence. Our study proposes a novel attention block that focuses the characteristics of atrial activations to improve sampling accuracy. In our experiment, accuracy (ACC) in the intra-patient evaluation for predicting the recurrence of AF was 99.38%. In the inter-patient evaluation, ACC of 3D-CNN was 81.48%, and the area under the curve (AUC) was 0.88. It can be concluded that the dynamic rendering of multiple isopotential maps can not only comprehensively display the conduction of cardiac electrical activity on the body surface but also successfully predict the recurrence of AF after CA by using 3D isopotential sequence maps.

Preoperative Prediction of Catheter Ablation Outcome in Persistent Atrial Fibrillation Patients through Spectral Organization Analysis of the Surface Fibrillatory Waves

Catheter ablation (CA) is a commonly used treatment for persistent atrial fibrillation (AF). Since its medium/long-term success rate remains limited, preoperative prediction of its outcome is gaining clinical interest to optimally select candidates for the procedure. Among predictors based on the surface electrocardiogram, the dominant frequency (DF) and harmonic exponential decay (γ) of the fibrillatory waves (f-waves) have reported promising but clinically insufficient results. Hence, the main goal of this work was to conduct a broader analysis of the f-wave harmonic spectral structure to improve CA outcome prediction through several entropy-based measures computed on different frequency bands. On a database of 151 persistent AF patients under radio-frequency CA and a follow-up of 9 months, the newly introduced parameters discriminated between patients who relapsed to AF and those who maintained SR at about 70%, which was statistically superior to the DF and approximately similar to γ. They also provided complementary information to γ through different combinations in multivariate models based on lineal discriminant analysis and report classification performance improvement of about 5%. These results suggest that the presence of larger harmonics and a proportionally smaller DF peak is associated with a decreased probability of AF recurrence after CA.

Simultaneous Whole-Chamber Non-contact Mapping of Highest Dominant Frequency Sites During Persistent Atrial Fibrillation: A Prospective Ablation Study

Purpose

Sites of highest dominant frequency (HDF) are implicated by many proposed mechanisms underlying persistent atrial fibrillation (persAF). We hypothesized that prospectively identifying and ablating dynamic left atrial HDF sites would favorably impact the electrophysiological substrate of persAF. We aim to assess the feasibility of prospectively identifying HDF sites by global simultaneous left atrial mapping.

Methods

PersAF patients with no prior ablation history underwent global simultaneous left atrial non-contact mapping. 30 s of electrograms recorded during AF were exported into a bespoke MATLAB interface to identify HDF regions, which were then targeted for ablation, prior to pulmonary vein isolation. Following ablation of each region, change in AF cycle length (AFCL) was documented (≥ 10 ms considered significant). Baseline isopotential maps of ablated regions were retrospectively analyzed looking for rotors and focal activation or extinction events.

Results

A total of 51 HDF regions were identified and ablated in 10 patients (median DF 5.8Hz, range 4.4-7.1Hz). An increase in AFCL of was seen in 20 of the 51 regions (39%), including AF termination in 4 patients. 5 out of 10 patients (including the 4 patients where AF termination occurred with HDF-guided ablation) were free from AF recurrence at 1 year. The proportion of HDF occurrences in an ablated region was not associated with change in AFCL (τ = 0.11, p = 0.24). Regions where AFCL decreased by 10 ms or more (i.e., AF disorganization) after ablation also showed lowest baseline spectral organization (p < 0.033 for any comparison). Considering all ablated regions, the average proportion of HDF events which were also HRI events was 8.0 ± 13%. Focal activations predominated (537/1253 events) in the ablated regions on isopotential maps, were modestly associated with the proportion of HDF occurrences represented by the ablated region (Kendall's τ = 0.40, p < 0.0001), and very strongly associated with focal extinction events (τ = 0.79, p < 0.0001). Rotors were rare (4/1253 events).

Conclusion

Targeting dynamic HDF sites is feasible and can be efficacious, but lacks specificity in identifying relevant human persAF substrate. Spectral organization may have an adjunctive role in preventing unnecessary substrate ablation. Dynamic HDF sites are not associated with observable rotational activity on isopotential mapping, but epi-endocardial breakthroughs could be contributory.

Patterned Illumination Techniques in Optogenetics: An Insight Into Decelerating Murine Hearts

Much has been reported about optogenetic based cardiac arrhythmia treatment and the corresponding characterization of photostimulation parameters, but still, our capacity to interact with the underlying spatiotemporal excitation patterns relies mainly on electrical and/or pharmacological approaches. However, these well-established treatments have always been an object of somehow heated discussions. Though being acutely life-saving, they often come with potential side-effects leading to a decreased functionality of the complex cardiac system. Recent optogenetic studies showed the feasibility of the usage of photostimulation as a defibrillation method with comparatively high success rates. Although, these studies mainly concentrated on the description as well as on the comparison of single photodefibrillation approaches, such as locally focused light application and global illumination, less effort was spent on the description of excitation patterns during actual photostimulation. In this study, the authors implemented a multi-site photodefibrillation technique in combination with Multi-Lead electrocardiograms (ECGs). The technical connection of real-time heart rhythm measurements and the arrhythmia counteracting light control provides a further step toward automated arrhythmia classification, which can lead to adaptive photodefibrillation methods. In order to show the power effectiveness of the new approach, transgenic murine hearts expressing channelrhodopsin-2 ex vivo were investigated using circumferential micro-LED and ECG arrays. Thus, combining the best of two methods by giving the possibility to illuminate either locally or globally with differing pulse parameters. The optical technique presented here addresses a number of challenges of technical cardiac optogenetics and is discussed in the context of arrhythmic development during photostimulation.

Critical appraisal of technologies to assess electrical activity during atrial fibrillation: a position paper from the European Heart Rhythm Association and European Society of Cardiology Working Group on eCardiology in collaboration with the Heart Rhythm Society, Asia Pacific Heart Rhythm Society, Latin American Heart Rhythm Society and Computing in Cardiology

We aim to provide a critical appraisal of basic concepts underlying signal recording and processing technologies applied for (i) atrial fibrillation (AF) mapping to unravel AF mechanisms and/or identifying target sites for AF therapy and (ii) AF detection, to optimize usage of technologies, stimulate research aimed at closing knowledge gaps, and developing ideal AF recording and processing technologies. Recording and processing techniques for assessment of electrical activity during AF essential for diagnosis and guiding ablative therapy including body surface electrocardiograms (ECG) and endo- or epicardial electrograms (EGM) are evaluated. Discussion of (i) differences in uni-, bi-, and multi-polar (omnipolar/Laplacian) recording modes, (ii) impact of recording technologies on EGM morphology, (iii) global or local mapping using various types of EGM involving signal processing techniques including isochronal-, voltage- fractionation-, dipole density-, and rotor mapping, enabling derivation of parameters like atrial rate, entropy, conduction velocity/direction, (iv) value of epicardial and optical mapping, (v) AF detection by cardiac implantable electronic devices containing various detection algorithms applicable to stored EGMs, (vi) contribution of machine learning (ML) to further improvement of signals processing technologies. Recording and processing of EGM (or ECG) are the cornerstones of (body surface) mapping of AF. Currently available AF recording and processing technologies are mainly restricted to specific applications or have technological limitations. Improvements in AF mapping by obtaining highest fidelity source signals (e.g. catheter-electrode combinations) for signal processing (e.g. filtering, digitization, and noise elimination) is of utmost importance. Novel acquisition instruments (multi-polar catheters combined with improved physical modelling and ML techniques) will enable enhanced and automated interpretation of EGM recordings in the near future.

The Effects of Fibrotic Cell Type and Its Density on Atrial Fibrillation Dynamics: An In Silico Study

Remodeling in atrial fibrillation (AF) underlines the electrical and structural changes in the atria, where fibrosis is a hallmark of arrhythmogenic structural alterations. Fibrosis is an important feature of the AF substrate and can lead to abnormal conduction and, consequently, mechanical dysfunction. The fibrotic process comprises the presence of fibrotic cells, including fibroblasts, myofibroblasts and fibrocytes, which play an important role during fibrillatory dynamics. This work assesses the effect of the diffuse fibrosis density and the intermingled presence of the three types of fibrotic cells on the dynamics of persistent AF. For this purpose, the three fibrotic cells were electrically coupled to cardiomyocytes in a 3D realistic model of human atria. Low (6.25%) and high (25%) fibrosis densities were implemented in the left atrium according to a diffuse fibrosis representation. We analyze the action potential duration, conduction velocity and fibrillatory conduction patterns. Additionally, frequency analysis was performed in 50 virtual electrograms. The tested fibrosis configurations generated a significant conduction velocity reduction, where the larger effect was observed at high fibrosis density (up to 82% reduction in the fibrocytes configuration). Increasing the fibrosis density intensifies the vulnerability to multiple re-entries, zigzag propagation, and chaotic activity in the fibrillatory conduction. The most complex propagation patterns were observed at high fibrosis densities and the fibrocytes are the cells with the largest proarrhythmic effect. Left-to-right dominant frequency gradients can be observed for all fibrosis configurations, where the fibrocytes configuration at high density generates the most significant gradients (up to 4.5 Hz). These results suggest the important role of different fibrotic cell types and their density in diffuse fibrosis on the chaotic propagation patterns during persistent AF.

Automatic Extraction of Recurrent Patterns of High Dominant Frequency Mapping During Human Persistent Atrial Fibrillation

Purpose: Identifying targets for catheter ablation remains challenging in persistent atrial fibrillation (persAF). The dominant frequency (DF) of atrial electrograms during atrial fibrillation (AF) is believed to primarily reflect local activation. Highest DF (HDF) might be responsible for the initiation and perpetuation of persAF. However, the spatiotemporal behavior of DF remains not fully understood. Some DFs during persAF were shown to lack spatiotemporal stability, while others exhibit recurrent behavior. We sought to develop a tool to automatically detect recurrent DF patterns in persAF patients. Methods: Non-contact mapping of the left atrium (LA) was performed in 10 patients undergoing persAF HDF ablation. 2,048 virtual electrograms (vEGMs, EnSite Array, Abbott Laboratories, USA) were collected for up to 5 min before and after ablation. Frequency spectrum was estimated using fast Fourier transform and DF was identified as the peak between 4 and 10 Hz and organization index (OI) was calculated. The HDF maps were identified per 4-s window and an automated pattern recognition algorithm was used to find recurring HDF spatial patterns. Dominant patterns (DPs) were defined as the HDF pattern with the highest recurrence. Results: DPs were found in all patients. Patients in atrial flutter after ablation had a single DP over the recorded time period. The time interval (median [IQR]) of DP recurrence for the patients in AF after ablation (7 patients) decreased from 21.1 s [11.8 49.7 s] to 15.7 s [6.5 18.2 s]. The DF inside the DPs presented lower temporal standard deviation (0.18 ± 0.06 Hz vs. 0.29 ± 0.12 Hz, p < 0.05) and higher OI (0.35 ± 0.03 vs. 0.31 ± 0.04, p < 0.05). The atrial regions with the highest proportion of HDF region were the septum and the left upper pulmonary vein. Conclusion: Multiple recurrent spatiotemporal HDF patterns exist during persAF. The proposed method can identify and quantify the spatiotemporal repetition of the HDFs, where the high recurrences of DP may suggest a more organized rhythm. DPs presented a more consistent DF and higher organization compared with non-DPs, suggesting that DF with higher OI might be more likely to recur. Recurring patterns offer a more comprehensive dynamic insight of persAF behavior, and ablation targeting such regions may be beneficial.

Pulmonary Delivery of Antiarrhythmic Drugs for Rapid Conversion of New-Onset Atrial Fibrillation

Pharmacologic management of atrial fibrillation (AF) is a pressing problem. This arrhythmia afflicts >5 million individuals in the United States and prevalence is estimated to rise to 12 million by 2050. Although the pill-in-the-pocket regimen for self-administered AF cardioversion introduced over a decade ago has proven useful, significant drawbacks exist. Among these are the relatively long latency of effects in the range of hours along with potential for hypotension and other adverse effects. This experience prompted development of a new strategy for increasing plasma concentrations of antiarrhythmic drugs rapidly and for a limited time, namely, pulmonary delivery. In preclinical studies in Yorkshire pigs, intratracheal administration of flecainide was shown to cause a rapid, reproducible increase in plasma drug levels. Moreover, pulmonary delivery of flecainide converted AF to normal sinus rhythm by prolonging atrial depolarization, which slows intra-atrial conduction and seems to be directly correlated with efficacy in converting AF. The rapid rise in plasma flecainide levels optimizes its anti-AF effects while minimizing adverse influences on ventricular depolarization and contractility. A more concentrated and soluble formulation of flecainide using a novel cyclodextrin complex excipient reduced net drug delivery for AF conversion when compared to the acetate formulation. Inhalation of the beta-adrenergic blocking agent metoprolol slows ventricular rate and can also terminate AF. In human subjects, oral inhalation of flecainide acetate with a hand-held, breath-actuated nebulizer results in signature prolongation of the QRS complex without serious adverse events. Thus, pulmonary delivery is a promising advance in pharmacologic approach to management of AF.

Region-specific parasympathetic nerve remodeling in the left atrium contributes to creation of a vulnerable substrate for atrial fibrillation

Atrial fibrillation (AF) is the most common heart rhythm disorder and a major cause of stroke. Unfortunately, current therapies for AF are suboptimal, largely because the molecular mechanisms underlying AF are poorly understood. Since the autonomic nervous system is thought to increase vulnerability to AF, we used a rapid atrial pacing (RAP) canine model to investigate the anatomic and electrophysiological characteristics of autonomic remodeling in different regions of the left atrium. RAP led to marked hypertrophy of parent nerve bundles in the posterior left atrium (PLA), resulting in a global increase in parasympathetic and sympathetic innervation throughout the left atrium. Parasympathetic fibers were more heterogeneously distributed in the PLA when compared with other left atrial regions; this led to greater fractionation and disorganization of AF electrograms in the PLA. Computational modeling revealed that heterogeneously distributed parasympathetic activity exacerbates sympathetic substrate for wave break and reentry. We further discovered that levels of nerve growth factor (NGF) were greatest in the left atrial appendage (LAA), where AF was most organized. Preferential NGF release by the LAA - likely a direct function of frequency and regularity of atrial stimulation - may have important implications for creation of a vulnerable AF substrate.

Rotors Detected by Phase Analysis of Filtered, Epicardial Atrial Fibrillation Electrograms Colocalize With Regions of Conduction Block

BACKGROUND

Several recent studies suggest rotors detected by phase mapping may act as main drivers of persistent atrial fibrillation. However, the electrophysiological nature of detected rotors remains unclear. We performed a direct, 1:1 comparison between phase and activation time mapping in high-density, epicardial, direct-contact mapping files of human atrial fibrillation.

METHODS

Thirty-eight unipolar electrogram files of 10 s duration were recorded in patients with atrial fibrillation (n=20 patients) using a 16×16 electrode array placed on the epicardial surface of the left atrial posterior wall or the right atrial free wall. Phase maps and isochrone wave maps were constructed for all recordings. For each detected phase singularity (PS) with a lifespan of >1 cycle length, the corresponding conduction pattern was investigated in the isochrone wave maps.

RESULTS

When using sinusoidal recomposition and Hilbert Transform, 138 PSs were detected. One hundred and four out of 138 PSs were detected within 1 electrode distance (1.5 mm) from a line of conduction block between nonrotating wavefronts detected by activation mapping. Far fewer rotating wavefronts were detected when rotating activity was identified based on wave mapping (18 out of 8219 detected waves). Fourteen out of these 18 cases were detected as PSs in phase mapping. Phase analysis of filtered electrograms produced by simulated wavefronts separated by conduction block also identified PSs on the line of conduction block.

CONCLUSIONS

PSs identified by phase analysis of filtered epicardial electrograms colocalize with conduction block lines identified by activation mapping. Detection of PSs using phase analysis has a low specificity for identifying rotating wavefronts during human atrial fibrillation using activation mapping.

The stationary wavelet transform as an efficient reductor of powerline interference for atrial bipolar electrograms in cardiac electrophysiology

OBJECTIVE

The most relevant source of signal contamination in the cardiac electrophysiology (EP) laboratory is the ubiquitous powerline interference (PLI). To reduce this perturbation, algorithms including common fixed-bandwidth and adaptive-notch filters have been proposed. Although such methods have proven to add artificial fractionation to intra-atrial electrograms (EGMs), they are still frequently used. However, such morphological alteration can conceal the accurate interpretation of EGMs, specially to evaluate the mechanisms supporting atrial fibrillation (AF), which is the most common cardiac arrhythmia. Given the clinical relevance of AF, a novel algorithm aimed at reducing PLI on highly contaminated bipolar EGMs and, simultaneously, preserving their morphology is proposed.

APPROACH

The method is based on the wavelet shrinkage and has been validated through customized indices on a set of synthesized EGMs to accurately quantify the achieved level of PLI reduction and signal morphology alteration. Visual validation of the algorithm's performance has also been included for some real EGM excerpts.

MAIN RESULTS

The method has outperformed common filtering-based and wavelet-based strategies in the analyzed scenario. Moreover, it possesses advantages such as insensitivity to amplitude and frequency variations in the PLI, and the capability of joint removal of several interferences.

SIGNIFICANCE

The use of this algorithm in routine cardiac EP studies may enable improved and truthful evaluation of AF mechanisms.

Systematic differences of non-invasive dominant frequency estimation compared to invasive dominant frequency estimation in atrial fibrillation

Non-invasive analysis of atrial fibrillation (AF) using body surface mapping (BSM) has gained significant interest, with attempts at interpreting atrial spectro-temporal parameters from body surface signals. As these body surface signals could be affected by properties of the torso volume conductor, this interpretation is not always straightforward. This paper highlights the volume conductor effects and influences of the algorithm parameters for identifying the dominant frequency (DF) from cardiac signals collected simultaneously on the torso and atrial surface. Bi-atrial virtual electrograms (VEGMs) and BSMs were recorded simultaneously for 5 min from 10 patients undergoing ablation for persistent AF. Frequency analysis was performed on 4 s segments. DF was defined as the frequency with highest power between 4 and 10 Hz with and without applying organization index (OI) thresholds. The volume conductor effect was assessed by analyzing the highest DF (HDF) difference of each VEGM HDF against its BSM counterpart. Significant differences in HDF values between intra-cardiac and torso signals could be observed, independent of OI threshold. This difference increases with increasing endocardial HDF (BSM-VEGM median difference from -0.13 Hz for VEGM HDF at 6.25 ± 0.25 Hz to -4.24 Hz at 9.75 ± 0.25 Hz), thereby confirming the theory of the volume conductor effect in real-life situations. Applying an OI threshold strongly affected the BSM HDF area size and location and atrial HDF area location. These results suggest that volume conductor and measurement algorithm effects must be considered for appropriate clinical interpretation.

Validation and Opportunities of Electrocardiographic Imaging: From Technical Achievements to Clinical Applications

Electrocardiographic imaging (ECGI) reconstructs the electrical activity of the heart from a dense array of body-surface electrocardiograms and a patient-specific heart-torso geometry. Depending on how it is formulated, ECGI allows the reconstruction of the activation and recovery sequence of the heart, the origin of premature beats or tachycardia, the anchors/hotspots of re-entrant arrhythmias and other electrophysiological quantities of interest. Importantly, these quantities are directly and non-invasively reconstructed in a digitized model of the patient's three-dimensional heart, which has led to clinical interest in ECGI's ability to personalize diagnosis and guide therapy. Despite considerable development over the last decades, validation of ECGI is challenging. Firstly, results depend considerably on implementation choices, which are necessary to deal with ECGI's ill-posed character. Secondly, it is challenging to obtain (invasive) ground truth data of high quality. In this review, we discuss the current status of ECGI validation as well as the major challenges remaining for complete adoption of ECGI in clinical practice. Specifically, showing clinical benefit is essential for the adoption of ECGI. Such benefit may lie in patient outcome improvement, workflow improvement, or cost reduction. Future studies should focus on these aspects to achieve broad adoption of ECGI, but only after the technical challenges have been solved for that specific application/pathology. We propose 'best' practices for technical validation and highlight collaborative efforts recently organized in this field. Continued interaction between engineers, basic scientists, and physicians remains essential to find a hybrid between technical achievements, pathological mechanisms insights, and clinical benefit, to evolve this powerful technique toward a useful role in clinical practice.

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.

Atrial Fibrillation Dynamics and Ionic Block Effects in Six Heterogeneous Human 3D Virtual Atria with Distinct Repolarization Dynamics

Atrial fibrillation (AF) usually manifests as reentrant circuits propagating through the whole atria creating chaotic activation patterns. Little is yet known about how differences in electrophysiological and ionic properties between patients modulate reentrant patterns in AF. The goal of this study is to quantify how variability in action potential duration (APD) at different stages of repolarization determines AF dynamics and their modulation by ionic block using a set of virtual whole-atria human models. Six human whole-atria models are constructed based on the same anatomical structure and fiber orientation, but with different electrophysiological phenotypes. Membrane kinetics for each whole-atria model are selected with distinct APD characteristics at 20, 50, and 90% repolarization, from an experimentally calibrated population of human atrial action potential models, including AF remodeling and acetylcholine parasympathetic effects. Our simulations show that in all whole-atria models, reentrant circuits tend to organize around the pulmonary veins and the right atrial appendage, thus leading to higher dominant frequency (DF) and more organized activation in the left atrium than in the right atrium. Differences in APD in all phases of repolarization (not only APD₉₀) yielded quantitative differences in fibrillation patterns with long APD associated with slower and more regular dynamics. Long APD₅₀ and APD₂₀ were associated with increased interatrial conduction block and interatrial differences in DF and organization index, creating reentry instability and self-termination in some cases. Specific inhibitions of I_K1, I_NaK, or I_Na reduce DF and organization of the arrhythmia by enlarging wave meandering, reducing the number of secondary wavelets, and promoting interatrial block in all six virtual patients, especially for the phenotypes with short APD at 20, 50, and/or 90% repolarization. This suggests that therapies aiming at prolonging the early phase of repolarization might constitute effective antiarrhythmic strategies for the pharmacological management of AF. In summary, simulations report significant differences in atrial fibrillatory dynamics resulting from differences in APD at all phases of repolarization.

Waveform Integrity in Atrial Fibrillation: The Forgotten Issue of Cardiac Electrophysiology

Atrial fibrillation (AF) is the most common arrhythmia in clinical practice with an increasing prevalence of about 15% in the elderly. Despite other alternatives, catheter ablation is currently considered as the first-line therapy for the treatment of AF. This strategy relies on cardiac electrophysiology systems, which use intracardiac electrograms (EGM) as the basis to determine the cardiac structures contributing to sustain the arrhythmia. However, the noise-free acquisition of these recordings is impossible and they are often contaminated by different perturbations. Although suppression of nuisance signals without affecting the original EGM pattern is essential for any other later analysis, not much attention has been paid to this issue, being frequently considered as trivial. The present work introduces the first thorough study on the significant fallout that regular filtering, aimed at reducing acquisition noise, provokes on EGM pattern morphology. This approach has been compared with more refined denoising strategies. Performance has been assessed both in time and frequency by well established parameters for EGM characterization. The study comprised synthesized and real EGMs with unipolar and bipolar recordings. Results reported that regular filtering altered substantially atrial waveform morphology and was unable to remove moderate amounts of noise, thus turning time and spectral characterization of the EGM notably inaccurate. Methods based on Wavelet transform provided the highest ability to preserve EGM morphology with improvements between 20 and beyond 40%, to minimize dominant atrial frequency estimation error with up to 25% reduction, as well as to reduce huge levels of noise with up to 10 dB better reduction. Consequently, these algorithms are recommended as a replacement of regular filtering to avoid significant alterations in the EGMs. This could lead to more accurate and truthful analyses of atrial activity dynamics aimed at understanding and locating the sources of AF.

An interactive platform to guide catheter ablation in human persistent atrial fibrillation using dominant frequency, organization and phase mapping

BACKGROUND AND OBJECTIVE

Optimal targets for persistent atrial fibrillation (persAF) ablation are still debated. Atrial regions hosting high dominant frequency (HDF) are believed to participate in the initiation and maintenance of persAF and hence are potential targets for ablation, while rotor ablation has shown promising initial results. Currently, no commercially available system offers the capability to automatically identify both these phenomena. This paper describes an integrated 3D software platform combining the mapping of both frequency spectrum and phase from atrial electrograms (AEGs) to help guide persAF ablation in clinical cardiac electrophysiological studies.

METHODS

30s of 2048 non-contact AEGs (EnSite Array, St. Jude Medical) were collected and analyzed per patient. After QRST removal, the AEGs were divided into 4s windows with a 50% overlap. Fast Fourier transform was used for DF identification. HDF areas were identified as the maximum DF to 0.25Hz below that, and their centers of gravity (CGs) were used to track their spatiotemporal movement. Spectral organization measurements were estimated. Hilbert transform was used to calculate instantaneous phase.

RESULTS

The system was successfully used to guide catheter ablation for 10 persAF patients. The mean processing time was 10.4 ± 1.5min, which is adequate comparing to the normal electrophysiological (EP) procedure time (120∼180min).

CONCLUSIONS

A customized software platform capable of measuring different forms of spatiotemporal AEG analysis was implemented and used in clinical environment to guide persAF ablation. The modular nature of the platform will help electrophysiological studies in understanding of the underlying AF mechanisms.

Electrocardiographic Spectral Features for Long-Term Outcome Prognosis of Atrial Fibrillation Catheter Ablation

Atrial fibrillation (AF) is the most common arrhythmia in routine clinical practice. Despite many years of research, its mechanisms still are not well understood, thus reducing the effectiveness of AF treatments. Nowadays, pulmonary vein isolation by catheter ablation is the treatment of choice for AF resistant either to pharmacological or electrical cardioversion. However, given that long-term recurrences are common, an appropriate patient selection before the procedure is of paramount relevance in the improvement of AF catheter ablation outcome. The present work studies how several spectral features of the atrial activity (AA) from a single lead of the surface electrocardiogram (ECG) can become potential pre-ablation predictors of long-term (>2 months) sinus rhythm maintenance. Among all the analyzed spectral features, results indicated that the most significant single predictor of paroxysmal AF ablation treatment outcome was related to the amplitude of the first harmonic of the dominant frequency, providing sensitivity (Se), specificity (Sp) and accuracy (Ac) values of 90%, 42.86 and 77.78%, respectively. On the other hand, the AA harmonic structure was the most significant single predictor for persistent AF, with Se, Sp and Ac values of 100%, 54.55 and 77.27%, respectively. A logistic regression analysis, mainly based on spectral amplitudes as well as on the harmonic structure of the AA, provided a higher predictive ability both for paroxysmal AF (Se = 100%, Sp = 57.14% and Ac = 88.89%) and persistent AF (Se = 90.91%, Sp = 72.73 and Ac = 81.82%). In conclusion, the study of key AA spectral features from the surface ECG can provide a significant preoperative prognosis of AF catheter ablation outcome at long-term follow-up.

Several insights into the preprocessing of electrograms in atrial fibrillation for dominant frequency analysis

Background

Dominant frequency (DF) analysis of atrial electrograms has become an important method in characterizing atrial fibrillation (AF). As a classic method, Botteron’s approach is widely used in the preprocessing of frequency analysis during AF. It includes three steps: (1) band-pass filtering at 40–250 Hz, (2) absolute value, and (3) low-pass filtering at 20 Hz. This paper aims to expound the necessity and adjustability of each step.

Methods and results

Unipolar epicardial mapping signals were recorded during AF from eight mongrel dogs with cholinergic AF model. Episodes of these data were randomly selected to evaluate the impact of different pass bands and the necessity of low-pass filtering with 20 Hz cutoff frequency. Each episode of AF signal is 5 s long with a sampling rate of 2 kHz. Simulated electrograms were adopted to discuss the role of taking absolute value. Furthermore, direct spectral analysis method (FFT et al.) is compared with Botteron’s preprocessing approach. According to our statistical analysis, the pass band of 40–250 Hz was not the best, while 20–100 Hz presented the high accuracy rate of DF. From the comparing result of direct FFT without Botteron’s approach we deduced that the rectification of absolute value was meaningful for the fundamental atrial frequency. The final step, 20 Hz low-pass filter can completely be omitted in DF analysis. In consideration of the demand for real-time distribution of DF in clinical or experimental situations, down-sampling method and the impact of ventricular artifacts on DF was also discussed.

Conclusion

In the actual application of the three preprocessing steps, the pass band selection of band-pass filter can be adjusted and the rectification of taking absolute value is important. Nevertheless, the final step of 20 Hz low-pass filter is totally unnecessary. In real-time signal processing situations, taking down-sampling method and ignoring the ventricular artifacts can also have high performance in DF analysis of atrial electrograms.

Developing a New Computer-Aided Clinical Decision Support System for Prediction of Successful Postcardioversion Patients with Persistent Atrial Fibrillation

We propose a new algorithm to predict the outcome of direct-current electric (DCE) cardioversion for atrial fibrillation (AF) patients. AF is the most common cardiac arrhythmia and DCE cardioversion is a noninvasive treatment to end AF and return the patient to sinus rhythm (SR). Unfortunately, there is a high risk of AF recurrence in persistent AF patients; hence clinically it is important to predict the DCE outcome in order to avoid the procedure's side effects. This study develops a feature extraction and classification framework to predict AF recurrence patients from the underlying structure of atrial activity (AA). A multiresolution signal decomposition technique, based on matching pursuit (MP), was used to project the AA over a dictionary of wavelets. Seven novel features were derived from the decompositions and were employed in a quadratic discrimination analysis classification to predict the success of post-DCE cardioversion in 40 patients with persistent AF. The proposed algorithm achieved 100% sensitivity and 95% specificity, indicating that the proposed computational approach captures detailed structural information about the underlying AA and could provide reliable information for effective management of AF.

The effects of remodeling with heart failure on mode of initiation of ventricular fibrillation and its spatiotemporal organization

PURPOSE

The effect of the heart failure substrate on the initiation of ventricular fibrillation (VF) and its resulting mechanism is not known. The objective of this study was to determine the effects of substrate on VF initiation and its spatiotemporal organization in the heart failure model.

METHODS

Optical action potentials were recorded from LV wedge preparations either from structurally normal hearts (control, n = 11) or from congestive heart failure (CHF; n = 7), at the epicardial surface, endocardial surface which included a papillary muscle, and a transmural cross section. Action potential duration (APD(80)) was determined, and VF was initiated. A fast Fourier transform was calculated, and the dominant frequency (DF) was determined.

RESULTS

The CHF group showed increased VF vulnerability (69 vs 26 %, p < 0.03), and every mapped surface showed an APD(80) gradient which included islands of higher APDs on the transmural surface (M cells) which was not observed in controls. VF in the CHF group was characterized by stable, discrete, high-DF areas that correlated to either foci or spiral waves located on the transmural surface at the site of the papillary muscle. Overall, the top 10 % of DFs correlated to an APD of 101 ms while the bottom 10 % of DFs correlated to an APD of 126 ms (p < 0.01).

CONCLUSIONS

In the CHF model, APD gradients correlated with an increased vulnerability to VF, and the highest stable DFs were located on the transmural surface which was not seen in controls. This indicates that the CHF substrate creates unique APD and DF characteristics.

Catheter ablation of persistent atrial fibrillation: The importance of substrate modification

Accumulating data have shown that elimination of atrial fibrillation (AF) sources should be the goal in persistent AF ablation. Pulmonary vein isolation, linear lesions and complex fractionated atrial electrograms (CFAEs) ablation have shown limited efficacy in patients with persistent AF. A combined approach using voltage, CFAEs and dominant frequency (DF) mapping may be helpful for the identification of AF sources and subsequent focal substrate modification. The fibrillatory activity is maintained by intramural reentry centered on fibrotic patches. Voltage mapping may assist in the identification of fibrotic areas. Stable rotors display the higher DF and possibly drive AF. Furthermore, the single rotor is usually consistent with organized AF electrograms without fractionation. It is therefore quite possible that rotors are located at relatively “healthy islands” within the patchy fibrosis. This is supported by the fact that high DF sites have been negatively correlated to the amount of fibrosis. CFAEs are located in areas adjacent to high DF. In conclusion, patchy fibrotic areas displaying the maximum DF along with high organization index and the lower fractionation index are potential targets of ablation. Prospective studies are required to validate the efficacy of substrate modification in left atrial ablation outcomes.

The logical operator map identifies novel candidate markers for critical sites in patients with atrial fibrillation

The identification of suitable markers for critical patterns during atrial fibrillation (AF) may be crucial to guide an effective ablation treatment. Single parameter maps, based on dominant frequency and complex fractionated electrograms, have been proposed as a tool for electrogram-guided ablation, however the specificity of these markers is debated. Experimental studies suggest that AF critical patterns may be identified on the basis of specific rate and organization features, where rapid organized and rapid fragmented activities characterize respectively localized sources and critical substrates. In this paper we introduce the logical operator map, a novel mapping tool for a point-by-point identification and localization of AF critical sites. Based on advanced signal and image processing techniques, the approach combines in a single map electrogram-derived rate and organization features with tomographic anatomical detail. The construction of the anatomically-detailed logical operator map is based on the time-domain estimation of atrial rate and organization in terms of cycle length and wave-similarity, the logical combination of these indexes to obtain suitable markers of critical sites, and the multimodal integration of electrophysiological and anatomical information by segmentation and registration techniques. Logical operator maps were constructed in 14 patients with persistent AF, showing the capability of the combined rate and organization markers to identify with high selectivity the subset of electrograms associated with localized sources and critical substrates. The precise anatomical localization of these critical sites revealed the confinement of rapid organized sources in the left atrium with organization and rate gradients towards the surrounding tissue, and the presence of rapid fragmented electrograms in proximity of the sources. By merging in a single map the most relevant electrophysiological and anatomical features of the AF process, the logical operator map may have significant clinical impact as a direct, comprehensive tool to understand arrhythmia mechanisms in the single patient and guide more conservative, step-wise ablation.

Mathematical approaches to understanding and imaging atrial fibrillation: significance for mechanisms and management

Atrial fibrillation (AF) is the most common sustained arrhythmia in humans. The mechanisms that govern AF initiation and persistence are highly complex, of dynamic nature, and involve interactions across multiple temporal and spatial scales in the atria. This article aims to review the mathematical modeling and computer simulation approaches to understanding AF mechanisms and aiding in its management. Various atrial modeling approaches are presented, with descriptions of the methodological basis and advancements in both lower-dimensional and realistic geometry models. A review of the most significant mechanistic insights made by atrial simulations is provided. The article showcases the contributions that atrial modeling and simulation have made not only to our understanding of the pathophysiology of atrial arrhythmias, but also to the development of AF management approaches. A summary of the future developments envisioned for the field of atrial simulation and modeling is also presented. The review contends that computational models of the atria assembled with data from clinical imaging modalities that incorporate electrophysiological and structural remodeling could become a first line of screening for new AF therapies and approaches, new diagnostic developments, and new methods for arrhythmia prevention.

Preoperative study of the surface ECG for the prognosis of atrial fibrillation maze surgery outcome at discharge

The Cox-maze surgery is an effective procedure for terminating atrial fibrillation (AF) in patients requiring open-heart surgery associated with another heart disease. After the intervention, regardless of the patient's rhythm, all are treated with oral anticoagulants and antiarrhythmic drugs prior to discharge. Furthermore, patients maintaining AF before discharge could also be treated with electrical cardioversion (ECV). In view of this, a preoperative prognosis of the patient's rhythm at discharge would be helpful for optimizing drug therapy planning as well as for advancing ECV therapy. This work analyzes 30 preoperative electrocardiograms (ECGs) from patients suffering from AF in order to predict the Cox-maze surgery outcome at discharge. Two different characteristics of the AF pattern have been studied. On the one hand, the atrial activity (AA) organization, which provides information about the number of propagating wavelets in the atria, was investigated. AA organization has been successfully used in previous studies related to spontaneous reversion of paroxysmal AF and to the outcome of ECV. To assess organization, the dominant atrial frequency (DAF) and sample entropy (SampEn) have been computed. On the other hand, the second characteristic studied was the fibrillatory wave (f-wave) amplitude, which has been demonstrated to be a valuable indicator of the Cox-maze surgery outcome in previous studies. Moreover, this parameter has been obtained through a new methodology, based on computing the f-wave average power (fWP). Finally, all the computed indices were combined in a decision tree in order to improve prediction capability. Results for the DAF yielded a sensitivity (Se), a specificity (Sp) and an accuracy (Acc) of 61.54%, 82.35% and 73.33%, respectively. For SampEn the values were 69.23%, 76.00% and 73.33%, respectively, and for fWP they were 92.31%, 82.35% and 86.67%, respectively. Finally, the decision tree combining the three parameters analyzed improved the preoperative prognosis of the Cox-maze outcome with values of Se, Sp and Acc of 100%, 82.35% and 90%, respectively. As a consequence, the analysis of parameters related to the f-wave pattern, extracted from the preoperative ECG, has provided a considerable ability to predict the outcome of AF Cox-maze surgery at discharge.

Evaluating epicardial mapping electrogram by the method of dominant frequency and Lorenz plot

The methods of dominant frequency and Lorenz plot are used in this study to evaluate the activation rate and the activation rate variability of cardiac signals during atrial fibrillation. An epicardial mapping system was applied to acquire the atrial electrogram of mongrel dogs. The dominant frequency and Lorenz plot of each signal from various myocardial regions of the atria were analyzed. Our results show that both a frequency gradient and a variability gradient exist in the atria and the roots of pulmonary veins. The dominant frequencies of the anterior atria are higher than the posterior ones and the activation variability of both atria was higher than those of the pulmonary veins. A combination of these two methods may provide a more comprehensive understanding of the electrophysiology mechanism associated with atrial fibrillation.