1
|
Jaffery OA, Melki L, Slabaugh G, Good WW, Roney CH. A Review of Personalised Cardiac Computational Modelling Using Electroanatomical Mapping Data. Arrhythm Electrophysiol Rev 2024; 13:e08. [PMID: 38807744 PMCID: PMC11131150 DOI: 10.15420/aer.2023.25] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 12/27/2023] [Indexed: 05/30/2024] Open
Abstract
Computational models of cardiac electrophysiology have gradually matured during the past few decades and are now being personalised to provide patient-specific therapy guidance for improving suboptimal treatment outcomes. The predictive features of these personalised electrophysiology models hold the promise of providing optimal treatment planning, which is currently limited in the clinic owing to reliance on a population-based or average patient approach. The generation of a personalised electrophysiology model entails a sequence of steps for which a range of activation mapping, calibration methods and therapy simulation pipelines have been suggested. However, the optimal methods that can potentially constitute a clinically relevant in silico treatment are still being investigated and face limitations, such as uncertainty of electroanatomical data recordings, generation and calibration of models within clinical timelines and requirements to validate or benchmark the recovered tissue parameters. This paper is aimed at reporting techniques on the personalisation of cardiac computational models, with a focus on calibrating cardiac tissue conductivity based on electroanatomical mapping data.
Collapse
Affiliation(s)
- Ovais A Jaffery
- School of Engineering and Materials Science, Queen Mary University of London London, UK
| | - Lea Melki
- R&D Algorithms, Acutus Medical Carlsbad, CA, US
| | - Gregory Slabaugh
- Digital Environment Research Institute, Queen Mary University of London London, UK
| | | | - Caroline H Roney
- School of Engineering and Materials Science, Queen Mary University of London London, UK
| |
Collapse
|
2
|
Strocchi M, Rodero C, Roney CH, Mendonca Costa C, Plank G, Lamata P, Niederer SA. A Semi-automatic Pipeline for Generation of Large Cohorts of Four-Chamber Heart Meshes. Methods Mol Biol 2024; 2735:117-127. [PMID: 38038846 DOI: 10.1007/978-1-0716-3527-8_7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2023]
Abstract
Computational models for cardiac electro-mechanics have been increasingly used to further understand heart function. Small cohort and single patient computational studies provide useful insight into cardiac pathophysiology and response to therapy. However, these smaller studies have limited capability to capture the high level of anatomical variability seen in cardiology patients. Larger cohort studies are, on the other hand, more representative of the study population, but building several patient-specific anatomical meshes can be time-consuming and requires access to larger datasets of imaging data, image processing software to label anatomical structures and tools to create high fidelity anatomical meshes. Limited access to these tools and data might limit advances in this area of research. In this chapter, we present our semi-automatic pipeline to build patient-specific four-chamber heart meshes from CT imaging datasets, including ventricular myofibers and a set of universal ventricular and atrial coordinates. This pipeline was applied to CT images from both heart failure patients and healthy controls to generate cohorts of tetrahedral meshes suitable for electro-mechanics simulations. Both cohorts were made publicly available in order to promote computational studies employing large virtual cohorts.
Collapse
Affiliation(s)
- Marina Strocchi
- Department of Biomedical Engineering, King's College London, London, UK
| | - Cristobal Rodero
- Department of Biomedical Engineering, King's College London, London, UK
| | - Caroline H Roney
- Department of Biomedical Engineering, King's College London, London, UK
| | | | | | - Pablo Lamata
- Department of Biomedical Engineering, King's College London, London, UK
| | - Steven A Niederer
- Department of Biomedical Engineering, King's College London, London, UK.
| |
Collapse
|
3
|
Bifulco SF, Boyle PM. Computational Modeling and Simulation of the Fibrotic Human Atria. Methods Mol Biol 2024; 2735:105-115. [PMID: 38038845 DOI: 10.1007/978-1-0716-3527-8_6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2023]
Abstract
Patient-specific modeling of atrial electrical activity enables the execution of simulations that can provide mechanistic insights and provide novel solutions to vexing clinical problems. The geometry and fibrotic remodeling of the heart can be reconstructed from clinical-grade medical scans and used to inform personalized models with detail incorporated at the cell- and tissue-scale to represent changes in image-identified diseased regions. Here, we provide a rubric for the reconstruction of realistic atrial models from pre-segmented 3D renderings of the left atrium with fibrotic tissue regions delineated, which are the output from clinical-grade systems for quantifying fibrosis. We then provide a roadmap for using those models to carry out patient-specific characterization of the fibrotic substrate in terms of its potential to harbor reentrant drivers via cardiac electrophysiology simulations.
Collapse
Affiliation(s)
- Savannah F Bifulco
- Department of Bioengineering, University of Washington, Seattle, WA, USA
| | - Patrick M Boyle
- Department of Bioengineering, University of Washington, Seattle, WA, USA.
- Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, WA, USA.
- Center for Cardiovascular Biology, University of Washington, Seattle, WA, USA.
| |
Collapse
|
4
|
Roney CH, Solis Lemus JA, Lopez Barrera C, Zolotarev A, Ulgen O, Kerfoot E, Bevis L, Misghina S, Vidal Horrach C, Jaffery OA, Ehnesh M, Rodero C, Dharmaprani D, Ríos-Muñoz GR, Ganesan A, Good WW, Neic A, Plank G, Hopman LHGA, Götte MJW, Honarbakhsh S, Narayan SM, Vigmond E, Niederer S. Constructing bilayer and volumetric atrial models at scale. Interface Focus 2023; 13:20230038. [PMID: 38106921 PMCID: PMC10722212 DOI: 10.1098/rsfs.2023.0038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Accepted: 11/15/2023] [Indexed: 12/19/2023] Open
Abstract
To enable large in silico trials and personalized model predictions on clinical timescales, it is imperative that models can be constructed quickly and reproducibly. First, we aimed to overcome the challenges of constructing cardiac models at scale through developing a robust, open-source pipeline for bilayer and volumetric atrial models. Second, we aimed to investigate the effects of fibres, fibrosis and model representation on fibrillatory dynamics. To construct bilayer and volumetric models, we extended our previously developed coordinate system to incorporate transmurality, atrial regions and fibres (rule-based or data driven diffusion tensor magnetic resonance imaging (MRI)). We created a cohort of 1000 biatrial bilayer and volumetric models derived from computed tomography (CT) data, as well as models from MRI, and electroanatomical mapping. Fibrillatory dynamics diverged between bilayer and volumetric simulations across the CT cohort (correlation coefficient for phase singularity maps: left atrial (LA) 0.27 ± 0.19, right atrial (RA) 0.41 ± 0.14). Adding fibrotic remodelling stabilized re-entries and reduced the impact of model type (LA: 0.52 ± 0.20, RA: 0.36 ± 0.18). The choice of fibre field has a small effect on paced activation data (less than 12 ms), but a larger effect on fibrillatory dynamics. Overall, we developed an open-source user-friendly pipeline for generating atrial models from imaging or electroanatomical mapping data enabling in silico clinical trials at scale (https://github.com/pcmlab/atrialmtk).
Collapse
Affiliation(s)
- Caroline H. Roney
- School of Engineering and Materials Science, Queen Mary University of London, London, UK
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, UK
| | - Jose Alonso Solis Lemus
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, UK
- National Heart and Lung Institute, Imperial College London, London, UK
| | - Carlos Lopez Barrera
- School of Engineering and Materials Science, Queen Mary University of London, London, UK
- Center for Research in Advanced Materials S.C (CIMAV), Chihuahua, Mexico
| | - Alexander Zolotarev
- School of Engineering and Materials Science, Queen Mary University of London, London, UK
| | - Onur Ulgen
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, UK
| | - Eric Kerfoot
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, UK
| | - Laura Bevis
- School of Engineering and Materials Science, Queen Mary University of London, London, UK
| | - Semhar Misghina
- School of Engineering and Materials Science, Queen Mary University of London, London, UK
| | - Caterina Vidal Horrach
- School of Engineering and Materials Science, Queen Mary University of London, London, UK
| | - Ovais A. Jaffery
- School of Engineering and Materials Science, Queen Mary University of London, London, UK
| | - Mahmoud Ehnesh
- School of Engineering and Materials Science, Queen Mary University of London, London, UK
| | - Cristobal Rodero
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, UK
- National Heart and Lung Institute, Imperial College London, London, UK
| | - Dhani Dharmaprani
- College of Medicine and Public Health, Flinders University, Adelaide, Australia
| | - Gonzalo R. Ríos-Muñoz
- Bioengineering Department, Universidad Carlos III de Madrid, Madrid 28911, Spain
- Department of Cardiology, Gregorio Marañón Health Research Institute (IiSGM), Hospital General Universitario Gregorio Marañón, Madrid 28007, Spain
- Center for Biomedical Research in Cardiovascular Disease Network (CIBERCV), Madrid 28029, Spain
| | - Anand Ganesan
- College of Medicine and Public Health, Flinders University, Adelaide, Australia
| | | | | | - Gernot Plank
- Gottfried Schatz Research Center-Biophysics, Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | | | | | - Shohreh Honarbakhsh
- Electrophysiology Department, Barts Heart Centre, Barts Health NHS Trust, London, UK
| | - Sanjiv M. Narayan
- Department of Medicine and Cardiovascular Institute, Stanford University, Palo Alto, CA, USA
| | - Edward Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Bordeaux, France
- IMB, UMR 5251, University Bordeaux, Talence 33400, France
| | - Steven Niederer
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, UK
- National Heart and Lung Institute, Imperial College London, London, UK
- Turing Research and Innovation Cluster in Digital Twins (TRIC: DT), The Alan Turing Institute, London, UK
| |
Collapse
|
5
|
Rodero C, Baptiste TMG, Barrows RK, Lewalle A, Niederer SA, Strocchi M. Advancing clinical translation of cardiac biomechanics models: a comprehensive review, applications and future pathways. FRONTIERS IN PHYSICS 2023; 11:1306210. [PMID: 38500690 PMCID: PMC7615748 DOI: 10.3389/fphy.2023.1306210] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
Cardiac mechanics models are developed to represent a high level of detail, including refined anatomies, accurate cell mechanics models, and platforms to link microscale physiology to whole-organ function. However, cardiac biomechanics models still have limited clinical translation. In this review, we provide a picture of cardiac mechanics models, focusing on their clinical translation. We review the main experimental and clinical data used in cardiac models, as well as the steps followed in the literature to generate anatomical meshes ready for simulations. We describe the main models in active and passive mechanics and the different lumped parameter models to represent the circulatory system. Lastly, we provide a summary of the state-of-the-art in terms of ventricular, atrial, and four-chamber cardiac biomechanics models. We discuss the steps that may facilitate clinical translation of the biomechanics models we describe. A well-established software to simulate cardiac biomechanics is lacking, with all available platforms involving different levels of documentation, learning curves, accessibility, and cost. Furthermore, there is no regulatory framework that clearly outlines the verification and validation requirements a model has to satisfy in order to be reliably used in applications. Finally, better integration with increasingly rich clinical and/or experimental datasets as well as machine learning techniques to reduce computational costs might increase model reliability at feasible resources. Cardiac biomechanics models provide excellent opportunities to be integrated into clinical workflows, but more refinement and careful validation against clinical data are needed to improve their credibility. In addition, in each context of use, model complexity must be balanced with the associated high computational cost of running these models.
Collapse
Affiliation(s)
- Cristobal Rodero
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
| | - Tiffany M. G. Baptiste
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, Faculty of Life Sciences and Medicine, King’s College London, London, United Kingdom
| | - Rosie K. Barrows
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, Faculty of Life Sciences and Medicine, King’s College London, London, United Kingdom
| | - Alexandre Lewalle
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
| | - Steven A. Niederer
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
- Turing Research and Innovation Cluster in Digital Twins (TRIC: DT), The Alan Turing Institute, London, United Kingdom
| | - Marina Strocchi
- Cardiac Electro-Mechanics Research Group (CEMRG), National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London, United Kingdom
| |
Collapse
|
6
|
Bifulco SF, Macheret F, Scott GD, Akoum N, Boyle PM. Explainable Machine Learning to Predict Anchored Reentry Substrate Created by Persistent Atrial Fibrillation Ablation in Computational Models. J Am Heart Assoc 2023; 12:e030500. [PMID: 37581387 PMCID: PMC10492949 DOI: 10.1161/jaha.123.030500] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/16/2023] [Accepted: 07/21/2023] [Indexed: 08/16/2023]
Abstract
Background Postablation arrhythmia recurrence occurs in ~40% of patients with persistent atrial fibrillation. Fibrotic remodeling exacerbates arrhythmic activity in persistent atrial fibrillation and can play a key role in reentrant arrhythmia, but emergent interaction between nonconductive ablation-induced scar and native fibrosis (ie, residual fibrosis) is poorly understood. Methods and Results We conducted computational simulations in pre- and postablation left atrial models reconstructed from late gadolinium enhanced magnetic resonance imaging scans to test the hypothesis that ablation in patients with persistent atrial fibrillation creates new substrate conducive to recurrent arrhythmia mediated by anchored reentry. We trained a random forest machine learning classifier to accurately pinpoint specific nonconductive tissue regions (ie, areas of ablation-delivered scar or vein/valve boundaries) with the capacity to serve as substrate for anchored reentry-driven recurrent arrhythmia (area under the curve: 0.91±0.03). Our analysis suggests there is a distinctive nonconductive tissue pattern prone to serving as arrhythmogenic substrate in postablation models, defined by a specific size and proximity to residual fibrosis. Conclusions Overall, this suggests persistent atrial fibrillation ablation transforms substrate that favors functional reentry (ie, rotors meandering in excitable tissue) into an arrhythmogenic milieu more conducive to anchored reentry. Our work also indicates that explainable machine learning and computational simulations can be combined to effectively probe mechanisms of recurrent arrhythmia.
Collapse
Affiliation(s)
| | - Fima Macheret
- Division of CardiologyUniversity of WashingtonSeattleWAUSA
| | - Griffin D. Scott
- Department of BioengineeringUniversity of WashingtonSeattleWAUSA
| | - Nazem Akoum
- Department of BioengineeringUniversity of WashingtonSeattleWAUSA
- Division of CardiologyUniversity of WashingtonSeattleWAUSA
| | - Patrick M. Boyle
- Department of BioengineeringUniversity of WashingtonSeattleWAUSA
- Institute for Stem Cell and Regenerative MedicineUniversity of WashingtonSeattleWAUSA
- Center for Cardiovascular BiologyUniversity of WashingtonSeattleWAUSA
| |
Collapse
|
7
|
Strocchi M, Longobardi S, Augustin CM, Gsell MAF, Petras A, Rinaldi CA, Vigmond EJ, Plank G, Oates CJ, Wilkinson RD, Niederer SA. Cell to whole organ global sensitivity analysis on a four-chamber heart electromechanics model using Gaussian processes emulators. PLoS Comput Biol 2023; 19:e1011257. [PMID: 37363928 DOI: 10.1371/journal.pcbi.1011257] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Accepted: 06/09/2023] [Indexed: 06/28/2023] Open
Abstract
Cardiac pump function arises from a series of highly orchestrated events across multiple scales. Computational electromechanics can encode these events in physics-constrained models. However, the large number of parameters in these models has made the systematic study of the link between cellular, tissue, and organ scale parameters to whole heart physiology challenging. A patient-specific anatomical heart model, or digital twin, was created. Cellular ionic dynamics and contraction were simulated with the Courtemanche-Land and the ToR-ORd-Land models for the atria and the ventricles, respectively. Whole heart contraction was coupled with the circulatory system, simulated with CircAdapt, while accounting for the effect of the pericardium on cardiac motion. The four-chamber electromechanics framework resulted in 117 parameters of interest. The model was broken into five hierarchical sub-models: tissue electrophysiology, ToR-ORd-Land model, Courtemanche-Land model, passive mechanics and CircAdapt. For each sub-model, we trained Gaussian processes emulators (GPEs) that were then used to perform a global sensitivity analysis (GSA) to retain parameters explaining 90% of the total sensitivity for subsequent analysis. We identified 45 out of 117 parameters that were important for whole heart function. We performed a GSA over these 45 parameters and identified the systemic and pulmonary peripheral resistance as being critical parameters for a wide range of volumetric and hemodynamic cardiac indexes across all four chambers. We have shown that GPEs provide a robust method for mapping between cellular properties and clinical measurements. This could be applied to identify parameters that can be calibrated in patient-specific models or digital twins, and to link cellular function to clinical indexes.
Collapse
Affiliation(s)
- Marina Strocchi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - Stefano Longobardi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | | | | | - Argyrios Petras
- Johann Radon Institute for Computational and Applied Mathematics (RICAM), Linz, Austria
| | - Christopher A Rinaldi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- Guy's and St Thomas' NHS Foundation Trust, London, United Kingdom
| | - Edward J Vigmond
- University of Bordeaux, CNRS, Bordeaux, Talence, France
- IHU Liryc, Bordeaux, Talence, France
| | - Gernot Plank
- Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | - Chris J Oates
- Newcastle University, Newcastle upon Tyne, United Kingdom
| | | | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
- Alan Turing Institute, London, United Kingdom
| |
Collapse
|
8
|
Rossi S, Abdala L, Woodward A, Vavalle JP, Henriquez CS, Griffith BE. Rule-based definition of muscle bundles in patient-specific models of the left atrium. Front Physiol 2022; 13:912947. [PMID: 36311246 PMCID: PMC9597256 DOI: 10.3389/fphys.2022.912947] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Accepted: 06/29/2022] [Indexed: 11/13/2022] Open
Abstract
Atrial fibrillation (AF) is the most common arrhythmia encountered clinically, and as the population ages, its prevalence is increasing. Although the CHA2DS2- VASc score is the most used risk-stratification system for stroke risk in AF, it lacks personalization. Patient-specific computer models of the atria can facilitate personalized risk assessment and treatment planning. However, a challenge faced in creating such models is the complexity of the atrial muscle arrangement and its influence on the atrial fiber architecture. This work proposes a semi-automated rule-based algorithm to generate the local fiber orientation in the left atrium (LA). We use the solutions of several harmonic equations to decompose the LA anatomy into subregions. Solution gradients define a two-layer fiber field in each subregion. The robustness of our approach is demonstrated by recreating the fiber orientation on nine models of the LA obtained from AF patients who underwent WATCHMAN device implantation. This cohort of patients encompasses a variety of morphology variants of the left atrium, both in terms of the left atrial appendages (LAAs) and the number of pulmonary veins (PVs). We test the fiber construction algorithm by performing electrophysiology (EP) simulations. Furthermore, this study is the first to compare its results with other rule-based algorithms for the LA fiber architecture definition available in the literature. This analysis suggests that a multi-layer fiber architecture is important to capture complex electrical activation patterns. A notable advantage of our approach is the ability to reconstruct the main LA fiber bundles in a variety of morphologies while solving for a small number of harmonic fields, leading to a comparatively straightforward and reproducible approach.
Collapse
Affiliation(s)
- Simone Rossi
- Department of Mathematics, UNC Chapel Hill, Chapel Hill, NC, United States
| | - Laryssa Abdala
- Department of Mathematics, UNC Chapel Hill, Chapel Hill, NC, United States
| | - Andrew Woodward
- Advanced Medical Imaging Lab, UNC Chapel Hill, Chapel Hill, NC, United States
| | - John P. Vavalle
- Department of Medicine, UNC Chapel Hill, Chapel Hill, NC, United States
| | - Craig S. Henriquez
- Department of Biomedical Engineering, Duke University, Durham, NC, United States
| | - Boyce E. Griffith
- Department of Mathematics, UNC Chapel Hill, Chapel Hill, NC, United States
- Department of Biomedical Engineering, UNC Chapel Hill, Chapel Hill, NC, United States
- McAllister Heart Institute, UNC Chapel Hill, Chapel Hill, NC, United States
| |
Collapse
|
9
|
Azzolin L, Eichenlaub M, Nagel C, Nairn D, Sanchez J, Unger L, Dössel O, Jadidi A, Loewe A. Personalized ablation vs. conventional ablation strategies to terminate atrial fibrillation and prevent recurrence. Europace 2022; 25:211-222. [PMID: 35943361 PMCID: PMC9907752 DOI: 10.1093/europace/euac116] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2022] [Accepted: 06/17/2022] [Indexed: 11/14/2022] Open
Abstract
AIMS The long-term success rate of ablation therapy is still sub-optimal in patients with persistent atrial fibrillation (AF), mostly due to arrhythmia recurrence originating from arrhythmogenic sites outside the pulmonary veins. Computational modelling provides a framework to integrate and augment clinical data, potentially enabling the patient-specific identification of AF mechanisms and of the optimal ablation sites. We developed a technology to tailor ablations in anatomical and functional digital atrial twins of patients with persistent AF aiming to identify the most successful ablation strategy. METHODS AND RESULTS Twenty-nine patient-specific computational models integrating clinical information from tomographic imaging and electro-anatomical activation time and voltage maps were generated. Areas sustaining AF were identified by a personalized induction protocol at multiple locations. State-of-the-art anatomical and substrate ablation strategies were compared with our proposed Personalized Ablation Lines (PersonAL) plan, which consists of iteratively targeting emergent high dominant frequency (HDF) regions, to identify the optimal ablation strategy. Localized ablations were connected to the closest non-conductive barrier to prevent recurrence of AF or atrial tachycardia. The first application of the HDF strategy had a success of >98% and isolated only 5-6% of the left atrial myocardium. In contrast, conventional ablation strategies targeting anatomical or structural substrate resulted in isolation of up to 20% of left atrial myocardium. After a second iteration of the HDF strategy, no further arrhythmia episode could be induced in any of the patient-specific models. CONCLUSION The novel PersonAL in silico technology allows to unveil all AF-perpetuating areas and personalize ablation by leveraging atrial digital twins.
Collapse
Affiliation(s)
- Luca Azzolin
- Corresponding author. Tel: +393381319986, E-mail address:
| | | | - Claudia Nagel
- Institute of Biomedical Engineering at Karlsruhe Institute of Technology, Building 30.33, Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | - Deborah Nairn
- Institute of Biomedical Engineering at Karlsruhe Institute of Technology, Building 30.33, Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | - Jorge Sanchez
- Institute of Biomedical Engineering at Karlsruhe Institute of Technology, Building 30.33, Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | - Laura Unger
- Institute of Biomedical Engineering at Karlsruhe Institute of Technology, Building 30.33, Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | - Olaf Dössel
- Institute of Biomedical Engineering at Karlsruhe Institute of Technology, Building 30.33, Fritz-Haber-Weg 1, 76131 Karlsruhe, Germany
| | | | | |
Collapse
|
10
|
Saliani A, Biswas S, Jacquemet V. Simulation of atrial fibrillation in a non-ohmic propagation model with dynamic gap junctions. CHAOS (WOODBURY, N.Y.) 2022; 32:043113. [PMID: 35489863 DOI: 10.1063/5.0082763] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Accepted: 03/24/2022] [Indexed: 06/14/2023]
Abstract
Gap junctions exhibit nonlinear electrical properties that have been hypothesized to be relevant to arrhythmogenicity in a structurally remodeled tissue. Large-scale implementation of gap junction dynamics in 3D propagation models remains challenging. We aim to quantify the impact of nonlinear diffusion during episodes of arrhythmias simulated in a left atrial model. Homogenization of conduction properties in the presence of nonlinear gap junctions was performed by generalizing a previously developed mathematical framework. A monodomain model was solved in which conductivities were time-varying and depended on transjunctional potentials. Gap junction conductances were derived from a simplified Vogel-Weingart model with first-order gating and adjustable time constant. A bilayer interconnected cable model of the left atrium with 100 μm resolution was used. The diffusion matrix was recomputed at each time step according to the state of the gap junctions. Sinus rhythm and atrial fibrillation episodes were simulated in remodeled tissue substrates. Slow conduction was induced by reduced coupling and by diffuse or stringy fibrosis. Simulations starting from the same initial conditions were repeated with linear and nonlinear gap junctions. The discrepancy in activation times between the linear and nonlinear diffusion models was quantified. The results largely validated the linear approximation for conduction velocities >20 cm/s. In very slow conduction substrates, the discrepancy accumulated over time during atrial fibrillation, eventually leading to qualitative differences in propagation patterns, while keeping the descriptive statistics, such as cycle lengths, unchanged. The discrepancy growth rate was increased by impaired conduction, fibrosis, conduction heterogeneity, lateral uncoupling, fast gap junction time constant, and steeper action potential duration restitution.
Collapse
Affiliation(s)
- Ariane Saliani
- Institute of Biomedical Engineering, Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, C.P. 6128, succ. Centre-ville, Montreal, Quebec H3C 3J7, Canada
| | - Subhamoy Biswas
- Institute of Biomedical Engineering, Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, C.P. 6128, succ. Centre-ville, Montreal, Quebec H3C 3J7, Canada
| | - Vincent Jacquemet
- Institute of Biomedical Engineering, Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, C.P. 6128, succ. Centre-ville, Montreal, Quebec H3C 3J7, Canada
| |
Collapse
|
11
|
Feng Y, Roney CH, Bayer JD, Niederer SA, Hocini M, Vigmond EJ. Detection of focal source and arrhythmogenic substrate from body surface potentials to guide atrial fibrillation ablation. PLoS Comput Biol 2022; 18:e1009893. [PMID: 35312675 PMCID: PMC8970486 DOI: 10.1371/journal.pcbi.1009893] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Revised: 03/31/2022] [Accepted: 02/02/2022] [Indexed: 11/18/2022] Open
Abstract
Focal sources (FS) are believed to be important triggers and a perpetuation mechanism for paroxysmal atrial fibrillation (AF). Detecting FS and determining AF sustainability in atrial tissue can help guide ablation targeting. We hypothesized that sustained rotors during FS-driven episodes indicate an arrhythmogenic substrate for sustained AF, and that non-invasive electrical recordings, like electrocardiograms (ECGs) or body surface potential maps (BSPMs), could be used to detect FS and AF sustainability. Computer simulations were performed on five bi-atrial geometries. FS were induced by pacing at cycle lengths of 120-270 ms from 32 atrial sites and four pulmonary veins. Self-sustained reentrant activities were also initiated around the same 32 atrial sites with inexcitable cores of radii of 0, 0.5 and 1 cm. FS fired for two seconds and then AF inducibility was tested by whether activation was sustained for another second. ECGs and BSPMs were simulated. Equivalent atrial sources were extracted using second-order blind source separation, and their cycle length, periodicity and contribution, were used as features for random forest classifiers. Longer rotor duration during FS-driven episodes indicates higher AF inducibility (area under ROC curve = 0.83). Our method had accuracy of 90.6±1.0% and 90.6±0.6% in detecting FS presence, and 93.1±0.6% and 94.2±1.2% in identifying AF sustainability, and 80.0±6.6% and 61.0±5.2% in determining the atrium of the focal site, from BSPMs and ECGs of five atria. The detection of FS presence and AF sustainability were insensitive to vest placement (±9.6%). On pre-operative BSPMs of 52 paroxysmal AF patients, patients classified with initiator-type FS on a single atrium resulted in improved two-to-three-year AF-free likelihoods (p-value < 0.01, logrank tests). Detection of FS and arrhythmogenic substrate can be performed from ECGs and BSPMs, enabling non-invasive mapping towards mechanism-targeted AF treatment, and malignant ectopic beat detection with likely AF progression.
Collapse
Affiliation(s)
- Yingjing Feng
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux, France
- Univ. Bordeaux, IMB, UMR 5251, Talence, France
| | - Caroline H. Roney
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Jason D. Bayer
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux, France
- Univ. Bordeaux, IMB, UMR 5251, Talence, France
| | - Steven A. Niederer
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Mélèze Hocini
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux, France
- Bordeaux University Hospital (CHU), Electrophysiology and Ablation Unit, Pessac, France
| | - Edward J. Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, Pessac-Bordeaux, France
- Univ. Bordeaux, IMB, UMR 5251, Talence, France
| |
Collapse
|
12
|
Roney CH, Sim I, Yu J, Beach M, Mehta A, Alonso Solis-Lemus J, Kotadia I, Whitaker J, Corrado C, Razeghi O, Vigmond E, Narayan SM, O’Neill M, Williams SE, Niederer SA. Predicting Atrial Fibrillation Recurrence by Combining Population Data and Virtual Cohorts of Patient-Specific Left Atrial Models. Circ Arrhythm Electrophysiol 2022; 15:e010253. [PMID: 35089057 PMCID: PMC8845531 DOI: 10.1161/circep.121.010253] [Citation(s) in RCA: 29] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Accepted: 01/03/2022] [Indexed: 01/14/2023]
Abstract
BACKGROUND Current ablation therapy for atrial fibrillation is suboptimal, and long-term response is challenging to predict. Clinical trials identify bedside properties that provide only modest prediction of long-term response in populations, while patient-specific models in small cohorts primarily explain acute response to ablation. We aimed to predict long-term atrial fibrillation recurrence after ablation in large cohorts, by using machine learning to complement biophysical simulations by encoding more interindividual variability. METHODS Patient-specific models were constructed for 100 atrial fibrillation patients (43 paroxysmal, 41 persistent, and 16 long-standing persistent), undergoing first ablation. Patients were followed for 1 year using ambulatory ECG monitoring. Each patient-specific biophysical model combined differing fibrosis patterns, fiber orientation maps, electrical properties, and ablation patterns to capture uncertainty in atrial properties and to test the ability of the tissue to sustain fibrillation. These simulation stress tests of different model variants were postprocessed to calculate atrial fibrillation simulation metrics. Machine learning classifiers were trained to predict atrial fibrillation recurrence using features from the patient history, imaging, and atrial fibrillation simulation metrics. RESULTS We performed 1100 atrial fibrillation ablation simulations across 100 patient-specific models. Models based on simulation stress tests alone showed a maximum accuracy of 0.63 for predicting long-term fibrillation recurrence. Classifiers trained to history, imaging, and simulation stress tests (average 10-fold cross-validation area under the curve, 0.85±0.09; recall, 0.80±0.13; precision, 0.74±0.13) outperformed those trained to history and imaging (area under the curve, 0.66±0.17) or history alone (area under the curve, 0.61±0.14). CONCLUSION A novel computational pipeline accurately predicted long-term atrial fibrillation recurrence in individual patients by combining outcome data with patient-specific acute simulation response. This technique could help to personalize selection for atrial fibrillation ablation.
Collapse
Affiliation(s)
- Caroline H. Roney
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
- School of Engineering and Materials Science, Queen Mary University of London, United Kingdom (C.H.R.)
| | - Iain Sim
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
| | - Jin Yu
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
| | - Marianne Beach
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
| | - Arihant Mehta
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
| | - Jose Alonso Solis-Lemus
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
| | - Irum Kotadia
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
| | - John Whitaker
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
- The Department of Internal Medicine, Cardiovascular Division, Brigham and Women’s Hospital, Boston, MA (J.W.)
| | - Cesare Corrado
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
| | - Orod Razeghi
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
| | - Edward Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, France (E.V.)
- Univ. Bordeaux, IMB, UMR 5251, F-33400 Talence, France (E.V.)
| | - Sanjiv M. Narayan
- Department of Medicine and Cardiovascular Institute, Stanford University, Palo Alto, CA (S.M.N.)
| | - Mark O’Neill
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
| | - Steven E. Williams
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
- Centre for Cardiovascular Science, College of Medicine and Veterinary Medicine, University of Edinburgh (S.E.W.)
| | - Steven A. Niederer
- School of Biomedical Engineering and Imaging Sciences, King’s College London, United Kingdom (C.H.R., I.S., J.Y., M.B., A.M., J.A.S.-L., I.K., J.W., C.C., O.R., M.O., S.E.W., S.A.N.)
| |
Collapse
|
13
|
Roney CH, Child N, Porter B, Sim I, Whitaker J, Clayton RH, Laughner JI, Shuros A, Neuzil P, Williams SE, Razavi RS, O'Neill M, Rinaldi CA, Taggart P, Wright M, Gill JS, Niederer SA. Time-Averaged Wavefront Analysis Demonstrates Preferential Pathways of Atrial Fibrillation, Predicting Pulmonary Vein Isolation Acute Response. Front Physiol 2021; 12:707189. [PMID: 34646149 PMCID: PMC8503618 DOI: 10.3389/fphys.2021.707189] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2021] [Accepted: 08/24/2021] [Indexed: 11/13/2022] Open
Abstract
Electrical activation during atrial fibrillation (AF) appears chaotic and disorganised, which impedes characterisation of the underlying substrate and treatment planning. While globally chaotic, there may be local preferential activation pathways that represent potential ablation targets. This study aimed to identify preferential activation pathways during AF and predict the acute ablation response when these are targeted by pulmonary vein isolation (PVI). In patients with persistent AF (n = 14), simultaneous biatrial contact mapping with basket catheters was performed pre-ablation and following each ablation strategy (PVI, roof, and mitral lines). Unipolar wavefront activation directions were averaged over 10 s to identify preferential activation pathways. Clinical cases were classified as responders or non-responders to PVI during the procedure. Clinical data were augmented with a virtual cohort of 100 models. In AF pre-ablation, pathways originated from the pulmonary vein (PV) antra in PVI responders (7/7) but not in PVI non-responders (6/6). We proposed a novel index that measured activation waves from the PV antra into the atrial body. This index was significantly higher in PVI responders than non-responders (clinical: 16.3 vs. 3.7%, p = 0.04; simulated: 21.1 vs. 14.1%, p = 0.02). Overall, this novel technique and proof of concept study demonstrated that preferential activation pathways exist during AF. Targeting patient-specific activation pathways that flowed from the PV antra to the left atrial body using PVI resulted in AF termination during the procedure. These PV activation flow pathways may correspond to the presence of drivers in the PV regions.
Collapse
Affiliation(s)
- Caroline H. Roney
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Nicholas Child
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Bradley Porter
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Iain Sim
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - John Whitaker
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Richard H. Clayton
- INSIGNEO Institute for In Silico Medicine and Department of Computer Science, University of Sheffield, Sheffield, United Kingdom
| | | | - Allan Shuros
- Boston Scientific Corp, St. Paul, MN, United States
| | - Petr Neuzil
- Department of Cardiology, Na Holmolce Hospital, Prague, Czechia
| | - Steven E. Williams
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Reza S. Razavi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Mark O'Neill
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | | | - Peter Taggart
- Institute of Cardiovascular Science, University College London, London, United Kingdom
| | - Matt Wright
- Department of Cardiology, Guy's and St Thomas' Hospital, London, United Kingdom
| | - Jaswinder S. Gill
- Department of Cardiology, Guy's and St Thomas' Hospital, London, United Kingdom
| | - Steven A. Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| |
Collapse
|
14
|
Saliani A, Irakoze É, Jacquemet V. Simulation of diffuse and stringy fibrosis in a bilayer interconnected cable model of the left atrium. Europace 2021; 23:i169-i177. [PMID: 33751082 DOI: 10.1093/europace/euab001] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2020] [Accepted: 01/04/2021] [Indexed: 11/13/2022] Open
Abstract
AIMS The aim of this study is to design a computer model of the left atrium for investigating fibre-orientation-dependent microstructure such as stringy fibrosis. METHODS AND RESULTS We developed an approach for automatic construction of bilayer interconnected cable models from left atrial geometry and epi- and endocardial fibre orientation. The model consisted of two layers (epi- and endocardium) of longitudinal and transverse cables intertwined-like fabric threads, with a spatial discretization of 100 µm. Model validation was performed by comparison with cubic volumetric models in normal conditions. Then, diffuse (n = 2904), stringy (n = 3600), and mixed fibrosis patterns (n = 6840) were randomly generated by uncoupling longitudinal and transverse connections in the interconnected cable model. Fibrosis density was varied from 0% to 40% and mean stringy obstacle length from 0.1 to 2 mm. Total activation time, apparent anisotropy ratio, and local activation time jitter were computed during normal rhythm in each pattern. Non-linear regression formulas were identified for expressing measured propagation parameters as a function of fibrosis density and obstacle length (stringy and mixed patterns). Longer obstacles (even below tissue space constant) were independently associated with prolonged activation times, increased anisotropy, and local fluctuations in activation times. This effect was increased by endo-epicardial dissociation and mitigated when fibrosis was limited to the epicardium. CONCLUSION Interconnected cable models enable the study of microstructure in organ-size models despite limitations in the description of transmural structures.
Collapse
Affiliation(s)
- Ariane Saliani
- Department of Pharmacology and Physiology, Institute of Biomedical Engineering, Université de Montréal, Montréal, QC H3T 1J4, Canada.,Research Center, Hôpital du Sacré-Cœur de Montréal, 5400 boul. Gouin Ouest, Montréal, QC H4J 1C5, Canada
| | - Éric Irakoze
- Department of Pharmacology and Physiology, Institute of Biomedical Engineering, Université de Montréal, Montréal, QC H3T 1J4, Canada.,Research Center, Hôpital du Sacré-Cœur de Montréal, 5400 boul. Gouin Ouest, Montréal, QC H4J 1C5, Canada
| | - Vincent Jacquemet
- Department of Pharmacology and Physiology, Institute of Biomedical Engineering, Université de Montréal, Montréal, QC H3T 1J4, Canada.,Research Center, Hôpital du Sacré-Cœur de Montréal, 5400 boul. Gouin Ouest, Montréal, QC H4J 1C5, Canada
| |
Collapse
|
15
|
Hwang I, Jin Z, Park JW, Kwon OS, Lim B, Hong M, Kim M, Yu HT, Kim TH, Uhm JS, Joung B, Lee MH, Pak HN. Computational Modeling for Antiarrhythmic Drugs for Atrial Fibrillation According to Genotype. Front Physiol 2021; 12:650449. [PMID: 34054570 PMCID: PMC8155488 DOI: 10.3389/fphys.2021.650449] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2021] [Accepted: 03/22/2021] [Indexed: 01/11/2023] Open
Abstract
Background: The efficacy of antiarrhythmic drugs (AAD) can vary in patients with atrial fibrillation (AF), and the PITX2 gene affects the responsiveness of AADs. We explored the virtual AAD (V-AAD) responses between wild-type and PITX2 +/--deficient AF conditions by realistic in silico AF modeling. Methods: We tested the V-AADs in AF modeling integrated with patients' 3D-computed tomography and 3D-electroanatomical mapping, acquired in 25 patients (68% male, 59.8 ± 9.8 years old, 32.0% paroxysmal type). The ion currents for the PITX2 +/- deficiency and each AAD (amiodarone, sotalol, dronedarone, flecainide, and propafenone) were defined based on previous publications. Results: We compared the wild-type and PITX2 +/- deficiency in terms of the action potential duration (APD90), conduction velocity (CV), maximal slope of restitution (Smax), and wave-dynamic parameters, such as the dominant frequency (DF), phase singularities (PS), and AF termination rates according to the V-AADs. The PITX2 +/--deficient model exhibited a shorter APD90 (p < 0.001), a lower Smax (p < 0.001), mean DF (p = 0.012), PS number (p < 0.001), and a longer AF cycle length (AFCL, p = 0.011). Five V-AADs changed the electrophysiology in a dose-dependent manner. AAD-induced AFCL lengthening (p < 0.001) and reductions in the CV (p = 0.033), peak DF (p < 0.001), and PS number (p < 0.001) were more significant in PITX2 +/--deficient than wild-type AF. PITX2 +/--deficient AF was easier to terminate with class IC AADs than the wild-type AF (p = 0.018). Conclusions: The computational modeling-guided AAD test was feasible for evaluating the efficacy of multiple AADs in patients with AF. AF wave-dynamic and electrophysiological characteristics are different among the PITX2-deficient and the wild-type genotype models.
Collapse
|
16
|
Bifulco SF, Scott GD, Sarairah S, Birjandian Z, Roney CH, Niederer SA, Mahnkopf C, Kuhnlein P, Mitlacher M, Tirschwell D, Longstreth WT, Akoum N, Boyle PM. Computational modeling identifies embolic stroke of undetermined source patients with potential arrhythmic substrate. eLife 2021; 10:e64213. [PMID: 33942719 PMCID: PMC8143793 DOI: 10.7554/elife.64213] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2020] [Accepted: 04/16/2021] [Indexed: 12/25/2022] Open
Abstract
Cardiac magnetic resonance imaging (MRI) has revealed fibrosis in embolic stroke of undetermined source (ESUS) patients comparable to levels seen in atrial fibrillation (AFib). We used computational modeling to understand the absence of arrhythmia in ESUS despite the presence of putatively pro-arrhythmic fibrosis. MRI-based atrial models were reconstructed for 45 ESUS and 45 AFib patients. The fibrotic substrate's arrhythmogenic capacity in each patient was assessed computationally. Reentrant drivers were induced in 24/45 (53%) ESUS and 22/45 (49%) AFib models. Inducible models had more fibrosis (16.7 ± 5.45%) than non-inducible models (11.07 ± 3.61%; p<0.0001); however, inducible subsets of ESUS and AFib models had similar fibrosis levels (p=0.90), meaning that the intrinsic pro-arrhythmic substrate properties of fibrosis in ESUS and AFib are indistinguishable. This suggests that some ESUS patients have latent pre-clinical fibrotic substrate that could be a future source of arrhythmogenicity. Thus, our work prompts the hypothesis that ESUS patients with fibrotic atria are spared from AFib due to an absence of arrhythmia triggers.
Collapse
Affiliation(s)
- Savannah F Bifulco
- Department of Bioengineering, University of WashingtonSeattleUnited States
| | - Griffin D Scott
- Department of Bioengineering, University of WashingtonSeattleUnited States
| | - Sakher Sarairah
- Division of Cardiology, University of WashingtonSeattleUnited States
| | - Zeinab Birjandian
- Division of Cardiology, University of WashingtonSeattleUnited States
- Department of Neurology, University of WashingtonSeattleUnited States
| | - Caroline H Roney
- School of Biomedical Engineering and Imaging Sciences, King’s College LondonLondonUnited Kingdom
| | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King’s College LondonLondonUnited Kingdom
| | | | | | | | - David Tirschwell
- Department of Neurology, University of WashingtonSeattleUnited States
| | - WT Longstreth
- Department of Neurology, University of WashingtonSeattleUnited States
- Department of Epidemiology, University of WashingtonSeattleUnited States
| | - Nazem Akoum
- Division of Cardiology, University of WashingtonSeattleUnited States
| | - Patrick M Boyle
- Department of Bioengineering, University of WashingtonSeattleUnited States
- Center for Cardiovascular Biology, University of WashingtonSeattleUnited States
- Institute for Stem Cell and Regenerative Medicine, University of WashingtonSeattleUnited States
| |
Collapse
|
17
|
Roney CH, Beach ML, Mehta AM, Sim I, Corrado C, Bendikas R, Solis-Lemus JA, Razeghi O, Whitaker J, O’Neill L, Plank G, Vigmond E, Williams SE, O’Neill MD, Niederer SA. In silico Comparison of Left Atrial Ablation Techniques That Target the Anatomical, Structural, and Electrical Substrates of Atrial Fibrillation. Front Physiol 2020; 11:1145. [PMID: 33041850 PMCID: PMC7526475 DOI: 10.3389/fphys.2020.572874] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 08/18/2020] [Indexed: 12/17/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Caroline H. Roney
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Marianne L. Beach
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Arihant M. Mehta
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Iain Sim
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Cesare Corrado
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Rokas Bendikas
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Jose A. Solis-Lemus
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Orod Razeghi
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - John Whitaker
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Louisa O’Neill
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Gernot Plank
- Department of Biophysics, Medical University of Graz, Graz, Austria
| | - Edward Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Bordeaux, France
| | - Steven E. Williams
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Mark D. O’Neill
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Steven A. Niederer
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| |
Collapse
|
18
|
Strocchi M, Augustin CM, Gsell MAF, Karabelas E, Neic A, Gillette K, Razeghi O, Prassl AJ, Vigmond EJ, Behar JM, Gould J, Sidhu B, Rinaldi CA, Bishop MJ, Plank G, Niederer SA. A publicly available virtual cohort of four-chamber heart meshes for cardiac electro-mechanics simulations. PLoS One 2020; 15:e0235145. [PMID: 32589679 PMCID: PMC7319311 DOI: 10.1371/journal.pone.0235145] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2020] [Accepted: 06/09/2020] [Indexed: 12/12/2022] Open
Abstract
Computational models of the heart are increasingly being used in the development of devices, patient diagnosis and therapy guidance. While software techniques have been developed for simulating single hearts, there remain significant challenges in simulating cohorts of virtual hearts from multiple patients. To facilitate the development of new simulation and model analysis techniques by groups without direct access to medical data, image analysis techniques and meshing tools, we have created the first publicly available virtual cohort of twenty-four four-chamber hearts. Our cohort was built from heart failure patients, age 67±14 years. We segmented four-chamber heart geometries from end-diastolic (ED) CT images and generated linear tetrahedral meshes with an average edge length of 1.1±0.2mm. Ventricular fibres were added in the ventricles with a rule-based method with an orientation of -60° and 80° at the epicardium and endocardium, respectively. We additionally refined the meshes to an average edge length of 0.39±0.10mm to show that all given meshes can be resampled to achieve an arbitrary desired resolution. We ran simulations for ventricular electrical activation and free mechanical contraction on all 1.1mm-resolution meshes to ensure that our meshes are suitable for electro-mechanical simulations. Simulations for electrical activation resulted in a total activation time of 149±16ms. Free mechanical contractions gave an average left ventricular (LV) and right ventricular (RV) ejection fraction (EF) of 35±1% and 30±2%, respectively, and a LV and RV stroke volume (SV) of 95±28mL and 65±11mL, respectively. By making the cohort publicly available, we hope to facilitate large cohort computational studies and to promote the development of cardiac computational electro-mechanics for clinical applications.
Collapse
Affiliation(s)
- Marina Strocchi
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, City of London, United Kingdom
| | | | | | - Elias Karabelas
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, City of London, United Kingdom
| | | | - Karli Gillette
- Institute of Biophysics, Medical University of Graz, Graz, Steiermark, Austria
| | - Orod Razeghi
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, City of London, United Kingdom
| | - Anton J. Prassl
- Institute of Biophysics, Medical University of Graz, Graz, Steiermark, Austria
| | - Edward J. Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, fondation Bordeaux Université, F-33600 Pessac- Bordeaux, France
- University of Bordeaux, IMB, UMR 5251, F-33400 Talence, France
| | - Jonathan M. Behar
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, City of London, United Kingdom
- Guy’s and St Thomas’ NHS Foundation Trust, London, City of London, United Kingdom
| | - Justin Gould
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, City of London, United Kingdom
- Guy’s and St Thomas’ NHS Foundation Trust, London, City of London, United Kingdom
| | - Baldeep Sidhu
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, City of London, United Kingdom
- Guy’s and St Thomas’ NHS Foundation Trust, London, City of London, United Kingdom
| | - Christopher A. Rinaldi
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, City of London, United Kingdom
- Guy’s and St Thomas’ NHS Foundation Trust, London, City of London, United Kingdom
| | - Martin J. Bishop
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, City of London, United Kingdom
| | - Gernot Plank
- Institute of Biophysics, Medical University of Graz, Graz, Steiermark, Austria
| | - Steven A. Niederer
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, City of London, United Kingdom
| |
Collapse
|
19
|
Abstract
Atrial anisotropy affects electrical propagation patterns, anchor locations of atrial reentrant drivers, and atrial mechanics. However, patient-specific atrial fibre fields and anisotropy measurements are not currently available, and consequently assigning fibre fields to atrial models is challenging. We aimed to construct an atrial fibre atlas from a high-resolution DTMRI dataset that optimally reproduces electrophysiology simulation predictions corresponding to patient-specific fibre fields, and to develop a methodology for automatically assigning fibres to patient-specific anatomies. We extended an atrial coordinate system to map the pulmonary veins, vena cava and appendages to standardised positions in the coordinate system corresponding to the average location across the anatomies. We then expressed each fibre field in this atrial coordinate system and calculated an average fibre field. To assess the effects of fibre field on patient-specific modelling predictions, we calculated paced activation time maps and electrical driver locations during AF. In total, 756 activation time maps were calculated (7 anatomies with 9 fibre maps and 2 pacing locations, for the endocardial, epicardial and bilayer surface models of the LA and RA). Patient-specific fibre fields had a relatively small effect on average paced activation maps (range of mean local activation time difference for LA fields: 2.67-3.60 ms, and for RA fields: 2.29-3.44 ms), but had a larger effect on maximum LAT differences (range for LA 12.7-16.6%; range for RA 11.9-15.0%). A total of 126 phase singularity density maps were calculated (7 anatomies with 9 fibre maps for the LA and RA bilayer models). The fibre field corresponding to anatomy 1 had the highest median PS density map correlation coefficient for LA bilayer simulations (0.44 compared to the other correlations, ranging from 0.14 to 0.39), while the average fibre field had the highest correlation for the RA bilayer simulations (0.61 compared to the other correlations, ranging from 0.37 to 0.56). For sinus rhythm simulations, average activation time is robust to fibre field direction; however, maximum differences can still be significant. Patient specific fibres are more important for arrhythmia simulations, particularly in the left atrium. We propose using the fibre field corresponding to DTMRI dataset 1 for LA simulations, and the average fibre field for RA simulations as these optimally predicted arrhythmia properties.
Collapse
|
20
|
Gharaviri A, Bidar E, Potse M, Zeemering S, Verheule S, Pezzuto S, Krause R, Maessen JG, Auricchio A, Schotten U. Epicardial Fibrosis Explains Increased Endo-Epicardial Dissociation and Epicardial Breakthroughs in Human Atrial Fibrillation. Front Physiol 2020; 11:68. [PMID: 32153419 PMCID: PMC7047215 DOI: 10.3389/fphys.2020.00068] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2019] [Accepted: 01/21/2020] [Indexed: 01/22/2023] Open
Abstract
Background Atrial fibrillation (AF) is accompanied by progressive epicardial fibrosis, dissociation of electrical activity between the epicardial layer and the endocardial bundle network, and transmural conduction (breakthroughs). However, causal relationships between these phenomena have not been demonstrated yet. Our goal was to test the hypothesis that epicardial fibrosis suffices to increase endo–epicardial dissociation (EED) and breakthroughs (BT) during AF. Methods We simulated the effect of fibrosis in the epicardial layer on EED and BT in a detailed, high-resolution, three-dimensional model of the human atria with realistic electrophysiology. The model results were compared with simultaneous endo–epicardial mapping in human atria. The model geometry, specifically built for this study, was based on MR images and histo-anatomical studies. Clinical data were obtained in four patients with longstanding persistent AF (persAF) and three patients without a history of AF. Results The AF cycle length (AFCL), conduction velocity (CV), and EED were comparable in the mapping studies and the simulations. EED increased from 24.1 ± 3.4 to 56.58 ± 6.2% (p < 0.05), and number of BTs per cycle from 0.89 ± 0.55 to 6.74 ± 2.11% (p < 0.05), in different degrees of fibrosis in the epicardial layer. In both mapping data and simulations, EED correlated with prevalence of BTs. Fibrosis also increased the number of fibrillation waves per cycle in the model. Conclusion A realistic 3D computer model of AF in which epicardial fibrosis was increased, in the absence of other pathological changes, showed increases in EED and epicardial BT comparable to those in longstanding persAF. Thus, epicardial fibrosis can explain both phenomena.
Collapse
Affiliation(s)
- Ali Gharaviri
- Department of Physiology, Maastricht University, Maastricht, Netherlands.,Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera Italiana, Lugano, Switzerland
| | - Elham Bidar
- Maastricht University Medical Centre, Maastricht, Netherlands
| | - Mark Potse
- Inria Bordeaux - Sud-Ouest Research Centre, Talence, France.,IMB, UMR 5251, Université de Bordeaux, Talence, France.,IHU Liryc, Electrophysiology and Heart Modeling Institute, Foundation Bordeaux Université, Bordeaux, France
| | - Stef Zeemering
- Department of Physiology, Maastricht University, Maastricht, Netherlands
| | - Sander Verheule
- Department of Physiology, Maastricht University, Maastricht, Netherlands
| | - Simone Pezzuto
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera Italiana, Lugano, Switzerland
| | - Rolf Krause
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera Italiana, Lugano, Switzerland
| | - Jos G Maessen
- Maastricht University Medical Centre, Maastricht, Netherlands
| | - Angelo Auricchio
- Center for Computational Medicine in Cardiology, Institute of Computational Science, Università della Svizzera Italiana, Lugano, Switzerland.,Fondazione Cardiocentro Ticino, Lugano, Switzerland
| | - Ulrich Schotten
- Department of Physiology, Maastricht University, Maastricht, Netherlands
| |
Collapse
|
21
|
Lim B, Kim J, Hwang M, Song JS, Lee JK, Yu HT, Kim TH, Uhm JS, Joung B, Lee MH, Pak HN. In situ procedure for high-efficiency computational modeling of atrial fibrillation reflecting personal anatomy, fiber orientation, fibrosis, and electrophysiology. Sci Rep 2020; 10:2417. [PMID: 32051487 PMCID: PMC7016008 DOI: 10.1038/s41598-020-59372-x] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2019] [Accepted: 11/06/2019] [Indexed: 12/22/2022] Open
Abstract
We previously reported the feasibility and efficacy of a simulation-guided clinical catheter ablation of atrial fibrillation (AF) in an in-silico AF model. We developed a highly efficient realistic AF model reflecting the patient endocardial voltage and local conduction and tested its clinical feasibility. We acquired > 500 endocardial bipolar electrograms during right atrial pacing at the beginning of the AF ablation procedures. Based on the clinical bipolar electrograms, we generated simulated voltage maps by applying fibrosis and local activation maps adjusted for the fiber orientation. The software's accuracy (CUVIA2.5) was retrospectively tested in 17 patients and feasibility prospectively in 10 during clinical AF ablation. Results: We found excellent correlations between the clinical and simulated voltage maps (R = 0.933, p < 0.001) and clinical and virtual local conduction (R = 0.958, p < 0.001). The proportion of virtual local fibrosis was 15.4, 22.2, and 36.9% in the paroxysmal AF, persistent AF, and post-pulmonary vein isolation (PVI) states, respectively. The reconstructed virtual bipolar electrogram exhibited a relatively good similarities of morphology to the local clinical bipolar electrogram (R = 0.60 ± 0.08, p < 0.001). Feasibility testing revealed an in situ procedural computing time from the clinical data acquisition to wave-dynamics analyses of 48.2 ± 4.9 min. All virtual analyses were successfully achieved during clinical PVI procedures. We developed a highly efficient, realistic, in situ procedural simulation model reflective of individual anatomy, fiber orientation, fibrosis, and electrophysiology that can be applied during AF ablation.
Collapse
Affiliation(s)
- Byounghyun Lim
- Yonsei University Health System, Seoul, Republic of Korea
| | - Jaehyeok Kim
- Yonsei University Health System, Seoul, Republic of Korea
| | - Minki Hwang
- Yonsei University Health System, Seoul, Republic of Korea
| | - Jun-Seop Song
- Yonsei University Health System, Seoul, Republic of Korea
| | - Jung Ki Lee
- Yonsei University Health System, Seoul, Republic of Korea
| | - Hee-Tae Yu
- Yonsei University Health System, Seoul, Republic of Korea
| | - Tae-Hoon Kim
- Yonsei University Health System, Seoul, Republic of Korea
| | - Jae-Sun Uhm
- Yonsei University Health System, Seoul, Republic of Korea
| | - Boyoung Joung
- Yonsei University Health System, Seoul, Republic of Korea
| | - Moon-Hyung Lee
- Yonsei University Health System, Seoul, Republic of Korea
| | - Hui-Nam Pak
- Yonsei University Health System, Seoul, Republic of Korea.
| |
Collapse
|
22
|
Strocchi M, Gsell MAF, Augustin CM, Razeghi O, Roney CH, Prassl AJ, Vigmond EJ, Behar JM, Gould JS, Rinaldi CA, Bishop MJ, Plank G, Niederer SA. Simulating ventricular systolic motion in a four-chamber heart model with spatially varying robin boundary conditions to model the effect of the pericardium. J Biomech 2020; 101:109645. [PMID: 32014305 PMCID: PMC7677892 DOI: 10.1016/j.jbiomech.2020.109645] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2019] [Revised: 01/15/2020] [Accepted: 01/15/2020] [Indexed: 12/11/2022]
Abstract
The pericardium affects cardiac motion by limiting epicardial displacement normal to the surface. In computational studies, it is important for the model to replicate realistic motion, as this affects the physiological fidelity of the model. Previous computational studies showed that accounting for the effect of the pericardium allows for a more realistic motion simulation. In this study, we describe the mechanism through which the pericardium causes improved cardiac motion. We simulated electrical activation and contraction of the ventricles on a four-chamber heart in the presence and absence of the effect of the pericardium. We simulated the mechanical constraints imposed by the pericardium by applying normal Robin boundary conditions on the ventricular epicardium. We defined a regional scaling of normal springs stiffness based on image-derived motion from CT images. The presence of the pericardium reduced the error between simulated and image-derived end-systolic configurations from 12.8±4.1 mm to 5.7±2.5 mm. First, the pericardium prevents the ventricles from spherising during isovolumic contraction, reducing the outward motion of the free walls normal to the surface and the upwards motion of the apex. Second, by restricting the inward motion of the free and apical walls of the ventricles the pericardium increases atrioventricular plane displacement by four folds during ejection. Our results provide a mechanistic explanation of the importance of the pericardium in physiological simulations of electromechanical cardiac function.
Collapse
Affiliation(s)
- Marina Strocchi
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK.
| | | | | | - Orod Razeghi
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Caroline H Roney
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Anton J Prassl
- Department of Biophysics, Medical University of Graz, Graz, Austria
| | - Edward J Vigmond
- University of Bordeaux, Talence, France; LIRYC Electrophysiology and Heart Modeling Institute, Campus Xavier Arnozan, Pessac, France
| | - Jonathan M Behar
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK; Cardiology Department, Guys and St Thomas' NHS Foundation Trust, London, UK
| | - Justin S Gould
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK; Cardiology Department, Guys and St Thomas' NHS Foundation Trust, London, UK
| | - Christopher A Rinaldi
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK; Cardiology Department, Guys and St Thomas' NHS Foundation Trust, London, UK
| | - Martin J Bishop
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Gernot Plank
- Department of Biophysics, Medical University of Graz, Graz, Austria
| | - Steven A Niederer
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| |
Collapse
|
23
|
Bayer JD, Boukens BJ, Krul SPJ, Roney CH, Driessen AHG, Berger WR, van den Berg NWE, Verkerk AO, Vigmond EJ, Coronel R, de Groot JR. Acetylcholine Delays Atrial Activation to Facilitate Atrial Fibrillation. Front Physiol 2019; 10:1105. [PMID: 31551802 PMCID: PMC6737394 DOI: 10.3389/fphys.2019.01105] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2018] [Accepted: 08/09/2019] [Indexed: 11/13/2022] Open
Abstract
Background Acetylcholine (ACh) shortens action potential duration (APD) in human atria. APD shortening facilitates atrial fibrillation (AF) by reducing the wavelength for reentry. However, the influence of ACh on electrical conduction in human atria and its contribution to AF are unclear, particularly when combined with impaired conduction from interstitial fibrosis. Objective To investigate the effect of ACh on human atrial conduction and its role in AF with computational, experimental, and clinical approaches. Methods S1S2 pacing (S1 = 600 ms and S2 = variable cycle lengths) was applied to the following human AF computer models: a left atrial appendage (LAA) myocyte to quantify the effects of ACh on APD, maximum upstroke velocity (V max ), and resting membrane potential (RMP); a monolayer of LAA myocytes to quantify the effects of ACh on conduction; and 3) an intact left atrium (LA) to determine the effects of ACh on arrhythmogenicity. Heterogeneous ACh and interstitial fibrosis were applied to the monolayer and LA models. To corroborate the simulations, APD and RMP from isolated human atrial myocytes were recorded before and after 0.1 μM ACh. At the tissue level, LAAs from AF patients were optically mapped ex vivo using Di-4-ANEPPS. The difference in total activation time (AT) was determined between AT initially recorded with S1 pacing, and AT recorded during subsequent S1 pacing without (n = 6) or with (n = 7) 100 μM ACh. Results In LAA myocyte simulations, S1 pacing with 0.1 μM ACh shortened APD by 41 ms, hyperpolarized RMP by 7 mV, and increased V max by 27 mV/ms. In human atrial myocytes, 0.1 μM ACh shortened APD by 48 ms, hyperpolarized RMP by 3 mV, and increased V max by 6 mV/ms. In LAA monolayer simulations, S1 pacing with ACh hyperpolarized RMP to delay total AT by 32 ms without and 35 ms with fibrosis. This led to unidirectional conduction block and sustained reentry in fibrotic LA with heterogeneous ACh during S2 pacing. In AF patient LAAs, S1 pacing with ACh increased total AT from 39.3 ± 26 ms to 71.4 ± 31.2 ms (p = 0.036) compared to no change without ACh (56.7 ± 29.3 ms to 50.0 ± 21.9 ms, p = 0.140). Conclusion In fibrotic atria with heterogeneous parasympathetic activation, ACh facilitates AF by shortening APD and slowing conduction to promote unidirectional conduction block and reentry.
Collapse
Affiliation(s)
- Jason D Bayer
- Electrophysiology and Heart Modeling Institute (IHU-LIRYC), Bordeaux University Foundation, Bordeaux, France.,Institute of Mathematics of Bordeaux (U5251), University of Bordeaux, Bordeaux, France
| | - Bastiaan J Boukens
- Department of Medical Biology, Academic Medical Center, Amsterdam, Netherlands
| | - Sébastien P J Krul
- Department of Cardiology, Academic Medical Center, Amsterdam, Netherlands
| | - Caroline H Roney
- Division of Imaging Sciences and Bioengineering, King's College London, London, United Kingdom
| | | | - Wouter R Berger
- Department of Cardiology, Academic Medical Center, Amsterdam, Netherlands.,Department of Cardiology, Heart Center, OLVG, Amsterdam, Netherlands
| | | | - Arie O Verkerk
- Department of Medical Biology, Academic Medical Center, Amsterdam, Netherlands.,Department of Experimental Cardiology, Academic Medical Center, Amsterdam, Netherlands
| | - Edward J Vigmond
- Electrophysiology and Heart Modeling Institute (IHU-LIRYC), Bordeaux University Foundation, Bordeaux, France.,Institute of Mathematics of Bordeaux (U5251), University of Bordeaux, Bordeaux, France
| | - Ruben Coronel
- Electrophysiology and Heart Modeling Institute (IHU-LIRYC), Bordeaux University Foundation, Bordeaux, France.,Department of Experimental Cardiology, Academic Medical Center, Amsterdam, Netherlands
| | - Joris R de Groot
- Department of Cardiology, Academic Medical Center, Amsterdam, Netherlands
| |
Collapse
|
24
|
Roney CH, Williams SE, Cochet H, Mukherjee RK, O'Neill L, Sim I, Whitaker J, Razeghi O, Klein GJ, Vigmond EJ, O'Neill M, Niederer SA. Patient-specific simulations predict efficacy of ablation of interatrial connections for treatment of persistent atrial fibrillation. Europace 2019; 20:iii55-iii68. [PMID: 30476055 PMCID: PMC6251187 DOI: 10.1093/europace/euy232] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2018] [Accepted: 10/12/2018] [Indexed: 11/23/2022] Open
Abstract
Aims Treatments for persistent atrial fibrillation (AF) offer limited efficacy. One potential strategy aims to return the right atrium (RA) to sinus rhythm (SR) by ablating interatrial connections (IAC) to isolate the atria, but there is limited clinical data to evaluate this ablation approach. We aimed to use simulation to evaluate and predict patient-specific suitability for ablation of IAC to treat AF. Methods and results Persistent AF was simulated in 12 patient-specific geometries, incorporating electrophysiological heterogeneity and fibres, with IAC at Bachmann’s bundle, the coronary sinus, and fossa ovalis. Simulations were performed to test the effect of left atrial (LA)-to-RA frequency gradient and fibrotic remodelling on IAC ablation efficacy. During AF, we simulated ablation of one, two, or all three IAC, with or without pulmonary vein isolation and determined if this altered or terminated the arrhythmia. For models without structural remodelling, ablating all IAC terminated RA arrhythmia in 83% of cases. Models with the LA-to-RA frequency gradient removed had an increased success rate (100% success). Ablation of IACs is less effective in cases with fibrotic remodelling (interstitial fibrosis 50% success rate; combination remodelling 67%). Mean number of phase singularities in the RA was higher pre-ablation for IAC failure (success 0.6 ± 0.8 vs. failure 3.2 ± 2.5, P < 0.001). Conclusion This simulation study predicts that IAC ablation is effective in returning the RA to SR for many cases. Patient-specific modelling approaches have the potential to stratify patients prior to ablation by predicting if drivers are located in the LA or RA. We present a platform for predicting efficacy and informing patient selection for speculative treatments.
Collapse
Affiliation(s)
- Caroline H Roney
- School of Biomedical Engineering & Imaging Sciences, King's College London, St. Thomas' Hospital, Westminster Bridge Road, UK
| | - Steven E Williams
- School of Biomedical Engineering & Imaging Sciences, King's College London, St. Thomas' Hospital, Westminster Bridge Road, UK
| | - Hubert Cochet
- LIRYC Electrophysiology and Heart Modeling Institute, Bordeaux Fondation, Avenue du Haut-Lévèque, Pessac, France
| | - Rahul K Mukherjee
- School of Biomedical Engineering & Imaging Sciences, King's College London, St. Thomas' Hospital, Westminster Bridge Road, UK
| | - Louisa O'Neill
- School of Biomedical Engineering & Imaging Sciences, King's College London, St. Thomas' Hospital, Westminster Bridge Road, UK
| | - Iain Sim
- School of Biomedical Engineering & Imaging Sciences, King's College London, St. Thomas' Hospital, Westminster Bridge Road, UK
| | - John Whitaker
- School of Biomedical Engineering & Imaging Sciences, King's College London, St. Thomas' Hospital, Westminster Bridge Road, UK
| | - Orod Razeghi
- School of Biomedical Engineering & Imaging Sciences, King's College London, St. Thomas' Hospital, Westminster Bridge Road, UK
| | | | - Edward J Vigmond
- LIRYC Electrophysiology and Heart Modeling Institute, Bordeaux Fondation, Avenue du Haut-Lévèque, Pessac, France.,IMB, Univ. Bordeaux, Talence, France
| | - Mark O'Neill
- School of Biomedical Engineering & Imaging Sciences, King's College London, St. Thomas' Hospital, Westminster Bridge Road, UK
| | - Steven A Niederer
- School of Biomedical Engineering & Imaging Sciences, King's College London, St. Thomas' Hospital, Westminster Bridge Road, UK
| |
Collapse
|
25
|
Computationally guided personalized targeted ablation of persistent atrial fibrillation. Nat Biomed Eng 2019; 3:870-879. [PMID: 31427780 PMCID: PMC6842421 DOI: 10.1038/s41551-019-0437-9] [Citation(s) in RCA: 147] [Impact Index Per Article: 29.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2018] [Accepted: 07/03/2019] [Indexed: 12/12/2022]
Abstract
Atrial fibrillation (AF) — the most common arrhythmia — significantly increases the risk of stroke and heart failure. Although catheter ablation can restore normal heart rhythms, patients with persistent AF who develop atrial fibrosis often undergo multiple failed ablations and thus increased procedural risks. Here, we present personalized computational modelling for the reliable predetermination of ablation targets, which are then used to guide the ablation procedure in patients with persistent AF and atrial fibrosis. We first show that a computational model of the atria of patients identifies fibrotic tissue that if ablated will not sustain AF. We then integrated the target-ablation sites in a clinical-mapping system, and tested its feasibility in 10 patients with persistent AF. The computational prediction of ablation targets avoids lengthy electrical mapping and could improve the accuracy and efficacy of targeted AF ablation in patients whilst eliminating the need for repeat procedures.
Collapse
|
26
|
Aronis KN, Ali RL, Liang JA, Zhou S, Trayanova NA. Understanding AF Mechanisms Through Computational Modelling and Simulations. Arrhythm Electrophysiol Rev 2019; 8:210-219. [PMID: 31463059 PMCID: PMC6702471 DOI: 10.15420/aer.2019.28.2] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/08/2019] [Accepted: 06/17/2019] [Indexed: 12/21/2022] Open
Abstract
AF is a progressive disease of the atria, involving complex mechanisms related to its initiation, maintenance and progression. Computational modelling provides a framework for integration of experimental and clinical findings, and has emerged as an essential part of mechanistic research in AF. The authors summarise recent advancements in development of multi-scale AF models and focus on the mechanistic links between alternations in atrial structure and electrophysiology with AF. Key AF mechanisms that have been explored using atrial modelling are pulmonary vein ectopy; atrial fibrosis and fibrosis distribution; atrial wall thickness heterogeneity; atrial adipose tissue infiltration; development of repolarisation alternans; cardiac ion channel mutations; and atrial stretch with mechano-electrical feedback. They review modelling approaches that capture variability at the cohort level and provide cohort-specific mechanistic insights. The authors conclude with a summary of future perspectives, as envisioned for the contributions of atrial modelling in the mechanistic understanding of AF.
Collapse
Affiliation(s)
- Konstantinos N Aronis
- Department of Biomedical Engineering and the Institute for Computational Medicine, Johns Hopkins UniversityBaltimore, MD, US
- Division of Cardiology, Johns Hopkins HospitalBaltimore, MD, US
| | - Rheeda L Ali
- Department of Biomedical Engineering and the Institute for Computational Medicine, Johns Hopkins UniversityBaltimore, MD, US
| | - Jialiu A Liang
- Department of Biomedical Engineering and the Institute for Computational Medicine, Johns Hopkins UniversityBaltimore, MD, US
| | - Shijie Zhou
- Department of Biomedical Engineering and the Institute for Computational Medicine, Johns Hopkins UniversityBaltimore, MD, US
| | - Natalia A Trayanova
- Department of Biomedical Engineering and the Institute for Computational Medicine, Johns Hopkins UniversityBaltimore, MD, US
| |
Collapse
|
27
|
Roney CH, Pashaei A, Meo M, Dubois R, Boyle PM, Trayanova NA, Cochet H, Niederer SA, Vigmond EJ. Universal atrial coordinates applied to visualisation, registration and construction of patient specific meshes. Med Image Anal 2019; 55:65-75. [PMID: 31026761 PMCID: PMC6543067 DOI: 10.1016/j.media.2019.04.004] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2018] [Revised: 03/07/2019] [Accepted: 04/15/2019] [Indexed: 11/26/2022]
Abstract
We introduce a coordinate system for the atria based on anatomical landmarks. We construct the coordinates from solutions to Laplace’s equation. We demonstrate the mapping of both scalar and vector data between different atria. The coordinate system was used for registration and 2D visualisation of multimodal data. Patient specific meshes with atrial structures and fibre direction were constructed using just five landmark points.
Integrating spatial information about atrial physiology and anatomy in a single patient from multimodal datasets, as well as generalizing these data across patients, requires a common coordinate system. In the atria, this is challenging due to the complexity and variability of the anatomy. We aimed to develop and validate a Universal Atrial Coordinate (UAC) system for the following applications: combination and assessment of multimodal data; comparison of spatial data across patients; 2D visualization; and construction of patient specific geometries to test mechanistic hypotheses. Left and right atrial LGE-MRI data were segmented and meshed. Two coordinates were calculated for each atrium by solving Laplace’s equation, with boundary conditions assigned using five landmark points. The coordinate system was used to map spatial information between atrial meshes, including scalar fields measured using different mapping modalities, and atrial anatomic structures and fibre directions from a reference geometry. Average error in point transfer from a source mesh to a destination mesh and back again was less than 0.1 mm for the left atrium and 0.02 mm for the right atrium. Patient specific meshes were constructed using the coordinate system and phase singularity density maps from arrhythmia simulations were visualised in 2D. In conclusion, we have developed a universal atrial coordinate system allowing automatic registration of imaging and electroanatomic mapping data, 2D visualisation, and patient specific model creation.
Collapse
Affiliation(s)
- Caroline H Roney
- School of Biomedical Engineering and Imaging Sciences, King's College London, United Kingdom.
| | - Ali Pashaei
- Institute of Electrophysiology and Heart Modeling (IHU Liryc), Foundation Bordeaux University, Pessac-Bordeaux, France; IMB Bordeaux Institute of Mathematics, University of Bordeaux, 351 cours de la Libération, Talence 33405, France
| | - Marianna Meo
- Institute of Electrophysiology and Heart Modeling (IHU Liryc), Foundation Bordeaux University, Pessac-Bordeaux, France; University of Bordeaux, CRCTB, U1045, Bordeaux, France; INSERM, CRCTB, U1045, Bordeaux, France
| | - Rémi Dubois
- Institute of Electrophysiology and Heart Modeling (IHU Liryc), Foundation Bordeaux University, Pessac-Bordeaux, France; University of Bordeaux, CRCTB, U1045, Bordeaux, France; INSERM, CRCTB, U1045, Bordeaux, France
| | | | | | - Hubert Cochet
- Institute of Electrophysiology and Heart Modeling (IHU Liryc), Foundation Bordeaux University, Pessac-Bordeaux, France
| | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, United Kingdom
| | - Edward J Vigmond
- Institute of Electrophysiology and Heart Modeling (IHU Liryc), Foundation Bordeaux University, Pessac-Bordeaux, France; IMB Bordeaux Institute of Mathematics, University of Bordeaux, 351 cours de la Libération, Talence 33405, France
| |
Collapse
|
28
|
Aronis KN, Ali R, Trayanova NA. The role of personalized atrial modeling in understanding atrial fibrillation mechanisms and improving treatment. Int J Cardiol 2019; 287:139-147. [PMID: 30755334 DOI: 10.1016/j.ijcard.2019.01.096] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/28/2018] [Revised: 01/24/2019] [Accepted: 01/28/2019] [Indexed: 12/13/2022]
Abstract
Atrial fibrillation is the most common arrhythmia in humans and is associated with high morbidity, mortality and health-related expenses. Computational approaches have been increasingly utilized in atrial electrophysiology. In this review we summarize the recent advancements in atrial fibrillation modeling at the organ scale. Multi-scale atrial models now incorporate high level detail of atrial anatomy, tissue ultrastructure and fibrosis distribution. We provide the state-of-the art methodologies in developing personalized atrial fibrillation models with realistic geometry and tissue properties. We then focus on the use of multi-scale atrial models to gain mechanistic insights in AF. Simulations using atrial models have provided important insight in the mechanisms underlying AF, showing the importance of the atrial fibrotic substrate and altered atrial electrophysiology in initiation and maintenance of AF. Last, we summarize the translational evidence that supports incorporation of computational modeling in clinical practice for development of personalized treatment strategies for patients with AF. In early-stages clinical studies, AF models successfully identify patients where pulmonary vein isolation alone is not adequate for treatment of AF and suggest novel targets for ablation. We conclude with a summary of the future developments envisioned for the field of atrial computational electrophysiology.
Collapse
Affiliation(s)
- Konstantinos N Aronis
- Department of Biomedical Engineering and the Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA; Division of Cardiology, Johns Hopkins Hospital, Baltimore, MD, USA
| | - Rheeda Ali
- Department of Biomedical Engineering and the Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA
| | - Natalia A Trayanova
- Department of Biomedical Engineering and the Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA.
| |
Collapse
|
29
|
Loewe A, Poremba E, Oesterlein T, Luik A, Schmitt C, Seemann G, Dössel O. Patient-Specific Identification of Atrial Flutter Vulnerability-A Computational Approach to Reveal Latent Reentry Pathways. Front Physiol 2019; 9:1910. [PMID: 30692934 PMCID: PMC6339942 DOI: 10.3389/fphys.2018.01910] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2018] [Accepted: 12/18/2018] [Indexed: 11/23/2022] Open
Abstract
Atypical atrial flutter (AFlut) is a reentrant arrhythmia which patients frequently develop after ablation for atrial fibrillation (AF). Indeed, substrate modifications during AF ablation can increase the likelihood to develop AFlut and it is clinically not feasible to reliably and sensitively test if a patient is vulnerable to AFlut. Here, we present a novel method based on personalized computational models to identify pathways along which AFlut can be sustained in an individual patient. We build a personalized model of atrial excitation propagation considering the anatomy as well as the spatial distribution of anisotropic conduction velocity and repolarization characteristics based on a combination of a priori knowledge on the population level and information derived from measurements performed in the individual patient. The fast marching scheme is employed to compute activation times for stimuli from all parts of the atria. Potential flutter pathways are then identified by tracing loops from wave front collision sites and constricting them using a geometric snake approach under consideration of the heterogeneous wavelength condition. In this way, all pathways along which AFlut can be sustained are identified. Flutter pathways can be instantiated by using an eikonal-diffusion phase extrapolation approach and a dynamic multifront fast marching simulation. In these dynamic simulations, the initial pattern eventually turns into the one driven by the dominant pathway, which is the only pathway that can be observed clinically. We assessed the sensitivity of the flutter pathway maps with respect to conduction velocity and its anisotropy. Moreover, we demonstrate the application of tailored models considering disease-specific repolarization properties (healthy, AF-remodeled, potassium channel mutations) as well as applicabiltiy on a clinical dataset. Finally, we tested how AFlut vulnerability of these substrates is modulated by exemplary antiarrhythmic drugs (amiodarone, dronedarone). Our novel method allows to assess the vulnerability of an individual patient to develop AFlut based on the personal anatomical, electrophysiological, and pharmacological characteristics. In contrast to clinical electrophysiological studies, our computational approach provides the means to identify all possible AFlut pathways and not just the currently dominant one. This allows to consider all relevant AFlut pathways when tailoring clinical ablation therapy in order to reduce the development and recurrence of AFlut.
Collapse
Affiliation(s)
- Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Emanuel Poremba
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Tobias Oesterlein
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Armin Luik
- Medizinische Klinik IV, Städtisches Klinikum Karlsruhe, Karlsruhe, Germany
| | - Claus Schmitt
- Medizinische Klinik IV, Städtisches Klinikum Karlsruhe, Karlsruhe, Germany
| | - Gunnar Seemann
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
- Institute for Experimental Cardiovascular Medicine, University Heart Center Freiburg Bad Krozingen, Freiburg, Germany
- Faculty of Medicine, Albert-Ludwigs University, Freiburg, Germany
| | - Olaf Dössel
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| |
Collapse
|
30
|
Kim IS, Lim B, Shim J, Hwang M, Yu HT, Kim TH, Uhm JS, Kim SH, Joung B, On YK, Oh S, Oh YS, Nam GB, Lee MH, Shim EB, Kim YH, Pak HN. Clinical Usefulness of Computational Modeling-Guided Persistent Atrial Fibrillation Ablation: Updated Outcome of Multicenter Randomized Study. Front Physiol 2019; 10:1512. [PMID: 31920716 PMCID: PMC6928133 DOI: 10.3389/fphys.2019.01512] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Accepted: 12/02/2019] [Indexed: 01/13/2023] Open
Abstract
OBJECTIVE Catheter ablation of persistent atrial fibrillation (AF) is still challenging, no optimal extra-pulmonary vein lesion set is known. We previously reported the clinical feasibility of computational modeling-guided AF catheter ablation. METHODS We randomly assigned 118 patients with persistent AF (77.8% men, age 60.8 ± 9.9 years) to the computational modeling-guided ablation group (53 patients) and the empirical ablation group (55 patients) based on the operators' experience. For virtual ablation, four virtual linear and one electrogram-guided lesion sets were tested on patient heart computed tomogram-based models, and the lesion set with the fastest termination time was reported to the operator in the modeling-guided ablation group. The primary outcome was freedom from atrial tachyarrhythmias lasting longer than 30 s after a single procedure. RESULTS During 31.5 ± 9.4 months, virtual ablation procedures were available in 95.2% of the patients (108/118). Clinical recurrence rate was significantly lower after a modeling-guided ablation than after an empirical ablation (20.8 vs. 40.0%, log-rank p = 0.042). Modeling-guided ablation was independently associated with a better long-term rhythm outcome of persistent AF ablation (HR = 0.29 [0.12-0.69], p = 0.005). The rhythm outcome of the modeling-guided ablation showed better trends in males, non-obese patients with a less remodeled atrium (left atrial dimension < 50 mm), ejection fraction ≥ 50%, and those without hypertension or diabetes (p < 0.01). There were no significant differences between the groups for the total procedure time (p = 0.403), ablation time (p = 0.510), and major complication rate (p = 0.900). CONCLUSION Among patients with persistent AF, the computational modeling-guided ablation was superior to the empirical catheter ablation regarding the rhythm outcome. CLINICAL TRIAL REGISTRATION This study was registered with the ClinicalTrials.gov, number NCT02171364.
Collapse
Affiliation(s)
- In-Soo Kim
- Yonsei University Health System, Seoul, South Korea
| | | | - Jaemin Shim
- Cardiovascular Center, Korea University, Seoul, South Korea
| | - Minki Hwang
- Yonsei University Health System, Seoul, South Korea
| | - Hee Tae Yu
- Yonsei University Health System, Seoul, South Korea
| | - Tae-Hoon Kim
- Yonsei University Health System, Seoul, South Korea
| | - Jae-Sun Uhm
- Yonsei University Health System, Seoul, South Korea
| | - Sung-Hwan Kim
- Division of Cardiology, Catholic University of Korea, Seoul, South Korea
| | | | - Young Keun On
- Samsung Medical Center, Sungkyunkwan University, Seoul, South Korea
| | - Seil Oh
- Division of Cardiology, Seoul National University, Seoul, South Korea
| | - Yong-Seog Oh
- Division of Cardiology, Catholic University of Korea, Seoul, South Korea
| | - Gi-Byung Nam
- Asan Medical Center, University of Ulsan, Ulsan, South Korea
| | | | - Eun Bo Shim
- Department of Mechanical and Biomedical Engineering, Kangwon National University, Chuncheon, South Korea
| | - Young-Hoon Kim
- Cardiovascular Center, Korea University, Seoul, South Korea
| | - Hui-Nam Pak
- Yonsei University Health System, Seoul, South Korea
- *Correspondence: Hui-Nam Pak,
| | | |
Collapse
|
31
|
Roney CH, Whitaker J, Sim I, O'Neill L, Mukherjee RK, Razeghi O, Vigmond EJ, Wright M, O'Neill MD, Williams SE, Niederer SA. A technique for measuring anisotropy in atrial conduction to estimate conduction velocity and atrial fibre direction. Comput Biol Med 2019; 104:278-290. [PMID: 30415767 PMCID: PMC6506689 DOI: 10.1016/j.compbiomed.2018.10.019] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2018] [Revised: 10/17/2018] [Accepted: 10/17/2018] [Indexed: 01/04/2023]
Abstract
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.
Collapse
Affiliation(s)
- Caroline H Roney
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom.
| | - John Whitaker
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Iain Sim
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Louisa O'Neill
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Rahul K Mukherjee
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Orod Razeghi
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Edward J Vigmond
- LIRYC Electrophysiology and Heart Modeling Institute, Campus Xavier Arnozan, Avenue du Haut Lévêque, 33600, Pessac, France; Univ. Bordeaux, IMB, UMR 5251, F-33400, Talence, France
| | - Matthew Wright
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Mark D O'Neill
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Steven E Williams
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| | - Steven A Niederer
- School of Imaging Sciences and Biomedical Engineering, King's College London, London, United Kingdom
| |
Collapse
|
32
|
Irakoze É, Jacquemet V. Simulated P wave morphology in the presence of endo-epicardial activation delay. Europace 2018; 20:iii16-iii25. [PMID: 30476058 DOI: 10.1093/europace/euy229] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2018] [Accepted: 09/18/2018] [Indexed: 11/14/2022] Open
Abstract
Aims Evidences of asynchrony between epicardial and endocardial activation in the atrial wall have been reported. We used a computer model of the atria and torso to investigate the consequences of such activation delay on P wave morphology, while controlling for P wave duration. Methods and results We created 390 models of the atria based on the same geometry. These models differed by atrial wall thickness (from 2 to 3 mm), transmural coupling, and tissue conductivity in the endocardial and epicardial layers. Among them, 18 were in baseline, 186 had slower conduction in the epicardium layer and 186 in the endocardial layer. Conduction properties were adjusted in such a way that total activation time was the same in all models. P waves on a 16-lead system were simulated during sinus rhythm. Activation maps were similar in all cases. Endo-epicardial delay varied between -5.5 and 5.5 ms vs. 0 ± 0.5 ms in baseline. All P waves had the same duration but variability in their morphology was observed. With slower epicardial conduction, P wave amplitude was reduced by an average of 20% on leads V3-V5 and P wave area decreased by 50% on leads V1-V2 and by 40% on lead V3. Reversed, lower magnitude effects were observed with slower endocardial conduction. Conclusion An endo-epicardial delay of a few milliseconds is sufficient to significantly alter P wave morphology, even if the activation map remains the same.
Collapse
Affiliation(s)
- Éric Irakoze
- Département de Pharmacologie et Physiologie, Institut de Génie Biomédical, Université de Montréal, Montréal, QC, Canada
- Centre de Recherche, Hôpital du Sacré-Cœur de Montréal, 5400 boul. Gouin Ouest, Montréal, QC, Canada
| | - Vincent Jacquemet
- Département de Pharmacologie et Physiologie, Institut de Génie Biomédical, Université de Montréal, Montréal, QC, Canada
- Centre de Recherche, Hôpital du Sacré-Cœur de Montréal, 5400 boul. Gouin Ouest, Montréal, QC, Canada
| |
Collapse
|
33
|
Roney CH, Ng FS, Debney MT, Eichhorn C, Nachiappan A, Chowdhury RA, Qureshi NA, Cantwell CD, Tweedy JH, Niederer SA, Peters NS, Vigmond EJ. Determinants of new wavefront locations in cholinergic atrial fibrillation. Europace 2018; 20:iii3-iii15. [PMID: 30476057 PMCID: PMC6251188 DOI: 10.1093/europace/euy235] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2018] [Accepted: 10/10/2018] [Indexed: 01/10/2023] Open
Abstract
AIMS Atrial fibrillation (AF) wavefront dynamics are complex and difficult to interpret, contributing to uncertainty about the mechanisms that maintain AF. We aimed to investigate the interplay between rotors, wavelets, and focal sources during fibrillation. METHODS AND RESULTS Arrhythmia wavefront dynamics were analysed for four optically mapped canine cholinergic AF preparations. A bilayer computer model was tuned to experimental preparations, and varied to have (i) fibrosis in both layers or the epicardium only, (ii) different spatial acetylcholine distributions, (iii) different intrinsic action potential duration between layers, and (iv) varied interlayer connectivity. Phase singularities (PSs) were identified and tracked over time to identify rotational drivers. New focal wavefronts were identified using phase contours. Phase singularity density and new wavefront locations were calculated during AF. There was a single dominant mechanism for sustaining AF in each of the preparations, either a rotational driver or repetitive new focal wavefronts. High-density PS sites existed preferentially around the pulmonary vein junctions. Three of the four preparations exhibited stable preferential sites of new wavefronts. Computational simulations predict that only a small number of connections are functionally important in sustaining AF, with new wavefront locations determined by the interplay between fibrosis distribution, acetylcholine concentration, and heterogeneity in repolarization within layers. CONCLUSION We were able to identify preferential sites of new wavefront initiation and rotational activity, in order to determine the mechanisms sustaining AF. Electrical measurements should be interpreted differently according to whether they are endocardial or epicardial recordings.
Collapse
Affiliation(s)
- Caroline H Roney
- ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London, UK
- LIRYC Electrophysiology and Heart Modeling Institute, Bordeaux Fondation, Avenue du Haut-Lévèque, Pessac, France
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Fu Siong Ng
- ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Michael T Debney
- ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Christian Eichhorn
- ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Arun Nachiappan
- ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Rasheda A Chowdhury
- ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Norman A Qureshi
- ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Chris D Cantwell
- ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Jennifer H Tweedy
- ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Nicholas S Peters
- ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London, UK
| | - Edward J Vigmond
- LIRYC Electrophysiology and Heart Modeling Institute, Bordeaux Fondation, Avenue du Haut-Lévèque, Pessac, France
- Univ. Bordeaux, IMB UMR 5251, F-33400 Talence, France
| |
Collapse
|
34
|
Rossi S, Gaeta S, Griffith BE, Henriquez CS. Muscle Thickness and Curvature Influence Atrial Conduction Velocities. Front Physiol 2018; 9:1344. [PMID: 30420809 PMCID: PMC6215968 DOI: 10.3389/fphys.2018.01344] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2018] [Accepted: 09/06/2018] [Indexed: 12/04/2022] Open
Abstract
Electroanatomical mapping is currently used to provide clinicians with information about the electrophysiological state of the heart and to guide interventions like ablation. These maps can be used to identify ectopic triggers of an arrhythmia such as atrial fibrillation (AF) or changes in the conduction velocity (CV) that have been associated with poor cell to cell coupling or fibrosis. Unfortunately, many factors are known to affect CV, including membrane excitability, pacing rate, wavefront curvature, and bath loading, making interpretation challenging. In this work, we show how endocardial conduction velocities are also affected by the geometrical factors of muscle thickness and wall curvature. Using an idealized three-dimensional strand, we show that transverse conductivities and boundary conditions can slow down or speed up signal propagation, depending on the curvature of the muscle tissue. In fact, a planar wavefront that is parallel to a straight line normal to the mid-surface does not remain normal to the mid-surface in a curved domain. We further demonstrate that the conclusions drawn from the idealized test case can be used to explain spatial changes in conduction velocities in a patient-specific reconstruction of the left atrial posterior wall. The simulations suggest that the widespread assumption of treating atrial muscle as a two-dimensional manifold for electrophysiological simulations will not accurately represent the endocardial conduction velocities in regions of the heart thicker than 0.5 mm with significant wall curvature.
Collapse
Affiliation(s)
- Simone Rossi
- Cardiovascular Modeling and Simulation Laboratory, Carolina Center for Interdisciplinary Applied Mathematics, University of North Carolina, Chapel Hill, NC, United States
| | - Stephen Gaeta
- Clinical Cardiac Electrophysiology/Cardiology Division, Duke University Medical Center, Durham, NC, United States
| | - Boyce E. Griffith
- Cardiovascular Modeling and Simulation Laboratory, Carolina Center for Interdisciplinary Applied Mathematics, University of North Carolina, Chapel Hill, NC, United States
- Departments of Mathematics, Applied Physical Sciences, and Biomedical Engineering, University of North Carolina, Chapel Hill, NC, United States
- McAllister Heart Institute, University of North Carolina, Chapel Hill, NC, United States
| | - Craig S. Henriquez
- Department of Biomedical Engineering, Pratt School of Engineering, Duke University, Durham, NC, United States
| |
Collapse
|
35
|
Saha M, Roney CH, Bayer JD, Meo M, Cochet H, Dubois R, Vigmond EJ. Wavelength and Fibrosis Affect Phase Singularity Locations During Atrial Fibrillation. Front Physiol 2018; 9:1207. [PMID: 30246796 PMCID: PMC6139329 DOI: 10.3389/fphys.2018.01207] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2018] [Accepted: 08/10/2018] [Indexed: 01/06/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Mirabeau Saha
- IMB, UMR 5251, University of Bordeaux, Pessac, France
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux University, Pessac, France
| | - Caroline H. Roney
- Department of Biomedical Engineering, King's College London, London, United Kingdom
| | - Jason D. Bayer
- IMB, UMR 5251, University of Bordeaux, Pessac, France
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux University, Pessac, France
| | - Marianna Meo
- IMB, UMR 5251, University of Bordeaux, Pessac, France
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux University, Pessac, France
| | - Hubert Cochet
- IMB, UMR 5251, University of Bordeaux, Pessac, France
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux University, Pessac, France
| | - Remi Dubois
- IMB, UMR 5251, University of Bordeaux, Pessac, France
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux University, Pessac, France
| | - Edward J. Vigmond
- IMB, UMR 5251, University of Bordeaux, Pessac, France
- IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux University, Pessac, France
| |
Collapse
|
36
|
Vidmar D, Alhusseini MI, Narayan SM, Rappel WJ. Characterizing Electrogram Signal Fidelity and the Effects of Signal Contamination on Mapping Human Persistent Atrial Fibrillation. Front Physiol 2018; 9:1232. [PMID: 30237766 PMCID: PMC6135945 DOI: 10.3389/fphys.2018.01232] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2018] [Accepted: 08/15/2018] [Indexed: 11/30/2022] Open
Abstract
Objective: Determining accurate intracardiac maps of atrial fibrillation (AF) in humans can be difficult, owing primarily to various sources of contamination in electrogram signals. The goal of this study is to develop a measure for signal fidelity and to develop methods to quantify robustness of observed rotational activity in phase maps subject to signal contamination. Methods: We identified rotational activity in phase maps of human persistent AF using the Hilbert transform of sinusoidally recomposed signals, where localized ablation at rotational sites terminated fibrillation. A novel measure of signal fidelity was developed to quantify signal quality. Contamination is then introduced to the underlying electrograms by removing signals at random, adding noise to computations of cycle length, and adding realistic far-field signals. Mean tip number N and tip density δ, defined as the proportion of time a region contains a tip, at the termination site are computed to compare the effects of contamination. Results: Domains of low signal fidelity correspond to the location of rotational cores. Removing signals and altering cycle length accounted for minor changes in tip density, while targeted removal of low fidelity electrograms can result in a significant increase in tip density and stability. Far-field contamination was found to obscure rotation at the termination site. Conclusion: Rotational activity in clinical AF can produce domains of low fidelity electrogram recordings at rotational cores. Observed rotational patterns in phase maps appear most sensitive to far-field activation. These results may inform novel methods to map AF in humans which can be tested directly in patients at electrophysiological study and ablation.
Collapse
Affiliation(s)
- David Vidmar
- Department of Physics, University of California, San Diego, San Diego, CA, United States
| | - Mahmood I. Alhusseini
- Division of Cardiology, Department of Medicine, Stanford University, Stanford, CA, United States
| | - Sanjiv M. Narayan
- Division of Cardiology, Department of Medicine, Stanford University, Stanford, CA, United States
| | - Wouter-Jan Rappel
- Department of Physics, University of California, San Diego, San Diego, CA, United States
- *Correspondence: Wouter-Jan Rappel
| |
Collapse
|
37
|
Roney CH, Bayer JD, Cochet H, Meo M, Dubois R, Jaïs P, Vigmond EJ. Variability in pulmonary vein electrophysiology and fibrosis determines arrhythmia susceptibility and dynamics. PLoS Comput Biol 2018; 14:e1006166. [PMID: 29795549 PMCID: PMC5997352 DOI: 10.1371/journal.pcbi.1006166] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2017] [Revised: 06/12/2018] [Accepted: 04/30/2018] [Indexed: 11/28/2022] Open
Abstract
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. Atrial fibrillation is the most commonly encountered cardiac arrhythmia, affecting a significant portion of the population. Currently, ablation is the most effective treatment but success rates are less than optimal, being 70% one-year post-treatment. There is a large effort to find better ablation strategies to permanently cure the condition. Pulmonary vein isolation by ablation is more or less the standard of care, but many questions remain since pulmonary vein ectopy by itself does not explain all of the clinical successes or failures. We used computer simulations to investigate how electrophysiological properties of the pulmonary veins can affect rotor formation and maintenance in patients suffering from atrial fibrillation. We used complex, biophysical representations of cellular electrophysiology in highly detailed geometries constructed from patient scans. We heterogeneously varied electrophysiological and structural properties to see their effects on rotor initiation and maintenance. Our study suggests a metric for indicating the likelihood of success of pulmonary vein isolation. Thus either measuring this clinically, or running patient-specific simulations to estimate this metric may suggest whether ablation in addition to pulmonary vein isolation should be performed. Our study provides motivation for a retrospective clinical study or experimental study into this metric.
Collapse
Affiliation(s)
- Caroline H. Roney
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France
- Division of Imaging Sciences and Biomedical Engineering, King’s College London, London, United Kingdom
| | - Jason D. Bayer
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France
- Univ. Bordeaux, IMB UMR 5251, CNRS, F-33400 Talence, France
| | - Hubert Cochet
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France
- Hôpital Cardiologique du Haut-L’évêque, Université de Bordeaux, LIRYC Institute: IHU LIRYC ANR-10-IAHU-04 and Equipex MUSIC ANR-11-EQPX-0030, Bordeaux, France
| | - Marianna Meo
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France
| | - Rémi Dubois
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France
| | - Pierre Jaïs
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France
- Hôpital Cardiologique du Haut-L’évêque, Université de Bordeaux, LIRYC Institute: IHU LIRYC ANR-10-IAHU-04 and Equipex MUSIC ANR-11-EQPX-0030, Bordeaux, France
| | - Edward J. Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France
- Univ. Bordeaux, IMB UMR 5251, CNRS, F-33400 Talence, France
- * E-mail:
| |
Collapse
|
38
|
Lim B, Hwang M, Song JS, Ryu AJ, Joung B, Shim EB, Ryu H, Pak HN. Effectiveness of atrial fibrillation rotor ablation is dependent on conduction velocity: An in-silico 3-dimensional modeling study. PLoS One 2017; 12:e0190398. [PMID: 29287119 PMCID: PMC5747478 DOI: 10.1371/journal.pone.0190398] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2017] [Accepted: 12/14/2017] [Indexed: 11/29/2022] Open
Abstract
BACKGROUND We previously reported that stable rotors are observed in in-silico human atrial fibrillation (AF) models, and are well represented by a dominant frequency (DF). In the current study, we hypothesized that the outcome of DF ablation is affected by conduction velocity (CV) conditions and examined this hypothesis using in-silico 3D-AF modeling. METHODS We integrated 3D CT images of left atrium obtained from 10 patients with persistent AF (80% male, 61.8±13.5 years old) into in-silico AF model. We compared AF maintenance durations (max 300s), spatiotemporal stabilities of DF, phase singularity (PS) number, life-span of PS, and AF termination or defragmentation rates after virtual DF ablation with 5 different CV conditions (0.2, 0.3, 0.4, 0.5, and 0.6m/s). RESULTS 1. AF maintenance duration (p<0.001), spatiotemporal mean variance of DF (p<0.001), and the number of PS (p = 0.023) showed CV dependent bimodal patterns (highest at CV0.4m/s and lowest at CV0.6m/s) consistently. 2. After 10% highest DF ablation, AF defragmentation rates were the lowest at CV0.4m/s (37.8%), but highest at CV0.5 and 0.6m/s (all 100%, p<0.001). 3. In the episodes with AF termination or defragmentation followed by 10% highest DF ablation, baseline AF maintenance duration was shorter (p<0.001), spatiotemporal mean variance of DF was lower (p = 0.014), and the number of PS was lower (p = 0.004) than those with failed AF defragmentation after DF ablation. CONCLUSION Virtual ablation of DF, which may indicate AF driver, was more likely to terminate or defragment AF with spatiotemporally stable DF, but not likely to do so in long-lasting and sustained AF conditions, depending on CV.
Collapse
Affiliation(s)
- Byounghyun Lim
- Yonsei University Health System, Seoul, Republic of Korea
| | - Minki Hwang
- Yonsei University Health System, Seoul, Republic of Korea
| | - Jun-Seop Song
- Yonsei University Health System, Seoul, Republic of Korea
| | - Ah-Jin Ryu
- Department of Mechanical and Biomedical Engineering, Kangwon National University, Chuncheon, Ganwon-do, Republic of Korea
| | - Boyoung Joung
- Yonsei University Health System, Seoul, Republic of Korea
| | - Eun Bo Shim
- Department of Mechanical and Biomedical Engineering, Kangwon National University, Chuncheon, Ganwon-do, Republic of Korea
| | - Hyungon Ryu
- NVIDIA, Yonsei University, Department of Mathematics, Seoul, Republic of Korea
| | - Hui-Nam Pak
- Yonsei University Health System, Seoul, Republic of Korea
| |
Collapse
|
39
|
Campana C, Akar FG. Commentary: Atrial Fibrillation Dynamics and Ionic Block Effects in Six Heterogeneous Human 3D Virtual Atria with Distinct Repolarization Dynamics. Front Bioeng Biotechnol 2017; 5:59. [PMID: 29057224 PMCID: PMC5635327 DOI: 10.3389/fbioe.2017.00059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2017] [Accepted: 09/20/2017] [Indexed: 11/26/2022] Open
Affiliation(s)
- Chiara Campana
- Icahn School of Medicine at Mount Sinai, The Cardiovascular Institute, New York, NY, United States
| | - Fadi G Akar
- Icahn School of Medicine at Mount Sinai, The Cardiovascular Institute, New York, NY, United States
| |
Collapse
|
40
|
Roney CH, Cantwell CD, Bayer JD, Qureshi NA, Lim PB, Tweedy JH, Kanagaratnam P, Peters NS, Vigmond EJ, Ng FS. Spatial Resolution Requirements for Accurate Identification of Drivers of Atrial Fibrillation. Circ Arrhythm Electrophysiol 2017; 10:e004899. [PMID: 28500175 PMCID: PMC5434962 DOI: 10.1161/circep.116.004899] [Citation(s) in RCA: 85] [Impact Index Per Article: 12.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/20/2016] [Accepted: 04/11/2017] [Indexed: 11/25/2022]
Abstract
Supplemental Digital Content is available in the text. 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.
Collapse
Affiliation(s)
- Caroline H Roney
- From the ElectroCardioMaths Programme (C.H.R., C.D.C., N.A.Q., P.B.L., P.K., N.S.P., F.S.N.), and the Department of Bioengineering (J.H.T.), Imperial College London, United Kingdom; IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac-Bordeaux, France (J.D.B., E.J.V.); and Université de Bordeaux, IMB, UMR 5251, Talence, France (J.D.B., E.J.V.)
| | - Chris D Cantwell
- From the ElectroCardioMaths Programme (C.H.R., C.D.C., N.A.Q., P.B.L., P.K., N.S.P., F.S.N.), and the Department of Bioengineering (J.H.T.), Imperial College London, United Kingdom; IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac-Bordeaux, France (J.D.B., E.J.V.); and Université de Bordeaux, IMB, UMR 5251, Talence, France (J.D.B., E.J.V.)
| | - Jason D Bayer
- From the ElectroCardioMaths Programme (C.H.R., C.D.C., N.A.Q., P.B.L., P.K., N.S.P., F.S.N.), and the Department of Bioengineering (J.H.T.), Imperial College London, United Kingdom; IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac-Bordeaux, France (J.D.B., E.J.V.); and Université de Bordeaux, IMB, UMR 5251, Talence, France (J.D.B., E.J.V.)
| | - Norman A Qureshi
- From the ElectroCardioMaths Programme (C.H.R., C.D.C., N.A.Q., P.B.L., P.K., N.S.P., F.S.N.), and the Department of Bioengineering (J.H.T.), Imperial College London, United Kingdom; IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac-Bordeaux, France (J.D.B., E.J.V.); and Université de Bordeaux, IMB, UMR 5251, Talence, France (J.D.B., E.J.V.)
| | - Phang Boon Lim
- From the ElectroCardioMaths Programme (C.H.R., C.D.C., N.A.Q., P.B.L., P.K., N.S.P., F.S.N.), and the Department of Bioengineering (J.H.T.), Imperial College London, United Kingdom; IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac-Bordeaux, France (J.D.B., E.J.V.); and Université de Bordeaux, IMB, UMR 5251, Talence, France (J.D.B., E.J.V.)
| | - Jennifer H Tweedy
- From the ElectroCardioMaths Programme (C.H.R., C.D.C., N.A.Q., P.B.L., P.K., N.S.P., F.S.N.), and the Department of Bioengineering (J.H.T.), Imperial College London, United Kingdom; IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac-Bordeaux, France (J.D.B., E.J.V.); and Université de Bordeaux, IMB, UMR 5251, Talence, France (J.D.B., E.J.V.)
| | - Prapa Kanagaratnam
- From the ElectroCardioMaths Programme (C.H.R., C.D.C., N.A.Q., P.B.L., P.K., N.S.P., F.S.N.), and the Department of Bioengineering (J.H.T.), Imperial College London, United Kingdom; IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac-Bordeaux, France (J.D.B., E.J.V.); and Université de Bordeaux, IMB, UMR 5251, Talence, France (J.D.B., E.J.V.)
| | - Nicholas S Peters
- From the ElectroCardioMaths Programme (C.H.R., C.D.C., N.A.Q., P.B.L., P.K., N.S.P., F.S.N.), and the Department of Bioengineering (J.H.T.), Imperial College London, United Kingdom; IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac-Bordeaux, France (J.D.B., E.J.V.); and Université de Bordeaux, IMB, UMR 5251, Talence, France (J.D.B., E.J.V.).
| | - Edward J Vigmond
- From the ElectroCardioMaths Programme (C.H.R., C.D.C., N.A.Q., P.B.L., P.K., N.S.P., F.S.N.), and the Department of Bioengineering (J.H.T.), Imperial College London, United Kingdom; IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac-Bordeaux, France (J.D.B., E.J.V.); and Université de Bordeaux, IMB, UMR 5251, Talence, France (J.D.B., E.J.V.)
| | - Fu Siong Ng
- From the ElectroCardioMaths Programme (C.H.R., C.D.C., N.A.Q., P.B.L., P.K., N.S.P., F.S.N.), and the Department of Bioengineering (J.H.T.), Imperial College London, United Kingdom; IHU Liryc, Electrophysiology and Heart Modeling Institute, Fondation Bordeaux Université, Pessac-Bordeaux, France (J.D.B., E.J.V.); and Université de Bordeaux, IMB, UMR 5251, Talence, France (J.D.B., E.J.V.)
| |
Collapse
|
41
|
Roney CH, Bayer JD, Zahid S, Meo M, Boyle PMJ, Trayanova NA, Haïssaguerre M, Dubois R, Cochet H, Vigmond EJ. Modelling methodology of atrial fibrosis affects rotor dynamics and electrograms. Europace 2017; 18:iv146-iv155. [PMID: 28011842 DOI: 10.1093/europace/euw365] [Citation(s) in RCA: 75] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2016] [Accepted: 06/28/2016] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Caroline H Roney
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France.,Univ. Bordeaux, IMB UMR 5251, CNRS, F-33400 Talence, France
| | - Jason D Bayer
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France.,Univ. Bordeaux, IMB UMR 5251, CNRS, F-33400 Talence, France
| | - Sohail Zahid
- Department of Biomedical Engineering and Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA
| | - Marianna Meo
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France.,Univ. Bordeaux, Centre de recherche Cardio-Thoracique de Bordeaux, U1045, F- 33000, Bordeaux, France
| | - Patrick M J Boyle
- Department of Biomedical Engineering and Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA
| | - Natalia A Trayanova
- Department of Biomedical Engineering and Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA
| | - Michel Haïssaguerre
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France.,Univ. Bordeaux, Centre de recherche Cardio-Thoracique de Bordeaux, U1045, F- 33000, Bordeaux, France.,Bordeaux University Hospital (CHU), Electrophysiology and Ablation Unit, F-33600 Pessac, France
| | - Rémi Dubois
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France.,Univ. Bordeaux, Centre de recherche Cardio-Thoracique de Bordeaux, U1045, F- 33000, Bordeaux, France
| | - Hubert Cochet
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France.,Univ. Bordeaux, Centre de recherche Cardio-Thoracique de Bordeaux, U1045, F- 33000, Bordeaux, France.,Department of Cardiac Imaging, Bordeaux University Hospital, Bordeaux, France
| | - Edward J Vigmond
- IHU Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac- Bordeaux, France.,Univ. Bordeaux, IMB UMR 5251, CNRS, F-33400 Talence, France
| |
Collapse
|
42
|
Green HD, Thomas G, Terry JR. Signal Reconstruction of Pulmonary Vein Recordings Using a Phenomenological Mathematical Model: Application to Pulmonary Vein Isolation Therapy. Front Physiol 2017; 8:496. [PMID: 28769816 PMCID: PMC5511835 DOI: 10.3389/fphys.2017.00496] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2016] [Accepted: 06/28/2017] [Indexed: 12/19/2022] Open
Abstract
Atrial fibrillation (AF), the most prevalent cardiac arrhythmia, is commonly initiated by ectopic beats originating from a small myocardial sleeve extending over the pulmonary veins. Pulmonary vein isolation therapy attempts to isolate the pulmonary veins from the left atrium by ablating tissue, commonly by using radiofrequency ablation. During this procedure, the cardiologist records electrical activity using a lasso catheter, and the activation pattern recorded is used as a guide toward which regions to ablate. However, poor contact between electrode and tissue can lead to important regions of electrical activity not being recorded in clinic. We reproduce these signals through the use of a phenomenological model of the cardiac action potential on a cylinder, which we fit to post-AF atrial cells, and model the bipolar electrodes of the lasso catheter by an approximation of the surface potential. The resulting activation pattern is validated by direct comparison with those of clinical recordings. A potential application of the model is to reconstruct the missing electrical activity, minimizing the impact of the information loss on the clinical procedure, and we present results to demonstrate this.
Collapse
Affiliation(s)
- Harry D Green
- College of Engineering, Mathematics and Physical Sciences, University of ExeterExeter, United Kingdom.,Wellcome Trust Centre for Biomedical Modelling and Analysis, University of ExeterExeter, United Kingdom.,Living Systems Institute, University of ExeterExeter, United Kingdom
| | - Glyn Thomas
- Bristol Heart InstituteBristol, United Kingdom
| | - John R Terry
- College of Engineering, Mathematics and Physical Sciences, University of ExeterExeter, United Kingdom.,Wellcome Trust Centre for Biomedical Modelling and Analysis, University of ExeterExeter, United Kingdom.,Living Systems Institute, University of ExeterExeter, United Kingdom.,EPSRC Centre for Predictive Modelling in Healthcare, University of ExeterExeter, United Kingdom
| |
Collapse
|
43
|
Fibrosis and Atrial Fibrillation: Computerized and Optical Mapping; A View into the Human Atria at Submillimeter Resolution. JACC Clin Electrophysiol 2017; 3:531-546. [PMID: 29159313 DOI: 10.1016/j.jacep.2017.05.002] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Abstract
Recent studies strongly suggest that the majority of atrial fibrillation (AF) patients with diagnosed or subclinical cardiac diseases have established or even pre-existing fibrotic structural remodeling, which may lead to conduction abnormalities and reentrant activity that sustain AF. As conventional treatments fail to treat AF in far too many cases, an urgent need exists to identify specific structural arrhythmogenic fibrosis patterns, which may maintain AF, in order to identify effective ablation targets for AF treatment. However, the existing challenge is to define what exact structural remodeling within the complex 3D human atrial wall is arrhythmogenic, as well as linking arrhythmogenic fibrosis to an underlying mechanism of AF maintenance in the clinical setting. This review is focused on the role of 3D fibrosis architecture in the mechanisms of AF maintenance revealed by submillimeter, high-resolution ex-vivo imaging modalities directly of human atria, as well as from in-silico 3D computational techniques that can be able to overcome in-vivo clinical limitations. The systematic integration of functional and structural imaging ex-vivo may inform the necessary integration of electrode and structural mapping in-vivo. A holistic view of AF driver mechanisms may begin to identify the defining characteristics or "fingerprints" of reentrant AF drivers, such as 3D fibrotic architecture, in order to design optimal patient-specific ablation strategies.
Collapse
|
44
|
Pashakhanloo F, Herzka DA, Ashikaga H, Mori S, Gai N, Bluemke DA, Trayanova NA, McVeigh ER. Myofiber Architecture of the Human Atria as Revealed by Submillimeter Diffusion Tensor Imaging. Circ Arrhythm Electrophysiol 2016; 9:e004133. [PMID: 27071829 DOI: 10.1161/circep.116.004133] [Citation(s) in RCA: 112] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/15/2016] [Accepted: 03/15/2016] [Indexed: 11/16/2022]
Abstract
BACKGROUND Accurate knowledge of the human atrial fibrous structure is paramount in understanding the mechanisms of atrial electric function in health and disease. Thus far, such knowledge has been acquired from destructive sectioning, and there is a paucity of data about atrial fiber architecture variability in the human population. METHODS AND RESULTS In this study, we have developed a customized 3-dimensional diffusion tensor magnetic resonance imaging sequence on a clinical scanner that makes it possible to image an entire intact human heart specimen ex vivo at submillimeter resolution. The data from 8 human atrial specimens obtained with this technique present complete maps of the fibrous organization of the human atria. The findings demonstrate that the main features of atrial anatomy are mostly preserved across subjects although the exact location and orientation of atrial bundles vary. Using the full tractography data, we were able to cluster, visualize, and characterize the distinct major bundles in the human atria. Furthermore, quantitative characterization of the fiber angles across the atrial wall revealed that the transmural fiber angle distribution is heterogeneous throughout different regions of the atria. CONCLUSIONS The application of submillimeter diffusion tensor magnetic resonance imaging provides an unprecedented level of information on both human atrial structure, as well as its intersubject variability. The high resolution and fidelity of this data could enhance our understanding of structural contributions to atrial rhythm and pump disorders and lead to improvements in their targeted treatment.
Collapse
Affiliation(s)
- Farhad Pashakhanloo
- From the Departments of Biomedical Engineering (F.P., D.A.H., N.A.T., E.R.M.), Medicine (H.A.), and Radiology (S.M., E.R.M), Johns Hopkins University, Baltimore, MD; Radiology and Imaging Sciences, National Institutes of Health Clinical Center, Bethesda, MD (N.G, D.A.B.); and Departments of Bioengineering, Medicine, and Radiology, University of California, San Diego (E.R.M.)
| | - Daniel A Herzka
- From the Departments of Biomedical Engineering (F.P., D.A.H., N.A.T., E.R.M.), Medicine (H.A.), and Radiology (S.M., E.R.M), Johns Hopkins University, Baltimore, MD; Radiology and Imaging Sciences, National Institutes of Health Clinical Center, Bethesda, MD (N.G, D.A.B.); and Departments of Bioengineering, Medicine, and Radiology, University of California, San Diego (E.R.M.)
| | - Hiroshi Ashikaga
- From the Departments of Biomedical Engineering (F.P., D.A.H., N.A.T., E.R.M.), Medicine (H.A.), and Radiology (S.M., E.R.M), Johns Hopkins University, Baltimore, MD; Radiology and Imaging Sciences, National Institutes of Health Clinical Center, Bethesda, MD (N.G, D.A.B.); and Departments of Bioengineering, Medicine, and Radiology, University of California, San Diego (E.R.M.)
| | - Susumu Mori
- From the Departments of Biomedical Engineering (F.P., D.A.H., N.A.T., E.R.M.), Medicine (H.A.), and Radiology (S.M., E.R.M), Johns Hopkins University, Baltimore, MD; Radiology and Imaging Sciences, National Institutes of Health Clinical Center, Bethesda, MD (N.G, D.A.B.); and Departments of Bioengineering, Medicine, and Radiology, University of California, San Diego (E.R.M.)
| | - Neville Gai
- From the Departments of Biomedical Engineering (F.P., D.A.H., N.A.T., E.R.M.), Medicine (H.A.), and Radiology (S.M., E.R.M), Johns Hopkins University, Baltimore, MD; Radiology and Imaging Sciences, National Institutes of Health Clinical Center, Bethesda, MD (N.G, D.A.B.); and Departments of Bioengineering, Medicine, and Radiology, University of California, San Diego (E.R.M.)
| | - David A Bluemke
- From the Departments of Biomedical Engineering (F.P., D.A.H., N.A.T., E.R.M.), Medicine (H.A.), and Radiology (S.M., E.R.M), Johns Hopkins University, Baltimore, MD; Radiology and Imaging Sciences, National Institutes of Health Clinical Center, Bethesda, MD (N.G, D.A.B.); and Departments of Bioengineering, Medicine, and Radiology, University of California, San Diego (E.R.M.)
| | - Natalia A Trayanova
- From the Departments of Biomedical Engineering (F.P., D.A.H., N.A.T., E.R.M.), Medicine (H.A.), and Radiology (S.M., E.R.M), Johns Hopkins University, Baltimore, MD; Radiology and Imaging Sciences, National Institutes of Health Clinical Center, Bethesda, MD (N.G, D.A.B.); and Departments of Bioengineering, Medicine, and Radiology, University of California, San Diego (E.R.M.)
| | - Elliot R McVeigh
- From the Departments of Biomedical Engineering (F.P., D.A.H., N.A.T., E.R.M.), Medicine (H.A.), and Radiology (S.M., E.R.M), Johns Hopkins University, Baltimore, MD; Radiology and Imaging Sciences, National Institutes of Health Clinical Center, Bethesda, MD (N.G, D.A.B.); and Departments of Bioengineering, Medicine, and Radiology, University of California, San Diego (E.R.M.).
| |
Collapse
|
45
|
Li C, Lim B, Hwang M, Song JS, Lee YS, Joung B, Pak HN. The Spatiotemporal Stability of Dominant Frequency Sites in In-Silico Modeling of 3-Dimensional Left Atrial Mapping of Atrial Fibrillation. PLoS One 2016; 11:e0160017. [PMID: 27459377 PMCID: PMC4961424 DOI: 10.1371/journal.pone.0160017] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2016] [Accepted: 07/12/2016] [Indexed: 12/19/2022] Open
Abstract
BACKGROUND We previously reported that stable rotors were observed in in-silico human atrial fibrillation (AF) models, and were well represented by dominant frequency (DF). We explored the spatiotemporal stability of DF sites in 3D-AF models imported from patient CT images of the left atrium (LA). METHODS We integrated 3-D CT images of the LA obtained from ten patients with persistent AF (male 80%, 61.8 ± 13.5 years old) into an in-silico AF model. After induction, we obtained 6 seconds of AF simulation data for DF analyses in 30 second intervals (T1-T9). The LA was divided into ten sections. Spatiotemporal changes and variations in the temporal consistency of DF were evaluated at each section of the LA. The high DF area was defined as the area with the highest 10% DF. RESULTS 1. There was no spatial consistency in the high DF distribution at each LA section during T1-T9 except in one patient (p = 0.027). 2. Coefficients of variation for the high DF area were highly different among the ten LA sections (p < 0.001), and they were significantly higher in the four pulmonary vein (PV) areas, the LA appendage, and the peri-mitral area than in the other LA sections (p < 0.001). 3. When we conducted virtual ablation of 10%, 15%, and 20% of the highest DF areas (n = 270 cases), AF was changed to atrial tachycardia (AT) or terminated at a rate of 40%, 57%, and 76%, respectively. CONCLUSIONS Spatiotemporal consistency of the DF area was observed in 10% of AF patients, and high DF areas were temporally variable. Virtual ablation of DF is moderately effective in AF termination and AF changing into AT.
Collapse
Affiliation(s)
- Changyong Li
- Yonsei University Health System, Seoul, Republic of Korea
| | - Byounghyun Lim
- Yonsei University Health System, Seoul, Republic of Korea
| | - Minki Hwang
- Yonsei University Health System, Seoul, Republic of Korea
| | - Jun-Seop Song
- Yonsei University Health System, Seoul, Republic of Korea
| | - Young-Seon Lee
- Yonsei University Health System, Seoul, Republic of Korea
| | - Boyoung Joung
- Yonsei University Health System, Seoul, Republic of Korea
| | - Hui-Nam Pak
- Yonsei University Health System, Seoul, Republic of Korea
| |
Collapse
|
46
|
Bayer JD, Roney CH, Pashaei A, Jaïs P, Vigmond EJ. Novel Radiofrequency Ablation Strategies for Terminating Atrial Fibrillation in the Left Atrium: A Simulation Study. Front Physiol 2016; 7:108. [PMID: 27148061 PMCID: PMC4828663 DOI: 10.3389/fphys.2016.00108] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2015] [Accepted: 03/07/2016] [Indexed: 12/19/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jason D Bayer
- Electrophysiology and Heart Modeling Institute (LIRYC), Bordeaux University FoundationPessac, France; Cardiothoracic Research Center of Bordeaux (Inserm U 1045), University of BordeauxBordeaux, France
| | - Caroline H Roney
- Electrophysiology and Heart Modeling Institute (LIRYC), Bordeaux University FoundationPessac, France; Institute of Mathematics of Bordeaux (IMB), University of BordeauxTalence, France
| | - Ali Pashaei
- Electrophysiology and Heart Modeling Institute (LIRYC), Bordeaux University FoundationPessac, France; Institute of Mathematics of Bordeaux (IMB), University of BordeauxTalence, France
| | - Pierre Jaïs
- Electrophysiology and Heart Modeling Institute (LIRYC), Bordeaux University FoundationPessac, France; Cardiothoracic Research Center of Bordeaux (Inserm U 1045), University of BordeauxBordeaux, France; Haut-Lévêque Cardiology Hospital, University Hospital Center (CHU) of BordeauxPessac, France
| | - Edward J Vigmond
- Electrophysiology and Heart Modeling Institute (LIRYC), Bordeaux University FoundationPessac, France; Institute of Mathematics of Bordeaux (IMB), University of BordeauxTalence, France
| |
Collapse
|
47
|
Hwang M, Song JS, Lee YS, Li C, Shim EB, Pak HN. Electrophysiological Rotor Ablation in In-Silico Modeling of Atrial Fibrillation: Comparisons with Dominant Frequency, Shannon Entropy, and Phase Singularity. PLoS One 2016; 11:e0149695. [PMID: 26909492 PMCID: PMC4766081 DOI: 10.1371/journal.pone.0149695] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2015] [Accepted: 02/02/2016] [Indexed: 11/18/2022] Open
Abstract
BACKGROUND Although rotors have been considered among the drivers of atrial fibrillation (AF), the rotor definition is inconsistent. We evaluated the nature of rotors in 2D and 3D in- silico models of persistent AF (PeAF) by analyzing phase singularity (PS), dominant frequency (DF), Shannon entropy (ShEn), and complex fractionated atrial electrogram cycle length (CFAE-CL) and their ablation. METHODS Mother rotor was spatiotemporally defined as stationary reentries with a meandering tip remaining within half the wavelength and lasting longer than 5 s. We generated 2D- and 3D-maps of the PS, DF, ShEn, and CFAE-CL during AF. The spatial correlations and ablation outcomes targeting each parameter were analyzed. RESULTS 1. In the 2D PeAF model, we observed a mother rotor that matched relatively well with DF (>9 Hz, 71.0%, p<0.001), ShEn (upper 2.5%, 33.2%, p<0.001), and CFAE-CL (lower 2.5%, 23.7%, p<0.001). 2. The 3D-PeAF model also showed mother rotors that had spatial correlations with DF (>5.5 Hz, 39.7%, p<0.001), ShEn (upper 8.5%, 15.1%, p <0.001), and CFAE (lower 8.5%, 8.0%, p = 0.002). 3. In both the 2D and 3D models, virtual ablation targeting the upper 5% of the DF terminated AF within 20 s, but not the ablations based on long-lasting PS, high ShEn area, or lower CFAE-CL area. CONCLUSION Mother rotors were observed in both 2D and 3D human AF models. Rotor locations were well represented by DF, and their virtual ablation altered wave dynamics and terminated AF.
Collapse
Affiliation(s)
- Minki Hwang
- Division of Cardiology, Yonsei University Health System, Seoul, Republic of Korea
| | - Jun-Seop Song
- Division of Cardiology, Yonsei University Health System, Seoul, Republic of Korea
| | - Young-Seon Lee
- Division of Cardiology, Yonsei University Health System, Seoul, Republic of Korea
| | - Changyong Li
- Division of Cardiology, Yonsei University Health System, Seoul, Republic of Korea
| | - Eun Bo Shim
- Department of Mechanical and Biomedical Engineering, Kangwon National University, Chuncheon, Kangwon-do, Republic of Korea
- * E-mail: (HNP); (EBS)
| | - Hui-Nam Pak
- Division of Cardiology, Yonsei University Health System, Seoul, Republic of Korea
- * E-mail: (HNP); (EBS)
| |
Collapse
|
48
|
Hwang M, Lee HS, Pak HN, Shim EB. Inducibility of human atrial fibrillation in an in silico model reflecting local acetylcholine distribution and concentration. THE KOREAN JOURNAL OF PHYSIOLOGY & PHARMACOLOGY : OFFICIAL JOURNAL OF THE KOREAN PHYSIOLOGICAL SOCIETY AND THE KOREAN SOCIETY OF PHARMACOLOGY 2015; 20:111-7. [PMID: 26807030 PMCID: PMC4722184 DOI: 10.4196/kjpp.2016.20.1.111] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2015] [Revised: 11/06/2015] [Accepted: 11/10/2015] [Indexed: 11/17/2022]
Abstract
Vagal nerve activity has been known to play a crucial role in the induction and maintenance of atrial fibrillation (AF). However, it is unclear how the distribution and concentration of local acetylcholine (ACh) promotes AF. In this study, we investigated the effect of the spatial distribution and concentration of ACh on fibrillation patterns in an in silico human atrial model. A human atrial action potential model with an ACh-dependent K+ current (IKAch) was used to examine the effect of vagal activation. A simulation of cardiac wave dynamics was performed in a realistic 3D model of the atrium. A model of the ganglionated plexus (GP) and nerve was developed based on the "octopus hypothesis". The pattern of cardiac wave dynamics was examined by applying vagal activation to the GP areas or randomly. AF inducibility in the octopus hypothesis-based GP and nerve model was tested. The effect of the ACh concentration level was also examined. In the single cell simulation, an increase in the ACh concentration shortened APD90 and increased the maximal slope of the restitution curve. In the 3D simulation, a random distribution of vagal activation promoted wavebreaks while ACh secretion limited to the GP areas did not induce a noticeable change in wave dynamics. The octopus hypothesis-based model of the GP and nerve exhibited AF inducibility at higher ACh concentrations. In conclusion, a 3D in silico model of the GP and parasympathetic nerve based on the octopus model exhibited higher AF inducibility with higher ACh concentrations.
Collapse
Affiliation(s)
- Minki Hwang
- Division of Cardiology, Yonsei University Health System, Seoul 03722, Korea
| | - Hyun-Seung Lee
- Department of Mechanical and Biomedical Engineering, Kangwon National University, Chuncheon 24341, Korea
| | - Hui-Nam Pak
- Division of Cardiology, Yonsei University Health System, Seoul 03722, Korea
| | - Eun Bo Shim
- Department of Mechanical and Biomedical Engineering, Kangwon National University, Chuncheon 24341, Korea
| |
Collapse
|
49
|
Pashaei A, Bayer J, Meillet V, Dubois R, Vigmond E. Computation and projection of spiral wave trajectories during atrial fibrillation: a computational study. Card Electrophysiol Clin 2015; 7:37-47. [PMID: 25784021 DOI: 10.1016/j.ccep.2014.11.001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Abstract
To show how atrial fibrillation rotor activity on the heart surface manifests as phase on the torso, fibrillation was induced on a geometrically accurate computer model of the human atria. The Hilbert transform, time embedding, and filament detection were compared. Electrical activity on the epicardium was used to compute potentials on different surfaces from the atria to the torso. The Hilbert transform produces erroneous phase when pacing for longer than the action potential duration. The number of phase singularities, frequency content, and the dominant frequency decreased with distance from the heart, except for the convex hull.
Collapse
Affiliation(s)
- Ali Pashaei
- LIRYC Electrophysiology and Heart Modelling Institute, University of Bordeaux, PTIB-Campus Xavier, Arnozan, Avenue du Haut Lévèque, Bordeaux 33600, France; Inserm U1045, Cardiothoracic Research Center, 146 rue Léo-Saignat, Bordeaux Cedex 33076, France.
| | - Jason Bayer
- LIRYC Electrophysiology and Heart Modelling Institute, University of Bordeaux, PTIB-Campus Xavier, Arnozan, Avenue du Haut Lévèque, Bordeaux 33600, France; Inserm U1045, Cardiothoracic Research Center, 146 rue Léo-Saignat, Bordeaux Cedex 33076, France
| | - Valentin Meillet
- LIRYC Electrophysiology and Heart Modelling Institute, University of Bordeaux, PTIB-Campus Xavier, Arnozan, Avenue du Haut Lévèque, Bordeaux 33600, France; Inserm U1045, Cardiothoracic Research Center, 146 rue Léo-Saignat, Bordeaux Cedex 33076, France
| | - Rémi Dubois
- LIRYC Electrophysiology and Heart Modelling Institute, University of Bordeaux, PTIB-Campus Xavier, Arnozan, Avenue du Haut Lévèque, Bordeaux 33600, France; Inserm U1045, Cardiothoracic Research Center, 146 rue Léo-Saignat, Bordeaux Cedex 33076, France
| | - Edward Vigmond
- LIRYC Electrophysiology and Heart Modelling Institute, University of Bordeaux, PTIB-Campus Xavier, Arnozan, Avenue du Haut Lévèque, Bordeaux 33600, France; Bordeaux Institute of Mathematics UMR 5251, University of Bordeaux, 351 cours de la Libération, Talence 33405, France
| |
Collapse
|
50
|
A two layers monodomain model of cardiac electrophysiology of the atria. J Math Biol 2015; 71:1607-41. [PMID: 25773466 DOI: 10.1007/s00285-015-0861-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2014] [Revised: 10/12/2014] [Indexed: 10/23/2022]
Abstract
Numerical simulations of the cardiac electrophysiology in the atria are often based on the standard bidomain or monodomain equations stated on a two-dimensional manifold. These simulations take advantage of the thinness of the atrial tissue, and their computational cost is reduced, as compared to three-dimensional simulations. However, these models do not take into account the heterogeneities located in the thickness of the tissue, like discontinuities of the fiber direction, although they can be a substrate for atrial arrhythmia (Hocini et al., Circulation 105(20):2442-2448, 2002; Ho et al., Cardiovasc Res 54(2):325-336, 2002; Nattel, Nature 415(6868):219-226, 2002). We investigate a two-dimensional model with two coupled, superimposed layers that allows to introduce three-dimensional heterogeneities, but retains a reasonable computational cost. We introduce the mathematical derivation of this model and error estimates with respect to the three-dimensional model. We give some numerical illustrations of its interest: we numerically show its convergence for vanishing thickness, introduce an optimization process of the coupling coefficient and assess its validity on physiologically relevant geometries. Our model would be an efficient tool to test the influence of three-dimensional fiber direction heterogeneities in reentries or atrial arrhythmia without using three-dimensional models.
Collapse
|