Quantifying the impact of shape uncertainty on predicted arrhythmias.

BACKGROUND: Personalised computer models are increasingly used to diagnose cardiac arrhythmias and tailor treatment. Patient-specific models of the left atrium are often derived from pre-procedural imaging of anatomy and fibrosis. These images contain noise that can affect simulation predictions. There are few computationally tractable methods for propagating uncertainties from images to clinical predictions. METHOD: We describe the left atrium anatomy using our Bayesian shape model that captures anatomical uncertainty in medical images and has been validated on 63 independent clinical images. This algorithm describes the left atrium anatomy using Nmodes=15 principal components, capturing 95% of the shape variance and calculated from 70 clinical cardiac magnetic resonance (CMR) images. Latent variables encode shape uncertainty: we evaluate their posterior distribution for each new anatomy. We assume a normally distributed prior. We use the unscented transform to sample from the posterior shape distribution. For each sample, we assign the local material properties of the tissue using the projection of late gadolinium enhancement CMR (LGE-CMR) onto the anatomy to estimate local fibrosis. To test which activation patterns an atrium can sustain, we perform an arrhythmia simulation for each sample. We consider 34 possible outcomes (31 macro-re-entries, functional re-entry, atrial fibrillation, and non-sustained arrhythmia). For each sample, we determine the outcome by comparing pre- and post-ablation activation patterns following a cross-field stimulus. RESULTS: We create patient-specific atrial electrophysiology models of ten patients. We validate the mean and standard deviation maps from the unscented transform with the same statistics obtained with 12,000 Monte Carlo (ground truth) samples. We found discrepancies <3% and <2% for the mean and standard deviation for fibrosis burden and activation time, respectively. For each patient case, we then compare the predicted outcome from a model built on the clinical data (deterministic approach) with the probability distribution obtained from the simulated samples. We found that the deterministic approach did not predict the most likely outcome in 80% of the cases. Finally, we estimate the influence of each source of uncertainty independently. Fixing the anatomy to the posterior mean and maintaining uncertainty in fibrosis reduced the prediction of self-terminating arrhythmias from ≃14% to ≃7%. Keeping the fibrosis fixed to the sample mean while retaining uncertainty in shape decreased the prediction of substrate-driven arrhythmias from ≃33% to ≃18% and increased the prediction of macro-re-entries from ≃54% to ≃68%. CONCLUSIONS: We presented a novel method for propagating shape uncertainty in atrial models through to uncertainty in numerical simulations. The algorithm takes advantage of the unscented transform to compute the output distribution of the outcomes. We validated the unscented transform as a viable sampling strategy to deal with anatomy uncertainty. We then showed that the prediction computed with a deterministic model does not always coincide with the most likely outcome. Finally, we found that shape uncertainty affects the predictions of macro-re-entries, while fibrosis uncertainty affects the predictions of functional re-entries.

Provocation and localization of atrial ectopy in patients with atrial septal defects.

BACKGROUND: Atrial fibrillation (AF) is associated with atrial septal defects (ASDs), but the mechanism of arrhythmia in these patients is poorly understood. We hypothesised that right-sided atrial ectopy may predominate in this cohort. Here, we aimed to localise the origin of spontaneous and provoked atrial ectopy in ASD patients. METHODS: Following invasive calibration of P-wave axes, 24-h Holter monitoring was used to determine the chamber of origin of spontaneous atrial ectopy. Simultaneous electrogram recording from multiple intra-cardiac catheters was used to determine the chamber of origin of isoprenaline-provoked ectopy. Comparison was made to a group of non-congenital heart disease AF patients. RESULTS: Amongst ASD patients, a right-sided origin for spontaneous atrial ectopy was significantly more prevalent than a left-sided origin (24/30 patients with right-sided ectopy vs. 14/30 with left-sided ectopy, P = 0.015). Amongst AF patients, there was no difference in the prevalence of spontaneous right vs. left-sided ectopy. For isoprenaline-provoked ectopy, there was no significant difference in the proportions of patients with right-sided or left-sided ectopy in either group. CONCLUSIONS: When spontaneous atrial ectopy occurs in ASD patients, it is significantly more prevalent from a right-sided than left-sided origin. Isoprenaline infusion did not reveal the predilection for right-sided ectopy during electrophysiology study.

Correction: Linking statistical shape models and simulated function in the healthy adult human heart.

[This corrects the article DOI: 10.1371/journal.pcbi.1008851.].

Regional left ventricle scar detection from routine cardiac computed tomography angiograms using latent space classification.

OBJECTIVES: The aim of this study is to develop an automated method of regional scar detection on clinically standard computed tomography angiography (CTA) using encoder-decoder networks with latent space classification. BACKGROUND: Localising scar in cardiac patients can assist in diagnosis and guide interventions. Magnetic resonance imaging (MRI) with late gadolinium enhancement (LGE) is the clinical gold standard for scar imaging; however, it is commonly contraindicated. CTA is an alternative imaging modality that has fewer contraindications and is widely used as a first-line imaging modality of cardiac applications. METHODS: A dataset of 79 patients with both clinically indicated MRI LGE and subsequent CTA scans was used to train and validate networks to classify septal and lateral scar presence within short axis left ventricle slices. Two designs of encoder-decoder networks were compared, with one encoding anatomical shape in the latent space. Ground truth was established by segmenting scar in MRI LGE and registering this to the CTA images. Short axis slices were taken from the CTA, which served as the input to the networks. An independent external set of 22 cases (27% the size of the cross-validation set) was used to test the best network. RESULTS: A network classifying lateral scar only achieved an area under ROC curve of 0.75, with a sensitivity of 0.79 and specificity of 0.62 on the independent test set. The results of septal scar classification were poor (AUC < 0.6) for all networks. This was likely due to a high class imbalance. The highest AUC network encoded anatomical shape information in the network latent space, indicating it was important for the successful classification of lateral scar. CONCLUSIONS: Automatic lateral wall scar detection can be performed from a routine cardiac CTA with reasonable accuracy, without any scar specific imaging. This requires only a single acquisition in the cardiac cycle. In a clinical setting, this could be useful for pre-procedure planning, especially where MRI is contraindicated. Further work with more septal scar present is warranted to improve the usefulness of this approach.

CArdiac MagnEtic resonance assessment of bi-Atrial fibrosis in secundum atrial septal defects patients: CAMERA-ASD study.

AIMS: Atrial septal defects (ASD) are associated with atrial arrhythmias, but the arrhythmia substrate in these patients is poorly defined. We hypothesized that bi-atrial fibrosis is present and that right atrial fibrosis is associated with atrial arrhythmias in ASD patients. We aimed to evaluate the extent of bi-atrial fibrosis in ASD patients and to investigate the relationships between bi-atrial fibrosis, atrial arrhythmias, shunt fraction, and age. METHODS AND RESULTS: Patients with uncorrected secundum ASDs (n = 36; 50.4 ± 13.6 years) underwent cardiac magnetic resonance imaging with atrial late gadolinium enhancement. Comparison was made to non-congenital heart disease patients (n = 36; 60.3 ± 10.5 years) with paroxysmal atrial fibrillation (AF). Cardiac magnetic resonance parameters associated with atrial arrhythmias were identified and the relationship between bi-atrial structure, age, and shunt fraction studied. Bi-atrial fibrosis burden was greater in ASD patients than paroxysmal AF patients (20.7 ± 14% vs. 10.1 ± 8.6% and 14.8 ± 8.5% vs. 8.6 ± 6.1% for right and left atria respectively, P = 0.001 for both). In ASD patients, right atrial fibrosis burden was greater in those with than without atrial arrhythmias (33.4 ± 18.7% vs. 16.8 ± 10.3%, P = 0.034). On receiver operating characteristic analysis, a right atrial fibrosis burden of 32% had a 92% specificity and 71% sensitivity for predicting the presence of atrial arrhythmias. Neither age nor shunt fraction was associated with bi-atrial fibrosis burden. CONCLUSION: Bi-atrial fibrosis burden is greater in ASD patients than non-congenital heart disease patients with paroxysmal AF. Right atrial fibrosis is associated with the presence of atrial arrhythmias in ASD patients. These findings highlight the importance of right atrial fibrosis to atrial arrhythmogenesis in ASD patients.

Automated Left Ventricle Ischemic Scar Detection in CT Using Deep Neural Networks.

Objectives: The aim of this study is to develop a scar detection method for routine computed tomography angiography (CTA) imaging using deep convolutional neural networks (CNN), which relies solely on anatomical information as input and is compatible with existing clinical workflows. Background: Identifying cardiac patients with scar tissue is important for assisting diagnosis and guiding interventions. Late gadolinium enhancement (LGE) magnetic resonance imaging (MRI) is the gold standard for scar imaging; however, there are common instances where it is contraindicated. CTA is an alternative imaging modality that has fewer contraindications and is faster than Cardiovascular magnetic resonance imaging but is unable to reliably image scar. Methods: A dataset of LGE MRI (200 patients, 83 with scar) was used to train and validate a CNN to detect ischemic scar slices using segmentation masks as input to the network. MRIs were segmented to produce 3D left ventricle meshes, which were sampled at points along the short axis to extract anatomical masks, with scar labels from LGE as ground truth. The trained CNN was tested with an independent CTA dataset (25 patients, with ground truth established with paired LGE MRI). Automated segmentation was performed to provide the same input format of anatomical masks for the network. The CNN was compared against manual reading of the CTA dataset by 3 experts. Results: Note that 84.7% cross-validated accuracy (AUC: 0.896) for detecting scar slices in the left ventricle on the MRI data was achieved. The trained network was tested against the CTA-derived data, with no further training, where it achieved an 88.3% accuracy (AUC: 0.901). The automated pipeline outperformed the manual reading by clinicians. Conclusion: Automatic ischemic scar detection can be performed from a routine cardiac CTA, without any scar-specific imaging or contrast agents. This requires only a single acquisition in the cardiac cycle. In a clinical setting, with near zero additional cost, scar presence could be detected to triage images, reduce reading times, and guide clinical decision-making.

Linking statistical shape models and simulated function in the healthy adult human heart.

Cardiac anatomy plays a crucial role in determining cardiac function. However, there is a poor understanding of how specific and localised anatomical changes affect different cardiac functional outputs. In this work, we test the hypothesis that in a statistical shape model (SSM), the modes that are most relevant for describing anatomy are also most important for determining the output of cardiac electromechanics simulations. We made patient-specific four-chamber heart meshes (n = 20) from cardiac CT images in asymptomatic subjects and created a SSM from 19 cases. Nine modes captured 90% of the anatomical variation in the SSM. Functional simulation outputs correlated best with modes 2, 3 and 9 on average (R = 0.49 ± 0.17, 0.37 ± 0.23 and 0.34 ± 0.17 respectively). We performed a global sensitivity analysis to identify the different modes responsible for different simulated electrical and mechanical measures of cardiac function. Modes 2 and 9 were the most important for determining simulated left ventricular mechanics and pressure-derived phenotypes. Mode 2 explained 28.56 ± 16.48% and 25.5 ± 20.85, and mode 9 explained 12.1 ± 8.74% and 13.54 ± 16.91% of the variances of mechanics and pressure-derived phenotypes, respectively. Electrophysiological biomarkers were explained by the interaction of 3 ± 1 modes. In the healthy adult human heart, shape modes that explain large portions of anatomical variance do not explain equivalent levels of electromechanical functional variation. As a result, in cardiac models, representing patient anatomy using a limited number of modes of anatomical variation can cause a loss in accuracy of simulated electromechanical function.

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.

Percutaneous secundum atrial septal defect closure for the treatment of atrial arrhythmia in the adult: A meta-analysis.

BACKGROUND: Atrial arrhythmias are common in patients with atrial septal defects (ASD) but the effects of percutaneous closure on atrial arrhythmia prevalence is unclear. We investigated the effects of ASD device closure and the impact of age at time of closure on prevalent atrial arrythmia. METHODS: Meta-analysis of studies reporting atrial arrhythmia prevalence in adult patients before and after percutaneous closure was performed. Primary outcomes were prevalence of 'all atrial arrhythmia' and atrial fibrillation alone post closure. Sub-group analysis examined the effects of closure according to age in patients; <40 years, ≥40 and ≥ 60 years. 25 studies were included. RESULTS: Meta-analysis of all studies demonstrated no reduction in all atrial arrhythmia or atrial fibrillation prevalence post-closure (OR 0.855, 95% CI 0.672 to 1.087, P = .201 and OR 0.818, 95% CI 0.645 to 1.038, P = .099, respectively). A weak reduction in all atrial arrhythmia and atrial fibrillation was seen in patients ≥40 years (OR 0.77, 95% CI 0.616 to 0.979, P = .032 and OR 0.760, 95% CI 0.6 to 0.964, P = .024, respectively) but not ≥60 years (OR 0.822, 95% CI 0.593 to 1.141, P = .242 and OR 0.83, 95% CI 0.598 to 1.152, P = .266, respectively). No data were available in patients <40 years. This, and other limitations, prevents conclusive assessment of the effect of age on arrhythmia prevalence. CONCLUSIONS: Overall, percutaneous ASD closure is not associated with a reduction in atrial arrhythmia prevalence in this meta-analysis. A weak benefit is seen in patients ≥40 years of age, not present in patients ≥60 years.

Gaussian process manifold interpolation for probabilistic atrial activation maps and uncertain conduction velocity.

In patients with atrial fibrillation, local activation time (LAT) maps are routinely used for characterizing patient pathophysiology. The gradient of LAT maps can be used to calculate conduction velocity (CV), which directly relates to material conductivity and may provide an important measure of atrial substrate properties. Including uncertainty in CV calculations would help with interpreting the reliability of these measurements. Here, we build upon a recent insight into reduced-rank Gaussian processes (GPs) to perform probabilistic interpolation of uncertain LAT directly on human atrial manifolds. Our Gaussian process manifold interpolation (GPMI) method accounts for the topology of the atrium, and allows for calculation of statistics for predicted CV. We demonstrate our method on two clinical cases, and perform validation against a simulated ground truth. CV uncertainty depends on data density, wave propagation direction and CV magnitude. GPMI is suitable for probabilistic interpolation of other uncertain quantities on non-Euclidean manifolds. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Probabilistic Interpolation of Uncertain Local Activation Times on Human Atrial Manifolds.

OBJECTIVE: Local activation time (LAT) mapping of the atria is important for targeted treatment of atrial arrhythmias, but current methods do not interpolate on the atrial manifold and neglect uncertainties associated with LAT observations. In this paper, we describe novel methods to, first, quantify uncertainties in LAT arising from bipolar electrogram analysis and assignment of electrode recordings to the anatomical mesh, second, interpolate uncertain LAT measurements directly on left atrial manifolds to obtain complete probabilistic activation maps, and finally, interpolate LAT jointly across both the manifold and different S1-S2 pacing protocols. METHODS: A modified center of mass approach was used to process bipolar electrograms, yielding a LAT estimate and error distribution from the electrogram morphology. An error distribution for assigning measurements to the anatomical mesh was estimated. Probabilistic LAT maps were produced by interpolating on a left atrial manifold using Gaussian Markov random fields, taking into account observation errors and characterizing LAT predictions by their mean and standard deviation. This approach was extended to interpolate across S1-S2 pacing protocols. RESULTS: We evaluated our approach using recordings from three patients undergoing atrial ablation. Cross-validation showed consistent and accurate prediction of LAT observations both at different locations on the left atrium and for different S1-S2 intervals. SIGNIFICANCE: Interpolation of scalar and vector fields across anatomical structures from point measurements is a challenging problem in biomedical engineering, compounded by uncertainties in measurements and meshes. New methods and approaches are required, and in this paper, we have demonstrated an effective method for probabilistic interpolation of uncertain LAT.

The impact of wall thickness and curvature on wall stress in patient-specific electromechanical models of the left atrium.

The left atrium (LA) has a complex anatomy with heterogeneous wall thickness and curvature. The anatomy plays an important role in determining local wall stress; however, the relative contribution of wall thickness and curvature in determining wall stress in the LA is unknown. We have developed electromechanical finite element (FE) models of the LA using patient-specific anatomical FE meshes with rule-based myofiber directions. The models of the LA were passively inflated to 10mmHg followed by simulation of the contraction phase of the atrial cardiac cycle. The FE models predicted maximum LA volumes of 156.5 mL, 99.3 mL and 83.4 mL and ejection fractions of 36.9%, 32.0% and 25.2%. The median wall thickness in the 3 cases was calculated as [Formula: see text] mm, [Formula: see text] mm, and [Formula: see text] mm. The median curvature was determined as [Formula: see text] [Formula: see text], [Formula: see text], and [Formula: see text]. Following passive inflation, the correlation of wall stress with the inverse of wall thickness and curvature was 0.55-0.62 and 0.20-0.25, respectively. At peak contraction, the correlation of wall stress with the inverse of wall thickness and curvature was 0.38-0.44 and 0.16-0.34, respectively. In the LA, the 1st principal Cauchy stress is more dependent on wall thickness than curvature during passive inflation and both correlations decrease during active contraction. This emphasizes the importance of including the heterogeneous wall thickness in electromechanical FE simulations of the LA. Overall, simulation results and sensitivity analyses show that in complex atrial anatomy it is unlikely that a simple anatomical-based law can be used to estimate local wall stress, demonstrating the importance of FE analyses.

Pulmonary vein encirclement using an Ablation Index-guided point-by-point workflow: cardiovascular magnetic resonance assessment of left atrial scar formation.

AIMS: A point-by-point workflow for pulmonary vein isolation (PVI) targeting pre-defined Ablation Index values (a composite of contact force, time, and power) and minimizing interlesion distance may optimize the creation of contiguous ablation lesions whilst minimizing scar formation. We aimed to compare ablation scar formation in patients undergoing PVI using this workflow to patients undergoing a continuous catheter drag workflow. METHODS AND RESULTS: Post-ablation cardiovascular magnetic resonance imaging was performed in patients undergoing 1st-time PVI using a parameter-guided point-by-point workflow (n = 26). Total left atrial scar burden and the width and continuity of the pulmonary vein encirclement were determined on analysis of atrial late gadolinium enhancement sequences. Comparison was made with a cohort of patients (n = 20) undergoing PVI using continuous drag lesions. Mean post-ablation scar burden and scar width were significantly lower in the point-by-point group than in the control group (6.6 ± 6.8% vs. 9.6 ± 5.0%, P = 0.03 and 7.9 ± 3.6 mm vs. 10.7 ± 2.3 mm, P = 0.003). More complete bilateral pulmonary vein encirclements were seen in the point-by-point group (P = 0.038). All patients achieved acute PVI. CONCLUSION: Pulmonary vein isolation using a point-by-point workflow is feasible and results in a lower scar burden and scar width with more complete pulmonary vein encirclements than a conventional drag lesion approach.

Advances in Real-Time MRI-Guided Electrophysiology.

PURPOSE OF REVIEW: Theoretical benefits of real-time MRI guidance over conventional electrophysiology include contemporaneous 3D substrate assessment and accurate intra-procedural guidance and evaluation of ablation lesions. We review the unique challenges inherent to MRI-guided electrophysiology and how to translate the potential benefits in the treatment of cardiac arrhythmias. RECENT FINDINGS: Over the last 5 years, there has been substantial progress, initially in animal models and more recently in clinical studies, to establish methods and develop workflows within the MR environment that resemble those of conventional electrophysiology laboratories. Real-time MRI-guided systems have been used to perform electroanatomic mapping and ablation in patients with atrial flutter, and there is interest in developing the technology to tackle more complex arrhythmias including atrial fibrillation and ventricular tachycardia. SUMMARY: Mainstream adoption of real-time MRI-guided electrophysiology will require demonstration of clinical benefit and will be aided by increased availability of devices suitable for use in the MRI environment.

Evaluation of a real-time magnetic resonance imaging-guided electrophysiology system for structural and electrophysiological ventricular tachycardia substrate assessment.

AIMS: Potential advantages of real-time magnetic resonance imaging (MRI)-guided electrophysiology (MR-EP) include contemporaneous three-dimensional substrate assessment at the time of intervention, improved procedural guidance, and ablation lesion assessment. We evaluated a novel real-time MR-EP system to perform endocardial voltage mapping and assessment of delayed conduction in a porcine ischaemia-reperfusion model. METHODS AND RESULTS: Sites of low voltage and slow conduction identified using the system were registered and compared to regions of late gadolinium enhancement (LGE) on MRI. The Sorensen-Dice similarity coefficient (DSC) between LGE scar maps and voltage maps was computed on a nodal basis. A total of 445 electrograms were recorded in sinus rhythm (range: 30-186) using the MR-EP system including 138 electrograms from LGE regions. Pacing captured at 103 sites; 47 (45.6%) sites had a stimulus-to-QRS (S-QRS) delay of ≥40 ms. Using conventional (0.5-1.5 mV) bipolar voltage thresholds, the sensitivity and specificity of voltage mapping using the MR-EP system to identify MR-derived LGE was 57% and 96%, respectively. Voltage mapping had a better predictive ability in detecting LGE compared to S-QRS measurements using this system (area under curve: 0.907 vs. 0.840). Using an electrical threshold of 1.5 mV to define abnormal myocardium, the total DSC, scar DSC, and normal myocardium DSC between voltage maps and LGE scar maps was 79.0 ± 6.0%, 35.0 ± 10.1%, and 90.4 ± 8.6%, respectively. CONCLUSION: Low-voltage zones and regions of delayed conduction determined using a real-time MR-EP system are moderately associated with LGE areas identified on MRI.

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.

Mind the gap: Quantification of incomplete ablation patterns after pulmonary vein isolation using minimum path search.

Pulmonary vein isolation (PVI) is a common procedure for the treatment of atrial fibrillation (AF) since the initial trigger for AF frequently originates in the pulmonary veins. A successful isolation produces a continuous lesion (scar) completely encircling the veins that stops activation waves from propagating to the atrial body. Unfortunately, the encircling lesion is often incomplete, becoming a combination of scar and gaps of healthy tissue. These gaps are potential causes of AF recurrence, which requires a redo of the isolation procedure. Late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) is a non-invasive method that may also be used to detect gaps, but it is currently a time-consuming process, prone to high inter-observer variability. In this paper, we present a method to semi-automatically identify and quantify ablation gaps. Gap quantification is performed through minimum path search in a graph where every node is a scar patch and the edges are the geodesic distances between patches. We propose the Relative Gap Measure (RGM) to estimate the percentage of gap around a vein, which is defined as the ratio of the overall gap length and the total length of the path that encircles the vein. Additionally, an advanced version of the RGM has been developed to integrate gap quantification estimates from different scar segmentation techniques into a single figure-of-merit. Population-based statistical and regional analysis of gap distribution was performed using a standardised parcellation of the left atrium. We have evaluated our method on synthetic and clinical data from 50 AF patients who underwent PVI with radiofrequency ablation. The population-based analysis concluded that the left superior PV is more prone to lesion gaps while the left inferior PV tends to have less gaps (p < .05 in both cases), in the processed data. This type of information can be very useful for the optimization and objective assessment of PVI interventions.

Atrial Fibrillation Ablation in Patients with Heart Failure: One Size Does Not Fit All.

Atrial fibrillation (AF) is common in patients with heart failure and is associated with poorer clinical outcomes compared with patients with heart failure alone. Recent evidence has challenged previous treatment paradigms in which rate control was considered equivalent to rhythm control in this population. Catheter ablation has emerged as a safe and effective treatment strategy in selected patients and overcomes the issues of limited efficacy and drug toxicities associated with pharmacological rhythm control. Numerous studies have explored the benefits of catheter ablation in patients with heart failure, but these have included heterogeneous patient cohorts and variable ablation strategies. This state-of-the-art review explores the evidence from these trials and examines the need for tailored, patient-specific strategies for AF ablation in patients with heart failure.