Layer-specific strain in patients with heart failure using cardiovascular magnetic resonance: not all layers are the same
Journal of Cardiovascular Magnetic Resonance volume 22, Article number: 81 (2020)
Global longitudinal strain (GLS), most commonly measured at the endocardium, has been shown to be superior to left ventricular (LV) ejection fraction (LVEF) for the identification of systolic dysfunction and prediction of outcomes in heart failure (HF). We hypothesized that strains measured at different myocardial layers (endocardium = ENDO, epicardium = EPI, average = AVE) will have distinct diagnostic and predictive performance for patients with HF.
Layer-specific GLS, layer-specific global circumferential strain (GCS) and global radial strain (GRS) were evaluated by cardiovascular magnetic resonance imaging (CMR) feature tracking in the Alberta HEART study. A total of 453 subjects consisted of healthy controls (controls, n = 77), at-risk for HF (at-risk, n = 143), HF with preserved ejection fraction (HFpEF, n = 87), HF with mid-range ejection fraction (HFmrEF, n = 88) and HF with reduced ejection fraction (HFrEF, n = 58). For outcomes analysis, CMR-derived imaging parameters were adjusted with a base model that included age and N-terminal prohormone of b-type natriuretic peptide (NT-proBNP) to test their independent association with 5-year all-cause mortality.
GLS_EPI distinguished all groups with preserved LVEF (controls − 16.5 ± 2.4% vs. at-risk − 15.5 ± 2.7% vs. HFpEF − 14.1 ± 3.0%, p < 0.001) while GLS_ENDO and all GCS (all layers) were similar among these groups. GRS was reduced in HFpEF (41.1 ± 13.8% versus 48.9 ± 10.7% in controls, p < 0.001) and the difference between GLS_EPI and GLS_ENDO were significantly larger in HFpEF as compared to controls. Within the preserved LVEF groups, reduced GRS and GLS_EPI were significantly associated with increased LV mass (LVM) and LVM/LV end-diastolic volume EDV (concentricity). In multivariable analysis, only GLS_AVE and GRS predicted 5-year all-cause mortality (all ps < 0.05), with the strongest association with 5-year all-cause mortality by Akaike Information Criterion analysis and significant incremental value for outcomes prediction beyond LVEF or GLS_ENDO by the likelihood ratio test.
Global strains measured on endocardium, epicardium or averaged across the wall thickness are not equivalent for the identification of systolic dysfunction or outcomes prediction in HF. The endocardium-specific strains were shown to have poorest all-around performance. GLS_AVE and GRS were the only CMR parameters to be significantly associated with 5-year all-cause mortality in multivariable analysis. GLS_EPI and GRS, as well as the difference in endocardial and epicardial strains, were sensitive to systolic dysfunction among HF patients with normal LVEF (> 55%), in whom lower strains were associated with increased concentricity.
Heart failure (HF) is a complex clinical syndrome with a wide array of characteristic cardiac structural and functional abnormalities . Patients are typically categorized according to left ventricular (LV) ejection fraction (LVEF), with HF and reduced ejection fraction defined as LVEF < 40% (HFrEF), HF with mid-range LVEF in the span of 40–49% (HFmrEF) and those with HF with preserved LVEF ≥ 50% (HFpEF) [2, 3]. However, LVEF has an inconsistent relationship with outcomes [4, 5] and does not distinguish those with preserved LVEF from those without cardiac disease . Global longitudinal strain (GLS) provides an alternate measure of systolic dysfunction that has been shown to be superior to LVEF in distinguishing HF patients from healthy subjects  and predicting adverse outcomes in acute HF .
However, myocardial structure and function are heterogeneous, with layer-specific fiber orientations ranging from largely circumferential at the mid-wall to more oblique at the endocardium (ENDO) and epicardium (EPI) . There is a gradient in myocardial strain across the wall, with decreasing values from ENDO to EPI [9, 10] which makes reported values dependent on measurement layer, with endocardial values being most commonly reported [6, 7, 11, 12]. The potential added value of layer-specific strain has been illustrated in various cardiovascular diseases. In those with suspected coronary artery disease, ENDO strains were shown to be superior to EPI strains for diagnosis . Other studies of ischemic heart disease have shown the superior performance of ENDO strain for the prediction of outcomes [14, 15]. Conversely, EPI longitudinal strain has been shown to provide incremental prognostic information for acute coronary syndrome  and hypertension , and improved diagnostic performance for myocarditis . EPI, mid-wall and endocardial GLS were shown to have similar predictive performance in population-based cohort . Layer-specific strains were recently reported in a small HF cohort .
The aims of the current study were to compare the utility of global strains assessed at different layers of the myocardium, measured with cardiovascular magnetic resonance imaging (CMR) feature-tracking, to characterize layer-dependence, to detect systolic dysfunction in HFpEF and to predict outcomes in all HF patients, regardless of LVEF.
The study was approved by the University of Alberta and University of Calgary Health Research Ethics Boards and all study participants provided written informed consent. We excluded those unable to provide informed consent or with a contraindication to CMR. Patient recruitment and testing has previously been reported . Briefly, we recruited patients with chronic HF and those at-risk for HF from ambulatory clinics and all subjects underwent comprehensive phenotyping that included a detailed history and physical examination, serum biomarkers and a multi-parametric CMR exam. Individuals at-risk for HF had a history of coronary artery disease, diabetes mellitus, hypertension, atrial fibrillation, and/or obesity without a diagnosis of HF (American Heart Association (AHA)/American College of Cardiology (ACC) class A and B) . HF patients (AHA/ACC Class C), were sub-grouped into those with preserved (HFpEF, LVEF ≥ 55%) [21, 22], midrange (HFmrEF, 40% ≤ LVEF < 55%) or reduced ejection fraction (HFrEF, LVEF < 40%) [2, 3]. Healthy controls were also recruited and underwent identical testing.
All subjects underwent a CMR examination on 1.5 T CMR scanners (Siemens Sonata or Avanto, Siemens Healthineers, Erlangen, Germany) at the University of Alberta or University of Calgary sites, respectively. To quantify cardiac structure and function, the CMR exam included standard balanced steady state free precession (bSSFP) cine imaging, with 10–14 short axis slices covering the entire ventricle, as well as two-chamber, three-chamber and four-chamber long axis views. Typical acquisition parameters: Repetition time/echo time (TR/TE) 2.8/1.4 ms, 50–70 degree flip angle, 8 mm slice thickness with a 2 mm gap for short axis slices, 256 × 192 matrix, 380 × 285 mm field of view, 10 views per segment with 25 or 30 reconstructed cardiac phases over the cardiac cycle. All cardiac images were acquired with electrocardiographic gating within an 8–12 s breath-hold per slice.
LV and right ventricular (RV) end-diastolic volume (LVEDV, RVEDV) and end-systolic volume (LVESV, RVESV), ejection fraction (LVEF, RVEF) and left ventricular mass (LVM) were measured using Syngo Argus, (Siemens Healthineers) or CVI42 (Circle Cardiovascular Imaging, Calgary, Canada). Both methods used only short axis slices for volume analysis and papillary muscle was included in the chamber volume. Volumes and LVM were normalized to body surface area (BSA), calculated using idealized weight for the given subject height . Relative wall thickness (RWT) was calculated from short-axis slices as average end-diastolic wall thickness divided by average LV end-diastolic diameter from two short-axis slices (mid and mid-basal slice locations). LVM/LVEDV (concentricity) was also calculated.
Strain was measured from bSSFP short and long axis cine images using a feature tracking approach develop in-house based on b-spline non-rigid registration . The strain evaluation procedure is summarized in Fig. 1. First, in-plane displacement fields were calculated from registration of all images in each cine image series to the reference end-diastolic image frame, separately for each short-axis and long-axis slice (Fig. 1a and b). Manual tracing was limited to tracing of endocardial (red) and epicardial (green) borders (Figs. 2 and 3) on the end-diastolic image frames (Fig. 1c). Tracings were completed by an experienced interpreter for all subjects (LX) who was blinded to clinical data and conventional CMR measures. Equally spaced contours between the endocardial and epicardial contours (blue) were generated at the end-diastolic frame (a subset of longitudinal contours are shown in Figs. 2 and 3 and 1d). A similar automated procedure was repeated for the radial direction. Subsequently, all end-diastolic contours were automatically propagated to all image frames over the full cardiac cycle using the previously calculated feature tracking displacement fields (Fig. 1d). Strain in each slice was calculated independently for all contours as the fractional change in length of the contour from end-diastole (L0) to end-systole (L) relative to end-diastolic length, reported as a percentage, (L − L0)/L0*100 (Lagrangian strain), as shown in Fig. 1d [25, 26]. GLS was calculated separately for the endocardium (GLS_ENDO), the epicardium (GLS_EPI) and the average of all contours (GLS_AVE). Peak systolic strain from all long axis slices were averaged to provide the reported global values (Fig. 2, right panels). The ratio, GLS_ENDO/GLS_EPI, as well as the absolute and relative differences between GLS_ENDO and GLS_EPI were also calculated for each subject. Similarly, global circumferential systolic strains (GCS_ENDO, GCS_EPI, GCS_AVE) were calculated as the average of the peak strains from two short-axis slices (mid and mid-basal) (Fig. 3). The ratio, GCS_ENDO/GCS_EPI, as well as the absolute and relative differences between GCS_ENDO and GCS_EPI were also calculated for each subject. Finally, global radial strains (GRS) were calculated separately from both long axis (GRS_LAX) and short axis slices (GRS_SAX). Contour lengths for the calculation of radial strains were measured as the distance between the endocardial and epicardial borders at end-diastole and end-systole (Fig. 1d). Radial strain for each slice was calculated as the average radial strain from all contours, and GRS as the average from all slices, for the two short-axis (GRS_SAX) and three long-axis slices (GRS_LAX), respectively.
GLS and GCS strains at all layers were also calculated using commercially available CMR feature tracking software (CVI42, Circle Cardiovascular Imaging) in a subset of 202 subjects, for a comparison of feature tracking methodology used in current study.
Continuous variables were expressed as mean ± standard deviation or median (25th, 75th percentile), as appropriate. Categorical variables were expressed as frequency and percentage. For missing data, the data was assumed to be missing at random. Our cohort had missing values for N-terminal prohormone of b-type natriuretic peptide (NT-proBNP, 11.0% in total, including 2.3% in healthy control, 3.5% in patients at risk for HF, 5.2% for patients with HF) and creatinine (8.6% % in total, including 3.1% in healthy control, 1.3% in patients at risk, 4.2% in HF patients). Multiple imputation by chained equations with 50 imputed data sets was used to generate missing data based on all candidate predictors and outcomes. We averaged results from the 50 imputations .
Two sample t-test (or Mann–Whitney U test) or one-way analysis of variances with Bonferroni post-hoc correction (or Dunn’s test) were used to compare continuous variables among groups of patients, as appropriate. Chi-square was used to compare the categorical baseline characteristics. The normal distribution of continuous variables was tested by the Shapiro–Wilk normality test. We applied logarithmic transformation to NT-proBNP and creatinine. Correlation between continuous variables was tested by Pearson correlation.
Clinical outcome was all-cause mortality over a 5-year follow-up. Univariable Cox proportional regression was performed for all demographic parameters, cardiovascular risk factors, cardiovascular disease history, concomitant diseases and CMR-derived imaging parameters. Sex and the parameters with univariable p value < 0.2, including age, systolic blood pressure, current smoker, presence of heart failure, coronary artery disease, atrial fibrillation or atrial flutter, chronic obstructive pulmonary disease, log (NT-proBNP) and log (creatinine) entered the forward selection approach based on Akaike Information Criterion (AIC) to build the optimal set of predictors as the base model for adjustment. In the multivariable Cox proportional hazard analysis conventional LV measurements, GRS, and layer-specific strains for GLS, GCS were each adjusted with the base model to test their independent association with the outcome. Additionally, the association of layer-specific strains and conventional LV measurements with clinical outcomes were quantified using the AIC values, where the lowest AIC score (AICminimum) indicates the best outcomes model. Each compared model included the base model and a single evaluated CMR parameter. For comparison of AIC values (between CMR parameters), Δi = AICi – AICminimum is the difference between Modeli and best fit model. If Δi ≤ 2, there is substantial support for similar predictive performance of the best fit model and the i-th model; for 4 ≤ Δi ≤ 7 there is considerably less support for i-th model; for Δi > 10 there is no support for i-th model . Finally, the incremental value of outcomes prediction by layer-specific strain over the other CMR variables was evaluated using the likelihood ratio test.
Intra-observer and inter-observer reproducibility of all strain parameters were evaluated by intra-class correlation coefficient and coefficient of variation in 20 randomly selected subjects with blinding to clinical data (the same subjects were used for both intra-observer and inter-observer calculations). Statistical analyses were performed using STATA (version 16.0, StataCorp LP, College Station, Texas, USA). A p-value less than 0.05 was considered significant for all tests.
From 466 subjects who underwent CMR exams, 453 subjects with analyzable images were identified: controls (n = 77), at-risk (n = 143), HFpEF (n = 87), HFmrEF (n = 88), and HFrEF (n = 58).
For statistical analysis of group differences, two subgroups of subjects were considered, those with preserved LVEF (> 55%), consisting of healthy controls, patients at-risk for HF and HFpEF (p-value 1 in Tables 1 and 2), and the three HF groups (p-value 2). Survival analysis included patients at-risk or with HF.
Those with HFpEF were older and had higher body mass index, higher medication use, more concomitant disease and higher serum NT-proBNP, as compared to healthy controls and those at-risk (all ps < 0.001, Table 1). Among the three HF groups, all had similar serum creatinine and concomitant disease, however those with HFpEF were slightly older with a larger proportion of females, slightly higher body mass index, and lower rate of beta blocker use (all ps < 0.05) and those with HFrEF had higher level of serum NT-proBNP with lower systolic blood pressure (both ps < 0.05) (Table 1).
Ventricular structure and function
Within the HF groups, there was significantly lower LVEF, lower strains and larger LVM and volumes from HFpEF to HFmrEF and to HFrEF groups, respectively (Table 2, all ps < 0.001). LVEF was similar within the preserved LVEF groups (controls, at-risk and HFpEF) by definition. All circumferential strain components, GLS_ENDO as well as LV and RV volumes were also similar in the preserved LVEF groups. GLS_EPI was the only strain parameter that distinguished all three groups in post-hoc analysis, with incrementally reduced systolic function from healthy controls, to those at-risk to HFpEF (− 16.5 ± 2.4% vs. − 15.5 ± 2.7% vs. − 14.1 ± 3.0%, p < 0.001). Other parameters, including GLS_AVE and GRS, also identified systolic dysfunction in the HFpEF group, with significantly reduced values versus healthy controls and at-risk groups (Table 2). Finally, HFpEF patients were also distinguished from healthy controls by significantly increased differences in endocardial and epicardial strains, measured as absolute or relative differences or a ratio (ENDO/EPI), for both GLS and GCS. Both LVM index and LVM/LVEDV also distinguished all three groups in post-hoc analysis, with a significant stepwise increase from controls, to those at-risk to HFpEF.
LVEF values were strongly overlapping in the three preserved LVEF groups and distinct in the three HF groups, by definition, as illustrated in Fig. 4a for all study subjects. GLS_EPI values, conversely, had overlapping values in all HF patients while also being reduced in the HFpEF group versus the other LVEF-matched groups (Fig. 4b). Using the 90th percentile in the control group as a cutoff for identification of systolic dysfunction, with consideration of sex differences, GLS_EPI identified 35/87 HFpEF patients (GLS_EPI > − 11.4% cutoff in males (Fig. 4d), GLS_EPI > -14.2% cutoff in females (Fig. 4f)), while GLS_ENDO identified a much smaller subgroup of 18/87 with reduced function (GLS_ENDO >− 16.0% cutoff in males, GLS_ENDO > − 18.7% cutoff in females) (not shown). Similar scatter plots including all study subjects are shown for all strain components (Fig. 5). All three long axis and two short axis slices were included in the strain analysis for all subjects.
GLS and GCS values were significantly reduced from ENDO to EPI layers in all groups. Additionally, the ratio of strains between the layers and the absolute and relative strain differences between the layers illustrate the relatively lower strains on the EPI (Table 2). Within the HF groups, the differences between ENDO and EPI strains (ratio, absolute and relative differences) were significantly reduced from HFpEF to HFmrEF to HFrEF (Table 2).
The significant linear relationship between strain and LV structure in the preserved LVEF groups is illustrated by comparison of GRS (GRS_SAX) and LVM/LVEDV (concentricity) values in HFpEF (top row) and all individuals with preserved LVEF (bottom row), with separate results by sex (Fig. 6). Similar significant linear correlations were observed when comparing GLS_EPI and LVM/LVEDV (not shown). However, all GCS components had no significant relationship with LVM/LVEDV (not shown).
Layer-specific strain in healthy controls
Among the 77 healthy controls subjects (29 male, median age 59 years), male subjects had lower strains than females at all locations except for GCS_ENDO which were similar (Table 3). For both GLS and GCS, values were lowest at the epicardium with a stepwise increase to the average value in myocardium and to the endocardium (all ps < 0.001). The magnitude of GRS values measured from long axis (GRS_LAX) and short axis slices (GRS_SAX) were similar (p > 0.05).
During a 5-year follow-up for the at-risk and HF subgroups, there were 33 events of all-cause mortality. The base model included age and log (NT-proBNP) for multivariable analysis. After adjustment with the base model, GLS_AVE and GRS (for both LAX and SAX-derived GRS) were the only independent predictors of mortality, with adjusted hazard ratio 1.10, 1.03 and 1.03, respectively, for 1% absolute decrease of strains, all ps < 0.05. LVEF and GLS_ENDO were unable to independently predict outcome (Table 4). Additionally, AIC values for the layer-specific GLSs, GCSs and conventional LV measurements were compared to evaluate the relative strength of their association with clinical outcomes. GLS_AVE and GRS had the lowest AIC values, which demonstrates the strongest association with the clinical outcome (Table 5). For the remaining parameters, GLS_ENDO, GLS_EPI, all GCS components and conventional LV measures (LVEF, indexed LVEDV (LVEDVI), indexed LVESV (LVESVI) and indexed LVM) had considerably less support for similar performance as GLS_AVE and GRS for outcomes prediction, with ΔAIC values of 3.5 to 6.1 . The use of GLS_AVE had significant incremental value for prediction of outcomes over LVEF and GLS_ENDO (Δχ2 = 5.1 and 7.4, respectively, both ps < 0.05), by the likelihood ratio test. GRS, from both at short axis and long axis images, had incremental value over LVEF (Δχ2 = 4.4 and 4.8, respectively, both ps < 0.05).
The intra-observer reproducibility of all layer-specific strains was excellent; CoV values ranged from 2.7 to 5.4%, with 4.3% for GLS_EPI and 2.7% for GLS_ENDO, and ICC values ranging from 0.97 to 0.99, with 0.98 for GLS_EPI and 0.99 for GLS_ENDO. The inter-observer reproducibility of all layer-specific strains was also excellent; CoV values ranged from 4.9 to 9%, with 4.9% for GLS_EPI and 8.7% for GLS_ENDO, and ICC ranged from 0.92 to 0.98, with 0.92 for GLS_EPI and 0.98 for GLS_ENDO. Comparison of the custom feature tracking software used in the current study and the commercially available CVI42 for all GLS and GCS components in 202 subjects illustrated good agreement between the methods. (Additional file 1: Fig. 1).
The main findings of the current study are that consideration of measurement layer for global strain is necessary for optimal identification of dysfunction and outcomes prediction in HF. The endocardium-specific strains were shown to have poorest performance both for detection of systolic dysfunction and outcomes prediction.
Layer-dependence of systolic dysfunction and association with remodeling in HFpEF
In subjects with preserved LVEF, GLS_EPI distinguished all three groups with preserved LVEF and identified ~ 40% of HFpEF patients as having reduced systolic function. In contrast, GCS components at all layers and GLS_ENDO performed similarly to LVEF (i.e. values were comparable in healthy controls, those with HF risk factors and HFpEF). Similar to recent studies in healthy subjects [10, 29, 30] and in HF , GLS and GCS values were incrementally decreased from ENDO to EPI in all groups. In the current study, it was also found that the difference between ENDO and EPI strains in individuals (relative and absolute differences) were significantly higher in HFpEF as compared to healthy controls. This is in agreement with the patterns of preserved ENDO strains and reduced GLS_EPI in this HF group.
Additionally, indexed LVM and LVM/LVEDV were significantly increased in the HFpEF group as compared to at-risk and healthy controls. Together, these findings describe a HFpEF phenotype with preserved ENDO function, paralleling LVEF, but with increased strain reduction across the wall, potentially associated with increased LVM and concentricity (mass/volume). GRS was also significantly reduced in the HFpEF group as compared to healthy controls, and was shown to be associated with increased concentricity. A similar relationship was observed for GLS_EPI, with more impaired function in those with increased concentricity. These observed associations between strains and structure also exist in comparison of the healthy men and women, where men had larger LVM, larger LVM/LVEDV, lower GLS_EPI, lower GRS and larger differences between ENDO and EPI strains.
In the current study, the use of a higher cutoff of 55% for preserved LVEF [21, 22], as compared to the more commonly used 50%, was used to reflect the larger number of HF patients in the 50%-55% range but very small proportion of at-risk or control subjects in this range. He et al. showed that mildly decreased LVEF, in the range of 40–55% (matching our HFmrEF group in the current study), is associated with eccentric remodeling and decreased chamber contractility, from invasive pressure–volume analysis, most comparable to subjects with overt systolic HF. The use of the lower 50% cutoff for preserved LVEF would have resulted in a significantly lower LVEF in our HFpEF group, compared to at-risk and healthy controls, and a large number of individuals in the HFpEF that are outside of the normal range of LVEF values from the control group and literature values.
For many of the reported structural and functional parameters that distinguish HFpEF from healthy controls, the values in the at-risk group were intermediate, between the control and HFpEF groups, with significant group differences on post-hoc analysis, suggesting early changes in these parameters may contribute to or be associated with future development of HF.
When comparing the three HF groups alone (preserved, mid-range and reduced LVEF), all CMR parameters were significantly different between all groups, with larger volumes, increased LV mass and reduced strain (all components) from HFpEF to HFrEF. The differences between endocardial and epicardial strains also followed this pattern, with reduced differences between the layers from HFpEF to HFrEF, paralleling the reduction in concentricity (mass/volume).
Heart failure outcomes
Among all strain components, only GLS_AVE and GRS were significantly predictive of mortality in HF patients and those at risk when including the key factors of age and NT-proBNP in the outcomes model. Both parameters had superior outcomes prediction performance as compared to commonly reported GLS_ENDO [6, 7, 11, 12, 31]. Not surprisingly, LVEF was also not predictive of outcomes, in agreement with large studies showing that HF outcomes are independent of LVEF, being similar for HFpEF, HFmrEF and HFrEF groups . Similar to our findings, a recent study of 463 patients with HFpEF showed no association between endocardial GLS and mortality or a composite of mortality or rehospitalization at 1 year . While reduced GLS_ENDO has previously been shown to have incremental value over LVEF in outcomes prediction in acute HF [7, 31], NT-proBNP, a robust and widely available predictor of mortality [33, 34], was not included in the reported outcomes models. Case in point, in the current study, GLS_ENDO, GLS_AVE and GLS_EPI and even LVEF were all significantly predictive of outcomes in univariable analysis, but of these only GLS_AVE remained independently associated with clinical outcomes when including NT-proBNP and age in the outcomes model. Similar to our findings, a recent CMR feature tracking study showed that GLS_AVE and GRS were predictive of outcomes in HF . Different from our study, GCS_AVE was also shown to be predictive of outcomes. This may reflect differences in the patient cohorts and the overall limited number of outcomes. However, the authors acknowledged that lack of NT-pro-BNP in their study, for multivariable analysis of risk, was a limitation given the very strong predictive power of NT-pro-BNP in heart failure.
The limited improvement in endocardial-specific strain over LVEF for detection of dysfunction or prognosis might be expected given the strong relationship between volumetric function, from which LVEF is defined, and deformations of the endocardial surface. Recently, Stokke et al. highlighted the direct relationship between endocardial strains and LVEF, and that LV wall thickness becomes a modulator of this relationship when considering average strains, measured across the wall thickness . Specifically, their mathematical model showed how thicker ventricles can have reduced average strains across the wall thickness in the presence of preserved LVEF. Similarly, MacIver et al. illustrated the mechanism by which mid-wall global strains can be reduced in the presence of normal LVEF as a result of increased wall thickness, based solely on geometry . Similar geometric models linking wall thickness, LVEF and strains, with similar conclusions, have previously been described . It is thus possible that the superior prognostic performance of GLS_AVE over endocardial strains reflects its dependence on both systolic strain and LV mass, where LV mass itself is predictive of cardiovascular outcomes [39, 40]. GRS had similar good outcomes prediction performance as compared to GLS_AVE. Like GLS_AVE, GRS directly incorporates geometric information from the full thickness of the myocardium. Additionally, reduced GRS was shown to be associated with increased concentricity (LV mass/LVEDV), so like GLS_AVE, may integrate structural remodeling and reduced function. GRS measured from short axis (GRS_SAX) or long axis (GRS_LAX) images had similar values in all groups and similar outcomes prediction performance. It is unclear why GLS_EPI was not significantly predictive of outcomes, given its sensitivity to systolic dysfunction in HF in the current study, and its association with structural remodeling, similar to GRS. It is possible that increased measurement variability at the epicardium has contributed to this finding. Similar to previous studies, reproducibility of ENDO strain was superior to those measured on the EPI .
It should be acknowledged that the conventional strain components, longitudinally or circumferentially, are not along the direction of fiber shortening. Fibers are helically orientated at the ENDO and EPI . It is possible that reduced strain along a single direction at the EPI, for example, could reflect changes in fiber orientations and function across the thickness of the heart wall .
Dependence on sex
Our study has a modest sample size and number of events for outcomes analysis, a reflection of the clinical stability of our HF cohort, which may limit the generalizability of the reported findings. However, based on the suggested 10 outcome events per predictor, our use of age, log (NT-proBNP) with each CMR parameter separately (i.e. 3 predictors in the composite model), our observed 33 events (mortality at 5 years) are sufficient. The base model, while only including age and NT-proBNP was established based on AIC forward selection including all cardiovascular disease risk factors and disease history with p-values < 0.2 and is thus robust and representative. Also, age and NT-proBNP have previously been reported to be the strongest predictors of mortality among clinical characteristics in patients with HF or without HF . An additional limitation is the generalizability of the reported findings given the heterogeneity of the underlying pathologies of HF, some of which may not follow the trends reported in the current study. For example, endocardial strains have been shown to have superior prognostic performance in ischemic disease [14, 15]. Nonetheless, the overall significant findings for GLS_EPI, GLS_AVE and GRS speak to the robustness of these metrics across a wide spectrum of functional abnormalities. The calculation of average strain across the wall is limited by the relatively poor contrast within the myocardium, which may lead to errors in this strain component. This is an intrinsic limitation of CMR feature tracking with conventional cine imaging, for which the myocardium has relatively uniform signal intensity. Two short axis slices at mid and mid-basal locations were used for assessment of circumferential and radial strains which would not be sensitive to apical wall motion abnormalities. The HFpEF group in the current study was slightly older than the healthy control and at-risk groups, which may contribute to the lower observed strains in this group, however, the relatively small expected decline in strains beyond 50 years of age  suggests these effects will be negligible. Age was also included as a co-factor in all statistical analyses to address potential age effects. CMR has lower availability than echocardiography for the measurement of strain. However, recent speckle tracking study of healthy subjects reported similar GLS layer-specific values to those reported in the controls in the current study [29, 30], suggesting that current echocardiographic methods are similar to CMR for layer specific stain evaluation. Also, direct comparison of strain measured with CMR feature tracking and speckle tracking echocardiography has shown good inter-technique agreement . Finally, the reported findings remain to be verified in an independent validation cohort.
Global strains measured on ENDO, EPI or average across the wall thickness are not equivalent for the identification of dysfunction or outcomes prediction in HF. The endocardium-specific strains were shown to have poorest all-around performance. GLS_AVE and GRS were the only CMR parameters to be significantly associated with 5-year all-cause mortality in multivariable analysis. GLS_EPI, GLS_AVE, GRS and the relative difference in endocardial and epicardial strains differentiated HFpEF patients from healthy controls, and increased LVM and LVM/LVEDV were generally associated with reduced strain in those with preserved LVEF.
Availability of data and materials
The data used in this study are available from the corresponding author on reasonable request.
American College of Cardiology
Angiotensin converting enzyme inhibitor
American Heart Association
Akaike information criterion
Angiotensin II receptor blocker
Body mass index
Body surface area
Balanced steady state free precession
Cardiovascular magnetic resonance
Global circumferential strain
Global longitudinal strain
Average global circumferential strain
Endocardial global circumferential strain
Epicardial global circumferential strain
Average global longitudinal strain
Endocardial global longitudinal strain
Epicardial global longitudinal strain
Global radial strain
Heart failure with midrange ejection fraction
Heart failure with preserved ejection fraction
Heart failure with reduced ejection fraction
Left ventricle/left ventricular
Left ventricular end-diastolic volume
Left ventricular end diastolic volume index
Left ventricular ejection fraction
Left ventricular end-systolic volume
Left ventricular end-systolic volume index
Left ventricular mass
N-terminal pro b-type natriuretic peptide
New York Heart Association Classification
Right ventricle/right ventricular
Right ventricular end-diastolic volume
Right ventricular end-diastolic volume index
Right ventricular end-systolic volume
Right ventricular end-systolic volume index
Relative wall thickness
Yancy CW, Jessup M, Bozkurt B, Butler J, Casey DE Jr, Drazner MH, et al. 2013 ACCF/AHA guideline for the management of heart failure: a report of the American College of Cardiology Foundation/American Heart Association Task Force on Practice Guidelines. J Am Coll Cardiol. 2013;62 (16):e147-239.
Ponikowski P, Voors AA, Anker SD, Bueno H, Cleland JG, Coats AJ, et al. ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure: The Task Force for the diagnosis and treatment of acute and chronic heart failure of the European Society of Cardiology (ESC). Developed with the special contribution of the Heart Failure Association (HFA) of the ESC. Eur J Heart Fail. 2016;18 (8):891–975.
Ezekowitz JA, O’Meara E, McDonald MA, Abrams H, Chan M, Ducharme A, et al. 2017 comprehensive update of the canadian cardiovascular society guidelines for the management of heart failure. Can J Cardiol. 2017;33 (11):1342–433.
Shah KS, Xu H, Matsouaka RA, Bhatt DL, Heidenreich PA, Hernandez AF, et al. Heart failure with preserved, borderline, and reduced ejection fraction: 5-year outcomes. J Am Coll Cardiol. 2017;70 (20):2476–86.
Aimo A, Januzzi JL Jr, Vergaro G, Petersen C, Pasanisi EM, Molinaro S, et al. Left ventricular ejection fraction for risk stratification in chronic systolic heart failure. Int J Cardiol. 2018;273:136–40.
Kraigher-Krainer E, Shah AM, Gupta DK, Santos A, Claggett B, Pieske B, et al. Impaired systolic function by strain imaging in heart failure with preserved ejection fraction. J Am Coll Cardiol. 2014;63 (5):447–56.
Park JJ, Park JB, Park JH, Cho GY. Global longitudinal strain to predict mortality in patients with acute heart failure. J Am Coll Cardiol. 2018;71 (18):1947–57.
Greenbaum RA, Ho SY, Gibson DG, Becker AE, Anderson RH. Left ventricular fibre architecture in man. Br Heart J. 1981;45 (3):248–63.
Shi J, Pan C, Kong D, Cheng L, Shu X. Left ventricular longitudinal and circumferential layer-specific myocardial strains and their determinants in healthy subjects. Echocardiography. 2016;33 (4):510–8.
Tanacli R, Hashemi D, Lapinskas T, Edelmann F, Gebker R, Pedrizzetti G, et al. Range variability in CMR feature tracking multilayer strain across different stages of heart failure. Sci Rep. 2019;9 (1):16478.
Shah AM, Claggett B, Sweitzer NK, Shah SJ, Anand IS, Liu L, et al. Prognostic importance of impaired systolic function in heart failure with preserved ejection fraction and the impact of spironolactone. Circulation. 2015;132 (5):402–14.
Iacoviello M, Puzzovivo A, Guida P, Forleo C, Monitillo F, Catanzaro R, et al. Independent role of left ventricular global longitudinal strain in predicting prognosis of chronic heart failure patients. Echocardiography. 2013;30 (7):803–11.
Sarvari SI, Haugaa KH, Zahid W, Bendz B, Aakhus S, Aaberge L, et al. Layer-specific quantification of myocardial deformation by strain echocardiography may reveal significant CAD in patients with non-ST-segment elevation acute coronary syndrome. JACC Cardiovasc Imaging. 2013;6 (5):535–44.
Scharrenbroich J, Hamada S, Keszei A, Schroder J, Napp A, Almalla M, et al. Use of two-dimensional speckle tracking echocardiography to predict cardiac events: Comparison of patients with acute myocardial infarction and chronic coronary artery disease. Clin Cardiol. 2018;41 (1):111–8.
Hamada S, Schroeder J, Hoffmann R, Altiok E, Keszei A, Almalla M, et al. Prediction of outcomes in patients with chronic ischemic cardiomyopathy by layer-specific strain echocardiography: a proof of concept. J Am Soc Echocardiogr. 2016;29 (5):412–20.
Skaarup KG, Iversen A, Jorgensen PG, Olsen FJ, Grove GL, Jensen JS, et al. Association between layer-specific global longitudinal strain and adverse outcomes following acute coronary syndrome. Eur Heart J Cardiovasc Imaging. 2018;19 (12):1334–42.
Lee WH, Liu YW, Yang LT, Tsai WC. Prognostic value of longitudinal strain of subepicardial myocardium in patients with hypertension. J Hypertens. 2016;34 (6):1195–200.
Caspar T, Fichot M, Ohana M, El Ghannudi S, Morel O, Ohlmann P. Late detection of left ventricular dysfunction using two-dimensional and three-dimensional speckle-tracking echocardiography in patients with history of nonsevere acute myocarditis. J Am Soc Echocardiogr. 2017;30 (8):756–62.
Kuznetsova T, Cauwenberghs N, Knez J, Yang WY, Herbots L, D’Hooge J, et al. Additive prognostic value of left ventricular systolic dysfunction in a population-based cohort. Circ Cardiovasc Imaging. 2016;9:7.
Ezekowitz JA, Becher H, Belenkie I, Clark AM, Duff HJ, Friedrich MG, et al. The Alberta Heart Failure Etiology and Analysis Research Team (HEART) study. BMC Cardiovasc Disord. 2014;14:91.
He KL, Burkhoff D, Leng WX, Liang ZR, Fan L, Wang J, et al. Comparison of ventricular structure and function in Chinese patients with heart failure and ejection fractions >55% versus 40% to 55% versus <40%. Am J Cardiol. 2009;103 (6):845–51.
Ueda T, Kawakami R, Nishida T, Onoue K, Soeda T, Okayama S, et al. Left ventricular ejection fraction (EF) of 55% as cutoff for late transition from heart failure (HF) with preserved EF to HF with mildly reduced EF. Circ J. 2015;79 (10):2209–15.
Robinson JD, Lupkiewicz SM, Palenik L, Lopez LM, Ariet M. Determination of ideal body weight for drug dosage calculations. Am J Hosp Pharm. 1983;40 (6):1016–9.
Klein S, Staring M, Murphy K, Viergever MA, Pluim JPW. elastix: a toolbox for intensity-based medical image registration. IEEE T Med Imaging. 2010;29 (1):196–205.
Pedrizzetti G, Claus P, Kilner PJ, Nagel E. Principles of cardiovascular magnetic resonance feature tracking and echocardiographic speckle tracking for informed clinical use. J Cardiovasc Magn Reson. 2016;18 (1):51.
Taylor RJ, Moody WE, Umar F, Edwards NC, Taylor TJ, Stegemann B, et al. Myocardial strain measurement with feature-tracking cardiovascular magnetic resonance: normal values. Eur Heart J Cardiovasc Imaging. 2015;16 (8):871–81.
Royston P, White IR. Multiple imputation by chained equations (MICE): implementation in stata. J Stat Softw. 2011;45 (4):1–20.
Burnham KP, Anderson DR. Multimodel inference - understanding AIC and BIC in model selection. Sociol Method Res. 2004;33 (2):261–304.
Alcidi GM, Esposito R, Evola V, Santoro C, Lembo M, Sorrentino R, et al. Normal reference values of multilayer longitudinal strain according to age decades in a healthy population: a single-centre experience. Eur Heart J Cardiovasc Imaging. 2018;19 (12):1390–6.
Nagata Y, Wu VC, Otsuji Y, Takeuchi M. Normal range of myocardial layer-specific strain using two-dimensional speckle tracking echocardiography. PLoS ONE. 2017;12 (6):e0180584.
Romano S, Judd RM, Kim RJ, Heitner JF, Shah DJ, Shenoy C, et al. Feature-tracking global longitudinal strain predicts mortality in patients with preserved ejection fraction: a multicenter study. JACC Cardiovasc Imaging. 2020;13 (4):940–7.
Buggey J, Alenezi F, Yoon HJ, Phelan M, DeVore AD, Khouri MG, et al. Left ventricular global longitudinal strain in patients with heart failure with preserved ejection fraction: outcomes following an acute heart failure hospitalization. ESC Heart Fail. 2017;4 (4):432–9.
Myhre PL, Vaduganathan M, Claggett BL, Anand IS, Sweitzer NK, Fang JC, et al. Association of natriuretic peptides with cardiovascular prognosis in heart failure with preserved ejection fraction: secondary analysis of the TOPCAT randomized clinical trial. JAMA Cardiol. 2018;3 (10):1000–5.
York MK, Gupta DK, Reynolds CF, Farber-Eger E, Wells QS, Bachmann KN, et al. B-type natriuretic peptide levels and mortality in patients with and without heart failure. J Am Coll Cardiol. 2018;71 (19):2079–88.
Sardana M, Konda P, Hashmath Z, Oldland G, Gaddam S, Miller R, et al. Usefulness of left ventricular strain by cardiac magnetic resonance feature-tracking to predict cardiovascular events in patients with and without heart failure. Am J Cardiol. 2019;123 (8):1301–8.
Stokke TM, Hasselberg NE, Smedsrud MK, Sarvari SI, Haugaa KH, Smiseth OA, et al. Geometry as a confounder when assessing ventricular systolic function: comparison between ejection fraction and strain. J Am Coll Cardiol. 2017;70 (8):942–54.
MacIver DH, Adeniran I, Zhang H. Left ventricular ejection fraction is determined by both global myocardial strain and wall thickness. Int J Cardiol Heart Vasc. 2015;7:113–8.
Aurigemma GP, Silver KH, Priest MA, Gaasch WH. Geometric changes allow normal ejection fraction despite depressed myocardial shortening in hypertensive left ventricular hypertrophy. J Am Coll Cardiol. 1995;26 (1):195–202.
Schillaci G, Verdecchia P, Porcellati C, Cuccurullo O, Cosco C, Perticone F. Continuous relation between left ventricular mass and cardiovascular risk in essential hypertension. Hypertension. 2000;35 (2):580–6.
Armstrong AC, Gidding S, Gjesdal O, Wu C, Bluemke DA, Lima JA. LV mass assessed by echocardiography and CMR, cardiovascular outcomes, and medical practice. JACC Cardiovasc Imaging. 2012;5 (8):837–48.
Kuetting DLR, Feisst A, Dabir D, Homsi R, Sprinkart AM, Luetkens J, et al. Comparison of magnetic resonance feature tracking with CSPAMM HARP for the assessment of global and regional layer specific strain. Int J Cardiol. 2017;244:340–6.
Maciver DH. The relative impact of circumferential and longitudinal shortening on left ventricular ejection fraction and stroke volume. Exp Clin Cardiol. 2012;17 (1):5–11.
Erley J, Genovese D, Tapaskar N, Alvi N, Rashedi N, Besser SA, et al. Echocardiography and cardiovascular magnetic resonance based evaluation of myocardial strain and relationship with late gadolinium enhancement. J Cardiovasc Magn Reson. 2019;21 (1):46.
The authors sincerely thank the patients and volunteers participating in our study. We gratefully acknowledge the program managers (Kelly Narine, PhD, Jennifer Beck BsAg, MBA, Lisa Tanguay, PhD, Beverly Armstrong, RN, MBA) and research coordinator (Marleen Irwin, RRT/C) for the Alberta HEART study.
Funding for this study was awarded by Alberta Innovates-Health Solutions (AI-HS, Grant #AHFMR ITG 200801018). Alberta Health Services and Alberta HEART investigators kindly contributed to this study. LX received a scholarship from the Sino-Li Ka Shing foundation. GYO receives salary awards from AI-HS, the Heart and Stroke Foundation of Canada and the Canadian Institutes of Health Research (CIHR). JRBD receives a salary awards from AI-HS and is a Canada Research Chair in Molecular Medicine. MJH is supported by the Moritz Chair in Geriatrics in the College of Nursing and Health Innovation at The University of Texas at Arlington.
Ethics approval and consent to participate
The study was approved by the University of Alberta and University of Calgary Health Research Ethics Boards and all study participants provided written informed consent.
Consent for publication
All authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Comparison of feature tracking software (custom software used in the current study versus CVI42) for calculation of layer-specific longitudinal and circumferential strains in a subset of 202 subjects (41 controls, 73 at-risk and 88 heart failure). See Table 2 for abbreviations.
About this article
Cite this article
Xu, L., Pagano, J.J., Haykowksy, M.J. et al. Layer-specific strain in patients with heart failure using cardiovascular magnetic resonance: not all layers are the same. J Cardiovasc Magn Reson 22, 81 (2020). https://doi.org/10.1186/s12968-020-00680-6