Author + information
- Received September 12, 2011
- Revision received November 16, 2011
- Accepted December 13, 2011
- Published online August 1, 2012.
- Loren P. Budge, MD⁎,
- Adam S. Helms, MD⁎,
- Michael Salerno, MD, PhD⁎,
- Christopher M. Kramer, MD⁎,
- Frederick H. Epstein, PhD† and
- Kenneth C. Bilchick, MD, MS⁎,⁎ ()
- ↵⁎Reprint requests and correspondence:
Dr. Kenneth Bilchick, University of Virginia Health System, Cardiovascular Division, Post Office Box 800158, Charlottesville, Virginia 22908
Objectives We sought to assess the effectiveness of automated mechanical dyssynchrony (MD) parameters based on regional heterogeneity of strain (circumferential [CURE], longitudinal [LURE], and radial uniformity ratio estimates) relative to parameters based on regional time to peak contraction with cardiac magnetic resonance (CMR) cine DENSE (Displacement Encoding with Stimulated Echoes) validated with myocardial tissue tagging (MTT) strain data.
Background Dyssynchrony measures based on the Fourier transformation (FT) of regional strain, such as CURE (previously evaluated in cardiac resynchronization therapy candidates), directly assess MD and yield straightforward global dyssynchrony indexes; however, performance relative to the 12-segment standard deviation of time to peak strain (SD12) or maximal regional delay in time to peak strain is unknown.
Methods Cine DENSE and MTT were obtained with CMR (1.5-T Siemens Avanto, Siemens, Erlangen, Germany) in 13 canines: 3 normal control subjects, 5 with tachycardia pacing-induced heart failure (HF) and left bundle branch ablation (LBBB-HF), and 5 with HF and narrow QRS (NQRS-HF). Strain and dyssynchrony parameters were determined with both CMR methods.
Results Both HF groups had reduced peak strains and left ventricular ejection fraction compared with normal cases. There was strong agreement between cine DENSE and MTT on the basis of intraclass correlation coefficients (CURE: 0.99, 95% CI: 0.96 to 1.00; LURE: 0.92, 95% CI: 0.77 to 0.98; circumferential strain [ECC]: 0.95, 95% CI: 0.72 to 0.99; longitudinal strain [ELL]: 0.82, 95% CI: 0.42 to 0.97). The FT-based metrics (scale 0 to 1), in particular CURE, discriminated highly between LBBB-HF and NQRS-HF groups (median difference): CURE: 0.60, 95% CI: 0.43 to 0.76; LURE: 0.39, 95% CI: 0.19 to 0.58; radial uniformity ratio estimate: 0.22, 95% CI: 0.04 to 0.40). In contrast, relative confidence intervals for group differences in time-to-peak parameters were wide, indicating less consistent discrimination (median difference): SD12-ECC: 52.5, 95% CI: −4.0 to 109.2; SD12-ELL: 40.9, 95% CI: −5.3 to 87.1; SD12-radial strain: 42.0, 95% CI: 0.4 to 83.6). Correlations between FT-based and time-to-peak parameters were significant (CURE/SD12-ECC: r = −0.62, p = 0.03; LURE/SD12-ELL: r = −0.76, p = 0.005) but not as tight as correlations between time-to-peak parameters.
Conclusions Automated FT-based circumferential, radial, and longitudinal dyssynchrony measures compare favorably with time-to-peak parameters. Cine DENSE was effective for this application and validated with MTT. Further clinical evaluation in cardiac resynchronization therapy candidates with CMR or other imaging modalities is warranted.
The overall benefit of cardiac resynchronization therapy (CRT) to selected patients with heart failure (HF) has been shown in clinical trials (1–4), but the nonresponse rate to CRT with current clinical indications is still significant at approximately 30% to 40% (5). Although other biological factors are involved in CRT response (6), cardiac imaging remains the best-studied tool to achieve more accurate identification of CRT responders. Even so, a number of early mechanical dyssynchrony parameters were unable to identify CRT responders effectively in a multicenter clinical trial (5). Although these results have led some to question whether cardiac imaging can identify patients most likely to respond to CRT, an alternative response has been to look for more effective imaging techniques and approaches for CRT patients. For example, there have been promising results shown for echo speckle tracking (7), 3-dimensional (3D) echo (8), and cardiac magnetic resonance (CMR) with myocardial tissue tagging (MTT) (9). Magnetic resonance (MR) cine Displacement Encoding with Stimulated Echoes (DENSE) is another imaging approach with high potential. DENSE encodes tissue displacement directly into the phase of the MR signal (10), provides high spatial resolution, accurately assesses circumferential, longitudinal, and radial strain, provides rapid strain analysis (11), and offers automatic contouring (12) without the need for tag detection or a minimum tag density.
In addition to the imaging modality used, the parameterization of dyssynchrony is also an important factor. Widely used dyssynchrony parameters are based on time to peak contraction, indexed as the delay in opposing walls or standard deviation of time to peak strain in 12 segments (SD12) (7,13–15). Although associated studies have shown positive results, identification in peaks in tissue velocity or strain curves might be time-consuming, subjective, and very difficult (or impossible) for segments with infarction, akinesis, dyskinesis, or very low maximum contractile strain. An alternative approach (not specific to CMR) is the use of automated global circumferential (circumferential uniformity ratio estimate [CURE]) (9,16), radial (radial uniformity ratio estimate [RURE]), and longitudinal (longitudinal uniformity ratio estimate [LURE]) dyssynchrony indexes, which are based on the Fourier transformation (FT) of regional strain (FT-based parameters) and do not require assessment of regional strain or velocity peaks.
In the present study, we aimed to assess the hypotheses that FT-based parameters would perform better than corresponding timing-based parameters and that cine DENSE would be feasible and valid for this application on the basis of comparison with MTT. We elected to compare these indexes in a validated animal model of HF with and without left bundle branch block (LBBB) rather than in a human CRT study. Because HF with dyssynchrony modeled this way has been shown to improve with resynchronization (16), it was an ideal model to test our hypotheses. Although the study was implemented with CMR, the findings of this study may be extended to any imaging technique that can generate regional strain curves from standard short-axis and long-axis slices.
Approval for the study was obtained through the University of Virginia Animal Care and Use Committee. As a validated model of dyssynchronous HF (17), we performed left bundle branch ablation (50 W, temperature ≤60°) before tachycardia pacing (LBBB-HF) with a 4-mm-tip ablation catheter in 5 canines, resulting in widening of the QRS duration from 55 to 60 ms to approximately 120 ms with a typical LBBB appearance. Persistence of LBBB was confirmed after a 30-min waiting period as well as at the time of final MR examination. Cardiomyopathy was induced in a total of 10 canines with right atrial tachy-pacing at 180 beats/min for 5 weeks (the 5 LBBB-HF animals, and 5 animals without left bundle ablation [NQRS-HF]). Three normal canines did not undergo ablation or tachy-pacing.
With a 1.5T Avanto scanner (Siemens Medical Solutions, Erlangen, Germany) with a 4-channel phased-array chest radiofrequency coil, we performed cine DENSE imaging in standard short-axis and long-axis planes in two orthogonal directions for each plane with the following parameters (10,18): interleaved spiral readout with 6 interleaves/image; repetition time/echo time 17/1.9 ms; section thickness 8 mm; field of view 350 × 350 mm; flip angle 15°; pixel size 2.8 × 2.8 mm; fat suppression; echo spacing 17 ms with view sharing; and displacement-encoding frequency 0.1 cycles/mm. The MTT images were obtained with a gradient echo grid tagging sequence: tag spacing 7 mm; repetition time/echo time 8.7/4.2 ms; slice thickness 8 mm; field of view 250 × 250 mm; flip angle 14°; pixel size 1.4 × 1.0 mm; and temporal resolution 17.5 ms with view sharing (19).
Image analysis and calculation of dyssynchrony metrics
After image acquisition, segmentation of the left ventricular (LV) myocardium was performed, a phase-unwrapping algorithm was applied to LV myocardium pixels, and displacements were calculated (11). Lagrangian strain was computed from displacements in 24 short-axis segments in multiple slices and then projected in both the radial (ERR) and circumferential (ECC) directions relative to the LV center of mass. Longitudinal strain (ELL) was determined by analysis of strain in standard long-axis planes in basal and mid-cavity segments. Left ventricular volumes and ejection fraction were calculated from cine steady-state free precession images with Argus software. Strains were calculated from cine DENSE images with an algorithm developed in Matlab (MathWorks, Natick, Massachusetts) and from MTT images with harmonic phase analysis software (Diagnosoft, Inc., Cary, North Carolina) (11).
Dyssynchrony was assessed by 2 types of analyses for each strain component: variability in time to peak strain and the FT-based uniformity ratio estimates (CURE, LURE, and RURE). Basic timing delays for each strain component were calculated as the difference in time from QRS onset to peak strain for opposing LV segments with earliest and latest peak strain (ECC-Delay, ELL-Delay, and ERR-Delay). The SD12-ECC, SD12-ELL, and SD12-ERR were constructed on the basis of the SD of time to peak strain in 6 basal and 6 mid-cavity segments (12 total segments) similar to the method Yu et al. (14) used for regional tissue velocity.
The CURE, LURE, and RURE were calculated as previously described (16). Briefly, this technique involves a simple and automated analysis on the basis of the FT of the spatial distribution of strain. Because synchronous contraction is represented primarily by zero-order FT terms and dyssynchronous contraction by first-order FT terms, a ratio measure involving these terms (the square root of the ratio of the zero-order term to the sum of zero and first-order terms) can be used to index dyssynchrony on a scale between 0 (dyssynchrony) and 1 (synchrony) (9).
The Wilcoxon rank sum test with exact p values was used to compare continuous variables, whereas the Hodges Lehmann estimate and exact 95% confidence interval (CI) were used to evaluate the median differences between groups. Boxplots were also constructed for a visual comparison of these differences between groups. For pairwise correlations, the Pearson correlation coefficient (r) and associated significance level are reported. Comparison between cine DENSE and MTT was performed with the intraclass correlation coefficients for absolute agreement between methods. Results are presented graphically with both scatter plots with regression lines and Bland-Altman plots. Statistical analysis was performed with SAS (version 9.2, SAS, Cary, North Carolina).
Successful left bundle branch ablations were completed in all LBBB-HF animals, with resulting QRS durations of 120 ms and characteristic LBBB morphology on the 12-lead electrocardiogram; QRS durations were 55 to 60 ms in the NQRS-HF group and control groups. After 5 weeks of tachy-pacing, all LBBB-HF and NQRS-HF group animals displayed LV dysfunction with left ventricular ejection fraction (LVEF) <0.30. The LBBB-HF group showed a trend for lower LVEF (median 0.15 vs. 0.24, p = 0.14) and had increased LV end-diastolic volumes (median 110 vs. 77 ml, p = 0.008) relative to the NQRS-HF group. Normal animals had a greater LVEF than HF animals (median 0.57 vs. 0.15, p = 0.017). Examples of CMR steady-state free precession and Lagrangian displacement movies for both an LBBB-HF case and NQRS-HF case may be found in Online Videos 1, 2, 3, and 4.
Comparison of cine DENSE and MTT
As shown in Figure 1, agreement between cine DENSE and MTT was excellent for peak ECC (R2 = 0.94, p < 0.01; intraclass correlation coefficient [ICC] for absolute agreement between methods: 0.95, 95% CI: 0.72 to 0.99) and CURE (R2 = 0.98, p < 0.001; ICC: 0.99, 95% CI: 0.96 to 1.00). There was also good agreement between DENSE and MTT for peak ELL (peak ELL: R2 = 0.86, p < 0.01; ICC: 0.82, 95% CI: 0.42 to 0.97) and LURE (R2 = 0.85, p < 0.001; ICC: 0.92, 95% CI: 0.77 to 0.98) (Online Appendix). As previously reported, MTT-based radial strain and dyssynchrony measures were of insufficient precision to include in our analysis (20).
Comparison between NQRS-HF and LBBB-HF groups
The NQRS-HF and LBBB-HF animals had marked heterogeneity in regional strain curves (see Online Appendix for examples of circumferential, longitudinal, and radial regional strain curves for each group). Although the regional strain polarity (positive or negative) was generally the same regardless of the strain component for NQRS-HF animals, typical regional strain versus time curves for the LBBB-HF group show simultaneous positive and negative strains in different sectors. This indicates early anteroseptal contraction and lateral wall stretch, followed by late lateral wall contraction and anteroseptal stretch. Statistically significant differences in all FT-based parameters (p < 0.01 for CURE and LURE; p = 0.016 for RURE) are evident from Table 1, although the magnitude of the difference between LBBB-HF and NQRS-HF groups was greatest for CURE, followed by LURE and RURE (CURE: 0.60, 95% CI: 0.43 to 0.76; LURE: 0.39, 95% CI: 0.19 to 0.58; RURE: 0.22, 95% CI: 0.04 to 0.40).
Boxplots comparing FT-based metrics with corresponding timing-based dyssynchrony parameters (SD12 and opposing wall delay) are shown in Figures 2, 3, and 4⇓⇓⇓ for circumferential, longitudinal, and radial parameters, respectively. From these boxplots, it is evident that although time-to-peak parameters indicate more dyssynchrony for LBBB-HF relative to NQRS-HF cases, the time-to-peak parameters are less consistent within each group (relative to FT-based parameters). The less consistent discrimination with time-to-peak parameters is also apparent from Table 1, which shows that even time-to-peak parameters based on 12 segments have much wider CIs (relative to the scale of the parameter) for the median differences between the LBBB-HF and NQRS-HF groups (SD12-ECC: 52.5, 95% CI: −4.0 to 109.2; SD12-ELL: 40.9, 95% CI: −5.3 to 87.1; SD12-ERR: 42.0, 95% CI: 0.4 to 83.6).
The 2 LBBB-HF cases in Figure 5 both have significant dyssynchrony and illustrate some of the advantages of FT-based parameters. In panel A, the circumferential and radial strain curves show absent early septal contraction (sometimes seen when LV systolic dysfunction is more severe or infarction is present). The CURE and RURE still both correctly indicate significant dyssynchrony, but time-to-peak measures do not because there are no early contractile strain peaks. There is better agreement in dyssynchrony assessments in panel B, which shows the classic pattern of early septal contraction and early free wall stretch followed by late free wall contraction and late septal stretch.
Comparison between normal and HF groups
Figures 2, 3, and 4 also show comparisons for all 9 dyssynchrony parameters between normal cases and those with HF. Dyssynchrony parameters in the normal control and NQRS-HF groups are similar, whereas the LBBB-HF group demonstrates the greatest amount of dyssynchrony. Median differences and CIs for dyssynchrony parameter differences between the LBBB-HF and normal groups are shown in Table 1.
Correlations among dyssynchrony metrics
Only moderate (statistically significant) correlations were present between the FT-based and time-to-peak parameters with ECC and ELL (CURE/SD12-ECC: r = −0.62, p = 0.03; LURE/SD12-ELL: r = −0.76, p = 0.005). As expected, the correlations were tighter between time-to-peak measures (SD12-ECC/ECC-Delay: r = 0.95, p < 0.001; SD12-ELL/ELL-Delay: r = 0.91, p < 0.001; SD12-ERR/ERR-Delay: r = 0.89, p < 0.001). A table with all the correlation data may be found in the Online Appendix.
The key findings of this study are: 1) cine DENSE provides high-quality dyssynchrony assessments on the basis of circumferential, longitudinal, or radial strain, with good correlation with MTT; 2) cine DENSE assessment of mechanical dyssynchrony in HF with FT-based parameters quantifying the extent of simultaneous myocardial stretch and contraction is feasible; and 3) FT-based dyssynchrony parameters (CURE, RURE, and LURE) derived from CMR effectively discriminate HF with and without LBBB and perform favorably relative to time-to-peak parameters.
Timing-based versus FT-based dyssynchrony parameters
Widely used dyssynchrony measures from echocardiography (7,13,14) and CMR published studies (15) are dependent on the time to peak contraction. Despite this, with both echocardiography and CMR, determination of time to peak contraction on the basis of regional tissue velocity or strain curves might be difficult when infarction, akinesis, or dyskinesis is present or when segmental velocity/strain curves have a predominant stretch pattern, low maximal amplitude, and/or multiple peaks. In contrast, FT-based measures such as CURE, LURE, and RURE measure the degree to which some segments contract while other segments stretch during the cardiac cycle and do not require manual identification of strain/velocity peaks, providing an automated and easily interpretable global measure of dyssynchrony. Because these FT-based dyssynchrony parameters can be easily determined with any imaging modality that generates strain curves from standard LV short-axis and long-axis views, this methodology has broad potential clinical applicability.
Clinical importance of dyssynchrony assessment by echocardiography
Early single-center echocardiographic studies with tissue Doppler imaging showed that a significant delay between time to peak systolic velocity at basal septal and basal lateral segments predicted a lower adverse clinical event rate after CRT (13). In addition, an alternative parameterization on the basis of the SD of time from QRS onset to peak systolic velocity for 12 LV segments (derived from tissue Doppler imaging) was also shown to be effective and associated with LV reverse remodeling after CRT (14). Subsequently, the multicenter PROSPECT (Predictors of Response to CRT) study “observed relatively low yield and high variability for TDI [tissue Doppler imaging] measures” (5). Although the methodology was criticized by some, the study seemed to have the positive effect of encouraging the development of more advanced echocardiographic dyssynchrony methods with assessments of radial and circumferential strain, such as echo speckle tracking and 3D echocardiography, and positive clinical studies were subsequently published on the basis of these methods (7,8). In addition, advanced strain analysis methods based on myocardial tagging (9) and cine DENSE (11,12) have also been developed for this purpose.
Clinical impact of refined dyssynchrony imaging
Despite a large number of CRT-related imaging studies, major guidelines for CRT patient selection in patients with a widened QRS do not include any criteria based on imaging (with the exception of the UK National Institute for Health and Clinical Excellence guidelines) (21). In spite of this, significant heterogeneity in CRT response still exists even in patients with a wide QRS, and there are compelling reasons to improve CRT selection criteria in these patients. These include an increased complication rate associated with CRT-D relative to the rate with standard ICDs. For example, in RAFT (Resynchronization/Defibrillation in Ambulatory Heart Failure Trial), the complication rate for CRT-D was 13%, which was significantly greater than the complication rate for standard ICDs (22). Procedure times and patient radiation exposure are also increased with CRT (23). In addition, CRT carries considerable costs related to the device, implantation, monitoring, and the need for multiple device changes over a lifetime. On the basis of these costs and current CRT nonresponse rates, if the medical community does not develop more effective ways to identify the best CRT candidates, it is possible that reimbursement could be restricted in the future, which could leave a number of potential CRT responders unable to receive this treatment.
In light of recent results for speckle tracking, 3D echocardiography, and CMR modalities such as DENSE and myocardial tagging, it is our opinion that cardiac imaging for mechanical dyssynchrony and viability holds the best promise for developing more meaningful CRT patient selection criteria, although a positive large multicenter randomized trial would be needed to influence guidelines. We believe that effective imaging for mechanical dyssynchrony requires both accurate assessment of cardiac mechanical function and optimal parameterization of the degree of dyssynchrony. Because the present study has addressed both these issues and demonstrated an advantage for FT-based parameters over time-to-peak parameters, we believe this work will have significant clinical impact. Of note, these FT-based dyssynchrony parameters could have potential broad application to strain data obtained not only from cine DENSE MR but also other imaging modalities capable of generating regional strain curves.
Limitations and rationale for translational approach
This is a translational rather than a clinical study, and our model of LBBB might not be representative of all types of mechanical dyssynchrony in human HF; however, this design provided a controlled and validated model to test our hypotheses. Further clinical investigation is indicated on the basis of these data. Of note, we chose a translational approach for several reasons. First, the model was necessary to demonstrate consistency of cine DENSE and MTT for dyssynchrony and strain assessment. Second, the pacing-induced HF model is a well-validated animal model of HF and produces LV dysfunction similar to the LV dysfunction present in patients referred for CRT (1,2). Third, induction of LBBB results in the prototypical model of mechanical dyssynchrony (17) that responds to resynchronization on the basis of prior published data (16). Fourth, the translational model has allowed us to determine that this technique has the potential to perform better than timing-based methods in human HF patients, which justifies further clinical study.
Automated FT-based circumferential, radial, and longitudinal dyssynchrony indexes based on heterogeneity of regional segmental strain can be easily determined without the need for time-to-peak contraction determinations and compare favorably with time-to-peak parameters. In the present study, we demonstrate the effectiveness of cine DENSE for this purpose, with validation based on strain obtained from harmonic phase analysis of MTT. Further clinical evaluation is warranted with cine DENSE CMR or other imaging modalities, because this methodology can be extended to any imaging modality capable of generating regional strain curves from standard LV short-axis and long-axis views.
The authors thank Xiaodong Zhong, PhD, Robert Janiczek, PhD, Andrew Gilliam, PhD, John Christopher, RT, and Craig Goodman, BS, for their assistance.
For supplementary figures and videos, please see the online version of this article.
Magnetic Resonance Cine Displacement Encoding With Stimulated Echoes Dyssynchrony Parameters for the Evaluation of Heart Failure
This work was supported by National Institutes of Health Grants K23 HL094761 (to Dr. Bilchick) and R01 EB001763 (to Dr. Epstein). Drs. Salerno, Kramer, and Epstein have received grant support from Siemens. Dr. Kramer is a consultant for St. Jude Medical. Pacemakers designated for animal use were donated by St. Jude Medical. All other authors have reported that they have no relationships relevant to the contents of this paper to disclose. Drs. Budge and Helms contributed equally to this work. Warren Manning, MD, served as Guest Editor for this paper.
- Abbreviations and Acronyms
- cardiac magnetic resonance
- cardiac resynchronization therapy (defibrillator)
- circumferential uniformity ratio estimate
- circumferential strain
- longitudinal strain
- radial strain
- Fourier transformation
- heart failure
- left bundle branch block
- longitudinal uniformity ratio estimate
- left ventricular
- left ventricular ejection fraction
- myocardial tissue tagging
- narrow QRS
- radial uniformity ratio estimate
- standard deviation of time to peak strain in 12 segments
- Received September 12, 2011.
- Revision received November 16, 2011.
- Accepted December 13, 2011.
- American College of Cardiology Foundation
- Chung E.S.,
- Leon A.R.,
- Tavazzi L.,
- et al.
- Aiba T.,
- Hesketh G.G.,
- Barth A.S.,
- et al.
- Tanaka H.,
- Nesser H.J.,
- Buck T.,
- et al.
- Bilchick K.C.,
- Dimaano V.,
- Wu K.C.,
- et al.
- Bax J.J.,
- Bleeker G.B.,
- Marwick T.H.,
- et al.
- Yu C.M.,
- Fung J.W.,
- Zhang Q.,
- et al.
- Westenberg J.J.,
- Lamb H.J.,
- van der Geest R.J.,
- et al.
- Leclercq C.,
- Faris O.,
- Tunin R.,
- et al.
- UK National Institute for Health and Clinical Excellence
- Perisinakis K.,
- Theocharopoulos N.,
- Damilakis J.,
- Manios E.,
- Vardas P.,
- Gourtsoyiannis N.