Patient-Specific Identification of Atrial Flutter Vulnerability–A Computational Approach to Reveal Latent Reentry Pathways

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

Personalizing Simulations of the Human Atria : Intracardiac Measurements, Tissue Conductivities, and Cellular Electrophysiology

This work addresses major challenges of heart model personalization. Analysis techniques for clinical intracardiac electrograms determine wave direction and conduction velocity from single beats. Electrophysiological measurements are simulated to validate the models. Uncertainties in tissue conductivities impact on simulated ECGs. A minimal model of cardiac myocytes is adapted to the atria. This makes personalized cardiac models a promising technique to improve treatment of atrial arrhythmias

Effects of Electrical and Structural Remodeling on Atrial Fibrillation Maintenance: A Simulation Study

Atrial fibrillation, a common cardiac arrhythmia, often progresses unfavourably: in patients with long-term atrial fibrillation, fibrillatory episodes are typically of increased duration and frequency of occurrence relative to healthy controls. This is due to electrical, structural, and contractile remodeling processes. We investigated mechanisms of how electrical and structural remodeling contribute to perpetuation of simulated atrial fibrillation, using a mathematical model of the human atrial action potential incorporated into an anatomically realistic three-dimensional structural model of the human atria. Electrical and structural remodeling both shortened the atrial wavelength - electrical remodeling primarily through a decrease in action potential duration, while structural remodeling primarily slowed conduction. The decrease in wavelength correlates with an increase in the average duration of atrial fibrillation/flutter episodes. The dependence of reentry duration on wavelength was the same for electrical vs. structural remodeling. However, the dynamics during atrial reentry varied between electrical, structural, and combined electrical and structural remodeling in several ways, including: (i) with structural remodeling there were more occurrences of fragmented wavefronts and hence more filaments than during electrical remodeling; (ii) dominant waves anchored around different anatomical obstacles in electrical vs. structural remodeling; (iii) dominant waves were often not anchored in combined electrical and structural remodeling. We conclude that, in simulated atrial fibrillation, the wavelength dependence of reentry duration is similar for electrical and structural remodeling, despite major differences in overall dynamics, including maximal number of filaments, wave fragmentation, restitution properties, and whether dominant waves are anchored to anatomical obstacles or spiralling freely

Simplified Cardiodynamic Tissue Electrophysiology Characterization, Reduced Order Modeling with Therapeutic Perspective

Atrial fibrillation (Afib) is the most common cardiac arrhythmia affecting millions of people around the world. Mapping and analysis of electrical activation patterns such as electric rotors during Afib is crucial in understanding arrhythmic mechanisms and assessment of diagnostic measures. To this end, there exists various mapping studies where textit{'quantitative'} features such as local activation time, dominant frequency, wave direction, and conduction velocity are extracted from recorded intracardiac electrograms (EGMs). However, obtaining quantitative features further adds to multiplicity of the data and henceforth does not help interpretation of measured signals as opposed to using a more compressed diagnostic terms such as linking the measurements to reentry mechanisms. Through some techniques it is possible to construct isopotential and phase mappings by the help of monophasic action potential recordings in higher spatial resolution. In those cases, however, both expensive mapping tools performing multi-site simultaneous recordings which are not available to most of electrophysiologists are required. On the other hand, the most commonly used catheters which provide high resolution but local measurements remain rather rudimentary in mapping a spatially more global arrhythmic behaviors in a simultaneous fashion. Spiral waves are tissue level phenomena observed in both clinical and experimental settings. They are the product of electrical rotors which are associated with reentry mechanisms during Afib. They can be reproduced using computer models of cardiac electrical activity. Current computer models vary in complexity, accuracy, and efficiency. One particular type is called biophysical models which are based on detailed ion channel interactions. Besides being computationally demanding, they are exceedingly complex and intractable preventing their use in a systems approach where multilevel events are generally considered together. Phenomenological models, on the other hand, include summarized details of ionic events yet preserve fundamental biophysical accuracy. A particular one of them, a minimal resistor model (MRM), was shown to reproduce relevant basic electrophysiological behaviors such as (action potential) AP and electrical restitution properties for human ventricular tissue. The objective in present thesis is to 'qualitatively' characterize fibrillatory wavefront propagation dynamics in cardiac tissue using simulated intracardiac EGMs obtained from most commonly used and lower cost catheter types providing high resolution but localized readings. Another purpose connected to the previous is to show adequacy of a phenomenological model, MRM, in reproducing biophysically related behaviors for human atria. In this respect, two category of problems are handled throughout the thesis: (1) parameter estimation of MRM and (2) discrimination of spiral wave behaviors through intracardiac EGMs simulated using MRM. In the first part, representativeness of MRM for human atrial electrophysiology is established through adaptation of it to a biophysically detailed model originated from experimental data. Specifically, a method is proposed for parameter estimation of the simple model, MRM, to match a targeted behavior such as AP and electrical restitutions first generated from a complex model, by using extended Kalman filter (EKF). In the second part, a method that receives intracardiac EGMs and returns corresponding wavefront propagation patterns classified in terms of electric rotor dynamics is introduced. The method incorporates an information theoretical distance which is called normalized compression distance (NCD) used for assessment of distance measure between simulated behaviors. Achieving outstanding performance together with robustness in discrimination through usage of simulated data enables a theoretical validation of the method. Proposed frameworks collectively yield (1) potential usability of a computationally efficient and easier in analysis model for tissue level cardiac events and (2) simplicity and practicality in clinics through a mapping from a multiple, complex EGM signals to electric rotor behaviors, symptoms more relevant to the diagnosis.Ph.D., Electrical Engineering -- Drexel University, 201

The contact electrogram and its architectural determinants in atrial fibrillation

The electrogram is the sine qua non of excitable tissues, yet classification in atrial fibrillation (AF) remains poorly related to substrate factors. The objective of this thesis was to establish the relationship between electrograms and two commonly implicated substrate factors, connexin 43 and fibrosis in AF. The substrates and methods chosen to achieve this ranged from human acutely induced AF using open chest surgical mapping (Chapter 6), ex vivo whole heart Langendorff (Chapter 7) with in vivo telemetry confirming spontaneous AF in a new species of rat, the Brown Norway and finally isolated atrial preparations from an older cohort of rats using orthogonal pacing and novel co-localisation methods at sub-millimetre resolution and in some atria, optical mapping (Chapter 8). In rodents, electrode size and spacing was varied (Chapters 5, 10) to study its effects on structure function correlations (Chapter 9). Novel indices of AF organisation and automated electrogram morphology were used to quantify function (Chapter 4). Key results include the discoveries that humans without any history of prior AF have sinus rhythm electrograms with high spectral frequency content, that wavefront propagation velocities correlated with fibrosis and connexin phosphorylation ratios, that AF heterogeneity of conduction correlates to fibrosis and that orthogonal pacing in heavily fibrosed atria causes anisotropy in electrogram-fibrosis correlations. Furthermore, fibrosis and connexin 43 have differing and distinct spatial resolutions in their relationship with AF organisational indices. In conclusion a new model of AF has been found, and structure function correlations shown on an unprecedented scale, but with caveats of electrode size and direction dependence. These findings impact structure function methods and prove the effect of substrate on AF organisation.Open Acces

Mechanistic Inquiry into the Role of Tissue Remodeling in Fibrotic Lesions in Human Atrial Fibrillation

AbstractAtrial fibrillation (AF), the most common arrhythmia in humans, is initiated when triggered activity from the pulmonary veins propagates into atrial tissue and degrades into reentrant activity. Although experimental and clinical findings show a correlation between atrial fibrosis and AF, the causal relationship between the two remains elusive. This study used an array of 3D computational models with different representations of fibrosis based on a patient-specific atrial geometry with accurate fibrotic distribution to determine the mechanisms by which fibrosis underlies the degradation of a pulmonary vein ectopic beat into AF. Fibrotic lesions in models were represented with combinations of: gap junction remodeling; collagen deposition; and myofibroblast proliferation with electrotonic or paracrine effects on neighboring myocytes. The study found that the occurrence of gap junction remodeling and the subsequent conduction slowing in the fibrotic lesions was a necessary but not sufficient condition for AF development, whereas myofibroblast proliferation and the subsequent electrophysiological effect on neighboring myocytes within the fibrotic lesions was the sufficient condition necessary for reentry formation. Collagen did not alter the arrhythmogenic outcome resulting from the other fibrosis components. Reentrant circuits formed throughout the noncontiguous fibrotic lesions, without anchoring to a specific fibrotic lesion

Modeling Human Atrial Patho-Electrophysiology from Ion Channels to ECG - Substrates, Pharmacology, Vulnerability, and P-Waves

Half of the patients suffering from atrial fibrillation (AF) cannot be treated adequately, today. This thesis presents multi-scale computational methods to advance our understanding of patho-mechanisms, to improve the diagnosis of patients harboring an arrhythmogenic substrate, and to tailor therapy. The modeling pipeline ranges from ion channels on the subcellular level up to the ECG on the body surface. The tailored therapeutic approaches carry the potential to reduce the burden of AF

Characterizing Atrial Fibrillation Substrate by Electrogram and Restitution Analysis

Vorhofflimmern ist die häufigste supraventrikuläre Arrhythmie in der klinischen Praxis. Es gibt Hinweise darauf, dass pathologisches Vorhofsubstrat (Fibrose) eine zentrale mechanistische Rolle bei der Aufrechterhaltung von Vorhofflimmern spielt. Die Behandlung von Vorhofflimmern erfolgt durch Ablation des fibrotischen Substrats. Der Nachweis eines solchen Substrats ist jedoch eine ungelöste Herausforderung, was durch die mangelnden positiven klinischen Ablationsergebnisse ersichtlich wird. Daher ist das Hauptthema dieser Arbeit die Charakterisierung des atrialen Substrats. Die Bestimmung von Signalmerkmalen an Stellen mit fibrotischem Substrat erleichtert die Erkennung und anschließende Ablation solcher Areale in Zukunft. Darüber hinaus kann das Verständnis der Art und Weise, wie diese Areale das Vorhofflimmern aufrechterhalten, die positiven Ergebnisse von Ablationseingriffen verbessern. Schließlich kann Restitutionsinformation ein weiteres Instrument zur Substratcharakterisierung sein, das bei der Unterscheidung zwischen pathologischen und nicht-pathologischen Arealen helfen und somit das Ablationsergebnis weiter verbessert. In dieser Arbeit werden zwei Ansätze zur Substratcharakterisierung vorgestellt: Zunächst wurde eine Charakterisierung des Substrats mit Hilfe des intraatrialen Elektrogramms vorgenommen. Dazu wurde eine Auswahl spezifischer Merkmale des Elektrogramms an Positionen evaluiert, die eine Terminierung von Vorhofflimmern nach Ablation zur Folge hatten. Die Studie beinhaltete 21 Patienten, bei denen eine Ablation nach Pulmonalvenenisolation das klinisch persistierende Vorhofflimmern beendete. Der klinisch vorgeschlagene Grenzwert der Spannungsamplitude von <0:5 mV wurde genutzt, um die Positionen der Ablation zu definieren. Die Bereiche, in denen das Vorhofflimmern erfolgreich terminiert wurde, wiesen ausgeprägte Elektrogramm-Muster auf. Diese waren gekennzeichnet durch kurze lokale Zykluslängen, die fraktionierte Potentiale und Niederspannungspotentiale enthielten. Gleichzeitig zeigten sie eine lokale Konsistenz und deckten einen Großteil der lokalen Vorhofflimmer-Zykluslänge ab. Die meisten dieser Bereiche wiesen auch im Sinusrhythmus pathologisch verzögerte atriale Spätpotentiale und fraktionierte Elektrogramme auf. Im zweiten Teil der Arbeit wurden Restitutionsdaten der lokalen Amplitude und der lokalen Leitungsgeschwindigkeit (CV) erfasst und genutzt, um daraus Informationen über das zugrunde liegende Substrat abzuleiten. Die Daten zur Restitution wurden von 22 Patienten mit Vorhofflimmern aus zwei Kliniken unter Verwendung eines S1S2-Protokolls mit Stimulationsintervallen von 180 ms bis 500 ms gewonnen. Um Restitutionsdaten der Patientengruppe zu erhalten, musste ein automatisierter Algorithmus entwickelt werden, der in der Lage ist, große Mengen an Stimulationsprotokolldaten zu lesen, zu segmentieren und auszuwerten. Dieser Algorithmus wurde in der vorliegenden Arbeit entwickelt und CVAR-Seg genannt. Der CVAR-Seg Algorithmus bietet eine rauschresistente Signalsegmentierung, die mit extremen Rauschpegeln getestet wurde, die weit über dem erwarteten klinischen Pegel lagen. CVAR-Seg wurde unter einer Open-Source-Lizenz für die Allgemeinheit bereitgestellt. Es ermöglicht aufgrund seines modularen Aufbaus den einfachen Austausch einzelner Verfahrensschritte durch alternative Methoden entsprechend den Bedürfnissen des Anwenders. Darüber hinaus wurde im Rahmen dieser Studie eine neuartige Methode, die sogenannte inverse Doppelellipsenmethode, zur Bestimmung der lokalen CV etabliert. Diese Methode schätzt die CV, die Faserorientierung und den Anisotropiefaktor bei beliebiger Elektrodenanordnung. In Simulationen reproduzierte die Doppelellipsenmethode die vorherrschende CV, Faserorientierung und Anisotropie genauer und robuster als die aktuell gängigste Methode. Zusätzlich erwies sich diese Methode als echtzeittauglich und könnte daher in klinischen Elektrophysiologiesystemen eingesetzt werden. Die Doppelellipsenmethode würde durch die lokalisierte Vermessung des Vorhofsubstrats ermöglichen während eines Kartierungsverfahrens gleichzeitig eine CV-Karte, eine Anisotropieverhältniskarte und eine Faserkarte zu erstellen. Die Restitutionsinformationen der Patientenkohorte wurden mit der CVARSeg-Pipeline und der inversen Doppelellipsenmethode ausgewertet, um Amplituden- und CV-Restitutionskurven zu erhalten. Zur Anpassung der Restitutionskurven wurde eine monoexponentielle Funktion verwendet. Die Parameter der angepassten Funktion, die die Restitutionskurven abbilden, wurden verwendet, um Unterschiede in den Restitutionseigenschaften zwischen pathologischem und nicht-pathologischem Substrat zu erkennen. Das Ergebnis zeigte, dass klinisch definierte pathologische Bereiche durch eine reduzierte Amplitudenasymptote und einen steilen Kurvenabfall bei erhöhter Stimulationsrate gekennzeichnet waren. CV-Kurven zeigten eine reduzierte Asymptote und eine große Variation im Parameter der den Kurvenabfall beschreibt. Darüber hinaus wurden die Restitutionsunterschiede innerhalb des Vorhofs an der posterioren und anterioren Wand verglichen, da die Literatur keine eindeutigen Ergebnisse lieferte. In dieser Arbeit wurde nachgewiesen, dass die posteriore Vorhofwand Amplituden- und CV-Restitutionskurven mit höherer Asymptote und moderaterer Krümmung verglichen mit der anterioren Vorhofwand aufweist. Um über den empirisch beschriebenen manuellen Schwellenwert hinauszugehen, wurde der Parameterraum, der von den Anpassungsparametern der Amplituden- und CV-Restitutionskurven aufgespannt wird, nach natürlich vorkommenden Clustern durchsucht. Obgleich Cluster vorhanden waren, deutete ihre unzureichende Trennung auf einen kontinuierlichen, sich mit dem Schweregrad der Substratpathologie verändernden Verlauf der Amplituden- und CV-Kurven hin. Schließlich wurde eine einfachere und schnellere Methode zur Erfassung von Restitutionsdaten vorgestellt, die einen vergleichbaren Informationsgehalt auf der Grundlage der maximalen Steigung anstelle einer vollständigen Restitutionskurve liefert. In dieser Arbeit werden zwei neue Methoden vorgestellt, der CVAR-Seg-Algorithmus und die inverse Doppelellipsenmethode, die eine Auswertung von S1S2 Stimulationsprotokollen und die Bestimmung der lokalen Leitungsgeschwindigkeit beschleunigen und verbessern. Darüber hinaus werden in dieser Arbeit Merkmale von pathologischem Gewebe definiert, die zur Identifizierung von Arrhythmiequellen beitragen. Somit trägt diese Arbeit dazu bei, die Therapie von Vorhofflimmern in Zukunft zu verbessern

Statistical and Graph-Based Signal Processing: Fundamental Results and Application to Cardiac Electrophysiology

The goal of cardiac electrophysiology is to obtain information about the mechanism, function, and performance of the electrical activities of the heart, the identification of deviation from normal pattern and the design of treatments. Offering a better insight into cardiac arrhythmias comprehension and management, signal processing can help the physician to enhance the treatment strategies, in particular in case of atrial fibrillation (AF), a very common atrial arrhythmia which is associated to significant morbidities, such as increased risk of mortality, heart failure, and thromboembolic events. Catheter ablation of AF is a therapeutic technique which uses radiofrequency energy to destroy atrial tissue involved in the arrhythmia sustenance, typically aiming at the electrical disconnection of the of the pulmonary veins triggers. However, recurrence rate is still very high, showing that the very complex and heterogeneous nature of AF still represents a challenging problem. Leveraging the tools of non-stationary and statistical signal processing, the first part of our work has a twofold focus: firstly, we compare the performance of two different ablation technologies, based on contact force sensing or remote magnetic controlled, using signal-based criteria as surrogates for lesion assessment. Furthermore, we investigate the role of ablation parameters in lesion formation using the late-gadolinium enhanced magnetic resonance imaging. Secondly, we hypothesized that in human atria the frequency content of the bipolar signal is directly related to the local conduction velocity (CV), a key parameter characterizing the substrate abnormality and influencing atrial arrhythmias. Comparing the degree of spectral compression among signals recorded at different points of the endocardial surface in response to decreasing pacing rate, our experimental data demonstrate a significant correlation between CV and the corresponding spectral centroids. However, complex spatio-temporal propagation pattern characterizing AF spurred the need for new signals acquisition and processing methods. Multi-electrode catheters allow whole-chamber panoramic mapping of electrical activity but produce an amount of data which need to be preprocessed and analyzed to provide clinically relevant support to the physician. Graph signal processing has shown its potential on a variety of applications involving high-dimensional data on irregular domains and complex network. Nevertheless, though state-of-the-art graph-based methods have been successful for many tasks, so far they predominantly ignore the time-dimension of data. To address this shortcoming, in the second part of this dissertation, we put forth a Time-Vertex Signal Processing Framework, as a particular case of the multi-dimensional graph signal processing. Linking together the time-domain signal processing techniques with the tools of GSP, the Time-Vertex Signal Processing facilitates the analysis of graph structured data which also evolve in time. We motivate our framework leveraging the notion of partial differential equations on graphs. We introduce joint operators, such as time-vertex localization and we present a novel approach to significantly improve the accuracy of fast joint filtering. We also illustrate how to build time-vertex dictionaries, providing conditions for efficient invertibility and examples of constructions. The experimental results on a variety of datasets suggest that the proposed tools can bring significant benefits in various signal processing and learning tasks involving time-series on graphs. We close the gap between the two parts illustrating the application of graph and time-vertex signal processing to the challenging case of multi-channels intracardiac signals

An iPS-derived in vitro model of human atrial conduction

Atrial fibrillation (AF) is the most common arrhythmia in the United States, affecting approximately 1 in 10 adults, and its prevalence is expected to rise as the population ages. Treatment options for AF are limited; moreover, the development of new treatments is hindered by limited (1) knowledge regarding human atrial electrophysiological endpoints (e.g., conduction velocity [CV]) and (2) accurate experimental models. Here, we measured the CV and refractory period, and subsequently calculated the conduction wavelength, in vivo (four subjects with AF and four controls), and ex vivo (atrial slices from human hearts). Then, we created an in vitro model of human atrial conduction using induced pluripotent stem (iPS) cells. This model consisted of iPS-derived human atrial cardiomyocytes plated onto a micropatterned linear 1D spiral design of Matrigel. The CV (34-41 cm/s) of the in vitro model was nearly five times faster than 2D controls (7-9 cm/s) and similar to in vivo (40-64 cm/s) and ex vivo (28-51 cm/s) measurements. Our iPS-derived in vitro model recapitulates key features of in vivo atrial conduction and may be a useful methodology to enhance our understanding of AF and model patient-specific disease