Bayesian Inference with Combined Dynamic and Sparsity Models: Application in 3D Electrophysiological Imaging

Data-driven inference is widely encountered in various scientific domains to convert the observed measurements into information that cannot be directly observed about a system. Despite the quickly-developing sensor and imaging technologies, in many domains, data collection remains an expensive endeavor due to financial and physical constraints. To overcome the limits in data and to reduce the demand on expensive data collection, it is important to incorporate prior information in order to place the data-driven inference in a domain-relevant context and to improve its accuracy. Two sources of assumptions have been used successfully in many inverse problem applications. One is the temporal dynamics of the system (dynamic structure). The other is the low-dimensional structure of a system (sparsity structure). In existing work, these two structures have often been explored separately, while in most high-dimensional dynamic system they are commonly co-existing and contain complementary information. In this work, our main focus is to build a robustness inference framework to combine dynamic and sparsity constraints. The driving application in this work is a biomedical inverse problem of electrophysiological (EP) imaging, which noninvasively and quantitatively reconstruct transmural action potentials from body-surface voltage data with the goal to improve cardiac disease prevention, diagnosis, and treatment. The general framework can be extended to a variety of applications that deal with the inference of high-dimensional dynamic systems

MRI-guided non-invasive epicardial mapping in patients with implanted pacing devices

Developing patient-specific 3D heart models to non-invasively localize arrhythmic foci. My research focussed on the application of MRI and complex electrocardiography in patients with MRI conditional pacemaker systems

When Cardiac Biophysics Meets Groupwise Statistics: Complementary Modelling Approaches for Patient-Specific Medicine

This habilitation manuscript contains research on biophysical and statistical modeling of the heart, as well as interactions between these two approaches

When Cardiac Biophysics Meets Groupwise Statistics: Complementary Modelling Approaches for Patient-Specific Medicine

This habilitation manuscript contains research on biophysical and statistical modeling of the heart, as well as interactions between these two approaches

Improving the Access of Cardiac Magnetic Resonance in Low-Middle Income Countries to Improve Cardiac Care: Rapid CMR

Non-communicable diseases– cancer and cardiovascular disease - are the leading causes of death in high-income countries (HICs). Cardiovascular disease, however, is increasing in Low-Middle Income Countries (LMICs) and is the emergent primary cause of mortality. Part of the reason is suboptimal therapies– from primary prevention to more advanced tertiary care. Not only are advanced therapies scarce but advanced diagnostic tests which apply to them are not fully available, and so diagnoses could be inaccurate and treatments poorly targeted. Within the portfolio and hierarchy of cardiovascular diagnostic testing, Cardiac Magnetic Resonance (CMR) is a crucial diagnostic imaging test that redefines diagnosis and enables targeted therapies, but is expensive with inadequate training and poor availability in LMICs countries. I demonstrated that CMR could be made fast, easy, and cheap – sufficient for delivery in five LMICs countries in three continents. To achieve this, I developed an abbreviated CMR protocol, focused on the core of CMR - volumes, function, and scar imaging (with selected additions like iron quantification), and by embedding the technical quality protocol within clinical care, training, and mentoring, so it proved to have diagnostic utility and change management, as well being a self-sustaining and essential service. I also used CMR as a research method in LMICs specifically to complement research in areas of a specific need to those countries, exploiting opportunities that were previously unavailable, with one chapter dedicated to evaluating early cardiovascular involvement in treated and non-treated people living with HIV in Peru, and a second chapter of the potential utility of CMR for screening cardiotoxicity and its comparison in precision with other cardiac imaging modalities in the UK, potentially extending the role of rapid CMR in HICs. Unlike traditional PhDs in medicine, my research involved technology adaptation, transfer, and collaboration. The project was multi-layered with political, social, educational, training, and partnership aspects, along with more traditional aspects such as clinical effectiveness and cost-effectiveness analysis. I showed the use of advanced cardiac imaging in LMICs by breaking down barriers, demonstrating that Rapid CMR can be possible in new clinical environments where much need exists

Computational modelling of the human heart and multiscale simulation of its electrophysiological activity aimed at the treatment of cardiac arrhythmias related to ischaemia and Infarction

[ES] Las enfermedades cardiovasculares constituyen la principal causa de morbilidad y mortalidad a nivel mundial, causando en torno a 18 millones de muertes cada año. De entre ellas, la más común es la enfermedad isquémica cardíaca, habitualmente denominada como infarto de miocardio (IM). Tras superar un IM, un considerable número de pacientes desarrollan taquicardias ventriculares (TV) potencialmente mortales durante la fase crónica del IM, es decir, semanas, meses o incluso años después la fase aguda inicial. Este tipo concreto de TV normalmente se origina por una reentrada a través de canales de conducción (CC), filamentos de miocardio superviviente que atraviesan la cicatriz del infarto fibrosa y no conductora. Cuando los fármacos anti-arrítmicos resultan incapaces de evitar episodios recurrentes de TV, la ablación por radiofrecuencia (ARF), un procedimiento mínimamente invasivo realizado mediante cateterismo en el laboratorio de electrofisiología (EF), se usa habitualmente para interrumpir de manera permanente la propagación eléctrica a través de los CCs responsables de la TV. Sin embargo, además de ser invasivo, arriesgado y requerir mucho tiempo, en casos de TVs relacionadas con IM crónico, hasta un 50% de los pacientes continúa padeciendo episodios recurrentes de TV tras el procedimiento de ARF. Por tanto, existe la necesidad de desarrollar nuevas estrategias pre-procedimiento para mejorar la planificación de la ARF y, de ese modo, aumentar esta tasa de éxito relativamente baja. En primer lugar, realizamos una revisión exhaustiva de la literatura referente a los modelos cardiacos 3D existentes, con el fin de obtener un profundo conocimiento de sus principales características y los métodos usados en su construcción, con especial atención sobre los modelos orientados a simulación de EF cardíaca. Luego, usando datos clínicos de un paciente con historial de TV relacionada con infarto, diseñamos e implementamos una serie de estrategias y metodologías para (1) generar modelos computacionales 3D específicos de paciente de ventrículos infartados que puedan usarse para realizar simulaciones de EF cardíaca a nivel de órgano, incluyendo la cicatriz del infarto y la región circundante conocida como zona de borde (ZB); (2) construir modelos 3D de torso que permitan la obtención del ECG simulado; y (3) llevar a cabo estudios in-silico de EF personalizados y pre-procedimiento, tratando de replicar los verdaderos estudios de EF realizados en el laboratorio de EF antes de la ablación. La finalidad de estas metodologías es la de localizar los CCs en el modelo ventricular 3D para ayudar a definir los objetivos de ablación óptimos para el procedimiento de ARF. Por último, realizamos el estudio retrospectivo por simulación de un caso, en el que logramos inducir la TV reentrante relacionada con el infarto usando diferentes configuraciones de modelado para la ZB. Validamos nuestros resultados mediante la reproducción, con una precisión razonable, del ECG del paciente en TV, así como en ritmo sinusal a partir de los mapas de activación endocárdica obtenidos invasivamente mediante sistemas de mapeado electroanatómico en este último caso. Esto permitió encontrar la ubicación y analizar las características del CC responsable de la TV clínica. Cabe destacar que dicho estudio in-silico de EF podría haberse efectuado antes del procedimiento de ARF, puesto que nuestro planteamiento está completamente basado en datos clínicos no invasivos adquiridos antes de la intervención real. Estos resultados confirman la viabilidad de la realización de estudios in-silico de EF personalizados y pre-procedimiento de utilidad, así como el potencial del abordaje propuesto para llegar a ser en un futuro una herramienta de apoyo para la planificación de la ARF en casos de TVs reentrantes relacionadas con infarto. No obstante, la metodología propuesta requiere de notables mejoras y validación por medio de es[CA] Les malalties cardiovasculars constitueixen la principal causa de morbiditat i mortalitat a nivell mundial, causant entorn a 18 milions de morts cada any. De elles, la més comuna és la malaltia isquèmica cardíaca, habitualment denominada infart de miocardi (IM). Després de superar un IM, un considerable nombre de pacients desenvolupen taquicàrdies ventriculars (TV) potencialment mortals durant la fase crònica de l'IM, és a dir, setmanes, mesos i fins i tot anys després de la fase aguda inicial. Aquest tipus concret de TV normalment s'origina per una reentrada a través dels canals de conducció (CC), filaments de miocardi supervivent que travessen la cicatriu de l'infart fibrosa i no conductora. Quan els fàrmacs anti-arítmics resulten incapaços d'evitar episodis recurrents de TV, l'ablació per radiofreqüència (ARF), un procediment mínimament invasiu realitzat mitjançant cateterisme en el laboratori de electrofisiologia (EF), s'usa habitualment per a interrompre de manera permanent la propagació elèctrica a través dels CCs responsables de la TV. No obstant això, a més de ser invasiu, arriscat i requerir molt de temps, en casos de TVs relacionades amb IM crònic fins a un 50% dels pacients continua patint episodis recurrents de TV després del procediment d'ARF. Per tant, existeix la necessitat de desenvolupar noves estratègies pre-procediment per a millorar la planificació de l'ARF i, d'aquesta manera, augmentar la taxa d'èxit, que es relativament baixa. En primer lloc, realitzem una revisió exhaustiva de la literatura referent als models cardíacs 3D existents, amb la finalitat d'obtindre un profund coneixement de les seues principals característiques i els mètodes usats en la seua construcció, amb especial atenció sobre els models orientats a simulació de EF cardíaca. Posteriorment, usant dades clíniques d'un pacient amb historial de TV relacionada amb infart, dissenyem i implementem una sèrie d'estratègies i metodologies per a (1) generar models computacionals 3D específics de pacient de ventricles infartats capaços de realitzar simulacions de EF cardíaca a nivell d'òrgan, incloent la cicatriu de l'infart i la regió circumdant coneguda com a zona de vora (ZV); (2) construir models 3D de tors que permeten l'obtenció del ECG simulat; i (3) dur a terme estudis in-silico de EF personalitzats i pre-procediment, tractant de replicar els vertaders estudis de EF realitzats en el laboratori de EF abans de l'ablació. La finalitat d'aquestes metodologies és la de localitzar els CCs en el model ventricular 3D per a ajudar a definir els objectius d'ablació òptims per al procediment d'ARF. Finalment, a manera de prova de concepte, realitzem l'estudi retrospectiu per simulació d'un cas, en el qual aconseguim induir la TV reentrant relacionada amb l'infart usant diferents configuracions de modelatge per a la ZV. Validem els nostres resultats mitjançant la reproducció, amb una precisió raonable, del ECG del pacient en TV, així com en ritme sinusal a partir dels mapes d'activació endocardíac obtinguts invasivament mitjançant sistemes de mapatge electro-anatòmic en aquest últim cas. Això va permetre trobar la ubicació i analitzar les característiques del CC responsable de la TV clínica. Cal destacar que aquest estudi in-silico de EF podria haver-se efectuat abans del procediment d'ARF, ja que el nostre plantejament està completament basat en dades clíniques no invasius adquirits abans de la intervenció real. Aquests resultats confirmen la viabilitat de la realització d'estudis in-silico de EF personalitzats i pre-procediment d'utilitat, així com el potencial de l'abordatge proposat per a arribar a ser en un futur una eina de suport per a la planificació de l'ARF en casos de TVs reentrants relacionades amb infart. No obstant això, la metodologia proposada requereix de notables millores i validació per mitjà d'estudis de simulació amb grans cohorts de pacients.[EN] Cardiovascular diseases represent the main cause of morbidity and mortality worldwide, causing around 18 million deaths every year. Among these diseases, the most common one is the ischaemic heart disease, usually referred to as myocardial infarction (MI). After surviving to a MI, a considerable number of patients develop life-threatening ventricular tachycardias (VT) during the chronic stage of the MI, that is, weeks, months or even years after the initial acute phase. This particular type of VT is typically sustained by reentry through slow conducting channels (CC), which are filaments of surviving myocardium that cross the non-conducting fibrotic infarct scar. When anti-arrhythmic drugs are unable to prevent recurrent VT episodes, radiofrequency ablation (RFA), a minimally invasive procedure performed by catheterization in the electrophysiology (EP) laboratory, is commonly used to interrupt the electrical conduction through the CCs responsible for the VT permanently. However, besides being invasive, risky and time-consuming, in the cases of VTs related to chronic MI, up to 50% of patients continue suffering from recurrent VT episodes after the RFA procedure. Therefore, there exists a need to develop novel pre-procedural strategies to improve RFA planning and, thereby, increase this relatively low success rate. First, we conducted an exhaustive review of the literature associated with the existing 3D cardiac models in order to gain a deep knowledge about their main features and the methods used for their construction, with special focus on those models oriented to simulation of cardiac EP. Later, using a clinical dataset of a chronically infarcted patient with a history of infarct-related VT, we designed and implemented a number of strategies and methodologies to (1) build patient-specific 3D computational models of infarcted ventricles that can be used to perform simulations of cardiac EP at the organ level, including the infarct scar and the surrounding region known as border zone (BZ); (2) construct 3D torso models that enable to compute the simulated ECG; and (3) carry out pre-procedural personalized in-silico EP studies, trying to replicate the actual EP studies conducted in the EP laboratory prior to the ablation. The goal of these methodologies is to allow locating the CCs into the 3D ventricular model in order to help in defining the optimal ablation targets for the RFA procedure. Lastly, as a proof-of-concept, we performed a retrospective simulation case study, in which we were able to induce an infarct-related reentrant VT using different modelling configurations for the BZ. We validated our results by reproducing with a reasonable accuracy the patient's ECG during VT, as well as in sinus rhythm from the endocardial activation maps invasively recorded via electroanatomical mapping systems in this latter case. This allowed us to find the location and analyse the features of the CC responsible for the clinical VT. Importantly, such in-silico EP study might have been conducted prior to the RFA procedure, since our approach is completely based on non-invasive clinical data acquired before the real intervention. These results confirm the feasibility of performing useful pre-procedural personalized in-silico EP studies, as well as the potential of the proposed approach to become a helpful tool for RFA planning in cases of infarct-related reentrant VTs in the future. Nevertheless, the developed methodology requires further improvements and validation by means of simulation studies including large cohorts of patients.During the carrying out of this doctoral thesis, the author Alejandro Daniel López Pérez was financially supported by the Ministerio de Economía, Industria y Competitividad of Spain through the program Ayudas para contratos predoctorales para la formación de doctores, with the grant number BES-2013-064089.López Pérez, AD. (2019). Computational modelling of the human heart and multiscale simulation of its electrophysiological activity aimed at the treatment of cardiac arrhythmias related to ischaemia and Infarction [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/124973TESI

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

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 book 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

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 book 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

Modified mass-spring system for physically based deformation modeling

Mass-spring systems are considered the simplest and most intuitive of all deformable models. They are computationally efficient, and can handle large deformations with ease. But they suffer several intrinsic limitations. In this book a modified mass-spring system for physically based deformation modeling that addresses the limitations and solves them elegantly is presented. Several implementations in modeling breast mechanics, heart mechanics and for elastic images registration are presented