Vibration-synchronized magnetic resonance imaging for the detection of myocardial elasticity changes

Cardiovascular magnetic resonance elastography: A review

Cardiovascular diseases are the leading cause of death worldwide. These cardiovascular diseases are associated with mechanical changes in the myocardium and aorta. It is known that stiffness is altered in many diseases, including the spectrum of ischemia, diastolic dysfunction, hypertension and hypertrophic cardiomyopathy. In addition, the stiffness of the aortic wall is altered in multiple diseases, including hypertension, coronary artery disease and aortic aneurysm formation. For example, in diastolic dysfunction in which the ejection fraction is preserved, stiffness can potentially be an important biomarker. Similarly, in aortic aneurysms, stiffness can provide valuable information with regard to rupture potential. A number of studies have addressed invasive and non-invasive approaches to test and measure the mechanical properties of the myocardium and aorta. One of the non-invasive approaches is magnetic resonance elastography (MRE). MRE is a phase-contrast magnetic resonance imaging technique that measures tissue stiffness non-invasively. This review article highlights the technical details and application of MRE in the quantification of myocardial and aortic stiffness in different disease states.

Relative identifiability of anisotropic properties from magnetic resonance elastography

Although magnetic resonance elastography (MRE) has been used to estimate isotropic stiffness in the heart, myocardium is known to have anisotropic properties. This study investigated the determinability of global transversely isotropic material parameters using MRE and finite-element modeling (FEM). A FEM-based material parameter identification method, using a displacement-matching objective function, was evaluated in a gel phantom and simulations of a left ventricular (LV) geometry with a histology-derived fiber field. Material parameter estimation was performed in the presence of Gaussian noise. Parameter sweeps were analyzed and characteristics of the Hessian matrix at the optimal solution were used to evaluate the determinability of each constitutive parameter. Four out of five material stiffness parameters (Young's modulii E1 and E3 , shear modulus G13 and damping coefficient s), which describe a transversely isotropic linear elastic material, were well determined from the MRE displacement field using an iterative FEM inversion method. However, the remaining parameter, Poisson's ratio, was less identifiable. In conclusion, Young's modulii, shear modulii and damping can theoretically be well determined from MRE data, but Poisson's ratio is not as well determined and could be set to a reasonable value for biological tissue (close to 0.5).

Estimation of transversely isotropic material properties from magnetic resonance elastography using the optimised virtual fields method

Magnetic resonance elastography (MRE) has been used to estimate isotropic myocardial stiffness. However, anisotropic stiffness estimates may give insight into structural changes that occur in the myocardium as a result of pathologies such as diastolic heart failure. The virtual fields method (VFM) has been proposed for estimating material stiffness from image data. This study applied the optimised VFM to identify transversely isotropic material properties from both simulated harmonic displacements in a left ventricular (LV) model with a fibre field measured from histology as well as isotropic phantom MRE data. Two material model formulations were implemented, estimating either 3 or 5 material properties. The 3-parameter formulation writes the transversely isotropic constitutive relation in a way that dissociates the bulk modulus from other parameters. Accurate identification of transversely isotropic material properties in the LV model was shown to be dependent on the loading condition applied, amount of Gaussian noise in the signal, and frequency of excitation. Parameter sensitivity values showed that shear moduli are less sensitive to noise than the other parameters. This preliminary investigation showed the feasibility and limitations of using the VFM to identify transversely isotropic material properties from MRE images of a phantom as well as simulated harmonic displacements in an LV geometry.

In vivo magnetic resonance elastography to estimate left ventricular stiffness in a myocardial infarction induced porcine model

PURPOSE

To estimate change in left ventricular (LV) end-systolic and end-diastolic myocardial stiffness (MS) in pigs induced with myocardial infarction (MI) with disease progression using cardiac magnetic resonance elastography (MRE) and to compare it against ex vivo mechanical testing, LV circumferential strain, and magnetic resonance imaging (MRI) relaxometry parameters (T₁ , T₂ , and extracellular volume fraction [ECV]).

MATERIALS AND METHODS

MRI (1.5T) was performed on seven pigs, before surgery (Bx), and 10 (D10), and 21 (D21) days after creating MI. Cardiac MRE-derived MS was measured in infarcted region (MIR) and remote region (RR), and validated against mechanical testing-derived MS obtained postsacrifice on D21. Circumferential strain and MRI relaxometry parameters (T₂ , T₁ , and ECV) were also obtained. Multiparametric analysis was performed to determine correlation between cardiac MRE-derived MS and 1) strain, 2) relaxometry parameters, and 3) mechanical testing.

RESULTS

Mean diastolic (D10: 5.09 ± 0.6 kPa; D21: 5.45 ± 0.7 kPa) and systolic (D10: 5.72 ± 0.8 kPa; D21: 6.34 ± 1.0 kPa) MS in MIR were significantly higher (P < 0.01) compared to mean diastolic (D10: 3.97 ± 0.4 kPa; D21: 4.12 ± 0.2 kPa) and systolic (D10: 5.08 ± 0.6 kPa; and D21: 5.16 ± 0.6 kPa) MS in RR. The increase in cardiac MRE-derived MS at D21 (MIR) was consistent and correlated strongly with mechanical testing-derived MS (r(diastolic) = 0.86; r(systolic) = 0.89). Diastolic MS in MIR demonstrated a negative correlation with strain (r = 0.58). Additionally, cardiac MRE-derived MS demonstrated good correlations with post-contrast T₁ (r(diastolic) = -0.549; r(systolic) = -0.741) and ECV (r(diastolic) = 0.548; r(systolic) = 0.703), and no correlation with T₂ .

CONCLUSION

As MI progressed, cardiac MRE-derived MS increased in MIR compared to RR, which significantly correlated with mechanical testing-derived MS, T₁ and ECV.

LEVEL OF EVIDENCE

1 J. Magn. Reson. Imaging 2017;45:1024-1033.

Investigation of the effects of myocardial anisotropy for shear wave elastography using impulsive force and harmonic vibration

The myocardium is known to be an anisotropic medium where the muscle fiber orientation changes through the thickness of the wall. Shear wave elastography methods use propagating waves which are measured by ultrasound or magnetic resonance imaging (MRI) techniques to characterize the mechanical properties of various tissues. Ultrasound- or MR-based methods have been used and the excitation frequency ranges for these various methods cover a large range from 24-500 Hz. Some of the ultrasound-based methods have been shown to be able to estimate the fiber direction. We constructed a model with layers of elastic, transversely isotropic materials that were oriented at different angles to simulate the heart wall in systole and diastole. We investigated the effect of frequency on the wave propagation and the estimation of fiber direction and wave speeds in the different layers of the assembled models. We found that waves propagating at low frequencies such as 30 or 50 Hz showed low sensitivity to the fiber direction but also had substantial bias in estimating the wave speeds in the layers. Using waves with higher frequency content (>200 Hz) allowed for more accurate fiber direction and wave speed estimation. These results have particular relevance for MR- and ultrasound-based elastography applications in the heart.

Measuring age-dependent myocardial stiffness across the cardiac cycle using MR elastography: A reproducibility study

PURPOSE

To assess reproducibility in measuring left ventricular (LV) myocardial stiffness in volunteers throughout the cardiac cycle using MR elastography (MRE) and to determine its correlation with age.

METHODS

Cardiac MRE (CMRE) was performed on 29 normal volunteers, with ages ranging from 21 to 73 years. For assessing reproducibility of CMRE-derived stiffness measurements, scans were repeated per volunteer. Wave images were acquired throughout the LV myocardium, and were analyzed to obtain mean stiffness during the cardiac cycle. CMRE-derived stiffness values were correlated to age.

RESULTS

Concordance correlation coefficient revealed good interscan agreement with rc of 0.77, with P-value < 0.0001. Significantly higher myocardial stiffness was observed during end-systole (ES) compared with end-diastole (ED) across all subjects. Additionally, increased deviation between ES and ED stiffness was observed with increased age.

CONCLUSION

CMRE-derived stiffness is reproducible, with myocardial stiffness changing cyclically across the cardiac cycle. Stiffness is significantly higher during ES compared with ED. With age, ES myocardial stiffness increases more than ED, giving rise to an increased deviation between the two.

Shear-wave amplitudes measured with cardiac MR elastography for diagnosis of diastolic dysfunction

PURPOSE

To test whether shear-wave amplitudes (SWAs) in the myocardium measured with cardiac magnetic resonance (MR) elastography enable diagnosis of myocardial relaxation abnormalities in patients with diastolic dysfunction.

MATERIALS AND METHODS

Each subject gave written informed consent to participate in this institutional review board-approved prospective study. Electrocardiographically triggered SWA-based cardiac MR elastography with 24.13-Hz external vibration frequency was performed in 50 subjects grouped into asymptomatic young (n = 10, 18-39 years) and asymptomatic old (n = 10, 40-68 years) subjects and patients with echocardiographically proved mild, moderate, or severe diastolic dysfunction (n = 30, 44-73 years). SWA images were analyzed in the left ventricular (LV) region and were normalized against reference SWA of the thoracic wall. Analysis of variance with Bonferroni-corrected pairwise comparison and Pearson correlation were used for statistical evaluation.

RESULTS

Young and old control subjects had normalized mean LV SWA of 0.67 ± 0.04 (standard error of the mean) and 0.56 ± 0.04 (P = .18, F test), respectively. Compared with the control groups, patients with mild, moderate, and severe diastolic dysfunction displayed significantly reduced normalized mean LV SWA of 0.37 ± 0.04, 0.34 ± 0.04, and 0.29 ± 0.04 (P < .001, F test), respectively, which was inversely correlated to the severity of diastolic dysfunction (R = -0.61, P < .001). The best cutoff value to differentiate between asymptomatic volunteers and patients was 0.43, yielding an area under the receiver operating characteristic curve of 0.92, with 90% sensitivity and 89.7% specificity.

CONCLUSION

LV SWA measured with cardiac MR elastography provides image contrast sensitive to myocardial relaxation abnormalities and shows significantly lower values in patients with diastolic dysfunction.

Isovolumetric elasticity alteration in the human heart detected by in vivo time-harmonic elastography

Time harmonic elastography (THE) has recently been introduced for measurement of the periodic alteration in myocardial shear modulus based on externally induced low-frequency acoustic vibrations produced by a loudspeaker. In this study, we propose further developments of cardiac THE toward a clinical modality including integration of the vibration source into the patient bed and automated parameter extraction from harmonic shear wave amplitudes, wall motion profiles and synchronized electrocardiographic records. This method has enabled us to evaluate the delay between wall motion and wave amplitude alteration for the measurement of isovolumetric times of elasticity alteration during contraction (τ(C)) and relaxation (τ(R)) in a group of 32 healthy volunteers. On average, the wave amplitudes changed between systole and diastole by a factor of 1.7 ± 0.3, with a τ(C) of 137 ± 61 ms and a τ(R) of 68 ± 73 ms, which agrees with results obtained with the more time-consuming and expensive cardiac magnetic resonance elastography. Furthermore, because of the high sampling rate, elasto-morphometric parameters such as transition times and the area of wave amplitude-cardiac motion cycles can be processed in an automated way for the future clinical detection of myocardial relaxation abnormalities.

Rapid acquisition of multifrequency, multislice and multidirectional MR elastography data with a fractionally encoded gradient echo sequence

In MR elastography (MRE), periodic tissue motion is phase encoded using motion-encoding gradients synchronized to an externally applied periodic mechanical excitation. Conventional methods result in extended scan time for quality phase images, thus limiting the broad application of MRE in the clinic. For practical scan times, researchers have been relying on one-dimensional or two-dimensional motion-encoding, low-phase sampling and a limited number of slices, and artifact-prone, single-shot, echo planar imaging (EPI) readout. Here, we introduce a rapid multislice pulse sequence capable of three-dimensional motion encoding that is also suitable for simultaneously encoding motion with multiple frequency components. This sequence is based on a gradient-recalled echo (GRE) sequence and exploits the principles of fractional encoding. This GRE MRE pulse sequence was validated as capable of acquiring full three-dimensional motion encoding of isotropic voxels in a large volume within less than a minute. This sequence is suitable for monofrequency and multifrequency MRE experiments. In homogeneous paraffin phantoms, the eXpresso sequence yielded similar storage modulus values as those obtained with conventional methods, although with markedly reduced variances (7.11 ± 0.26 kPa for GRE MRE versus 7.16 ± 1.33 kPa for the conventional spin-echo EPI sequence). The GRE MRE sequence obtained better phase-to-noise ratios than the equivalent spin-echo EPI sequence (matched for identical acquisition time) in both paraffin phantoms and in vivo data in the liver (59.62 ± 11.89 versus 27.86 ± 3.81, 61.49 ± 14.16 versus 24.78 ± 2.48 and 58.23 ± 10.39 versus 23.48 ± 2.91 in the X, Y and Z components, respectively, in the case of liver experiments). Phase-to-noise ratios were similar between GRE MRE used in monofrequency or multifrequency experiments (75.39 ± 14.93 versus 86.13 ± 18.25 at 28 Hz, 71.52 ± 24.74 versus 86.96 ± 30.53 at 56 Hz and 95.60 ± 36.96 versus 61.35 ± 26.25 at 84Hz, respectively).