Caroline H Roney

Techniques for automated local activation time annotation and conduction velocity estimation in cardiac mapping

Measurements of cardiac conduction velocity provide valuable functional and structural insight into the initiation and perpetuation of cardiac arrhythmias, in both a clinical and laboratory context. The interpretation of activation wavefronts and their propagation can identify mechanistic properties of a broad range of electrophysiological pathologies. However, the sparsity, distribution and uncertainty of recorded data make accurate conduction velocity calculation difficult. A wide range of mathematical approaches have been proposed for addressing this challenge, often targeted towards specific data modalities, species or recording environments. Many of these algorithms require identification of activation times from electrogram recordings which themselves may have complex morphology or low signal-to-noise ratio. This paper surveys algorithms designed for identifying local activation times and computing conduction direction and speed. Their suitability for use in different recording contexts and applications is assessed.

Spatial Resolution Requirements for Accurate Identification of Drivers of Atrial Fibrillation

Background— Recent studies have demonstrated conflicting mechanisms underlying atrial fibrillation (AF), with the spatial resolution of data often cited as a potential reason for the disagreement. The purpose of this study was to investigate whether the variation in spatial resolution of mapping may lead to misinterpretation of the underlying mechanism in persistent AF. Methods and Results— Simulations of rotors and focal sources were performed to estimate the minimum number of recording points required to correctly identify the underlying AF mechanism. The effects of different data types (action potentials and unipolar or bipolar electrograms) and rotor stability on resolution requirements were investigated. We also determined the ability of clinically used endocardial catheters to identify AF mechanisms using clinically recorded and simulated data. The spatial resolution required for correct identification of rotors and focal sources is a linear function of spatial wavelength (the distance between wavefronts) of the arrhythmia. Rotor localization errors are larger for electrogram data than for action potential data. Stationary rotors are more reliably identified compared with meandering trajectories, for any given spatial resolution. All clinical high-resolution multipolar catheters are of sufficient resolution to accurately detect and track rotors when placed over the rotor core although the low-resolution basket catheter is prone to false detections and may incorrectly identify rotors that are not present. Conclusions— The spatial resolution of AF data can significantly affect the interpretation of the underlying AF mechanism. Therefore, the interpretation of human AF data must be taken in the context of the spatial resolution of the recordings.

An automated algorithm for determining conduction velocity, wavefront direction and origin of focal cardiac arrhythmias using a multipolar catheter.

Determining locations of focal arrhythmia sources and quantifying myocardial conduction velocity (CV) are two major challenges in clinical catheter ablation cases. CV, wave-front direction and focal source location can be estimated from multipolar catheter data, but currently available methods are time-consuming, limited to specific electrode configurations, and can be inaccurate. We developed automated algorithms to rapidly identify CV from multipolar catheter data with any arrangement of electrodes, whilst providing estimates of wavefront direction and focal source position, which can guide the catheter towards a focal arrhythmic source. We validated our methods using simulations on realistic human left atrial geometry. We subsequently applied them to clinically-acquired intracardiac electrogram data, where CV and wavefront direction were accurately determined in all cases, whilst focal source locations were correctly identified in 2/3 cases. Our novel automated algorithms can potentially be used to guide ablation of focal arrhythmias in real-time in cardiac catheter laboratories.

ST-Elevation Magnitude Correlates With Right Ventricular Outflow Tract Conduction Delay in Type I Brugada ECG

Background: The substrate location and underlying electrophysiological mechanisms that contribute to the characteristic ECG pattern of Brugada syndrome (BrS) are still debated. Using noninvasive electrocardiographical imaging, we studied whole heart conduction and repolarization patterns during ajmaline challenge in BrS individuals. Methods and Results: A total of 13 participants (mean age, 44±12 years; 8 men), 11 concealed patients with type I BrS and 2 healthy controls, underwent an ajmaline infusion with electrocardiographical imaging and ECG recordings. Electrocardiographical imaging activation recovery intervals and activation timings across the right ventricle (RV) body, outflow tract (RVOT), and left ventricle were calculated and analyzed at baseline and when type I BrS pattern manifested after ajmaline infusion. Peak J-ST point elevation was calculated from the surface ECG and compared with the electrocardiographical imaging–derived parameters at the same time point. After ajmaline infusion, the RVOT had the greatest increase in conduction delay (5.4±2.8 versus 2.0±2.8 versus 1.1±1.6 ms; P=0.007) and activation recovery intervals prolongation (69±32 versus 39±29 versus 21±12 ms; P=0.0005) compared with RV or left ventricle. In controls, there was minimal change in J-ST point elevation, conduction delay, or activation recovery intervals at all sites with ajmaline. In patients with BrS, conduction delay in RVOT, but not RV or left ventricle, correlated to the degree of J-ST point elevation (Pearson R, 0.81; P<0.001). No correlation was found between J-ST point elevation and activation recovery intervals prolongation in the RVOT, RV, or left ventricle. Conclusions: Magnitude of ST (J point) elevation in the type I BrS pattern is attributed to degree of conduction delay in the RVOT and not prolongation in repolarization time.

A technique for visualising three-dimensional left atrial cardiac activation data in two dimensions with minimal distance distortion.

Electro-anatomic mapping and medical imaging systems, used during clinical procedures for treatment of atrial arrhythmias, frequently record and display measurements on an anatomical surface of the left atrium. As such, obtaining a complete picture of activation necessitates simultaneous views from multiple angles. In addition, post-processing of three-dimensional surface data is challenging, since algorithms are typically applicable to planar or volumetric data. We applied a surface flattening methodology to medical imaging data and electro-anatomic mapping data to generate a two-dimensional representation that best preserves distances, since the calculation of many clinically relevant metrics, including conduction velocity and rotor trajectory identification require an accurate representation of distance. Distance distortions were small and improved upon exclusion of the pulmonary veins. The technique is demonstrated using maps of local activation time, based on clinical data, and plotting rotor-core trajectories, using simulated data.

Influence of left atrial geometry on rotor core trajectories in a model of atrial fibrillation

Left atrial anatomy and myocardial architecture are known to influence rotor initiation and maintenance. However, identifying their relative contribution clinically is challenging. The present study aims to investigate in silico the effect of left atrial geometry in isolation on rotor generation and evolution through the spatiotemporal tracking of phase singularities. After meandering for a short period of time, rotors are attracted to specific areas of the chamber where there is high curvature, primarily near the base of the left atrial appendage and the junctions of the pulmonary veins. This suggests that the left atrial anatomy could play a key role in the perpetuation of fibrillatory activity.

Repolarization abnormalities unmasked with exercise in sudden cardiac death survivors with structurally normal hearts

Models of cardiac arrhythmogenesis predict that nonuniformity in repolarization and/or depolarization promotes ventricular fibrillation and is modulated by autonomic tone, but this is difficult to evaluate in patients. We hypothesize that such spatial heterogeneities would be detected by noninvasive ECG imaging (ECGi) in sudden cardiac death (SCD) survivors with structurally normal hearts under physiological stress.

Automated fiducial point selection for reducing registration error in the co-localisation of left atrium electroanatomic and imaging data

Registration of electroanatomic surfaces and segmented images for the co-localisation of structural and functional data typically requires the manual selection of fiducial points, which are used to initialise automated surface registration. The identification of equivalent points on geometric features by the human eye is heavily subjective, and error in their selection may lead to distortion of the transformed surface and subsequently limit the accuracy of data co-localisation. We propose that the manual trimming of the pulmonary veins through the region of greatest geometrical curvature, coupled with an automated angle-based fiducial-point selection algorithm, significantly reduces target registration error compared with direct manual selection of fiducial points.

A software platform for the comparative analysis of electroanatomic and imaging data including conduction velocity mapping

Electroanatomic mapping systems collect increasingly large quantities of spatially-distributed electrical data which may be potentially further scrutinized post-operatively to expose mechanistic properties which sustain and perpetuate atrial fibrillation. We describe a modular software platform, developed to post-process and rapidly analyse data exported from electroanatomic mapping systems using a range of existing and novel algorithms. Imaging data highlighting regions of scar can also be overlaid for comparison. In particular, we describe the conduction velocity (CV) mapping algorithm used to highlight wavefront behaviour. CV was found to be particularly sensitive to the spatial distribution of the triangulation points and corresponding activation times. A set of geometric conditions were devised for selecting suitable triangulations of the electrogram set for generating CV maps.

Variability in pulmonary vein electrophysiology and fibrosis determines arrhythmia susceptibility and dynamics

Success rates for catheter ablation of persistent atrial fibrillation patients are currently low; however, there is a subset of patients for whom electrical isolation of the pulmonary veins alone is a successful treatment strategy. It is difficult to identify these patients because there are a multitude of factors affecting arrhythmia susceptibility and maintenance, and the individual contributions of these factors are difficult to determine clinically. We hypothesised that the combination of pulmonary vein (PV) electrophysiology and atrial body fibrosis determine driver location and effectiveness of pulmonary vein isolation (PVI). We used bilayer biatrial computer models based on patient geometries to investigate the effects of PV properties and atrial fibrosis on arrhythmia inducibility, maintenance mechanisms, and the outcome of PVI. Short PV action potential duration (APD) increased arrhythmia susceptibility, while longer PV APD was found to be protective. Arrhythmia inducibility increased with slower conduction velocity (CV) at the LA/PV junction, but not for cases with homogeneous CV changes or slower CV at the distal PV. Phase singularity (PS) density in the PV region for cases with PV fibrosis was increased. Arrhythmia dynamics depend on both PV properties and fibrosis distribution, varying from meandering rotors to PV reentry (in cases with baseline or long APD), to stable rotors at regions of high fibrosis density. Measurement of fibrosis and PV properties may indicate patient specific susceptibility to AF initiation and maintenance. PV PS density before PVI was higher for cases in which AF terminated or converted to a macroreentry; thus, high PV PS density may indicate likelihood of PVI success.