Elasticity-based determination of isovolumetric phases in the human heart
© Elgeti et al; licensee BioMed Central Ltd. 2010
Received: 24 March 2010
Accepted: 27 October 2010
Published: 27 October 2010
To directly determine isovolumetric cardiac time intervals by magnetic resonance elastography (MRE) using the magnitude of the complex signal for deducing morphological information combined with the phase of the complex signal for tension-relaxation measurements.
Thirty-five healthy volunteers and 11 patients with relaxation abnormalities were subjected to transthoracic wave stimulation using vibrations of approximately 25 Hz. A k-space-segmented, ECG-gated gradient-recalled echo steady-state sequence with a 500-Hz bipolar motion-encoding gradient was used for acquiring a series of 360 complex images of a short-axis view of the heart at a frame rate of less than 5.2 ms. Magnitude images were employed for measuring the cross-sectional area of the left ventricle, while phase images were used for analyzing the amplitudes of the externally induced waves. The delay between the decrease in amplitude and onset of ventricular contraction was determined in all subjects and assigned to the time of isovolumetric tension. Conversely, the delay between the increase in wave amplitude and ventricular dilatation was used for measuring the time of isovolumetric elasticity relaxation.
Wave amplitudes decreased during systole and increased during diastole. The variation in wave amplitude occurred ahead of morphological changes. In healthy volunteers the time of isovolumetric elasticity relaxation was 75 ± 31 ms, which is significantly shorter than the time of isovolumetric tension of 136 ± 36 ms (P < 0.01). In patients with relaxation abnormalities (mild diastolic dysfunction, n = 11) isovolumetric elasticity relaxation was significantly prolonged, with 133 ± 57 ms (P < 0.01), whereas isovolumetric tension time was in the range of healthy controls (161 ± 45 ms; P = 0.053).
The complex MRE signal conveys complementary information on cardiac morphology and elasticity, which can be combined for directly measuring isovolumetric tension and elasticity relaxation in the human heart.
Times of isovolumetric contraction and relaxation and the resulting indices (e.g., Tei index) are frequently used to describe overall cardiac function . The standard method for measuring cardiac time intervals is echocardiography, which exploits the opening and closing times of aortic and mitral valves . Isovolumetric cardiac phases are characterized by changes in myocardial elasticity, while ventricular volume remains constant. Thus, the most direct determination of isovolumetric times combines measurement of both ventricular volume and myocardial elasticity. Cardiac volumes can be measured using ultrasound or magnetic resonance (MR) imaging. Myocardial elasticity can be measured noninvasively by elastography using either ultrasound or MR imaging and various mechanical stimuli such as time-harmonic vibrations [3–6], focused ultrasound pulses [7, 8] or transient intrinsic cardiac waves . In general, dynamic elastography relies on the evaluation of shear waves for recovering mechanical tissue parameters from amplitudes, wave lengths, or propagation speed. Classically, the wave equation is solved for elastic moduli using for example image-based algebraic Helmholtz inversion  or profile-based phase-gradient methods. Inversion techniques are well applicable in large organs such as the liver, brain, or larger groups of skeletal muscle. However, the complex anatomy of the heart requires consideration of finite geometries as proposed by Kolipaka et al.  and demonstrated by MR elastography (MRE) on a spherical rubber phantom. This approach assumes a lossless, isotropic, and linear elastic heart with a spherical thin-shell geometry that is vibrated solely by in-plane harmonic shear motion. Fewer assumptions are made in cardiac MRE described by Sack et al. , where no wave inversion is needed, and relative myocardial elasticity changes are deducible across the cardiac cycle. This wave-amplitude-sensitive MRE technique cannot measure absolute elasticity values; however, relative elasticity ratios are provided with an excellent time resolution of approximately 5 ms, which we will exploit in this study for deriving elasticity-based isovolumetric phases in the human heart. We will evaluate solely the relative timing between heart geometry and wave amplitudes based on one important observation made previously in cardiac MRE on healthy volunteers and animals: there is a distinct delay in the dynamics of wave amplitudes and ventricular volume such that the alteration of wave amplitudes begins 75 to 160 milliseconds earlier than changes in heart geometry [4, 12, 13]. Analyzing this timing we aim to develop wave amplitudes as a surrogate marker for alterations in myocardial elasticity during the cardiac cycle. We do not account for mechanical properties beyond shear elasticity such as viscosity, compressibility, and anisotropy or nonlinearity of elasticity. Therefore all these properties will be subsumed under the label of apparent elasticity, and the terms myocardial tension and relaxation will respectively refer to the increase and decrease in apparent elasticity.
In the following, the dynamics of ventricular geometry is estimated from the cross-sectional area of the left ventricle in a short-axis view of the heart. Wave amplitudes and cross-sectional areas are derived from the phase and the magnitudes of the complex MRE signal, respectively. Thus, all temporal information on the ventricular elasticity-volume cycle is available from the MRE signal. The pathophysiological significance of this information is demonstrated by applying the method to healthy volunteers and to patients with mild left ventricular relaxation abnormalities. Our hypothesis was that impaired ventricular relaxation can be directly deduced from a longer delay between the increase in wave amplitude and left ventricular dilatation at early diastole.
This study was approved by the local ethics committee (EA 1/055/07-1-4), and written informed consent was obtained from all subjects. The group of healthy volunteers (n = 35, aged between 18 and 59 years, mean 35 ± 8.5; 4 women) had no history of cardiac disease. The group of patients (n = 11, aged 47-70 years, mean 60 ± 8.5; 2 women) had echocardiographically proven relaxation abnormalities (E/A < 0.75 and E/E' = 8-12).
Cardiac MRE was performed using clinical 1.5 T scanners (Siemens Magnetom Sonata and Avanto) and 16-channel phased-array chest coils. An ECG-gated gradient-recalled echo steady-state sequence (repetition time, TR, 5.18 ms; time to echo, TE, 3.29 ms; flip angle, α,25°; twofold GRAPPA acceleration; 128 × 96 matrix; typical FoV 320 × 250 mm; 7 mm slice thickness; 48 single-line k-space segments corresponding to 48 phase-encoding steps; 1.3 kHz bandwidth) was sensitized to motion by applying a sinusoidal motion-encoding gradient (MEG) of 2 ms duration, corresponding to an MEG frequency of 500 Hz. Three MRE experiments were performed on each subject while consecutively alternating the MEG orientation between read-out, phase encoding, and slice selection. The amplitude of the MEG was 25 mT/m in each direction. The start of image acquisition was controlled by the R-wave of the ECG. A total of 360 separate images were acquired across approximately two cardiac cycles (1.86 seconds = TR × 360). During this interval, a single phase-encoding step was performed for each image. Acquisition of all k-space data for each image required 48 × 1.86 seconds. During each of the 1.86-second intervals of data acquisition, the patient held their breath in expiration, followed by a short period of 2.5 seconds of free breathing (one inhalation and one exhalation) before the next k-line acquisition started.
A customer-built actuator produced low-frequency acoustic vibrations, which were applied to the patient's chest using a rigid piston . A sinusoidal burst of 24.13 Hz frequency (1/8TR) was fed into the actuator after every eight TRs, resulting in 45 (360/8) continuous harmonic vibration cycles for each of the 48 phase-encoding steps. Depending on heart rate, one MRE scan for each vibration direction took approximately 3.5-4.0 min. Image slice orientation was set to a short cardiac axis in the anterior third of the left ventricle.
MRE data evaluation
Images were reconstructed from 96 single-line k-space segments after 48 phase-encoding steps (128 × 96 matrix, 75% phase sampling, and GRAPPA factor of 2 with 12 extra lines for image reconstruction). The magnitude and phase of the complex MRE signal, M(x,t) and φ(x,t), were exploited for analyzing the change in left ventricular (LV) cross-sectional areas, a LV (t), and wave amplitudes, U(t), respectively. a LV (t) was determined by manual segmentation of anteroposterior and septal-lateral diameters from 90 magnitude images averaged over four consecutive images M(x, nTR), M(x, [n+1]TR), M(x, [n+2]TR), and M(x, [n+3]TR) with n = 1,5,..360. Averaging was performed since manual segmentation of a LV required discrimination of morphological features of the heart based on the contrast of MRE magnitude images. As this contrast was inherently impaired due to the short TR, this type of image processing improved image quality. To further reduce physiological noise related to high k-space segmentation was averaged across MRE experiments of varying motion component sensitivity. The resulting time resolution of and a LV (t) was 4TR = 20.7 ms. Characteristic times of the cardiac cycle given by the duration of diastole (τdia) and systole (τsys) were determined from the a LV (t) function. Furthermore, the shortening fraction was calculated from the difference in a LV between end-systole and end-diastole, i.e., (a LV [end-diastole] - aLV [end-systole)])/a LV (end-diastole), corresponding to the shortening fraction derived from echocardiography .
Calculation of the unwrapped phase velocity with , as proposed in .
ii) Transformation of the resulting phase velocity to deflections with , where ω is the angular driving frequency. q (in the dimension of microns per rad) denotes the inverse encoding efficiency that scales phases to wave deflections u .
Data is given as mean ± standard deviation. Student's t-test was used to test the time intervals within each group for statistical differences. The two-sample Kolmogorov-Smirnov test was used to test the time intervals between volunteers and patients for statistical differences. Pearson's correlation test was performed to test for possible linear correlation of the measured variables. A P-value of < 0.05 was considered statistically significant.
In all experiments, wave amplitudes decreased due to myocardial contraction. The change in wave amplitude occurred before that in the cross-sectional area of the left ventricle.
τ A (ms)
136 ± 36
161 ± 46
P = 0.053, k = 0.44
τ A ′ (ms)
138 ± 37
173 ± 50
P = 0.06, k = 0.42
τ sys (ms)
297 ± 64
365 ± 87
P < 0.001, k = 1
τ A /τ sys
0.25 ± 0.06
0.36 ± 0.13
P = 0.01, k = 0.54
τ B (ms)
75 ± 31
133 ± 58
P = 0.005, k = 0.55
τ B ′ (ms)
76 ± 30
142 ± 68
P = 0.006, k = 0.52
τ dia (ms)
682 ± 108
517 ± 91
P = 0.0035, k = 0.61
RR interval (ms)
982 ± 102
888 ± 139
P = 0.08, k = 0.40
56.0 ± 5.5
44.8 ± 12.6
P = 0.008, k = 0.53
Cardiac time intervals give insight into the hemodynamics and mechanical function of the heart and are thus valuable indicators of cardiac health. This study for the first time analyzed cardiac intervals based on changes in elastic wave amplitudes in volunteers and in patients with mild diastolic dysfunction. The evaluation of the complex signal in MRE enabled us to derive the timing of elasticity variation in relation to that of changes in LV geometry from a single MRE examination. This paves the way for the development of software capable of automatically combining magnitude-based and phase-based information from cardiac MRE. On-site implementation would provide a new modality of cardiac MRI for the instantaneous assessment of isovolumetric times based on changes in myocardial elasticity. In this study, externally induced elastic waves were used for probing the inherent cohesion of the mechanical matrix of myocardium, which determines its macroscopic shear elasticity. Changes in shear elasticity of biological tissue are known to be sensitive to disease progression [18–20]. In particular, cardiac function is synonymous with variation in myocardial elasticity (i.e. shear modulus). A variety of heart diseases are associated with dysfunctional elasticity variations of the myocardium. An important cardiac disease pattern is diastolic dysfunction, which is characterized by an elevated elastic modulus of the ventricular wall muscle at relaxation [21–23]. This has driven the development of imaging modalities sensitive to shear modulus variations in the heart [3–9].
It is important to mention that the cardiac time intervals, τ A and τ B , analyzed in this study are not directly comparable to isovolumetric times measured by echocardiography or other more traditional methods relying on auscultation, ECG, phonocardiogram, or carotid pulse tracing . Morphology-based methods allow us to assess the dynamics of heart geometry such as valve position, ventricular volumes, tissue strain, and blood flow or tissue velocity. In M-mode echocardiography the isovolumetric time intervals have been calculated from the opening and closure times of aortic and mitral valve relative to the R-wave of the ECG . However, the well-known electromechanical delay limits the reliability of the time intervals deduced in this way. In Doppler echocardiography time intervals can be determined more reliably using the velocity of tissue motion [2, 26]. Since both methods evaluate the change in cardiac geometry, no elasticity-based information is obtained. Hence, the time constants determined with these methods are different from τ A and τ B proposed in this study.
Isovolumetric relaxation time (IVR) measured by tissue doppler imaging is 70 ± 22 ms (healthy volunteers 30-40 years of age) , showing good agreement with the time of elasticity relaxation, τ B , measured in our study (75 ± 31 ms). It is known that this time interval is prolonged in mild diastolic dysfunction . Surprisingly, the MRE-derived time of tension, τ A , in both groups investigated in our study (136 ± 36 ms and 161 ± 46 ms) is considerable longer than the echocardiographic isovolumetric contraction time of 63 ± 14 ms . This might be attributable to the fact that the elastic properties of muscle tissue are nonlinear under large deformation, rendering the measured MRE wave response stretch-dependent . Thus, the nonlinearity of myocardium results in an increase in elasticity during ventricular expansion even though the muscle is intrinsically relaxed. Upon electrical activation, the hyperelastic tension of stretched myocardium is superposed with active contraction, while the stretch-related hyperelasticity of the muscle decreases. U(t) might therefore decrease prior to systolic contraction, which would cause a longer τ A time than observed by echocardiography. Since our τ A and τ B parameters reflect the full stress-strain hysteresis of myocardium, they represent a novel measure of cardiac dynamics, whose clinical significance remains to be investigated in a large group of patients. In this study, the diagnostic significance of elasticity-based cardiac time intervals is suggested by the findings in 11 patients suffering from mild diastolic dysfunction. The fact that in this small group of patients, elasticity relaxation of the left ventricle was significantly slowed down (indicated by increases in τ B and τ A /τ sys ) provides further motivation for improving the applicability and accuracy of cardiac MRE and for determining its diagnostic relevance.
It is known that isovolumetric relaxation time is, to a certain amount, age dependent, increasing from 74 ± 7 ms between 41 and 60 years to 87 ± 7 ms after age 60 , whereas the isovolumetric contraction time remains widely unchanged . Therefore, further investigations are necessary to study the influence of age on our results. It is also known that, with severe diastolic dysfunction, isovolumetric relaxation time shortens, so that this parameter cannot diagnose severe diastolic dysfunction in the absence of further information .
Multiple factors can influence the wave amplitudes measured in cardiac MRE. We cannot discriminate shear elasticity from viscosity effects, from the reorientation of myocardial fibers that is associated with a change in the apparent elasticity (see discussion in ) or nonlinear effects as discussed above. A reciprocal correlation of U(t) with ventricular pressure was shown in . Any further interpretation of amplitude dynamics in terms of myocardial shear elasticity is still model-based and needs experimental validation.
Cardiac MRE provides temporally resolved information on both morphology and elasticity dynamics of the heart, which can be combined for determining isovolumetric times. The resulting elasticity-based timing parameters comprise tissue constitution and mechanical functioning of the heart and are thus relevant for the assessment of cardiac health. Future cardiac MRE experiments should include the evaluation of phase and magnitude information in order to further evaluate the extent to which time intervals in MRE can provide useful diagnostic information.
- Tei C, Ling LH, Hodge DO, Bailey KR, Oh JK, Rodeheffer RJ, Tajik AJ, Seward JB: New index of combined systolic and diastolic myocardial performance: a simple and reproducible measure of cardiac function--a study in normals and dilated cardiomyopathy. J Cardiol. 1995, 26: 357-366.PubMedGoogle Scholar
- Cui W, Roberson DA, Chen Z, Madronero LF, Cuneo BF: Systolic and diastolic time intervals measured from Doppler tissue imaging: normal values and Z-score tables, and effects of age, heart rate, and body surface area. J Am Soc Echocardiogr. 2008, 21: 361-370. 10.1016/j.echo.2007.05.034.View ArticlePubMedGoogle Scholar
- Rump J, Klatt D, Braun J, Warmuth C, Sack I: Fractional encoding of harmonic motions in MR elastography. Magn Reson Med. 2007, 57: 388-395. 10.1002/mrm.21152.View ArticlePubMedGoogle Scholar
- Sack I, Rump J, Elgeti T, Samani A, Braun J: MR elastography of the human heart: noninvasive assessment of myocardial elasticity changes by shear wave amplitude variations. Magn Reson Med. 2009, 61: 668-677. 10.1002/mrm.21878.View ArticlePubMedGoogle Scholar
- Kolipaka A, Mcgee KP, Araoz PA, Glaser KJ, Manduca A, Romano AJ, Ehman RL: MR elastography as a method for the assessment of myocardial stiffness: comparison with an established pressure-volume model in a left ventricular model of the heart. Magn Reson Med. 2009, 62: 135-140. 10.1002/mrm.21991.PubMed CentralView ArticlePubMedGoogle Scholar
- Robert B, Sinkus R, Gennisson JL, Fink M: Application of DENSE-MR-elastography to the human heart. Magn Reson Med. 2009, 62: 1155-1163. 10.1002/mrm.22124.View ArticlePubMedGoogle Scholar
- Bouchard RR, Hsu SJ, Wolf PD, Trahey GE: In vivo cardiac, acoustic-radiation-force-driven, shear wave velocimetry. Ultrason Imaging. 2009, 31: 201-213.PubMed CentralView ArticlePubMedGoogle Scholar
- Hsu SJ, Bouchard RR, Dumont DM, Ong CW, Wolf PD, Trahey GE: Novel acoustic radiation force impulse imaging methods for visualization of rapidly moving tissue. Ultrason Imaging. 2009, 31: 183-200.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanai H: Propagation of spontaneously actuated pulsive vibration in human heart wall and in vivo viscoelasticity estimation. IEEE Trans Ultrason Ferroelectr Freq Control. 2005, 52: 1931-1942. 10.1109/TUFFC.2005.1561662.View ArticlePubMedGoogle Scholar
- Oliphant TE, Manduca A, Ehman RL, Greenleaf JF: Complex-valued stiffness reconstruction for magnetic resonance elastography by algebraic inversion of the differential equation. Magn Reson Med. 2001, 45: 299-310. 10.1002/1522-2594(200102)45:2<299::AID-MRM1039>3.0.CO;2-O.View ArticlePubMedGoogle Scholar
- Kolipaka A, McGee KP, Araoz PA, Glaser KJ, Manduca A, Romano AJ, Ehman RL: MR elastography as a method for the assessment of myocardial stiffness: Comparison with an established pressure-volume model in a left ventricular model of the heart. Magn Reson Med. 2009, 62: 135-140. 10.1002/mrm.21991.PubMed CentralView ArticlePubMedGoogle Scholar
- Elgeti T, Rump J, Hamhaber U, Papazoglou S, Hamm B, Braun J, Sack I: Cardiac Magnetic Resonance Elastography - Initial Results. Invest Radiol. 2008, 43: 762-772. 10.1097/RLI.0b013e3181822085.View ArticlePubMedGoogle Scholar
- Elgeti T, Laule M, Kaufels N, Schnorr J, Hamm B, Samani A, Braun J, Sack I: Cardiac MR elastography: comparison with left ventricular pressure measurement. J Cardiovasc Magn Reson. 2009, 11-44.Google Scholar
- Slama M, Maizel J: Echocardiographic measurement of ventricular function. Curr Opin Crit Care. 2006, 12: 241-248. 10.1097/01.ccx.0000224869.86205.1a.View ArticlePubMedGoogle Scholar
- Asbach P, Klatt D, Hamhaber U, Braun J, Somasundaram R, Hamm B, Sack I: Assessment of liver viscoelasticity using multifrequency MR elastography. Magn Reson Med. 2008, 60: 373-379. 10.1002/mrm.21636.View ArticlePubMedGoogle Scholar
- Woodrum DA, Romano AJ, Lerman A, Pandya UH, Brosh D, Rossman PJ, Lerman LO, Ehman RL: Vascular wall elasticity measurement by magnetic resonance imaging. Magn Reson Med. 2006, 56: 593-600. 10.1002/mrm.20991.View ArticlePubMedGoogle Scholar
- Weissler AM, Harris WS, Schoenfeld CD: Systolic time intervals in heart failure in man. Circulation. 1968, 37: 149-159.View ArticlePubMedGoogle Scholar
- Yin M, Talwalkar JA, Glaser KJ, Manduca A, Grimm RC, Rossman PJ, Fidler JL, Ehman RL: Assessment of hepatic fibrosis with magnetic resonance elastography. Clin Gastroenterol Hepatol. 2007, 5: 1207-1213. 10.1016/j.cgh.2007.06.012. e1202PubMed CentralView ArticlePubMedGoogle Scholar
- Sinkus R, Siegmann K, Xydeas T, Tanter M, Claussen C, Fink M: MR elastography of breast lesions: understanding the solid/liquid duality can improve the specificity of contrast-enhanced MR mammography. Magn Reson Med. 2007, 58: 1135-1144. 10.1002/mrm.21404.View ArticlePubMedGoogle Scholar
- Wuerfel J, Paul F, Beierbach B, Hamhaber U, Klatt D, Papazoglou S, Zipp F, Martus P, Braun J, Sack I: MR-elastography reveals degradation of tissue integrity in multiple sclerosis. Neuroimage. 2010, 49: 2520-2525. 10.1016/j.neuroimage.2009.06.018.View ArticlePubMedGoogle Scholar
- Zile MR, Brutsaert DL: New concepts in diastolic dysfunction and diastolic heart failure: Part I: diagnosis, prognosis, and measurements of diastolic function. Circulation. 2002, 105: 1387-1393. 10.1161/hc1102.105289.View ArticlePubMedGoogle Scholar
- van Heerebeek L, Borbély A, Niessen HWM, Bronzwaer JGF, van der Velden J, Stienen GJM, Linke WA, Laarman GJ, Paulus WJ: Myocardial structure and function differ in systolic and diastolic heart failure. Circulation. 2006, 113: 1966-1973. 10.1161/CIRCULATIONAHA.105.587519.View ArticlePubMedGoogle Scholar
- Fomovsky G, Thomopoulos S, Holmes J: Contribution of extracellular matrix to the mechanical properties of the heart. J Mol Cell Cardiol. 2009Google Scholar
- Oh JK, Tajik J: The return of cardiac time intervals: the phoenix is rising. J Am Coll Cardiol. 2003, 42: 1471-1474. 10.1016/S0735-1097(03)01036-2.View ArticlePubMedGoogle Scholar
- Lauboeck H: Echocardiographic study of the isovolumetric contraction time. J Biomed Eng. 1980, 2: 281-284. 10.1016/0141-5425(80)90121-1.View ArticlePubMedGoogle Scholar
- Spencer KT, Kirkpatrick JN, Mor-Avi V, Decara JM, Lang RM: Age dependency of the Tei index of myocardial performance. J Am Soc Echocardiogr. 2004, 17: 350-352. 10.1016/j.echo.2004.01.003.View ArticlePubMedGoogle Scholar
- Garcia MJ, Thomas JD, Klein AL: New Doppler echocardiographic applications for the study of diastolic function. J Am Coll Cardiol. 1998, 32: 865-875. 10.1016/S0735-1097(98)00345-3.View ArticlePubMedGoogle Scholar
- Fung YC: Heart Muscle. In Biomechanics - Mechanical Properties of Living Tissue. 1993, New York: Springer-Verlag, 427-465.Google Scholar
- Nagueh SF, Appleton CP, Gillebert TC, Marino PN, Oh JK, Smiseth OA, Waggoner AD, Flachskampf FA, Pellikka PA, Evangelista A: Recommendations for the evaluation of left ventricular diastolic function by echocardiography. J Am Soc Echocardiogr. 2009, 22: 107-133. 10.1016/j.echo.2008.11.023.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.