Skip to main content

Pilot study using 3D–longitudinal strain computation in a multi-parametric approach for best selecting responders to cardiac resynchronization therapy



Almost all attempts to improve patient selection for cardiac resynchronization therapy (CRT) using echo-derived indices have failed so far. We sought to assess: the performance of homemade software for the automatic quantification of integral 3D regional longitudinal strain curves exploring left ventricular (LV) mechanics and the potential value of this tool to predict CRT response.


Forty-eight heart failure patients in sinus rhythm, referred for CRT-implantation (mean age: 65 years; LV-ejection fraction: 26%; QRS-duration: 160 milliseconds) were prospectively explored. Thirty-four patients (71%) had positive responses, defined as an LV end-systolic volume decrease ≥15% at 6-months. 3D–longitudinal strain curves were exported for analysis using custom-made algorithms. The integrals of the longitudinal strain signals (I L,peak) were automatically measured and calculated for all 17 LV-segments.


The standard deviation of longitudinal strain peak (SDI L,peak ) for all 17 LV-segments was greater in CRT responders than non-responders (1.18% s−1 [0.96; 1.35] versus 0.83% s−1 [0.55; 0.99], p = 0.007). The optimal cut-off value of SDI L,peak to predict response was 1.037%.s−1. In the 18-patients without septal flash, SDI L,peak was significantly higher in the CRT-responders.


This new automatic software for analyzing 3D longitudinal strain curves is avoiding previous limitations of imaging techniques for assessing dyssynchrony and then its value will have to be tested in a large group of patients.

Peer Review reports


Cardiac resynchronization therapy (CRT) has emerged as a relevant therapeutic intervention for the treatment of chronic heart failure [1,2,3,4,5]. Based on current guidelines, patient selection for this costly therapy relies mainly on heart failure clinical status, ejection fraction (EF), and QRS characteristics (width and morphology). The proportion of non-responders remains relatively high, estimated at up to 30% [6, 7].

Several 2D echocardiographic indices of mechanical dyssynchrony have been proposed to better identify therapy responders, yet the lack of reproducibility and the non-optimal quantification of LV mechanical dyssynchrony (LVMD) achievable with these 2D echocardiographic indices in multicenter trials have cast some doubt on the techniques’ clinical applicability [8,9,10].

In a further advancement, 3D speckle-tracking echocardiography (STE) has been proposed as an alternative and potentially more accurate method for quantifying LVMD and for identifying patients suitable for CRT [11,12,13].

These different approaches are typically based on analyzing differences in either myocardial velocity timing, by means of tissue Doppler imaging (TDI), or in myocardial deformation using 2D/3D STE. To describe the complexity of LV mechanics, we believe it is essential to perform a combined assessment of LV dyssynchrony and LV regional contractility using STE, particularly by means of longitudinal strain analysis.

We hypothesized that a new approach, based on automatic quantification of the integrals pertaining to 3D regional longitudinal strain signals, could provide valuable additional information about regional LV mechanics and function prior to any CRT procedure. The aims of this pilot study were to describe LV mechanics using 3D echocardiography integral-derived longitudinal strain parameters in patients eligible for CRT and to test the relevance of this new tool for predicting CRT response.


Study population

48 patients referred for CRT device implantation at two institutions, the Saint Philibert Catholic University Hospital (Lille, France) and the Rennes University Hospital (Rennes, France), were included in the study. Indications for CRT implantation were based on the 2010 ESC guidelines for CRT device use in heart failure [14]. The patients had no previous pacemaker or cardioverter-defibrillator implantation. Patients with a poor acoustic window were excluded (n = 5).

Ischemic etiology was defined by a history of previous myocardial infarction or prior coronary revascularization or if a > 75% stenosis was observed in ≥1 of the major epicardial coronary arteries [15]. The NYHA functional class reported was the highest reached by the patient. 12-lead surface electrocardiograms (ECGs) were recorded at 25 and 50 mm/s during intrinsic conduction before CRT-device implantation and then were analyzed by Rennes University’s ECG Core Center. The morphology was classified as either LBBB or non-LBBB (non-specific intraventricular conduction delay) [16]. Only patients with a right bundle branch block were excluded.

The devices were implanted by a standard procedure. The electrophysiologist was blinded with respect to the localization of scar and the main aim during the implantation was to obtain the narrowest QRS at the end of the procedure.

All patients provided informed consent to participate in this study, which was performed in accordance with the principles outlined in the Declaration of Helsinki on research in human subjects (CNIL declaration no. 1620030 V. 0).

Two-dimensional echocardiography and speckle-tracking echocardiography

Baseline echocardiography was performed prior to CRT implantation (ViVid e9; GE Healthcare, Horten, Norway). Digital, routine, gray-scale, 2D Doppler TDI cineloops were obtained from three consecutive cardiac cycles, along with speckle tracking echocardiographic cineloops from one cardiac cycle, all from the apical view (gray-scale frame rate ≥ 50 Hz; color frame rate > 100 Hz; and 2-, 3-, 4-chamber apical views). Off-line analyses were performed on digitally stored images (BT12-EchoPAC PC; GE Healthcare).

The echocardiography examination was conducted according to the American Society of Echocardiography guidelines [17]. The LV volumes and LV ejection fraction (LVEF) were calculated using the biplane modified Simpson’s rule. All of the measurements were averaged for three cardiac cycles. The LV pre-ejection interval and aortic valve closure values were determined using aortic Doppler profiles. LV global longitudinal strain (LVGLS) was measured off-line using automated functional imaging (AFI). After manual tracing of the endocardial LV border in the 4-, 2-, and 3-chamber views over one frame, the endocardial borders were automatically tracked throughout the cardiac cycle. LVGLS was averaged from all of the analyzable segments in all apical views.

Atrio-ventricular delay was calculated as the ratio between LV filling time and the RR interval. Atrio-ventricular dyssynchrony (AVD) was considered significant when the duration of LV filling time resulted <40% of the RR interval [18]. Interventricular dyssynchrony was defined as a left ventricular pre-ejection interval (LVPEI) >140 ms, with or without an interventricular mechanical delay (IVMD) >40 ms [18]. Intraventricular dyssynchrony was defined by the presence of one of the following: a septal to lateral wall delay by color TDI [19] ≥65 ms [20], or the presence of septal flash (SF). SF was defined as an early septal thickening or thinning within the isovolumetric contraction period, as detected both visually from the gray-scale short axis (SAX) and 4-chamber (4CH) views and from the parasternal long axis, SAX, and 4CH views obtained by M-mode [21].

Three-dimensional echocardiography and speckle-tracking echocardiography

Baseline 3D–echocardiography was performed in each patient prior to CRT implantation using a commercially available echocardiographic system (ViVid E9; GE Healthcare, Horten, Norway), equipped with a 4 V phased-array matrix transducer. Consecutive 6-beat ECG-gated sub-volumes were acquired from the apical approach, using second-harmonic imaging during end-expiratory apnea, to generate the full-volume data set. We paid particular attention to encompassing the entire LV cavity within the data set, which was digitally stored in a raw-data format and was exported to a separate workstation equipped with the 4D–AutoLVQ package (EchoPAC V.110.1.3, GE Healthcare, Horten, Norway) for off-line analysis of STE LV myocardial longitudinal deformation.

The end-diastolic frames required for contour detection were automatically displayed in quad view. Manual alignment, achieved by pivoting and translating the four-chamber plane, was undertaken to align the three apical views so that the corresponding intersection line of all of the planes was placed in the middle of the LV cavity, crossing the LV apex and center of the mitral valve in each view. We subsequently used the semi-automated option to identify a fitting geometric model. Importantly, the software required only two single points to be input manually (one at the apex and the other at the tip of the mitral leaflet) on the end-diastolic and end-systolic frames of the four-chamber view slice. The software automatically detected the LV cavity endocardial border in 3D and provided the measured LV volumes. If the endocardial border detection was judged inadequate by the examiner, the LV endocardial borders were manually adjusted in multiplanar layout (three apical and three transverse planes) with a point-click method, immediately followed by secondary automated refinement of boundary detection according to the results. Following assessment of the LV volumes and ejection fraction, an automatic trace of the epicardial border was displayed to identify the region of interest required for LV mass and myocardial deformation measurements by means of 3D STE. This epicardial trace was manually adjusted in order to include the entire LV wall thickness using the same point-click method. The longitudinal deformation parameters were reported as global (both peak and end-systolic) and regional (only end-systolic) and were presented as color-coded polar maps and time-strain traces of an LV 17-segment model (Fig. 1). The time required for all this process is close to the time required to record the acquisitions and the export them on a computer where the computation will be performed automatically.

Fig. 1
figure 1

LV dataset display with 3D speckle-tracking analysis of longitudinal myocardial deformation, using the 4D–AutoLVQ package (EchoPAC version 110.1.3, GE Healthcare, Horten, Norway). Microsoft Excel files of 3D longitudinal strain analyses were exported for dedicated analysis performed, with Matlab software (Mathworks Inc., USA)

Automatic longitudinal strain integrals – new application

The 3D longitudinal strain signals were exported in text format to undergo dedicated analysis using custom-made methods and algorithms developed in Matlab software (Mathworks Inc., Natik, Massasuchetts, USA). This analysis included a regional characterization of the longitudinal strain integrals and extended to the 3D case of a set of methods initially developed by our group for 2D strain analysis [22]. The strain integrals represented the accumulated strain during different time intervals of the cardiac cycle. In this study, two particular markers were extracted for each of the 17 LV segments: 1) strain integrals from the beginning of the cardiac cycle (QRS onset) to the instant of the corresponding longitudinal strain peak (I L,peak ); and 2), strain integrals from the QRS onset to the instant of aortic valve closure (I L,avc ) (Fig. 2). All values exceeding −5% were considered to be noise (irrelevant information) and were thus not considered when calculating the integral.

Fig. 2
figure 2

Longitudinal strain curve of one LV segment, analyzed using custom-made algorithms. The pink area represents the integral of the longitudinal strain signal from the beginning of the cardiac cycle (QRS onset) to the instant of the aortic valve closure (I L,avc ). The blue-shaded area represents the integral of longitudinal strain signal from the beginning of the cardiac cycle to the instant of the corresponding longitudinal strain peak (I L,peak). All values greater than −5% were considered noise and were thus not considered in calculating integrals. t peak: time to strain peak

We tested the integral-based indicators of the regional longitudinal strain signals, all of which were automatically calculated, revealing standard parameters such as “peak strain” (amplitude), “mean strain peak”, and “SD t,peak ” (standard deviation of time to strain peak), along with novel measurements as detailed below.

  • mean I L,peak and mean I L,avc represented the mean of I L,peak and I L,avc of all 17 LV segments, respectively.

  • SDI L,peak and SDI L,avc were the standard deviation of the integrals of the strain signals I L,peak and I L,avc of all 17 segments, respectively. SDI L,peak and SDI L,avc corresponded to the energy dispersion for all 17 segments at the longitudinal strain peak and at the instant of aortic valve closure, respectively.

  • DiffInt was calculated as the average of all 17 LV segments in terms of the difference between I L,avc and I L,peak for each segment. This value could be considered an indicator of the wasted energy developed by the ventricle after aortic valve closure.

  • MSDI (maximal difference between strain peak instants) was calculated as the ratio of the time difference between the last and first strain peaks occurring during the cardiac cycle and the duration of the cardiac cycle.

Echocardiographic response

At 6-months post-implant, all of the patients were reassessed with 2D echocardiography. Response to CRT was defined as a ≥ 15% reduction in LV end-systolic volume, compared with baseline [23, 24].

Statistical analysis

All of the normally distributed data are displayed as the means and standard deviations, with non-normally distributed data expressed as medians and interquartile ranges (IQRs). Normality was evaluated by the Shapiro-Wilk test. Comparisons between the groups were performed using Student’s t-test or the Welch two-sample t-test, with the Mann–Whitney U test applied for normally and non-abnormally distributed data. Categorical variables, expressed as counts and percentages, were compared using the chi-square test or Fisher’s exact test, as appropriate. A p-value <0.05 was considered statistically significant. Sex, non-ischemic etiology, LBBB morphology, QRS ≥ 150 ms, GLS, atrioventricular dyssynchrony, IVMD >40 ms, LVPEI, DTI septo-lateral delay, septal flash, SD t,peak , SDI L,peak , DiffInt, and MSDI were covariates entered into the univariate model.

All variables showing a p < 0.1 were inserted in the multivariable logistic regression analysis (stepwise entry method). Variables with a p-value <0.1 in the multivariate model were considered possible contributors (according to the fact that this is a pilot study). For all of the variables in the multivariate model, the net odds ratio (OR) was reported, along with its 95% confidence interval (CI) and p-value. Receiver-operating characteristics curves were individually constructed for dyssynchrony parameters to determine the optimal threshold (closest to the top-left corner), sensitivity, specificity, positive and negative predictive values, and diagnostic accuracy. To explore the value in predicting response to CRT, of the tested 3D–strain parameters over the known variables (i.e. QRS width, or septal flash) the ROC curves were traced and a comparison was done. In addition, and for assessing the complementarity or the concordance of the new potential 3D–strain tool vs the previously proposed predictive parameters, a Kappa test was performed. The analyses were performed using R software, version 3.0.3 (R Foundation for Statistical Computing, Vienna, Austria. URL:


Study population

A total of 48 heart-failure patients (mean age: 65 ± 10 years, 30 men) referred for CRT device implantation were included in this study. 15 patients (31%) had ischemic cardiomyopathy. The mean intrinsic QRS duration was 160 ms (IQR: 160–170), and a typical LBBB morphology was observed in 39 (79%) patients. The mean LVEF was 26 ± 6%.

At 6-months follow-up, 34 34 patients (71%) had a LV end-systolic volume reduction ≥15% (responders). Compared to non-responders, the responders were more often female (47 versus 14%, p = 0.033), had and had better LV performance at baseline as indicated by LVEF (28 ± 5% versus 23 ± 5%, p = 0.002), and GLS values (−9.8 ± 3.4% versus 6.5 ± 3.1%, p = 0.003). The main clinical and echocardiographic characteristics are shown in Table 1.

Table 1 Baseline characteristics of patients

Classical dyssynchrony echocardiographic parameters

The classical dyssynchrony parameters were analyzed manually, with the results displayed in Table 2. No significant difference was observed in atrioventricular dyssynchrony between CRT-responders and non-responders. The prevalence of IVMD >40 msec and of septal flash were significantly higher in responders than in non-responders (91 versus 50%, p = 0.003 and 79 versus 21%, p = 0.001, respectively).

Table 2 Classical dyssynchrony 2D–echocardiographic parameters and 3D- echocardiographic integral-based indicators of longitudinal strain

Automatic analyses and 3D integral-derived longitudinal strain parameters

In the overall population, the median strain peak was −10.2% (−11.6; −9.2), with no significant difference between CRT responders and non-responders (−10.6% [−11.7;-9.6] versus −9.7% [−11.3;-8.4], p = 0.302). The same relationship was observed for SD t,peak (101 ms [80; 123] versus 107 ms [66; 121], p = 0.626).

In the study population, mean I L,avc was lower than mean I L,peak (the mean of the differences: 0.97% s−1, 95% CI: 0.82–1.13, p < 0.0001), and only Mean I L,peak was significantly higher in CRT responders than in non-responders (1.80 ± 0.62% s−1 versus 1.39 ± 0.41% s−1, p = 0.029).

The SDI L,peak of all 17 LV segments differed significantly between CRT responders and non-responders (1.18% s−1 [0.96; 1.35] versus 0.83% s−1 [0.55; 0.99], p = 0.007). DiffInt and MSDI parameters were comparable between CRT responders and non-responders (Table 2).

Predictors of echocardiographic response

From all of the clinical, electrocardiographic, and echocardiographic variables entered into the model, the univariate regression analysis identified six variables with a p-value <0.1(Table 3). The multivariate regression analyses identified three variables as potentials predictors (Table 3): septal flash (OR: 14.1; 95% CI: 3.08–64.9, p = 0.001), SDI L,peak (OR: 12.1; 95% CI: 0.81–180, p = 0.078), and non-ischemic etiology (OR: 5.33; 95% CI: 0.92–31.1, p = 0.063). Among the 18 patients without septal flash, the SDI L,peak values were higher in CRT responders than in non-responders (1.12 ± 0.26% s−1versus 0.77 ± 0.34% s−1, p = 0.03) (Table 4).

Table 3 Factors associated with good response to cardiac resynchronization therapy (univariate and multivariate regression analyses)
Table 4 3D–echocardiographic integral-based indicators of longitudinal strain in patients without Septal flash

The receiver operator characteristic curve analysis for SDI L,peak values identified an optimal cut-off value of 1.037% s−1, with a sensitivity (Se) of 70.6% and specificity (Sp) of 78.6% (Fig. 3). The positive predictive value (PPV), negative predictive value (NPV), and diagnostic accuracy for a cut-off of 1.037% s−1 were 89%, 52%, and 0.73, respectively (Table 5). For the septal flash, the Se, Sp, PPV, NPV, and diagnostic accuracy were 79%, 79%, 90%, 61%, and 0.79, respectively.

Fig. 3
figure 3

Receiver operator characteristic curve analyses to predict reduction in LVESV ≥15% after CRT for SDI L,peak

Table 5 Sensitivity, specificity, positive predictive predictive value, negative predictive value, diagnostic accuracy in monoparametric and multiparametric analyses for reverse remodeling induced by cardiac resynchronization therapy

All of the multiparametric associations were tested, and they are displayed in Table 5. Better diagnostic accuracy was achieved with a combination of an SDI L,peak value >1.037% s−1 and interventricular dyssynchrony. The association of four parameters (SDI L,peak  > 1.037% s−1, atrioventricular dyssynchrony, interventricular dyssynchrony, and septal flash) increased the test specificity to 100%, although it decreased its sensitivity and NPV. The Kappa test testing the concordance between septal flash and SDI L,peak was 0.26 confirming that SDI L,peak is exploring another kind of mechanical dyssynchrony than septal flash. In addition combining the septal flash and the SDI L,peak  > 1.037% s−1 lead to an area under the curve: AUC = 0.86.


Our study investigated a new quantitative computation method for longitudinal strain curves. Integral-derived parameters were generated and described for the first time, and they were applied in heart failure patients eligible for CRT, producing the following primary findings: 1) automatic analysis of longitudinal strain curves provided new complementary data on LV mechanics by combining information on timing and LV regional performance; 2) the mean I L,peak and SDI L,peak of all 17 LV segments were higher in CRT responders than in non-responders; 3) While marginally significant in this study, SDI L,peak could have an independent value for predicting CRT response; and 4) combining a SDI L,peak value >1.037% s−1 and interventricular dyssynchrony appeared to be a promising multi-parametric approach for best predicting CRT response, and it was feasible and robust in patients with sufficient acoustic windows and a sinus rhythm. And perhaps the most important finding was that 3D strain-derived parameter could differentiate responders from non-responders among 18 patients without septal flash.

Automatic versus mechanical manual dyssynchrony analysis

Since the publication of the disappointing PROSPECT study results [8], several questions have remained unanswered regarding the real discriminatory value of mechanical dyssynchrony in CRT candidates.

The first source of error is undoubtedly related to the lack of standardized data-processing methods. Numerous parameters have been proposed [10]. After disappointing results and thanks to effort of scientists, it has been best understood how complex is cardiac mechanical, how complex is electromechanical coupling but also how valuable could be the imaging approach to best select patients for CRT [25,26,27]. It has been for instance, clearly shown that, patients with the same left-bundle branch block could exhibit completely different mechanical dyssynchrony patterns [28]. But, there are also, technical challenges to consider. The analysis of strain peaks can be difficult because there are different patterns of strain, which often have multiple peaks. Determining the relevant peak or peaks, along with the reproducibility of the timing of these peaks, remains a challenge. To overcome this difficulty and also to consider not only peaks and dyssynchrony but regional myocardial function as well, a number of authors have proposed novel operator-independent methods that automatically assess function and dyssynchrony, using predefined algorithms [28,29,30,31]. In our study, heart failure patients with dilated or ischemic cardiomyopathy were considered eligible for CRT. All of the patients presented with severe LV dysfunction, an altered LVGLS as detectable on 2D echocardiography, an altered mean strain peak with a very low amplitude of strain curves, and significant dispersion of time to strain peak (SD t,peak ). We thus proposed a new 3D automatic assessment of regional LV mechanics, avoiding as much as possible potential human error.

Toward a new step by step approach for assessing left ventricular dyssynchrony assessment and for CRT response prediction

Mechanical dyssynchrony assessment could have 2 steps. As the first step, septal flash should be searched for using a simple visual and/or M-mode approach. If septal flash (and/or apical rocking) is not found, 3D echo should be performed in the next step as an attempt to detect a novel predictor of volumetric response. In the present study, we confirmed such results but also that the septal flash is highly relevant [21]. That is emphasizing the potential value of multiparametric scores like the L2ANDS2 score [32]. However, septal flash was found in 79% of CRT responders and in only approximately 60% of all of the patients, regardless of cardiomyopathy etiology. In patients without any septal flash, SDI L,peak was significantly higher in CRT responders, suggesting that this new predictor could provide additional information for the optimal selection of CRT patients. After the performing of a multiparametric evaluation of mechanical dyssynchrony that considers the three levels of dyssynchrony (atrio-ventricular, interventricular, and intra-ventricular), it appeared rather clear that this imaging approach could be of additive value to ECG [33]. A combination of very simple dyssynchrony parameters, such as atrioventricular, interventricular and septal flash, which are easy to measure and are reproducible, with automated strain-derived parameters, should probably be tested in larger groups of patients based on the current, first validation. The combination of an SDI L,peak  > 1.037% s−1 and interventricular dyssynchrony appeared to be an interesting multiparametric approach for predicting a good CRT response.

This “second step” three-dimensional STE appeared promising for several reasons [34]: (I) all 17 segments of the LV were evaluated in their 3D motion, along with the relationships among them, thus avoiding the ‘out-of-plane’ phenomenon inherent to 2D imaging; (II) the full LV volume was assessed during a 6-beat acquisition, allowing for a rapid evaluation of the global and regional all-directional contractions; and (III) all of the 3D echocardiographic strain markers (longitudinal, radial, and circumferential area) exhibited good reproducibility [35]. In our study, the feasibility rate of 3D STE was 83%, with the scientific literature reporting feasibility rates ranging from 63 to 83% [36]. This feasibility rate is likely to increase further with the advent of new transducer technologies. Until now, the principal 3D echocardiography parameter has been the standard deviation of time to minimal systolic volume [11]. This parameter has proved a feasible and reliable parameter of LV mechanical dyssynchrony, which might even provide additional value compared to the current selection criteria for accurate CRT response prediction [11]. This approach is only looking at the analysis of the differences in timings. The assessment of remaining LV regional contractility is lacking [11]. In a small study, 3D strain dyssynchrony index, (area tracking approach using the average difference between peak and end-systolic area strain, derived from 16 LV segments), was proposed for predicting CRT-response [13]. An integrative approach of dyssynchrony and function, using a new longitudinal strain integral-derived method, has appeared more relevant [29, 37]. This approach is likely to be of particular interest in ischemic diseases for distinguishing the passive movement of scarred tissue segments [38]. The mean I L,avc was lower than the mean I L,peak, in our patients, indicating that most segments reached their maximal deformation after aortic valve closure. SDI L,peak , corresponding to the energy dispersion for all 17 segments at the longitudinal strain peak, appeared to be promising for the assessment of LV dyssynchrony and the prediction of LV reverse remodeling following CRT.

Study limitations

That is a first pilot study, using one kind of echo-machine in a limited number of patients but in two centers. The goal was to ‘validate’ a tool that was planned to be use in largest population. Since a large intervendor variability has been demonstrated in 3D strain systems, the cut-off values identified are valid only for GE technology [39]. Only the integral-derived longitudinal strain parameters were measured in this study [40]. One limitation to the image acquisition for 3D speckle-tracking was the relatively slow volume rate of 32 ± 7 volumes/s. The volume rate and the image definition will improved in next generation of 3D–probes, it will improve the value of our automatic approach that will be applicable in atrial fibrillation patients thanks to the single beat capabilities.

Another limitation is related to fact that it is a pilot study with the use of an endpoint that is questionable. LV-remodeling is not a perfect surrogate marker of the hardest clinical endpoints. The Endpoint proposed by Packer will have to be considered [41].


This new automatic analysis of 3D longitudinal strain curves, using integral-derived parameters, provided original information on LV mechanics by combining timings and LV regional contractility data. This approach could be of value for improving patient selection for CRT.







Automated function imaging


Cardiac resynchronization therapy


The average of 17 LV segments of the difference between I L,avc and I L,peak for each segment

I L,avc :

The integrals of longitudinal strain signals for all 17 LV segments from the beginning of the cardiac cycle (QRS onset) to the moment of the aortic valve closure

I L,peak :

The integrals of longitudinal strain signals for all 17 LV segments from the beginning of the cardiac cycle (QRS onset) to the moment of the corresponding longitudinal strain peak


Interquartile range


Interventricular mechanical delay


Left bundle branch block


Left ventricular


Left ventricular ejection fraction


Left ventricular global longitudinal strain


Left ventricular pre-ejection interval


Maximal difference between strain peak instants


New York Heart Association


Positive predictive value

SD t,peak :

Standard deviation of time to strain peak


Speckle-tracking echocardiography


Tissue doppler imaging


  1. Cazeau S, Leclercq C, Lavergne T, Walker S, Varma C, Linde C, et al. Effects of multisite biventricular pacing in patients with heart failure and intraventricular conduction delay. N Engl J Med. 2001;344(12):873–80.

    Article  CAS  PubMed  Google Scholar 

  2. Cleland JG, Daubert JC, Erdmann E, Freemantle N, Gras D, Kappenberger L, et al. The effect of cardiac resynchronization on morbidity and mortality in heart failure. N Engl J Med. 2005;352(15):1539–49.

    Article  CAS  PubMed  Google Scholar 

  3. Linde C, Abraham WT, Gold MR, Daubert C, Group RS. Cardiac resynchronization therapy in asymptomatic or mildly symptomatic heart failure patients in relation to etiology: results from the REVERSE (REsynchronization reVErses Remodeling in systolic left vEntricular dysfunction) study. J Am Coll Cardiol. 2010;56(22):1826–31.

    Article  PubMed  Google Scholar 

  4. Moss AJ, Hall WJ, Cannom DS, Klein H, Brown MW, Daubert JP, et al. Cardiac-resynchronization therapy for the prevention of heart-failure events. N Engl J Med. 2009;361(14):1329–38.

    Article  PubMed  Google Scholar 

  5. Tang AS, Wells GA, Talajic M, Arnold MO, Sheldon R, Connolly S, et al. Cardiac-resynchronization therapy for mild-to-moderate heart failure. N Engl J Med. 2010;363(25):2385–95.

    Article  CAS  PubMed  Google Scholar 

  6. European Heart Rhythm A, European Society of C, Heart Rhythm S, Heart Failure Society of A, American Society of E, American Heart A, et al. 2012 EHRA/HRS expert consensus statement on cardiac resynchronization therapy in heart failure: implant and follow-up recommendations and management. Heart Rhythm. 2012;9(9):1524–76.

    Article  Google Scholar 

  7. European Society of C, European Heart Rhythm A, Brignole M, Auricchio A, Baron-Esquivias G, Bordachar P, et al. 2013 ESC guidelines on cardiac pacing and cardiac resynchronization therapy: the task force on cardiac pacing and resynchronization therapy of the European Society of Cardiology (ESC). Developed in collaboration with the European heart rhythm association (EHRA). Europace. 2013;15(8):1070–118.

    Article  Google Scholar 

  8. Chung ES, Leon AR, Tavazzi L, Sun JP, Nihoyannopoulos P, Merlino J, et al. Results of the predictors of response to CRT (PROSPECT) trial. Circulation. 2008;117(20):2608–16.

    Article  PubMed  Google Scholar 

  9. Ruschitzka F, Abraham WT, Singh JP, Bax JJ, Borer JS, Brugada J, et al. Cardiac-resynchronization therapy in heart failure with a narrow QRS complex. N Engl J Med. 2013;369(15):1395–405.

    Article  CAS  PubMed  Google Scholar 

  10. Gorcsan J 3rd, Abraham T, Agler DA, Bax JJ, Derumeaux G, Grimm RA, et al. Echocardiography for cardiac resynchronization therapy: recommendations for performance and reporting--a report from the American Society of Echocardiography Dyssynchrony writing group endorsed by the Heart Rhythm Society. J Am Soc Echocardiogr. 2008;21(3):191–213.

    Article  PubMed  Google Scholar 

  11. 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(7):992–1000.

    Article  CAS  PubMed  Google Scholar 

  12. Thebault C, Donal E, Bernard A, Moreau O, Schnell F, Mabo P, et al. Real-time three-dimensional speckle tracking echocardiography: a novel technique to quantify global left ventricular mechanical dyssynchrony. Eur J Echocardiogr. 2011;12(1):26–32.

    Article  PubMed  Google Scholar 

  13. Tatsumi K, Tanaka H, Matsumoto K, Hiraishi M, Miyoshi T, Tsuji T, et al. Mechanical left ventricular dyssynchrony in heart failure patients with narrow QRS duration as assessed by three-dimensional speckle area tracking strain. Am J Cardiol. 2011;108(6):867–72.

    Article  PubMed  Google Scholar 

  14. Dickstein K, Vardas PE, Auricchio A, Daubert JC, Linde C, McMurray J, et al. 2010 focused update of ESC guidelines on device therapy in heart failure: an update of the 2008 ESC guidelines for the diagnosis and treatment of acute and chronic heart failure and the 2007 ESC guidelines for cardiac and resynchronization therapy. Developed with the special contribution of the heart failure association and the European heart rhythm association. Eur Heart J. 2010;31(21):2677–87.

    Article  PubMed  Google Scholar 

  15. Felker GM, Shaw LK, O'Connor CM. A standardized definition of ischemic cardiomyopathy for use in clinical research. J Am Coll Cardiol. 2002;39(2):210–8.

    Article  PubMed  Google Scholar 

  16. Zareba W, Klein H, Cygankiewicz I, Hall WJ, McNitt S, Brown M, et al. Effectiveness of cardiac resynchronization therapy by QRS morphology in the multicenter automatic defibrillator implantation trial-cardiac resynchronization therapy (MADIT-CRT). Circulation. 2011;123(10):1061–72.

    Article  PubMed  Google Scholar 

  17. Lang RM, Badano LP, Mor-Avi V, Afilalo J, Armstrong A, Ernande L, et al. Recommendations for cardiac chamber quantification by echocardiography in adults: an update from the american society of echocardiography and the European association of cardiovascular imaging. Eur Heart J Cardiovasc Imaging. 2015;16(3):233–71.

    Article  PubMed  Google Scholar 

  18. Cazeau S, Bordachar P, Jauvert G, Lazarus A, Alonso C, Vandrell MC, et al. Echocardiographic modeling of cardiac dyssynchrony before and during multisite stimulation: a prospective study. Pacing Clin Electrophysiol. 2003;26(1 Pt 2):137–43.

    Article  CAS  PubMed  Google Scholar 

  19. Sogaard P, Egeblad H, Kim WY, Jensen HK, Pedersen AK, Kristensen BO, et al. Tissue Doppler imaging predicts improved systolic performance and reversed left ventricular remodeling during long-term cardiac resynchronization therapy. J Am Coll Cardiol. 2002;40(4):723–30.

    Article  PubMed  Google Scholar 

  20. Bax JJ, Bleeker GB, Marwick TH, Molhoek SG, Boersma E, Steendijk P, et al. Left ventricular dyssynchrony predicts response and prognosis after cardiac resynchronization therapy. J Am Coll Cardiol. 2004;44(9):1834–40.

    Article  PubMed  Google Scholar 

  21. Parsai C, Bijnens B, Sutherland GR, Baltabaeva A, Claus P, Marciniak M, et al. Toward understanding response to cardiac resynchronization therapy: left ventricular dyssynchrony is only one of multiple mechanisms. Eur Heart J. 2009;30(8):940–9.

    Article  PubMed  Google Scholar 

  22. Bernard A, Donal E, Leclercq C, Schnell F, Fournet M, Reynaud A, et al. Impact of Cardiac Resynchronization Therapy on Left Ventricular Mechanics: Understanding the Response through a New Quantitative Approach Based on Longitudinal Strain Integrals. J Am Soc Echocardiogr. 2015;28(6):700–708.

  23. Yu CM, Fung WH, Lin H, Zhang Q, Sanderson JE, Lau CP. Predictors of left ventricular reverse remodeling after cardiac resynchronization therapy for heart failure secondary to idiopathic dilated or ischemic cardiomyopathy. Am J Cardiol. 2003;91(6):684–8.

    Article  PubMed  Google Scholar 

  24. Bertini M, Hoke U, van Bommel RJ, Ng AC, Shanks M, Nucifora G, et al. Impact of clinical and echocardiographic response to cardiac resynchronization therapy on long-term survival. Eur Heart J Cardiovasc Imaging. 2013;14(8):774–81.

    Article  PubMed  Google Scholar 

  25. Galli E, Leclercq C, Donal E. Mechanical dyssynchrony in heart failure: still a valid concept for optimizing treatment? Arch Cardiovasc Dis. 2017;110(1):60–8.

    Article  PubMed  Google Scholar 

  26. Bernard A, Donal E, Leclercq C, Schnell F, Fournet M, Reynaud A, et al. Impact of cardiac resynchronization therapy on left ventricular mechanics: understanding the response through a new quantitative approach based on longitudinal strain integrals. J Am Soc Echocardiogr. 2015;28(6):700–8.

    Article  PubMed  Google Scholar 

  27. Brunet-Bernard A, Leclercq C, Donal E. Defining patients at-risk of non-response to cardiac resynchronization therapy. Value of rest and exercise echocardiography. Int J Cardiol. 2014;171(2):279–81.

    Article  PubMed  Google Scholar 

  28. Leenders GE, Lumens J, Cramer MJ, De Boeck BW, Doevendans PA, Delhaas T, et al. Septal deformation patterns delineate mechanical dyssynchrony and regional differences in contractility: analysis of patient data using a computer model. Circ Heart Fail. 2012;5(1):87–96.

    Article  PubMed  Google Scholar 

  29. Prinzen FW, Vernooy K, De Boeck BW, Delhaas T. Mechano-energetics of the asynchronous and resynchronized heart. Heart Fail Rev. 2011;16(3):215–24.

    Article  PubMed  Google Scholar 

  30. Russell K, Eriksen M, Aaberge L, Wilhelmsen N, Skulstad H, Gjesdal O, et al. Assessment of wasted myocardial work: a novel method to quantify energy loss due to uncoordinated left ventricular contractions. Am J Physiol Heart Circ Physiol. 2013;305(7):H996–1003.

    Article  CAS  PubMed  Google Scholar 

  31. Szulik M, Tillekaerts M, Vangeel V, Ganame J, Willems R, Lenarczyk R, et al. Assessment of apical rocking: a new, integrative approach for selection of candidates for cardiac resynchronization therapy. Eur J Echocardiogr. 2010;11(10):863–9.

    Article  PubMed  Google Scholar 

  32. Brunet-Bernard A, Marechaux S, Fauchier L, Guiot A, Fournet M, Reynaud A, et al. Combined score using clinical, electrocardiographic, and echocardiographic parameters to predict left ventricular remodeling in patients having had cardiac resynchronization therapy six months earlier. Am J Cardiol. 2014;113(12):2045–51.

    Article  PubMed  Google Scholar 

  33. Lafitte S, Reant P, Zaroui A, Donal E, Mignot A, Bougted H, et al. Validation of an echocardiographic multiparametric strategy to increase responders patients after cardiac resynchronization: a multicentre study. Eur Heart J. 2009;30(23):2880–7.

    Article  PubMed  Google Scholar 

  34. Tanaka H, Hara H, Saba S, Gorcsan J 3rd. Usefulness of three-dimensional speckle tracking strain to quantify dyssynchrony and the site of latest mechanical activation. Am J Cardiol. 2010;105(2):235–42.

    Article  PubMed  Google Scholar 

  35. Reant P, Barbot L, Touche C, Dijos M, Arsac F, Pillois X, et al. Evaluation of global left ventricular systolic function using three-dimensional echocardiography speckle-tracking strain parameters. J Am Soc Echocardiogr. 2012;25(1):68–79.

    Article  PubMed  Google Scholar 

  36. Jasaityte R, Heyde B, D'Hooge J. Current state of three-dimensional myocardial strain estimation using echocardiography. J Am Soc Echocardiogr. 2013;26(1):15–28.

    Article  PubMed  Google Scholar 

  37. Lim P, Donal E, Lafitte S, Derumeaux G, Habib G, Reant P, et al. Multicentre study using strain delay index for predicting response to cardiac resynchronization therapy (MUSIC study). Eur J Heart Fail. 2011;13(9):984–91.

    Article  PubMed  Google Scholar 

  38. Lim P, Pasquet A, Gerber B, D'Hondt AM, Vancraeynest D, Gueret P, et al. Is postsystolic shortening a marker of viability in chronic left ventricular ischemic dysfunction? Comparison with late enhancement contrast magnetic resonance imaging. J Am Soc Echocardiogr. 2008;21(5):452–7.

    Article  PubMed  Google Scholar 

  39. Badano LP, Cucchini U, Muraru D, Al Nono O, Sarais C, Iliceto S. Use of three-dimensional speckle tracking to assess left ventricular myocardial mechanics: inter-vendor consistency and reproducibility of strain measurements. Eur Heart J Cardiovasc Imaging. 2013;14(3):285–93.

    Article  PubMed  Google Scholar 

  40. Mor-Avi V, Lang RM, Badano LP, Belohlavek M, Cardim NM, Derumeaux G, et al. Current and evolving echocardiographic techniques for the quantitative evaluation of cardiac mechanics: ASE/EAE consensus statement on methodology and indications endorsed by the Japanese Society of Echocardiography. Eur J Echocardiogr. 2011;12(3):167–205.

    Article  PubMed  Google Scholar 

  41. Packer M. Development and evolution of a hierarchical clinical composite end point for the evaluation of drugs and devices for acute and chronic heart failure: a 20-years perspective. Packer M Circulation. 2016;134(21):1664–78.

    Article  PubMed  Google Scholar 

Download references


To the whole team and especially the research nurses at the CIC-IT 1414


For performing the study, we receive a grant from the ‘fondation de l’Avenir’ (Fondation de France). Grant ET2–654.

Availability of data and materials

All the data a recorded and saved.

Authors’ contributions

All the authors validated and provided correction to the manuscript. The concept was imagined by Alfredo Hernandez, Maxime Fournet, Anne Bernard and Erwan Donal. Elena Galli provided a substantial help. Christophe Leclercq helps as well and implanted all the patients. Sylvestre Marechaux provided patients and data. Philippe Mabo lead the research lab and agreed and reviewed the work as well as Jean-Claude Daubert and Raphael Martins. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

All the authors contributed substantially to the project and agreed for sending the manuscript to Cardiovascular Ultrasound. There is no competing interest and the data are available.

Ethics approval and consent to participate

All patients provided informed consent to participate in this study, which was performed in accordance with the principles outlined in the Declaration of Helsinki on research in human subjects (CNIL declaration no. 1620030 V. 0).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Erwan Donal.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fournet, M., Bernard, A., Marechaux, S. et al. Pilot study using 3D–longitudinal strain computation in a multi-parametric approach for best selecting responders to cardiac resynchronization therapy. Cardiovasc Ultrasound 15, 15 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: