Modeling cardiac muscle fibers in ventricular and atrial electrophysiology simulations

Since myocardial fibers drive the electric signal propagation throughout the myocardium, accurately modeling their arrangement is essential for simulating heart electrophysiology (EP). Rule-Based-Methods (RBMs) represent a commonly used strategy to include cardiac fibers in computational models. A particular class of such methods is known as Laplace-Dirichlet-Rule-Based-Methods (LDRBMs) since they rely on the solution of Laplace problems. In this work we provide a unified framework, based on LDRBMs, for generating full heart muscle fibers. First, we review existing ventricular LDRBMs providing a communal mathematical description and introducing also some modeling improvements with respect to the existing literature. We then carry out a systematic comparison of LDRBMs based on meaningful biomarkers produced by numerical EP simulations. Next we propose, for the first time, a LDRBM to be used for generating atrial fibers. The new method, tested both on idealized and realistic atrial models, can be applied to any arbitrary geometries. Finally, we present numerical results obtained in a realistic whole heart where fibers are included for all the four chambers using the discussed LDRBMs

Characterization and modeling of the human left atrium using optical coherence tomography

With current needs to better understand the interaction between atrial tissue microstructure and atrial fibrillation dynamics, micrometer scale imaging with optical coherence tomography has significant potential to provide further insight on arrhythmia mechanisms and improve treatment guidance. However, optical coherence tomography imaging of cardiac tissue in humans is largely unexplored, and the ability of optical coherence tomography to identify the structural substrate of atrial fibrillation has not yet been investigated. Therefore, the objective of this thesis was to develop an optical coherence tomography imaging atlas of the human heart, study the utility of optical coherence tomography in providing useful features of human left atrial tissues, and develop a framework for optical coherence tomography-informed cardiac modeling that could be used to probe dynamics between electrophysiology and tissue structure. Human left atrial tissues were comprehensively imaged by optical coherence tomography for the first time, providing an imaging atlas that can guide identification of left atrial tissue features from optical coherence tomography imaging. Optical coherence tomography image features corresponding to myofiber and collagen fiber orientation, adipose tissue, endocardial thickness and composition, and venous media were established. Varying collagen fiber distributions in the myocardial sleeves were identified within the pulmonary veins. A scheme for mapping optical coherence tomography data of dissected left atrial tissues to a three-dimensional, anatomical model of the human left atrium was also developed, enabling the mapping of distributions of imaged adipose tissue and fiber orientation to the whole left atrial geometry. These results inform future applications of structural substrate mapping in the human left atrium using optical coherence tomography-integrated catheters, as well as potential directions of ex vivo optical coherence tomography atrial imaging studies. Additionally, we developed a workflow for creating optical mapping models of atrial tissue as informed by optical coherence tomography. Tissue geometry, fiber orientation, ablation lesion geometry, and heterogeneous tissue types were extracted from optical coherence tomography images and incorporated into tissue-specific meshes. Electrophysiological propagation was simulated and combined with photon scattering simulations to evaluate the influence of tissue-specific structure on electrical and optical mapping signals. Through tissue-specific modeling of myofiber orientation, ablation lesions, and heterogeneous tissue types, the influence of myofiber orientation on transmural activation, the relationship between fluorescent signals and lesion geometry, and the blurring of optical mapping signals in the presence of heterogeneous tissue types were investigated. By providing a comprehensive optical coherence tomography image database of the human left atrium and a workflow for developing optical coherence tomography-informed cardiac tissue models, this work establishes the foundation for utilizing optical coherence tomography to improve the structural substrate characterization of atrial fibrillation. Future developments include analysis of optical coherence tomography imaged tissue structure with respect to clinical presentation, development of automated processing to better leverage the large amount of imaging data, enhancements and validation of the modeling scheme, and in vivo evaluation of the left atrial structural substrate through optical coherence tomography-integrated catheter

Circ Arrhythm Electrophysiol

BackgroundAccurate knowledge of the human atrial fibrous structure is paramount in understanding the mechanisms of atrial electrical function in health and disease. Thus far such knowledge has been acquired from destructive sectioning, and there is a paucity of data regarding atrial fiber architecture variability in the human population.Methods and ResultsIn this study, we have developed a customized 3D diffusion tensor magnetic resonance imaging (DTMRI) sequence on a clinical scanner that makes it possible to image an entire intact human heart specimen ex vivo at sub-millimeter resolution. The data from eight human atrial specimens obtained with this technique present complete maps of the fibrous organization of the human atria. The findings demonstrate that the main features of atrial anatomy are mostly preserved across subjects, although the exact location and orientation of atrial bundles vary. Using the full tractography data, we were able to cluster, visualize, and characterize the distinct major bundles in the human atria. Further, quantitative characterization of the fiber angles across the atrial wall revealed that the transmural fiber angle distribution is heterogeneous throughout different regions of the atria.ConclusionsThe application of sub-millimeter DTMRI provides an unprecedented level of information regarding both human atrial structure as well as its inter-subject variability. The high resolution and fidelity of this data could enhance our understanding of structural contributions to atrial rhythm and pump disorders, and lead to improvements in their targeted treatment.DP1HL123271/DP/NCCDPHP CDC HHS/United StatesR01 HL142893/HL/NHLBI NIH HHS/United StatesDP1 HL123271/HL/NHLBI NIH HHS/United StatesImNIH/Intramural NIH HHS/United StatesR01 HL142496/HL/NHLBI NIH HHS/United StatesR01 HL126802/HL/NHLBI NIH HHS/United States2020-02-22T00:00:00Z27071829PMC70358847696vault:3467

lifex-fiber: an open tool for myofibers generation in cardiac computational models

Background: Modeling the whole cardiac function involves the solution of several complex multi-physics and multi-scale models that are highly computationally demanding, which call for simpler yet accurate, high-performance computational tools. Despite the efforts made by several research groups, no software for whole-heart fully coupled cardiac simulations in the scientific community has reached full maturity yet.Results: In this work we present life(x)-fiber, an innovative tool for the generation of myocardial fibers based on Laplace-Dirichlet Rule-Based Methods, which are the essential building blocks for modeling the electrophysiological, mechanical and electromechanical cardiac function, from single-chamber to whole-heart simulations. life(x)-fiber is the first publicly released module for cardiac simulations based on life(x), an open-source, high-performance Finite Element solver for multi-physics, multi-scale and multi-domain problems developed in the framework of the iHEART project, which aims at making in silico experiments easily reproducible and accessible to a wide community of users, including those with a background in medicine or bio-engineering.Conclusions: The tool presented in this document is intended to provide the scientific community with a computational tool that incorporates general state of the art models and solvers for simulating the cardiac function within a high-performance framework that exposes a user-and developer-friendly interface. This report comes with an extensive technical and mathematical documentation to welcome new users to the core structure of life(x)-fiber and to provide them with a possible approach to include the generated cardiac fibers into more sophisticated computational pipelines. In the near future, more modules will be successively published either as pre-compiled binaries for x86-64 Linux systems or as open source software

Fiber Organization has Little Effect on Electrical Activation Patterns during Focal Arrhythmias in the Left Atrium

Over the past two decades there has been a steady trend towards the development of realistic models of cardiac conduction with increasing levels of detail. However, making models more realistic complicates their personalization and use in clinical practice due to limited availability of tissue and cellular scale data. One such limitation is obtaining information about myocardial fiber organization in the clinical setting. In this study, we investigated a chimeric model of the left atrium utilizing clinically derived patient-specific atrial geometry and a realistic, yet foreign for a given patient fiber organization. We discovered that even significant variability of fiber organization had a relatively small effect on the spatio-temporal activation pattern during regular pacing. For a given pacing site, the activation maps were very similar across all fiber organizations tested

lifex-ep: a robust and efficient software for cardiac electrophysiology simulations

Background: Simulating the cardiac function requires the numerical solution of multi-physics and multi-scale mathematical models. This underscores the need for streamlined, accurate, and high-performance computational tools. Despite the dedicated endeavors of various research teams, comprehensive and user-friendly software programs for cardiac simulations, capable of accurately replicating both normal and pathological conditions, are still in the process of achieving full maturity within the scientific community. Results: This work introduces lifex-ep, a publicly available software for numerical simulations of the electrophysiology activity of the cardiac muscle, under both normal and pathological conditions. lifex-ep employs the monodomain equation to model the heart's electrical activity. It incorporates both phenomenological and second-generation ionic models. These models are discretized using the Finite Element method on tetrahedral or hexahedral meshes. Additionally, lifex-ep integrates the generation of myocardial fibers based on Laplace-Dirichlet Rule-Based Methods, previously released in Africa et al., 2023, within lifex-fiber. As an alternative, users can also choose to import myofibers from a file. This paper provides a concise overview of the mathematical models and numerical methods underlying lifex-ep, along with comprehensive implementation details and instructions for users. lifex-ep features exceptional parallel speedup, scaling efficiently when using up to thousands of cores, and its implementation has been verified against an established benchmark problem for computational electrophysiology. We showcase the key features of lifex-ep through various idealized and realistic simulations conducted in both normal and pathological scenarios. Furthermore, the software offers a user-friendly and flexible interface, simplifying the setup of simulations using self-documenting parameter files. Conclusions: lifex-ep provides easy access to cardiac electrophysiology simulations for a wide user community. It offers a computational tool that integrates models and accurate methods for simulating cardiac electrophysiology within a high-performance framework, while maintaining a user-friendly interface. lifex-ep represents a valuable tool for conducting in silico patient-specific simulations

Submillimeter diffusion tensor imaging and late gadolinium enhancement cardiovascular magnetic resonance of chronic myocardial infarction.

BackgroundKnowledge of the three-dimensional (3D) infarct structure and fiber orientation remodeling is essential for complete understanding of infarct pathophysiology and post-infarction electromechanical functioning of the heart. Accurate imaging of infarct microstructure necessitates imaging techniques that produce high image spatial resolution and high signal-to-noise ratio (SNR). The aim of this study is to provide detailed reconstruction of 3D chronic infarcts in order to characterize the infarct microstructural remodeling in porcine and human hearts.MethodsWe employed a customized diffusion tensor imaging (DTI) technique in conjunction with late gadolinium enhancement (LGE) cardiovascular magnetic resonance (CMR) on a 3T clinical scanner to image, at submillimeter resolution, myofiber orientation and scar structure in eight chronically infarcted porcine hearts ex vivo. Systematic quantification of local microstructure was performed and the chronic infarct remodeling was characterized at different levels of wall thickness and scar transmurality. Further, a human heart with myocardial infarction was imaged using the same DTI sequence.ResultsThe SNR of non-diffusion-weighted images was >100 in the infarcted and control hearts. Mean diffusivity and fractional anisotropy (FA) demonstrated a 43% increase, and a 35% decrease respectively, inside the scar tissue. Despite this, the majority of the scar showed anisotropic structure with FA higher than an isotropic liquid. The analysis revealed that the primary eigenvector orientation at the infarcted wall on average followed the pattern of original fiber orientation (imbrication angle mean: 1.96 ± 11.03° vs. 0.84 ± 1.47°, p = 0.61, and inclination angle range: 111.0 ± 10.7° vs. 112.5 ± 6.8°, p = 0.61, infarcted/control wall), but at a higher transmural gradient of inclination angle that increased with scar transmurality (r = 0.36) and the inverse of wall thickness (r = 0.59). Further, the infarcted wall exhibited a significant increase in both the proportion of left-handed epicardial eigenvectors, and in the angle incoherency. The infarcted human heart demonstrated preservation of primary eigenvector orientation at the thinned region of infarct, consistent with the findings in the porcine hearts.ConclusionsThe application of high-resolution DTI and LGE-CMR revealed the detailed organization of anisotropic infarct structure at a chronic state. This information enhances our understanding of chronic post-infarction remodeling in large animal and human hearts

High-Resolution Whole-Heart Imaging and Modeling for Studying Cardiac Arrhythmia

Cardiac arrhythmia is a life-threatening heart rhythm disorder affecting millions of people worldwide. The underlying structure of the heart plays an important role in cardiac activity and could promote rhythm disorders. Accurate knowledge of whole-heart cardiac geometry and microstructure in normal and disease hearts is essential for a complete understanding of the mechanisms of arrhythmias. This dissertation presents novel structural data at the whole-heart level aimed at advancing knowledge of cardiac structure in normal and infarcted hearts, and at constructing whole-heart computational models. A 3D diffusion tensor MRI (DTMRI) technique was implemented on a clinical scanner to image intact large animal and human hearts with high image quality and spatial resolution ex vivo. This method was first applied to reconstruct the 3D myofiber organization in 8 human atria nondestructively and at submillimeter resolution. The findings showed that the main features of atrial anatomy are mostly preserved across subjects despite variability in the exact location and orientation of the bundles. Further, we were able to cluster, visualize, and characterize the distinct major bundles in the human atria. Quantitative analysis of the fiber angles across the atrial wall revealed that the transmural fiber angle distribution is heterogeneous throughout the atria. We next studied microstructural remodeling in infarcted porcine and human hearts by combining DTMRI with high-resolution Late Gadolinium Enhancement imaging. This enabled us to provide reconstructions of both fiber architecture and scar distribution in infarcted hearts with an unprecedented level of detail, and to systematically quantify the transmural pattern of diffusion eigenvector orientation. Our results demonstrated that the fiber orientation is generally preserved inside the scar but at a higher transmural gradient of inclination angle. Lastly, we employed the obtained data to generate whole-heart computational models of infarcted hearts with detailed scar geometry and subject-specific fiber orientation. We used these models in simulations to investigate the contribution of the infarct microarchitecture to ventricular tachycardia. The simulation results showed that the reentry circuits traverse thin viable tissues with complex geometries located inside of the infarct. The high resolution of the images enabled 3D reconstruction and characterization of such structures

Relationship of Transmural Variations in Myofiber Contractility to Left Ventricular Ejection Fraction: Implications for Modeling Heart Failure Phenotype With Preserved Ejection Fraction

The pathophysiological mechanisms underlying preserved left ventricular (LV) ejection fraction (EF) in patients with heart failure and preserved ejection fraction (HFpEF) remain incompletely understood. We hypothesized that transmural variations in myofiber contractility with existence of subendocardial dysfunction and compensatory increased subepicardial contractility may underlie preservation of LVEF in patients with HFpEF. We quantified alterations in myocardial function in a mathematical model of the human LV that is based on the finite element method. The fiber-reinforced material formulation of the myocardium included passive and active properties. The passive material properties were determined such that the diastolic pressure-volume behavior of the LV was similar to that shown in published clinical studies of pressure-volume curves. To examine changes in active properties, we considered six scenarios: (1) normal properties throughout the LV wall; (2) decreased myocardial contractility in the subendocardium; (3) increased myocardial contractility in the subepicardium; (4) myocardial contractility decreased equally in all layers, (5) myocardial contractility decreased in the midmyocardium and subepicardium, (6) myocardial contractility decreased in the subepicardium. Our results indicate that decreased subendocardial contractility reduced LVEF from 53.2 to 40.5%. Increased contractility in the subepicardium recovered LVEF from 40.5 to 53.2%. Decreased contractility transmurally reduced LVEF and could not be recovered if subepicardial and midmyocardial contractility remained depressed. The computational results simulating the effects of transmural alterations in the ventricular tissue replicate the phenotypic patterns of LV dysfunction observed in clinical practice. In particular, data for LVEF, strain and displacement are consistent with previous clinical observations in patients with HFpEF, and substantiate the hypothesis that increased subepicardial contractility may compensate for subendocardial dysfunction and play a vital role in maintaining LVEF