A new algorithm to diagnose atrial ectopic origin from multi lead ECG systems - insights from 3D virtual human atria and torso

Rapid atrial arrhythmias such as atrial fibrillation (AF) predispose to ventricular arrhythmias, sudden cardiac death and stroke. Identifying the origin of atrial ectopic activity from the electrocardiogram (ECG) can help to diagnose the early onset of AF in a cost-effective manner. The complex and rapid atrial electrical activity during AF makes it difficult to obtain detailed information on atrial activation using the standard 12-lead ECG alone. Compared to conventional 12-lead ECG, more detailed ECG lead configurations may provide further information about spatio-temporal dynamics of the body surface potential (BSP) during atrial excitation. We apply a recently developed 3D human atrial model to simulate electrical activity during normal sinus rhythm and ectopic pacing. The atrial model is placed into a newly developed torso model which considers the presence of the lungs, liver and spinal cord. A boundary element method is used to compute the BSP resulting from atrial excitation. Elements of the torso mesh corresponding to the locations of the placement of the electrodes in the standard 12-lead and a more detailed 64-lead ECG configuration were selected. The ectopic focal activity was simulated at various origins across all the different regions of the atria. Simulated BSP maps during normal atrial excitation (i.e. sinoatrial node excitation) were compared to those observed experimentally (obtained from the 64-lead ECG system), showing a strong agreement between the evolution in time of the simulated and experimental data in the P-wave morphology of the ECG and dipole evolution. An algorithm to obtain the location of the stimulus from a 64-lead ECG system was developed. The algorithm presented had a success rate of 93%, meaning that it correctly identified the origin of atrial focus in 75/80 simulations, and involved a general approach relevant to any multi-lead ECG system. This represents a significant improvement over previously developed algorithms

Non-invasive localization of atrial ectopic beats by using simulated body surface P-wave integral maps

Non-invasive localization of continuous atrial ectopic beats remains a cornerstone for the treatment of atrial arrhythmias. The lack of accurate tools to guide electrophysiologists leads to an increase in the recurrence rate of ablation procedures. Existing approaches are based on the analysis of the P-waves main characteristics and the forward body surface potential maps (BSPMs) or on the inverse estimation of the electric activity of the heart from those BSPMs. These methods have not provided an efficient and systematic tool to localize ectopic triggers. In this work, we propose the use of machine learning techniques to spatially cluster and classify ectopic atrial foci into clearly differentiated atrial regions by using the body surface P-wave integral map (BSPiM) as a biomarker. Our simulated results show that ectopic foci with similar BSPiM naturally cluster into differentiated non-intersected atrial regions and that new patterns could be correctly classified with an accuracy of 97% when considering 2 clusters and 96% for 4 clusters. Our results also suggest that an increase in the number of clusters is feasible at the cost of decreasing accuracy.This work was partially supported by The "Programa Prometeu" from Conselleria d'Educacio Formacio I Ocupacio, Generalitat Valenciana (www.edu.gva.es/fio/index_es.asp) Award Number: PROMETEU/2016/088 to JS; The "Plan Estatal de Investigacion Cientifica y Tecnica y de Innovacion 2013-2016" from the Ministerio de Economia, Industria y Competitividad of Spain, Agencia Estatal de Investigacion (www.mineco.gob.es) and the European Commission (European Regional Development Funds - ERDF -FEDER) (ec.europa.eu/regional_policy/es/funding/erdf/) Award Number: DPI2016-75799-R to JS and The "Programa Estatal de Investigacion, Desarrollo e Innovacion Orientado a los Retos de la Sociedad" from the Ministerio de Economia y Competitividad of Spain, Agencia Estatal de Investigacion (www.mineco.gob.es) and the European Commission (European Regional Development Funds - ERDF -FEDER) (ec.europa.eu/regional_policy/es/funding/erdf/) Award Number: TIN2014-59932-JIN to AFA and RS. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Ferrer Albero, A.; Godoy, EJ.; Lozano, M.; Martínez Mateu, L.; Alonso Atienza, F.; Saiz Rodríguez, FJ.; Sebastián Aguilar, R. (2017). Non-invasive localization of atrial ectopic beats by using simulated body surface P-wave integral maps. PLoS ONE. 12(7):1-23. https://doi.org/10.1371/journal.pone.0181263S12312

Novel non-invasive algorithm to identify the origins of re-entry and ectopic foci in the atria from 64-lead ECGs: A computational study.

Atrial tachy-arrhytmias, such as atrial fibrillation (AF), are characterised by irregular electrical activity in the atria, generally associated with erratic excitation underlain by re-entrant scroll waves, fibrillatory conduction of multiple wavelets or rapid focal activity. Epidemiological studies have shown an increase in AF prevalence in the developed world associated with an ageing society, highlighting the need for effective treatment options. Catheter ablation therapy, commonly used in the treatment of AF, requires spatial information on atrial electrical excitation. The standard 12-lead electrocardiogram (ECG) provides a method for non-invasive identification of the presence of arrhythmia, due to irregularity in the ECG signal associated with atrial activation compared to sinus rhythm, but has limitations in providing specific spatial information. There is therefore a pressing need to develop novel methods to identify and locate the origin of arrhythmic excitation. Invasive methods provide direct information on atrial activity, but may induce clinical complications. Non-invasive methods avoid such complications, but their development presents a greater challenge due to the non-direct nature of monitoring. Algorithms based on the ECG signals in multiple leads (e.g. a 64-lead vest) may provide a viable approach. In this study, we used a biophysically detailed model of the human atria and torso to investigate the correlation between the morphology of the ECG signals from a 64-lead vest and the location of the origin of rapid atrial excitation arising from rapid focal activity and/or re-entrant scroll waves. A focus-location algorithm was then constructed from this correlation. The algorithm had success rates of 93% and 76% for correctly identifying the origin of focal and re-entrant excitation with a spatial resolution of 40 mm, respectively. The general approach allows its application to any multi-lead ECG system. This represents a significant extension to our previously developed algorithms to predict the AF origins in association with focal activities

Three-dimensional Multiscale Modelling and Simulation of Atria and Torso Electrophysiology

A better understanding of the electrical activity of the heart under physiological and pathological conditions has always been key for clinicians and researchers. Over the last years, the information in the P-wave signals has been extensively analysed to un-cover the mechanisms underlying atrial arrhythmias by localizing ectopic foci or high-frequency rotors. However, the relationship between the activation of the different areas of the atria and the characteristics of the P-wave signals or body surface poten-tial maps are still far from being completely understood. Multiscale anatomical and functional models of the heart are a new technological framework that can enable the investigation of the heart as a complex system. This thesis is centred in the construction of a multiscale framework that allows the realistic simulation of atrial and torso electrophysiology and integrates all the anatom-ical and functional descriptions described in the literature. The construction of such model involves the development of heterogeneous cellular and tissue electrophysiolo-gy models fitted to empirical data. It also requires an accurate 3D representation of the atrial anatomy, including tissue fibre arrangement, and preferential conduction axes. This multiscale model aims to reproduce faithfully the activation of the atria under physiological and pathological conditions. We use the model for two main applica-tions. First, to study the relationship between atrial activation and surface signals in sinus rhythm. This study should reveal the best places for recording P-waves signals in the torso, and which are the regions of the atria that make the most significant contri-bution to the body surface potential maps and determine the main P-wave characteris-tics. Second, to spatially cluster and classify ectopic atrial foci into clearly differenti-ated atrial regions by using the body surface P-wave integral map (BSPiM) as a bi-omarker. We develop a machine-learning pipeline trained from simulations obtained from the atria-torso model aiming to validate whether ectopic foci with similar BSPiM naturally cluster into differentiated non-intersected atrial regions, and whether new BSPiM could be correctly classified with high accuracy.En la actualidad, una mejor compresión de la actividad eléctrica del corazón en condi-ciones fisiológicas y patológicas es clave para médicos e investigadores. A lo largo de los últimos años, la información derivada de la onda P se ha utilizado para intentar descubrir los mecanismos subyacentes a las arritmias auriculares mediante la localiza-ción de focos ectópicos y rotores de alta frecuencia. Sin embargo, la relación entre la activación de distintas regiones auriculares y las características tanto de las ondas P como de la distribución de potencial en la superficie del torso está lejos de entenderse completamente. Los modelos cardíacos funcionales y anatómicos son una nueva he-rramienta que puede facilitar la investigación relativa al corazón entendido como sis-tema complejo. La presente tesis se centra en la construcción de un modelo multiescala para la simula-ción realista de la electrofisiología cardíaca tanto a nivel auricular como de torso, integrando toda la información anatómica y funcional disponible en la literatura. La construcción de este modelo implica el desarrollo, en base a datos experimentales, de modelos electrofisiológicos heterogéneos tanto celulares como tisulares. Así mismo, es imprescindible una representación tridimensional precisa de la anatomía auricular, incluyendo la dirección de fibras y los haces de conducción preferentes. Este modelo multiescala busca reproducir fielmente la activación auricular en condiciones fisiológi-cas y patológicas. Su uso se ha centrado fundamentalmente en dos aplicaciones. En primer lugar, estudiar la relación entre la activación auricular en ritmo sinusal y las señales en la superficie del torso. Este estudio busca definir la mejor ubicación para el registro de las ondas P en el torso así como determinar aquellas regiones auriculares que contribuyen fundamentalmente a la formación y distribución de potenciales super-ficiales así como a las características de las ondas P. En segundo lugar, agrupar y cla-sificar espacialmente los focos ectópicos en regiones auriculares claramente diferen-ciables empleando como biomarcador los mapas superficiales de integral de la onda P (BSPiM). Se ha desarrollado para ello una metodología de aprendizaje automático en la que las simulaciones obtenidas con el modelo multiescala aurícula-torso sirven de entrenamiento, permitiendo validar si los focos ectópicos cuyos BSPiMs son similares se agrupan de forma natural en regiones auriculares no intersectadas y si BSPiMs nue-vos podrían ser clasificados prospectivamente con gran precisión.Avui en dia, una millor comprenssió de l'activitat elèctrica del cor en condicions fisio-lògiques i patològiques és clau per a metges i investigadors. Al llarg dels últims anys, la informació derivada de l'ona P s'ha utilitzat per intentar descobrir els mecanismes subjacents a les arítmies auriculars mitjançant la localització de focus ectòpics i rotors d'alta freqüència. No obstant això, la relació entre l'activació de diferents regions auri-culars i les característiques tant de les ones P com de la distribució de potencial en la superfície del tors està lluny d'entendre's completament. Els models cardíacs funcionals i anatòmics són una nova eina que pot facilitar la recerca relativa al cor entès com a sistema complex. La present tesi es centra en la construcció d'un model multiescala per a la simulació realista de la electrofisiologia cardíaca tant a nivell auricular com de tors, integrant tota la informació anatòmica i funcional disponible en la literatura. La construcció d'aquest model implica el desenvolupament, sobre la base de dades experimentals, de models electrofisiològics heterogenis, tant cel·lulars com tissulars. Així mateix, és imprescindible una representació tridimensional precisa de l'anatomia auricular, in-cloent la direcció de fibres i els feixos de conducció preferents. Aquest model multies-cala busca reproduir fidelment l'activació auricular en condicions fisiològiques i pa-tològiques. El seu ús s'ha centrat fonamentalment en dues aplicacions. En primer lloc, estudiar la relació entre l'activació auricular en ritme sinusal i els senyals en la superfí-cie del tors. A més a més, amb aquest estudi també es busca definir la millor ubicació per al registre de les ones P en el tors, així com, determinar aquelles regions auriculars que contribueixen fonamentalment a la formació i distribució de potencials superfi-cials a l'hora que es caracteritzen les ones P. En segon lloc, agrupar i classificar espa-cialment els focus ectòpics en regions auriculars clarament diferenciables emprant com a biomarcador els mapes superficials d'integral de l'ona P (BSPiM). És per això que s'ha desenvolupat una metodologia d'aprenentatge automàtic en la qual les simulacions obtingudes amb el model multiescala aurícula-tors serveixen d'entrenament, la qual cosa permet validar si els focus ectòpics, llurs BSPiMs són similars, s'agrupen de for-ma natural en regions auriculars no intersectades i si BSPiMs nous podrien ser classifi-cats de manera prospectiva amb precisió.Ferrer Albero, A. (2017). Three-dimensional Multiscale Modelling and Simulation of Atria and Torso Electrophysiology [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/88402TESI

Contributions To The Methodology Of Electrocardiographic Imaging (ECGI) And Application Of ECGI To Study Mechanisms Of Atrial Arrhythmia, Post Myocardial Infarction Electrophysiological Substrate, And Ventricular Tachycardia In Patients

ABSTRACT OF THE DISSERTATION Contributions to the Methodology of Electrocardiographic Imaging: ECGI) and Application of ECGI to Study Mechanisms of Atrial Arrhythmia, Post Myocardial Infarction Electrophysiological Substrate, and Ventricular Tachycardia in Patients by Yong Wang Doctor of Philosophy in Biomedical Engineering Washington University in St. Louis, 2009 Professor Yoram Rudy, Chair Electrocardiographic Imaging: ECGI) is a noninvasive imaging modality for cardiac electrophysiology and arrhythmia. ECGI reconstructs epicardial potentials, electrograms and isochrones from body-surface electrocardiograms combined with heart-torso geometry from computed tomography: CT). The application of a new meshless method, the Method of Fundamental Solutions: MFS) is introduced to ECGI with the following major advantages: 1. Elimination of meshing and manual mesh optimization processes, thereby enhancing automation and speeding the ECGI procedure. 2. Elimination of mesh-induced artifacts. 3. Simpler implementation. These properties of MFS enhance the practical application of ECGI as a clinical diagnostic tool. The current ECGI mode of operation is offline with generation of epicardial potential maps delayed to data acquisition. A real time ECGI procedure is proposed, by which the epicardial potentials can be reconstructed while the body surface potential data are acquired: \u3c 1msec/frame) during a clinical procedure. This development enables real-time monitoring, diagnosis, and interactive guidance of intervention for arrhythmia therapy. ECGI is applied to map noninvasively the electrophysiological substrate in eight post-MI patients during sinus rhythm: SR). Contrast-enhanced MRI: ceMRI) is conducted to determine anatomical scar. ECGI imaged regions of electrical scar corresponded closely in location, extent, and morphology to the anatomical scars. In three patients, late diastolic potentials are imaged in the scar epicardial border zone during SR. Scar-related ventricular tachycardia: VT) in two patients are imaged, showing the VT activation sequence in relation to the abnormal electrophysiological substrate. ECGI imaging the substrate in a beat-by-beat fashion could potentially help in noninvasive risk stratification for post-MI arrhythmias and facilitate substrate-based catheter ablation of these arrhythmias. ECGI is applied to eleven consecutive patients referred for VT catheter ablation procedure. ECGI is performed either before: 8 patients) or during: 3 patients) the ablation procedure. Blinded ECGI and invasive electrophysiology: EP) study results are compared. Over a wide range of VT types and locations, ECGI results are consistent with EP data regarding localization of the arrhythmia origin: including myocardial depth) and mechanism: focal, reentrant, fascicular). ECGI also provides mechanistic electrophysiological insights, relating arrhythmia patterns to the myocardial substrate. The study shows ECGI has unique potential clinical advantages, especially for hemodynamically intolerant VT or VT that is difficult to induce. Because it provides local cardiac information, ECGI may aid in better understanding of mechanisms of ventricular arrhythmia. Further prospective trials of ECGI with clinical endpoints are warranted. Many mechanisms for the initiation and perpetuation of atrial fibrillation: AF) have been demonstrated over the last several decades. The tools to study these mechanisms in humans have limitations, the most common being invasiveness of a mapping procedure. In this paper, we present simultaneous noninvasive biatrial epicardial activation sequences of AF in humans, obtained using the Electrocardiographic Imaging: ECGI) system, and analyzed in terms of mechanisms and complexity of activation patterns. We performed ECGI in 36 patients with a diagnosis of AF. To determine ECGI atrial accuracy, atrial pacing from different sites was performed in six patients: 37 pacing events), and ECGI was compared to registered CARTO images. Then, ECGI was performed on all 36 patients during AF and ECGI epicardial maps were analyzed for mechanisms and complexity. ECGI noninvasively imaged the low-amplitude signals of AF in a wide range of patients: 97% procedural success). The spatial accuracy in determining initiation sites as simulated by atrial pacing was ~ 6mm. ECGI imaged many activation patterns of AF, most commonly multiple wavelets: 92%), with pulmonary vein: 69%) and non-pulmonary vein: 62%) trigger sites. Rotor activity was seen rarely: 15%). AF complexity increased with longer clinical history of AF, though the degree of complexity of nonparoxysmal AF varied and overlapped. ECGI offers a way to identify unique epicardial activation patterns of AF in a patient-specific manner. The results are consistent with contemporary animal models of AF mechanisms and highlight the coexistence of a variety of mechanisms among patients

Reconstruction of atrial ectopic focal and Re-entrant excitations from body surface potentials. Insights from 3D virtual human atria and torso

Non-invasive electrocardiographic imaging has been seen as a painless and economic method to map the electrical functions of the heart. However, it is still a great challenge to obtain accurate reconstruction of cardiac electrical activity from body surface potentials (BSP) due to the ill-posed behaviour of the cardiac inverse-problem. Though some advances have been made in solving the inverse-problem, few studies have been conducted for the atria, which have dramatic differences to the ventricles in their anatomical structures and electrophysiological properties. It is unclear either how the spatial resolution of electrodes on the BSP and rapid excitation rates of atrial activation during atrial fibrillation affect the accuracy of the inverse-problem. In this study, we used a biophysically detailed model of the human atria and torso to investigate effects of multi-lead ECG on the accuracy of reconstructed atrial excitation pattern on the epicardiac surface during the time courses of atrial fibrillation induced by electrical remodelling. It was shown that the solution of the atrial inverse-problem was dependent on the spatial resolution of electrodes on the body surface. The solution was also influenced by the morphology of the AP, rate and types of atrial excitation as well as the implantation of variant orders of the Tikhonov regularization method

Mapping of the electrical activity of human atria. Multiscale modelling and simulations

La fibrilación auricular es una de las arritmias cardíacas más comunes observadas en la práctica clínica. Por lo tanto, es de vital importancia desarrollar nuevas tecnologías destinadas a diagnosticar y acabar con este tipo de arritmia, para mejorar la calidad de vida de los pacientes y reducir los costes de los sistemas nacionales de salud. En los últimos años ha aumentado el uso de las nuevas técnicas de mapeo auricular, basadas en sistemas multi-electrodo para mapear la actividad eléctrica en humanos. Dichas técnicas permiten localizar y ablacionar los impulsores de la fibrilación auricular, como son las fuentes focales o los rotores. Sin embargo, todavía existe incertidumbre sobre su precisión y los procedimientos experimentales para su análisis están limitados debido a su carácter invasivo. Por lo tanto, las simulaciones computacionales son una herramienta muy útil para superar estas limitaciones, al permitir reproducir con fidelidad las observaciones experimentales, dividir el problema bajo estudio en sub-estudios más simples, y realizar investigaciones preliminares imposibles de llevar a cabo en el práctica clínica. Esta tesis doctoral se centra en el análisis de la precisión de los sistemas de mapeo multi-electrodo a través de modelos y simulaciones computacionales. Para ello, desarrollamos modelos realistas multi-escala con el objetivo de simular actividad eléctrica auricular reentrante, en primer lugar en una lámina de tejido auricular, y en segundo lugar en las aurículas completas. Posteriormente, analizamos los efectos de las configuraciones geométricas multi-electrodo en la precisión de la localización de los rotores, mediante el uso de agrupaciones multi-electrodo con distancias inter-electrodo equidistantes, así como a través de catéteres de tipo basket con distancias inter-electrodo no equidistantes. Después de calcular los electrogramas unipolares intracavitarios, realizamos mapas de fase, detecciones de singularidad de fase para rastrear los rotores, y mapas de frecuencia dominantes. Finalmente, descubrimos que la precisión de los sistemas de mapeo multi-electrodo depende de su posición dentro de la cavidad auricular, de la distancia entre los electrodos y el tejido, de la distancia inter-electrodo, y de la contribución de las fuentes de campo lejano. Además, como consecuencia de estos factores que pueden afectar a la precisión de los sistemas de mapeo multi-electrodo, observamos la aparición de rotores falsos que podrían contribuir al fracaso de los procesos de ablación de la fibrilación auricular.Atrial fibrillation is one of the most common cardiac arrhythmias seen in clinical practice. Therefore, it is of vital importance to develop new technologies aimed at diagnosing and terminating this kind of arrhythmia, to improve the quality of life of patients and to reduce costs to national health systems. In the last years, new atrial mapping techniques based on multi-electrode systems are increasingly being used to map the atrial electrical activity in humans and localise and target atrial fibrillation drivers in the form of focal sources or rotors. However, significant concerns remain about their accuracy and experimental approaches to analyse them are limited due to their invasive character. Therefore, computer simulations are a helpful tool to overcome these limitations since they can reproduce with fidelity experimental observations, permit to split the problem to treat into more simple substudies, and allow the possibility of performing preliminary investigations impossible to carry out in the clinical practice. This PhD thesis is focused on the analysis for accuracy of the multielectrode mapping systems through computational models and simulations. For this purpose, we developed realistic multiscale models in order to simulate atrial electrical reentrant activity, first in a sheet of atrial tissue and, then, in the whole atria. Then, we analysed the effects of the multi-electrode geometrical configurations on the accuracy of localizing rotors, by using multi-electrode arrays with equidistant inter-electrode distances, as well as multi-electrode basket catheters with non-equidistant inter-electrode distances. After computing the intracavitary unipolar electrograms, we performed phase maps, phase singularity detections to track rotors, and dominant frequency maps. We finally found out that the accuracy of multi-electrode mapping systems depends on their position inside the atrial cavity, the electrode-to-tissue distance, the inter-electrode distance, and the contribution of far field sources. Furthermore, as a consequence of these factors, false rotors might appear and could contribute to failure of atrial fibrillation ablation procedures.La fibril·lació auricular és una de les arítmies cardíaques més comuns observades en la pràctica clínica. Per tant, és de vital importància desenvolupar noves tecnologies destinades a diagnosticar i acabar amb aquest tipus d'arítmia, per tal de millorar la qualitat de vida dels pacients i reduir els costos dels sistemes nacionals de salut. En els últims anys, ha augmentat l'ús de les noves tècniques de mapeig auricular, basades en sistemes multielèctrode per a mapejar l'activitat elèctrica auricular en humans. Aquestes tècniques permeten localitzar i ablacionar els impulsors de la fibril·lació auricular, com són les fonts focals o els rotors. No obstant això, encara hi ha incertesa sobre la seua precisió i els procediments experimentals per al seu anàlisi estan limitats a causa del seu caràcter invasiu. Per tant, les simulacions computacionals són una eina molt útil per a superar aquestes limitacions, en permetre reproduir amb fidelitat les observacions experimentals, dividir el problema sota estudi en subestudis més simples, i realitzar investigacions preliminars impossibles de dur a terme en el pràctica clínica. Aquesta tesi doctoral es centra en l'anàlisi de la precisió del sistemes de mapeig multielèctrode mitjançant els models i les simulacions computacionals. Per a això, desenvolupàrem models realistes multiescala per tal de simular activitat elèctrica auricular reentrant, en primer lloc en una làmina de teixit auricular, i en segon lloc a les aurícules completes. Posteriorment, analitzàrem els efectes de les configuracions geomètriques multielèctrode en la precisió de la localització dels rotors, mitjançant l'ús d'agrupacions multielèctrode amb distàncies interelèctrode equidistants, així com catèters de tipus basket amb distàncies interelèctrode no equidistants. Després de calcular els electrogrames unipolars intracavitaris, vam realitzar mapes de fase, deteccions de singularitat de fase per a rastrejar els rotors, i mapes de freqüència dominants. Finalment, vam descobrir que la precisió dels sistemes de mapeig multielèctrode depèn de la seua posició dins de la cavitat auricular, de la distància entre els elèctrodes i el teixit, de la distància interelèctrode, i de la contribució de les fonts de camp llunyà. A més, com a conseqüència d'aquests factors, es va observar l'aparició de rotors falsos que podrien contribuir al fracàs de l'ablació de la fibril·lació auricular.Martínez Mateu, L. (2018). Mapping of the electrical activity of human atria. Multiscale modelling and simulations [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/104604TESI

Mathematical modeling approaches for the diagnosis and treatment of reentrant atrial tachyarrhythmias

[EN] Atrial tachyarrhythmias present a high prevalence in the developed world, and several studies predict that in the coming decades it will be increased. Micro or macro-reentrant mechanisms of the electrical wavefronts that govern the mechanical behavior of the heart are one of the main responsibles for the maintenance of these arrhythmias. Atrial flutter is maintained by a macro-reentry around an anatomical or functional obstacle located in the atria. In the case of atrial fibrillation, the hypothesis which describes high frequency rotors as dominant sources of the fibrillation and responsible for the maintenance of the arrhythmia, has been gaining relevance in the last years. However, the therapies that target high frequency sources have a limited efficacy with current techniques. Radiofrequency ablation allows the destruction of parts of the cardiac tissue resulting in the interruption of the reentrant circuit in case of macro-reentries or the isolation of micro-reentrant circuits. The non-invasive location of reentrant circuits would increment the efficacy of these therapies and would shorten surgery interventions. In parallel, pharmacological therapies modify ionic expressions associated to the excitability and electrical refractoriness of the cardiac tissue with the objective of hindering the maintenance of reentrant behaviors. These therapies require a deep knowledge of the ionic mechanisms underlying the reentrant behavior and its properties in order to be effective. The research in these mechanisms allows the evaluation of new targets for the treatment and thus may improve the efficacy in atrial fibrillation termination. In this thesis, mathematical modeling is used to go forward in the minimization of the limitations associated to these treatments. Body surface potential mapping has been evaluated, both clinically and by means of mathematical simulations for the diagnosis and location of macro-reentrant circuits. The analysis of phase maps obtained from multiple lead electrocardiographic recordings distributed in the whole torso allowed the discrimination between different reentrant circuits. It is the reason why this technique is presented as a tool for the non-invasive location of macro and micro-reentrant circuits. A population of mathematical models designed in this thesis based on the action potentials recordings of atrial cardiomyocites from 149 patients, allowed the evaluation of the ionic mechanisms defining the properties of reentrant behaviors. This study has allowed us defining the blockade of ICaL as a target for the pharmacological treatment. The blockade of this current is associated with the increase of the movement in the core of the rotor which easies the collision of the rotor with other wavefronts or anatomical obstacles promoting the extinction of the reentry. The variability observed between patients modeled in our population has allowed showing and explaining the mechanisms promoting divergent results of a single treatment. This is why the introduction of populations of models will allow the prevention of side effects associated to inter-subject variability and to go forward in the development of individualized therapies. These works are built through a simulation platform of cardiac electrophysiology based in Graphic Processing Units (GPUs) and developed in this thesis. The platform allows the simulation of cellular models, tissues and organs with a realistic geometry and shows features comparable to that of the platforms used by the most relevant electrophysiology research groups at the moment.[ES] Las taquiarritmias auriculares tienen una alta prevalencia en el mundo desarrollado, además diversos estudios poblacionales indican que en las próximas décadas ésta se verá incrementada. Los mecanismos de micro o macro-reentrada de los frentes de onda eléctricos que rigen el comportamiento mecánico del corazón, se presentan como una de las principales causas del mantenimiento de estas arritmias. El flutter auricular es mantenido por un macro-reentrada alrededor de un obstáculo anatómico o funcional en las aurículas, mientras que en el caso de la fibrilación auricular la hipótesis que define a los rotores de alta frecuencia como elementos dominantes y responsables del mantenimiento de la arritmia se ha ido imponiendo al resto en los últimos años. Sin embargo, las terapias que tienen como objetivo finalizar o aislar estas reentradas tienen todavía una eficacia limitada. La ablación por radiofrecuencia permite eliminar zonas del tejido cardiaco resultando en la interrupción del circuito de reentrada en el caso de macro-reentradas o el aislamiento de comportamientos micro-reentrantes. La localización no invasiva de los circuitos reentrantes incrementaría la eficacia de estas terapias y reduciría la duración de las intervenciones quirúrgicas. Por otro lado, las terapias farmacológicas alteran las expresiones iónicas asociadas a la excitabilidad y la refractoriedad del tejido con el fin de dificultar el mantenimiento de comportamientos reentrantes. Este tipo de terapias exigen incrementar el conocimiento de los mecanismos subyacentes que explican el proceso de reentrada y sus propiedades, la investigación de estos mecanismos permite definir las dianas terapéuticas que mejoran la eficacia en la extinción de estos comportamientos. En esta tesis el modelado matemático se utiliza para dar un paso importante en la minimización de las limitaciones asociadas a estos tratamientos. La cartografía eléctrica de superficie ha sido testada, clínicamente y con simulaciones matemática,s como técnica de diagnóstico y localización de circuitos macro-reentrantes. El análisis de mapas de fase obtenidos a partir de los registros multicanal de derivaciones electrocardiográficas distribuidas en la superficie del torso permite diferenciar distintos circuitos de reentrada. Es por ello que esta técnica de registro y análisis se presenta como una herramienta para la localización no invasiva de circuitos macro y micro-reentrantes. Una población de modelos matemáticos, diseñada en esta tesis a partir de los registros de los potenciales de acción de 149 pacientes, ha permitido evaluar los mecanismos iónicos que definen las propiedades asociadas a los procesos de reentrada. Esto ha permitido apuntar al bloqueo de la corriente ICaL como diana terapéutica. Ésta se asocia al incremento del movimiento del núcleo que facilita el impacto del rotor con otros frentes de onda u obstáculos extinguiéndose así el comportamiento reentrante. La variabilidad entre pacientes reflejada en la población de modelos ha permitido además mostrar los mecanismos por los cuales un mismo tratamiento puede mostrar efectos divergentes, así el uso de poblaciones de modelos matemáticos permitirá prevenir efectos secundarios asociados a la variabilidad entre pacientes y profundizar en el desarrollo de terapias individualizadas. Estos trabajos se cimientan sobre una plataforma de simulación de electrofisiología cardiaca de basado en Unidades de Procesado Gráfico (GPUs) y desarrollada en esta tesis. La plataforma permite la simulación de modelos celulares cardiacos así como de tejidos u órganos con geometría realista, mostrando unas prestaciones comparables con las de las utilizadas por los grupos de investigación más potentes en el campo de la electrofisiología.[CA] Les taquiarítmies auriculars tenen una alta prevalença en el món desenvolupat, a més diversos estudis poblacionals indiquen que en les pròximes dècades aquesta es veurà incrementada. Els mecanismes de micro o macro-reentrada dels fronts d'ona elèctrics que regeixen el comportament mecànic del cor, es presenten com una de les principals causes del manteniment d'aquestes arítmies. El flutter auricular és mantingut per una macro-reentrada al voltant d'un obstacle anatòmic o funcional en les aurícules, mentre que en el cas de la fibril·lació auricular la hipòtesi que defineix als rotors d'alta freqüència com a elements dominants i responsables del manteniment de l'arítmia s'ha anat imposant a la resta en els últims anys. No obstant això, les teràpies que tenen com a objectiu finalitzar o aïllar aquestes reentrades tenen encara una eficàcia limitada. L'ablació per radiofreqüència permet eliminar zones del teixit cardíac resultant en la interrupció del circuit de reentrada en el cas de macro-reentrades o l'aïllament de comportaments micro-reentrants. La localització no invasiva dels circuits reentrants incrementaria l'eficàcia d'aquestes teràpies i reduiria la durada de les intervencions quirúrgiques. D'altra banda, les teràpies farmacològiques alteren les expressions iòniques associades a la excitabilitat i la refractaritat del teixit amb la finalitat de dificultar el manteniment de comportaments reentrants. Aquest tipus de teràpies exigeixen incrementar el coneixement dels mecanismes subjacents que expliquen el procés de reentrada i les seues propietats, la recerca d'aquests mecanismes permet definir les dianes terapèutiques que milloren l'eficàcia en l'extinció d'aquests comportaments. En aquesta tesi el modelatge matemàtic s'utilitza per a fer un pas important en la minimització de les limitacions associades a aquests tractaments. La cartografia elèctrica de superfície ha sigut testada, clínicament i amb simulacions matemàtiques com a tècnica de diagnòstic i localització de circuits macro-reentrants. L'anàlisi de mapes de fase obtinguts a partir dels registres multicanal de derivacions electrocardiogràfiques distribuïdes en la superfície del tors permet diferenciar diferents circuits de reentrada. És per açò que aquesta tècnica de registre i anàlisi es presenta com una eina per a la localització no invasiva de circuits macro i micro-reentrants. Una població de models matemàtics, dissenyada en aquesta tesi a partir dels registres dels potencials d'acció de 149 pacients, ha permès avaluar els mecanismes iònics que defineixen les propietats associades als processos de reentrada. Açò ha permès apuntar al bloqueig del corrent ICaL com a diana terapèutica. Aquesta s'associa a l'increment del moviment del nucli que facilita l'impacte del rotor amb altres fronts d'ona o obstacles extingint-se així el comportament reentrant. La variabilitat entre pacients reflectida en la població de models ha permès a més mostrar els mecanismes pels quals un mateix tractament pot mostrar efectes divergents, així l'ús de poblacions de models matemàtics permetrà prevenir efectes secundaris associats a la variabilitat entre pacients i aprofundir en el desenvolupament de teràpies individualitzades. Aquests treballs es fonamenten sobre una plataforma de simulació de electrofisiologia cardíaca basat en Unitats de Processament Gràfic (GPUs) i desenvolupada en aquesta tesi. La plataforma permet la simulació de models cel·lulars cardíacs així com de teixits o òrgans amb geometria realista, mostrant unes prestacions comparables amb les de les utilitzades per els grups de recerca més importants en aquesta área.Liberos Mascarell, A. (2016). Mathematical modeling approaches for the diagnosis and treatment of reentrant atrial tachyarrhythmias [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/62166TESI

A Computational Based Approach for Non-invasive Localization of Atrial ectopic foci

Las arritmias auriculares son las arritmias cardı́acas más comunes, afectan a seis millones de personas en Europa e imponen una enorme carga sanitaria en la sociedad. Las nuevas tecnologı́as médicas están ayudando a los electrofisiólogos a adaptar el tratamiento a cada paciente de diferentes maneras. Por ejemplo, la resonancia magnética (MRI) permite evaluar la distribución espacial de la fibrosis auricular; los mapas electroanatómicos (EAM) permiten obtener una caracterización eléctrica de los tejidos en tiempo real; Las imágenes electrocardiográficas (ECGI) permiten estudiar la actividad eléctrica cardı́aca de forma no invasiva; y la ablación por radiofrecuencia (RFA), permite eliminar el tejido patológico en el corazón que desencadena o mantiene una arritmia. A pesar del acceso a tecnologı́as avanzadas y de la existencia de guı́as clı́nicas bien desarrolladas para el tratamiento de las arritmias auriculares, las tasas de éxito del tratamiento a largo plazo siguen siendo bajas, debido a la complejidad de la enfermedad. Por lo tanto, existe una necesidad imperiosa de mejorar los resultados clı́nicos en beneficio de los pacientes y el sistema de salud. Se podrı́an emplear modelos biofı́sicos detallados de las aurı́culas y el torso para integrar todos los datos del paciente en un solo modelo 3D de referencia capaz de reproducir los complejos patrones de activación eléctrica observados en experimentos y la clı́nica. Sin embargo, existen algunas limitaciones relacionadas con la dificultad de construir tales modelos para cada paciente o realizar un número considerable de simulaciones para planificar la terapia óptima de RFA. Teniendo en cuenta todas esas limitaciones, proponemos utilizar modelos biofı́sicos detallados y simulaciones como una herramienta para entrenar sistemas de aprendizaje automático, para lo cual dispondrı́amos de todos los datos y variables del problema, que serı́an imposibles de obtener en un entorno clı́nico real. Por lo tanto, podemos realizar cientos de simulaciones electrofisiológicas, considerando una variedad de escenarios y patologı́as comunes, y entrenar un sistema que deberı́a ser capaz de reconocerlos a partir de un conjunto limitado de datos no invasivos del paciente, como un electrocardiograma (ECG), o mapa de potencial de superficie corporal (BSPM).Abstract Atrial arrhythmias are the most common cardiac arrhythmia, affecting six million people in Europe and imposing a huge healthcare bur- den on society. New technologies are helping electrophysiologists to tailor the treatment to each patient in different ways. For instance, magnetic resonance imaging (MRI) allows to assess the spatial distribution of atrial fibrosis; electro-anatomical maps (EAM) permit to obtain an electrical char- acterization of tissue in real-time; electrocardiographic imaging (ECGI) al- lows to study cardiac electrical activity non-invasively; and radiofrequency ablation (RFA), allows to eliminate pathological tissue in the heart that is triggering or sustaining an arrhythmia. Despite the access to advanced technologies and well-developed clinical guidelines for the management of atrial arrhythmia, long-term treatment success rates remain low, due to the complexity of the disease. Therefore, there is a compelling need to improve clinical outcomes for the benefit of patients and the healthcare system. Detailed biophysical models of the atria and torso could be employed to integrate all the patient data into a single reference 3D model able to re- produce the complex electrical activation patterns observed in experiments and clinics. However, there are some limitations related to the difficulty of building such models for each patient, or performing a substantial number of simulations to plan the optimal RFA therapy. Considering all those lim- itations, we propose to use detailed biophysical models and simulations as a tool to train machine learning systems, for which we have all the ground- truth data which would be impossible to obtain in a real clinical setting. Therefore, we can perform hundreds of electrophysiology simulations, con- sidering a variety of common scenarios and pathologies, and train a system that should be able to recognize them from a limited set of non-invasive pa- tient data, such as an electrocardiogram (ECG), or a body surface potential map (BSPM)