Feasibility and initial experience of assessment of mechanical dyssynchrony using cardiovascular magnetic resonance and semi-automatic border detection
Journal of Cardiovascular Magnetic Resonance volume 10, Article number: 49 (2008)
The systolic dyssynchrony index (SDI) has been introduced as a measure of mechanical dyssynchrony using three-dimensional echocardiography to select patients who may benefit from cardiac resynchronization therapy (CRT). However, three-dimensional echocardiography may be inadequate in a number of patients with suboptimal acoustic window and no single echocardiographic measure of dyssynchrony has proven to be of value in selecting patients for CRT. Thus, the aim of this study was to determine the value of cardiovascular magnetic resonance (CMR) for the assessment of the SDI in patients with reduced LV function as well as in healthy controls using semi-automatic border tracking.
We investigated a total of 45 patients including 35 patients (65 ± 8 years) with reduced LV function (EF 30 ± 11%) and a wide QRS complex as well as 10 control subjects (42 ± 21 years, EF 70 ± 11%). For cine imaging a standard SSFP imaging sequence was used with a temporal resolution of 40 frames per RR-interval. Quantitative analysis was performed off-line using a software prototype for semi-automatic border detection. Global volumes, ejection fraction and the SDI were calculated in each subject. SDI was compared with standard echocardiographic parameters of dyssynchrony.
The mean SDI differed significantly between patients (14 ± 5%) and controls (5 ± 2%, p < 0.001). An exponential correlation between the EF and the SDI was observed (r = -0.84; p < 0.001). In addition, a significant association between the SDI and the standard deviation of time to peak systolic motion of 12 LV segments (Ts-SD) determined by echocardiography was observed (r = 0.66, p = 0.002).
The results of this preliminary study suggest that CMR with semi-automatic border detection may be useful for the assessment of mechanical dyssynchrony in patients with reduced LV function.
No trial registration due to recruitment period between October 2004 and November 2006
Several studies have shown that the presence of intra-left ventricular (LV) dyssynchrony is one of the main predictors of response to cardiac resynchronization therapy (CRT) [1–3]. The surface electrocardiogram may not be optimal for patient selection because electrical dyssynchrony (wide QRS complex) does not accurately reflect intra-LV dyssynchrony [4, 5]. Tissue Doppler imaging (TDI) can assess LV dyssynchrony, i.e. the mechanical delay between septum and lateral wall contraction [2, 3, 6, 7]. In single center studies TDI as well as strain and strain rate imaging have been shown to be potentially useful in the selection of candidates for CRT [2, 3, 6–8]. However, these techniques have important technical limitations . In case of TDI regional systolic velocities may reflect passive motion of segments due to heart motion or tethering by adjacent segments. In addition, it is difficult to measure myocardial velocities in the apical segments of the LV.
More recently volumetric real time three-dimensional echocardiography has been used in conjunction with semi-automatic border detection algorithms to evaluate heart-failure patients for the presence of mechanical asynchrony [10, 11]. Myocardial wall motion has been shown to be accurately quantified by volume-time curve analysis with 3D analysis software . The systolic dyssynchrony index (SDI) defined as the standard deviation of the time to minimal systolic volume of all 16 segments was introduced as a measure of mechanical dyssynchrony . This index has been shown to be potentially useful for selection of patients who may benefit from CRT . However, three dimensional echocardiography may be inadequate in a number of patients due to limited temporal resolution and suboptimal acoustic window. In addition, no single echocardiographic measure of dyssynchrony has proved to be of value in selecting patients for CRT  and has been validated against mortality. Therefore, novel techniques to quantify mechanical dyssynchrony using cardiovascular magnetic resonance (CMR) may be of potential value.
Thus, the primary aim of this study was to compare the SDI in patients with reduced LV function and in healthy controls using cine CMR and semi-automatic border detection algorithms.
Patients and controls
The study group of this feasibility study consisted of 45 subjects including 35 heart failure (HF) patients (65 ± 8 years, 22 men) with reduced systolic LV function (EF 30 ± 11%) and wide QRS complex ≥ 120 ms (mean QRS duration 158 ± 22 ms). Ten normal subjects without cardiac disease (42 ± 21 years, EF 70 ± 11%) with normal QRS duration (83 ± 5 ms on electrocardiogram) were evaluated for comparison. The inclusion criteria for the HF patients were as follows: 1) heart failure New York Heart Association (NYHA) class III or IV 2) known ischemic (n = 13) or non-ischemic cardiomyopathy (n = 22) with at least a moderately impaired LV function (EF ≤ 40%) as assessed by CMR and 3) wide QRS complex ≥ 120 ms with left bundle branch block or interventricular conduction delay. Patients with acute coronary syndromes were excluded from the study. The patients were recruited between October 2004 and November 2006.
The ischemic or non-ischemic etiology of the cardiomyopathy was assessed using invasive cardiac catheterization. We performed cine CMR in all patients and control subjects. In addition, we retrospectively analysed echocardiographic data of the HF patients.
Informed consent was obtained from each patient prior to the investigation in accordance with local Ethics Committee requirements.
CMR protocol and generation of volume-time curves
Cine CMR was performed on a 1.5 T whole-body MR scanner (Intera, Philips Medical Systems, Best, The Netherlands) with electrocardiographic gating. Images were obtained using a five element cardiac synergy coil during repeated breath-holds (approximately 7 to 15 seconds). For cine imaging a standard SSFP imaging sequence (imaging parameters: echo time (TE) 1.8 ms, repetition time (TR) 3.6 ms, spatial resolution 1.8 × 1.5 mm2, slice thickness 8 mm) was used with a temporal resolution of 40 frames per RR-interval. Three long-axis images and a stack of 10 to 12 short axis images covering the complete LV were acquired.
The data set was transferred to a separate workstation and reconstructed off-line to a 4D volume data set. Quantitative analysis was performed using a software prototype for semi-automatic border detection (TomTec, Unterschleissheim, Germany). Endocardial borders were manually outlined in the three long-axis images at end-diastole and end-systole. Based on these initial contours LV casts are automatically calculated for each of the 40 time frames per RR interval providing a global volume-time-curve. Manual editing of endocardial borders in selected short-axis slices was possible. These changes in 2D images were subsequently adopted in the 3D data set. The cast were automatically subdivided into 17 subsegments resulting in 17 segmental volume-time curves (Fig. 1). The SDI was calculated from the time taken to reach minimum regional systolic volume for each segment as a percentage of the cardiac cycle. The SDI is defined as the standard deviation of these timings . Higher SDI denotes increasing intraventricular dyssynchrony. In addition, global volumes and ejection fraction were calculated in each patient.
We retrospectively analysed echocardiographic data off-line using commercial software (Echopac 6.1, General Electric-Vingmed) in 32 of 35 (91%) HF patients. These patients received a comprehensive echocardiographic study using a commercially available ultrasonographic system (GE Vingmed, Vivid 7, Horten, Norway). Interventricular mechanical dyssynchrony was assessed by determination of aortic and pulmonary pre-ejection times. A difference of more than 40 ms between the pulmonary and the aortic pre-ejection times (dPEP) was considered as significant delay . For assessment of intraventricular mechanical asynchrony, measurement of septal to posterior wall thickening delay (SPWMD) by M-mode echocardiography was performed using a cut-off of ≥ 130 ms indicating asynchrony . In addition, in 20 patients (57%) tissue Doppler images from three apical views (2, 3 and 4 chamber views) were available and the standard deviation of the time to peak systolic motion in 12 segments (Ts-SD) according to Yu et al was determined . A cut-off of Ts-SD = 34.4 ms was considered as measure for the presence of asynchrony .
Continuous variables are expressed as mean values ± standard deviation. Categorical data are presented as frequencies for the study population as well as for the control group and were compared using the Chi-Square test. The unpaired Student's t test was used for comparisons of continuous variables, i.e. for comparison of the ejection fraction, LV volume and SDI between patients and control subjects. The degree of association between SDI and EF as well as between SDI and echocardiographic parameters of asynchrony (Ts-SD, SPWMD, dPEP) was investigated by computing the Pearson's correlation coefficient. In addition, intraobserver and interobserver correlations were evaluated by using the Pearson's correlation coefficient.
Intraobserver variability of SDI was derived from manually outlined endocardial borders in the three long-axis images at end-diastole and end-systole from CMR of 10 randomly selected patients. Studies were performed by the same observer on 2 occasions 3 weeks apart. Interobserver variability was derived from the same subjects, performed by two blinded observers. Both intra- and interobserver variabilities were calculated as the standard deviation of the absolute difference between 2 measurements divided by the mean of both measurements, and expressed as a percentage. In Bland Altman analyses, the intraobserver and interobserver agreements, expressed in terms of the mean difference ± 2 standard deviation were evaluated.
The global significance level for all statistical tests conducted was chosen to α = 0.05. Statistical analysis was performed with the use of SPSS 10.0 (SPSS Inc., Chicago, USA).
Clinical characteristics and CMR results of patients and controls
The clinical baseline characteristics of the patient population with heart failure and the healthy control subjects are demonstrated in Table 1.
Data analysis using the semi-automatic software was feasible in all patients. Manual correction of endocardial contours had to be performed in most of the patients (n = 39, 87%). The mean time required for complete analysis was 10 minutes. However, major corrections defined by an additional time required for manual editing of endocardial borders in selected short-axis slices > 5 minutes (range 6–15 minutes) were needed in only a smaller proportion of patients (n = 10, 22%). Eight (80%) of these patients belonged to the HF group. In patients with HF the EF was significantly lower (EF 30 ± 11%) compared to the control subjects (EF 70 ± 11%, p < 0.001). In addition, LVEDV and LVESV measured by CMR were larger in HF patients than in controls (Table 2). The mean SDI was significantly increased in patients with HF (14 ± 5%, range 6.0%–23.4%) compared to control subjects (4.9 ± 1.5%, range 3.2%–7.8%; p < 0.001). An example of segmental volume-time-curves for a patient with marked dyssynchrony and a normal control subject is shown in Fig. 2. A significant exponential correlation between the mean SDI and EF was observed (r = -0.84, p < 0.001; Fig. 3). The correlation between mean SDI and EF for patients with ischemic cardiomyopathy and non-ischemic cardiomyopathy was r = -0.59 and r = -0.76, respectively. In addition, a strong correlation between the QRS complex duration and SDI was noted (r = 0.67, p < 0.001).
Intraobserver and interobserver variability for SDI was 2.2% and 3.4%, respectively. In addition, intraobserver and interobserver correlation was r = 0.99 and r = 0.94, respectively. In Bland Altman analyses, the intraobserver and interobserver agreements, expressed in terms of the mean difference ± 2 standard deviation (upper and lower limits of agreement) were 0.02 (-1.87 to 1,91%) and 0.16 (-3.45 to 3.77%), respectively.
A significant association between the SDI assessed by CMR and the standard deviation of time to peak systolic motion of 12 LV segments (Ts-SD) determined by Tissue Doppler imaging was observed (r = 0.66, p = 0.002, n = 20, Fig. 4). However, no significant associations between the SDI and SPWMD (r = 0.31, p = 0.11, n = 29) as well as between the SDI and dPEP (r = 0.04, p = 0.82, n = 32) were noted.
The results of this preliminary study show that myocardial wall motion can be quantified by volume-time curve analysis for assessment of mechanical dyssynchrony using CMR with semi-automatic border detection in analogy to 3D echocardiography. Using 3D imaging techniques it is possible to integrate the function of all 17 segments into the SDI which can be used to quantify the severity of LV dyssynchrony.
There is currently no gold-standard for the assessment of LV dyssynchrony. As electrical dyssynchrony does not accurately reflect intra-LV dyssynchrony, QRS-duration was found to be a poor predictor of CRT long term response [7, 15]. Ventricular mechanical dyssynchrony can be demonstrated by numerous echocardiographic parameters. However, these techniques have important technical limitations . In case of SPWMD, this parameter may not be useful in many patients, especially those with ischemic cardiomyopathy, in whom the inter-ventricular septum may be akinetic. This may partially contribute to the worse association between SDI and SPWMD observed in this study. The lack of correlation between SDI and dPEP may at least in part result from the fact that dPEP is a measure of interventricular dyssynchrony whereas the SDI reflect intra-LV dyssynchrony.
Real time 3D echocardiography has been shown to provide a precise evaluation of LV volumes even in conditions of altered load . Kapetanakis et al  used 3D echocardiography for guiding CRT and introduced the concept of the SDI. Despite the advantages of 3D echocardiography for evaluation of LV asynchrony, there are also limitations. LV volume analysis is image sector and quality dependent i.e. clear endocardial borders and complete LV cavity acquisition is required, which may be difficult in severely dilated left ventricles. Moreover, 3D echocardiography has a relatively low temporal resolution compared to other echo modalities such as M-Mode or TDI.
As no single echocardiographic measure of dyssynchrony has proven to be of value in selecting patients for CRT , CMR may be an alternative approach to select patients for CRT. In addition, late gadolinium enhancement (LGE) is able to determine the extent and transmurality of scar tissue. Recent data show that the total scar burden is an important factor influencing response to CRT . Thus, evaluation for viability, scar tissue and mechanical dyssynchrony at the same time using CMR may be important in the selection process for CRT candidates. Moreover, CMR is a powerful imaging modality for quantitative assessment of 3D myocardial structure und function. Measurements of end diastolic, end systolic volumes as well as ejection fraction are accurate and highly reproducible . A quantitative analysis of regional LV function has been performed by MR tagging and strain encoded imaging . However, the length of the post processing analysis procedure may often pose limitations to the clinical application.
In this CMR study a semi-automatic border detection prototype for quantitative analysis of LV dyssynchrony was used. The SDI values for normal subjects as well as for patients with HF in this study are in excellent agreement to previously published data using real-time 3D echocardiography . In addition, in HF patients an association between the SDI and the "echocardiographic dyssynchrony index", the Ts-SD, was shown in our study. Moreover, our study confirms the findings of Kapetanakis et al  that mechanical dyssynchrony is increasingly prevalent with worsening LV function which is supported by the reported correlation between the SDI and the EF. As 3D imaging not only allows to accurately evaluate EF but also demonstrates changes in regional myocardial motion by calculating regional volume-time curves, this may be important for preoperative planning before CRT. A recent CMR study evaluated a tissue synchronization index showing CMR to be a powerful method for assessing LV dyssynchrony and also for predicting mortality after CRT . In contrast to this study using semi-automatic border tracking, left ventricular volumes were quantified using manual planimetry of all short axis cine images in the former study . As manual planimetry of short axis images is possibly afflicted with inaccuracies and is more time consuming this poses limitations to its application in clinical routine. The ability of 3D imaging with CMR to discriminate between patients with heart failure and healthy control patients, to measure EF accurately, to demonstrate scar tissue and to assess SDI with high reproducibility from a 3D data set reflects its capability to image the entire left ventricle with superior spatial resolution. Thus, CMR may be a powerful method in the selection process for CRT candidates with the ability to identify LV dyssynchrony and to demonstrate scar tissue, which may prohibit response to CRT, at the same time.
The study population was relatively small. In addition, we did not assess a potential influence of infarct related wall thinning or motion abnormality on the calculation of the SDI. Thus, further prospective studies in larger populations also in comparisons to other imaging modalities are required to confirm this methodology to assess LV dyssynchrony and to compare with patient outcome.
The results of this feasibility study suggest that CMR with semi-automatic border detection may be useful for the assessment of mechanical dyssynchrony in patients with reduced LV function.
Yu CM, Chau E, Sanderson JE, Fan K, Tang MO, Fung WH, Lin H, Kong SL, Lam YM, Hill MR, Lau CP: Tissue Doppler echocardiographic evidence of reverse remodeling and improved synchronicity by simultaneously delaying regional contraction after biventricular pacing therapy in heart failure. Circulation. 2002, 105: 438-445. 10.1161/hc0402.102623.
Sogaard P, Egeblad H, Kim WY, Jensen HK, Pedersen AK, Kristensen BO, Mortensen PT: Tissue Doppler imaging predicts improved systolic performance and reversed left ventricular remodeling during long-term cardiac resynchronization therapy. J Am Coll Cardiol. 2002, 40: 723-730. 10.1016/S0735-1097(02)02010-7.
Bax JJ, Bleeker GB, Marwick TH, Molhoek SG, Boersma E, Steendijk P, Wall van der EE, Schalij MJ: Left ventricular dyssynchrony predicts response and prognosis after cardiac resynchronization therapy. J Am Coll Cardiol. 2004, 44: 1834-1840. 10.1016/j.jacc.2004.08.016.
Leclercq C, Faris O, Tunin R, Johnson J, Kato R, Evans F, Spinelli J, Halperin H, McVeigh E, Kass DA: Systolic improvement and mechanical resynchronization does not require electrical synchrony in the dilated failing heart with left bundle-brunch block. Circulation. 2002, 106: 1760-1763. 10.1161/01.CIR.0000035037.11968.5C.
Kass DA: Predicting cardiac resynchronization response by QRS duration: the long and short of it. J Am Coll Cardiol. 2003, 42: 2125-2127. 10.1016/j.jacc.2003.09.021.
Breithardt OA, Stellbrink C, Kramer AP, Sinha AM, Franke A, Salo R, Schiffgens B, Huvelle E, Auricchio A, PATH-CHF Study group. Pacing Therapies for congestive heart failure: Echocardiographic quantification of left ventricular asynchrony predicts an acute hemodynamic benefit of cardiac resynchronization therapy. J Am Coll Cardiol. 2002, 40: 536-545. 10.1016/S0735-1097(02)01987-3.
Penicka M, Bartunek J, De Bruyne B, Vanderheyden M, Goethals M, De Zutter M, Brugada P, Geelen P: Improvement of left ventricular function after cardiac resynchronization therapy is predicted by tissue Doppler imaging echocardiography. Circulation. 2004, 109: 978-983. 10.1161/01.CIR.0000116765.43251.D7.
Breithardt OA, Stellbrink C, Herbots L, Claus P, Sinha AM, Bijnens B, Hanrath P, Sutherlands GR: Cardiac resynchronization therapy can reverse abnormal myocardial strain distribution in patients with heart failure and left bundle branch block. J Am Coll Cardiol. 2003, 42: 486-494. 10.1016/S0735-1097(03)00709-5.
Rosen BD, Lardo AC, Berger RD: Imaging of myocardial dyssynchrony in congestive heart failure. Heart Fail Rev. 2006, 11: 289-303. 10.1007/s10741-006-0230-1.
Kapetanakis S, Kearney MT, Siva A, Gall N, Cooklin M, Monaghan MJ: Real-time three-dimensional echocardiography: a novel technique to quantify global left ventricular mechanical dyssynchrony. Circulation. 2005, 112: 992-1000. 10.1161/CIRCULATIONAHA.104.474445.
Park SM, Kim KC, Jeon MJ, Lee CK, Kim DH, Park KS, Lee WH, Kwan J: Assessment of left ventricular asynchrony using volume-time curves of 16 segments by real-time 3 dimensional echocardiography: comparison with tissue Doppler imaging. Eur J Heart Fail. 2007, 9: 62-67. 10.1016/j.ejheart.2006.04.009.
Chung ES, Leon AR, Tavazzi L, Sun JP, Nihoyannopoulos P, Merlino J, Abraham WT, Ghio S, Leclercq C, Bax JJ, Yu CM, Gorcsan J, St John Sutton M, DeSutter J, Murillo J: Results of the Predictors of Response to CRT (PROSPECT) trial. Circulation. 2008, 117: 2608-2616. 10.1161/CIRCULATIONAHA.107.743120.
Bax JJ, Ansalone G, Breithardt OA, Derumeaux G, Leclercq C, Schalij MJ, Sogaard P, St John Sutton M, Nihoyannopoulos P: Echocardiographic evaluation of cardiac resynchronization therapy: ready for routine clinical use? A critical appraisal. J Am Coll Cardiol. 2004, 44: 1-9. 10.1016/j.jacc.2004.02.055.
Pitzalis MV, Iacoviello M, Romito R, Massari F, Rizzon B, Luzzi G, Guida P, Andriani A, Mastropasqua F, Rizzon P: Cardiac resynchronization therapy tailored by echocardiographic evaluation of ventricular asynchrony. J Am Coll Cardiol. 2002, 40: 1615-1622. 10.1016/S0735-1097(02)02337-9.
Abraham WT, Fisher WG, Smith AL, Delurgio DB, Leon AR, Loh E, Kocovic DZ, Packer M, Clavell AL, Hayes DL, Ellestad M, Trupp RJ, Underwood J, Pickering F, Truex C, McAtee P, Messenger J, MIRACLE Study Group. Multicenter InSync Randomized Clinical Evaluation: Cardiac resynchronization in chronic heart failure. N Engl J Med. 2002, 346: 1845-1853. 10.1056/NEJMoa013168.
Atkins BZ, Kuo J, Shah AS, Hutcheson KA, Glower DD, von Ramm OT: Real-time three-dimensional echocardiography to construct clinically ready, load-independent indices of myocardial contractile performance. J Am Soc Echocardiogr. 2003, 16 (9): 922-930. 10.1016/S0894-7317(03)00515-7.
Ypenburg C, Roes SD, Bleeker GB, Kaandorp TA, de Roos A, Schalij MJ, Wall van der EE, Bax JJ: Effect of total scar burden on contrast-enhanced magnetic resonance imaging on response to cardiac resynchronization therapy. Am J Cardiol. 2007, 99: 657-660. 10.1016/j.amjcard.2006.09.115.
Doherty NE, Seelos KC, Suzuki J, Caputo GR, O'Sullivan M, Sobol SM, Cavero P, Chatterjee K, Parmley WW, Higgins CB: Application of cine nuclear magnetic resonance imaging for sequential evaluation of response to angiotensin-converting enzyme inhibitor therapy in dilated cardiomyopathy. J Am Coll Cardiol. 1992, 19: 1294-1302.
Lardo AC, Abraham TP, Kass DA: Magnetic resonance imaging assessment of ventricular dyssynchrony: current and emerging concepts. J Am Coll Cardiol. 2005, 46: 2223-2228. 10.1016/j.jacc.2005.09.015.
Chalil S, Stegemann B, Muhyaldeen S, Khadjooi K, Smith RE, Jordan PJ, Leyva F: Intraventricular dyssynchrony predicts mortality and morbidity after cardiac resynchronization therapy: a study using cardiovascular magnetic resonance tissue synchronization imaging. J Am Coll Cardiol. 2007, 50: 243-252. 10.1016/j.jacc.2007.03.035.
The authors thank Dr. Ole A. Breithardt, MD for his assistance in the preparation of the manuscript.
The authors declare that they have no competing interests.
RK conceived and designed the study, analyzed the data and drafted the manuscript. MN operated the CMR scanner and helped draft the manuscript. GS stored the data and participated in analyzing the CMR data. GAK operated the CMR scanner and helped draft the manuscript. SS performed the statistical analysis. RWG and MK participated in study design, coordination and scientific input. HPK participated in analysing the CMR data as well as in scientific input. All authors read and approved the final manuscript.
About this article
Cite this article
Koos, R., Neizel, M., Schummers, G. et al. Feasibility and initial experience of assessment of mechanical dyssynchrony using cardiovascular magnetic resonance and semi-automatic border detection. J Cardiovasc Magn Reson 10, 49 (2008). https://doi.org/10.1186/1532-429X-10-49