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

MR perfusion source mapping depicts venous territories and reveals perfusion modulation during neural activation

The cerebral venous system plays a crucial role in neurological and vascular conditions, yet its hemodynamics remain underexplored due to its complexity and variability across individuals. To address this, we develop a venous perfusion source mapping method using Displacement Spectrum MRI, a non-contrast technique that leverages blood water as an endogenous tracer. Our technique encodes spatial information into the magnetization of blood water spins during tagging and detects it once the tagged blood reaches the brain's surface, where the signal-to-noise ratio is 3-4 times higher. We resolve the sources of blood entering the imaging slice across short (10 ms) to long (3 s) evolution times, effectively capturing perfusion sources in reverse. This approach enables the measurement of slow venous blood flow, including potential contributions from capillary beds and surrounding tissue. We demonstrate perfusion source mapping in the superior cerebral veins, verify its sensitivity to global perfusion modulation induced by caffeine, and establish its specificity by showing repeatable local perfusion modulation during neural activation. From all blood within the imaging slice, our method localizes the portion originating from an activated region upstream.

In-silico heart model phantom to validate cardiac strain imaging

The quantification of cardiac strains as structural indices of cardiac function has a growing prevalence in clinical diagnosis. However, the highly heterogeneous four-dimensional (4D) cardiac motion challenges accurate "regional" strain quantification and leads to sizable differences in the estimated strains depending on the imaging modality and post-processing algorithm, limiting the translational potential of strains as incremental biomarkers of cardiac dysfunction. There remains a crucial need for a feasible benchmark that successfully replicates complex 4D cardiac kinematics to determine the reliability of strain calculation algorithms. In this study, we propose an in-silico heart phantom derived from finite element (FE) simulations to validate the quantification of 4D regional strains. First, as a proof-of-concept exercise, we created synthetic magnetic resonance (MR) images for a hollow thick-walled cylinder under pure torsion with an exact solution and demonstrated that "ground-truth" values can be recovered for the twist angle, which is also a key kinematic index in the heart. Next, we used mouse-specific FE simulations of cardiac kinematics to synthesize dynamic MR images by sampling various sectional planes of the left ventricle (LV). Strains were calculated using our recently developed non-rigid image registration (NRIR) framework in both problems. Moreover, we studied the effects of image quality on distorting regional strain calculations by conducting in-silico experiments for various LV configurations. Our studies offer a rigorous and feasible tool to standardize regional strain calculations to improve their clinical impact as incremental biomarkers.

In-silico heart model phantom to validate cardiac strain imaging

The quantification of cardiac strains as structural indices of cardiac function has a growing prevalence in clinical diagnosis. However, the highly heterogeneous four-dimensional (4D) cardiac motion challenges accurate "regional" strain quantification and leads to sizable differences in the estimated strains depending on the imaging modality and post-processing algorithm, limiting the translational potential of strains as incremental biomarkers of cardiac dysfunction. There remains a crucial need for a feasible benchmark that successfully replicates complex 4D cardiac kinematics to determine the reliability of strain calculation algorithms. In this study, we propose an in-silico heart phantom derived from finite element (FE) simulations to validate the quantification of 4D regional strains. First, as a proof-of-concept exercise, we created synthetic magnetic resonance (MR) images for a hollow thick-walled cylinder under pure torsion with an exact solution and demonstrated that "ground-truth" values can be recovered for the twist angle, which is also a key kinematic index in the heart. Next, we used mouse-specific FE simulations of cardiac kinematics to synthesize dynamic MR images by sampling various sectional planes of the left ventricle (LV). Strains were calculated using our recently developed non-rigid image registration (NRIR) framework in both problems. Moreover, we studied the effects of image quality on distorting regional strain calculations by conducting in-silico experiments for various LV configurations. Our studies offer a rigorous and feasible tool to standardize regional strain calculations to improve their clinical impact as incremental biomarkers.

Myocardial mesostructure and mesofunction

The complex and highly organized structural arrangement of some five billion cardiomyocytes directs the coordinated electrical activity and mechanical contraction of the human heart. The characteristic transmural change in cardiomyocyte orientation underlies base-to-apex shortening, circumferential shortening, and left ventricular torsion during contraction. Individual cardiomyocytes shorten ∼15% and increase in diameter ∼8%. Remarkably, however, the left ventricular wall thickens by up to 30-40%. To accommodate this, the myocardium must undergo significant structural rearrangement during contraction. At the mesoscale, collections of cardiomyocytes are organized into sheetlets, and sheetlet shear is the fundamental mechanism of rearrangement that produces wall thickening. Herein, we review the histological and physiological studies of myocardial mesostructure that have established the sheetlet shear model of wall thickening. Recent developments in tissue clearing techniques allow for imaging of whole hearts at the cellular scale, whereas magnetic resonance imaging (MRI) and computed tomography (CT) can image the myocardium at the mesoscale (100 µm to 1 mm) to resolve cardiomyocyte orientation and organization. Through histology, cardiac diffusion tensor imaging (DTI), and other modalities, mesostructural sheetlets have been confirmed in both animal and human hearts. Recent in vivo cardiac DTI methods have measured reorientation of sheetlets during the cardiac cycle. We also examine the role of pathological cardiac remodeling on sheetlet organization and reorientation, and the impact this has on ventricular function and dysfunction. We also review the unresolved mesostructural questions and challenges that may direct future work in the field.

Comparison of interpolation methods of predominant cardiomyocyte orientation from in vivo and ex vivo cardiac diffusion tensor imaging data

Cardiac electrophysiology and cardiac mechanics both depend on the average cardiomyocyte long-axis orientation. In the realm of personalized medicine, knowledge of the patient-specific changes in cardiac microstructure plays a crucial role. Patient-specific computational modelling has emerged as a tool to better understand disease progression. In vivo cardiac diffusion tensor imaging (cDTI) is a vital tool to non-destructively measure the average cardiomyocyte long-axis orientation in the heart. However, cDTI suffers from long scan times, rendering volumetric, high-resolution acquisitions challenging. Consequently, interpolation techniques are needed to populate bio-mechanical models with patient-specific average cardiomyocyte long-axis orientations. In this work, we compare five interpolation techniques applied to in vivo and ex vivo porcine input data. We compare two tensor interpolation approaches, one rule-based approximation, and two data-driven, low-rank models. We demonstrate the advantage of tensor interpolation techniques, resulting in lower interpolation errors than do low-rank models and rule-based methods adapted to cDTI data. In an ex vivo comparison, we study the influence of three imaging parameters that can be traded off against acquisition time: in-plane resolution, signal to noise ratio, and number of acquired short-axis imaging slices.

Reproducibility of global and segmental myocardial strain using cine DENSE at 3 T: a multicenter cardiovascular magnetic resonance study in healthy subjects and patients with heart disease

BACKGROUND

While multiple cardiovascular magnetic resonance (CMR) methods provide excellent reproducibility of global circumferential and global longitudinal strain, achieving highly reproducible segmental strain is more challenging. Previous single-center studies have demonstrated excellent reproducibility of displacement encoding with stimulated echoes (DENSE) segmental circumferential strain. The present study evaluated the reproducibility of DENSE for measurement of whole-slice or global circumferential (E_cc), longitudinal (E_ll) and radial (E_rr) strain, torsion, and segmental E_cc at multiple centers.

METHODS

Six centers participated and a total of 81 subjects were studied, including 60 healthy subjects and 21 patients with various types of heart disease. CMR utilized 3 T scanners, and cine DENSE images were acquired in three short-axis planes and in the four-chamber long-axis view. During one imaging session, each subject underwent two separate DENSE scans to assess inter-scan reproducibility. Each subject was taken out of the scanner and repositioned between the scans. Intra-user, inter-user-same-site, inter-user-different-site, and inter-user-Human-Deep-Learning (DL) comparisons assessed the reproducibility of different users analyzing the same data. Inter-scan comparisons assessed the reproducibility of DENSE from scan to scan. The reproducibility of whole-slice or global E_cc, E_ll and E_rr, torsion, and segmental E_cc were quantified using Bland-Altman analysis, the coefficient of variation (CV), and the intraclass correlation coefficient (ICC). CV was considered excellent for CV ≤ 10%, good for 10% < CV ≤ 20%, fair for 20% < CV ≤ 40%, and poor for CV > 40. ICC values were considered excellent for ICC > 0.74, good for ICC 0.6 < ICC ≤ 0.74, fair for ICC 0.4 < ICC ≤ 0.59, poor for ICC < 0.4.

RESULTS

Based on CV and ICC, segmental E_cc provided excellent intra-user, inter-user-same-site, inter-user-different-site, inter-user-Human-DL reproducibility and good-excellent inter-scan reproducibility. Whole-slice E_cc and global E_ll provided excellent intra-user, inter-user-same-site, inter-user-different-site, inter-user-Human-DL and inter-scan reproducibility. The reproducibility of torsion was good-excellent for all comparisons. For whole-slice E_rr, CV was in the fair-good range, and ICC was in the good-excellent range.

CONCLUSIONS

Multicenter data show that 3 T CMR DENSE provides highly reproducible whole-slice and segmental E_cc, global E_ll, and torsion measurements in healthy subjects and heart disease patients.

Reproducibility of global and segmental myocardial strain using cine DENSE at 3 T: a multicenter cardiovascular magnetic resonance study in healthy subjects and patients with heart disease

BACKGROUND

While multiple cardiovascular magnetic resonance (CMR) methods provide excellent reproducibility of global circumferential and global longitudinal strain, achieving highly reproducible segmental strain is more challenging. Previous single-center studies have demonstrated excellent reproducibility of displacement encoding with stimulated echoes (DENSE) segmental circumferential strain. The present study evaluated the reproducibility of DENSE for measurement of whole-slice or global circumferential (Ecc), longitudinal (Ell) and radial (Err) strain, torsion, and segmental Ecc at multiple centers.

METHODS

Six centers participated and a total of 81 subjects were studied, including 60 healthy subjects and 21 patients with various types of heart disease. CMR utilized 3 T scanners, and cine DENSE images were acquired in three short-axis planes and in the four-chamber long-axis view. During one imaging session, each subject underwent two separate DENSE scans to assess inter-scan reproducibility. Each subject was taken out of the scanner and repositioned between the scans. Intra-user, inter-user-same-site, inter-user-different-site, and inter-user-Human-Deep-Learning (DL) comparisons assessed the reproducibility of different users analyzing the same data. Inter-scan comparisons assessed the reproducibility of DENSE from scan to scan. The reproducibility of whole-slice or global Ecc, Ell and Err, torsion, and segmental Ecc were quantified using Bland-Altman analysis, the coefficient of variation (CV), and the intraclass correlation coefficient (ICC). CV was considered excellent for CV ≤ 10%, good for 10% < CV ≤ 20%, fair for 20% < CV ≤ 40%, and poor for CV > 40. ICC values were considered excellent for ICC > 0.74, good for ICC 0.6 < ICC ≤ 0.74, fair for ICC 0.4 < ICC ≤ 0.59, poor for ICC < 0.4.

RESULTS

Based on CV and ICC, segmental Ecc provided excellent intra-user, inter-user-same-site, inter-user-different-site, inter-user-Human-DL reproducibility and good-excellent inter-scan reproducibility. Whole-slice Ecc and global Ell provided excellent intra-user, inter-user-same-site, inter-user-different-site, inter-user-Human-DL and inter-scan reproducibility. The reproducibility of torsion was good-excellent for all comparisons. For whole-slice Err, CV was in the fair-good range, and ICC was in the good-excellent range.

CONCLUSIONS

Multicenter data show that 3 T CMR DENSE provides highly reproducible whole-slice and segmental Ecc, global Ell, and torsion measurements in healthy subjects and heart disease patients.

Cardiac Diffusion Tensor Biomarkers of Chronic Infarction Based on In Vivo Data

In vivo cardiac diffusion tensor imaging (cDTI) data were acquired in swine subjects six to ten weeks post-myocardial infarction (MI) to identify microstructural-based biomarkers of MI. Diffusion tensor invariants, diffusion tensor eigenvalues, and radial diffusivity (RD) are evaluated in the infarct, border, and remote myocardium, and compared with extracellular volume fraction (ECV) and native T1 values. Additionally, to aid the interpretation of the experimental results, the diffusion of water molecules was numerically simulated as a function of ECV. Finally, findings based on in vivo measures were confirmed using higher-resolution and higher signal-to-noise data acquired ex vivo in the same subjects. Mean diffusivity, diffusion tensor eigenvalues, and RD increased in the infarct and border regions compared to remote myocardium, while fractional anisotropy decreased. Secondary (e₂) and tertiary (e₃) eigenvalues increased more significantly than the primary eigenvalue in the infarct and border regions. These findings were confirmed by the diffusion simulations. Although ECV presented the largest increase in infarct and border regions, e₂, e₃, and RD increased the most among non-contrast-based biomarkers. RD is of special interest as it summarizes the changes occurring in the radial direction and may be more robust than e₂ or e₃ alone.

Personalization of biomechanical simulations of the left ventricle by in-vivo cardiac DTI data: Impact of fiber interpolation methods

Simulations of cardiac electrophysiology and mechanics have been reported to be sensitive to the microstructural anisotropy of the myocardium. Consequently, a personalized representation of cardiac microstructure is a crucial component of accurate, personalized cardiac biomechanical models. In-vivo cardiac Diffusion Tensor Imaging (cDTI) is a non-invasive magnetic resonance imaging technique capable of probing the heart's microstructure. Being a rather novel technique, issues such as low resolution, signal-to noise ratio, and spatial coverage are currently limiting factors. We outline four interpolation techniques with varying degrees of data fidelity, different amounts of smoothing strength, and varying representation error to bridge the gap between the sparse in-vivo data and the model, requiring a 3D representation of microstructure across the myocardium. We provide a workflow to incorporate in-vivo myofiber orientation into a left ventricular model and demonstrate that personalized modelling based on fiber orientations from in-vivo cDTI data is feasible. The interpolation error is correlated with a trend in personalized parameters and simulated physiological parameters, strains, and ventricular twist. This trend in simulation results is consistent across material parameter settings and therefore corresponds to a bias introduced by the interpolation method. This study suggests that using a tensor interpolation approach to personalize microstructure with in-vivo cDTI data, reduces the fiber uncertainty and thereby the bias in the simulation results.

Using synthetic data generation to train a cardiac motion tag tracking neural network

A CNN based method for cardiac MRI tag tracking was developed and validated. A synthetic data simulator was created to generate large amounts of training data using natural images, a Bloch equation simulation, a broad range of tissue properties, and programmed ground-truth motion. The method was validated using both an analytical deforming cardiac phantom and in vivo data with manually tracked reference motion paths. In the analytical phantom, error was investigated relative to SNR, and accurate results were seen for SNR>10 (displacement error <0.3 mm). Excellent agreement was seen in vivo for tag locations (mean displacement difference = -0.02 pixels, 95% CI [-0.73, 0.69]) and calculated cardiac circumferential strain (mean difference = 0.006, 95% CI [-0.012, 0.024]). Automated tag tracking with a CNN trained on synthetic data is both accurate and precise.

DiSpect: Displacement spectrum imaging of flow and tissue perfusion using spin-labeling and stimulated echoes

PURPOSE

We propose a new method, displacement spectrum (DiSpect) imaging, for probing in vivo complex tissue dynamics such as motion, flow, diffusion, and perfusion. Based on stimulated echoes and image phase, our flexible approach enables observations of the spin dynamics over short (milliseconds) to long (seconds) evolution times.

METHODS

The DiSpect method is a Fourier-encoded variant of displacement encoding with stimulated echoes, which encodes bulk displacement of spins that occurs between tagging and imaging in the image phase. However, this method fails to capture partial volume effects as well as blood flow. The DiSpect variant mitigates this by performing multiple scans with increasing displacement-encoding steps. Fourier analysis can then resolve the multidimensional spectrum of displacements that spins exhibit over the mixing time. In addition, repeated imaging following tagging can capture dynamic displacement spectra with increasing mixing times.

RESULTS

We demonstrate properties of DiSpect MRI using flow phantom experiments as well as in vivo brain scans. Specifically, the ability of DiSpect to perform retrospective vessel-selective perfusion imaging at multiple mixing times is highlighted.

CONCLUSION

The DiSpect variant is a new tool in the arsenal of MRI techniques for probing complex tissue dynamics. The flexibility and the rich information it provides open the possibility of alternative ways to quantitatively measure numerous complex spin dynamics, such as flow and perfusion within a single exam.

A kinematic model-based analysis framework for 3D Cine-DENSE-validation with an axially compressed gel phantom and application in sheep before and after antero-apical myocardial infarction

PURPOSE

Myocardial strain is increasingly used to assess left ventricular (LV) function. Incorporation of LV deformation into finite element (FE) modeling environment with subsequent strain calculation will allow analysis to reach its full potential. We describe a new kinematic model-based analysis framework (KMAF) to calculate strain from 3D cine-DENSE (displacement encoding with stimulated echoes) MRI.

METHODS

Cine-DENSE allows measurement of 3D myocardial displacement with high spatial accuracy. The KMAF framework uses cine cardiovascular magnetic resonance (CMR) to facilitate cine-DENSE segmentation, interpolates cine-DENSE displacement, and kinematically deforms an FE model to calculate strain. This framework was validated in an axially compressed gel phantom and applied in 10 healthy sheep and 5 sheep after myocardial infarction (MI).

RESULTS

Excellent Bland-Altman agreement of peak circumferential (E_cc ) and longitudinal (E_ll ) strain (mean difference = 0.021 ± 0.04 and -0.006 ± 0.03, respectively), was found between KMAF estimates and idealized FE simulation. E_rr had a mean difference of -0.014 but larger variation (±0.12). Cine-DENSE estimated end-systolic (ES) E_cc , E_ll and E_rr exhibited significant spatial variation for healthy sheep. Displacement magnitude was reduced on average by 27%, 42%, and 56% after MI in the remote, adjacent and MI regions, respectively.

CONCLUSIONS

The KMAF framework allows accurate calculation of 3D LV E_cc and E_ll from cine-DENSE.

Arbitrary Point Tracking with Machine Learning to Measure Cardiac Strains in Tagged MRI

Cardiac tagged MR images allow for deformation fields to be measured in the heart by tracking the motion of tag lines throughout the cardiac cycle. Machine learning (ML) algorithms enable accurate and robust tracking of tag lines. Herein, the use of a massive synthetic physics-driven training dataset with known ground truth was used to train an ML network to enable tracking any number of points at arbitrary positions rather than anchored to the tag lines themselves. The tag tracking and strain calculation methods were investigated in a computational deforming cardiac phantom with known (ground truth) strain values. This enabled both tag tracking and strain accuracy to be characterized for a range of image acquisition and tag tracking parameters. The methods were also tested on in vivo volunteer data. Median tracking error was <0.26mm in the computational phantom, and strain measurements were improved in vivo when using the arbitrary point tracking for a standard clinical protocol.

Fully-automated global and segmental strain analysis of DENSE cardiovascular magnetic resonance using deep learning for segmentation and phase unwrapping

BACKGROUND

Cardiovascular magnetic resonance (CMR) cine displacement encoding with stimulated echoes (DENSE) measures heart motion by encoding myocardial displacement into the signal phase, facilitating high accuracy and reproducibility of global and segmental myocardial strain and providing benefits in clinical performance. While conventional methods for strain analysis of DENSE images are faster than those for myocardial tagging, they still require manual user assistance. The present study developed and evaluated deep learning methods for fully-automatic DENSE strain analysis.

METHODS

Convolutional neural networks (CNNs) were developed and trained to (a) identify the left-ventricular (LV) epicardial and endocardial borders, (b) identify the anterior right-ventricular (RV)-LV insertion point, and (c) perform phase unwrapping. Subsequent conventional automatic steps were employed to compute strain. The networks were trained using 12,415 short-axis DENSE images from 45 healthy subjects and 19 heart disease patients and were tested using 10,510 images from 25 healthy subjects and 19 patients. Each individual CNN was evaluated, and the end-to-end fully-automatic deep learning pipeline was compared to conventional user-assisted DENSE analysis using linear correlation and Bland Altman analysis of circumferential strain.

RESULTS

LV myocardial segmentation U-Nets achieved a DICE similarity coefficient of 0.87 ± 0.04, a Hausdorff distance of 2.7 ± 1.0 pixels, and a mean surface distance of 0.41 ± 0.29 pixels in comparison with manual LV myocardial segmentation by an expert. The anterior RV-LV insertion point was detected within 1.38 ± 0.9 pixels compared to manually annotated data. The phase-unwrapping U-Net had similar or lower mean squared error vs. ground-truth data compared to the conventional path-following method for images with typical signal-to-noise ratio (SNR) or low SNR (p < 0.05), respectively. Bland-Altman analyses showed biases of 0.00 ± 0.03 and limits of agreement of - 0.04 to 0.05 or better for deep learning-based fully-automatic global and segmental end-systolic circumferential strain vs. conventional user-assisted methods.

CONCLUSIONS

Deep learning enables fully-automatic global and segmental circumferential strain analysis of DENSE CMR providing excellent agreement with conventional user-assisted methods. Deep learning-based automatic strain analysis may facilitate greater clinical use of DENSE for the quantification of global and segmental strain in patients with cardiac disease.

Myofiber strain in healthy humans using DENSE and cDTI

PURPOSE

Myofiber strain, E_ff , 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 E_ff 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 E_ff . The precision of the post-processing pipeline was assessed using a computational phantom study. Transmural myofiber strain was compared to circumferential strain, E_cc , in healthy volunteers using a Wilcoxon sign rank test.

RESULTS

In vivo, computed E_ff was found uniform transmurally compared to E_cc (-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 E_{cc_C} = -0.14 [-0.15, -0.11], P = .002, E_{ff_C} vs E_{cc_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.

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.