Introduction

Cardiac resynchronization therapy (CRT) is an established treatment for patients with heart failure, reduced left ventricular (LV) ejection fraction, and a prolonged QRS caused by a left bundle branch block (LBBB) or nonspecific intraventricular conduction delay [1]. CRT aims to restore LV mechanics and improve hemodynamic by resynchronization of LV electrical activation [2]. Unfortunately, the effect of CRT is limited in 30–40% of the patients, partly due to a lack of optimal criteria for patient selection [3, 4]. In current international guidelines the selection criteria for CRT are limited to clinical parameters, ECG parameters and LV ejection fraction ≤ 35% [1]. Patient selection for CRT may be improved with additional parameters reflecting mechanical dyssynchrony or discoordination obtained with strain analysis on imaging [4,5,6,7]. These parameters reflect the LV mechanical consequences caused by an inhomogeneous electrical activation. Mechanical dyssynchrony parameters are based on timing differences between particular LV segments [8, 9]. However, these mechanical dyssynchrony parameters showed disappointing results in large multi-centre trials [9]. More promising parameters focus on discoordination, reflecting a percentage or fraction of opposing (i.e. inefficient) deformation [6, 10,11,12]. These parameters are determined using myocardial strain analysis, which can be obtained with several cardiac imaging techniques, including cardiac magnetic resonance imaging (CMR) with tagging (CMR-TAG), CMR cine images and a post-processing technique named feature tracking (CMR-FT), and speckle tracking echocardiography (STE) [13, 14]. Although CMR-TAG is regarded as the non-invasive ‘gold-standard’, it is generally limited to scientific applications, requiring specific imaging protocols, sequences, and dedicated post-processing software. Clinical application of CMR-FT and STE is more feasible compared to CMR-TAG, as both techniques are applicable to images obtained during standard clinical imaging protocols [14,15,16]. Nevertheless, both techniques (i.e. CMR-FT and STE) lack validation on strain parameters reflecting mechanical dyssynchrony and discoordination. Thus, no study has yet compared results obtained with all three techniques (i.e. CMR-TAG, CMR-FT and STE) in patients eligible for CRT. This study aims to compare circumferential strain parameters obtained with CMR-FT and STE versus gold-standard CMR-TAG in patients eligible for CRT. The comparison of indices reflecting mechanical dyssynchrony and discoordination are of specific interest.

Materials and methods

This sub study is part of the Markers of Acute Response to CRT (MARC) study (Cohfar, CTMM, The Netherlands, clinicaltrials.gov: NCT01519908), which was designed to investigate predictors for response on CRT. The MARC study included 240 patients planned for CRT implantation in six medical centres in the Netherlands, using previously published in- and exclusion criteria [17]. Twenty-seven of the 240 patients were included in this sub-study, as these patients gave consent for an additional CMR examination including myocardial tagging in the VU university medical centre (Amsterdam, The Netherlands). All subjects gave written informed consent and the local medical ethics committees approved data collection and management. The investigation conforms to the principles outlined in the Declaration of Helsinki.

Echocardiographic examination

Echocardiographic examinations were performed on either GE Vivid7, GE Vivid9 (General Electric Healthcare, Chicago, Illinois, USA), or Philips iE33 (Philips Medical Systems, Best, The Netherlands) ultrasound machines prior to CRT implantation by all participating centres and analysed by the echocardiographic core lab (WE and MC, UMC Utrecht, Utrecht, The Netherlands).

Acquisition—standard echocardiographic images

Standard echocardiographic images were obtained, including a parasternal short axis view at the papillary muscle and at the mitral valve level [18]. Image quality and frame rate (50–100 Hz) were optimized for offline speckle tracking analysis. Pulsed-wave Doppler images of the mitral valve inlet and LV outflow tract were obtained of mitral valve and aortic valve closure (AVC) to define systole.

Offline analysis—speckle tracking echocardiography

Echocardiographic images were exported as DICOM-files for vendor-independent strain analysis (TomTec 2D Cardiac Performance Analysis (2DCPA) version 1.2.1.2, TomTec Imaging Systems, Unterschleissheim, Germany). A region of interest was placed by user defined markers at the endocardial border. The epicardial border was excluded, as it often lacked a clear border zone. The region of interest was automatically separated into six segments. Segments were excluded if, even after repeated adjustment of the region of interest, adequate tracking was not achievable. The marker for reference length was placed at QRS onset.

STE results were exported for analysis with author written scripts for Matlab 2014b (Mathworks, Natick, MA, USA). Segmental strain curves were discarded in case of low signal-to-noise ratio as judged by two independent investigators (WE and AZ). At least two segments needed to be analysable per wall. Results of strain parameters of the septum were based on averages of maximal four septal segments (i.e. basal- and mid-level of inferoseptal and anteroseptal segments) while the lateral wall parameters were based on averages of maximal four lateral wall segments (i.e. basal- and mid-level of inferolateral and anterolateral segments). The post-processing and selection of analysable segments and averaging into one septal and one lateral wall strain curve, was similar for STE, CMR-TAG and CMR-FT.

Cardiac magnetic resonance imaging

CMR examinations were performed on a 1.5T system (Magnetom Avanto, Siemens, Erlangen, Germany) with the use of a phased array cardiac receiver coil. Although performed on a different moment compared to STE, both standard CMR cine images for the CMR-FT analysis and CMR tagging images were obtained in the same examination.

Acquisition—standard CMR images

Standard CMR cine images were acquired using a retrospectively ECG-gated balanced steady-state free-precession (SSFP) sequence during end-expiratory breath holding. A stack of 8–12 consecutive short axis cine images was acquired covering the entire LV. Typical image acquisition parameters were: slice thickness 5 mm, slice gap 5 mm, echo time (TE) 1.6 ms, repetition time (TR) 3.2 ms, temporal resolution < 50 ms, in-plane spatial resolution 1.5 by 2.1 mm, flip angle 60°. The number of reconstructed temporal phases within the cardiac cycle was set at 20. Subsequently, high temporal resolution (TE 1.7 ms, TR 3.4 ms, temporal resolution ~ 15 ms) cine imaging of the LV in the three-chamber view was performed to assess the opening and closure times of the mitral and aortic valve.

Acquisition—CMR tagging images

Before contrast injection, tagged images were acquired at three short-axis slices (basal, mid, apical) using a complementary spatial modulation of magnetization (CSPAMM) line tagging sequence with segmented ECG-gated acquisitions and serial breath holds [16]. Typical image acquisition parameters were: slice thickness 6 mm, TE 1.7 ms, TR 3.6 ms, temporal resolution < 15 ms, in-plane spatial resolution 1.3 by 4.3 mm, flip angle 20°, tag spacing 7 mm. The number of reconstructed temporal phases within the cardiac cycle was set at 55.

Offline analysis—CMR tagging

Tagged CMR images were exported and analysed (AZ, RN) with the SinMod technique (inTag, v2.0, CREATIS lab, Lyon, France, run as a plug-in for OsiriX Imaging Software v6.5, Pixmeo, Switzerland) [19]. Apical slices were discarded, while basal and mid short-axis slices were used for analysis in order to match the STE slice positions. After selecting the area of interest, endocardial and epicardial contours were manually drawn in the end-systolic phase and automatically propagated. A template was placed dividing the LV in six equally sized regions, similar to STE. The myocardium was divided in three layers (i.e. endo-, mid-, epi-wall layer). Results of the mid-wall layer were used, as these results are independent of contour placement.

Offline analysis—CMR feature tracking

Semi-automated FT analysis software (QStrain Research Edition evaluation version 1.3.0.10, Medis, Leiden, The Netherlands) was used to analyse short-axis cine images corresponding with the mid and basal slice-location of the CMR-TAG images (AZ and RN). Apical slices were discarded to match STE. First, endo- and epicardial contours were manually drawn in both end-diastolic and end-systolic frames and propagated automatically. Both endocardial and epicardial features were included for strain analysis, resulting in myocardial strain. The LV was divided in six regions, similar to the other techniques.

Basic strain parameters

The following parameters were obtained for the septal and lateral wall (Fig. 1). (1) Peak strain was the maximal negative peak strain during the cardiac cycle. (2) AVC strain was defined as the strain value at aortic valve closure. (3) Time to maximal peak (TTPmax) was the time difference between the start of the strain curve to most negative peak strain. Furthermore, (4) average systolic strain rate (i.e. average strain rate between mitral valve closure and AVC) and (5) average diastolic strain rate (i.e. average strain rate after AVC) were obtained.

Fig. 1
figure 1

Overview of imaging techniques and corresponding myocardial strain analysis. Examples of imaging techniques (top row) and resulting strain signals (bottom row) of one specific patient. Each column represents a single technique with the corresponding strain results. Examples of derived parameters are shown per graph. Basic strain parameters are indicated with a number, dyssynchrony and discoordination parameters are indicated with a character. Strain signals of the septum (black line) and lateral wall (grey line) are given, with the aortic valve closure (grey vertical line) as end of systole. CMR cardiac magnetic resonance imaging, AVC aortic valve closure, AVC strain strain value at aortic valve closure, TTP max time to maximal peak shortening, onset-delay time delay between onset of shortening of septal and lateral wall, peak-delay septal to lateral wall delay of TTPmax, SRS sept systolic rebound stretch of the septum, SSI systolic stretch index, ISF sep–lat internal stretch fraction of septal and lateral wall

Dyssynchrony parameters

Three parameters of dyssynchrony were analysed. (a) Onset-delay was determined as the absolute time delay between onset of shortening of the septal and lateral wall. (b) Peak delay was calculated as the absolute difference between lateral and septal wall TTPmax. (c) The TTPSD was calculated as the standard deviation of TTPmax of all analysable segments of the total LV.

Regional discoordination parameters

Three regional discoordination parameters were analysed. (d) Systolic rebound stretch of the septum (SRSsept) was defined as the total amount of systolic stretch after initial shortening of the septum (Fig. 1). (e) Systolic stretch index (SSI) was calculated by adding SRSsept to all systolic stretch of the lateral wall [11]. (f) Internal stretch factor (ISF) was calculated as the fraction of all systolic stretch compared to cumulative systolic shortening for the septal and lateral wall (ISFsep–lat). (g) Septal strain curves were categorized in three types, determined by their shape, LBBB-1: double-peaked systolic stretch, LBBB-2: early pre-ejection shortening peak followed by prominent systolic stretching and LBBB-3: pseudo normal shortening with a late-systolic shortening peak, followed by less pronounced end-systolic stretch (Fig. 2) [12].

Fig. 2
figure 2

LBBB pattern categorization. Septal strain pattern categorization and distribution of strain patterns found by the three imaging techniques. The distribution per imaging technique is given vertical in the upper panel. The cross-over of patients from CMR tagging to speckle tracking echocardiography and CMR feature tracking is displayed by arrows. The thickness of the arrows matches the number of patients crossing over. The number of patients crossing over is also given by a number in each arrow. Specific examples of the three patterns are given in the lower panel. Black curve: septal strain, grey dashed curved: lateral wall strain. CMR cardiac magnetic resonance imaging, LBBB-1 double peak shortening, LBBB-2 predominant stretch, LBBB-3 pseudo-normal shortening, n number of patients

Discoordination parameters of the total LV

Finally, two discoordination parameters reflecting the total LV were analysed. (h) The internal stretch factor of the total LV (ISFLV) was determined using all analysable segments. ISFLV was determined as the total amount of stretch divided by the total amount of shortening during systole (supplemental Fig. 1) [20]. (i) Lastly, the circumferential uniformity ratio estimates (CURE) was calculated, ranging from 0 (i.e. total dyssynchrony) to 1 (i.e. perfectly synchronous) [21].

Statistical analysis

Statistical analysis was performed (BG and MR) using R version 3.3.2 (The R foundation for Statistical Computing), and the R-packages psych version 1.5.8 (for calculation of Cohen’s kappa coefficients, ICCs and their associated p values). Results obtained with the three techniques were compared using the intra-class correlation coefficient (ICC) for absolute agreement between techniques (ICC2 according to Shrout and Fleiss) [22] and Spearman rank or Pearson correlation coefficient (R) depending on normality of data. An ICC ≥ 0.75 was classified as excellent, 0.60–0.74 as good, 0.40–0.59 as fair, and < 0.40 as poor [23]. Bland–Altman plots were made to observe the agreement between modalities. The mean difference and limits of agreement (± 1.96 standard deviation) of the Bland–Altman plot were used a reference of agreement. Lastly, Cohen’s kappa coefficient was calculated as the level of agreement between modalities on septal strain pattern categorization. A statistical result with a p value < 0.05 was deemed significant.

Results

Study population

Twenty-seven patients with CMR tagging images were included, of which a detailed description is given in Table 1. In these patients 94% of all segments were analysable with CMR-TAG, 87% with CMR-FT and 89% with STE. Frame rate of echocardiographic images was on average 65 ± 11 Hz, which corresponds to a temporal resolution of ~ 15 ms. Temporal resolution of CMR-TAG was ~ 14 ms, while it was ~ 40 ms for CMR-FT.

Table 1 Baseline characteristics

Basic strain parameters

Overall, agreement of CMR-TAG and CMR-FT was higher compared to agreement of CMR-TAG and STE for basic strain parameters. This applied for ICC values, Bland–Altman characteristics and the correlation coefficient (R) (Table 3). (1) For CMR-FT AVC strain of the septum was fair (ICC 0.55, R 0.67), while it was poor for STE (ICC 0.23, R 0.47). The ICC of AVC strain of the lateral wall was fair for CMR-FT (ICC 0.50, R 0.50) and poor for STE (ICC 0.08, R 0.10). (2) Peak strain of the septum had a fair ICC for CMR-FT (ICC 0.58, R 0.55) and a poor ICC for STE (ICC 0.155, 0.42). Lateral wall peak strain also had a fair ICC for CMR-FT (ICC 0.54, R 0.59) and a poor ICC for STE (ICC 0.01, R 0.02). (3) TTPmax of the septal and lateral wall showed an apparent wide distribution in the Bland–Altman plots for both comparisons (Fig. 3). Septal TTPmax caused a large spread in results, while the lateral wall TTPmax were more similar for both comparisons. Although ICC for TTPmax was poor for all comparisons, CMR-FT showed better agreement with CMR-TAG compared to STE for the septum (CMR-FT: ICC 0.17, R 0.11 and STE: ICC 0.00, R −0.16) and the lateral wall (CMR-FT: ICC 0.34, R 0.40 and STE: ICC 0.13, R 0.23). (4) Systolic strain rate showed comparable results to AVC strain. For CMR-FT systolic strain rate of the septum was fair (ICC 0.56, R 0.66), while it was poor for STE (ICC 0.25, R 0.45). ICC of the systolic strain rate of the lateral wall was fair for CMR-FT (ICC 0.575, R 0.58) and poor for STE (ICC 0.05, R 0.05). (5) Diastolic strain rate showed good ICC for the septal (ICC 0.64, R 0.66) and excellent ICC for the lateral wall (ICC 0.82, 0.82) for CMR-TAG versus CMR-FT. ICC of diastolic strain rate was poor for both walls comparing CMR-TAG and STE (septum: ICC 0.34, R 0.50, lateral wall: ICC 0.23, R 0.38).

Fig. 3
figure 3

Bland–Altman plots of basic strain parameters. Bland–Altman plots for CMR-TAG versus CMR-FT and CMR-TAG versus STE of three basic strain parameters. The mean of two techniques is plotted on the x-axis and the difference on the y-axis. The mean difference is displayed as a solid red line, while the limits of agreement are displayed as dotted red lines. Septal values are given as dots, while lateral wall values are given as crosses. AVC strain strain at aortic valve closure time, Peak strain highest negative peak strain value, TTP max time to maximal peak strain, CMR cardiac magnetic resonance imaging, TAG tagging, FT feature tracking, STE speckle tracking echocardiography

Table 2 Strain parameters of each myocardial strain analysis modality

Dyssynchrony parameters

(a) Onset delay was quite similar for CMR-TAG and CMR-FT, with a mean difference in the Bland–Altman plot of −2.5 ms (Supplemental Fig. 2). The corresponding ICC was fair (ICC 0.42, R 0.23). CMR-TAG versus STE also had a low mean difference of -1.9 ms, although the limits of agreement where larger, combined with a poor ICC (ICC 0.024, R −0.08). (b) Peak delay of CMR-TAG was overall larger compared to CMR-FT and STE (Supplemental Fig. 2). ICC was fair for CMR-FT (ICC 0.45, R 0.46), and poor for STE (ICC 0.23, R 0.27). TTPSD (c) showed a fair ICC for CMR-FT (ICC 0.46, R 0.49), and poor ICC for STE (ICC 0.20, R 0.19).

Regional discoordination parameters

(d) SRSsept showed a fair ICC for CMR-FT (ICC 0.41, R 0.65), while agreement was poor for STE (ICC 0.30, R 0.41). CMR-TAG showed overall higher values for SRSsept compared to both other imaging techniques. The difference of CMR-TAG to CMR-FT and STE were mostly positive, indicating an underestimation by CMR-FT and STE (Fig. 4). (e) SSI also showed an overall underestimation by CMR-FT and STE compared to CMR-TAG. Agreement on SSI was fair for CMR-FT (ICC 0.58, R 0.68) and STE (ICC 0.55, R 0.70). (f) ISFsep–lat was comparable between techniques, ICC’s of both CMR-FT (ICC 0.53, R 0.45) and STE (ICC 0.46, R 0.69) were fair. Overall values were still lower by CMR-FT and STE compared to CMR-TAG (Fig. 4). For septal strain patterns (g) the kappa value of CMR-TAG versus CMR-FT (0.465 p < 0.001) was higher compared to the kappa of CMR-TAG versus STE (0.265, p < 0.001). The number of patients crossing over from LBBB-1 or LBBB-2 on the one hand, and LBBB-3 on the other, using CMR-TAG and CMR-FT is rather low (n = 4, 15%), especially compared to CMR-TAG and STE (n = 8, 30%) (Fig. 2).

Fig. 4
figure 4

Bland–Altman plots of regional discoordination parameters. Bland–Altman plots for CMR-TAG versus CMR-FT and CMR-TAG versus STE of regional discoordination parameters (i.e. SRSsept, SSI and ISFsep–lat). The mean value of one patient analysed with the two techniques is plotted on the x-axis and the difference on the y-axis. The mean difference is displayed as a solid red line, while the limits of agreement are displayed as dotted red lines. SRS sept septal systolic rebound stretch, SSI systolic stretch index, ISF sep–lat internal stretch factor of septum and lateral wall, CMR cardiac magnetic resonance imaging, TAG tagging, FT feature tracking, STE speckle tracking echocardiography

Discoordination parameters of the total LV

ICC of ISFLV (h) of CMR-FT (ICC 0.55, R 0.66) was the highest off all dyssynchrony and discoordination parameters. Both ICC and R values were lower for STE (ICC 0.32, R 0.42). The CURE index (i) showed rather comparable values between techniques (Fig. 5) with relative narrow limits of agreement in the Bland–Altman plot. Both CMR-FT (ICC 0.485, R 0.37) and STE (ICC 0.41, R 0.36) resulted in a fair ICC value for CURE compared to CMR-TAG (Table 3).

Fig. 5
figure 5

Bland–Altman plots of discoordination parameters of the total LV. Bland–Altman plots for CMR-TAG versus CMR-FT and CMR-TAG versus STE of two discoordination parameters, obtained from the total LV. The mean value of one patient analysed with the two techniques is plotted on the x-axis and the difference on the y-axis. The mean difference is displayed as a solid red line, while the limits of agreement are displayed as dotted red lines. ISF LV internal stretch factor of the total LV, CURE circumferential uniformity ratio estimates, LV left ventricle, CMR cardiac magnetic resonance imaging, TAG tagging, FT feature tracking, STE speckle tracking echocardiography

Table 3 Intra-class correlation and correlation of CMR tagging versus CMR feature tracking and CMR tagging versus speckle tracking echocardiography

Discussion

This study explores the comparison of strain parameters in CRT candidates of two widely available strain analysis techniques, speckle tracking echocardiography and CMR feature tracking, with gold-standard CMR tagging. While most basic strain and dyssynchrony parameters differed substantially between techniques, there were apparent similarities found for discoordination parameters. This finding is promising, as discoordination parameters are potential predictors for CRT response [10, 11, 20]. The CMR-based techniques (i.e. CMR-TAG and CMR-FT) showed the highest agreement, shown in fair ICC values, higher R values, and relative narrow limits of agreement of the Bland–Altman plots. STE mostly had a poor agreement with CMR-TAG. CMR-FT may therefore be a valuable alternative to CMR-TAG for analysis of discoordination parameters in patients eligible for CRT.

Comparison of imaging techniques

To the best of our knowledge, this is the first study to compare strain parameters between different strain analysis techniques in a CRT patient population. The overall agreement between CMR-FT and CMR-TAG was higher compared to the agreement between STE and CMR-TAG. We would like to discuss three considerations to ascribe this difference. Firstly, STE uses a different imaging source, while both CMR-FT and CMR-TAG are obtained with the same imaging modality. Second, as part of the protocol, echocardiographic examinations and CMR scans were not performed on the same day. Therefore, physiological differences, such as loading conditions and heart rate, may have interfered with agreement of STE and CMR-TAG. Third, the imaging plane used for CMR and echocardiography is possibly different. Echocardiographic parasternal short-axis views were obtained from a single intercostal position, angulating the echo probe to the mitral valve annulus plane and the papillary muscle plane. These imaging planes may thereby be partly oblique, while CMR imaging planes were ‘true’ short-axis views. Furthermore, CMR-FT and CMR-TAG images were acquired on the almost exact same slice position, while the anatomical plane of STE images may be different. Another factor causing discrepancies between techniques is the specific manufacturer used for strain analysis with either CMR-TAG, CMR-FT, or STE [24]. Results of CMR-TAG, CMR-FT, and STE are contemporary, as they are dependent on specific analysis algorithms which are constantly under development. Although earlier studies show less favourable agreement of CMR-TAG and CMR-FT,[25, 26] recent developments are more promising [27, 28]. This trend is in accordance with our results, as we found that CMR-FT had fair agreement with CMR-TAG. However, further improvements are necessary, as results obtained with different imaging techniques can still differ largely for the individual patient. These differences may have underestimated the agreement between STE and CMR-TAG, compared to CMR-TAG versus CMR-FT. Nevertheless, echocardiography has its known limitations. High quality images are required for reliable strain analysis with STE,[29, 30] but can be difficult in this selection of patients. Frame rate is directly related to the temporal resolution, which is often high in echocardiographic images, especially compared to the relative low frame rate of standard cine images used for CMR-FT. A low frame rate causes under sampling and may lead to misinterpretation of peak and time-to-peak values in strain signals [29]. The frame rate of CMR-TAG was relatively high and comparable to STE in our study. Therefore, CMR-TAG may be considered a true gold-standard technique in this study, as imaging quality and frame rate of the implemented tagging protocol were optimized.

Assessment of strain parameters

Peak strain parameters showed fair correlation, especially between CMR techniques, except for timing indices of the septum. The maximal peak of septal strain can shift easily in case of dyssynchrony, as there are often multiple peaks (e.g. LBBB-1 and LBBB-2 patterns). Changes in absolute strain values of these peaks can drastically change TTPmax. In previous studies, most dyssynchrony and discoordination parameters have been primarily analysed with a single imaging technique. While some (i.e. CURE and ISFLV) are predominantly used in CMR-based studies,[20] others (i.e. SRSsept, peak-delay and septal strain patterns) are primarily derived with STE [12]. In our study, basic strain parameters, and more complex parameters of mechanical dyssynchrony showed apparent variations among the three techniques. However, the three techniques did show fair agreement on discoordination parameters. This indicates that these parameters adequately reflect mechanical discoordination and that they are detectable by multiple modalities. Discoordination parameters are promising as predictors for CRT response [10, 11, 20]. The predictive value of discoordination parameters is even known in combination with electrocardiographic parameters [6, 11]. ISFLV and CURE are predictors of CRT response and use information of all available LV segments,[31] therefore reflecting total LV discoordination [20, 21]. These parameters are also less susceptible to outliers compared to basic strain parameters, as they contain information on all segments [20, 21]. Parameters being calculated using averages of multiple segments (i.e. SRSsept, SSI) also showed fair agreement between modalities. Obtaining deformation characteristics using averages of multiple segments may therefore reduce noise and measurement variability. Specific pre-specified septal strain patterns are known to predict CRT response, as LBBB-1 and LBBB-2 patterns are associated with volumetric response after CRT, while LBBB-3 is not [12, 32]. The relative high agreement between CMR-TAG and CMR-FT on LBBB-1 and LBBB-2 on the one hand, and LBBB-3 on the other is therefore promising for further implementation of septal strain pattern categorization using CMR.

Myocardial strain orientation

STE parameters are mainly validated with longitudinal strain,[33, 34] while CMR is predominantly based on circumferential strain [20, 21]. Circumferential strain is more intuitive, as mid-myocardial fibres are orientated in the circumferential direction and short-axis images represent all segments distributed around the LV at each level (i.e. basal, mid or apical) [35]. The method of determining circumferential strain calculation differs between the three methods. Both the CMR-FT and STE software track specific myocardial details, respectively ‘features’ and ‘speckles’, of the endo- and epicardial border [14]. The specific wall layer used for strain analysis differed between techniques. The results of the endocardial layer were used for STE, as the epicardial layer often lacked an appropriate border zone. Strain values of CMR-FT were a product of endocardial and epicardial strain. This is in contrast to CMR-TAG, of which strain of the mid-wall layer was used [19]. The difference between the approaches may have biased the overall level of agreement. Endocardial strain is known to give higher peak values compared to epicardial strain [36], and might also be higher than midmyocardial values, which can be appreciated in the positive mean difference between CMR-TAG and STE on peak strain and AVC-strain in the Bland–Altman results. This difference may have also affected the agreement of dyssynchrony and discoordination parameters.

Limitations

CMR imaging with myocardial tagging was performed in a small subset of patients from the MARC study, which may have given outliers a relatively large effect on results. The patient population was moreover limited to patients eligible for CRT, reducing variability in measurements. These results should therefore be validated in a larger cohort. However, strain measurements are of particular interest in this specific population to improve patient selection for CRT. As mentioned, the study protocol has also influenced results, as echocardiographic and CMR examination were not performed on the same day. Moreover, differences in imaging plane between CMR and STE are possible and strain analysis was not performed on the same wall layers. ECG triggering differs between imaging techniques, as ECG electrodes were repositioned between examinations and a different lead may have been used. Moreover, ECG triggering of STE was placed at QRS onset, while the top of the R wave is used for CMR. ECG triggering affects the reference value and may have affected subsequent values of timing and absolute changes. While STE relied on end-diastolic region of interest placement, CMR-FT used both end-systolic and end-diastolic region of interests to determine myocardial strain. The reliability of CMR-FT may therefore be higher. Echocardiography was moreover obtained with ultrasound machines from two vendors, possibly introducing differences in source data. The overall lower agreement of CMR-TAG and STE should therefore be appreciated carefully. STE was performed with circumferential strain obtained from short axis images, for a more direct comparison between techniques. While circumferential strain is widely used in scientific publications, standardization of algorithms of STE has also mainly been done for longitudinal strain [5, 34]. Longitudinal strain assessed with STE may therefore have a higher reliability and reproducibility compared to circumferential strain. The effect of using longitudinal or circumferential strain derived with STE for prediction of CRT response deserves attention in future work.

Clinical application

The overall reasonable agreement between CMR-TAG and CMR-FT is promising for clinical application. CMR-FT might be a reasonable alternative for CMR-TAG and STE, as suitable CMR cine images are more easily available in clinical practice, compared to the highly specialized CMR-TAG protocols. Detection of mechanical discoordination with CMR-FT is a valuable addition to CMR, which already constitutes an important imaging tool in CRT-candidates for accurate determination of the LV ejection fraction and scar tissue localization [31]. On the other hand, a portable and bed-side tool like STE might have the highest clinical applicability, of which most discoordination parameters also showed fair agreement compared to the CMR-TAG. The reasonable agreement of the three techniques on mechanical discoordination parameters is moreover promising for the prediction of response to CRT. The implemented discoordination parameters were previously associated with CRT response in single centre studies [10, 11, 20]. However, previous markers of CRT response failed to take the final step to clinical application, partly because validation to gold-standard techniques was missing [9]. As the specific methods and modality may slightly differ from previous publications, further studies are needed for implementation into clinical practice. Future studies will focus on the predictive value of these parameters using follow-up data in this specific population.

Conclusions

In conclusion, comparison of strain analysis techniques showed that CMR-FT had an overall fair agreement with gold-standard CMR-TAG. Although agreement between STE and CMR-TAG was overall lower, direct comparison was limited by technical and methodological differences. The agreement was highest for parameters of mechanical discoordination, compared to basic strain or dyssynchrony parameters. CMR-FT is therefore a potentially valuable clinical alternative for CMR-TAG and STE, especially in the evaluation of mechanical discoordination in CRT-candidates.