TIME RESOLVED DISPLACEMENT-BASED REGISTRATION OF IN VIVO CDTI CARDIOMYOCYTE ORIENTATIONS

Characterizing variability in passive myocardial stiffness in healthy human left ventricles using personalized MRI and finite element modeling

Abnormal passive stiffness of the heart muscle (myocardium) is evident in the pathophysiology of several cardiovascular diseases, making it an important indicator of heart health. Recent advancements in cardiac imaging and biophysical modeling now enable more effective evaluation of this biomarker. Estimating passive myocardial stiffness can be accomplished through an MRI-based approach that requires comprehensive subject-specific input data. This includes the gross cardiac geometry (e.g. from conventional cine imaging), regional diastolic kinematics (e.g. from tagged MRI), microstructural configuration (e.g. from diffusion tensor imaging), and ventricular diastolic pressure, whether invasively measured or non-invasively estimated. Despite the progress in cardiac biomechanics simulations, developing a framework to integrate multiphase and multimodal cardiac MRI data for estimating passive myocardial stiffness has remained a challenge. Moreover, the sensitivity of estimated passive myocardial stiffness to input data has not been fully explored. This study aims to: (1) develop a framework for integrating subject-specific in vivo MRI data into in silico left ventricular finite element models to estimate passive myocardial stiffness, (2) apply the framework to estimate the passive myocardial stiffness of multiple healthy subjects under assumed filling pressure, and (3) assess the sensitivity of these estimates to loading conditions and myofiber orientations. This work contributes toward the establishment of a range of reference values for material parameters of passive myocardium in healthy human subjects. Notably, in this study, beat-to-beat variation in left ventricular end-diastolic pressure was found to have a greater influence on passive myocardial material parameter estimation than variation in fiber orientation.

Estimating cardiomyofiber strain in vivo by solving a computational model

Since heart contraction results from the electrically activated contraction of millions of cardiomyocytes, a measure of cardiomyocyte shortening mechanistically underlies cardiac contraction. In this work we aim to measure preferential aggregate cardiomyocyte ("myofiber") strains based on Magnetic Resonance Imaging (MRI) data acquired to measure both voxel-wise displacements through systole and myofiber orientation. In order to reduce the effect of experimental noise on the computed myofiber strains, we recast the strains calculation as the solution of a boundary value problem (BVP). This approach does not require a calibrated material model, and consequently is independent of specific myocardial material properties. The solution to this auxiliary BVP is the displacement field corresponding to assigned values of myofiber strains. The actual myofiber strains are then determined by minimizing the difference between computed and measured displacements. The approach is validated using an analytical phantom, for which the ground-truth solution is known. The method is applied to compute myofiber strains using in vivo displacement and myofiber MRI data acquired in a mid-ventricular left ventricle section in N=8 swine subjects. The proposed method shows a more physiological distribution of myofiber strains compared to standard approaches that directly differentiate the displacement field.

Estimating Aggregate Cardiomyocyte Strain Using In Vivo Diffusion and Displacement Encoded MRI

Changes in left ventricular (LV) aggregate cardiomyocyte orientation and deformation underlie cardiac function and dysfunction. As such, in vivo aggregate cardiomyocyte "myofiber" strain ( [Formula: see text]) has mechanistic significance, but currently there exists no established technique to measure in vivo [Formula: see text]. The objective of this work is to describe and validate a pipeline to compute in vivo [Formula: see text] from magnetic resonance imaging (MRI) data. Our pipeline integrates LV motion from multi-slice Displacement ENcoding with Stimulated Echoes (DENSE) MRI with in vivo LV microstructure from cardiac Diffusion Tensor Imaging (cDTI) data. The proposed pipeline is validated using an analytical deforming heart-like phantom. The phantom is used to evaluate 3D cardiac strains computed from a widely available, open-source DENSE Image Analysis Tool. Phantom evaluation showed that a DENSE MRI signal-to-noise ratio (SNR) ≥20 is required to compute [Formula: see text] with near-zero median strain bias and within a strain tolerance of 0.06. Circumferential and longitudinal strains are also accurately measured under the same SNR requirements, however, radial strain exhibits a median epicardial bias of -0.10 even in noise-free DENSE data. The validated framework is applied to experimental DENSE MRI and cDTI data acquired in eight ( N=8 ) healthy swine. The experimental study demonstrated that [Formula: see text] has decreased transmural variability compared to radial and circumferential strains. The spatial uniformity and mechanistic significance of in vivo [Formula: see text] make it a compelling candidate for characterization and early detection of cardiac dysfunction.

High-Resolution Ex Vivo Microstructural MRI After Restoring Ventricular Geometry via 3D Printing

Computational modeling of the heart requires accurately incorporating both gross anatomical detail and local microstructural information. Together, these provide the necessary data to build 3D meshes for simulation of cardiac mechanics and electrophysiology. Recent MRI advances make it possible to measure detailed heart motion in vivo, but in vivo microstructural imaging of the heart remains challenging. Consequently, the most detailed measurements of microstructural organization and microanatomical infarct details are obtained ex vivo. The objective of this work was to develop and evaluate a new method for restoring ex vivo ventricular geometry to match the in vivo configuration. This approach aids the integration of high-resolution ex vivo microstructural information with in vivo motion measurements. The method uses in vivo cine imaging to generate surface meshes, then creates a 3D printed left ventricular (LV) blood pool cast and a pericardial mold to restore the ex vivo cardiac geometry to a mid-diastasis reference configuration. The method was evaluated in healthy (N = 7) and infarcted (N = 3) swine. Dice similarity coefficients were calculated between in vivo and ex vivo images for the LV cavity (0.93 ± 0.01), right ventricle (RV) cavity (0.80 ± 0.05), and the myocardium (0.72 ± 0.04). The R ² coefficient between in vivo and ex vivo LV and RV cavity volumes were 0.95 and 0.91, respectively. These results suggest that this method adequately restores ex vivo geometry to match in vivo geometry. This approach permits a more precise incorporation of high-resolution ex vivo anatomical and microstructural data into computational models that use in vivo data for simulation of cardiac mechanics and electrophysiology.