Accelerated whole-heart 3D CSPAMM for myocardial motion quantification

Highly accelerated free-breathing real-time myocardial tagging for exercise cardiovascular magnetic resonance

BACKGROUND

Exercise cardiovascular magnetic resonance (Ex-CMR) myocardial tagging would enable quantification of myocardial deformation after exercise. However, current electrocardiogram (ECG)-segmented sequences are limited for Ex-CMR.

METHODS

We developed a highly accelerated balanced steady-state free-precession real-time tagging technique for 3 T. A 12-fold acceleration was achieved using incoherent sixfold random Cartesian sampling, twofold truncated outer phase encoding, and a deep learning resolution enhancement model. The technique was tested in two prospective studies. In a rest study of 27 patients referred for clinical CMR and 19 healthy subjects, a set of ECG-segmented for comparison and two sets of real-time tagging images for repeatability assessment were collected in 2-chamber and short-axis views with spatiotemporal resolution 2.0 × 2.0 mm2 and 29 ms. In an Ex-CMR study of 26 patients with known or suspected cardiac disease and 23 healthy subjects, real-time images were collected before and after exercise. Deformation was quantified using measures of short-axis global circumferential strain (GCS). Two experienced CMR readers evaluated the image quality of all real-time data pooled from both studies using a 4-point Likert scale for tagline quality (1-excellent; 2-good; 3-moderate; 4-poor) and artifact level (1-none; 2-minimal; 3-moderate; 4-significant). Statistical evaluation included Pearson correlation coefficient (r), intraclass correlation coefficient (ICC), and coefficient of variation (CoV).

RESULTS

In the rest study, deformation was successfully quantified in 90% of cases. There was a good correlation (r = 0.71) between ECG-segmented and real-time measures of GCS, and repeatability was good to excellent (ICC = 0.86 [0.71, 0.94]) with a CoV of 4.7%. In the Ex-CMR study, deformation was successfully quantified in 96% of subjects pre-exercise and 84% of subjects post-exercise. Short-axis and 2-chamber tagline quality were 1.6 ± 0.7 and 1.9 ± 0.8 at rest and 1.9 ± 0.7 and 2.5 ± 0.8 after exercise, respectively. Short-axis and 2-chamber artifact level was 1.2 ± 0.5 and 1.4 ± 0.7 at rest and 1.3 ± 0.6 and 1.5 ± 0.8 post-exercise, respectively.

CONCLUSION

We developed a highly accelerated real-time tagging technique and demonstrated its potential for Ex-CMR quantification of myocardial deformation. Further studies are needed to assess the clinical utility of our technique.

Cardiac Magnetic Resonance Imaging in Appraising Myocardial Strain and Biomechanics: A Current Overview

Subclinical alterations in myocardial structure and function occur early during the natural disease course. In contrast, clinically overt signs and symptoms occur during late phases, being associated with worse outcomes. Identification of such subclinical changes is critical for timely diagnosis and accurate management. Hence, implementing cost-effective imaging techniques with accuracy and reproducibility may improve long-term prognosis. A growing body of evidence supports using cardiac magnetic resonance (CMR) to quantify deformation parameters. Tissue-tagging (TT-CMR) and feature-tracking CMR (FT-CMR) can measure longitudinal, circumferential, and radial strains and recent research emphasize their diagnostic and prognostic roles in ischemic heart disease and primary myocardial illnesses. Additionally, these methods can accurately determine LV wringing and functional dynamic geometry parameters, such as LV torsion, twist/untwist, LV sphericity index, and long-axis strain, and several studies have proved their utility in prognostic prediction in various cardiovascular patients. More recently, few yet important studies have suggested the superiority of fast strain-encoded imaging CMR-derived myocardial strain in terms of accuracy and significantly reduced acquisition time, however, more studies need to be carried out to establish its clinical impact. Herein, the current review aims to provide an overview of currently available data regarding the role of CMR in evaluating myocardial strain and biomechanics.

Quantification of left ventricular strain and torsion by joint analysis of 3D tagging and cine MR images

Cardiovascular magnetic resonance (CMR) imaging is the gold standard for the non-invasive assessment of left-ventricular (LV) function. Prognostic value of deformation metrics extracted directly from regular SSFP CMR images has been shown by numerous studies in the clinical setting, but with some limitations to detect torsion of the myocardium. Tagged CMR introduces trackable features in the myocardium that allow for the assessment of local myocardial deformation, including torsion; it is, however, limited in the quantification of radial strain, which is a decisive metric for assessing the contractility of the heart. In order to improve SSFP-only and tagged-only approaches, we propose to combine the advantages of both image types by fusing global shape motion obtained from SSFP images with the local deformation obtained from tagged images. To this end, tracking is first performed on SSFP images, and subsequently, the resulting motion is utilized to mask and track tagged data. Our implementation is based on a recent finite element-based motion tracking tool with mechanical regularization. Joint SSFP and tagged images registration performance is assessed based on deformation metrics including LV strain and twist using human and in-house porcine datasets. Results show that joint analysis of SSFP and 3DTAG images provides better quantification of LV strain and twist as either data source alone.

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

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

In-silico study of accuracy and precision of left-ventricular strain quantification from 3D tagged MRI

Cardiac Magnetic Resonance Imaging (MRI) allows quantifying myocardial tissue deformation and strain based on the tagging principle. In this work, we investigate accuracy and precision of strain quantification from synthetic 3D tagged MRI using equilibrated warping. To this end, synthetic biomechanical left-ventricular tagged MRI data with varying tag distance, spatial resolution and signal-to-noise ratio (SNR) were generated and processed to quantify errors in radial, circumferential and longitudinal strains relative to ground truth. Results reveal that radial strain is more sensitive to image resolution and noise than the other strain components. The study also shows robustness of quantifying circumferential and longitudinal strain in the presence of geometrical inconsistencies of 3D tagged data. In conclusion, our study points to the need for higher-resolution 3D tagged MRI than currently available in practice in order to achieve sufficient accuracy of radial strain quantification.

Cardiovascular magnetic resonance imaging of functional and microstructural changes of the heart in a longitudinal pig model of acute to chronic myocardial infarction

BACKGROUND

We examined the dynamic response of the myocardium to infarction in a longitudinal porcine study using relaxometry, functional as well as diffusion cardiovascular magnetic resonance (CMR). We sought to compare non contrast CMR methods like relaxometry and in-vivo diffusion to contrast enhanced imaging and investigate the link of microstructural and functional changes in the acute and chronically infarcted heart.

METHODS

CMR was performed on five myocardial infarction pigs and four healthy controls. In the infarction group, measurements were obtained 2 weeks before 90 min occlusion of the left circumflex artery, 6 days after ischemia and at 5 as well as 9 weeks as chronic follow-up. The timing of measurements was replicated in the control cohort. Imaging consisted of functional cine imaging, 3D tagging, T2 mapping, native as well as gadolinium enhanced T1 mapping, cardiac diffusion tensor imaging, and late gadolinium enhancement imaging.

RESULTS

Native T1, extracellular volume (ECV) and mean diffusivity (MD) were significantly elevated in the infarcted region while fractional anisotropy (FA) was significantly reduced. During the transition from acute to chronic stages, native T1 presented minor changes (< 3%). ECV as well as MD increased from acute to the chronic stages compared to baseline: ECV: 125 ± 24% (day 6) 157 ± 24% (week 5) 146 ± 60% (week 9), MD: 17 ± 7% (day 6) 33 ± 14% (week 5) 29 ± 15% (week 9) and FA was further reduced: - 31 ± 10% (day 6) - 38 ± 8% (week 5) - 36 ± 14% (week 9). T2 as marker for myocardial edema was significantly increased in the ischemic area only during the acute stage (83 ± 3 ms infarction vs. 58 ± 2 ms control p < 0.001 and 61 ± 2 ms in the remote area p < 0.001). The analysis of functional imaging revealed reduced left ventricular ejection fraction, global longitudinal strain and torsion in the infarct group. At the same time the transmural helix angle (HA) gradient was steeper in the chronic follow-up and a correlation between longitudinal strain and transmural HA gradient was detected (r = 0.59 with p < 0.05). Comparing non-gadolinium enhanced data T2 mapping showed the largest relative change between infarct and remote during the acute stage (+ 33 ± 4% day 6, with p = 0.013 T2 vs. MD, p = 0.009 T2 vs. FA and p = 0.01 T2 vs. T1) while FA exhibited the largest relative change between infarct and remote during the chronic follow-up (+ 31 ± 2% week 5, with p = N.S. FA vs. MD, p = 0.03 FA vs. T2 and p = 0.003 FA vs. T1). Overall, diffusion parameters provided a higher contrast (> 23% for MD and > 27% for FA) during follow-up compared to relaxometry (T1 17-18%/T2 10-20%).

CONCLUSION

During chronic follow-up after myocardial infarction, cardiac diffusion tensor imaging provides a higher sensitivity for mapping microstructural alterations when compared to non-contrast enhanced relaxometry with the added benefit of providing directional tensor information to assess remodelling of myocyte aggregate orientations, which cannot be otherwise assessed.

Segmental strain analysis for the detection of chronic ischemic scars in non-contrast cardiac MRI cine images

Cardiac magnetic resonance imaging (MRI) with late gadolinium enhancement (LGE) is considered the gold standard for scar detection after myocardial infarction. In times of increasing skepticism about gadolinium depositions in brain tissue and contraindications of gadolinium administration in some patient groups, tissue strain-based techniques for detecting ischemic scars should be further developed as part of clinical protocols. Therefore, the objective of the present work was to investigate whether segmental strain is noticeably affected in chronic infarcts and thus can be potentially used for infarct detection based on routinely acquired non-contrast cine images in patients with known coronary artery disease (CAD). Forty-six patients with known CAD and chronic scars in LGE images (5 female, mean age 52 ± 19 years) and 24 gender- and age-matched controls with normal cardiac MRI (2 female, mean age 47 ± 13 years) were retrospectively enrolled. Global (global peak circumferential [GPCS], global peak longitudinal [GPLS], global peak radial strain [GPRS]) and segmental (segmental peak circumferential [SPCS], segmental peak longitudinal [SPLS], segmental peak radial strain [SPRS]) strain parameters were calculated from standard non-contrast balanced SSFP cine sequences using commercially available software (Segment CMR, Medviso, Sweden). Visual wall motion assessment of short axis cine images as well as segmental circumferential strain calculations (endo-/epicardially contoured short axis cine and resulting polar plot strain map) of every patient and control were presented in random order to two independent blinded readers, which should localize potentially infarcted segments in those datasets blinded to LGE images and patient information. Global strain values were impaired in patients compared to controls (GPCS p = 0.02; GPLS p = 0.04; GPRS p = 0.01). Patients with preserved ejection fraction showed also impeded GPCS compared to healthy individuals (p = 0.04). In patients, mean SPCS was significantly impaired in subendocardially (-  5.4% ± 2) and in transmurally infarcted segments (- 1.2% ± 3) compared to remote myocardium (- 12.9% ± 3, p = 0.02 and 0.03, respectively). ROC analysis revealed an optimal cut-off value for SPCS for discriminating infarcted from remote myocardium of - 7.2% with a sensitivity of 89.4% and specificity of 85.7%. Mean SPRS was impeded in transmurally infarcted segments (15.9% ± 6) compared to SPRS of remote myocardium (31.4% ± 5; p = 0.02). The optimal cut-off value for SPRS for discriminating scar tissue from remote myocardium was 16.6% with a sensitivity of 83.3% and specificity of 76.5%. 80.3% of all in LGE infarcted segments (118/147) were correctly localized in segmental circumferential strain calculations based on non-contrast cine images compared to 53.7% (79/147) of infarcted segments detected by visual wall motion assessment (p > 0.01). Global strain parameters are impaired in patients with chronic infarcts compared to controls. Mean SPCS and SPRS in scar tissue is impeded compared to remote myocardium in infarcts patients. Blinded to LGE images, two readers correctly localized 80% of infarcted segments in segmental circumferential strain calculations based on non-contrast cine images, in contrast to only 54% of infarcted segments detected due to wall motion abnormalities in visual wall motion assessment. Analysis of segmental circumferential strain shows a promising method for detection of chronic scars in routinely acquired, non-contrast cine images for patients who cannot receive or decline gadolinium.

HARP-I: A Harmonic Phase Interpolation Method for the Estimation of Motion From Tagged MR Images

We proposed a novel method called HARP-I, which enhances the estimation of motion from tagged Magnetic Resonance Imaging (MRI). The harmonic phase of the images is unwrapped and treated as noisy measurements of reference coordinates on a deformed domain, obtaining motion with high accuracy using Radial Basis Functions interpolations. Results were compared against Shortest Path HARP Refinement (SP-HR) and Sine-wave Modeling (SinMod), two harmonic image-based techniques for motion estimation from tagged images. HARP-I showed a favorable similarity with both methods under noise-free conditions, whereas a more robust performance was found in the presence of noise. Cardiac strain was better estimated using HARP-I at almost any motion level, giving strain maps with less artifacts. Additionally, HARP-I showed better temporal consistency as a new method was developed to fix phase jumps between frames. In conclusion, HARP-I showed to be a robust method for the estimation of motion and strain under ideal and non-ideal conditions.

A miniaturized endocardial electromagnetic energy harvester for leadless cardiac pacemakers

Life expectancy of contemporary cardiac pacemakers is limited due to the use of an internal primary battery. Repeated device replacement interventions are necessary, which leads to an elevated risk for patients and an increase of health care costs. The aim of our study is to investigate the feasibility of powering an endocardial pacemaker by converting a minimal amount of the heart's kinetic energy into electric energy. The intrinsic cardiac muscle activity makes it an ideal candidate as continuous source of energy for endocardial pacemakers. For this reason, we developed a prototype able to generate enough power to supply a pacing circuit at different heart rates. The prototype consists of a mass imbalance that drives an electromagnetic generator while oscillating. We developed a mathematical model to estimate the amount of energy harvested from the right ventricle. Finally, the implemented prototype was successfully tested during in-vitro and in-vivo experiments.

Statistical Characterization of Human Brain Deformation During Mild Angular Acceleration Measured In Vivo by Tagged Magnetic Resonance Imaging

Understanding of in vivo brain biomechanical behavior is critical in the study of traumatic brain injury (TBI) mechanisms and prevention. Using tagged magnetic resonance imaging, we measured spatiotemporal brain deformations in 34 healthy human volunteers under mild angular accelerations of the head. Two-dimensional (2D) Lagrangian strains were examined throughout the brain in each subject. Strain metrics peaked shortly after contact with a padded stop, corresponding to the inertial response of the brain after head deceleration. Maximum shear strain of at least 3% was experienced at peak deformation by an area fraction (median±standard error) of 23.5±1.8% of cortical gray matter, 15.9±1.4% of white matter, and 4.0±1.5% of deep gray matter. Cortical gray matter strains were greater in the temporal cortex on the side of the initial contact with the padded stop and also in the contralateral temporal, frontal, and parietal cortex. These tissue-level deformations from a population of healthy volunteers provide the first in vivo measurements of full-volume brain deformation in response to known kinematics. Although strains differed in different tissue type and cortical lobes, no significant differences between male and female head accelerations or strain metrics were found. These cumulative results highlight important kinematic features of the brain's mechanical response and can be used to facilitate the evaluation of computational simulations of TBI.

Hierarchical Template Matching for 3D Myocardial Tracking and Cardiac Strain Estimation

Myocardial tracking and strain estimation can non-invasively assess cardiac functioning using subject-specific MRI. As the left-ventricle does not have a uniform shape and functioning from base to apex, the development of 3D MRI has provided opportunities for simultaneous 3D tracking, and 3D strain estimation. We have extended a Local Weighted Mean (LWM) transformation function for 3D, and incorporated in a Hierarchical Template Matching model to solve 3D myocardial tracking and strain estimation problem. The LWM does not need to solve a large system of equations, provides smooth displacement of myocardial points, and adapt local geometric differences in images. Hence, 3D myocardial tracking can be performed with 1.49 mm median error, and without large error outliers. The maximum error of tracking is up to 24% reduced compared to benchmark methods. Moreover, the estimated strain can be insightful to improve 3D imaging protocols, and the computer code of LWM could also be useful for geo-spatial and manufacturing image analysis researchers.

Equilibrated warping: Finite element image registration with finite strain equilibrium gap regularization

In this paper, we propose a novel continuum finite strain formulation of the equilibrium gap regularization for image registration. The equilibrium gap regularization essentially penalizes any deviation from the solution of a hyperelastic body in equilibrium with arbitrary loads prescribed at the boundary. It thus represents a regularization with strong mechanical basis, especially suited for cardiac image analysis. We describe the consistent linearization and discretization of the regularized image registration problem, in the framework of the finite elements method. The method is implemented using FEniCS & VTK, and distributed as a freely available python library. We show that the equilibrated warping method is effective and robust: regularization strength and image noise have minimal impact on motion tracking, especially when compared to strain-based regularization methods such as hyperelastic warping. We also show that equilibrated warping is able to extract main deformation features on both tagged and untagged cardiac magnetic resonance images.

A Framework for the Generation of Realistic Synthetic Cardiac Ultrasound and Magnetic Resonance Imaging Sequences From the Same Virtual Patients

The use of synthetic sequences is one of the most promising tools for advanced in silico evaluation of the quantification of cardiac deformation and strain through 3-D ultrasound (US) and magnetic resonance (MR) imaging. In this paper, we propose the first simulation framework which allows the generation of realistic 3-D synthetic cardiac US and MR (both cine and tagging) image sequences from the same virtual patient. A state-of-the-art electromechanical (E/M) model was exploited for simulating groundtruth cardiac motion fields ranging from healthy to various pathological cases, including both ventricular dyssynchrony and myocardial ischemia. The E/M groundtruth along with template MR/US images and physical simulators were combined in a unified framework for generating synthetic data. We efficiently merged several warping strategies to keep the full control of myocardial deformations while preserving realistic image texture. In total, we generated 18 virtual patients, each with synthetic 3-D US, cine MR, and tagged MR sequences. The simulated images were evaluated both qualitatively by showing realistic textures and quantitatively by observing myocardial intensity distributions similar to real data. In particular, the US simulation showed a smoother myocardium/background interface than the state-of-the-art. We also assessed the mechanical properties. The pathological subjects were discriminated from the healthy ones by both global indexes (ejection fraction and the global circumferential strain) and regional strain curves. The synthetic database is comprehensive in terms of both pathology and modality, and has a level of realism sufficient for validation purposes. All the 90 sequences are made publicly available to the research community via an open-access database.

Studying Dynamic Myofiber Aggregate Reorientation in Dilated Cardiomyopathy Using In Vivo Magnetic Resonance Diffusion Tensor Imaging

Supplemental Digital Content is available in the text.

Background—

The objective of this study is to assess the dynamic alterations of myocardial microstructure and strain between diastole and systole in patients with dilated cardiomyopathy relative to healthy controls using the magnetic resonance diffusion tensor imaging, myocardial tagging, and biomechanical modeling.

Methods and Results—

Dual heart-phase diffusion tensor imaging was successfully performed in 9 patients and 9 controls. Tagging data were acquired for the diffusion tensor strain correction and cardiac motion analysis. Mean diffusivity, fractional anisotropy, and myocyte aggregate orientations were compared between both cohorts. Cardiac function was assessed by left ventricular ejection fraction, torsion, and strain. Computational modeling was used to study the impact of cardiac shape on fiber reorientation and how fiber orientations affect strain. In patients with dilated cardiomyopathy, a more longitudinal orientation of diastolic myofiber aggregates was measured compared with controls. Although a significant steepening of helix angles (HAs) during contraction was found in the controls, consistent change in HAs during contraction was absent in patients. Left ventricular ejection fraction, cardiac torsion, and strain were significantly lower in the patients compared with controls. Computational modeling revealed that the dilated heart results in reduced HA changes compared with a normal heart. Reduced torsion was found to be exacerbated by steeper HAs.

Conclusions—

Diffusion tensor imaging revealed reduced reorientation of myofiber aggregates during cardiac contraction in patients with dilated cardiomyopathy relative to controls. Left ventricular remodeling seems to be an important factor in the changes to myocyte orientation. Steeper HAs are coupled with a worsening in strain and torsion. Overall, the findings provide new insights into the structural alterations in patients with dilated cardiomyopathy.

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.

Estimation of myocardial deformation using correlation image velocimetry

Background

Tagged Magnetic Resonance (tMR) imaging is a powerful technique for determining cardiovascular abnormalities. One of the reasons for tMR not being used in routine clinical practice is the lack of easy-to-use tools for image analysis and strain mapping. In this paper, we introduce a novel interdisciplinary method based on correlation image velocimetry (CIV) to estimate cardiac deformation and strain maps from tMR images.

Methods

CIV, a cross-correlation based pattern matching algorithm, analyses a pair of images to obtain the displacement field at sub-pixel accuracy with any desired spatial resolution. This first time application of CIV to tMR image analysis is implemented using an existing open source Matlab-based software called UVMAT. The method, which requires two main input parameters namely correlation box size (C_B) and search box size (S_B), is first validated using a synthetic grid image with grid sizes representative of typical tMR images. Phantom and patient images obtained from a Medical Imaging grand challenge dataset (http://stacom.cardiacatlas.org/motion-tracking-challenge/) were then analysed to obtain cardiac displacement fields and strain maps. The results were then compared with estimates from Harmonic Phase analysis (HARP) technique.

Results

For a known displacement field imposed on both the synthetic grid image and the phantom image, CIV is accurate for 3-pixel and larger displacements on a 512 × 512 image with (C_B,S_B)=(25,55) pixels. Further validation of our method is achieved by showing that our estimated landmark positions on patient images fall within the inter-observer variability in the ground truth. The effectiveness of our approach to analyse patient images is then established by calculating dense displacement fields throughout a cardiac cycle, and were found to be physiologically consistent. Circumferential strains were estimated at the apical, mid and basal slices of the heart, and were shown to compare favorably with those of HARP over the entire cardiac cycle, except in a few (∼4) of the segments in the 17-segment AHA model. The radial strains, however, are underestimated by our method in most segments when compared with HARP.

Conclusions

In summary, we have demonstrated the capability of CIV to accurately and efficiently quantify cardiac deformation from tMR images. Furthermore, physiologically consistent displacement fields and circumferential strain curves in most regions of the heart indicate that our approach, upon automating some pre-processing steps and testing in clinical trials, can potentially be implemented in a clinical setting.

Electronic supplementary material

The online version of this article (doi:10.1186/s12880-017-0195-7) contains supplementary material, which is available to authorized users.

Three-dimensional Cardiac MR Imaging: Related Techniques and Clinical Applications

Three-dimensional (3D) cardiac magnetic resonance (MR) imaging has several advantages, including the easy coverage of the entire heart without misregistration, reduction of breath-holding times, and availability for postprocessing reconstruction. These advantages are associated with some techniques such as breath-hold or navigator gating and parallel imaging. However, the image quality of 3D cardiac MR images is compromised by the use of a shorter repetition time and parallel imaging. Thus, a steady-state free precession sequence, contrast agent administration, and presaturation pulses are used to maintain the image quality. In this review, we introduce the MR imaging techniques used in 3D cardiac MR imaging and demonstrate the typical 3D cardiac MR images, followed by discussion about their advantages and disadvantages.

Estimation of passive and active properties in the human heart using 3D tagged MRI

Advances in medical imaging and image processing are paving the way for personalised cardiac biomechanical modelling. Models provide the capacity to relate kinematics to dynamics and-through patient-specific modelling-derived material parameters to underlying cardiac muscle pathologies. However, for clinical utility to be achieved, model-based analyses mandate robust model selection and parameterisation. In this paper, we introduce a patient-specific biomechanical model for the left ventricle aiming to balance model fidelity with parameter identifiability. Using non-invasive data and common clinical surrogates, we illustrate unique identifiability of passive and active parameters over the full cardiac cycle. Identifiability and accuracy of the estimates in the presence of controlled noise are verified with a number of in silico datasets. Unique parametrisation is then obtained for three datasets acquired in vivo. The model predictions show good agreement with the data extracted from the images providing a pipeline for personalised biomechanical analysis.

Multiphysics and multiscale modelling, data-model fusion and integration of organ physiology in the clinic: ventricular cardiac mechanics

With heart and cardiovascular diseases continually challenging healthcare systems worldwide, translating basic research on cardiac (patho)physiology into clinical care is essential. Exacerbating this already extensive challenge is the complexity of the heart, relying on its hierarchical structure and function to maintain cardiovascular flow. Computational modelling has been proposed and actively pursued as a tool for accelerating research and translation. Allowing exploration of the relationships between physics, multiscale mechanisms and function, computational modelling provides a platform for improving our understanding of the heart. Further integration of experimental and clinical data through data assimilation and parameter estimation techniques is bringing computational models closer to use in routine clinical practice. This article reviews developments in computational cardiac modelling and how their integration with medical imaging data is providing new pathways for translational cardiac modelling.

Fast 3-Breath-Hold 3-Dimensional Tagging Cardiac Magnetic Resonance in Patients with Hypertrophic Myocardial Diseases: A Feasibility Study

Tagging CMR has been established as the standard reference for measurement of myocardial strain. The current 2D tagging technique requires multiple breath-holds to cover the whole heart and cannot show the 3D motions of the left ventricle. We performed fast 3-breath-hold 3D tagging with localized tagging preparation and complementary spatial modulation of magnetization in 10 patients with hypertrophic myocardial diseases and 6 normal volunteers. The left wall motion was observed at any view angle, which allowed for the identification of regional and global hypokinesis using the fast 3D tagging. Although a decrease in the circumferential strain and LGE were observed at the basal septum in hypertrophic cardiomyopathy, they were not located together in each patient. In hypertensive heart disease, the decrease in circumferential strain was observed more widely than LGE, and the summed strain of all segments was significantly decreased. The decrease in strain and LGE were observed diffusely in cardiac amyloidosis. In conclusion, fast 3-breath-hold 3D tagging is feasible for the regional and global strain analysis. The location of reduced circumferential strain is not necessarily the same as that of LGE and is related to the global cardiac function in patients with hypertrophic myocardial diseases.

Second-order motion compensated PRESS for cardiac spectroscopy

PURPOSE

Second-order motion compensation for point-resolved spectroscopy (PRESS) is proposed to allow for robust single-voxel cardiac spectroscopy throughout the entire cardiac cycle and at various heart rates.

METHODS

Bipolar FID spoiling gradient pairs compensating for first and second-order motion were designed and implemented into a cardiac-triggered PRESS sequence on a clinical MR system. A numerical three-dimensional model of cardiac motion was used to optimize and validate the gradient waveforms. In vivo measurements in healthy volunteers were obtained to assess the signal-to-noise ratio (SNR) and triglyceride-to-water ratio (TG/W). SNR gains and variability of TG/W of the proposed approach were evaluated against a conventional PRESS sequence with optimized gradients.

RESULTS

The proposed sequence increases the mean SNR by 32% (W) and 23% (TG) on average with significantly lower variability for different trigger delays. The variability of TG/W quantification over the cardiac cycle is significantly decreased with second-order motion compensated PRESS when compared with conventional PRESS with reduced-spoiler gradients (coefficient of variation: 0.1 ± 0.02 versus 0.37 ± 0.26).

CONCLUSION

Second-order motion compensated PRESS effectively reduces cardiac motion-induced signal degradation during FID spoiling, providing higher SNR and less variability for TG/W quantification. The sequence is considered promising to assess the TG/W modulation during various interventions including pharmacologically induced stress. Magn Reson Med 77:57-64, 2017. © 2016 Wiley Periodicals, Inc.

Multi-oriented windowed harmonic phase reconstruction for robust cardiac strain imaging

The purpose of this paper is to develop a method for direct estimation of the cardiac strain tensor by extending the harmonic phase reconstruction on tagged magnetic resonance images to obtain more precise and robust measurements. The extension relies on the reconstruction of the local phase of the image by means of the windowed Fourier transform and the acquisition of an overdetermined set of stripe orientations in order to avoid the phase interferences from structures outside the myocardium and the instabilities arising from the application of a gradient operator. Results have shown that increasing the number of acquired orientations provides a significant improvement in the reproducibility of the strain measurements and that the acquisition of an extended set of orientations also improves the reproducibility when compared with acquiring repeated samples from a smaller set of orientations. Additionally, biases in local phase estimation when using the original harmonic phase formulation are greatly diminished by the one here proposed. The ideas here presented allow the design of new methods for motion sensitive magnetic resonance imaging, which could simultaneously improve the resolution, robustness and accuracy of motion estimates.

Estimation of passive and active properties in the human heart using 3D tagged MRI

Advances in medical imaging and image processing are paving the way for personalised cardiac biomechanical modelling. Models provide the capacity to relate kinematics to dynamics and-through patient-specific modelling-derived material parameters to underlying cardiac muscle pathologies. However, for clinical utility to be achieved, model-based analyses mandate robust model selection and parameterisation. In this paper, we introduce a patient-specific biomechanical model for the left ventricle aiming to balance model fidelity with parameter identifiability. Using non-invasive data and common clinical surrogates, we illustrate unique identifiability of passive and active parameters over the full cardiac cycle. Identifiability and accuracy of the estimates in the presence of controlled noise are verified with a number of in silico datasets. Unique parametrisation is then obtained for three datasets acquired in vivo. The model predictions show good agreement with the data extracted from the images providing a pipeline for personalised biomechanical analysis.

An integrated electromechanical-growth heart model for simulating cardiac therapies

An emerging class of models has been developed in recent years to predict cardiac growth and remodeling (G&R). We recently developed a cardiac G&R constitutive model that predicts remodeling in response to elevated hemodynamics loading, and a subsequent reversal of the remodeling process when the loading is reduced. Here, we describe the integration of this G&R model to an existing strongly coupled electromechanical model of the heart. A separation of timescale between growth deformation and elastic deformation was invoked in this integrated electromechanical-growth heart model. To test our model, we applied the G&R scheme to simulate the effects of myocardial infarction in a realistic left ventricular (LV) geometry using the finite element method. We also simulate the effects of a novel therapy that is based on alteration of the infarct mechanical properties. We show that our proposed model is able to predict key features that are consistent with experiments. Specifically, we show that the presence of a non-contractile infarct leads to a dilation of the left ventricle that results in a rightward shift of the pressure volume loop. Our model also predicts that G&R is attenuated by a reduction in LV dilation when the infarct stiffness is increased.

3D harmonic phase tracking with anatomical regularization

This paper presents a novel algorithm that extends HARP to handle 3D tagged MRI images. HARP results were regularized by an original regularization framework defined in an anatomical space of coordinates. In the meantime, myocardium incompressibility was integrated in order to correct the radial strain which is reported to be more challenging to recover. Both the tracking and regularization of LV displacements were done on a volumetric mesh to be computationally efficient. Also, a window-weighted regression method was extended to cardiac motion tracking which helps maintain a low complexity even at finer scales. On healthy volunteers, the tracking accuracy was found to be as accurate as the best candidates of a recent benchmark. Strain accuracy was evaluated on synthetic data, showing low bias and strain errors under 5% (excluding outliers) for longitudinal and circumferential strains, while the second and third quartiles of the radial strain errors are in the (-5%,5%) range. In clinical data, strain dispersion was shown to correlate with the extent of transmural fibrosis. Also, reduced deformation values were found inside infarcted segments.

Effects of incremental beta-blocker dosing on myocardial mechanics of the human left ventricle: MRI 3D-tagging insight into pharmacodynamics supports theory of inner antagonism

Beta-blockers contribute to treatment of heart failure. Their mechanism of action, however, is incompletely understood. Gradients in beta-blocker sensitivity of helically aligned cardiomyocytes compared with counteracting transversely intruding cardiomyocytes seem crucial. We hypothesize that selective blockade of transversely intruding cardiomyocytes by low-dose beta-blockade unloads ventricular performance. Cardiac magnetic resonance imaging (MRI) 3D tagging delivers parameters of myocardial performance. We studied 13 healthy volunteers by MRI 3D tagging during escalated intravenous administration of esmolol. The circumferential, longitudinal, and radial myocardial shortening was determined for each dose. The curves were analyzed for peak value, time-to-peak, upslope, and area-under-the-curve. At low doses, from 5 to 25 μg·kg−1·min−1, peak contraction increased while time-to-peak decreased yielding a steeper upslope. Combining the values revealed a left shift of the curves at low doses compared with baseline without esmolol. At doses of 50 to 150 μg·kg−1·min−1, a right shift with flattening occurred. In healthy volunteers we found more pronounced myocardial shortening at low compared with clinical dosage of beta-blockers. In patients with ventricular hypertrophy and higher prevalence of transversely intruding cardiomyocytes selective low-dose beta-blockade could be even more effective. MRI 3D tagging could help to determine optimal individual beta-blocker dosing avoiding undesirable side effects.

Analysis of passive cardiac constitutive laws for parameter estimation using 3D tagged MRI

An unresolved issue in patient-specific models of cardiac mechanics is the choice of an appropriate constitutive law, able to accurately capture the passive behavior of the myocardium, while still having uniquely identifiable parameters tunable from available clinical data. In this paper, we aim to facilitate this choice by examining the practical identifiability and model fidelity of constitutive laws often used in cardiac mechanics. Our analysis focuses on the use of novel 3D tagged MRI, providing detailed displacement information in three dimensions. The practical identifiability of each law is examined by generating synthetic 3D tags from in silico simulations, allowing mapping of the objective function landscape over parameter space and comparison of minimizing parameter values with original ground truth values. Model fidelity was tested by comparing these laws with the more complex transversely isotropic Guccione law, by characterizing their passive end-diastolic pressure–volume relation behavior, as well as by considering the in vivo case of a healthy volunteer. These results show that a reduced form of the Holzapfel–Ogden law provides the best balance between identifiability and model fidelity across the tests considered.

Correlation of electrocardiogram and regional cardiac magnetic resonance imaging findings in ST-elevation myocardial infarction: a literature review

BACKGROUND

Patients with acute ST-elevation myocardial infarction (STEMI) benefit substantially from emergent coronary reperfusion. The principal mechanism is to open the occluded coronary artery to minimize myocardial injury. Thus the size of the area at risk is a critical determinant of the patient outcome, although other factors, such as reperfusion injury, have major impact on the final infarct size. Acute coronary occlusion almost immediately induces metabolic changes within the myocardium, which can be assessed with both the electrocardiogram (ECG) and cardiac magnetic resonance (CMR) imaging.

METHODS

The 12-lead ECG is the principal diagnostic method to detect and risk-stratify acute STEMI. However, to achieve a correct diagnosis, it is paramount to compare different ECG parameters with golden standards in imaging, such as CMR. In this review, we discuss aspects of ECG and CMR in the assessment of acute regional ischemic changes in the myocardium using the 17 segment model of the left ventricle presented by American Heart Association (AHA), and their relation to coronary artery anatomy.

RESULTS

Using the 17 segment model of AHA, the segments 12 and 16 remain controversial. There is an important overlap in myocardial blood supply at the antero-lateral region between LAD and LCx territories concerning these two segments.

CONCLUSION

No all-encompassing correlation can be found between ECG and CMR findings in acute ischemia with respect to coronary anatomy.

Dual-phase cardiac diffusion tensor imaging with strain correction

Purpose

In this work we present a dual-phase diffusion tensor imaging (DTI) technique that incorporates a correction scheme for the cardiac material strain, based on 3D myocardial tagging.

Methods

In vivo dual-phase cardiac DTI with a stimulated echo approach and 3D tagging was performed in 10 healthy volunteers. The time course of material strain was estimated from the tagging data and used to correct for strain effects in the diffusion weighted acquisition. Mean diffusivity, fractional anisotropy, helix, transverse and sheet angles were calculated and compared between systole and diastole, with and without strain correction. Data acquired at the systolic sweet spot, where the effects of strain are eliminated, served as a reference.

Results

The impact of strain correction on helix angle was small. However, large differences were observed in the transverse and sheet angle values, with and without strain correction. The standard deviation of systolic transverse angles was significantly reduced from 35.9±3.9° to 27.8°±3.5° (p<0.001) upon strain-correction indicating more coherent fiber tracks after correction. Myocyte aggregate structure was aligned more longitudinally in systole compared to diastole as reflected by an increased transmural range of helix angles (71.8°±3.9° systole vs. 55.6°±5.6°, p<0.001 diastole). While diastolic sheet angle histograms had dominant counts at high sheet angle values, systolic histograms showed lower sheet angle values indicating a reorientation of myocyte sheets during contraction.

Conclusion

An approach for dual-phase cardiac DTI with correction for material strain has been successfully implemented. This technique allows assessing dynamic changes in myofiber architecture between systole and diastole, and emphasizes the need for strain correction when sheet architecture in the heart is imaged with a stimulated echo approach.

A prospective evaluation of cardiovascular magnetic resonance measures of dyssynchrony in the prediction of response to cardiac resynchronization therapy

BACKGROUND

Many patients with electrical dyssynchrony who undergo cardiac resynchronization therapy (CRT) do not obtain substantial benefit. Assessing mechanical dyssynchrony may improve patient selection. Results from studies using echocardiographic imaging to measure dyssynchrony have ultimately proved disappointing. We sought to evaluate cardiac motion in patients with heart failure and electrical dyssynchrony using cardiovascular magnetic resonance (CMR). We developed a framework for comparing measures of myocardial mechanics and evaluated how well they predicted response to CRT.

METHODS

CMR was performed at 1.5 Tesla prior to CRT. Steady-state free precession (SSFP) cine images and complementary modulation of magnetization (CSPAMM) tagged cine images were acquired. Images were processed using a novel framework to extract regional ventricular volume-change, thickening and deformation fields (strain). A systolic dyssynchrony index (SDI) for all parameters within a 16-segment model of the ventricle was computed with high SDI denoting more dyssynchrony. Once identified, the optimal measure was applied to a second patient population to determine its utility as a predictor of CRT response compared to current accepted predictors (QRS duration, LBBB morphology and scar burden).

RESULTS

Forty-four patients were recruited in the first phase (91% male, 63.3 ± 14.1 years; 80% NYHA class III) with mean QRSd 154 ± 24 ms. Twenty-one out of 44 (48%) patients showed reverse remodelling (RR) with a decrease in end systolic volume (ESV) ≥ 15% at 6 months. Volume-change SDI was the strongest predictor of RR (PR 5.67; 95% CI 1.95-16.5; P = 0.003). SDI derived from myocardial strain was least predictive. Volume-change SDI was applied as a predictor of RR to a second population of 50 patients (70% male, mean age 68.6 ± 12.2 years, 76% NYHA class III) with mean QRSd 146 ± 21 ms. When compared to QRSd, LBBB morphology and scar burden, volume-change SDI was the only statistically significant predictor of RR in this group.

CONCLUSION

A systolic dyssynchrony index derived from volume-change is a highly reproducible measurement that can be derived from routinely acquired SSFP cine images and predicts RR following CRT whilst an SDI of regional strain does not.

Left ventricular pressure gating in ovine cardiac studies: a software-based method

Cardiac imaging using magnetic resonance requires a gating signal in order to compensate for motion. Human patients are routinely scanned using an electrocardiogram (ECG) as a gating signal during imaging. However, we found that in sheep the ECG is not a reliable method for gating. We developed a software based method that allowed us to use the left ventricular pressure (LVP) as a reliable gating signal. By taking the time derivative of the LVP (dP/dt), we were able to start imaging at both end-diastole for systolic phase images, and end-systole for diastolic phase images. We also used MR tissue tagging to calculate 3D strain information during diastole. Using the LVP in combination with our digital circuit provided a reliable and time efficient method for ovine cardiac imaging. Unlike the ECG signal the left ventricular pressure was a clean signal and allowed for accurate, nondelay based triggering during systole and diastole.

A computational model that predicts reverse growth in response to mechanical unloading

Ventricular growth is widely considered to be an important feature in the adverse progression of heart diseases, whereas reverse ventricular growth (or reverse remodeling) is often considered to be a favorable response to clinical intervention. In recent years, a number of theoretical models have been proposed to model the process of ventricular growth while little has been done to model its reverse. Based on the framework of volumetric strain-driven finite growth with a homeostatic equilibrium range for the elastic myofiber stretch, we propose here a reversible growth model capable of describing both ventricular growth and its reversal. We used this model to construct a semi-analytical solution based on an idealized cylindrical tube model, as well as numerical solutions based on a truncated ellipsoidal model and a human left ventricular model that was reconstructed from magnetic resonance images. We show that our model is able to predict key features in the end-diastolic pressure-volume relationship that were observed experimentally and clinically during ventricular growth and reverse growth. We also show that the residual stress fields generated as a result of differential growth in the cylindrical tube model are similar to those in other nonidentical models utilizing the same geometry.

Analysis of 3D cardiac deformations with 3D SinMod

In this paper, we propose a novel 3D sine wave modeling (3D SinMod) approach to automatic analysis of 3D cardiac deformations. An accelerated 3D complementary spatial modulation of magnetization (CSPAMM) tagging technique was used to modulate the myocardial tissue and to acquire 3D MR data sets of the whole-heart including three orthog- onal tags within three breath-holds. Each tag set is able to assess the motion along a direction perpendicular to the tag lines. With the application of CSPAMM, the effect of tag fading encountered in SPAMM tagging due to T1 relaxation is mitigated and tag deformations can be visualized for the entire cardiac cycle, including diastolic phases. In the proposed approach, the environment around each voxel in the 3D volume is modeled as a moving sine wavefront with local frequency and amplitude. The biggest advantage of the proposed technique is that the entire framework, from data acquisition to data analysis is in the 3D domain, which permits quantification of both the in-plane and through-plane motion components. The accuracy and the effectiveness of the proposed method has been validated using both simulated and in vivo tag data.

Benchmarking framework for myocardial tracking and deformation algorithms: an open access database

In this paper we present a benchmarking framework for the validation of cardiac motion analysis algorithms. The reported methods are the response to an open challenge that was issued to the medical imaging community through a MICCAI workshop. The database included magnetic resonance (MR) and 3D ultrasound (3DUS) datasets from a dynamic phantom and 15 healthy volunteers. Participants processed 3D tagged MR datasets (3DTAG), cine steady state free precession MR datasets (SSFP) and 3DUS datasets, amounting to 1158 image volumes. Ground-truth for motion tracking was based on 12 landmarks (4 walls at 3 ventricular levels). They were manually tracked by two observers in the 3DTAG data over the whole cardiac cycle, using an in-house application with 4D visualization capabilities. The median of the inter-observer variability was computed for the phantom dataset (0.77 mm) and for the volunteer datasets (0.84 mm). The ground-truth was registered to 3DUS coordinates using a point based similarity transform. Four institutions responded to the challenge by providing motion estimates for the data: Fraunhofer MEVIS (MEVIS), Bremen, Germany; Imperial College London - University College London (IUCL), UK; Universitat Pompeu Fabra (UPF), Barcelona, Spain; Inria-Asclepios project (INRIA), France. Details on the implementation and evaluation of the four methodologies are presented in this manuscript. The manually tracked landmarks were used to evaluate tracking accuracy of all methodologies. For 3DTAG, median values were computed over all time frames for the phantom dataset (MEVIS=1.20mm, IUCL=0.73 mm, UPF=1.10mm, INRIA=1.09 mm) and for the volunteer datasets (MEVIS=1.33 mm, IUCL=1.52 mm, UPF=1.09 mm, INRIA=1.32 mm). For 3DUS, median values were computed at end diastole and end systole for the phantom dataset (MEVIS=4.40 mm, UPF=3.48 mm, INRIA=4.78 mm) and for the volunteer datasets (MEVIS=3.51 mm, UPF=3.71 mm, INRIA=4.07 mm). For SSFP, median values were computed at end diastole and end systole for the phantom dataset(UPF=6.18 mm, INRIA=3.93 mm) and for the volunteer datasets (UPF=3.09 mm, INRIA=4.78 mm). Finally, strain curves were generated and qualitatively compared. Good agreement was found between the different modalities and methodologies, except for radial strain that showed a high variability in cases of lower image quality.

Understanding the mechanisms amenable to CRT response: from pre-operative multimodal image data to patient-specific computational models

This manuscript describes our recent developments towards better understanding of the mechanisms amenable to cardiac resynchronization therapy response. We report the results from a full multimodal dataset corresponding to eight patients from the euHeart project. The datasets include echocardiography, MRI and electrophysiological studies. We investigate two aspects. The first one focuses on pre-operative multimodal image data. From 2D echocardiography and 3D tagged MRI images, we compute atlas based dyssynchrony indices. We complement these indices with presence and extent of scar tissue and correlate them with CRT response. The second one focuses on computational models. We use pre-operative imaging to generate a patient-specific computational model. We show results of a fully automatic personalized electromechanical simulation. By case-per-case discussion of the results, we highlight the potential and key issues of this multimodal pipeline for the understanding of the mechanisms of CRT response and a better patient selection.

Generation of synthetic but visually realistic time series of cardiac images combining a biophysical model and clinical images

We propose a new approach for the generation of synthetic but visually realistic time series of cardiac images based on an electromechanical model of the heart and real clinical 4-D image sequences. This is achieved by combining three steps. The first step is the simulation of a cardiac motion using an electromechanical model of the heart and the segmentation of the end diastolic image of a cardiac sequence. We use biophysical parameters related to the desired condition of the simulated subject. The second step extracts the cardiac motion from the real sequence using nonrigid image registration. Finally, a synthetic time series of cardiac images corresponding to the simulated motion is generated in the third step by combining the motion estimated by image registration and the simulated one. With this approach, image processing algorithms can be evaluated as we know the ground-truth motion underlying the image sequence. Moreover, databases of visually realistic images of controls and patients can be generated for which the underlying cardiac motion and some biophysical parameters are known. Such databases can open new avenues for machine learning approaches.

Three-dimensional regional strain analysis in porcine myocardial infarction: a 3T magnetic resonance tagging study

BACKGROUND

Previous studies of mechanical strain anomalies in myocardial infarction (MI) have been largely limited to analysis of one-dimensional (1D) and two-dimensional (2D) strain parameters. Advances in cardiovascular magnetic resonance (CMR) methods now permit a complete three-dimensional (3D) interrogation of myocardial regional strain. The aim of this study was to investigate the incremental value of CMR-based 3D strain and to test the hypothesis that 3D strain is superior to 1D or 2D strain analysis in the assessment of viability using a porcine model of infarction.

METHODS

Infarction was induced surgically in 20 farm pigs. Cine, late gadolinium enhancement, and CMR tagging images were acquired at 11 days before (baseline), and 11 days (early) and 1 month (late) after induction of infarct. Harmonic phase analysis was performed to measure circumferential, longitudinal, and radial strains in myocardial segments, which were defined based on the transmurality of delayed enhancement. Univariate, bivariate, and multivariate logistic regression models of strain parameters were created and analyzed to compare the overall diagnostic accuracy of 3D strain analysis with 1D and 2D analyses in identifying the infarct and its adjacent regions from healthy myocardium.

RESULTS

3D strain differed significantly in infarct, adjacent, and remote segments (p<0.05) at early and late post-MI. In univariate, bivariate, and multivariate analyses, circumferential, longitudinal, and radial strains were significant factors (p<0.001) in differentiation of infarct and adjacent segments from baseline values. In identification of adjacent segments, receiver operating characteristic analysis using the 3D strain multivariate model demonstrated a significant improvement (p<0.01) in overall diagnostic accuracy in comparison with 2D (circumferential and radial) and 1D (circumferential) models (3D: 96%, 2D: 81%, and 1D: 71%). A similar trend was observed in identification of infarct segments.

CONCLUSIONS

Cumulative 3D strain information accurately identifies infarcts and their neighboring regions from healthy myocardium. The 3D interrogation of myocardial contractility provides incremental diagnostic accuracy in delineating the dysfunctional and nonviable myocardium in comparison with 1D or 2D quantification of strain. The infarct neighboring regions are the major beneficiaries of the 3D assessment of regional strain.

Undersampled cine 3D tagging for rapid assessment of cardiac motion

BACKGROUND

CMR allows investigating cardiac contraction, rotation and torsion non-invasively by the use of tagging sequences. Three-dimensional tagging has been proposed to cover the whole-heart but data acquisition requires three consecutive breath holds and hence demands considerable patient cooperation. In this study we have implemented and studied k-t undersampled cine 3D tagging in conjunction with k-t PCA reconstruction to potentially permit for single breath-hold acquisitions.

METHODS

The performance of undersampled cine 3D tagging was investigated using computer simulations and in-vivo measurements in 8 healthy subjects and 5 patients with myocardial infarction. Fully sampled data was obtained and compared to retrospectively and prospectively undersampled acquisitions. Fully sampled data was acquired in three consecutive breath holds. Prospectively undersampled data was obtained within a single breath hold. Based on harmonic phase (HARP) analysis, circumferential shortening, rotation and torsion were compared between fully sampled and undersampled data using Bland-Altman and linear regression analysis.

RESULTS

In computer simulations, the error for circumferential shortening was 2.8 ± 2.3% and 2.7 ± 2.1% for undersampling rates of R = 3 and 4 respectively. Errors in ventricular rotation were 2.5 ± 1.9% and 3.0 ± 2.2% for R = 3 and 4. Comparison of results from fully sampled in-vivo data acquired with prospectively undersampled acquisitions showed a mean difference in circumferential shortening of -0.14 ± 5.18% and 0.71 ± 6.16% for R = 3 and 4. The mean differences in rotation were 0.44 ± 1.8° and 0.73 ± 1.67° for R = 3 and 4, respectively. In patients peak, circumferential shortening was significantly reduced (p < 0.002 for all patients) in regions with late gadolinium enhancement.

CONCLUSION

Undersampled cine 3D tagging enables significant reduction in scan time of whole-heart tagging and facilitates quantification of shortening, rotation and torsion of the left ventricle without adding significant errors compared to previous 3D tagging approaches.

Energy harvesting from the beating heart by a mass imbalance oscillation generator

Energy-harvesting devices attract wide interest as power supplies of today's medical implants. Their long lifetime will spare patients from repeated surgical interventions. They also offer the opportunity to further miniaturize existing implants such as pacemakers, defibrillators or recorders of bio signals. A mass imbalance oscillation generator, which consists of a clockwork from a commercially available automatic wrist watch, was used as energy harvesting device to convert the kinetic energy from the cardiac wall motion to electrical energy. An MRI-based motion analysis of the left ventricle revealed basal regions to be energetically most favorable for the rotating unbalance of our harvester. A mathematical model was developed as a tool for optimizing the device's configuration. The model was validated by an in vitro experiment where an arm robot accelerated the harvesting device by reproducing the cardiac motion. Furthermore, in an in vivo experiment, the device was affixed onto a sheep heart for 1 h. The generated power in both experiments-in vitro (30 μW) and in vivo (16.7 μW)-is sufficient to power modern pacemakers.

A CMR study of the effects of tissue edema and necrosis on left ventricular dyssynchrony in acute myocardial infarction: implications for cardiac resynchronization therapy

BACKGROUND

In acute myocardial infarction (AMI), both tissue necrosis and edema are present and both might be implicated in the development of intraventricular dyssynchrony. However, their relative contribution to transient dyssynchrony is not known. Cardiovascular magnetic resonance (CMR) can detect necrosis and edema with high spatial resolution and it can quantify dyssynchrony by tagging techniques.

METHODS

Patients with a first AMI underwent percutaneous coronary interventions (PCI) of the infarct-related artery within 24 h of onset of chest pain. Within 5-7 days after the event and at 4 months, CMR was performed. The CMR protocol included the evaluation of intraventricular dyssynchrony by applying a novel 3D-tagging sequence to the left ventricle (LV) yielding the CURE index (circumferential uniformity ratio estimate; 1 = complete synchrony). On T2-weighted images, edema was measured as high-signal (> 2 SD above remote tissue) along the LV mid-myocardial circumference on 3 short-axis images (% of circumference corresponding to the area-at-risk). In analogy, on late-gadolinium enhancement (LGE) images, necrosis was quantified manually as percentage of LV mid-myocardial circumference on 3 short-axis images. Necrosis was also quantified on LGE images covering the entire LV (expressed as % LV mass). Finally, salvaged myocardium was calculated as the area-at-risk minus necrosis (expressed as % of LV circumference).

RESULTS

After successful PCI (n = 22, 2 female, mean age: 57 ± 12y), peak troponin T was 20 ± 36ug/l and the LV ejection fraction on CMR was 41 ± 8%. Necrosis mass was 30 ± 10% and CURE was 0.91 ± 0.05. Edema was measured as 58 ± 14% of the LV circumference. In the acute phase, the extent of edema correlated with dyssynchrony (r2 = -0.63, p < 0.01), while extent of necrosis showed borderline correlation (r2 = -0.19, p = 0.05). PCI resulted in salvaged myocardium of 27 ± 14%. LV dyssynchrony (=CURE) decreased at 4 months from 0.91 ± 0.05 to 0.94 ± 0.03 (p < 0.004, paired t-test). At 4 months, edema was absent and scar %LV slightly shrunk to 23.7 ± 10.0% (p < 0.002 vs baseline). Regression of LV dyssynchrony during the 4 months follow-up period was predicted by both, the extent of edema and its necrosis component in the acute phase.

CONCLUSIONS

In the acute phase of infarction, LV dyssynchrony is closely related to the extent of edema, while necrosis is a poor predictor of acute LV dyssynchrony. Conversely, regression of intraventricular LV dyssynchrony during infarct healing is predicted by the extent of necrosis in the acute phase.

A comprehensive cardiac motion estimation framework using both untagged and 3-D tagged MR images based on nonrigid registration

In this paper, we present a novel technique based on nonrigid image registration for myocardial motion estimation using both untagged and 3-D tagged MR images. The novel aspect of our technique is its simultaneous usage of complementary information from both untagged and 3-D tagged MR images. To estimate the motion within the myocardium, we register a sequence of tagged and untagged MR images during the cardiac cycle to a set of reference tagged and untagged MR images at end-diastole. The similarity measure is spatially weighted to maximize the utility of information from both images. In addition, the proposed approach integrates a valve plane tracker and adaptive incompressibility into the framework. We have evaluated the proposed approach on 12 subjects. Our results show a clear improvement in terms of accuracy compared to approaches that use either 3-D tagged or untagged MR image information alone. The relative error compared to manually tracked landmarks is less than 15% throughout the cardiac cycle. Finally, we demonstrate the automatic analysis of cardiac function from the myocardial deformation fields.

Heterogeneous distribution of left ventricular contractile injury in chronic aortic insufficiency

BACKGROUND

Global systolic strain has been described previously in patients with chronic aortic insufficiency (AI). This study explored regional differences in contractile injury.

METHODS

Tagged magnetic resonance images of the left ventricle (LV) were acquired and analyzed to calculate systolic strain in 42 patients with chronic AI. Multiparametric systolic strain analysis was applied to relate cardiac function in AI patients to a normal strain database (N = 60). AI patients were classified as having normal or poor function based on their results. A two-way repeated-measures analysis of variance was applied to analyze regional differences in injury.

RESULTS

The mean and standard deviation of raw strain values (circumferential strain, longitudinal strain, and minimum principal strain angle) are presented over the entire LV in our normal strain database. Of the 42 patients with AI, 15 could be defined as having poor function by multiparametric systolic strain analysis. In AI patients with poor function, statistical analysis showed significant differences in injury between standard LV regions (F(3.69,44.33) = 3.47, p = 0.017) and levels (F(1.49,17.88) = 4.41, p = 0.037) of the LV, whereas no significant differences were seen in the group with normal cardiac function.

CONCLUSIONS

Patients with poor function, as defined by multiparametric systolic strain z scores, exhibit a consistent, heterogeneous pattern of contractile injury in which the septum and posterior regions at the base are most injured.

First evidence of depressed contractility in the border zone of a human myocardial infarction

BACKGROUND

The temporal progression in extent and severity of regional myofiber contractile dysfunction in normally perfused border zone (BZ) myocardium adjacent to a myocardial infarction (MI) has been shown to be an important pathophysiologic feature of the adverse remodeling process in large animal models. We sought, for the first time, to document the presence of impaired contractility of the myofibers in the human BZ myocardium.

METHODS

A 62-year-old man who experienced an MI in 1985 and had recently had complete revascularization was studied. Myofiber systolic contractile stress developed in the normally perfused BZ adjacent to the MI (T(max_B)) and that developed in regions remote from the MI (T(max_R)) were quantified using cardiac catheterization, magnetic resonance imaging, and mathematical modeling.

RESULTS

The resulting finite element model of the patient's beating left ventricle was able to simulate the reduced systolic strains measured using magnetic resonance imaging at matching left ventricular pressures and volumes. The T(max_B) (73.1 kPa) was found to be greatly reduced relative to T(max_R) (109.5 kPa). These results were found to be independent of assumptions relating to BZ myofiber orientation.

CONCLUSIONS

The results of this study document the presence of impaired contractility of the myofibers in the BZ myocardium and support its role in the post-MI remodeling process in patients. To fully establish this important conclusion serial evaluations beginning at the time of the index MI will need to be performed in a cohort of patients. The current study supports the importance and demonstrates the feasibility of larger and longer-term studies.