Skip to main content

Layer-specific strain in patients with heart failure using cardiovascular magnetic resonance: not all layers are the same



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 [1]. 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 [6]. 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 [6] and predicting adverse outcomes in acute HF [7].

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) [8]. 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 [13]. 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 [16] and hypertension [17], and improved diagnostic performance for myocarditis [18]. EPI, mid-wall and endocardial GLS were shown to have similar predictive performance in population-based cohort [19]. Layer-specific strains were recently reported in a small HF cohort [10].

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.


Study population

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 [20]. 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) [20]. 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.

CMR protocol

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.

Image analysis

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 [23]. 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 [24]. 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.

Fig. 1
figure 1

Strain Method. a Registration of cine images is used to calculate pixel-wise displacement fields between each image frame and the reference frames. b The process of registration yields displacement images (∆xi, ∆yi) that contain the translation of each pixel from the reference frame to the target ith image frame. c The user specifies endocardial (red) and epicardial (green) points on the end-diastolic frame. A uniform grid of points is generated at end-diastole consisting of 20 equally spaced contours between the user-specified endocardial and epicardial tracings from, with 100 equally spaced points along each contour. d Each point, or equivalently, each series of points in a contour, is propagated to all other cardiac phases using the previously calculated displacement fields. Four sample contours are shown. Strain for each contour is calculated as the fraction change in the contour length to end-systole. For regions that contain multiple contours, as shown here for GLS_AVE and GRS, average strains are calculated as the average of the strains from the individual contours. GLS_AVE, average global longitudinal strain; GRS, global radial strain

Fig. 2
figure 2

Layer-specific contours in long axis cine images, at endocardial (red), epicardial (green) and equally spaced intramyocardial contours (blue) at end-diastole and end-systole used for the calculation of layer-specific strain. A subset of the 20 total intramyocardial contours used for calculation of average strain are shown. Sample tracings for a four-chamber view are shown for a health control (top) and a patient with heart failure (HF) with preserved ejection fraction (HFpEF) (bottom). Strain values at each layer over the full cardiac cycle are shown on the right

Fig. 3
figure 3

Layer-specific contours in short axis cine images, at endocardial (red), epicardial (green) and equally spaced intramyocardial contours (blue) at end-diastole and end-systole used for the calculation of layer-specific strain. A subset of the 20 total intramyocardial contours used for calculation of average strain are shown. Sample tracings for a mid-ventricular short-axis view are shown for a health control (top) and a patient with HFpEF (bottom). Strain values at each layer over the full cardiac cycle are shown on the right

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.

Statistical approach

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 [27].

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 [28]. 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.

Table 1 Baseline Characteristics in the 5 Subgroups of Subjects
Table 2 Cardiac structure and functions in the 5 subgroups of subjects

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.

Fig. 4
figure 4

Scatter and box plots of LV ejection fraction (LVEF) (a) and epicardial global longitudinal strain (GLS_EPI) (b) in all five groups. For each group, box plots show the median value, 25th and 75th percentiles and the full extent of the data. Outliers are identified as unfilled circles in the box plot. Scatter and box plots for LVEF (c and e) and GLS_EPI (d and f) in the preserved LVEF groups are shown with grouping by sex. Individuals with reduced GLS_EPI values are contained within the dashed boxes in d (men: 11% of controls, 17% of the at risk group and 32% of HFpEF are within the box) and f (women: 9% of healthy controls, 14% of the at risk group and 45% of HFpEF are within the box). All subjects in the preserved LVEF groups have normal LVEF by definition

Fig. 5
figure 5

Scatter and box plots for LVEF and all strain components for all five groups. Box plots show the median value, 25th and 75th percentiles and the full extent of the data. Outliers are identified as unfilled circles in the box plot

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).

Fig. 6
figure 6

Relationship between GRS and LV Mass / LVEDV (concentricity) in HFpEF patients (a and b) and in all subjects with preserved LVEF (c and d). See Table 2 for abbreviations

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).

Table 3 Comparison of layer-specific strain, LV mass and relative wall thickness between two sexes in healthy controls

Clinical outcomes

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 [28]. 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).

Table 4 Cox proportional hazard regression analysis for the outcome of 5-year all-cause mortality
Table 5 Discrimination performance of layer-specific Strains by Akaike information Criterion (AIC) for death at 5 years


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 [19], 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 [4]. 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 [32]. 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 [35]. 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 [36]. 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 [37]. Similar geometric models linking wall thickness, LVEF and strains, with similar conclusions, have previously been described [38]. 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 [41].

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 [8]. 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 [42].

Dependence on sex

Similar to previous studies [9, 30], men were shown to have lower absolute strain values than women for most strain components, necessitating the definition of sex-specific normal values.

Study limitations

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 [34]. 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 [29] 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 [43]. 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


Atrial fibrillation/flutter


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


Heart failure with midrange ejection fraction


Heart failure with preserved ejection fraction


Heart failure with reduced ejection fraction


Long axis


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


Short axis


Echo time


Repetition time


  1. 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.

    Article  PubMed  Google Scholar 

  2. 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.

    Article  PubMed  Google Scholar 

  3. 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.

    Article  PubMed  Google Scholar 

  4. 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.

    Article  PubMed  Google Scholar 

  5. 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.

    Article  PubMed  Google Scholar 

  6. 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.

    Article  PubMed  Google Scholar 

  7. 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.

    Article  PubMed  Google Scholar 

  8. Greenbaum RA, Ho SY, Gibson DG, Becker AE, Anderson RH. Left ventricular fibre architecture in man. Br Heart J. 1981;45 (3):248–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. 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.

    Article  CAS  PubMed  Google Scholar 

  10. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. 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.

    Article  PubMed  Google Scholar 

  13. 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.

    Article  PubMed  Google Scholar 

  14. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 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.

    Article  PubMed  Google Scholar 

  16. 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.

    Article  PubMed  Google Scholar 

  17. 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.

    Article  CAS  PubMed  Google Scholar 

  18. 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.

    Article  PubMed  Google Scholar 

  19. 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.

    Article  Google Scholar 

  20. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 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.

    Article  PubMed  Google Scholar 

  23. 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.

    CAS  PubMed  Google Scholar 

  24. 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.

    Article  Google Scholar 

  25. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  26. 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.

    Article  PubMed  Google Scholar 

  27. Royston P, White IR. Multiple imputation by chained equations (MICE): implementation in stata. J Stat Softw. 2011;45 (4):1–20.

    Article  Google Scholar 

  28. Burnham KP, Anderson DR. Multimodel inference - understanding AIC and BIC in model selection. Sociol Method Res. 2004;33 (2):261–304.

    Article  Google Scholar 

  29. 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.

    PubMed  Google Scholar 

  30. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. 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.

    Article  PubMed  Google Scholar 

  32. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  36. 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.

    Article  PubMed  Google Scholar 

  37. 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.

    PubMed  PubMed Central  Google Scholar 

  38. 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.

    Article  CAS  PubMed  Google Scholar 

  39. 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.

    Article  CAS  PubMed  Google Scholar 

  40. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 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.

    Article  PubMed  Google Scholar 

  42. 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.

    PubMed  PubMed Central  Google Scholar 

  43. 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.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


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.

Author information

Authors and Affiliations




MJH, JAE, GYO, JRBD, TA, IP and RBT contributed to the study design. LX was responsible for study design, data analysis, statistical analysis, and drafting the manuscript. JP contributed to data acquisition. IP and YM contributed to data analysis. RBT contributed to data acquisition and analysis, development of image analysis software and drafting the manuscript. All authors critically revised and approved the final manuscript.

Corresponding author

Correspondence to Richard B. Thompson.

Ethics declarations

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

Not applicable.

Competing interests

All authors declare no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1: Fig. 1

. 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.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: