Deformation analysis of 3D tagged cardiac images using an optical flow method

Computational Analysis of Cardiac Contractile Function

PURPOSE OF REVIEW

Heart failure results in the high incidence and mortality all over the world. Mechanical properties of myocardium are critical determinants of cardiac function, with regional variations in myocardial contractility demonstrated within infarcted ventricles. Quantitative assessment of cardiac contractile function is therefore critical to identify myocardial infarction for the early diagnosis and therapeutic intervention.

RECENT FINDINGS

Current advancement of cardiac functional assessments is in pace with the development of imaging techniques. The methods tailored to advanced imaging have been widely used in cardiac magnetic resonance, echocardiography, and optical microscopy. In this review, we introduce fundamental concepts and applications of representative methods for each imaging modality used in both fundamental research and clinical investigations. All these methods have been designed or developed to quantify time-dependent 2-dimensional (2D) or 3D cardiac mechanics, holding great potential to unravel global or regional myocardial deformation and contractile function from end-systole to end-diastole. Computational methods to assess cardiac contractile function provide a quantitative insight into the analysis of myocardial mechanics during cardiac development, injury, and remodeling.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Neural-Network-Based Diagnosis Using 3-Dimensional Myocardial Architecture and Deformation: Demonstration for the Differentiation of Hypertrophic Cardiomyopathy

The diagnosis of cardiomyopathy states may benefit from machine-learning (ML) based approaches, particularly to distinguish those states with similar phenotypic characteristics. Three-dimensional myocardial deformation analysis (3D-MDA) has been validated to provide standardized descriptors of myocardial architecture and deformation, and may therefore offer appropriate features for the training of ML-based diagnostic tools. We aimed to assess the feasibility of automated disease diagnosis using a neural network trained using 3D-MDA to discriminate hypertrophic cardiomyopathy (HCM) from its mimic states: cardiac amyloidosis (CA), Anderson–Fabry disease (AFD), and hypertensive cardiomyopathy (HTNcm). 3D-MDA data from 163 patients (mean age 53.1 ± 14.8 years; 68 females) with left ventricular hypertrophy (LVH) of known etiology was provided. Source imaging data was from cardiac magnetic resonance (CMR). Clinical diagnoses were as follows: 85 HCM, 30 HTNcm, 30 AFD, and 18 CA. A fully-connected-layer feed-forward neural was trained to distinguish HCM vs. other mimic states. Diagnostic performance was compared to threshold-based assessments of volumetric and strain-based CMR markers, in addition to baseline clinical patient characteristics. Threshold-based measures provided modest performance, the greatest area under the curve (AUC) being 0.70. Global strain parameters exhibited reduced performance, with AUC under 0.64. A neural network trained exclusively from 3D-MDA data achieved an AUC of 0.94 (sensitivity 0.92, specificity 0.90) when performing the same task. This study demonstrates that ML-based diagnosis of cardiomyopathy states performed exclusively from 3D-MDA is feasible and can distinguish HCM from mimic disease states. These findings suggest strong potential for computer-assisted diagnosis in clinical practice.

Temporally coherent cardiac motion tracking from cine MRI: Traditional registration method and modern CNN method

Purpose

Cardiac motion tracking enables quantitative evaluation of myocardial strain, which is clinically interesting in cardiovascular disease research. However, motion tracking is difficult to perform manually. In this paper, we aim to develop and compare two fully automated motion tracking methods for the steady state free precession (SSFP) cine magnetic resonance imaging (MRI), and explore their use in real clinical scenario with different patient groups.

Methods

We proposed two automated cardiac motion tracking method: (a) a traditional registration‐based method, named full cardiac cycle registration, which simultaneously tracks all cine frames within a full cardiac cycle by joint registration of all frames; and (b) a modern convolutional neural network (CNN)‐based method, named Groupwise MotionNet, which enhances the temporal coherence by fusing motion along a continuous time scale. Both methods were evaluated on the healthy volunteer data from the MICCAI 2011 STACOM Challenge, as well as on patient data including hypertrophic cardiomyopathy (HCM) and myocardial infarction (MI).

Results

The full cardiac cycle registration method achieved an average end‐point error (EPE) 2.89 ± 1.57 mm for cardiac motion tracking, with computation time of around 9 min per short‐axis cine MRI (size 128 × 128, 30 cardiac phases). In comparison, the Groupwise MotionNet achieved an average EPE of 0.94 ± 1.59 mm, taking < 1 s for a full cardiac phases. Further experiments showed that registration method had stable performance, independent of patient cohort and MRI machine, while the CNN‐based method relied on the training data to deliver consistently accurate results.

Conclusion

Both registration‐based and CNN‐based method can track the cardiac motion from SSFP cine MRI in a fully automated manner, while taking temporal coherence into account. The registration method is generic, robust, but relatively slow; the CNN‐based method trained with heterogeneous data was able to achieve high tracking accuracy with real‐time performance.

A gradient-based optical-flow cardiac motion estimation method for cine and tagged MR images

A new method is proposed to quantify the myocardial motion from both 2D C(ine)-MRI and T(agged)-MRI sequences. The tag pattern offers natural landmarks within the image that makes it possible to accurately quantify the motion within the myocardial wall. Therefore, several methods have been proposed for T-MRI. However, the lack of salient features within the cardiac wall in C-MRI hampers local motion estimation. Our method aims to ensure the local intensity and shape features invariance during motion through the iterative minimization of a cost function via a random walk scheme. The proposed approach is evaluated on realistic simulated C-MRI and T-MRI sequences. The results show more than 53% improvements on displacement estimation, and more than 24% on strain estimation for both C-MRI and T-MRI sequences, as compared to state-of-the-art cardiac motion estimators. Preliminary experiments on clinical data have shown a good ability of the proposed method to detect abnormal motion patterns related to pathology. If those results are confirmed on large databases, this would open up the possibility for more accurate diagnosis of cardiac function from standard C-MRI examinations and also the retrospective study of prior studies.

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.

Effects of hydrogel injection on borderzone contractility post-myocardial infarction

Injectable hydrogels are a potential therapy for mitigating adverse left ventricular (LV) remodeling after myocardial infarction (MI). Previous studies using magnetic resonance imaging (MRI) have shown that hydrogel treatment improves systolic strain in the borderzone (BZ) region surrounding the infarct. However, the corresponding contractile properties of the BZ myocardium are still unknown. The goal of the current study was to quantify the in vivo contractile properties of the BZ myocardium post-MI in an ovine model treated with an injectable hydrogel. Contractile properties were determined 8 weeks following posterolateral MI by minimizing the difference between in vivo strains and volume calculated from MRI and finite element model predicted strains and volume. This was accomplished by using a combination of MRI, catheterization, finite element modeling, and numerical optimization. Results show contractility in the BZ of animals treated with hydrogel injection was significantly higher than untreated controls. End-systolic (ES) fiber stress was also greatly reduced in the BZ of treated animals. The passive stiffness of the treated infarct region was found to be greater than the untreated control. Additionally, the wall thickness in the infarct and BZ regions was found to be significantly higher in the treated animals. Treatment with hydrogel injection significantly improved BZ function and reduced LV remodeling, via altered MI properties. These changes are linked to a reduction in the ES fiber stress in the BZ myocardium surrounding the infarct. The current results imply that injectable hydrogels could be a viable therapy for maintaining LV function post-MI.

Myocardial motion analysis based on an optical flow method using tagged MR images

We developed a method of velocimetry based on an optical flow method using quantitative analyses of tagged magnetic resonance (MR) images (tagged MR-optical flow velocimetry, tMR-O velocimetry). The purpose of our study was to examine the accuracy of measurement of the proposed tMR-O velocimetry. We performed retrospective pseudo-electrocardiogram (ECG) gating tagged cine MR imaging on a rotating phantom. We optimized imaging parameters for tagged MR imaging, and validated the accuracy of tMR-O velocimetry. Our results indicated that the difference between the reference velocities and the computed velocities measured using optimal imaging parameters was less than 1%. In addition, we performed tMR-O velocimetry and echocardiography on 10 healthy volunteers, for four sections of the heart (apical, midventricular, and basal sections aligned with the short-axis, and a four-chamber section aligned with the long-axis), and obtained radial and longitudinal myocardial velocities in these sections. We compared the myocardial velocities obtained using tMR-O velocimetry with those obtained using echocardiography. Our results showed good agreement between tMR-O velocimetry and echocardiography in the radial myocardial velocities in three short-axial sections and longitudinal myocardial velocities on the midventricular portion of the four-chamber section in the long-axis. In the study conducted on the rotating phantom, tMR-O velocimetry showed high accuracy; moreover, in the healthy volunteers, the myocardial velocities obtained using tMR-O velocimetry were relatively similar to those obtained using echocardiography. In conclusion, tMR-O velocimetry is a potentially feasible method for analyzing myocardial motion in the human heart.

Computational Investigation of Transmural Differences in Left Ventricular Contractility

Myocardial contractility of the left ventricle (LV) plays an essential role in maintaining normal pump function. A recent ex vivo experimental study showed that cardiomyocyte force generation varies across the three myocardial layers of the LV wall. However, the in vivo distribution of myocardial contractile force is still unclear. The current study was designed to investigate the in vivo transmural distribution of myocardial contractility using a noninvasive computational approach. For this purpose, four cases with different transmural distributions of maximum isometric tension (Tmax) and/or reference sarcomere length (lR) were tested with animal-specific finite element (FE) models, in combination with magnetic resonance imaging (MRI), pressure catheterization, and numerical optimization. Results of the current study showed that the best fit with in vivo MRI-derived deformation was obtained when Tmax assumed different values in the subendocardium, midmyocardium, and subepicardium with transmurally varying lR. These results are consistent with recent ex vivo experimental studies, which showed that the midmyocardium produces more contractile force than the other transmural layers. The systolic strain calculated from the best-fit FE model was in good agreement with MRI data. Therefore, the proposed noninvasive approach has the capability to predict the transmural distribution of myocardial contractility. Moreover, FE models with a nonuniform distribution of myocardial contractility could provide a better representation of LV function and be used to investigate the effects of transmural changes due to heart disease.

Fully automatic left ventricular myocardial strain estimation in 2D short-axis tagged magnetic resonance imaging

Cardiovascular diseases are among the leading causes of death and frequently result in local myocardial dysfunction. Among the numerous imaging modalities available to detect these dysfunctional regions, cardiac deformation imaging through tagged magnetic resonance imaging (t-MRI) has been an attractive approach. Nevertheless, fully automatic analysis of these data sets is still challenging. In this work, we present a fully automatic framework to estimate left ventricular myocardial deformation from t-MRI. This strategy performs automatic myocardial segmentation based on B-spline explicit active surfaces, which are initialized using an annular model. A non-rigid image-registration technique is then used to assess myocardial deformation. Three experiments were set up to validate the proposed framework using a clinical database of 75 patients. First, automatic segmentation accuracy was evaluated by comparing against manual delineations at one specific cardiac phase. The proposed solution showed an average perpendicular distance error of 2.35  ±  1.21 mm and 2.27  ±  1.02 mm for the endo- and epicardium, respectively. Second, starting from either manual or automatic segmentation, myocardial tracking was performed and the resulting strain curves were compared. It is shown that the automatic segmentation adds negligible differences during the strain-estimation stage, corroborating its accuracy. Finally, segmental strain was compared with scar tissue extent determined by delay-enhanced MRI. The results proved that both strain components were able to distinguish between normal and infarct regions. Overall, the proposed framework was shown to be accurate, robust, and attractive for clinical practice, as it overcomes several limitations of a manual analysis.

Clinical feasibility and validation of 3D principal strain analysis from cine MRI: comparison to 2D strain by MRI and 3D speckle tracking echocardiography

Two-dimensional (2D) strain analysis is constrained by geometry-dependent reference directions of deformation (i.e. radial, circumferential, and longitudinal) following the assumption of cylindrical chamber architecture. Three-dimensional (3D) principal strain analysis may overcome such limitations by referencing intrinsic (i.e. principal) directions of deformation. This study aimed to demonstrate clinical feasibility of 3D principal strain analysis from routine 2D cine MRI with validation to strain from 2D tagged cine analysis and 3D speckle tracking echocardiography. Thirty-one patients undergoing cardiac MRI were studied. 3D strain was measured from routine, multi-planar 2D cine SSFP images using custom software designed to apply 4D deformation fields to 3D cardiac models to derive principal strain. Comparisons of strain estimates versus those by 2D tagged cine, 2D non-tagged cine (feature tracking), and 3D speckle tracking echocardiography (STE) were performed. Mean age was 51 ± 14 (36% female). Mean LV ejection fraction was 66 ± 10% (range 37–80%). 3D principal strain analysis was feasible in all subjects and showed high inter- and intra-observer reproducibility (ICC range 0.83–0.97 and 0.83–0.98, respectively—p < 0.001 for all directions). Strong correlations of minimum and maximum principal strain were respectively observed versus the following: 3D STE estimates of longitudinal (r = 0.81 and r = −0.64), circumferential (r = 0.76 and r = −0.58) and radial (r = −0.80 and r = 0.63) strain (p < 0.001 for all); 2D tagged cine estimates of longitudinal (r = 0.81 and r = −0.81), circumferential (r = 0.87 and r = −0.85), and radial (r = −0.76 and r = 0.81) strain (p < 0.0001 for all); and 2D cine (feature tracking) estimates of longitudinal (r = 0.85 and −0.83), circumferential (r = 0.88 and r = −0.87), and radial strain (r = −0.79 and r = 0.84, p < 0.0001 for all). 3D principal strain analysis is feasible using routine, multi-planar 2D cine MRI and shows high reproducibility with strong correlations to 2D conventional strain analysis and 3D STE-based analysis. Given its independence from geometry-related directions of deformation this technique may offer unique benefit for the detection and prognostication of myocardial disease, and warrants expanded investigation.

Passive diastolic modelling of human ventricles: Effects of base movement and geometrical heterogeneity

Left-ventricular (LV) remodelling, associated with diastolic heart failure, is driven by an increase in myocardial stress. Therefore, normalisation of LV wall stress is the cornerstone of many therapeutic treatments. However, information regarding such regional stress-strain for human LV is still limited. Thus, the objectives of our study were to determine local diastolic stress-strain field in healthy LVs, and consequently, to identify the regional variations amongst them due to geometric heterogeneity. Effects of LV base movement on diastolic model predictions, which were ignored in the literature, were further explored. Personalised finite-element modelling of five normal human bi-ventricles was carried out using subject-specific myocardium properties. Model prediction was validated individually through comparison with end-diastolic volume and a new shape-volume based measurement of LV cavity, extracted from magnetic resonance imaging. Results indicated that incorporation of LV base movement improved the model predictions (shape-volume relevancy of LV cavity), and therefore, it should be considered in future studies. The LV endocardium always experienced higher fibre stress compared to the epicardium for all five subjects. The LV wall near base experienced higher stress compared to equatorial and apical locations. The lateral LV wall underwent greater stress distribution (fibre and sheet stress) compared to other three regions. In addition, normal ranges of different stress-strain components in different regions of LV wall were reported for five healthy ventricles. This information could be used as targets for future computational studies to optimise diastolic heart failure treatments or design new therapeutic interventions/devices.

A graph theoretic approach for computing 3D+time biventricular cardiac strain from tagged MRI data

Tagged magnetic resonance imaging (tMRI) is a well-established method for evaluating regional mechanical function of the heart. Many techniques have been developed to compute 2D or 3D cardiac deformation and strain from tMRI images. In this paper, we present a new method for measuring 3D plus time biventricular myocardial strain from tMRI data. The method is composed of two parts. First, we use a Gabor filter bank to extract tag points along tag lines. Second, each tag point is classified to one of a set of indexed reference tag lines using a point classification with graph cuts (PCGC) algorithm and a motion compensation technique. 3D biventricular deformation and strain is computed at each image time frame from the classified tag points using a previously published finite difference method. The strain computation is fully automatic after myocardial contours are defined near end-diastole and end-systole. An in-vivo dataset composed of 30 human imaging studies with a range of pathologies was used for validation. Strains computed with the PCGC method with no manual corrections were compared to strains computed from both manually placed tag points and a manually-corrected unwrapped phase method. A typical cardiac imaging study with 10 short-axis slices and 6 long-axis slices required 30 min for contouring followed by 44 min of automated processing. The results demonstrate that the proposed method can reconstruct accurate 3D plus time cardiac strain maps with minimal user intervention.

Computational Modeling of Healthy Myocardium in Diastole

In order to better understand the mechanics of the heart and its disorders, engineers increasingly make use of the finite element method (FEM) to investigate healthy and diseased cardiac tissue. However, FEM is only as good as the underlying constitutive model, which remains a major challenge to the biomechanics community. In this study, a recently developed structurally based constitutive model was implemented to model healthy left ventricular myocardium during passive diastolic filling. This model takes into account the orthotropic response of the heart under loading. In-vivo strains were measured from magnetic resonance images (MRI) of porcine hearts, along with synchronous catheterization pressure data, and used for parameter identification of the passive constitutive model. Optimization was performed by minimizing the difference between MRI measured and FE predicted strains and cavity volumes. A similar approach was followed for the parameter identification of a widely used phenomenological constitutive law, which is based on a transversely isotropic material response. Results indicate that the parameter identification with the structurally based constitutive law is more sensitive to the assigned fiber architecture and the fit between the measured and predicted strains is improved with more realistic sheet angles. In addition, the structurally based model is capable of generating a more physiological end-diastolic pressure-volume relationship in the ventricle.

Myocardial motion estimation of tagged cardiac magnetic resonance images using tag motion constraints and multi-level b-splines interpolation

Myocardial motion estimation of tagged cardiac magnetic resonance (TCMR) images is of great significance in clinical diagnosis and the treatment of heart disease. Currently, the harmonic phase analysis method (HARP) and the local sine-wave modeling method (SinMod) have been proven as two state-of-the-art motion estimation methods for TCMR images, since they can directly obtain the inter-frame motion displacement vector field (MDVF) with high accuracy and fast speed. By comparison, SinMod has better performance over HARP in terms of displacement detection, noise and artifacts reduction. However, the SinMod method has some drawbacks: 1) it is unable to estimate local displacements larger than half of the tag spacing; 2) it has observable errors in tracking of tag motion; and 3) the estimated MDVF usually has large local errors. To overcome these problems, we present a novel motion estimation method in this study. The proposed method tracks the motion of tags and then estimates the dense MDVF by using the interpolation. In this new method, a parameter estimation procedure for global motion is applied to match tag intersections between different frames, ensuring specific kinds of large displacements being correctly estimated. In addition, a strategy of tag motion constraints is applied to eliminate most of errors produced by inter-frame tracking of tags and the multi-level b-splines approximation algorithm is utilized, so as to enhance the local continuity and accuracy of the final MDVF. In the estimation of the motion displacement, our proposed method can obtain a more accurate MDVF compared with the SinMod method and our method can overcome the drawbacks of the SinMod method. However, the motion estimation accuracy of our method depends on the accuracy of tag lines detection and our method has a higher time complexity.

MRI evaluation of injectable hyaluronic acid-based hydrogel therapy to limit ventricular remodeling after myocardial infarction

Injectable biomaterials are an attractive therapy to attenuate left ventricular (LV) remodeling after myocardial infarction (MI). Although studies have shown that injectable hydrogels improve cardiac structure and function in vivo, temporal changes in infarct material properties after treatment have not been assessed. Emerging imaging and modeling techniques now allow for serial, non-invasive estimation of infarct material properties. Specifically, cine magnetic resonance imaging (MRI) assesses global LV structure and function, late-gadolinium enhancement (LGE) MRI enables visualization of infarcted tissue to quantify infarct expansion, and spatial modulation of magnetization (SPAMM) tagging provides passive wall motion assessment as a measure of tissue strain, which can all be used to evaluate infarct properties when combined with finite element (FE) models. In this work, we investigated the temporal effects of degradable hyaluronic acid (HA) hydrogels on global LV remodeling, infarct thinning and expansion, and infarct stiffness in a porcine infarct model for 12 weeks post-MI using MRI and FE modeling. Hydrogel treatment led to decreased LV volumes, improved ejection fraction, and increased wall thickness when compared to controls. FE model simulations demonstrated that hydrogel therapy increased infarct stiffness for 12 weeks post-MI. Thus, evaluation of myocardial tissue properties through MRI and FE modeling provides insight into the influence of injectable hydrogel therapies on myocardial structure and function post-MI.

Temporal Changes in Infarct Material Properties: An In Vivo Assessment Using Magnetic Resonance Imaging and Finite Element Simulations

BACKGROUND

Infarct expansion initiates and sustains adverse left ventricular (LV) remodeling after myocardial infarction (MI) and is influenced by temporal changes in infarct material properties. Data from ex vivo biaxial extension testing support this hypothesis; however, infarct material properties have never been measured in vivo. The goal of the current study was to serially quantify the in vivo material properties and fiber orientation of infarcted myocardium over a 12-week period in a porcine model of MI.

METHODS

A combination of magnetic resonance imaging (MRI), catheterization, finite element modeling, and numeric optimization was used to analyze posterolateral MI. Specifically, properties were determined by minimizing the difference between in vivo strains and volume calculated from MRI and strains and volume predicted by finite element modeling.

RESULTS

In 1 week after MI, the infarct region was found to be approximately 20 times stiffer than normal diastolic myocardium. Over the course of 12 weeks, the infarct region became progressively less stiff as the LV dilated and ejection fraction decreased. The infarct thinned by nearly half during the remodeling period, and infarct fiber angles became more circumferentially oriented.

CONCLUSIONS

The results reported here are consistent with previously described ex vivo biaxial extension studies of infarct material properties and the circumferential change of collagen orientation in posterolateral infarcts. The current study represents a significant advance in that the method used allows for the serial assessment of an individual infarct in vivo over time and avoids the inherent limitations related to the testing of excised tissues.

Injectable microsphere gel progressively improves global ventricular function, regional contractile strain, and mitral regurgitation after myocardial infarction

BACKGROUND

There is continued need for therapies which reverse or abate the remodeling process after myocardial infarction (MI). In this study, we evaluate the longitudinal effects of calcium hydroxyapatite microsphere gel on regional strain, global ventricular function, and mitral regurgitation (MR) in a porcine MI model.

METHODS

Twenty-five Yorkshire swine were enrolled. Five were dedicated weight-matched controls. Twenty underwent posterolateral infarction by direct ligation of the circumflex artery and its branches. Infarcted animals were randomly divided into the following 4 groups: 1-week treatment; 1-week control; 4-week treatment; and 4-week control. After infarction, animals received either twenty 150 μL calcium hydroxyapatite gel or saline injections within the infarct. At their respective time points, echocardiograms, cardiac magnetic resonance imaging, and tissue were collected for evaluation of MR, regional and global left ventricular function, wall thickness, and collagen content.

RESULTS

Global and regional left ventricular functions were depressed in all infarcted subjects at 1 week compared with healthy controls. By 4-weeks post-infarction, global function had significantly improved in the calcium hydroxyapatite group compared with infarcted controls (ejection fraction 0.485 ± 0.019 vs 0.38 ± 0.017, p < 0.01). Similarly, regional borderzone radial contractile strain (16.3% ± 1.5% vs 11.2% ± 1.5%, p = 0.04), MR grade (0.4 ± 0.2 vs 1.2 ± 0.2, p = 0.04), and infarct thickness (7.8 ± 0.5 mm vs 4.5 ± 0.2 mm, p < 0.01) were improved at this time point in the treatment group compared with infarct controls.

CONCLUSIONS

Calcium hydroxyapatite injection after MI progressively improves global left ventricular function, borderzone function, and mitral regurgitation. Using novel biomaterials to augment infarct material properties is a viable alternative in the current management of heart failure.

Regional myocardial three-dimensional principal strains during postinfarction remodeling

BACKGROUND

The purpose of this study was to quantify myocardial three-dimensional (3D) principal strains as the left ventricle (LV) remodels after myocardial infarction (MI). Serial quantification of myocardial strains is important for understanding the mechanical response of the LV to MI. Principal strains convert the 3D LV wall-based strain matrix with three normal and three shear elements, to a matrix with three nonzero normal elements, thereby eliminating the shear elements, which are difficult to physically interpret.

METHODS

The study was designed to measure principal strains of the remote, border zone, and infarct regions in a porcine model of post-MI LV remodeling. Magnetic resonance imaging was used to measure function and strain at baseline, 1 week, and 4 weeks after infarct. Principal strain was measured using 3D acquisition and the optical flow method for displacement tracking.

RESULTS

Principal strains were altered as the LV remodeled. Maximum principal strain magnitude decreased in all regions, including the noninfarcted remote, while maximum principal strain angles rotated away from the radial direction in the border zone and infarct. Minimum principal strain magnitude followed a similar pattern; however, strain angles were altered in all regions. Evolution of principal strains correlated with adverse LV remodeling.

CONCLUSIONS

Using a state-of-the-art imaging and optical flow method technique, 3D principal strains can be measured serially after MI in pigs. Results are consistent with progressive infarct stretching as well as with decreased contractile function in the border zone and remote myocardial regions.

Estimating passive mechanical properties in a myocardial infarction using MRI and finite element simulations

Myocardial infarction (MI) triggers a series of maladaptive events that lead to structural and functional changes in the left ventricle. It is crucial to better understand the progression of adverse remodeling, in order to develop effective treatment. In addition, being able to assess changes in vivo would be a powerful tool in the clinic. The goal of the current study is to quantify the in vivo material properties of infarcted and remote myocardium 1 week after MI, as well as the orientation of collagen fibers in the infarct. This will be accomplished by using a combination of magnetic resonance imaging (MRI), catheterization, finite element modeling, and numerical optimization to analyze a porcine model ([Formula: see text]) of posterolateral myocardial infarction. Specifically, properties will be determined by minimizing the difference between in vivo strains and volume calculated from MRI and finite element model predicted strains and volume. The results indicate that the infarct region is stiffer than the remote region and that the infarct collagen fibers become more circumferentially oriented 1 week post-MI. These findings are consistent with previous studies, which employed ex vivo techniques. The proposed methodology will ultimately provide a means of predicting remote and infarct mechanical properties in vivo at any time point post-MI.

Estimating passive mechanical properties in a myocardial infarction using MRI and finite element simulations

Myocardial infarction (MI) triggers a series of maladaptive events that lead to structural and functional changes in the left ventricle. It is crucial to better understand the progression of adverse remodeling, in order to develop effective treatment. In addition, being able to assess changes in vivo would be a powerful tool in the clinic. The goal of the current study is to quantify the in vivo material properties of infarcted and remote myocardium 1 week after MI, as well as the orientation of collagen fibers in the infarct. This will be accomplished by using a combination of magnetic resonance imaging (MRI), catheterization, finite element modeling, and numerical optimization to analyze a porcine model ([Formula: see text]) of posterolateral myocardial infarction. Specifically, properties will be determined by minimizing the difference between in vivo strains and volume calculated from MRI and finite element model predicted strains and volume. The results indicate that the infarct region is stiffer than the remote region and that the infarct collagen fibers become more circumferentially oriented 1 week post-MI. These findings are consistent with previous studies, which employed ex vivo techniques. The proposed methodology will ultimately provide a means of predicting remote and infarct mechanical properties in vivo at any time point post-MI.

Left ventricular strain and its pattern estimated from cine CMR and validation with DENSE

Measurement of local strain provides insight into the biomechanical significance of viable myocardium. We attempted to estimate myocardial strain from cine cardiovascular magnetic resonance (CMR) images by using a b-spline deformable image registration method. Three healthy volunteers and 41 patients with either recent or chronic myocardial infarction (MI) were studied at 1.5 Tesla with both cine and DENSE CMR. Regional circumferential and radial left ventricular strains were estimated from cine and DENSE acquisitions. In all healthy volunteers, there was no difference for peak circumferential strain (- 0.18 ± 0.04 versus - 0.18 ± 0.03, p = 0.76) between cine and DENSE CMR, however peak radial strain was overestimated from cine (0.84 ± 0.37 versus 0.49 ± 0.2, p < 0.01). In the patient study, the peak strain patterns predicted by cine were similar to the patterns from DENSE, including the strain evolution related to recovery time and strain patterns related to MI scar extent. Furthermore, cine-derived strain disclosed different strain patterns in MI and non-MI regions, and regions with transmural and non-transmural MI as DENSE. Although there were large variations with radial strain measurements from cine CMR images, useful circumferential strain information can be obtained from routine clinical CMR imaging. Cine strain analysis has potential to improve the diagnostic yield from routine CMR imaging in clinical practice.

Preclinical evaluation of the engineered stem cell chemokine stromal cell-derived factor 1α analog in a translational ovine myocardial infarction model

RATIONALE

After myocardial infarction, there is an inadequate blood supply to the myocardium, and the surrounding borderzone becomes hypocontractile.

OBJECTIVE

To develop a clinically translatable therapy, we hypothesized that in a preclinical ovine model of myocardial infarction, the modified endothelial progenitor stem cell chemokine, engineered stromal cell-derived factor 1α analog (ESA), would induce endothelial progenitor stem cell chemotaxis, limit adverse ventricular remodeling, and preserve borderzone contractility.

METHODS AND RESULTS

Thirty-six adult male Dorset sheep underwent permanent ligation of the left anterior descending coronary artery, inducing an anteroapical infarction, and were randomized to borderzone injection of saline (n=18) or ESA (n=18). Ventricular function, geometry, and regional strain were assessed using cardiac MRI and pressure-volume catheter transduction. Bone marrow was harvested for in vitro analysis, and myocardial biopsies were taken for mRNA, protein, and immunohistochemical analysis. ESA induced greater chemotaxis of endothelial progenitor stem cells compared with saline (P<0.01) and was equivalent to recombinant stromal cell-derived factor 1α (P=0.27). Analysis of mRNA expression and protein levels in ESA-treated animals revealed reduced matrix metalloproteinase 2 in the borderzone (P<0.05), with elevated levels of tissue inhibitor of matrix metalloproteinase 1 and elastin in the infarct (P<0.05), whereas immunohistochemical analysis of borderzone myocardium showed increased capillary and arteriolar density in the ESA group (P<0.01). Animals in the ESA treatment group also had significant reductions in infarct size (P<0.01), increased maximal principle strain in the borderzone (P<0.01), and a steeper slope of the end-systolic pressure-volume relationship (P=0.01).

CONCLUSIONS

The novel, biomolecularly designed peptide ESA induces chemotaxis of endothelial progenitor stem cells, stimulates neovasculogenesis, limits infarct expansion, and preserves contractility in an ovine model of myocardial infarction.

Method to create regional mechanical dyssynchrony maps from short-axis cine steady-state free-precession images

PURPOSE

To develop a robust method to assess regional mechanical dyssynchrony from cine short-axis MR images. Cardiac resynchronization therapy (CRT) is an effective treatment for patients with heart failure and evidence of left-ventricular (LV) dyssynchrony. Patient response to CRT is greatest when the LV pacing lead is placed in the most dyssynchronous segment. Existing techniques for assessing regional dyssynchrony require difficult acquisition and/or postprocessing. Our goal was to develop a widely applicable and robust method to assess regional mechanical dyssynchrony.

MATERIALS AND METHODS

Using the endocardial boundary, radial displacement curves (RDCs) were generated throughout the LV. Cross-correlation was used to determine the delay time between each RDC and a patient-specific reference. Delay times were projected onto the American Heart Association 17-segment model creating a regional dyssynchrony map. Our method was tested in 10 normal individuals and 10 patients enrolled for CRT (QRS > 120 ms, NYHA III-IV, EF < 35%).

RESULTS

Delay times over the LV were 23.9 ± 33.8 ms and 93.1 ± 99.9 ms (P < 0.001) in normal subjects and patients, respectively. Interobserver reproducibility for segment averages was 6.8 ± 39.3 ms and there was 70% agreement in identifying the latest contracting segment.

CONCLUSION

We have developed a method that can reliably calculate regional delay times from cine steady-state free-precession (SSFP) images. Maps of regional dyssynchrony could be used to identify the latest-contracting segment to assist in CRT lead implantation.

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.

Review of Journal of Cardiovascular Magnetic Resonance 2011

There were 83 articles published in the Journal of Cardiovascular Magnetic Resonance (JCMR) in 2011, which is an 11% increase in the number of articles since 2010. The quality of the submissions continues to increase. The editors had been delighted with the 2010 JCMR Impact Factor of 4.33, although this fell modestly to 3.72 for 2011. The impact factor undergoes natural variation according to citation rates of papers in the 2 years following publication, and is significantly influenced by highly cited papers such as official reports. However, we remain very pleased with the progress of the journal's impact over the last 5 years. Our acceptance rate is approximately 25%, and has been falling as the number of articles being submitted has been increasing. In accordance with Open-Access publishing, the JCMR articles go on-line as they are accepted with no collating of the articles into sections or special thematic issues. For this reason, the Editors feel it is useful to summarize the papers for the readership into broad areas of interest or theme, which we feel would be useful, so that areas of interest from the previous year can be reviewed in a single article in relation to each other and other recent JCMR articles. The papers are presented in broad themes and set in context with related literature and previously published JCMR papers to guide continuity of thought in the journal. We hope that you find the open-access system increases wider reading and citation of your papers, and that you will continue to send your quality manuscripts to JCMR for publication.

Optimized local infarct restraint improves left ventricular function and limits remodeling

BACKGROUND

Preventing expansion and dyskinetic movement of a myocardial infarction (MI) can limit left ventricular (LV) remodeling. Using a device designed to produce variable alteration of infarct stiffness and geometry, we sought to understand how these parameters affect LV function and remodeling early after MI.

METHODS

Ten pigs had posterolateral infarctions. An unexpanded device was placed in 5 animals at the time of infarction and 5 animals served as untreated controls. One week after MI animals underwent magnetic resonance imaging to assess LV size and regional function. In the treatment group, after initial imaging, the device was expanded with 2, 4, 6, 8, and 10 mL of saline. The optimal degree of inflation was defined as that which maximized stroke volume (SV). The device was left optimally inflated in the treatment animals for 3 additional weeks.

RESULTS

One week after MI, device inflation to 6 mL or greater significantly (p < 0.05) decreased end-systolic volume (0 mL, 59.9 mL ± 3.8; 6 mL, 54.0 mL ± 3.1; 8 mL, 50.5 mL ± 4.8; and 10 mL, 46.1 mL ± 2.2) and increased ejection fraction (EF) (0 mL, 0.346 ± 0.016; 6 mL, 0.0397 ± 0.009; 8 mL, 0.431 ± 0.027; and 10 mL, 0.441 ± 0.009). Systolic volume significantly (p < 0.05) improved for the 6 mL and 8 mL volumes (0 mL, 31.2 mL ± 2.6; 6 mL, 35.7 mL ± 2.0; and 8 mL, 37.5 mL ± 1.9) but trended downward for 10 mL (36.6 mL ± 2.8). At 4 weeks after MI, end-diastolic volume and end-systolic volume were unchanged from 1-week values in the treatment group while the control group continued to dilate. Systolic volume (38.2 ± 4.4 mL vs 34.0.1 ± 4.8 mL, p = 0.08) and EF (0.360 ± 0.026 vs 0.276 ± 0.014, p = 0.04) were also better in the treatment animals.

CONCLUSIONS

Optimized isolated infarct restraint can limit adverse LV remodeling after MI. The tested device affords the potential for a patient-specific therapy to preserve cardiac function after MI.

MR assessment of regional myocardial mechanics

Regional myocardial function can be measured by several MR techniques including tissue tagging, phase velocity mapping, and more recently, displacement encoding with stimulated echoes (DENSE) and strain encoding (SENC). Each of these techniques was developed separately and has undergone significant change since its original implementation. As a result, in the current literature, the common features and the differences between the techniques and what they measure are often unclear and confusing. This review article delivers an extensively referenced introductory text which clarifies the current methodology from the starting point of the Bloch equations. By doing this in a consistent way for each method, the similarities and differences between them are highlighted. In addition, their capabilities and limitations are discussed, together with their relative advantages and disadvantages. While the focus is on sequence design and development, the principal parameters measured by each technique are also summarized, together with brief results, with the reader being directed to the extensive literature on data processing and clinical applications for more detail.

Validation of continuously tagged MRI for the measurement of dynamic 3D skeletal muscle tissue deformation

PURPOSE

Typically spatial modulation of the magnetization (SPAMM) tagged magnetic resonance imaging (MRI) requires many repeated motion cycles limiting the applicability to highly repeatable tissue motions only. This paper describes the validation of a novel SPAMM tagged MRI and post-processing framework for the measurement of complex and dynamic 3D soft tissue deformation following just three motion cycles. Techniques are applied to indentation induced deformation measurement of the upper arm and a silicone gel phantom.

METHODS

A SPAMM tagged MRI methodology is presented allowing continuous (3.3-3.6 Hz) sampling of 3D dynamic soft tissue deformation using non segmented 3D acquisitions. The 3D deformation is reconstructed by the combination of three mutually orthogonal tagging directions, thus requiring only three repeated motion cycles. In addition a fully automatic post-processing framework is presented employing Gabor scale-space and filter-bank analysis for tag extrema segmentation and triangulated surface fitting aided by Gabor filter bank derived surface normals. Deformation is derived following tracking of tag surface triplet triangle intersections. The dynamic deformation measurements were validated using indentation tests (∼20 mm deep at 12 mm/s) on a silicone gel soft tissue phantom containing contrasting markers which provide a reference measure of deformation. In addition, the techniques were evaluated in vivo for dynamic skeletal muscle tissue deformation measurement during indentation of the biceps region of the upper arm in a volunteer.

RESULTS

For the phantom and volunteer tag point location precision were 44 and 92 μm, respectively resulting in individual displacements precisions of 61 and 91 μm, respectively. For both the phantom and volunteer data cumulative displacement measurement accuracy could be evaluated and the difference between initial and final locations showed a mean and standard deviation of 0.44 and 0.59 mm for the phantom and 0.40 and 0.73 mm for the human data. Finally accuracy of (cumulative) displacement was evaluated using marker tracking in the silicone gel phantom. Differences between true and predicted marker locations showed a mean of 0.35 mm and a standard deviation of 0.63 mm.

CONCLUSIONS

A novel SPAMM tagged MRI and fully automatic post-processing framework for the measurement of complex 3D dynamic soft tissue deformation following just three repeated motion cycles was presented. The techniques demonstrate dynamic measurement of complex 3D soft tissue deformation at subvoxel accuracy and precision and were validated for 3.3-3.6 Hz sampling of deformation speeds up to 12 mm/s.

Cardiac motion and deformation recovery from MRI: a review

Magnetic resonance imaging (MRI) is a highly advanced and sophisticated imaging modality for cardiac motion tracking and analysis, capable of providing 3D analysis of global and regional cardiac function with great accuracy and reproducibility. In the past few years, numerous efforts have been devoted to cardiac motion recovery and deformation analysis from MR image sequences. Many approaches have been proposed for tracking cardiac motion and for computing deformation parameters and mechanical properties of the heart from a variety of cardiac MR imaging techniques. In this paper, an updated and critical review of cardiac motion tracking methods including major references and those proposed in the past ten years is provided. The MR imaging and analysis techniques surveyed are based on cine MRI, tagged MRI, phase contrast MRI, DENSE, and SENC. This paper can serve as a tutorial for new researchers entering the field.

[Three-dimensional imaging in cardiac resynchronization therapy]

Cardiac resynchronization therapy improves clinical symptoms and prognosis of heart failure patients. However, it has been shown that up to 40% of patients do not respond to this therapy. Three main determinants of cardiac resynchronization therapy response have been identified: left ventricular dyssynchrony, left ventricular lead position, and extent and location of myocardial scar tissue. Two-dimensional echocardiography is the first imaging technique to evaluate patients who may be candidates for cardiac resynchronization therapy. However, a multimodality approach based on 3-dimensional imaging techniques may provide a more comprehensive evaluation of these patients by combining the assessment of the aforementioned pathophysiological determinants of cardiac resynchronization therapy response.

Review of journal of cardiovascular magnetic resonance 2010

There were 75 articles published in the Journal of Cardiovascular Magnetic Resonance (JCMR) in 2010, which is a 34% increase in the number of articles since 2009. The quality of the submissions continues to increase, and the editors were delighted with the recent announcement of the JCMR Impact Factor of 4.33 which showed a 90% increase since last year. Our acceptance rate is approximately 30%, but has been falling as the number of articles being submitted has been increasing. In accordance with Open-Access publishing, the JCMR articles go on-line as they are accepted with no collating of the articles into sections or special thematic issues. Last year for the first time, the Editors summarized the papers for the readership into broad areas of interest or theme, which we felt would be useful to practitioners of cardiovascular magnetic resonance (CMR) so that you could review areas of interest from the previous year in a single article in relation to each other and other recent JCMR articles 1. This experiment proved very popular with a very high rate of downloading, and therefore we intend to continue this review annually. The papers are presented in themes and comparison is drawn with previously published JCMR papers to identify the continuity of thought and publication in the journal. We hope that you find the open-access system increases wider reading and citation of your papers, and that you will continue to send your quality manuscripts to JCMR for publication.

Regional left ventricular myocardial contractility and stress in a finite element model of posterobasal myocardial infarction

Recently, a noninvasive method for determining regional myocardial contractility, using an animal-specific finite element (FE) model-based optimization, was developed to study a sheep with anteroapical infarction (Sun et al., 2009, "A Computationally Efficient Formal Optimization of Regional Myocardial Contractility in a Sheep With Left Ventricular Aneurysm," ASME J. Biomech. Eng., 131(11), p. 111001). Using the methodology developed in the previous study (Sun et al., 2009, "A Computationally Efficient Formal Optimization of Regional Myocardial Contractility in a Sheep With Left Ventricular Aneurysm," ASME J. Biomech. Eng., 131(11), p. 111001), which incorporates tagged magnetic resonance images, three-dimensional myocardial strains, left ventricular (LV) volumes, and LV cardiac catheterization pressures, the regional myocardial contractility and stress distribution of a sheep with posterobasal infarction were investigated. Active material parameters in the noninfarcted border zone (BZ) myocardium adjacent to the infarct (T(max_B)), in the myocardium remote from the infarct (T(max_R)), and in the infarct (T(max_I)) were estimated by minimizing the errors between FE model-predicted and experimentally measured systolic strains and LV volumes using the previously developed optimization scheme. The optimized T(max_B) was found to be significantly depressed relative to T(max_R), while T(max_I) was found to be zero. The myofiber stress in the BZ was found to be elevated, relative to the remote region. This could cause further damage to the contracting myocytes, leading to heart failure.

Validation of SPAMM tagged MRI based measurement of 3D soft tissue deformation

PURPOSE

This study presents and validates a novel (non-ECG-triggered) MRI sequence based on spatial modulation of the magnetization (SPAMM) to noninvasively measure 3D (quasistatic) soft tissue deformations using only six acquisitions (three static and three indentations). In the current SPAMM tagged MRI approaches, data are typically constructed from many repeated motion cycles. This has so far restricted its application to the measurement of highly repeatable and periodic movements (e.g., cardiac deformation). In biomechanical applications where soft tissue deformation is artificially induced, often by indentation, significant repeatability constraints exist, and for clinical applications, discomfort and health issues generally preclude a large number of repetitions.

METHODS

A novel (non-ECG-triggered) SPAMM tagged MRI sequence is presented, whereby a single 1-1 (first order) SPAMM set is acquired following a 3D transient field echo acquisition. Full 3D deformation measurement is achieved through the combination of only six acquisitions (three static and three motion cycles). The 3D deformation measurements were validated using quasistatic indentation tests and marker tracking in a silicone gel soft tissue phantom. In addition, the technique's ability to measure 3D soft tissue deformation in vivo was evaluated using indentation of the biceps region of the upper arm in a volunteer.

RESULTS

Following comparison to marker tracking in the silicone gel phantom, the SPAMM tagged MRI based displacement measurement demonstrated subvoxel accuracy with a mean displacement difference of 72 microm and a standard deviation of 289 microm. In addition, precision of displacement magnitude was evaluated for both the phantom and the volunteer data. The standard deviations of the displacement magnitude with respect to the average displacement magnitude were 75 and 169 microm for the phantom and volunteer data, respectively.

CONCLUSIONS

The subvoxel accuracy and precision demonstrated in the phantom in combination with the precision comparison between the phantom and the volunteer data provide confidence in the methods presented for measurement of soft tissue deformation in vivo. To the author's knowledge, since only six acquisitions are required, the presented methodology is the fastest SPAMM tagged MRI method currently available for the noninvasive measurement of quasistatic 3D soft tissue deformation.

On the mechanisms enabling myocardial edema contrast in bSSFP-based imaging approaches

The biophysical mechanisms influencing balanced steady-state free precession (bSSFP) based edema imaging in the setting of acute myocardial infarction are not well understood. To assess the various mechanisms that enable the detection of myocardial edema on bSSFP-based imaging approaches (cine bSSFP and T(2)-prepared bSSFP), experiments were conducted in canine models subjected to ischemia-reperfusion injury. Results showed that in addition to relaxation effects, the alteration in thermal equilibrium (M(0)) (including magnetization transfer) has a significant contribution to the image contrast between edematous and healthy myocardium. The relative signal-intensity ratios between edematous and healthy myocardium were: 1.51 ± 0.18 (cine bSSFP) and 1.58 ± 0.20 (T(2)-prepared bSSFP); the theoretically estimated relative relaxation and M(0) effects were: 1.17 ± 0.09 and 1.30 ± 0.19, respectively (cine bSSFP), and 1.49 ± 0.23 and 1.06 ± 0.07, respectively (T(2)-prepared bSSFP). There were no significant difference between cine bSSFP and T(2)-prep bSSFP relative signal-intensity ratios. However, the relative relaxation effect in cine bSSFP was significantly lower than in T(2)-prep bSSFP (P < 0.05), and the M(0) effect in cine bSSFP was significantly higher than in T(2)-prep bSSFP (P < 0.05). Hence the acquisition strategies that wish to maximize myocardial edema contrast in cine bSSFP imaging should take both relaxation and M(0) effects into account.

Cardiac resynchronization therapy guided by cardiovascular magnetic resonance

Cardiac resynchronization therapy (CRT) is an established treatment for patients with symptomatic heart failure, severely impaired left ventricular (LV) systolic dysfunction and a wide (> 120 ms) complex. As with any other treatment, the response to CRT is variable. The degree of pre-implant mechanical dyssynchrony, scar burden and scar localization to the vicinity of the LV pacing stimulus are known to influence response and outcome. In addition to its recognized role in the assessment of LV structure and function as well as myocardial scar, cardiovascular magnetic resonance (CMR) can be used to quantify global and regional LV dyssynchrony. This review focuses on the role of CMR in the assessment of patients undergoing CRT, with emphasis on risk stratification and LV lead deployment.