Search

Purpose: Myofiber strain, Eff , is a mechanistically relevant metric of cardiac cell shortening and is expected to be spatially uniform in healthy populations, making it a prime candidate for the evaluation of local cardiomyocyte contractility. In this study, a new, efficient pipeline was proposed to combine microstructural cDTI and functional DENSE data in order to estimate Eff in vivo.Methods: Thirty healthy volunteers were scanned with three long-axis (LA) and three short-axis (SA) DENSE slices using 2D displacement encoding and one SA slice of cDTI. The total acquisition time was 11 minutes ± 3 minutes across volunteers. The pipeline first generates 3D SA displacements from all DENSE slices which are then combined with cDTI data to generate a cine of myofiber orientations and compute Eff . The precision of the post-processing pipeline was assessed using a computational phantom study. Transmural myofiber strain was compared to circumferential strain, Ecc , in healthy volunteers using a Wilcoxon sign rank test.Results: In vivo, computed Eff was found uniform transmurally compared to Ecc (-0.14[-0.15, -0.12] vs -0.18 [-0.20, -0.16], P < .001, -0.14 [-0.16, -0.12] vs -0.16 [-0.17, -0.13], P < .001 and -0.14 [-0.16, -0.12] vs Ecc_C = -0.14 [-0.15, -0.11], P = .002, Eff_C vs Ecc_C in the endo, mid, and epi layers, respectively).Conclusion: We demonstrate that it is possible to measure in vivo myofiber strain in a healthy human population in 10 minutes per subject. Myofiber strain was observed to be spatially uniform in healthy volunteers making it a potential biomarker for the evaluation of local cardiomyocyte contractility in assessing cardiovascular dysfunction.

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.

Computational models of cardiac contraction can provide critical insight into cardiac function and dysfunction. A necessary step before employing these computational models is their validation. Here we propose a series of validation criteria based on left ventricular (LV) global (ejection fraction and twist) and local (strains in a cylindrical coordinate system, aggregate cardiomyocyte shortening, and low myocardial compressibility) MRI measures to characterize LV motion and deformation during contraction. These validation criteria are used to evaluate an LV finite element model built from subject-specific anatomy and aggregate cardiomyocyte orientations reconstructed from diffusion tensor MRI. We emphasize the key role of the simulation boundary conditions in approaching the physiologically correct motion and strains during contraction. We conclude by comparing the global and local validation criteria measures obtained using two different boundary conditions: the first constraining the LV base and the second taking into account the presence of the pericardium, which leads to greatly improved motion and deformation.

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\pm 0.01\)), right ventricle (RV) cavity (\(0.80\pm 0.05\)), and the myocardium (\(0.72\pm 0.04\)). The \(R^2\) 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.

In vivo cardiac microstructure acquired using cardiac diffusion tensor imaging (cDTI) is a critical component of patient-specific models of cardiac electrophysiology and mechanics. In order to limit bulk motion artifacts and acquisition time, cDTI microstructural data is acquired at a single cardiac phase necessitating registration to the reference configuration on which the patient-specific computational models are based. Herein, we propose a method to register subject-specific microstructural data to an arbitrary cardiac phase using measured cardiac displacements. We validate our approach using a subject-specific computational phantom based on data from human subjects. Compared to a geometry-based non-rigid registration method, the displacement-based registration leads to improved accuracy (less than 1° versus 10° average median error in cardiomyocyte angular differences) and tighter confidence interval (3° versus 65° average upper limit of the 95% confidence interval).

Metrics of regional myocardial function can detect the onset of cardiovascular disease, evaluate the response to therapy, and provide mechanistic insight into cardiac dysfunction. Knowledge of local myocardial microstructure is necessary to distinguish between isotropic and anisotropic contributions of local deformation and to quantify myofiber kinematics, a microstructurally anchored measure of cardiac function. Using a computational model we combine in vivo cardiac displacement and diffusion tensor data to evaluate pointwise the deformation gradient tensor and isotropic and anisotropic deformation invariants. In discussing the imaging methods and the model construction, we identify potential improvements to increase measurement accuracy. We conclude by demonstrating the applicability of our method to compute myofiber strain in five healthy volunteers.