Cardiovascular magnetic resonance myocardial feature tracking detects quantitative wall motion during dobutamine stress
Journal of Cardiovascular Magnetic Resonance volume 13, Article number: 58 (2011)
Dobutamine stress cardiovascular magnetic resonance (DS-CMR) is an established tool to assess hibernating myocardium and ischemia. Analysis is typically based on visual assessment with considerable operator dependency. CMR myocardial feature tracking (CMR-FT) is a recently introduced technique for tissue voxel motion tracking on standard steady-state free precession (SSFP) images to derive circumferential and radial myocardial mechanics.
We sought to determine the feasibility and reproducibility of CMR-FT for quantitative wall motion assessment during intermediate dose DS-CMR.
10 healthy subjects were studied at 1.5 Tesla. Myocardial strain parameters were derived from SSFP cine images using dedicated CMR-FT software (Diogenes MRI prototype; Tomtec; Germany). Right ventricular (RV) and left ventricular (LV) longitudinal strain (EllRV and EllLV) and LV long-axis radial strain (ErrLAX) were derived from a 4-chamber view at rest. LV short-axis circumferential strain (EccSAX) and ErrSAX; LV ejection fraction (EF) and volumes were analyzed at rest and during dobutamine stress (10 and 20 μg · kg-1· min-1).
In all volunteers strain parameters could be derived from the SSFP images at rest and stress. EccSAX values showed significantly increased contraction with DSMR (rest: -24.1 ± 6.7; 10 μg: -32.7 ± 11.4; 20 μg: -39.2 ± 15.2; p < 0.05). ErrSAX increased significantly with dobutamine (rest: 19.6 ± 14.6; 10 μg: 31.8 ± 20.9; 20 μg: 42.4 ± 25.5; p < 0.05). In parallel with these changes; EF increased significantly with dobutamine (rest: 56.9 ± 4.4%; 10 μg: 70.7 ± 8.1; 20 μg: 76.8 ± 4.6; p < 0.05). Observer variability was best for LV circumferential strain (EccSAX ) and worst for RV longitudinal strain (EllRV) as determined by 95% confidence intervals of the difference.
CMR-FT reliably detects quantitative wall motion and strain derived from SSFP cine imaging that corresponds to inotropic stimulation. The current implementation may need improvement to reduce observer-induced variance. Within a given CMR lab; this novel technique holds promise of easy and fast quantification of wall mechanics and strain.
Cardiovascular magnetic resonance (CMR) plays an increasingly important role in the diagnosis and assessment of coronary artery disease (CAD). It has evolved into a comprehensive clinical tool with the unique capability of assessing myocardial function; viability and perfusion in a single examination .
Wall motion analysis with CMR has a pivotal role in clinical practice. It is considered the gold standard for visualizing left ventricular (LV) endocardial wall motion at rest; as well as during low and high dose dobutamine stress to assess myocardial hibernation and ischemia. At the present time; image analysis is most commonly performed qualitatively. However diagnostic accuracy of qualitative assessment has been shown to be considerably operator dependant .
Deformation assessment of tagged lines within the myocardium may overcome these limitations however requires acquisition of additional tagging sequences and post processing . Recently CMR myocardial feature tracking (FT); a technique analogous to echocardiographic speckle tracking; has been introduced . CMR-FT allows tracking of tissue voxel motion of cine-CMR images with a potential to assess longitudinal; circumferential and radial myocardial strain as well as velocity; displacement and torsion independent of additional sequences. A good agreement of CMR-FT versus myocardial tagging with harmonic phase imaging (HARP) as a reference standard has been demonstrated . However it is unclear; whether this approach would allow the response to dobutamine stress to be quantified . The aim of the current study was to determine the ability of CMR-FT for quantitative wall motion assessment at rest and during intermediate dose dobutamine stress in healthy volunteers.
Ten healthy volunteers underwent CMR on a 1.5 Tesla scanner (Intera R 18.104.22.168; Philips Medical Systems; Best; The Netherlands). The study protocol was approved by the Institutional Review Board at the University of Nebraska Medical Center. All participants gave written informed consent.
Cardiovascular magnetic resonance
All CMR measurements were performed in the supine position using a 5-channel cardiac surface coil. LV dimensions and function were assessed with an ECG-gated steady state free-precession cine sequence during brief periods of breath-holding in the following planes: ventricular 2-chamber; 4-chamber; and 12 to 14 equidistant short-axis planes (slice thickness 6-8 mm; gap 0-2 mm) completely covering both ventricles. The field of view was 360 × 480 mm and matrix size 196 × 172. Dobutamine stress imaging was performed as previously described . Repeat short-axis stacks were acquired with 10 and 20 μg · kg-1· min-1 of dobutamine; respectively.
Ventricular volumes and function
End-diastolic (EDV) and end-systolic volumes (ESV); stroke volume (SV); and ejection fraction (EF) were measured as previously described using commercially available software packages (View Forum; Philips) . Ventricular volumes were adjusted to body surface area. All parameters were analysed at rest; 10 and 20 μg · kg-1· min-1 of dobutamine stress.
CMR-FT analysis of strain was performed using a dedicated software prototype (Diogenes MRI; Tomtec; Germany). The 4-chamber view was used to calculate right ventricular (RV) and LV longitudinal strain and LV radial strain (EllRV and EllLV and ErrLAX) at rest. LV short axis circumferential (EccSAX) and radial strains (ErrSAX) were derived from a mid-ventricular short-axis view containing both papillary muscles. The RV upper septal insertion point of the LV was manually detected to allow accurate segmentation according to a recognized standard model . Endocardial contours were manually drawn in all analyzed slices by one skilled observer (AS; 7 years of experience). EccSAX and ErrSAX were analysed at rest; 10 and 20 μg · kg-1· min-1 of dobutamine stress. A second observer (SK; 4 years of experience) re-analysed the images to assess inter-observer variability. The mid-ventricular short axis images analyzed by the second observer were at exactly the same slice position as for the first observer. The first observer repeated the measurements after a period of 4 weeks to assess intra-observer variability. Figure 1 shows a representative example of the tracking of LV and RV in the respective views.
Comparison with natural radial strain
Natural radial strain values were obtained as an external reference standard and compared to the respective CMR-FT ErrSAX values . In brief end-diastolic and end-systolic wall-thickness (EDWT and ESWT) were quantified in identical segments as analysed for ErrSAX using commercially available software (Philips View Forum; The Netherlands) . Natural radial strain values were calculated according to the following equation: log e (ESWT/EDWT) as previously validated . 95% confidence intervalls of the difference and p-values were calculated to compare the 2 techniques .
We have applied a paired t-test followed by a Bonferroni-Holm correction as a multiple test procedure to compare measurements at rest and with dobutamine stress after proving a normal distribution of the sample. Intra- and inter-observer variability analysis were performed using the method proposed by Bland and Altman . A p-value of <0.05 was considered statistically significant. All data analysis was performed with PASW statistics for Mac 18.0.0 (SPSS Inc.; Chicago; Illinois; USA).
The image quality was sufficient to perform strain analysis in all segments for all subjects. Gender was equally distributed and LV and RV volumes were within normal limits . Participant demographics are shown in table 1. There were no side effects to dobutamine exposure. There was significant (p < 0.05) increase of heart rate; mean blood pressure and cardiac output between rest and both levels of dobutamine as well as between 10 and 20 μg · kg-1· min-1 of dobutamine
Strain parameters at rest
Results at rest are displayed in table 1.
Dobutamine stress cardiovascular magnetic resonance
Changes in EccSAX and ErrSAX were significant between rest and both levels of dobutamine as well as between 10 and 20 μg · kg-1· min-1 of dobutamine (table 2; table 3; figure 2). In parallel with these changes LV-EF increased significantly with 10 and 20 μg · kg-1· min-1 of dobutamine (table 2; figure 2).
Intra- and inter-observer variability
All parameters were reproducible on an intra- and inter-observer level. Bland Altman plots are displayed in figure 3 and table 4 shows the 95% confidence intervals of the difference between the repeated measurements. Observer variability at rest was best for EccSAX and worst for EllRV as determined by 95% confidence intervals of the difference. Observer variability did not significantly increase with dobutamine stress (table 5).
Comparison with natural radial strain
There was reasonable agreement between mean ErrSAX and natural radial strain as demonstrated in figure 4.
The current study includes a unique population of healthy volunteers studied at rest and with DS-CMR and demonstrates several important findings.
Firstly; CMR-FT can quantify wall motion changes between rest and dobutamine stress. Secondly; we noted considerable intra- and inter-observer variability for all parameters; which was most pronounced for RV longitudinal strain and smallest for LV circumferential strain. Thirdly; normal values of CMR-FT cover a large range with considerable overlap between rest and stress between different subjects (Figure 2) and closely correlate with hemodynamic changes secondary to changed LV function.
In essence; our data demonstrate that CMR-FT strain parameters can be derived from routine SFFP cine sequences at varying levels of DS-CMR. The current reference standard for quantitative wall motion assessment with CMR is myocardial tagging [11, 14, 15]. Myocardial tagging based strain assessment has been demonstrated to improve diagnostic accuracy of DS-CMR in patients with suspected CAD as well as in patients with myocardial hibernation . Adding quantitative analysis to DS-CMR in CAD may not only increase diagnostic accuracy as compared to visual analysis for the detection of ischemia during high-dose dobutamine stress  but may also detect quantitative changes in myocardial strain already detectable at lower stress levels; thereby increasing feasibility of the test .
There is evidence that CMR-FT may also be of clinical utility. Hor et al have recently shown that CMR-FT provides similar results compared to HARP myocardial tagging in a patient population with Duchenne muscular dystrophy . Maret et al demonstrated that CMR-FT can be used in CAD patients to accurately detect strain in the radial and longitudinal direction correlating with the presence of myocardial scarring . CMR-FT is also useful for the assessment of myocardial viability . Detection of quantitative contractile reserve in patients with myocardial hibernation using low dose dobutamine has the potential to predict functional recovery after revascularisation with higher accuracy in the future . In this context; CMR-FT may serve as an additional tool alongside conventional visual analysis; thus facilitating the detection of subtle wall motion abnormalities and the identification of contractile reserve; particularly for the less experienced observer.
The importance of the intra- and interobserver variability documented in the current study needs to be taken into consideration. This has also been reported by echocardiography based speckle tracking studies and our results are similar . The parameter with highest variability at rest in the current study was longitudinal strain of the RV indicating that the analysis of the thin-walled RV with CMR-FT is not yet adequately accurate. This might also be explained by difficulties in endocardial tracking due to difficulty in accurately following the tricuspid valve annulus motion with the CMR-FT software; and RV trabeculations that also lead to greater variability in RV volumetric assessment . The most robust parameter in our study was circumferential strain of the LV; which might be clinically valuable. CMR- FT algorithm allows reliable and easy border tracking; the frame-to-frame displacement of features tracked is equivalent to evaluating the local velocity (ratio between displacement and time interval); allowing automatic evaluation of tissue motion during the cardiac cycle. The tracking results from this algorithm may be more reliable due to the inherently high image quality with CMR. However echocardiographic speckle tracking has better temporal resolution than CMR-FT. In addition our collective consisted of good breath holders and as a consequence we obtained good image quality at an intermediate stress level. It is therefore not surprising that observer-induced variance did not significantly increase with stress. Whether potentially degraded image quality at higher stress levels in patients who are not able to hold their breath would substantially obscure CMR-FT results needs to be prospectively assessed. Interestingly results in radial strain from matching segments from the short-axis and long-axis orientation were not equal. Whether this could be explained by more extensive through-plane motion in the short-axis orientation or increased susceptibility of the 4-chamber view to the breath-holding position of the diaphragm needs to be investigated in healthy volunteers and in patients with scarred areas and segments with wall motion abnormalities. In particular future studies need to investigate whether DS-CMR accuracy could be improved with CMR-FT information that is available with any DS-CMR stress study. There is evidence to suggest that these quantitative parameters have prognostic implications. Stanton and colleagues demonstrated that automated echocardiography speckle-tracking derived global EllLV is a superior predictor of outcome compared to either EF or wall motion score index; and suggested that EllLV may even become the optimal method to assess global left ventricular systolic function . As CMR-FT is a relatively new method such evidence is not yet available and future studies need to investigate whether CMR-FT could also provide such assessment.
The sample size of the current study was relatively small. Future studies will need to reassess these parameters in a larger cohort of volunteers and patients. Global Ell; which has been previously shown to be an important; prognostic echocardiographic parameter was only assessed at rest. This was due to time constraints as a whole stack of short axis images had been acquired at each stage of dobutamine for volumetry. Also we did not perform any echocardiographic measurements to compare with CMR-FT data; which needs to be addressed in future studies. Finally the current work aimed to determine the feasibility of CMR-FT during dobutamine stress in a collective of healthy volunteers. Future research needs to prospectively validate this novel technique in pathologies such as coronary artery disease; valvular disease or congenital disorders.
CMR-FT allows derivation of strain mechanics from SSFP cine images at rest and during dobutamine stress CMR corresponding to global hemodynamic changes. The current analysis algorithm requires improvement to reduce observer-induced variance; which at present is comparable to data reported from 2D strain by echocardiography. Within a given CMR lab; this novel CMR-FT technique holds promise for easy and fast quantification of wall mechanics and strain.
acute myocardial infarction
cardiovascular magnetic resonance
dobutamine stress cardiovascular magnetic resonance
left ventricular short-axis circumferential strain
left ventricular longitudinal strain
right ventricular longitudinal strain
end-diastolic wall thickness
left ventricular long-axis radial strain
left ventricular short-axis radial strain
myocardial feature tracking
Morton G, Schuster A, Perera D, Nagel E: Cardiac magnetic resonance imaging to guide complex revascularization in stable coronary artery disease. European heart journal. 2010, 31: 2209-2215. 10.1093/eurheartj/ehq256.
Paetsch I, Jahnke C, Ferrari VA, Rademakers FE, Pellikka PA, Hundley WG, Poldermans D, Bax JJ, Wegscheider K, Fleck E, Nagel E: Determination of interobserver variability for identifying inducible left ventricular wall motion abnormalities during dobutamine stress magnetic resonance imaging. Eur Heart J. 2006, 27: 1459-1464.
Zerhouni EA, Parish DM, Rogers WJ, Yang A, Shapiro EP: Human heart: tagging with MR imaging--a method for noninvasive assessment of myocardial motion. Radiology. 1988, 169: 59-63.
Maret E, Todt T, Brudin L, Nylander E, Swahn E, Ohlsson JL, Engvall JE: Functional measurements based on feature tracking of cine magnetic resonance images identify left ventricular segments with myocardial scar. Cardiovasc Ultrasound. 2009, 7: 53-10.1186/1476-7120-7-53.
Hor KN, Gottliebson WM, Carson C, Wash E, Cnota J, Fleck R, Wansapura J, Klimeczek P, Al-Khalidi HR, Chung ES, et al: Comparison of magnetic resonance feature tracking for strain calculation with harmonic phase imaging analysis. JACC Cardiovasc Imaging. 2010, 3: 144-151. 10.1016/j.jcmg.2009.11.006.
Schuster A, Nagel E: Toward Full Quantification of Wall Motion with MRI. curr cardiovasc imaging rep. 2011, 85-86. 10.1007/s12410-011-9072-x.
Nagel E, Lehmkuhl HB, Bocksch W, Klein C, Vogel U, Frantz E, Ellmer A, Dreysse S, Fleck E: Noninvasive diagnosis of ischemia-induced wall motion abnormalities with the use of high-dose dobutamine stress MRI: comparison with dobutamine stress echocardiography. Circulation. 1999, 99: 763-770.
Thiele H, Paetsch I, Schnackenburg B, Bornstedt A, Grebe O, Wellnhofer E, Schuler G, Fleck E, Nagel E: Improved accuracy of quantitative assessment of left ventricular volume and ejection fraction by geometric models with steady-state free precession. Journal of cardiovascular magnetic resonance: official journal of the Society for Cardiovascular Magnetic Resonance. 2002, 4: 327-339. 10.1081/JCMR-120013298.
Cerqueira MD, Weissman NJ, Dilsizian V, Jacobs AK, Kaul S, Laskey WK, Pennell DJ, Rumberger JA, Ryan T, Verani MS, Imaging AHAWGoMSaRfC: Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart: a statement for healthcare professionals from the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association. Circulation. 2002, 105: 539-542. 10.1161/hc0402.102975.
Mirsky I, Parmley WW: Assessment of passive elastic stiffness for isolated heart muscle and the intact heart. Circulation Research. 1973, 33: 233-243.
Attili AK, Schuster A, Nagel E, Reiber JHC, van der Geest RJ: Quantification in cardiac MRI: advances in image acquisition and processing. The international journal of cardiovascular imaging. 2010, 26 (Suppl 1): 27-40.
Bland JM, Altman DG: Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986, 1: 307-310.
Hudsmith LE, Petersen SE, Francis JM, Robson MD, Neubauer S: Normal human left and right ventricular and left atrial dimensions using steady state free precession magnetic resonance imaging. Journal of cardiovascular magnetic resonance: official journal of the Society for Cardiovascular Magnetic Resonance. 2005, 7: 775-782. 10.1080/10976640500295516.
Shehata ML, Cheng S, Osman NF, Bluemke DA, Lima JAC: Myocardial tissue tagging with cardiovascular magnetic resonance. Journal of cardiovascular magnetic resonance: official journal of the Society for Cardiovascular Magnetic Resonance. 2009, 11: 55-10.1186/1532-429X-11-55.
Ibrahim E-SH: Myocardial tagging by Cardiovascular Magnetic Resonance: evolution of techniques--pulse sequences, analysis algorithms, and applications. Journal of cardiovascular magnetic resonance: official journal of the Society for Cardiovascular Magnetic Resonance. 2011, 13: 36-10.1186/1532-429X-13-36.
Nagel E, Schuster A: Shortening without contraction: new insights into hibernating myocardium. JACC Cardiovasc Imaging. 2010, 3: 731-733. 10.1016/j.jcmg.2010.05.002.
Kuijpers D, Ho KY, van Dijkman PR, Vliegenthart R, Oudkerk M: Dobutamine cardiovascular magnetic resonance for the detection of myocardial ischemia with the use of myocardial tagging. Circulation. 2003, 107: 1592-1597. 10.1161/01.CIR.0000060544.41744.7C.
Korosoglou G, Lehrke S, Wochele A, Hoerig B, Lossnitzer D, Steen H, Giannitsis E, Osman NF, Katus HA: Strain-encoded CMR for the detection of inducible ischemia during intermediate stress. JACC Cardiovasc Imaging. 2010, 3: 361-371. 10.1016/j.jcmg.2009.11.015.
Schuster A, Paul M, Bettencourt N, Morton G, Chiribiri A, Ishida M, Hussain S, Jogiya R, Kutty S, Bigalke B, et al: Cardiovascular magnetic resonance myocardial feature tracking for quantitative viability assessment in ischemic cardiomyopathy. International Journal of Cardiology. 2011, [E-pub ahead of print], doi:10.1016/j.ijcard.2011.10.137
Bree D, Wollmuth JR, Cupps BP, Krock MD, Howells A, Rogers J, Moazami N, Pasque MK: Low-dose dobutamine tissue-tagged magnetic resonance imaging with 3-dimensional strain analysis allows assessment of myocardial viability in patients with ischemic cardiomyopathy. Circulation. 2006, 114: I33-36. 10.1161/CIRCULATIONAHA.105.600858.
Bansal M, Jeffriess L, Leano R, Mundy J, Marwick TH: Assessment of myocardial viability at dobutamine echocardiography by deformation analysis using tissue velocity and speckle-tracking. JACC Cardiovasc Imaging. 2010, 3: 121-131. 10.1016/j.jcmg.2009.09.025.
Beerbaum P, Barth P, Kropf S, Sarikouch S, Kelter-Kloepping A, Franke D, Gutberlet M, Kuehne T: Cardiac function by MRI in congenital heart disease: impact of consensus training on interinstitutional variance. J Magn Reson Imaging. 2009, 30: 956-966. 10.1002/jmri.21948.
Stanton T, Leano R, Marwick TH: Prediction of all-cause mortality from global longitudinal speckle strain: comparison with ejection fraction and wall motion scoring. Circ Cardiovasc Imaging. 2009, 2: 356-364.
AS receives grant support from the British Heart Foundation (BHF) (RE/08/003 and FS/10/029/28253) and the Biomedical Research Centre (BRC-CTF 196). SK receives grant support from the American College of Cardiology Foundation; the Edna Ittner Pediatric Foundation; and the Children's Hospital and Medical Center Foundation.
The authors declare that they have no competing interests.
AS and SK designed the study protocol; analyzed the data and drafted the manuscript. AP; PG and DAD performed the CMR studies and helped to draft the manuscript. VP and MRM participated in the study design and helped to draft the manuscript. BB participated in the study design and helped with the statistical analysis. PB and EN designed the study protocol and drafted the manuscript. All authors read and approved the final manuscript.
Andreas Schuster, Shelby Kutty, Philipp Beerbaum and Eike Nagel contributed equally to this work.
About this article
Cite this article
Schuster, A., Kutty, S., Padiyath, A. et al. Cardiovascular magnetic resonance myocardial feature tracking detects quantitative wall motion during dobutamine stress. J Cardiovasc Magn Reson 13, 58 (2011). https://doi.org/10.1186/1532-429X-13-58