Personalized noninvasive imaging of volumetric cardiac electrophysiology

Three-dimensionally distributed electrical functioning is the trigger of mechanical contraction of the heart. Disturbance of this electrical flow is known to predispose to mechanical catastrophe but, due to its amenability to certain intervention techniques, a detailed understanding of subject-specific cardiac electrophysiological conditions is of great medical interest. In current clinical practice, body surface potential recording is the standard tool for diagnosing cardiac electrical dysfunctions. However, successful treatments normally require invasive catheter mapping for a more detailed observation of these dysfunctions. In this dissertation, we take a system approach to pursue personalized noninvasive imaging of volumetric cardiac electrophysiology. Under the guidance of existing scientific knowledge of the cardiac electrophysiological system, we extract the subject specific cardiac electrical information from noninvasive body surface potential mapping and tomographic imaging data of individual subjects. In this way, a priori knowledge of system physiology leads the physiologically meaningful interpretation of personal data; at the same time, subject-specific information contained in the data identifies parameters in individual systems that differ from prior knowledge. Based on this perspective, we develop a physiological model-constrained statistical framework for the quantitative reconstruction of the electrical dynamics and inherent electrophysiological property of each individual cardiac system. To accomplish this, we first develop a coupled meshfree-BE (boundary element) modeling approach to represent existing physiological knowledge of the cardiac electrophysiological system on personalized heart-torso structures. Through a state space system approach and sequential data assimilation techniques, we then develop statistical model-data coupling algorithms for quantitative reconstruction of volumetric transmembrane potential dynamics and tissue property of 3D myocardium from body surface potential recoding of individual subjects. We also introduce a data integration component to build personalized cardiac electrophysiology by fusing tomographic image and BSP sequence of the same subject. In addition, we develop a computational reduction strategy that improves the efficiency and stability of the framework. Phantom experiments and real-data human studies are performed for validating each of the framework’s major components. These experiments demonstrate the potential of our framework in providing quantitative understanding of volumetric cardiac electrophysiology for individual subjects and in identifying latent threats in individual’s heart. This may aid in personalized diagnose, treatment planning, and fundamentally, prevention of fatal cardiac arrhythmia

A meshless fragile points method for rule-based definition of myocardial fiber orientation

Background and objective: Rule-based methods are commonly used to estimate the arrangement of myocardial fibers by solving the Laplace problem with appropriate Dirichlet boundary conditions. Existing algorithms are using the Finite Element Method (FEM) to solve the Laplace–Dirichlet problem. However, meshless methods are under development for cardiac electrophysiology simulation. The objective of this work is to propose a meshless rule based method for the determination of myocardial fiber arrangement without requiring a mesh discretization as it is required by FEM. Methods: The proposed method employs the Fragile Points Method (FPM) for the solution of the Laplace–Dirichlet problem. FPM uses simple discontinuous trial functions and single-point exact integration for linear trial functions that set it as a promising alternative to the Finite Element Method. We derive the FPM formulation of the Laplace–Dirichlet and we estimate ventricular and atrial fiber arrangements according to rules based on histology findings for four different geometries. The obtained fiber arrangements from FPM are compared with the ones obtained from FEM by calculating the angle between the fiber vector fields of the two methods for three different directions (i.e., longitudinal, sheet, transverse). Results:The fiber arrangements that were generated with FPM were in close agreement with the generated arrangements from FEM for all three directions. The mean angle difference between the FPM and FEM vector fields were lower than for the ventricular fiber arrangements and lower than for the atrial fiber arrangements. Discussion:The proposed meshless rule-based method was proven to generate myocardial fiber arrangements with very close agreement with FEM while alleviates the requirement for a mesh of the latter. This is of great value for cardiac electrophysiology solvers that are based on meshless methods since they require a well defined myocardial fiber arrangement to simulate accurately the propagation of electrical signals in the heart. Combining a meshless solution for both the determination of the fibers and the electrical signal propagation can allow for solution that do not require the definition of a mesh. To our knowledge, this work is the first one to propose a meshless rule-based method for myocardial fiber arrangement determination

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

Meshfree and Particle Methods in Biomechanics: Prospects and Challenges

The use of meshfree and particle methods in the field of bioengineering and biomechanics has significantly increased. This may be attributed to their unique abilities to overcome most of the inherent limitations of mesh-based methods in dealing with problems involving large deformation and complex geometry that are common in bioengineering and computational biomechanics in particular. This review article is intended to identify, highlight and summarize research works on topics that are of substantial interest in the field of computational biomechanics in which meshfree or particle methods have been employed for analysis, simulation or/and modeling of biological systems such as soft matters, cells, biological soft and hard tissues and organs. We also anticipate that this review will serve as a useful resource and guide to researchers who intend to extend their work into these research areas. This review article includes 333 references

Simulation of action potential propagation based on the ghost structure method

In this paper, a ghost structure (GS) method is proposed to simulate the monodomain model in irregular computational domains using finite difference without regenerating body-fitted grids. In order to verify the validity of the GS method, it is first used to solve the Fitzhugh-Nagumo monodomain model in rectangular and circular regions at different states (the stationary and moving states). Then, the GS method is used to simulate the propagation of the action potential (AP) in transverse and longitudinal sections of a healthy human heart, and with left bundle branch block (LBBB). Finally, we analyze the AP and calcium concentration under healthy and LBBB conditions. Our numerical results show that the GS method can accurately simulate AP propagation with different computational domains either stationary or moving, and we also find that LBBB will cause the left ventricle to contract later than the right ventricle, which in turn affects synchronized contraction of the two ventricles

Analysis of age-related left ventricular collagen remodeling in living donors: Implications in arrhythmogenesis

Age-related fibrosis in the left ventricle (LV) has been mainly studied in animals by assessing collagen content. Using second-harmonic generation microscopy and image processing, we evaluated amount, aggregation and spatial distribution of LV collagen in young to old pigs, and middle-age and elder living donors. All collagen features increased when comparing adult and old pigs with young ones, but not when comparing adult with old pigs or middle-age with elder individuals. Remarkably, all collagen parameters strongly correlated with lipofuscin, a biological age marker, in humans. By building patient-specific models of human ventricular tissue electrophysiology, we confirmed that amount and organization of fibrosis modulated arrhythmia vulnerability, and that distribution should be accounted for arrhythmia risk assessment. In conclusion, we characterize the age-associated changes in LV collagen and its potential implications for ventricular arrhythmia development. Consistency between pig and human results substantiate the pig as a relevant model of age-related LV collagen dynamics. © 2022 The Author(s

Steady-state and transient effects of SK channel block and adrenergic stimulation to counteract acetylcholine-induced arrhythmogenic effects in the human atria: a computational study

Hyperactivity of the parasympathetic nervous system has been linked to the development of paroxysmal atrial fibrillation (AF). The parasympathetic neurotransmitter acetylcholine (ACh) causes a reduction in action potential (AP) duration (APD) and an increase in resting membrane potential (RMP), both of which contribute to enhance the risk for reentry. Research suggests that small-conductance calcium activated potassium (SK) channels may be an effective target for treating AF. Therapies targeting the autonomic nervous system, either alone or in combination with other drugs, have been explored and have been shown to decrease the incidence of atrial arrhythmias. This study uses computational modeling and simulation to examine the impact of SK channel block (SKb) and-adrenergic stimulation through Isoproterenol (Iso) on countering the negative effects of cholinergic activity in human atrial cell and 2D tissue models. The steady-state effects of Iso and/or SKb on AP shape, APD at 90% repolarization (APD90) and RMP were evaluated. The ability to terminate stable rotational activity in cholinergically-stimulated 2D tissue models of AF was also investigated. A range of SKb and Iso application kinetics, which reflect varying drug binding rates, were taken into consideration. The results showed that SKb alone prolonged APD90 and was able to stop sustained rotors in the presence of ACh concentrations up to 0.01 M. Iso terminated rotors under all tested ACh concentrations, but resulted in highly-variable steady-state outcomes depending on baseline AP morphology. Importantly, the combination of SKb and Iso resulted in greater APD90 prolongation and showed promising anti-arrhythmic potential by stopping stable rotors and preventing re-inducibility

Steady-state and transient effects of SK channel block and adrenergic stimulation to counteract acetylcholine-induced arrhythmogenic effects in the human atria: A computational study

Hyperactivity of the parasympathetic nervous system has been linked to the development of paroxysmal atrial fibrillation (AF). The parasympathetic neurotransmitter acetylcholine (ACh) causes a reduction in action potential (AP) duration (APD) and an increase in resting membrane potential (RMP), both of which contribute to enhance the risk for reentry. Research suggests that small-conductance calcium activated potassium (SK) channels may be an effective target for treating AF. Therapies targeting the autonomic nervous system, either alone or in combination with other drugs, have been explored and have been shown to decrease the incidence of atrial arrhythmias. This study uses computational modeling and simulation to examine the impact of SK channel block (SKb) and β-adrenergic stimulation through Isoproterenol (Iso) on countering the negative effects of cholinergic activity in human atrial cell and 2D tissue models. The steady-state effects of Iso and/or SKb on AP shape, APD at 90% repolarization (APD90) and RMP were evaluated. The ability to terminate stable rotational activity in cholinergically-stimulated 2D tissue models of AF was also investigated. A range of SKb and Iso application kinetics, which reflect varying drug binding rates, were taken into consideration. The results showed that SKb alone prolonged APD90 and was able to stop sustained rotors in the presence of ACh concentrations up to 0.01 μM. Iso terminated rotors under all tested ACh concentrations, but resulted in highly-variable steady-state outcomes depending on baseline AP morphology. Importantly, the combination of SKb and Iso resulted in greater APD90 prolongation and showed promising anti-arrhythmic potential by stopping stable rotors and preventing re-inducibility

Transferring Generalized Knowledge from Physics-based Simulation to Clinical Domain

A primary factor for the success of machine learning is the quality of labeled training data. However, in many fields, labeled data can be costly, difficult, or even impossible to acquire. In comparison, computer simulation data can now be generated at a much higher abundance with a much lower cost. These simulation data could potentially solve the problem of data deficiency in many machine learning tasks. Nevertheless, due to model assumptions, simplifications and possible errors, there is always a discrepancy between simulated and real data. This discrepancy needs to be addressed when transferring the knowledge from simulation to real data. Furthermore, simulation data is always tied to specific settings of models parameters, many of which have a considerable range of variations yet not necessarily relevant to the machine learning task of interest. The knowledge extracted from simulation data must thus be generalizable across these parameter variations before being transferred. In this dissertation, we address the two outlined challenges in leveraging simulation data to overcome the shortage of labeled real data, . We do so in a clinical task of localizing the origin of ventricular activation from 12 lead electrocardiograms (ECGs), where the clinical ECG data with labeled sites of origin in the heart can only be invasively available. By adopting the concept of domain adaptation, we address the discrepancy between simulated and clinical ECG data by learning the shift between the two domains using a large amount of simulation data and a small amount of clinical data. By adopting the concept of domain generalization, we then address the reliance of simulated ECG data on patient-specific geometrical models by learning to generalize simulated ECG data across subjects, before transferring them to clinical data. Evaluated on in-vivo premature ventricular contraction (PVC) patients, we demonstrate the feasibility of utilizing a large number of offline simulated ECG datasets to enable the prediction of the origin of arrhythmia with only a small number of clinical ECG data on a new patient

Location of parasympathetic innervation regions from electrograms to guide atrial fibrillation ablation therapy: an in silico modeling study

The autonomic nervous system (ANS) plays an essential role in the generation and maintenance of cardiac arrhythmias. The cardiac ANS can be divided into its extrinsic and intrinsic components, with the latter being organized in an epicardial neural network of interconnecting axons and clusters of autonomic ganglia called ganglionated plexi (GPs). GP ablation has been associated with a decreased risk of atrial fibrillation (AF) recurrence, but the accurate location of GPs is required for ablation to be effective. Although GP stimulation triggers both sympathetic and parasympathetic ANS branches, a predominance of parasympathetic activity has been shown. This study aims was to develop a method to locate atrial parasympathetic innervation sites based on measurements from a grid of electrograms (EGMs). Electrophysiological models representative of non-AF, paroxysmal AF (PxAF), and persistent AF (PsAF) tissues were developed. Parasympathetic effects were modeled by increasing the concentration of the neurotransmitter acetylcholine (ACh) in randomly distributed circles across the tissue. Different circle sizes of ACh and fibrosis geometries were considered, accounting for both uniform diffuse and non-uniform diffuse fibrosis. Computational simulations were performed, from which unipolar EGMs were computed in a 16 × 1 6 electrode mesh. Different distances of the electrodes to the tissue (0.5, 1, and 2 mm) and noise levels with signal-to-noise ratio (SNR) values of 0, 5, 10, 15, and 20 dB were tested. The amplitude of the atrial EGM repolarization wave was found to be representative of the presence or absence of ACh release sites, with larger positive amplitudes indicating that the electrode was placed over an ACh region. Statistical analysis was performed to identify the optimal thresholds for the identification of ACh sites. In all non-AF, PxAF, and PsAF tissues, the repolarization amplitude rendered successful identification. The algorithm performed better in the absence of fibrosis or when fibrosis was uniformly diffuse, with a mean accuracy of 0.94 in contrast with a mean accuracy of 0.89 for non-uniform diffuse fibrotic cases. The algorithm was robust against noise and worked for the tested ranges of electrode-to-tissue distance. In conclusion, the results from this study support the feasibility to locate atrial parasympathetic innervation sites from the amplitude of repolarization wave