Modelling methodology of atrial fibrosis affects rotor dynamics and electrograms.

AIMS: Catheter ablation is an effective technique for terminating atrial arrhythmia. However, given a high atrial fibrillation (AF) recurrence rate, optimal ablation strategies have yet to be defined. Computer modelling can be a powerful aid but modelling of fibrosis, a major factor associated with AF, is an open question. Several groups have proposed methodologies based on imaging data, but no comparison to determine which methodology best corroborates clinically observed reentrant behaviour has been performed. We examined several methodologies to determine the best method for capturing fibrillation dynamics. METHODS AND RESULTS: Patient late gadolinium-enhanced magnetic resonance imaging data were transferred onto a bilayer atrial computer model and used to assign fibrosis distributions. Fibrosis was modelled as conduction disturbances (lower conductivity, edge splitting, or percolation), transforming growth factor-ß1 ionic channel effects, myocyte-fibroblast coupling, and combinations of the preceding. Reentry was induced through pulmonary vein ectopy and the ensuing rotor dynamics characterized. Non-invasive electrocardiographic imaging data of the patients in AF was used for comparison. Electrograms were computed and the fractionation durations measured over the surface. Edge splitting produced more phase singularities from wavebreaks than the other representations. The number of phase singularities seen with percolation was closer to the clinical values. Addition of fibroblast coupling had an organizing effect on rotor dynamics. Simple tissue conductivity changes with ionic changes localized rotors over fibrosis which was not observed with clinical data. CONCLUSION: The specific representation of fibrosis has a large effect on rotor dynamics and needs to be carefully considered for patient specific modelling.

In silico Comparison of Left Atrial Ablation Techniques That Target the Anatomical, Structural, and Electrical Substrates of Atrial Fibrillation.

Catheter ablation therapy for persistent atrial fibrillation (AF) typically includes pulmonary vein isolation (PVI) and may include additional ablation lesions that target patient-specific anatomical, electrical, or structural features. Clinical centers employ different ablation strategies, which use imaging data together with electroanatomic mapping data, depending on data availability. The aim of this study was to compare ablation techniques across a virtual cohort of AF patients. We constructed 20 paroxysmal and 30 persistent AF patient-specific left atrial (LA) bilayer models incorporating fibrotic remodeling from late-gadolinium enhancement (LGE) MRI scans. AF was simulated and post-processed using phase mapping to determine electrical driver locations over 15 s. Six different ablation approaches were tested: (i) PVI alone, modeled as wide-area encirclement of the pulmonary veins; PVI together with: (ii) roof and inferior lines to model posterior wall box isolation; (iii) isolating the largest fibrotic area (identified by LGE-MRI); (iv) isolating all fibrotic areas; (v) isolating the largest driver hotspot region [identified as high simulated phase singularity (PS) density]; and (vi) isolating all driver hotspot regions. Ablation efficacy was assessed to predict optimal ablation therapies for individual patients. We subsequently trained a random forest classifier to predict ablation response using (a) imaging metrics alone, (b) imaging and electrical metrics, or (c) imaging, electrical, and ablation lesion metrics. The optimal ablation approach resulting in termination, or if not possible atrial tachycardia (AT), varied among the virtual patient cohort: (i) 20% PVI alone, (ii) 6% box ablation, (iii) 2% largest fibrosis area, (iv) 4% all fibrosis areas, (v) 2% largest driver hotspot, and (vi) 46% all driver hotspots. Around 20% of cases remained in AF for all ablation strategies. The addition of patient-specific and ablation pattern specific lesion metrics to the trained random forest classifier improved predictive capability from an accuracy of 0.73 to 0.83. The trained classifier results demonstrate that the surface areas of pre-ablation driver regions and of fibrotic tissue not isolated by the proposed ablation strategy are both important for predicting ablation outcome. Overall, our study demonstrates the need to select the optimal ablation strategy for each patient. It suggests that both patient-specific fibrosis properties and driver locations are important for planning ablation approaches, and the distribution of lesions is important for predicting an acute response.

Wavelength and Fibrosis Affect Phase Singularity Locations During Atrial Fibrillation.

The mechanisms underlying atrial fibrillation (AF), the most common sustained cardiac rhythm disturbance, remain elusive. Atrial fibrosis plays an important role in the development of AF and rotor dynamics. Both electrical wavelength (WL) and the degree of atrial fibrosis change as AF progresses. However, their combined effect on rotor core location remains unknown. The aim of this study was to analyze the effects of WL change on rotor core location in both fibrotic and non-fibrotic atria. Three patient specific fibrosis distributions (total fibrosis content: 16.6, 22.8, and 19.2%) obtained from clinical imaging data of persistent AF patients were incorporated in a bilayer atrial computational model. Fibrotic effects were modeled as myocyte-fibroblast coupling + conductivity remodeling; structural remodeling; ionic current changes + conductivity remodeling; and combinations of these methods. To change WL, action potential duration (APD) was varied from 120 to 240ms, representing the range of clinically observed AF cycle length, by modifying the inward rectifier potassium current (IK1) conductance between 80 and 140% of the original value. Phase singularities (PSs) were computed to identify rotor core locations. Our results show that IK1 conductance variation resulted in a decrease of APD and WL across the atria. For large WL in the absence of fibrosis, PSs anchored to regions with high APD gradient at the center of the left atrium (LA) anterior wall and near the junctions of the inferior pulmonary veins (PVs) with the LA. Decreasing the WL induced more PSs, whose distribution became less clustered. With fibrosis, PS locations depended on the fibrosis distribution and the fibrosis implementation method. The proportion of PSs in fibrotic areas and along the borders varied with both WL and fibrosis modeling method: for patient one, this was 4.2-14.9% as IK1 varied for the structural remodeling representation, but 12.3-88.4% using the combination of structural remodeling with myocyte-fibroblast coupling. The degree and distribution of fibrosis and the choice of implementation technique had a larger effect on PS locations than the WL variation. Thus, distinguishing the fibrotic mechanisms present in a patient is important for interpreting clinical fibrosis maps to create personalized models.

Rotor Tracking Using Phase of Electrograms Recorded During Atrial Fibrillation.

Extracellular electrograms recorded during atrial fibrillation (AF) are challenging to interpret due to the inherent beat-to-beat variability in amplitude and duration. Phase mapping represents these voltage signals in terms of relative position within the cycle, and has been widely applied to action potential and unipolar electrogram data of myocardial fibrillation. To date, however, it has not been applied to bipolar recordings, which are commonly acquired clinically. The purpose of this study is to present a novel algorithm for calculating phase from both unipolar and bipolar electrograms recorded during AF. A sequence of signal filters and processing steps are used to calculate phase from simulated, experimental, and clinical, unipolar and bipolar electrograms. The algorithm is validated against action potential phase using simulated data (trajectory centre error <0.8 mm); between experimental multi-electrode array unipolar and bipolar phase; and for wavefront identification in clinical atrial tachycardia. For clinical AF, similar rotational content (R 2 = 0.79) and propagation maps (median correlation 0.73) were measured using either unipolar or bipolar recordings. The algorithm is robust, uses standard signal processing techniques, and accurately quantifies AF wavefronts and sources. Identifying critical sources, such as rotors, in AF, may allow for more accurate targeting of ablation therapy and improved patient outcomes.

Novel Radiofrequency Ablation Strategies for Terminating Atrial Fibrillation in the Left Atrium: A Simulation Study.

Pulmonary vein isolation (PVI) with radiofrequency ablation (RFA) is the cornerstone of atrial fibrillation (AF) therapy, but few strategies exist for when it fails. To guide RFA, phase singularity (PS) mapping locates reentrant electrical waves (rotors) that perpetuate AF. The goal of this study was to test existing and develop new RFA strategies for terminating rotors identified with PS mapping. It is unsafe to test experimental RFA strategies in patients, so they were evaluated in silico using a bilayer computer model of the human atria with persistent AF (pAF) electrical (ionic) and structural (fibrosis) remodeling. pAF was initiated by rapidly pacing the right (RSPV) and left (LSPV) superior pulmonary veins during sinus rhythm, and rotor dynamics quantified by PS analysis. Three RFA strategies were studied: (i) PVI, roof, and mitral lines; (ii) circles, perforated circles, lines, and crosses 0.5-1.5 cm in diameter/length administered near rotor locations/pathways identified by PS mapping; and (iii) 4-8 lines streamlining the sequence of electrical activation during sinus rhythm. As in pAF patients, 2 ± 1 rotors with cycle length 185 ± 4 ms and short PS duration 452 ± 401 ms perpetuated simulated pAF. Spatially, PS density had weak to moderate positive correlations with fibrosis density (RSPV: r = 0.38, p = 0.35, LSPV: r = 0.77, p = 0.02). RFA PVI, mitral, and roof lines failed to terminate pAF, but RFA perforated circles and lines 1.5 cm in diameter/length terminated meandering rotors from RSPV pacing when placed at locations with high PS density. Similarly, RFA circles, perforated circles, and crosses 1.5 cm in diameter/length terminated stationary rotors from LSPV pacing. The most effective strategy for terminating pAF was to streamline the sequence of activation during sinus rhythm with >4 RFA lines. These results demonstrate that co-localizing 1.5 cm RFA lesions with locations of high PS density is a promising strategy for terminating pAF rotors. For patients immune to PVI, roof, mitral, and PS guided RFA strategies, streamlining patient-specific activation sequences during sinus rhythm is a robust but challenging alternative.