Tracking myocardial motion from cine DENSE images using spatiotemporal phase unwrapping and temporal fitting

Fully automated and comprehensive MRI-based left-ventricular contractility analysis in post-chemotherapy breast cancer patients

OBJECTIVE

This study investigated the occurrence of cardiotoxicity-related left-ventricular (LV) contractile dysfunction in breast cancer patients following treatment with antineoplastic chemotherapy agents.

METHODS

A validated and automated MRI-based LV contractility analysis tool consisting of quantization-based boundary detection, unwrapping of image phases and the meshfree Radial Point Interpolation Method was used toward measuring LV chamber quantifications (LVCQ), three-dimensional strains and torsions in patients and healthy subjects. Data were acquired with the Displacement Encoding with Stimulated Echoes (DENSE) sequence on 21 female patients and 21 age-matched healthy females. Estimates of patient LVCQs from DENSE acquisitions were validated in comparison to similar steady-state free precession measurements and their strain results validated via Bland-Altman interobserver agreements. The occurrence of LV abnormalities was investigated via significant differences in contractility measurements (LVCQs, strains and torsions) between patients and healthy subjects.

RESULTS

Repeated measures analysis showed similarities between LVCQ measurements from DENSE and steady-state free precession, including cardiac output (4.7 ± 0.4 L, 4.6 ± 0.4 L, p = 0.8), and LV ejection fractions (59±6%, 58±5%, p = 0.2). Differences found between patients and healthy subjects included enlarged basal diameter (5.0 ± 0.5 cm vs 4.4 ± 0.5 cm, p < 0.01), apical torsion (6.0 ± 1.1° vs 9.7 ± 1.4°, p < 0.001) and global longitudinal strain (-0.15 ± 0.02 vs. -0.21 ± 0.04, p < 0.001), but not LV ejection fraction (59±6% vs. 63±6%, p = 0.1).

CONCLUSION

The results from the statistical analysis reveal the possibility of LV abnormalities in the post-chemotherapy patients via enlarged basal diameter and reduced longitudinal strain and torsion, in comparison to healthy subjects.

ADVANCES IN KNOWLEDGE

This study shows that subclinical LV abnormalities in post-chemotherapy breast cancer patients can be detected with an automated technique for the comprehensive analysis of contractile parameters.

Quantification of regional right ventricular strain in healthy rats using 3D spiral cine dense MRI

Statistical data from clinical studies suggests that right ventricular (RV) circumferential strain (Ecc) and longitudinal strain (Ell) are significant biomarkers for many cardiovascular diseases. However, a detailed and regional characterization of these strains in the RV is very limited. In the current study, RV images were obtained with 3D spiral cine DENSE MRI in healthy rats. An algorithm for surface growing was proposed in order to fit irregular topology. Specifically, a new custom plugin for the DENSEanalysis program, called 3D DENSE Plugin for Crescent Organ, was developed for surface reconstruction and precise segmentation of organs with sharp curvature, such as the murine RV. The RV free wall (RVFW) was divided into three longitudinal thirds (i.e., basal, middle, and apical) with each one partitioned into circumferential fourths (i.e., anterior, anteriorlateral, inferiorlateral and inferior). Peak systolic strains were quantified for each segment and comparisons were performed statistically. The inclusion of a new plugin was able to generate global values for Ecc and Ell that are in good agreement with previous findings using MRI. Despite no regional variation found in the peak Ecc, the peak Ell exhibited regional variation at the anterior side of the RV, which is potentially due to differences in biventricular torsion at the RV insertion point and fiber architecture. These results provide fundamental insights into the regional contractile function of the RV in healthy rat and could act as a normative baseline for future studies on regional changes induced by disease or treatment.

Assessment of the regional distribution of normalized circumferential strain in the thoracic and abdominal aorta using DENSE cardiovascular magnetic resonance

BACKGROUND

Displacement Encoding with Stimulated Echoes (DENSE) cardiovascular magnetic resonance (CMR) of the aortic wall offers the potential to improve patient-specific diagnostics and prognostics of diverse aortopathies by quantifying regionally heterogeneous aortic wall strain in vivo. However, before regional mapping of strain can be used to clinically assess aortic pathology, an evaluation of the natural variation of normal regional aortic kinematics is required.

METHOD

Aortic spiral cine DENSE CMR was performed at 3 T in 30 healthy adult subjects (range 18 to 65 years) at one or more axial locations that are at high risk for aortic aneurysm or dissection: the infrarenal abdominal aorta (IAA, n = 11), mid-descending thoracic aorta (DTA, n = 17), and/or distal aortic arch (DAA, n = 11). After implementing custom noise-reduction techniques, regional circumferential Green strain of the aortic wall was calculated across 16 sectors around the aortic circumference at each location and normalized by the mean circumferential strain for comparison between individuals.

RESULTS

The distribution of normalized circumferential strain (NCS) was heterogeneous for all locations evaluated. Despite large differences in mean strain between subjects, comparisons of NCS revealed consistent patterns of strain distribution for similar groupings of patients by axial location, age, and/or mean displacement angle. NCS at local systole was greatest in the lateral/posterolateral walls in the IAAs (1.47 ± 0.27), medial wall in anteriorly displacing DTAs (1.28 ± 0.20), lateral wall in posteriorly displacing DTAs (1.29 ± 0.29), superior curvature in DAAs < 50 years-old (1.93 ± 0.22), and medial wall in DAAs > 50 years (2.29 ± 0.58). The distribution of strain was strongly influenced by the location of the vertebra and other surrounding structures unique to each location.

CONCLUSIONS

Regional in vivo circumferential strain in the adult aorta is unique to each axial location and heterogeneous around its circumference, but can be grouped into consistent patterns defined by basic patient-specific metrics following normalization. The heterogeneous strain distributions unique to each group may be due to local peri-aortic constraints (particularly at the aorto-vertebral interface), heterogeneous material properties, and/or heterogeneous flow patterns. These results must be carefully considered in future studies seeking to clinically interpret or computationally model patient-specific aortic kinematics.

Can post-chemotherapy cardiotoxicity be detected in long-term survivors of breast cancer via comprehensive 3D left-ventricular contractility (strain) analysis?

PURPOSE

This study applied a novel and automated contractility analysis tool to investigate possible cardiotoxicity-related left-ventricular (LV) dysfunction in breast cancer patients following treatment with anti-neoplastic chemotherapy agents (CTA). Subclinical dysfunction otherwise undetected via LV ejection fraction (LVEF) was determined.

METHODS

Deformation data were acquired with the Displacement Encoding with Stimulated Echoes (DENSE) MRI sequence on 16 female patients who had CTA-based treatment. The contractility analysis tool consisting of image quantization-based boundary detection and the meshfree Radial Point Interpolation Method was used to compare chamber quantifications, 3D regional strains and torsion between patients and healthy subjects (N = 26 females with N = 14 age-matched). Quantifications of patient LVEFs from DENSE and Steady-State Free Precession (SSFP) acquisitions were compared, Bland-Altman interobserver agreements measured on their strain results and differences in contractile parameters with healthy subjects determined via Student's t-tests.

RESULTS

A significant difference was not found between DENSE and SSFP-based patient LVEFs at 58 ± 7% vs 57 ± 9%, p = 0.6. Bland-Altman agreements were - 0.01 ± 0.05 for longitudinal strain and 0.1 ± 1.3° for torsion. Differences in basal diameter indicating enlargement, 5.2 ± 0.5 cm vs 4.5 ± 0.5 cm, p < 0.01, and torsion, 4.7 ± 1.0° vs 8.1 ± 1.1°, p < 0.001 in the mid-ventricle and 5.9 ± 1.2° vs 10.2 ± 0.9°, p < 0.001 apically, were seen between patients and age-matched healthy subjects and similarly in longitudinal strain, but not in LVEF.

CONCLUSIONS

Results from the statistical analysis reveal the likelihood of LV remodeling in this patient subpopulation otherwise not indicated by LVEF measurements.

Cardiac Magnetic Resonance as a Tool to Assess Dyssynchrony

Cardiac resynchronization therapy (CRT) improves cardiac mechanics and quality of life in many patients with evidence of electromechanical cardiac dyssynchrony. However, up to 30% of patients receiving CRT do not respond to therapy. The mediator for poor response likely varies among patients; however, careful evaluation of mechanical dyssynchrony may inform management strategies. In this article, some of the methods and supporting evidence for dyssynchrony assessment with MRI as a predictor for CRT response are presented. The case is made for pre-implant assessment with MRI because of its ability to characterize scar, coronary venous distribution, and regional strain patterns.

Quantifying cardiac-induced brain tissue expansion using DENSE

Brain tissue undergoes viscoelastic deformation and volumetric strain as it expands over the cardiac cycle due to blood volume changes within the underlying microvasculature. Volumetric strain measurements may therefore provide insights into small vessel function and tissue viscoelastic properties. Displacement encoding via stimulated echoes (DENSE) is an MRI technique that can quantify the submillimetre displacements associated with brain tissue motion. Despite previous studies reporting brain tissue displacements using DENSE and other MRI techniques, a complete picture of brain tissue volumetric strain over the cardiac cycle has not yet been obtained. To address this need we implemented 3D cine-DENSE at 7 T and 3 T to investigate the feasibility of measuring cardiac-induced volumetric strain as a marker for small vessel blood volume changes. Volumetric strain over the entire cardiac cycle was computed for the whole brain and for grey and white matter tissue separately in six healthy human subjects. Signal-to-noise ratio (SNR) measurements were used to determine the voxel-wise volumetric strain noise. Mean peak whole brain volumetric strain at 7 T (mean ± SD) was (4.5 ± 1.0) × 10^-4 (corresponding to a volume expansion of 0.48 ± 0.1 mL), which is in agreement with literature values for cerebrospinal fluid that is displaced into the spinal canal to maintain a stable intracranial pressure. The peak volumetric strain ratio of grey to white matter was 4.4 ± 2.8, reflecting blood volume and tissue stiffness differences between these tissue types. The mean peak volumetric strains of grey and white matter tissue were found to be significantly different (p < 0.001). The mean SNR at 7 T and 3 T of the DENSE measurements was 22.0 ± 7.3 and 7.0 ± 2.8 respectively, which currently limits a voxel-wise strain analysis at both field strengths. We demonstrate that tissue specific quantification of volumetric strain is feasible with DENSE. This metric holds potential for studying blood volume pulsations in the ageing brain in healthy and diseased states.

The feasibility of a novel limited field of view spiral cine DENSE sequence to assess myocardial strain in dilated cardiomyopathy

Objective

Develop an accelerated cine displacement encoding with stimulated echoes (DENSE) cardiovascular magnetic resonance (CMR) sequence to enable clinically feasible myocardial strain evaluation in patients with dilated cardiomyopathy (DCM).

Materials and methods

A spiral cine DENSE sequence was modified by limiting the field of view in two dimensions using in-plane slice-selective pulses in the stimulated echo. This reduced breath hold duration from 20RR to 14RR intervals. Following phantom and pilot studies, the feasibility of the sequence to assess peak radial, circumferential, and longitudinal strain was tested in control subjects (n = 18) and then applied in DCM patients (n = 29).

Results

DENSE acquisition was possible in all participants. Elements of the data were not analysable in 1 control (6%) and 4 DCM r(14%) subjects due to off-resonance or susceptibility artefacts and low signal-to-noise ratio. Peak radial, circumferential, short-axis contour strain and longitudinal strain was reduced in DCM patients (p < 0.001 vs. controls) and strain measurements correlated with left ventricular ejection fraction (with circumferential strain r = − 0.79, p < 0.0001; with vertical long-axis strain r = − 0.76, p < 0.0001). All strain measurements had good inter-observer agreement (ICC > 0.80), except peak radial strain.

Discussion

We demonstrate the feasibility of CMR strain assessment in healthy controls and DCM patients using an accelerated cine DENSE technique. This may facilitate integration of strain assessment into routine CMR studies.

Electronic supplementary material

The online version of this article (10.1007/s10334-019-00735-5) contains supplementary material, which is available to authorized users.

Predictors of segmental myocardial functional recovery in patients after an acute ST-Elevation myocardial infarction

Objective

We hypothesized that Displacement Encoding with Stimulated Echoes (DENSE) and feature-tracking derived circumferential strain would provide incremental prognostic value over the extent of infarction for recovery of segmental myocardial function.

Methods

Two hundred and sixty-one patients (mean age 59 years, 73% male) underwent MRI 2 days post-ST elevation myocardial infarction (STEMI) and 241 (92%) underwent repeat imaging 6 months later.

The MRI protocol included cine, 2D-cine DENSE, T2 mapping and late enhancement.

Wall motion scoring was assessed by 2-blinded observers and adjudicated by a third. (WMS: 1=normal, 2=hypokinetic, 3=akinetic, 4=dyskinetic). WMS improvement was defined as a decrease in WMS ≥ 1, and normalization where WMS = 1 on follow-up. Segmental circumferential strain was derived utilizing DENSE and feature-tracking.

A generalized linear mixed model with random effect of subject was constructed and used to account for repeated sampling when investigating predictors of segmental myocardial improvement or normalization

Results

At baseline and follow-up, 1416 segments had evaluable data for all parameters. Circumferential strain by DENSE (p < 0.001) and feature-tracking (p < 0.001), extent of oedema (p < 0.001), infarct size (p < 0.001), and microvascular obstruction (p < 0.001) were associates of both improvement and normalization of WMS. Circumferential strain provided incremental predictive value even after accounting for infarct size, extent of oedema and microvascular obstruction, for segmental improvement (DENSE: odds ratio, 95% confidence intervals: 1.08 per −1% peak strain, 1.05–1.12, p < 0.001, feature-tracking: odds ratio, 95% confidence intervals: 1.05 per −1% peak strain, 1.03–1.07, p < 0.001) and segmental normalization (DENSE: 1.08 per −1% peak strain, 1.04–1.12, p < 0.001, feature-tracking: 1.06 per −1% peak strain, 1.04–1.08, p < 0.001).

Conclusions

Circumferential strain provides incremental prognostic value over segmental infarct size in patients post STEMI for predicting segmental improvement or normalization by wall-motion scoring.

Cardiac MRI biomarkers for Duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is a fatal inherited genetic disorder that results in progressive muscle weakness and ultimately loss of ambulation, respiratory failure and heart failure. Cardiac MRI (MRI) plays an increasingly important role in the diagnosis and clinical care of boys with DMD and associated cardiomyopathies. Conventional cardiac MRI biomarkers permit measurements of global cardiac function and presence of fibrosis, but changes in these measures are late manifestations. Emerging MRI biomarkers of myocardial function and structure include the estimation of rotational mechanics and regional strain using MRI tagging; T1-mapping; and T2-mapping, a marker of inflammation, edema and fat. These emerging biomarkers provide earlier insights into cardiac involvement in DMD, improving patient care and aiding the evaluation of emerging therapies.

Magnetic Resonance Elastography of Rodent Brain

Magnetic resonance elastography (MRE) is a non-invasive imaging technique, using the propagation of mechanical waves as a probe to palpate biological tissues. It consists in three main steps: production of shear waves within the tissue; encoding subsequent tissue displacement in magnetic resonance images; and extraction of mechanical parameters based on dedicated reconstruction methods. These three steps require an acoustic-frequency mechanical actuator, magnetic resonance imaging acquisition, and a post-processing tool for which no turnkey technology is available. The aim of the present review is to outline the state of the art of reported set-ups to investigate rodent brain mechanical properties. The impact of experimental conditions in dimensioning the set-up (wavelength and amplitude of the propagated wave, spatial resolution, and signal-to-noise ratio of the acquisition) on the accuracy and precision of the extracted parameters is discussed, as well as the influence of different imaging sequences, scanners, electromagnetic coils, and reconstruction algorithms. Finally, the performance of MRE in demonstrating viscoelastic differences between structures constituting the physiological rodent brain, and the changes in brain parameters under pathological conditions, are summarized. The recently established link between biomechanical properties of the brain as obtained on MRE and structural factors assessed by histology is also studied. This review intends to give an accessible outline on how to conduct an elastography experiment, and on the potential of the technique in providing valuable information for neuroscientists.

Free-breathing cine DENSE MRI using phase cycling with matchmaking and stimulated-echo image-based navigators

PURPOSE

This study aimed to develop a self-navigated method for free-breathing spiral cine displacement encoding with stimulated echoes (DENSE), a myocardial strain imaging technique that uses phase-cycling for artifact suppression. The method needed to address 2 consequences of motion for DENSE: striping artifacts from incomplete suppression of the T1 -relaxation echo and blurring.

METHODS

The method identifies phase-cycled spiral interleaves at matched respiratory phases by minimizing the residual signal due to T1 relaxation after phase-cycling subtraction. Next, the method reconstructs image-based navigators from matched phase-cycled interleaves that are comprised of the stimulated echo (ste-iNAVs). Ste-iNAVs are used for motion estimation and compensation of k-space data. The method was demonstrated in phantoms and compared to diaphragm-based navigator (dNAV) and conventional iNAV (c-iNAV) methods for the reconstruction of free-breathing volunteer data sets (N = 10).

RESULTS

Phantom experiments demonstrated that the proposed method removes striping artifacts and blurring due to motion. Volunteer results showed that respiratory motion measured by ste-iNAVs was better correlated than c-iNAVs to dNAV data (R2  = 0.82 ± 0.03 vs. 0.70 ± 0.05, P < 0.05). Match-making reconstructions of free-breathing data sets achieved lower residual T1 -relaxation echo energy (1.04 ± 0.01 vs. 1.18 ± 0.04 for dNAV and 1.18 ± 0.03 for c-iNAV, P < 0.05), higher apparent SNR (11.93 ± 1.05 vs. 10.68 ± 1.06 for dNAV and 10.66 ± 0.99 for c-iNAV, P < 0.05), and better phase quality (0.147 ± 0.012 vs. 0.166 ± 0.017 for dNAV, P = 0.06, and 0.168 ± 0.015 for c-iNAV, P < 0.05) than dNAV and c-iNAV methods.

CONCLUSION

For free-breathing cine DENSE, the proposed method addresses both types of breathing-induced artifacts and provides better quality images than conventional dNAV and iNAV methods.

3-D Measurements of Acceleration-Induced Brain Deformation via Harmonic Phase Analysis and Finite-Element Models

OBJECTIVE

To obtain dense spatiotemporal measurements of brain deformation from two distinct but complementary head motion experiments: linear and rotational accelerations.

METHODS

This study introduces a strategy for integrating harmonic phase analysis of tagged magnetic resonance imaging (MRI) and finite-element models to extract mechanically representative deformation measurements. The method was calibrated using simulated as well as experimental data, demonstrated in a phantom including data with image artifacts, and used to measure brain deformation in human volunteers undergoing rotational and linear acceleration.

RESULTS

Evaluation methods yielded a displacement error of 1.1 mm compared to human observers and strain errors between [Formula: see text] for linear acceleration and [Formula: see text] for rotational acceleration. This study also demonstrates an approach that can reduce error by 86% in the presence of corrupted data. Analysis of results shows consistency with 2-D motion estimation, agreement with external sensors, and the expected physical behavior of the brain.

CONCLUSION

Mechanical regularization is useful for obtaining dense spatiotemporal measurements of in vivo brain deformation under different loading regimes.

SIGNIFICANCE

The measurements suggest that the brain's 3-D response to mild accelerations includes distinct patterns observable using practical MRI resolutions. This type of measurement can provide validation data for computer models for the study of traumatic brain injury.

Demonstration of circumferential heterogeneity in displacement and strain in the abdominal aortic wall by spiral cine DENSE MRI

BACKGROUND

Knowledge of tissue properties of the abdominal aorta can improve understanding of vascular disease and guide interventional approaches. Existing MRI methods to quantify aortic wall displacement and strain are unable to discern circumferential heterogeneity.

PURPOSE

To assess regional variation in abdominal aortic wall displacement and strain as a function of circumferential position using spiral cine displacement encoding with stimulated echoes (DENSE).

STUDY TYPE

Prospective.

POPULATION

Cardiovascular disease-free men (n = 8) and women (n = 9) ages 30-42.

SEQUENCES

Prospective electrocardiogram (ECG)-gated and navigator echo-gated spiral, cine 2D DENSE and retrospective ECG-gated phase contrast MR (PCMR) sequences at 3T.

ASSESSMENT

In-plane displacement values of the aortic wall acquired with DENSE were used to determine radial and circumferential aortic wall motion. A quadrilateral-based 2D strain calculation method was implemented to determine strain from the displacement field. Peak displacement and its radial and circumferential contributions as well as peak circumferential strain were compared among eight circumferential wall segments. Distensibility was calculated using PCMR and compared with homogenized circumferential strain.

STATISTICAL TESTS

To account for repeated measurements in volunteers, linear mixed models for mean sector values were created for displacement magnitude, circumferential displacement, radial displacement, and circumferential strain. Comparisons were made between sectors. Calculated distensibility and homogenized circumferential strain were compared using Bland-Altman analysis. Statistical significance was defined as P < 0.05.

RESULTS

Displacement was highest in the anterior wall (1.5 ± 0.7 mm) and was primarily in the radial as compared with circumferential direction (1.04 ± 0.05 mm vs. 0.81 ± 0.42 mm). Circumferential strain was highest in the lateral walls (left 0.16 ± 0.05 and right 0.21 ± 0.12) with homogenized circumferential strain of 0.14 ± 0.05.

DATA CONCLUSION

DENSE imaging in the abdominal aortic wall demonstrated that the anterior aortic wall exhibits the greatest displacement, while the lateral wall experiences the largest circumferential strain.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;49:731-743.

Encoding and readout strategies in magnetic resonance elastography

Magnetic resonance elastography (MRE) has evolved significantly since its inception. Advances in motion-encoding gradient design and readout strategies have led to improved encoding and signal-to-noise ratio (SNR) efficiencies, which in turn allow for higher spatial resolution, increased coverage, and/or shorter scan times. The purpose of this review is to summarize MRE wave-encoding and readout approaches in a unified mathematical framework to allow for a comparative assessment of encoding and SNR efficiency of the various methods available. Besides standard full- and fractional-wave-encoding approaches, advanced techniques including flow compensation, sample interval modulation and multi-shot encoding are considered. Signal readout using fast k-space trajectories, reduced field of view, multi-slice, and undersampling techniques are summarized and put into perspective. The review is concluded with a foray into displacement and diffusion encoding as alternative and/or complementary techniques.

Comparison of left ventricular strains and torsion derived from feature tracking and DENSE CMR

BACKGROUND

Cardiovascular magnetic resonance (CMR) feature tracking is increasingly used to quantify cardiac mechanics from cine CMR imaging, although validation against reference standard techniques has been limited. Furthermore, studies have suggested that commonly-derived metrics, such as peak global strain (reported in 63% of feature tracking studies), can be quantified using contours from just two frames - end-diastole (ED) and end-systole (ES) - without requiring tracking software. We hypothesized that mechanics derived from feature tracking would not agree with those derived from a reference standard (displacement-encoding with stimulated echoes (DENSE) imaging), and that peak strain from feature tracking would agree with that derived using simple processing of only ED and ES contours.

METHODS

We retrospectively identified 88 participants with 186 pairs of DENSE and balanced steady state free precession (bSSFP) image slices acquired at the same locations across two institutions. Left ventricular (LV) strains, torsion, and dyssynchrony were quantified from both feature tracking (TomTec Imaging Systems, Circle Cardiovascular Imaging) and DENSE. Contour-based strains from bSSFP images were derived from ED and ES contours. Agreement was assessed with Bland-Altman analyses and coefficients of variation (CoV). All biases are reported in absolute percentage.

RESULTS

Comparison results were similar for both vendor packages (TomTec and Circle), and thus only TomTec Imaging System data are reported in the abstract for simplicity. Compared to DENSE, mid-ventricular circumferential strain (Ecc) from feature tracking had acceptable agreement (bias: - 0.4%, p = 0.36, CoV: 11%). However, feature tracking significantly overestimated the magnitude of Ecc at the base (bias: - 4.0% absolute, p < 0.001, CoV: 18%) and apex (bias: - 2.4% absolute, p = 0.01, CoV: 15%), underestimated torsion (bias: - 1.4 deg/cm, p < 0.001, CoV: 41%), and overestimated dyssynchrony (bias: 26 ms, p < 0.001, CoV: 76%). Longitudinal strain (Ell) had borderline-acceptable agreement (bias: - 0.2%, p = 0.77, CoV: 19%). Contour-based strains had excellent agreement with feature tracking (biases: - 1.3-0.2%, CoVs: 3-7%).

CONCLUSION

Compared to DENSE as a reference standard, feature tracking was inaccurate for quantification of apical and basal LV circumferential strains, longitudinal strain, torsion, and dyssynchrony. Feature tracking was only accurate for quantification of mid LV circumferential strain. Moreover, feature tracking is unnecessary for quantification of whole-slice strains (e.g. base, apex), since simplified processing of only ED and ES contours yields very similar results to those derived from feature tracking. Current feature tracking technology therefore has limited utility for quantification of cardiac mechanics.

Typical readout durations in spiral cine DENSE yield blurred images and underestimate cardiac strains at both 3.0 T and 1.5 T

INTRODUCTION

Displacement encoding with stimulated echoes (DENSE) is a phase contrast technique that encodes tissue displacement into phase images, which are typically processed into measures of cardiac function such as strains. For improved signal to noise ratio and spatiotemporal resolution, DENSE is often acquired with a spiral readout using an 11.1 ms readout duration. However, long spiral readout durations are prone to blurring due to common phenomena such as off-resonance and T2* decay, which may alter the resulting quantifications of strain. We hypothesized that longer readout durations would reduce image quality and underestimate cardiac strains at both 3.0 T and 1.5 T and that using short readout durations could overcome these limitations.

MATERIAL AND METHODS

Computational simulations were performed to investigate the relationship between off-resonance and T2* decay, the spiral cine DENSE readout duration, and measured radial and circumferential strain. Five healthy participants subsequently underwent 2D spiral cine DENSE at both 3.0 T and 1.5 T with several different readout durations 11.1 ms and shorter. Pearson correlations were used to assess the relationship between cardiac strains and the spiral readout duration.

RESULTS

Simulations demonstrated that long readout durations combined with off-resonance and T2* decay yield blurred images and underestimate strains. With the typical 11.1 ms DENSE readout, blurring was present in the anterior and lateral left ventricular segments of participants and was markedly improved with shorter readout durations. Radial and circumferential strains from those segments were significantly correlated with the readout duration. Compared to the 1.9 ms readout, the 11.1 ms readout underestimated radial and circumferential strains in those segments at both field strengths by up to 19.6% and 1.5% (absolute), or 42% and 7% (relative), respectively.

CONCLUSIONS

Blurring is present in spiral cine DENSE images acquired at both 3.0 T and 1.5 T using the typical 11.1 ms readout duration, which yielded substantially reduced radial strains and mildly reduced circumferential strains. Clinical studies using spiral cine DENSE should consider these limitations, while future technical advances may need to leverage accelerated techniques to improve the robustness and accuracy of the DENSE acquisition rather than focusing solely on reduced acquisition time.

Differential Effects of Isoproterenol on Regional Myocardial Mechanics in Rat using 3D cine DENSE Cardiovascular Magnetic Resonance

The present study assessed the acute effects of isoproterenol on left ventricular (LV) mechanics in healthy rats with the hypothesis that ß-adrenergic stimulation influences the mechanics of different myocardial regions of the LV wall in different ways. To accomplish this, magnetic resonance images were obtained in the LV of healthy rats with or without isoproterenol infusion. The LV contours were divided into basal, mid-ventricular, and apical regions. Additionally, the mid-ventricular myocardium was divided into three transmural layers with each layer partitioned into four segments (i.e., septal, inferior, lateral, and anterior). Peak systolic strains and torsion were quantified for each region. Isoproterenol significantly increased peak systolic radial strain and circumferential-longitudinal shear strain, as well as ventricular torsion, throughout the basal, mid-ventricle, and apical regions. In the mid-ventricle, isoproterenol significantly increased peak systolic radial strain, and induced significant increases in peak systolic circumferential strain and longitudinal strain in the septum. Isoproterenol consistently increased peak systolic circumferential-longitudinal shear strain in all mid-ventricular segments. Ventricular torsion was significantly increased in nearly all segments except the inferior sub-endocardium. The effects of isoproterenol on LV systolic mechanics (i.e., 3D strains and torsion) in healthy rats depend on the region. This region-dependency is also strain component-specific. These results provide insight into the regional response of LV mechanics to ß-adrenergic stimulation in rats, and could act as a baseline for future studies on subclinical abnormalities associated with the inotropic response in heart disease.

Regional Quantification of Brain Tissue Strain Using Displacement-Encoding With Stimulated Echoes Magnetic Resonance Imaging

Intrinsic cardiac-induced deformation of brain tissue is thought to be important in the pathophysiology of various neurological disorders. In this study, we evaluated the feasibility of utilizing displacement encoding with stimulated echoes (DENSE) magnetic resonance imaging (MRI) to quantify two-dimensional (2D) neural tissue strain using cardiac-driven brain pulsations. We examined eight adult healthy volunteers with an electrocardiogram-gated spiral DENSE sequence performed at the midsagittal plane on a 3 Tesla MRI scanner. Displacement, pixel-wise trajectories, and principal strains were determined in seven regions of interest (ROI): the brain stem, cerebellum, corpus callosum, and four cerebral lobes. Quantification of small neural tissue motion and strain along with their spatial and temporal variations in different brain regions was found to be feasible using DENSE. The medial and inferior brain structures (brain stem, cerebellum, and corpus callosum) had significantly larger motion and strain compared to structures located more peripherally. The brain stem had the largest peak mean displacement (PMD) (187 ± 50 μm, mean ± SD). The largest mean principal strains in compression and extension were observed in the brain stem (0.38 ± 0.08%) and the corpus callosum (0.37 ± 0.08%), respectively. Measured values in percent strain were altered by as much as 0.1 between repeated scans. This study showed that DENSE can quantify regional variations in brain tissue motion and strain and has the potential to be utilized as a tool to evaluate the changes in brain tissue dynamics resulting from alterations in biomechanical stresses and tissue properties.

In vivo quantification of regional circumferential Green strain in the thoracic and abdominal aorta by 2D spiral cine DENSE MRI

INTRODUCTION

Regional tissue mechanics play a fundamental role in patient-specific cardiovascular function. Nevertheless, regional assessments of aortic kinematics remain lacking due to the challenge of imaging the thin aortic wall. Herein, we present a novel application of DENSE (Displacement Encoding with Stimulated Echoes) MRI to quantify the circumferential Green strain of the thoracic and abdominal aorta.

METHODS

2D spiral cine DENSE and steady-state free procession (SSFP) cine images were acquired at 3T at the infrarenal aorta (IAA), descending thoracic aorta (DTA), or distal aortic arch (DAA) in a pilot study of 6 healthy volunteers. DENSE data was processed with multiple custom noise-reduction techniques to calculate circumferential Green strain across 16 equispaced sectors around the aorta. Each volunteer was scanned twice to evaluate interstudy repeatability.

RESULTS

Circumferential strain was heterogeneously distributed in all volunteers and locations. Spatial heterogeneity index by location was 0.37 (IAA), 0.28 (DTA), and 0.59 (DAA). Mean peak strain by DENSE for each cross-section was consistent with the homogenized linearized strain estimated from SSFP cine. The mean difference in peak strain across all sectors following repeat imaging was -0.1±2.2%, with a mean absolute difference of 1.7%.

CONCLUSIONS

Aortic cine DENSE MRI is a viable non-invasive technique for quantifying heterogeneous regional aortic wall strain and has significant potential to improve patient-specific clinical assessments of numerous aortopathies, as well as to provide the lacking spatiotemporal data required to refine computational models of aortic growth and remodeling.

Introduction to a mechanism for automated myocardium boundary detection with displacement encoding with stimulated echoes (DENSE)

Objective:

Displacement ENcoding with Stimulated Echoes (DENSE) is an MRI technique developed to encode phase related to myocardial tissue displacements, and the displacement information directly applied towards detecting left-ventricular (LV) myocardial motion during the cardiac cycle. The purpose of this study is to present a novel, three-dimensional (3D) DENSE displacement-based and magnitude image quantization-based, semi-automated detection technique for myocardial wall motion, whose boundaries are used for rapid and automated computation of 3D myocardial strain.

Methods:

The architecture of this boundary detection algorithm is primarily based on pixelwise spatiotemporal increments in LV tissue displacements during the cardiac cycle and further reinforced by radially searching for pixel-based image gradients in multithreshold quantized magnitude images. This spatiotemporal edge detection methodology was applied to all LV partitions and their subsequent timeframes that lead to full 3D LV reconstructions. It was followed by quantifications of 3D chamber dimensions and myocardial strains, whose rapid computation was the primary motivation behind developing this algorithm. A pre-existing two-dimensional (2D) semi-automated contouring technique was used in parallel to validate the accuracy of the algorithm and both methods tested on DENSE data acquired in (N = 14) healthy subjects. Chamber quantifications between methods were compared using paired t-tests and Bland–Altman analysis established regional strain agreements.

Results:

There were no significant differences in the results of chamber quantifications between the 3D semi-automated and existing 2D boundary detection techniques. This included comparisons of ejection fractions, which were 0.62 ± 0.04 vs 0.60 ± 0.06 (p = 0.23) for apical, 0.60 ± 0.04 vs 0.59 ± 0.05 (p = 0.76) for midventricular and 0.56 ± 0.04 vs 0.58 ± 0.05 (p = 0.07) for basal segments, that were quantified using the 3D semi-automated and 2D pre-existing methodologies, respectively. Bland–Altman agreement between regional strains generated biases of 0.01 ± 0.06, –0.01 ± 0.01 and 0.0 ± 0.06 for the radial, circumferential and longitudinal directions, respectively.

Conclusion:

A new, 3D semi-automated methodology for contouring the entire LV and rapidly generating chamber quantifications and regional strains is presented that was validated in relation to an existing 2D contouring technique.

Advances in knowledge:

This study introduced a scientific tool for rapid, semi-automated generation of clinical information regarding shape and function in the 3D LV.

Evaluation of a Novel Finite Element Model of Active Contraction in the Heart

Finite element (FE) modeling is becoming a widely used approach for the investigation of global heart function. In the present study, a novel model of cellular-level systolic contraction, which includes both length- and velocity-dependence, was implemented into a 3D non-linear FE code. To validate this new FE implementation, an optimization procedure was used to determine the contractile parameters, associated with sarcomeric function, by comparing FE-predicted pressure and strain to experimental measures collected with magnetic resonance imaging and catheterization in the ventricles of five healthy rats. The pressure-volume relationship generated by the FE models matched well with the experimental data. Additionally, the regional distribution of end-systolic strains and circumferential-longitudinal shear angle exhibited good agreement with experimental results overall, with the main deviation occurring in the septal region. Moreover, the FE model predicted a heterogeneous distribution of sarcomere re-lengthening after ventricular ejection, which is consistent with previous in vivo studies. In conclusion, the new FE active contraction model was able to predict the global performance and regional mechanical behaviors of the LV during the entire cardiac cycle. By including more accurate cellular-level mechanisms, this model could provide a better representation of the LV and enhance cardiac research related to both systolic and diastolic dysfunction.

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

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).

An in-vivo comparison of stimulated-echo and motion compensated spin-echo sequences for 3 T diffusion tensor cardiovascular magnetic resonance at multiple cardiac phases

BACKGROUND

Stimulated-echo (STEAM) and, more recently, motion-compensated spin-echo (M2-SE) techniques have been used for in-vivo diffusion tensor cardiovascular magnetic resonance (DT-CMR) assessment of cardiac microstructure. The two techniques differ in the length scales of diffusion interrogated, their signal-to-noise ratio efficiency and sensitivity to both motion and strain. Previous comparisons of the techniques have used high performance gradients at 1.5 T in a single cardiac phase. However, recent work using STEAM has demonstrated novel findings of microscopic dysfunction in cardiomyopathy patients, when DT-CMR was performed at multiple cardiac phases. We compare STEAM and M2-SE using a clinical 3 T scanner in three potentially clinically interesting cardiac phases.

METHODS

Breath hold mid-ventricular short-axis DT-CMR was performed in 15 subjects using M2-SE and STEAM at end-systole, systolic sweet-spot and diastasis. Success was defined by ≥50% of the myocardium demonstrating normal helix angles. From successful acquisitions DT-CMR results relating to tensor orientation, size and shape were compared between sequences and cardiac phases using non-parametric statistics. Strain information was obtained using cine spiral displacement encoding with stimulated echoes for comparison with DT-CMR results.

RESULTS

Acquisitions were successful in 98% of STEAM and 76% of M2-SE cases and visual helix angle (HA) map scores were higher for STEAM at the sweet-spot and diastasis. There were significant differences between sequences (p < 0.05) in mean diffusivity (MD), fractional anisotropy (FA), tensor mode, transmural HA gradient and absolute second eigenvector angle (E2A). Differences in E2A between systole and diastole correlated with peak radial strain for both sequences (p ≤ 0.01).

CONCLUSION

M2-SE and STEAM can be performed equally well at peak systole at 3 T using standard gradients, but at the sweet-spot and diastole STEAM is more reliable and image quality scores are higher. Differences in DT-CMR results are potentially due to differences in motion sensitivity and the longer diffusion time of STEAM, although the latter appears to be the dominant factor. The benefits of both sequences should be considered when planning future studies and sequence and cardiac phase specific normal ranges should be used for comparison.

Association between left ventricular mechanics and diffuse myocardial fibrosis in patients with repaired Tetralogy of Fallot: a cross-sectional study

BACKGROUND

Patients with repaired tetralogy of Fallot (TOF) have progressive, adverse biventricular remodeling, leading to abnormal contractile mechanics. Defining the mechanisms underlying this dysfunction, such as diffuse myocardial fibrosis, may provide insights into poor long-term outcomes. We hypothesized that left ventricular (LV) diffuse fibrosis is related to impaired LV mechanics.

METHODS

Patients with TOF were evaluated with cardiac magnetic resonance in which modified Look-Locker (MOLLI) T1-mapping and spiral cine Displacement encoding (DENSE) sequences were acquired at three LV short-axis positions. Linear mixed modeling was used to define the association between regional LV mechanics from DENSE based on regional T1-derived diffuse fibrosis measures, such as extracellular volume fraction (ECV).

RESULTS

Forty patients (26 ± 11 years) were included. LV ECV was generally within normal range (0.24 ± 0.05). For LV mechanics, peak circumferential strains (-15 ± 3%) and dyssynchrony indices (16 ± 8 ms) were moderately impaired, while peak radial strains (29 ± 8%) were generally normal. After adjusting for patient age, sex, and regional LV differences, ECV was associated with log-adjusted LV dyssynchrony index (β = 0.67) and peak LV radial strain (β = -0.36), but not LV circumferential strain. Moreover, post-contrast T1 was associated with log-adjusted LV diastolic circumferential strain rate (β = 0.37).

CONCLUSIONS

We observed several moderate associations between measures of fibrosis and impaired mechanics, particularly the LV dyssynchrony index and peak radial strain. Diffuse fibrosis may therefore be a causal factor in some ventricular dysfunction in TOF.

Regional quantification of myocardial mechanics in rat using 3D cine DENSE cardiovascular magnetic resonance

Rat models have assumed an increasingly important role in cardiac research. However, a detailed profile of regional cardiac mechanics, such as strains and torsion, is lacking for rats. We hypothesized that healthy rat left ventricles (LVs) exhibit regional differences in cardiac mechanics, which are part of normal function. In this study, images of the LV were obtained with 3D cine displacement encoding with stimulated echoes (DENSE) cardiovascular magnetic resonance in 10 healthy rats. To evaluate regional cardiac mechanics, the LV was divided into basal, mid-ventricular, and apical regions. The myocardium at the mid-LV was further partitioned into four wall segments (i.e. septal, inferior, lateral, and anterior) and three transmural layers (i.e. sub-endocardium, mid-myocardium, and sub-epicardium). The six Lagrangian strain components (i.e. Err , Ecc , Ell , Ecl , Erl , and Ecr ) were computed from the 3D displacement field and averaged within each region of interest. Torsion was quantified using the circumferential-longitudinal shear angle. While peak systolic Ecl differed between the mid-ventricle and apex, the other five components of peak systolic strain were similar across the base, mid-ventricle, and apex. In the mid-LV myocardium, Ecc decreased gradually from the sub-endocardial to the sub-epicardial layer. Ell demonstrated significant differences between the four wall segments, with the largest magnitude in the inferior segment. Err was uniform among the four wall segments. Ecl varied along the transmural direction and among wall segments, whereas Erl differed only among the wall segments. Erc was not associated with significant variations. Torsion also varied along the transmural direction and among wall segments. These results provide fundamental insights into the regional contractile function of healthy rat hearts, and form the foundation for future studies on regional changes induced by disease or treatments.

Magnetic Resonance Imaging of Myocardial Strain After Acute ST-Segment-Elevation Myocardial Infarction: A Systematic Review

The purpose of this systematic review is to provide a clinically relevant, disease-based perspective on myocardial strain imaging in patients with acute myocardial infarction or stable ischemic heart disease. Cardiac magnetic resonance imaging uniquely integrates myocardial function with pathology. Therefore, this review focuses on strain imaging with cardiac magnetic resonance. We have specifically considered the relationships between left ventricular (LV) strain, infarct pathologies, and their associations with prognosis. A comprehensive literature review was conducted in accordance with the PRISMA guidelines. Publications were identified that (1) described the relationship between strain and infarct pathologies, (2) assessed the relationship between strain and subsequent LV outcomes, and (3) assessed the relationship between strain and health outcomes. In patients with acute myocardial infarction, circumferential strain predicts the recovery of LV systolic function in the longer term. The prognostic value of longitudinal strain is less certain. Strain differentiates between infarcted versus noninfarcted myocardium, even in patients with stable ischemic heart disease with preserved LV ejection fraction. Strain recovery is impaired in infarcted segments with intramyocardial hemorrhage or microvascular obstruction. There are practical limitations to measuring strain with cardiac magnetic resonance in the acute setting, and knowledge gaps, including the lack of data showing incremental value in clinical practice. Critically, studies of cardiac magnetic resonance strain imaging in patients with ischemic heart disease have been limited by sample size and design. Strain imaging has potential as a tool to assess for early or subclinical changes in LV function, and strain is now being included as a surrogate measure of outcome in therapeutic trials.

Evaluation of the impact of strain correction on the orientation of cardiac diffusion tensors with in vivo and ex vivo porcine hearts

PURPOSE

To evaluate the importance of strain-correcting stimulated echo acquisition mode echo-planar imaging cardiac diffusion tensor imaging.

METHODS

Healthy pigs (n = 11) were successfully scanned with a 3D cine displacement-encoded imaging with stimulated echoes and a monopolar-stimulated echo-planar imaging diffusion tensor imaging sequence at 3 T during diastasis, peak systole, and strain sweet spots in a midventricular short-axis slice. The same diffusion tensor imaging sequence was repeated ex vivo after arresting the hearts in either a relaxed (KCl-induced) or contracted (BaCl₂ -induced) state. The displacement-encoded imaging with stimulated echoes data were used to strain-correct the in vivo cardiac diffusion tensor imaging in diastole and systole. The orientation of the primary (helix angles) and secondary (E2A) diffusion eigenvectors was compared with and without strain correction and to the strain-free ex vivo data.

RESULTS

Strain correction reduces systolic E2A significantly when compared without strain correction and ex vivo (median absolute E2A = 34.3° versus E2A = 57.1° (P = 0.01), E2A = 60.5° (P = 0.006), respectively). The systolic distribution of E2A without strain correction is closer to the contracted ex vivo distribution than with strain correction, root mean square deviation of 0.027 versus 0.038.

CONCLUSIONS

The current strain-correction model amplifies the contribution of microscopic strain to diffusion resulting in an overcorrection of E2A. Results show that a new model that considers cellular rearrangement is required. Magn Reson Med 79:2205-2215, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Impaired right ventricular contractile function in childhood obesity and its association with right and left ventricular changes: a cine DENSE cardiac magnetic resonance study

BACKGROUND

Pediatric obesity is a growing public health problem, which is associated with increased risk of cardiovascular disease and premature death. Left ventricular (LV) remodeling (increased myocardial mass and thickness) and contractile dysfunction (impaired longitudinal strain) have been documented in obese children, but little attention has been paid to the right ventricle (RV). We hypothesized that obese/overweight children would have evidence of RV remodeling and contractile dysfunction.

METHODS

One hundred and three children, ages 8-18 years, were prospectively recruited and underwent cardiovascular magnetic resonance (CMR), including both standard cine imaging and displacement encoding with stimulated echoes (DENSE) imaging, which allowed for quantification of RV geometry and function/mechanics. RV free wall longitudinal strain was quantified from the end-systolic four-chamber DENSE image. Linear regression was used to quantify correlations of RV strain with LV strain and measurements of body composition (adjusted for sex and height). Analysis of variance was used to study the relationship between RV strain and LV remodeling types (concentric remodeling, eccentric/concentric hypertrophy).

RESULTS

The RV was sufficiently visualized with DENSE in 70 (68%) subjects, comprising 36 healthy weight (13.6 ± 2.7 years) and 34 (12.1 ± 2.9 years) obese/overweight children. Obese/overweight children had a 22% larger RV mass index (8.2 ± 0.9 vs 6.7 ± 1.1 g/m2.7, p < 0.001) compared to healthy controls. RV free wall longitudinal strain was impaired in obese/overweight children (-16 ± 4% vs -19 ± 5%, p = 0.02). Ten (14%) out of 70 children had LV concentric hypertrophy, and these children had the most impaired RV longitudinal strain compared to those with normal LV geometry (-13 ± 4% vs -19 ± 5%, p = 0.002). RV longitudinal strain was correlated with LV longitudinal strain (r = 0.34, p = 0.004), systolic blood pressure (r = 0.33, p = 0.006), as well as BMI z-score (r = 0.28, p = 0.02), waist (r = 0.31, p = 0.01), hip (r = 0.40, p = 0.004) and abdominal (r = 0.38, p = 0.002) circumference, height and sex adjusted.

CONCLUSIONS

Obese/overweight children have evidence of RV remodeling (increased RV mass) and RV contractile dysfunction (impaired free wall longitudinal strain). Moreover, RV longitudinal strain correlates with LV longitudinal strain, and children with LV concentric hypertrophy show the most impaired RV function. These results suggest there may be a common mechanism underlying both remodeling and dysfunction of the left and right ventricles in obese/overweight children.

Microstructurally Anchored Cardiac Kinematics by Combining In Vivo DENSE MRI and cDTI

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.

Magnetic resonance imaging of myocardial strain: A review of current approaches

Contraction of the heart is central to its purpose of pumping blood around the body. While simple global function measures (such as the ejection fraction) are most commonly used in the clinical assessment of cardiac function, MRI also provides a range of approaches for quantitatively characterizing regional cardiac function, including the local deformation (or strain) within the heart wall. While they have been around for some years, these methods are still undergoing further technical development, and they have had relatively little clinical evaluation. However, they can provide potentially useful new ways to assess cardiac function, which may be able to contribute to better classification and treatment of heart disease. This article provides some basic background on the physical and physiological factors that determine the motion of the heart, in health and disease and then reviews some of the ways that MRI methods are being developed to image and quantify strain within the myocardium.

LEVEL OF EVIDENCE

4 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2017;46:1263-1280.

Motion Estimation with Finite-Element Biomechanical Models and Tracking Constraints from Tagged MRI

Noninvasive measurements of tissue deformation provide bio-mechanical insights of an organ, which can be used as clinical functional biomarkers or experimental data for validating computational simulations. However, acquisition of 3D displacement information is susceptible to experimental inconsistency and limited scan time. In this research, we describe the process of tracking tagged magnetic resonance imaging (MRI) as enforcing harmonic phase conservation in finite-element (FE) models. This concept is demonstrated as a tool for motion estimation in a brain motion phantom, the heart, and the tongue. Our results demonstrate that the new methodology offers robustness to edge and large-displacement artifacts, and that it can be seamlessly coupled with numerical simulations for estimating fiber stretch in residually stressed tissue, or for inverse identification of muscle activation.

Right Ventricular Strain, Torsion, and Dyssynchrony in Healthy Subjects Using 3D Spiral Cine DENSE Magnetic Resonance Imaging

Mechanics of the left ventricle (LV) are important indicators of cardiac function. The role of right ventricular (RV) mechanics is largely unknown due to the technical limitations of imaging its thin wall and complex geometry and motion. By combining 3D Displacement Encoding with Stimulated Echoes (DENSE) with a post-processing pipeline that includes a local coordinate system, it is possible to quantify RV strain, torsion, and synchrony. In this study, we sought to characterize RV mechanics in 50 healthy individuals and compare these values to their LV counterparts. For each cardiac frame, 3D displacements were fit to continuous and differentiable radial basis functions, allowing for the computation of the 3D Cartesian Lagrangian strain tensor at any myocardial point. The geometry of the RV was extracted via a surface fit to manually delineated endocardial contours. Throughout the RV, a local coordinate system was used to transform from a Cartesian strain tensor to a polar strain tensor. It was then possible to compute peak RV torsion as well as peak longitudinal and circumferential strain. A comparable analysis was performed for the LV. Dyssynchrony was computed from the standard deviation of regional activation times. Global circumferential strain was comparable between the RV and LV (-18.0% for both) while longitudinal strain was greater in the RV (-18.1% vs. -15.7%). RV torsion was comparable to LV torsion (6.2 vs. 7.1 degrees, respectively). Regional activation times indicated that the RV contracted later but more synchronously than the LV. 3D spiral cine DENSE combined with a post-processing pipeline that includes a local coordinate system can resolve both the complex geometry and 3D motion of the RV.

Multiparametric CMR imaging of infarct remodeling in a percutaneous reperfused Yucatan mini-pig model

To further understanding of the temporal evolution and pathophysiology of adverse ventricular remodeling over the first 60 days following a myocardial infarction (MI) in both the infarcted and remote myocardium, we performed multi-parametric cardiac magnetic resonance (CMR) imaging in a closed-chest chronic Yucatan mini-pig model of reperfused MI. Ten animals underwent 90 min left anterior descending artery occlusion and reperfusion. Three animals served as controls. Multiparametric CMR (1.5T) was performed at baseline, Day 2, Day 30 and in four animals on Day 60 after MI. Left ventricular (LV) volumes and infarct size were measured. T₁ and T₂ mapping sequences were performed to measure values in the infarct and remote regions. Remote region collagen fractions were compared between infarcted animals and controls. Procedure success was 80%. The model created large infarcts (28 ± 5% of LV mass on Day 2), which led to significant adverse myocardial remodeling that stabilized beyond 30 days. Native T₁ values did not reliably differentiate remote and infarct regions acutely. There was no evidence of remote fibrosis as indicated by partition coefficient and collagen fraction analyses. The infarct T₂ values remained elevated up to 60 days after MI. Multiparametric CMR in this model showed significant adverse ventricular remodeling 30 days after MI similar to that seen in humans. In addition, this study demonstrated that remote fibrosis is absent and that infarct T₂ signal remains chronically elevated in this model. These findings need to be considered when designing preclinical trials using CMR endpoints.

Using a respiratory navigator significantly reduces variability when quantifying left ventricular torsion with cardiovascular magnetic resonance

BACKGROUND

Left ventricular (LV) torsion is an important indicator of cardiac function that is limited by high inter-test variability (50% of the mean value). We hypothesized that this high inter-test variability is partly due to inconsistent breath-hold positions during serial image acquisitions, which could be significantly improved by using a respiratory navigator for cardiovascular magnetic resonance (CMR) based quantification of LV torsion.

METHODS

We assessed respiratory-related variability in measured LV torsion with two distinct experimental protocols. First, 17 volunteers were recruited for CMR with cine displacement encoding with stimulated echoes (DENSE) in which a respiratory navigator was used to measure and then enforce variability in end-expiratory position between all LV basal and apical acquisitions. From these data, we quantified the inter-test variability of torsion in the absence and presence of enforced end-expiratory position variability, which established an upper bound for the expected torsion variability. For the second experiment (in 20 new, healthy volunteers), 10 pairs of cine DENSE basal and apical images were each acquired from consecutive breath-holds and consecutive navigator-gated scans (with a single acceptance position). Inter-test variability of torsion was compared between the breath-hold and navigator-gated scans to quantify the variability due to natural breath-hold variation. To demonstrate the importance of these variability reductions, we quantified the reduction in sample size required to detect a clinically meaningful change in LV torsion with the use of a respiratory navigator.

RESULTS

The mean torsion was 3.4 ± 0.2°/cm. From the first experiment, enforced variability in end-expiratory position translated to considerable variability in measured torsion (0.56 ± 0.34°/cm), whereas inter-test variability with consistent end-expiratory position was 57% lower (0.24 ± 0.16°/cm, p < 0.001). From the second experiment, natural respiratory variability from consecutive breath-holds translated to a variability in torsion of 0.24 ± 0.10°/cm, which was significantly higher than the variability from navigator-gated scans (0.18 ± 0.06°/cm, p = 0.02). By using a respiratory navigator with DENSE, theoretical sample sizes were reduced from 66 to 16 and 26 to 15 as calculated from the two experiments.

CONCLUSIONS

A substantial portion (22-57%) of the inter-test variability of LV torsion can be reduced by using a respiratory navigator to ensure a consistent breath-hold position between image acquisitions.

Imaging left-ventricular mechanical activation in heart failure patients using cine DENSE MRI: Validation and implications for cardiac resynchronization therapy

PURPOSE

To image late mechanical activation and identify effective left-ventricular (LV) pacing sites for cardiac resynchronization therapy (CRT). There is variability in defining mechanical activation time, with some studies using the time to peak strain (TPS) and some using the time to the onset of circumferential shortening (TOS). We developed improved methods for imaging mechanical activation and evaluated them in heart failure (HF) patients undergoing CRT.

MATERIALS AND METHODS

We applied active contours to cine displacement encoding with stimulated echoes (DENSE) strain images to detect TOS. Six healthy volunteers underwent magnetic resonance imaging (MRI) at 1.5T, and 50 patients underwent pre-CRT MRI (strain, scar, volumes) and echocardiography, assessment of the electrical activation time (Q-LV) at the LV pacing site, and echocardiography assessment of LV reverse remodeling 6 months after CRT. TPS at the LV pacing site was also measured by DENSE.

RESULTS

The latest TOS was greater in HF patients vs. healthy subjects (112 ± 28 msec vs. 61 ± 7 msec, P < 0.01). The correlation between TOS and Q-LV was strong (r > 0.75; P < 0.001) and better than between TPS and Q-LV (r < 0.62; P ≥ 0.006). Twenty-three of 50 patients had the latest activating segment in a region other than the mid-ventricular lateral wall, the most common site for the CRT LV lead. Using a multivariable model, TOS/QRS was significantly associated with LV reverse remodeling even after adjustment for overall dyssynchrony and scar (P < 0.05), whereas TPS was not (P = 0.49).

CONCLUSION

Late activation by cine DENSE TOS analysis is associated with improved LV reverse remodeling with CRT and deserves further study as a tool to achieve optimal LV lead placement in CRT.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2017;46:887-896.

A Novel Method for Estimating Myocardial Strain: Assessment of Deformation Tracking Against Reference Magnetic Resonance Methods in Healthy Volunteers

We developed a novel method for tracking myocardial deformation using cardiac magnetic resonance (CMR) cine imaging. We hypothesised that circumferential strain using deformation-tracking has comparable diagnostic performance to a validated method (Displacement Encoding with Stimulated Echoes- DENSE) and potentially diagnostically superior to an established cine-strain method (feature-tracking). 81 healthy adults (44.6 ± 17.7 years old, 47% male), without any history of cardiovascular disease, underwent CMR at 1.5 T including cine, DENSE, and late gadolinium enhancement in subjects >45 years. Acquisitions were divided into 6 segments, and global and segmental peak circumferential strain were derived and analysed by age and sex. Peak circumferential strain differed between the 3 groups (DENSE: −19.4 ± 4.8%; deformation-tracking: −16.8 ± 2.4%; feature-tracking: −28.7 ± 4.8%) (ANOVA with Tukey post-hoc, F-value 279.93, p < 0.01). DENSE and deformation-tracking had better reproducibility than feature-tracking. Intra-class correlation co-efficient was >0.90. Larger magnitudes of strain were detected in women using deformation-tracking and DENSE, but not feature-tracking. Compared with a reference method (DENSE), deformation-tracking using cine imaging has similar diagnostic performance for circumferential strain assessment in healthy individuals. Deformation-tracking could potentially obviate the need for bespoke strain sequences, reducing scanning time and is more reproducible than feature-tracking.

Preliminary investigation of multiparametric strain Z-score (MPZS) computation using displacement encoding with simulated echoes (DENSE) and radial point interpretation method (RPIM)

PURPOSE

To describe and assess an automated normalization method for identifying sentinel (septal) regions of myocardial dysfunction in nonischemic, nonvalvular dilated cardiomyopathy (DCM), using an unprecedented combination of the navigator-gated 3D spiral displacement encoding with stimulated echoes (DENSE) magnetic resonance imaging (MRI), radial point interpolation (RPIM) and multiparametric strain z-score (MPZS).

MATERIALS AND METHODS

Navigator-gated 3D spiral DENSE, in a 1.5T MRI machine, was used for acquiring the displacement encoded complex images, MR Analytical Software System (MASS) for automated boundary detection and automated meshfree RPIM for left-ventricular (LV) myocardial strain computation to analyze MPZS in 36 subjects (with n = 17 DCM patients). Pearson's r correlation established relations between global/sentinel MPZS and ejection fraction (EF). The time taken for combined RPIM-MPZS computations was recorded.

RESULTS

Maximum MPZS differences were seen between anteroseptal and posterolateral regions in the base (2.0 ± 0.3 vs. 0.9 ± 0.5) and the mid-wall (2.1 ± 0.4 vs. 1.0 ± 0.4). These regional differences were found to be consistent with historically documented septal injury in nonischemic DCM. Correlations were 0.6 between global MPZS and EF, and 0.7 between sentinel MPZS and EF. The time taken for combined RPIM-MPZS computations per subject was 18.9 ± 5.9 seconds.

CONCLUSION

Heterogeneous contractility found in the sentinel regions with the current automated MPZS computation scheme and the correlation found between MPZS and EF may lead to the creation of a new clinical metric in LV DCM surveillance. J. MAGN. RESON. IMAGING 2016;44:993-1002.

Cardiovascular magnetic resonance detects the progression of impaired myocardial perfusion reserve and increased left-ventricular mass in mice fed a high-fat diet

BACKGROUND

Impaired myocardial perfusion reserve (MPR) is prevalent in obesity and diabetes, even in the absence of obstructive coronary artery disease (CAD), and is prognostic of adverse events. We sought to establish the time course of reduced MPR and to investigate associated vascular and tissue properties in mice fed a high-fat diet (HFD), as they are an emerging model of human obesity, diabetes, and reduced MPR without obstructive CAD.

METHODS

C57Bl/6 mice fed a HFD or a low-fat diet (control) were imaged at 6, 12, 18 and 24 weeks post-diet. The cardiovascular magnetic resonance (CMR) protocol included multi-slice cine imaging to assess ejection fraction (EF), left-ventricular (LV) mass, LV wall thickness (LVWT), and LV volumes, and first-pass perfusion CMR to quantify MPR. Coronary vascular reactivity, aortic atherosclerosis, myocardial capillary density and tissue fibrosis were also assessed.

RESULTS

Body weight was increased in HFD mice at 6-24 weeks post-diet (p < 0.05 vs. control). MPR in HFD mice was reduced and LV mass and LVWT were increased in HFD mice at 18 and 24 weeks post-diet (p < 0.05 vs. control). Coronary arteriolar vascular reactivity to adenosine and acetylcholine were reduced in HFD mice (p < 0.05 vs. control). There were no significant differences in cardiac volumes, EF, or capillary density measurements between the two groups. Histology showed interstitial fibrosis in HFD and no aortic atherosclerosis in either group.

CONCLUSIONS

C57Bl/6 mice fed a HFD for 18-24 weeks have progressively increased LV mass and impaired MPR with fibrosis, normal capillary density and no aortic plaque. These results establish C57Bl/6 mice fed a HFD for 18-24 weeks as a model of impaired MPR without obstructive CAD due to obesity and diabetes.

An interactive videogame designed to improve respiratory navigator efficiency in children undergoing cardiovascular magnetic resonance

BACKGROUND

Advanced cardiovascular magnetic resonance (CMR) acquisitions often require long scan durations that necessitate respiratory navigator gating. The tradeoff of navigator gating is reduced scan efficiency, particularly when the patient's breathing patterns are inconsistent, as is commonly seen in children. We hypothesized that engaging pediatric participants with a navigator-controlled videogame to help control breathing patterns would improve navigator efficiency and maintain image quality.

METHODS

We developed custom software that processed the Siemens respiratory navigator image in real-time during CMR and represented diaphragm position using a cartoon avatar, which was projected to the participant in the scanner as visual feedback. The game incentivized children to breathe such that the avatar was positioned within the navigator acceptance window (±3 mm) throughout image acquisition. Using a 3T Siemens Tim Trio, 50 children (Age: 14 ± 3 years, 48 % female) with no significant past medical history underwent a respiratory navigator-gated 2D spiral cine displacement encoding with stimulated echoes (DENSE) CMR acquisition first with no feedback (NF) and then with the feedback game (FG). Thirty of the 50 children were randomized to undergo extensive off-scanner training with the FG using a MRI simulator, or no off-scanner training. Navigator efficiency, signal-to-noise ratio (SNR), and global left-ventricular strains were determined for each participant and compared.

RESULTS

Using the FG improved average navigator efficiency from 33 ± 15 to 58 ± 13 % (p < 0.001) and improved SNR by 5 % (p = 0.01) compared to acquisitions with NF. There was no difference in navigator efficiency (p = 0.90) or SNR (p = 0.77) between untrained and trained participants for FG acquisitions. Circumferential and radial strains derived from FG acquisitions were slightly reduced compared to NF acquisitions (-16 ± 2 % vs -17 ± 2 %, p < 0.001; 40 ± 10 % vs 44 ± 11 %, p = 0.005, respectively). There were no differences in longitudinal strain (p = 0.38).

CONCLUSIONS

Use of a respiratory navigator feedback game during navigator-gated CMR improved navigator efficiency in children from 33 to 58 %. This improved efficiency was associated with a 5 % increase in SNR for spiral cine DENSE. Extensive off-scanner training was not required to achieve the improvement in navigator efficiency.

Reproducibility of cine displacement encoding with stimulated echoes (DENSE) in human subjects

OBJECTIVES

To test the hypothesis that two-dimensional (2D) displacement encoding via stimulated echoes (DENSE) is a reproducible technique for the depiction of segmental myocardial motion in human subjects.

MATERIALS AND METHODS

Following the approval of the institutional review board (IRB), 17 healthy volunteers without documented history of cardiovascular disease were recruited. For each participant, 2D DENSE were performed twice (at different days) and the images were obtained at basal, midventricular and apical levels of the left ventricle (LV) with a short-axis view. The radial thickening strain (Err), circumferential strain (Ecc), twist and torsion were calculated. The intra-, inter-observer and inter-study variations of DENSE-derived myocardial motion indices were evaluated using coefficient of variation (CoV) and intra-class correlation coefficient (ICC).

RESULTS

In total, there are 272 pairs of myocardial segments (data points) for comparison. There is good intra- and inter-observer reproducibility for all DENSE-derived measures in 17 participants. There is good inter-study reproducibility for peak Ecc (CoV=19.64%, ICC=0.8896, p<0.001), twist (CoV=33.11%, ICC=0.9135, p<0.001) and torsion (CoV=13.96%, ICC=0.8684, p<0.001). There is moderate inter-study reproducibility for Err (CoV=38.89%, ICC=0.7022, p<0.001).

CONCLUSION

DENSE is a reproducible technique for characterizing LV regional systolic myocardial motion on a per-segment basis in healthy volunteers.

Optimal configuration of respiratory navigator gating for the quantification of left ventricular strain using spiral cine displacement encoding with stimulated echoes (DENSE) MRI

PURPOSE

To determine the optimal respiratory navigator gating configuration for the quantification of left ventricular strain using spiral cine displacement encoding with stimulated echoes (DENSE) MRI.

MATERIALS AND METHODS

Two-dimensional spiral cine DENSE was performed on a 3 Tesla MRI using two single-navigator configurations (retrospective, prospective) and a combined "dual-navigator" configuration in 10 healthy adults and 20 healthy children. The adults also underwent breathhold DENSE as a reference standard for comparisons. Peak left ventricular strains, signal-to-noise ratio (SNR), and navigator efficiency were compared. Subjects also underwent dual-navigator gating with and without visual feedback to determine the effect on navigator efficiency.

RESULTS

There were no differences in circumferential, radial, and longitudinal strains between navigator-gated and breathhold DENSE (P = 0.09-0.95) (as confidence intervals, retrospective: [-1.0%-1.1%], [-7.4%-2.0%], [-1.0%-1.2%]; prospective: [-0.6%-2.7%], [-2.8%-8.3%], [-0.3%-2.9%]; dual: [-1.6%-0.5%], [-8.3%-3.2%], [-0.8%-1.9%], respectively). The dual configuration maintained SNR compared with breathhold acquisitions (16 versus 18, P = 0.06). SNR for the prospective configuration was lower than for the dual navigator in adults (P = 0.004) and children (P < 0.001). Navigator efficiency was higher (P < 0.001) for both retrospective (54%) and prospective (56%) configurations compared with the dual configuration (35%). Visual feedback improved the dual configuration navigator efficiency to 55% (P < 0.001).

CONCLUSION

When quantifying left ventricular strains using spiral cine DENSE MRI, a dual navigator configuration results in the highest SNR in adults and children. In adults, a retrospective configuration has good navigator efficiency without a substantial drop in SNR. Prospective gating should be avoided because it has the lowest SNR. Visual feedback represents an effective option to maintain navigator efficiency while using a dual navigator configuration.

LEVEL OF EVIDENCE

2 J. Magn. Reson. Imaging 2017;45:786-794.

Quantifying the effect of tissue deformation on diffusion-weighted MRI: a mathematical model and an efficient simulation framework applied to cardiac diffusion imaging

Cardiac motion presents a major challenge in diffusion weighted MRI, often leading to large signal losses that necessitate repeated measurements. The diffusion process in the myocardium is difficult to investigate because of the unqualified sensitivity of diffusion measurements to cardiac motion. A rigorous mathematical formalism is introduced to quantify the effect of tissue motion in diffusion imaging. The presented mathematical model, based on the Bloch-Torrey equations, takes into account deformations according to the laws of continuum mechanics. Approximating this mathematical model by using finite elements method, numerical simulations can predict the sensitivity of the diffusion signal to cardiac motion. Different diffusion encoding schemes are considered and the diffusion weighted MR signals, computed numerically, are compared to available results in literature. Our numerical model can identify the existence of two time points in the cardiac cycle, at which the diffusion is unaffected by myocardial strain and cardiac motion. Of course, these time points depend on the type of diffusion encoding scheme. Our numerical results also show that the motion sensitivity of the diffusion sequence can be reduced by using either spin echo technique with acceleration motion compensation diffusion gradients or stimulated echo acquisition mode with unipolar and bipolar diffusion gradients.

Accelerated two-dimensional cine DENSE cardiovascular magnetic resonance using compressed sensing and parallel imaging

BACKGROUND

Cine Displacement Encoding with Stimulated Echoes (DENSE) provides accurate quantitative imaging of cardiac mechanics with rapid displacement and strain analysis; however, image acquisition times are relatively long. Compressed sensing (CS) with parallel imaging (PI) can generally provide high-quality images recovered from data sampled below the Nyquist rate. The purposes of the present study were to develop CS-PI-accelerated acquisition and reconstruction methods for cine DENSE, to assess their accuracy for cardiac imaging using retrospective undersampling, and to demonstrate their feasibility for prospectively-accelerated 2D cine DENSE imaging in a single breathhold.

METHODS

An accelerated cine DENSE sequence with variable-density spiral k-space sampling and golden angle rotations through time was implemented. A CS method, Block LOw-rank Sparsity with Motion-guidance (BLOSM), was combined with sensitivity encoding (SENSE) for the reconstruction of under-sampled multi-coil spiral data. Seven healthy volunteers and 7 patients underwent 2D cine DENSE imaging with fully-sampled acquisitions (14-26 heartbeats in duration) and with prospectively rate-2 and rate-4 accelerated acquisitions (14 and 8 heartbeats in duration). Retrospectively- and prospectively-accelerated data were reconstructed using BLOSM-SENSE and SENSE. Image quality of retrospectively-undersampled data was quantified using the relative root mean square error (rRMSE). Myocardial displacement and circumferential strain were computed for functional assessment, and linear correlation and Bland-Altman analyses were used to compare accelerated acquisitions to fully-sampled reference datasets.

RESULTS

For retrospectively-undersampled data, BLOSM-SENSE provided similar or lower rRMSE at rate-2 and lower rRMSE at rate-4 acceleration compared to SENSE (p < 0.05, ANOVA). Similarly, for retrospective undersampling, BLOSM-SENSE provided similar or better correlation with reference displacement and strain data at rate-2 and better correlation at rate-4 acceleration compared to SENSE. Bland-Altman analyses showed similar or better agreement for displacement and strain data at rate-2 and better agreement at rate-4 using BLOSM-SENSE compared to SENSE for retrospectively-undersampled data. Rate-2 and rate-4 prospectively-accelerated cine DENSE provided good image quality and expected values of displacement and strain.

CONCLUSIONS

BLOSM-SENSE-accelerated spiral cine DENSE imaging with 2D displacement encoding can be acquired in a single breathhold of 8-14 heartbeats with high image quality and accurate assessment of myocardial displacement and circumferential strain.

Myocardial strain in healthy adults across a broad age range as revealed by cardiac magnetic resonance imaging at 1.5 and 3.0T: Associations of myocardial strain with myocardial region, age, and sex

Purpose

To assess myocardial strain using cine displacement encoding with stimulated echoes (DENSE) using 1.5T and 3.0T MRI in healthy adults.

Materials and Methods

Healthy adults without any history of cardiovascular disease underwent magnetic resonance imaging (MRI) at 1.5T and 3.0T within 2 days. The MRI protocol included balanced steady‐state free‐precession (b‐SSFP), 2D cine‐echo planar imaging (EPI)‐DENSE, and late gadolinium enhancement in subjects >45 years. Acquisitions were divided into six segments; global and segmental peak longitudinal and circumferential strain were derived and analyzed by field strength, age, and gender.

Results

In all, 89 volunteers (mean age 44.8 ± 18.0 years, range: 18–87 years) underwent MRI at 1.5T, and 88 of these subjects underwent MRI at 3.0T (1.4 ± 1.4 days between the scans). Compared with 3.0T, the magnitudes of global circumferential (–19.5 ± 2.6% vs. –18.47 ± 2.6%; P = 0.001) and longitudinal (–12.47 ± 3.2% vs. –10.53 ± 3.1%; P = 0.004) strain were greater at 1.5T. At 1.5T, longitudinal strain was greater in females than in males: –10.17 ± 3.4% vs. –13.67 ± 2.4%; P = 0.001. Similar observations occurred for circumferential strain at 1.5T (–18.72 ± 2.2% vs. –20.10 ± 2.7%; P = 0.014) and at 3.0T (–17.92 ± 1.8% vs. –19.1 ± 3.1%; P = 0.047). At 1.5T, longitudinal and circumferential strain were not associated with age after accounting for sex (longitudinal strain P = 0.178, circumferential strain P = 0.733). At 3.0T, longitudinal and circumferential strain were associated with age (P < 0.05). Longitudinal strain values were greater in the apico‐septal, basal‐lateral, and mid‐lateral segments and circumferential strain in the inferior, infero‐lateral, and antero‐lateral LV segments.

Conclusion

Myocardial strain parameters as revealed by cine‐DENSE at different MRI field strengths were associated with myocardial region, age, and sex. J. Magn. Reson. Imaging 2016;44:1197–1205.

Analyzing myocardial torsion based on tissue phase mapping cardiovascular magnetic resonance

BACKGROUND

The purpose of this work is to analyze differences in left ventricular torsion between volunteers and patients with non-ischemic cardiomyopathy based on tissue phase mapping (TPM) cardiovascular magnetic resonance (CMR).

METHODS

TPM was performed on 27 patients with non-ischemic cardiomyopathy and 14 normal volunteers. Patients underwent a standard CMR including late gadolinium enhancement (LGE) for the assessment of myocardial scar and ECG-gated cine CMR for global cardiac function. TPM was acquired in short-axis orientation at base, mid, and apex for all subjects. After evaluation by experienced observers, the patients were divided in subgroups according to the presence or absence of LGE (LGE+/LGE-), local wall motion abnormalities (WM+/WM-), and having a preserved (≥50%) or reduced (<50%) ejection fraction (EF+/EF-). TPM data was semi-automatically segmented and global LV torsion was computed for each cardiac time frame for endocardial and epicardial layers, and for the entire myocardium.

RESULTS

Maximum myocardial torsion was significantly lower for patients with reduced EF compared to controls (0.21 ± 0.15°/mm vs. 0.36 ± 0.11°/mm, p = 0.018), but also for patients with wall motion abnormalities (0.21 ± 0.13°/mm vs. 0.36 ± 0.11°/mm, p = 0.004). Global myocardial torsion showed a positive correlation (r = 0.54, p < 0.001) with EF. Moreover, endocardial torsion was significantly higher than epicardial torsion for EF+ subjects (0.56 ± 0.33°/mm vs. 0.34 ± 0.18°/mm, p = 0.039) and for volunteers (0.46 ± 0.16°/mm vs. 0.30 ± 0.09°/mm, p = 0.004). The difference in maximum torsion between endo- and epicardial layers was positively correlated with EF (r = 0.47, p = 0.002) and age (r = 0.37, p = 0.016) for all subjects.

CONCLUSIONS

TPM can be used to detect significant differences in LV torsion in patients with reduced EF and in the presence of local wall motion abnormalities. We were able to quantify torsion differences between the endocardium and epicardium, which vary between patient subgroups and are correlated to age and EF.

2D cine DENSE with low encoding frequencies accurately quantifies cardiac mechanics with improved image characteristics

BACKGROUND

Displacement Encoding with Stimulated Echoes (DENSE) encodes displacement into the phase of the magnetic resonance signal. The encoding frequency (ke) maps the measured phase to tissue displacement while the strength of the encoding gradients affects image quality. 2D cine DENSE studies have used a ke of 0.10 cycles/mm, which is high enough to remove an artifact-generating echo from k-space, provide high sensitivity to tissue displacements, and dephase the blood pool. However, through-plane dephasing can remove the unwanted echo and dephase the blood pool without relying on high ke. Additionally, the high sensitivity comes with the costs of increased phase wrapping and intra-voxel dephasing. We hypothesized that ke below 0.10 cycles/mm can be used to improve image characteristics and provide accurate measures of cardiac mechanics.

METHODS

Spiral cine DENSE images were obtained for 10 healthy subjects and 10 patients with a history of heart disease on a 3 T Siemens Trio. A mid-ventricular short-axis image was acquired with different ke: 0.02, 0.04, 0.06, 0.08, and 0.10 cycles/mm. Peak twist, circumferential strain, and radial strain were compared between acquisitions employing different ke using Bland-Altman analyses and coefficients of variation. The percentage of wrapped pixels in the phase images at end-systole was calculated for each ke. The dephasing of the blood signal and signal to noise ratio (SNR) were also calculated and compared.

RESULTS

Negligible differences were seen in strains and twist for all ke between 0.04 and 0.10 cycles/mm. These differences were of the same magnitude as inter-test differences. Specifically, the acquisitions with 0.04 cycles/mm accurately quantified cardiac mechanics and had zero phase wrapping. Compared to 0.10 cycles/mm, the acquisitions with 0.04 cycles/mm had 9 % greater SNR and negligible differences in blood pool dephasing.

CONCLUSIONS

For 2D cine DENSE with through-plane dephasing, the encoding frequency can be lowered to 0.04 cycles/mm without compromising the quantification of twist or strain. The amount of wrapping can be reduced with this lower value to greatly simplify the input to unwrapping algorithms. The strain and twist results from studies using different encoding frequencies can be directly compared.

Left ventricular mechanical dysfunction in diet-induced obese mice is exacerbated during inotropic stress: a cine DENSE cardiovascular magnetic resonance study

BACKGROUND

Obesity is a risk factor for cardiovascular disease. There is evidence of impaired left ventricular (LV) function associated with obesity, which may relate to cardiovascular mortality, but some studies have reported no dysfunction. Ventricular function data are generally acquired under resting conditions, which could mask subtle differences and potentially contribute to these contradictory findings. Furthermore, abnormal ventricular mechanics (strains, strain rates, and torsion) may manifest prior to global changes in cardiac function (i.e., ejection fraction) and may therefore represent more sensitive markers of cardiovascular disease. This study evaluated LV mechanics under both resting and stress conditions with the hypothesis that the LV mechanical dysfunction associated with obesity is exacerbated with stress and manifested at earlier stages of disease compared to baseline.

METHODS

C57BL/6J mice were randomized to a high-fat or control diet (60 %, 10 % kcal from fat, respectively) for varying time intervals (n = 7 - 10 subjects per group per time point, 100 total; 4 - 55 weeks on diet). LV mechanics were quantified under baseline (resting) and/or stress conditions (40 μg/kg/min continuous infusion of dobutamine) using cine displacement encoding with stimulated echoes (DENSE) with 7.4 ms temporal resolution on a 7 T Bruker ClinScan. Peak strain, systolic strain rates, and torsion were quantified. A linear mixed model was used with Benjamini-Hochberg adjustments for multiple comparisons.

RESULTS

Reductions in LV peak longitudinal strain at baseline were first observed in the obese group after 42 weeks, with no differences in systolic strain rates or torsion. Conversely, reductions in longitudinal strain and circumferential and radial strain rates were seen under inotropic stress conditions after only 22 weeks on diet. Furthermore, stress cardiovascular magnetic resonance (CMR) evaluation revealed supranormal values of LV radial strain and torsion in the obese group early on diet, followed by later deficits.

CONCLUSIONS

Differences in left ventricular mechanics in obese mice are exacerbated under stress conditions. Stress CMR demonstrated a broader array of mechanical dysfunction and revealed these differences at earlier time points. Thus, it may be important to evaluate cardiac function in the setting of obesity under stress conditions to fully elucidate the presence of ventricular dysfunction.

Influence of material property variability on the mechanical behaviour of carotid atherosclerotic plaques: a 3D fluid-structure interaction analysis

Mechanical analysis has been shown to be complementary to luminal stenosis in assessing atherosclerotic plaque vulnerability. However, patient-specific material properties are not available and the effect of material properties variability has not been fully quantified. Media and fibrous cap (FC) strips from carotid endarterectomy samples were classified into hard, intermediate and soft according to their incremental Young's modulus. Lipid and intraplaque haemorrhage/thrombus strips were classified as hard and soft. Idealised geometry-based 3D fluid-structure interaction analyses were performed to assess the impact of material property variability in predicting maximum principal stress (Stress-P1 ) and stretch (Stretch-P1 ). When FC was thick (1000 or 600 µm), Stress-P1 at the shoulder was insensitive to changes in material stiffness, whereas Stress-P1 at mid FC changed significantly. When FC was thin (200 or 65 µm), high stress concentrations shifted from the shoulder region to mid FC, and Stress-P1 became increasingly sensitive to changes in material properties, in particular at mid FC. Regardless of FC thickness, Stretch-P1 at these locations was sensitive to changes in material properties. Variability in tissue material properties influences both the location and overall stress/stretch value. This variability needs to be accounted for when interpreting the results of mechanical modelling.