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.

High-resolution mapping of reentrant atrial tachycardias: Relevance of low bipolar voltage

BACKGROUND

Bipolar voltage is widely used to characterize the atrial substrate but has been poorly validated, particularly during clinical tachycardias.

OBJECTIVE

The purpose of this study was to evaluate the diagnostic performance of voltage thresholds for identifying regions of slow conduction during reentrant atrial tachycardias (ATs).

METHODS

Thirty bipolar voltage and activation maps created during reentrant ATs were analyzed to (1) examine the relationship between voltage amplitude and conduction velocity (CV), (2) measure the diagnostic ability of voltage thresholds to predict CV, and (3) identify determinants of AT circuit dimensions. Voltage amplitude was categorized as "normal" (>0.50 mV), "abnormal" (0.05-0.50 mV), or "scar" (<0.05 mV); slow conduction was defined as <30 cm/s.

RESULTS

A total of 266,457 corresponding voltage and CV data points were included for analysis. Voltage and CV were moderately correlated (r = 0.407; P < .001). Bipolar voltage predicted regions of slow conduction with an area under the receiver operating characteristic curve of 0.733 (95% confidence interval 0.731-0.735). A threshold of 0.50 mV had 91% sensitivity and 35% specificity for identifying slow conduction, whereas 0.05 mV had 36% sensitivity and 87% specificity, with an optimal voltage threshold of 0.15 mV. Analyses restricted to the AT circuits identified weaker associations between voltage and CV and an optimal voltage threshold of 0.25 mV.

CONCLUSION

Widely used bipolar voltage amplitude thresholds to define "abnormal" and "scar" tissue in the atria are, respectively, sensitive and specific for identifying regions of slow conduction during reentrant ATs. However, overall, the association of voltage with CV is modest. No clinical predictors of AT circuit dimensions were identified.

Atrial conduction velocity mapping: clinical tools, algorithms and approaches for understanding the arrhythmogenic substrate

Characterizing patient-specific atrial conduction properties is important for understanding arrhythmia drivers, for predicting potential arrhythmia pathways, and for personalising treatment approaches. One metric that characterizes the health of the myocardial substrate is atrial conduction velocity, which describes the speed and direction of propagation of the electrical wavefront through the myocardium. Atrial conduction velocity mapping algorithms are under continuous development in research laboratories and in industry. In this review article, we give a broad overview of different categories of currently published methods for calculating CV, and give insight into their different advantages and disadvantages overall. We classify techniques into local, global, and inverse methods, and discuss these techniques with respect to their faithfulness to the biophysics, incorporation of uncertainty quantification, and their ability to take account of the atrial manifold.

Slow Conduction Corridors and Pivot Sites Characterize the Electrical Remodeling in Atrial Fibrillation

OBJECTIVES

This study aimed to evaluate the progression of electrophysiological phenomena in a cohort of patients with paroxysmal atrial fibrillation (PAF) and persistent atrial fibrillation (PsAF).

BACKGROUND

Electrical remodeling has been conjectured to determine atrial fibrillation (AF) progression.

METHODS

High-density electroanatomic maps during sinus rhythm of 20 patients with AF (10 PAF, 10 PsAF) were compared with 5 healthy control subjects (subjects undergoing ablation of a left-sided accessory pathway). A computational postprocessing of electroanatomic maps was performed to identify specific electrophysiological phenomena: slow conductions corridors, defined as discrete areas of conduction velocity <50 cm/s, and pivot points, defined as sites showing high wave-front curvature documented by a curl module >2.5 1/s.

RESULTS

A progressive decrease of mean conduction velocity was recorded across the groups (111.6 ± 55.5 cm/s control subjects, 97.1 ± 56.3 cm/s PAF, and 84.7 ± 55.7 cm/s PsAF). The number and density of slow conduction corridors increase in parallel with the progression of AF (8.6 ± 2.2 control subjects, 13.3 ± 3.2 PAF, and 20.5 ± 4.5 PsAF). In PsAF the atrial substrate is characterized by a higher curvature of wave-front propagation (0.86 ± 0.71 1/s PsAF vs 0.74 ± 0.63 1/s PAF; P = 0.003) and higher number of pivot points (25.1 ± 13.8 PsAF vs 9.5 ± 6.7 PAF; P < 0.0001). Slow conductions: corridors were mostly associated with pivot sites tending to cluster around pulmonary veins antra.

CONCLUSIONS

The electrical remodeling hinges mainly on corridors of slow conduction and higher curvature of wave-front propagation. Pivot points associated to SC corridors may be the major determinants for functional localized re-entrant circuits creating the substrate for maintenance of AF.

Manifold Approximating Graph Interpolation of Cardiac Local Activation Time

Objective:

Local activation time (LAT) mapping of cardiac chambers is vital for targeted treatment of cardiac arrhythmias in catheter ablation procedures. Current methods require too many LAT observations for an accurate interpolation of the necessarily sparse LAT signal extracted from intracardiac electrograms (EGMs). Additionally, conventional performance metrics for LAT interpolation algorithms do not accurately measure the quality of interpolated maps. We propose, first, a novel method for spatial interpolation of the LAT signal which requires relatively few observations; second, a realistic sub-sampling protocol for LAT interpolation testing; and third, a new color-based metric for evaluation of interpolation quality that quantifies perceived differences in LAT maps.

Methods:

We utilize a graph signal processing framework to reformulate the irregular spatial interpolation problem into a semi-supervised learning problem on the manifold with a closed-form solution. The metric proposed uses a color difference equation and color theory to quantify visual differences in generated LAT maps.

Results:

We evaluate our approach on a dataset consisting of seven LAT maps from four patients obtained by the CARTO electroanatomic mapping system during premature ventricular complex (PVC) ablation procedures. Random sub-sampling and re-interpolation of the LAT observations show excellent accuracy for relatively few observations, achieving on average 6% lower error than state-of-the-art techniques for only 100 observations.

Conclusion:

Our study suggests that graph signal processing methods can improve LAT mapping for cardiac ablation procedures.

Significance:

The proposed method can reduce patient time in surgery by decreasing the number of LAT observations needed for an accurate LAT map.

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.

A Computational Study of the Electrophysiological Substrate in Patients Suffering From Atrial Fibrillation

In the context of cardiac electrophysiology, we propose a novel computational approach to highlight and explain the long-debated mechanisms behind atrial fibrillation (AF) and to reliably numerically predict its induction and sustainment. A key role is played, in this respect, by a new way of setting a parametrization of electrophysiological mathematical models based on conduction velocities; these latter are estimated from high-density mapping data, which provide a detailed characterization of patients' electrophysiological substrate during sinus rhythm. We integrate numerically approximated conduction velocities into a mathematical model consisting of a coupled system of partial and ordinary differential equations, formed by the monodomain equation and the Courtemanche-Ramirez-Nattel model. Our new model parametrization is then adopted to predict the formation and self-sustainment of localized reentries characterizing atrial fibrillation, by numerically simulating the onset of ectopic beats from the pulmonary veins. We investigate the paroxysmal and the persistent form of AF starting from electro-anatomical maps of two patients. The model's response to stimulation shows how substrate characteristics play a key role in inducing and sustaining these arrhythmias. Localized reentries are less frequent and less stable in case of paroxysmal AF, while they tend to anchor themselves in areas affected by severe slow conduction in case of persistent AF.

Comparison of clockwise and counterclockwise right atrial flutter using high-resolution mapping and automated velocity measurements

BACKGROUND

Only few studies have been performed that explore the electrophysiological differences between clockwise (CW) and counterclockwise (CCW) right atrial (RA) cavotricuspid isthmus (CTI)-dependent atrial flutter (AFL) using the high-resolution Rhythmia mapping system.

OBJECTIVES

We sought to compare CW and CCW CTI-dependent AFL in pure right AFL patients (pts) using the ultra-high-definition (ultra-HD) Rhythmia mapping system and we mathematically developed a cartography model based on automatic velocity RA measurements to identify electrophysiological AFL specificities.

METHODS AND RESULTS

Thirty-three pts were recruited. The mean age was 71 ± 13 years old. The sinus venosus (SV) block line was present in 32/33 of cases (97%) and no significant difference was found between CCW and CW CTI AFL (100% vs. 91%; p = .7). No line was localized in the region of the crista terminalis (CT). A superior gap was present in the posterior line in 14/31 (45.2%) but this was similarly present in CCW AFL, when compared to CW AFL (10/22 [45.5%] vs. 4/10 [40%]; p = .9). When present, the extension of the posterior line of block was observed in 18/31 pts (58%) without significant differences between CCW and CW CI AFL (12/22 [54.5%] vs. 6/10 [60%]; p = .9) The Eustachian ridge line of block was similarly present in both groups (82% [18/22] vs. 45.5% [5/11]; p = .2). The absence of the Eustachian ridge line of block led to significantly slowed velocity in this area (28 ± 10 cm/s; n = 8), and the velocities were similarly altered between both groups (26 ± 10 [4/22] vs. 29.8 ± 11 cm/s [4/11]; p = .6). We created mathematical, three-dimensional RA reconstruction-velocity model measurements. In each block localization, when the block line was absent, velocity was significantly slowed (≤20 cm/s). A systematic slowdown in conduction velocity was observed at the entrance and exit of the CTI in 100% of cases. This alteration to the conduction entrance was localized at the lateral side of the CTI for the CCW AFL and at the septal side of the CTI for CW AFL. The exit-conduction alteration was localized at the CTI septal side for the CCW AFL and at the CTI lateral side for the CW AFL.

CONCLUSION

The ultra-HD Rhythmia mapping system confirmed the absence of significant electrophysiological differences between CCW and CW AFL. The mechanistic posterior SV and Eustachian ridge block lines were confirmed in each arrhythmia. A systematic slowing down at the entrance and exit of the CTI was demonstrated in both CCW and CW AFL, but in reverse positions.

Electrocardiographic predictors of successful resynchronization of left bundle branch block by His bundle pacing

BACKGROUND

His bundle pacing (HBP) is an alternative to biventricular pacing (BVP) for delivering cardiac resynchronization therapy (CRT) in patients with heart failure and left bundle branch block (LBBB). It is not known whether ventricular activation times and patterns achieved by HBP are equivalent to intact conduction systems and not all patients with LBBB are resynchronized by HBP.

OBJECTIVE

To compare activation times and patterns of His-CRT with BVP-CRT, LBBB and intact conduction systems.

METHODS

In patients with LBBB, noninvasive epicardial mapping (ECG imaging) was performed during BVP and temporary HBP. Intrinsic activation was mapped in all subjects. Left ventricular activation times (LVAT) were measured and epicardial propagation mapping (EPM) was performed, to visualize epicardial wavefronts. Normal activation pattern and a normal LVAT range were determined from normal subjects.

RESULTS

Forty-five patients were included, 24 with LBBB and LV impairment, and 21 with normal 12-lead ECG and LV function. In 87.5% of patients with LBBB, His-CRT successfully shortened LVAT by ≥10 ms. In 33.3%, His-CRT resulted in complete ventricular resynchronization, with activation times and patterns indistinguishable from normal subjects. EPM identified propagation discontinuity artifacts in 83% of patients with LBBB. This was the best predictor of whether successful resynchronization was achieved by HBP (logarithmic odds ratio, 2.19; 95% confidence interval, 0.07-4.31; p = .04).

CONCLUSION

Noninvasive electrocardiographic mapping appears to identify patients whose LBBB can be resynchronized by HBP. In contrast to BVP, His-CRT may deliver the maximum potential ventricular resynchronization, returning activation times, and patterns to those seen in normal hearts.

Granger Causality-Based Analysis for Classification of Fibrillation Mechanisms and Localization of Rotational Drivers

BACKGROUND

The mechanisms sustaining myocardial fibrillation remain disputed, partly due to a lack of mapping tools that can accurately identify the mechanism with low spatial resolution clinical recordings. Granger causality (GC) analysis, an econometric tool for quantifying causal relationships between complex time-series, was developed as a novel fibrillation mapping tool and adapted to low spatial resolution sequentially acquired data.

METHODS

Ventricular fibrillation (VF) optical mapping was performed in Langendorff-perfused Sprague-Dawley rat hearts (n=18), where novel algorithms were developed using GC-based analysis to (1) quantify causal dependence of neighboring signals and plot GC vectors, (2) quantify global organization with the causality pairing index, a measure of neighboring causal signal pairs, and (3) localize rotational drivers (RDs) by quantifying the circular interdependence of neighboring signals with the circular interdependence value. GC-based mapping tools were optimized for low spatial resolution from downsampled optical mapping data, validated against high-resolution phase analysis and further tested in previous VF optical mapping recordings of coronary perfused donor heart left ventricular wedge preparations (n=12), and adapted for sequentially acquired intracardiac electrograms during human persistent atrial fibrillation mapping (n=16).

RESULTS

Global VF organization quantified by causality pairing index showed a negative correlation at progressively lower resolutions (50% resolution: P=0.006, R2=0.38, 12.5% resolution, P=0.004, R2=0.41) with a phase analysis derived measure of disorganization, locations occupied by phase singularities. In organized VF with high causality pairing index values, GC vector mapping characterized dominant propagating patterns and localized stable RDs, with the circular interdependence value showing a significant difference in driver versus nondriver regions (0.91±0.05 versus 0.35±0.06, P=0.0002). These findings were further confirmed in human VF. In persistent atrial fibrillation, a positive correlation was found between the causality pairing index and presence of stable RDs (P=0.0005,R2=0.56). Fifty percent of patients had RDs, with a low incidence of 0.9±0.3 RDs per patient.

CONCLUSIONS

GC-based fibrillation analysis can measure global fibrillation organization, characterize dominant propagating patterns, and map RDs using low spatial resolution sequentially acquired data.

Challenges Associated with Interpreting Mechanisms of AF

Determining optimal treatment strategies for complex arrhythmogenesis in AF is confounded by the lack of consensus regarding the mechanisms causing AF. Studies report different mechanisms for AF, ranging from hierarchical drivers to anarchical multiple activation wavelets. Differences in the assessment of AF mechanisms are likely due to AF being recorded across diverse models using different investigational tools, spatial scales and clinical populations. The authors review different AF mechanisms, including anatomical and functional re-entry, hierarchical drivers and anarchical multiple wavelets. They then describe different cardiac mapping techniques and analysis tools, including activation mapping, phase mapping and fibrosis identification. They explain and review different data challenges, including differences between recording devices in spatial and temporal resolutions, spatial coverage and recording surface, and report clinical outcomes using different data modalities. They suggest future research directions for investigating the mechanisms underlying human AF.

A technique for measuring anisotropy in atrial conduction to estimate conduction velocity and atrial fibre direction

BACKGROUND

Cardiac conduction properties exhibit large variability, and affect patient-specific arrhythmia mechanisms. However, it is challenging to clinically measure conduction velocity (CV), anisotropy and fibre direction. Our aim is to develop a technique to estimate conduction anisotropy and fibre direction from clinically available electrical recordings.

METHODS

We developed and validated automated algorithms for estimating cardiac CV anisotropy, from any distribution of recording locations on the atrial surface. The first algorithm is for elliptical wavefront fitting to a single activation map (method 1), which works well close to the pacing location, but decreases in accuracy further from the pacing location (due to spatial heterogeneity in the conductivity and fibre fields). As such, we developed a second methodology for measuring local conduction anisotropy, using data from two or three activation maps (method 2: ellipse fitting to wavefront propagation velocity vectors from multiple activation maps).

RESULTS

Ellipse fitting to CV vectors from two activation maps (method 2) leads to an improved estimation of longitudinal and transverse CV compared to method 1, but fibre direction estimation is still relatively poor. Using three activation maps with method 2 provides accurate estimation, with approximately 70% of atrial fibres estimated within 20∘. We applied the technique to clinical activation maps to demonstrate the presence of heterogeneous conduction anisotropy, and then tested the effects of this conduction anisotropy on predicted arrhythmia dynamics using computational simulation.

CONCLUSIONS

We have developed novel algorithms for calculating CV and measuring the direction dependency of atrial activation to estimate atrial fibre direction, without the need for specialised pacing protocols, using clinically available electrical recordings.

Characterizing localized reentry with high-resolution mapping: Evidence for multiple slow conducting isthmuses within the circuit

BACKGROUND

Reentrant circuits are considered to be critically dependent on a single protected slow conducting isthmus.

OBJECTIVE

The purpose of this study was to investigate conduction properties and electrogram (EGM) characteristics of the entire circuit in localized atrial reentrant circuits using high-resolution mapping.

METHODS

Fifteen localized reentrant atrial tachycardias were studied with high-resolution mapping (Rhythmia). EGMs along the entire circuit were analyzed offline for fractionation, duration, and amplitude. Maps were exported to MATLAB (MathWorks) to measure bipolar voltage and conduction velocities (CVs) within the circuit. Slow conduction was defined as <30 cm/s.

RESULTS

Fifteen localized re-entrant circuits (12 left atrial, 3 right atrial) with mean cycle length 273 ± 40 ms were analyzed using high-resolution maps (22,389 ± 13,375 EGMs). A mean of 4.5 ± 1.6 slow conduction corridors were identified per circuit. Although the entire circuit was of low voltage, the bipolar voltage in slow conducting corridors was significantly lower than the rest of the circuit (0.22 ± 0.20 mV vs 0.50 ± 0.48 mV; P <.001). The mean conduction velocity of the circuit, excluding slow conduction areas, was 90.3 ± 34.3 cm/s vs 13.9 ± 3.5 cm/s (P <.001) in the slow conduction corridors. EGM analysis at the slowest conduction corridors demonstrated fractionation (100%) with longer EGM duration compared to the other slow conduction corridors along the circuit (99 ± 9 ms vs 74 ± 11 ms; P = .003).

CONCLUSION

In contrast to current understanding, localized atrial reentrant circuits have multiple sequential "corridors" of very slow conduction (2-7) that contribute to maintenance of arrhythmia. The localized reentry occurs in low-voltage areas, with voltage further reduced in these multiple slow conducting corridors.