Abstract
Background Heart failure (HF) and diabetes epidemiologically share a bidirectional relationship, with each increasing the incidence and worsening the prognosis of the other. Global longitudinal strain (GLS) measured by feature tracking cardiovascular magnetic resonance (CMR) has prognostic value in HF. However, this has not been established in HF patients with diabetes.
Objectives This study aimed to assess the prognostic value of GLS in HF patients with diabetes.
Methods Consecutive patients (n=315) with HF from six tertiary cardiac centres in Singapore underwent CMR at 3T, including GLS, late gadolinium enhancement (LGE), native T1 mapping, and extracellular volume fraction (ECV) mapping. The primary outcome was a composite of all-cause mortality or HF hospitalisation.
Results Compared to those without diabetes (n=156), the diabetes group (n=159) had a higher prevalence of LGE (76 vs 60%, p<0.05), higher T1 (1285±42 vs 1269±42ms, p<0.001) and higher ECV (30.5±3.5 vs 28.8±4.1%, p<0.001). Compared to those without diabetes, diabetes conferred a worse prognosis (log rank chi-squared 12, p<0.001, 74 events, median [interquartile range] 23 [18-24] months follow-up). In multivariable Cox regression models in the diabetes group, GLS was prognostic (hazard ratio 1.15 [95% confidence interval 1.06-1.25], p<0.001). When stratified by both diabetes status and median absolute GLS (9.9%), patients with diabetes and poorer GLS had the worst prognosis (log-rank chi-squared 31, p<0.001).
Conclusion Compared to heart failure patients without diabetes, those with diabetes have worse CMR markers of fibrosis and adverse prognosis. GLS is a powerful and independent prognostic marker in heart failure patients with diabetes.
Introduction
The diabetic heart is characterised by a number of structural abnormalities including diffuse interstitial myocardial fibrosis, myocyte hypertrophy, and impaired coronary microvascular perfusion, which have all been implicated in the development of both diastolic and systolic dysfunction [1, 2]. Patients with heart failure (HF) and diabetes have consistently worse clinical outcomes, including higher risk of hospitalisation for HF and death, compared to those without diabetes [3-5]. These findings appear to hold regardless of whether the HF is ischaemic or non-ischaemic in etiology, and regardless of left ventricular (LV) ejection fraction.
Cardiac magnetic resonance (CMR) imaging has become the non-invasive reference standard for evaluating the etiology of heart failure due to its ability to accurately assess cardiac morphology, function, and myocardial tissue characteristics. In particular, late gadolinium enhancement (LGE) permits visualisation of focal replacement myocardial fibrosis, while T1 mapping pre- and post-gadolinium contrast enables non-invasive measurement of the myocardial extracellular volume fraction (ECV), a quantitative marker of myocardial diffuse interstitial fibrosis. Additionally, global longitudinal strain (GLS), defined as the change in the LV myocardial length between diastole and systole divided by the original end-diastolic length, provides a measure of LV systolic function by CMR that is effectively the same as GLS measured by echocardiography [6, 7].
Patients with diabetes frequently have impaired GLS and have a higher degree of myocardial fibrosis as assessed by ECV and histology [8, 9]. GLS by CMR has been shown to have prognostic value in HF regardless of ejection fraction [10-12]. This has not been confirmed in the subgroup of patients with HF and diabetes. The current study therefore aimed to assess the prognostic significance of GLS by CMR in a cohort of patients with HF and diabetes across the spectrum of LV ejection fraction. We hypothesized that GLS would have an incremental prognostic association in this group, beyond LV ejection fraction and LGE.
Methods
Patients with HF were recruited prospectively across six tertiary cardiac centres in Singapore (Asian neTwork for Translational Research and Cardiovascular Trials [ATTRaCT], ClinicalTrials.gov NCT02791009). Patients were included if they presented to hospital with a primary diagnosis of HF, or if they attended a hospital clinic within 6 months of an episode of decompensated HF (requiring hospitalization or treatment in an out-patient setting). In all cases, a trained cardiologist adjudicated the clinical diagnosis of HF. The exclusion criteria were: HF primarily due to severe valve disease, HF due to acute coronary syndrome resulting in a transient episode of acute pulmonary oedema, severe renal failure (estimated glomerular filtration rate < 15 mL/min per 1.73m2), specific causes of HF (constrictive pericarditis, complex adult congenital heart disease, hypertrophic cardiomyopathy, eosinophilic myocarditis, cardiac amyloidosis, and acute chemotherapy-induced cardiomyopathy), isolated right heart failure, and life threatening non-cardiac co-morbidity with life expectancy of <1 year.
All patients underwent clinical assessment at baseline. Blood was drawn for biochemistry and circulating biochemical markers (N-terminal pro b-type natriuretic peptide [NT-proBNP] and high-sensitivity troponin T [hs-TnT]).
Ethics approval was obtained from the local Centralized Institutional Review Board in Singapore, and all participants provided written informed consent. The study was conducted in accordance with the principles of the Declaration of Helsinki.
Clinical Outcomes
The primary outcome was a composite of time to either first hospitalization for HF (regardless of prior history of hospitalization for HF) or all-cause mortality. Follow-up was conducted through a clinic visit at 6 months after baseline assessment and structured phone interviews with the participants at the 1- and 2-year timepoints. Data in patients who were lost to follow-up were censored at the date when the patient was last known to be alive and had not experienced an event.
Cardiovascular magnetic resonance image acquisition
All patients in the ATTRaCT cohort were invited and assessed for suitability for CMR. Those who agreed and were eligible underwent a standardized CMR protocol with a 3 Tesla MRI scanner (Ingenia, Philips Healthcare, Best, The Netherlands). Balanced steady-state free precision cines were acquired in the standard long-axis views and a short-axis stack from base to apex, as described previously [13]. LGE images were acquired at 10 min after 0.1 mmol/kg of gadobutrol (Gadovist®, Bayer Pharma AG, Germany) with a phase sensitive inversion-recovery gradient-echo imaging sequence. Typical parameters were: repetition time (TR) = 6.1 ms; echo time (TE) = 3 ms; time to inversion (TI) = 320–340 ms, flip angle 25°, voxel size = 1.5×1.7×8 mm, SENSE factor = 2.4, slice thickness 8mm with 2mm gap to match short-axis cine slice positions. The inversion time for optimal myocardial nulling was selected from an inversion time scout sequence. T1 maps were acquired at the basal and mid-ventricular short-axis levels, pre- and 15-min post-contrast with modified Look-Locker Inversion-recovery (MOLLI) 5s(3s)3s and 4s(1s)3s(1s)2s acquisition schemes, respectively [14].
CMR analysis
Image analysis was performed using CVI42 software (Circle Cardiovascular Imaging, Calgary, Canada) by trained imaging fellows at the National Heart Research Institute of Singapore CMR Core Laboratory, who were blinded to the clinical information of the patients. Ventricular volumes, mass and ejection fraction were measured from the short-axis cine stack, using manual contouring of the left ventricle in end-diastole and end-systole, excluding papillary muscles, as detailed previously [13]. LV volumes and mass data were indexed to body surface area. The presence of LGE was assessed qualitatively by two readers according to the recommendations by the Society of CMR [15]. Average native and post-contrast myocardial T1 values were measured by placing a region of interest (ROI) within the middle third of the short-axis myocardial wall at the basal- and mid-ventricular levels, while avoiding regions of focal LGE. The myocardium-blood pool interface was carefully avoided in order to prevent partial volume effects. Pre- and post-contrast blood T1 values were measured in a ROI drawn within the blood pool. ECV was calculated from the pre- and post-contrast average blood and myocardial T1 values, as described previously [16, 17]. Myocardial strain was analysed in the cine images using the Tissue Tracking Plugin [18].
Statistics
Normality was assessed for continuous variables using the Shapiro-Wilk test. Normally distributed data are presented as mean ± standard deviation. Non-normally distributed data are presented as median [interquartile range]. Comparisons were performed for continuous variables using the parametric Student t-test or the non-parametric Mann-Whitney U test. Categorical variables are presented as number (percentage) and compared using the χ2 test.
Univariable Cox regression analysis was performed to identify prognostic variables in the entire cohort and also in the diabetes subgroup. Clinically relevant variables (age, sex, body mass index, diabetes, smoking, coronary artery disease, hypertension, New York Heart Association functional class, systolic blood pressure, NT-proBNP and hs-TnT) and CMR markers of function, remodelling and fibrosis (LV end diastolic volume index, mass index, ejection fraction, GLS, LGE, Native T1, ECV) were tested in the Cox models. Covariates with a p-value <0.05 in univariable analyses were entered into the multivariable Cox model to identify independently prognostic variables, using forward stepwise selection. Event-free survival curves were examined using the Kaplan-Meier method and compared with the log-rank test. Statistical analyses were performed using SPSS Version 28 (Statistical Package for the Social Sciences, International Business Machines, Inc., Armonk, New York, USA) and GraphPad Prism 9.4.1 (GraphPad Software, Inc., San Diego, California, USA). A two-sided p-value <0.05 was considered as statistically significant.
Results
Figure 1 shows a flow chart describing patient inclusion. Of the 623 patients enrolled in the ATTRaCT study, 523 patients underwent a baseline CMR study. We excluded 168 studies performed without contrast. An additional 40 patients were excluded due to incomplete baseline clinical or CMR data sets, alternative diagnosis on the basis of CMR, or loss to follow-up. In total, 315 subjects (diabetes, n=159; without diabetes, n= 156) were included in the study cohort.
Abbreviations: HF: heart failure; ECG: 12-lead electrocardiogram; CMR: Cardiovascular magnetic resonance; HCM: Hypertrophic cardiomyopathy.
The baseline clinical characteristics are shown in Table 1. Compared to patients without diabetes, patients with diabetes were older (60±10 vs 56±12 years, p<0.001), more likely to have a history of hypertension (76 vs 52%, p<0.001), coronary artery disease (72 vs 52%, p<0.001) and had worse NYHA Functional Class (median Class II vs I, p=0.04). Patients with diabetes had higher creatinine (96 [81-118] vs 91 [77-105] μmol/L, p=0.02), elevated cardiac biomarkers of NT-proBNP (1091 [326-2272] vs 579 [232-1136] pg/mL, p<0.001) and hs-TnT (27 [14-41] versus 15 [9-22] ng/L, p<0.001). There was no difference in LV ejection fraction and GLS between the groups. The diabetes group had a higher prevalence of LGE (76 vs 60%, p=0.002, driven by higher rates of ischaemic LGE (51 vs 33%, p=0.003). Prevalence of non-ischaemic LGE did not differ between the groups. The diabetes group had higher native T1 (1285±42 vs 1269±42 ms, p<0.001) and ECV (30.5±3.5 vs 28.8±4.1%, p<0.001). Diabetes was independently associated with the presence of LGE and increased ECV, after adjustment for potential confounders (p<0.005 for both analyses).
Baseline clinical characteristics of the cohort according to presence or absence of diabetes.
Predictors of Primary Outcome in All Patients with Heart Failure
Over a median follow-up of 23 [18-24] months, 74 patients experienced the primary outcome (52 hospitalisations for HF, 22 deaths). In univariable Cox regression analyses for the entire cohort, clinical and plasma biomarkers associated with the primary outcome included diabetes, NYHA functional class, systolic blood pressure, NT-proBNP. CMR markers of adverse remodelling (LV mass index and end diastolic volume index), contractile function (LV ejection fraction and GLS) and myocardial fibrosis (presence of LGE, native T1 and ECV) were predictors of worse outcomes. Diabetes, NYHA functional class and GLS remained independently associated with outcomes in the multivariable analysis (Table 2 and Figure 2).
Univariable and multivariable Cox regression models (stepwise selection) in the entire cohort for the primary outcome of all-cause mortality or heart failure hospitalisation.
(A) Adverse prognosis in patients with heart failure and diabetes; (B) Patients with diabetes and GLS worse than cohort median (−9.9%) had worst prognosis. Patients with either diabetes or GLS worse than median had similar outcomes. Abbreviations: DM: diabetes mellitus; GLS: global longitudinal strain.
Predictors of Primary Outcome in Patients with Heart Failure and Diabetes
In the diabetes group, 50 patients experienced the primary outcome (35 hospitalisations for HF, 15 deaths). In univariable analyses, systolic blood pressure, NT-proBNP, LV ejection fraction, GLS and native T1 were associated with the primary outcome (Table 3). GLS remained associated with outcomes in the multivariable analysis (hazard ratio 1.15 [95% confidence interval 1.06-1.25], p<0.001). In other words, each percentage point of worsening in GLS was associated with a 15% increase in risk of hospitalisation for HF or death. Patients with diabetes and worse GLS (GLS > median value of -9.9%) had the worst prognosis (log-rank p<0.001, Figure 3). Of note, patients with diabetes and GLS ≤ -9.9% and patients without diabetes and GLS > -9.9% had similar outcomes (Figure 2B; p=0.70).
Univariable and multivariable Cox regression models (stepwise selection) in the diabetes group for the primary outcome of all-cause mortality or heart failure hospitalisation.
Top panel shows a patient with heart failure, diabetes and GLS ≤ -9.9%. Bottom panel shows a patient with heart failure, diabetes and GLS > -9.9%. Centre panel: event-free survival curves showing adverse prognosis in patients with heart failure, diabetes and GLS worse than cohort median (−9.9%). Abbreviations: DM: diabetes mellitus; GLS: global longitudinal strain.
Discussion
The main finding of this study is that GLS is a powerful independent predictor of adverse outcomes in patients with HF and diabetes, providing incremental prognostic information beyond several CMR markers of cardiac remodelling, function, and fibrosis. In the diabetes group, worse GLS (> median -9.9%) was associated with adverse prognosis. GLS is known to have prognostic value in HF, regardless of etiology or whether the EF is reduced or preserved [11, 12]. In this study, we have extended those findings to the diabetes subgroup, who are known to have a higher risk of adverse cardiovascular outcomes. Our findings are in agreement with a prior study which showed that speckle tracking echocardiography (STE) GLS has prognostic value in diabetes patients with dilated cardiomyopathy [19]. The results from the current study raise the possibility of GLS-guided risk stratification and management in patients with HF and diabetes. Indeed, there is emerging evidence of improvement in GLS with antidiabetic medications which may enable this in future [19-21].
GLS by CMR is well-suited for routine clinical use. The technique relies on software packages which track the endocardial and epicardial borders, and reference values are specific for the software used for analysis. STE has a number of limitations, including dependence on high quality 2D images, and results are often affected by foreshortening, dropout of the apical and anterolateral segments on apical views, and/or suboptimal acoustic windows [12, 22, 23]. GLS by CMR overcomes these limitations and can be measured from routinely acquired bright-blood, steady-state free precession long-axis cine imaging.
Disadvantages of GLS by CMR include its susceptibility to through-plane motion artefacts, limited temporal resolution for quantification of strain rate, and limitations in patients with contraindications to CMR [12]. Furthermore, GLS by CMR values may be affected by inter-vendor differences in algorithms, similar to STE, which has resulted in difficulties establishing reference values [23]. Nonetheless, GLS has shown close correlation with STE and has become an important component of the CMR examination alongside volumetric analysis and tissue characterisation in the assessment of HF [24]. A number of the limitations for GLS may be overcome by the mitral annular plane systolic excursion (MAPSE), a measure of left ventricular longitudinal function measured as the left-sided atrioventricular plane displacement. A recent study showed that artificial intelligence-derived MAPSE remained robust across vendors and demonstrated stronger associations with adverse outcomes in a large patient cohort compared to GLS [25]. Unfortunately, MAPSE was not measured in our study.
Diabetes was a strong independent predictor of the composite outcome of hospitalisation for HF or all-cause mortality in this HF cohort. This is in agreement with prior studies which have consistently shown worse cardiovascular outcomes in HF patients with diabetes, regardless of HF etiology and whether the ejection fraction is reduced or preserved [5, 26, 27]. Furthermore, Kaplan-Meier survival analysis showed similar event rates in the diabetes group with above median GLS (≤ -9.9%) compared to those without diabetes and below median GLS. These findings are similar to data from the Candesartan in Heart failure: Assessment of Reduction in Mortality and morbidity (CHARM) study, which demonstrated that patients with HF, preserved ejection fraction, and diabetes had a greater rate of HF hospitalisation than those with low ejection fraction and no diabetes [5]. These findings highlight the urgent need for therapeutic advances in patients with HF and concomitant diabetes.
In this study, HF patients with diabetes had worse CMR markers of myocardial injury, inflammation, and fibrosis. The association between diabetes and elevated ECV remained after adjusting for potential confounders that differed between the diabetic and non-diabetic groups including age, hypertension, coronary artery disease and presence of LGE. Furthermore, diabetes remained associated with the presence of LGE in logistic regression models, even accounting for age, sex, hypertension and coronary artery disease. These findings are in agreement with prior studies, including a recent meta-analysis, which showed an association between diabetes and a higher degree of myocardial fibrosis as assessed by histology as well as ECV by CMR [8].
Both focal and diffuse myocardial fibrotic processes are known to occur in patients with diabetes independently of co-morbid conditions, including coronary atherosclerosis and hypertension. Diffuse interstitial and perivascular fibrosis are structural hallmarks of diabetic cardiomyopathy, but focal replacement fibrosis can also be seen, even in the absence of coronary artery disease [28]. Several mechanisms may explain the fibrosis burden in diabetes. Hyperglycaemia is thought to upregulate the expression of profibrotic factors such as transforming growth factor beta 1 and down-regulate the activity of the matrix metalloproteinases [29]. Hyperglycaemia is also known to promote the formation of advanced glycation end products (AGEs) which can cross-link collagen in the interstitium, increasing their resistance to degradation. AGEs can result in generation of reactive oxygen species and oxidative stress which further promotes a pro-fibrotic state [8]. Pro-inflammatory cytokines and chemokines, as well as increased renin-angiotensin-aldosterone system activity in diabetes have also been implicated in the development of myocardial fibrosis.
Native T1 values were higher in the diabetes group and may reflect interstitial expansion due to myocardial fibrosis as well as myocardial oedema affecting the cellular and interstitial compartments [30]. There is some disagreement in the literature regarding whether native T1 values are increased in diabetes. Several studies have shown an association between diabetes and increased T1 [31-34]. Although a positive association was found between diabetes and native T1 in a recent meta-analysis, the result was not statistically significant [8]. The authors of that study suggested that the lack of statistical significance may have resulted from limited sample size of the included studies. Indeed, the diabetes cohort in the present study was larger than any of those in the included studies. Furthermore, T1 values are known to be dependent on a variety of factors, including field strength, pulse sequence, and region of measurement within the myocardium. In the present study, T1 maps were obtained at 3T using the same MOLLI sequence for all patients, and with experienced observers performing standardized analyses that may improve reproducibility and eliminate technical cofounders [35, 36]. Whilst it does appear that diabetes is associated with higher T1, larger studies controlling for the variability in T1 are required to more conclusively answer this question.
Limitations
One limitation of this study is the absence of measures of glycaemic control, which are known to be prognostic in patients with HF and diabetes [37]. However, STE GLS is known to have prognostic value independent of glycaemic control in diabetes cohorts with preserved and reduced ejection fraction, and therefore this is unlikely to have altered the prognostic associations for GLS by CMR [19, 38]. Our HF cohort included patients with both ischaemic and non-ischaemic etiologies, as well as both preserved and reduced ejection fraction. Unfortunately, the study was not powered for analyses of these subgroups. Documented history of coronary artery disease did not associate with outcomes in the univariable analysis and therefore etiology of HF is unlikely to have affected the results. Finally, this was a single centre study using feature tracking software from a single vendor. A larger, multicentre study would be required to determine whether these results are more broadly applicable across different analysis platforms.
Conclusion
Patients with HF and diabetes had worse CMR markers of injury, inflammation, and fibrosis, and an adverse prognosis. GLS is an important and independent prognostic marker in this group. Future studies should explore whether GLS-guided risk stratification and management can improve outcomes in this group of patients.
Data Availability
The data that support the findings of this study are included in this article or available from the corresponding author upon reasonable request.
Acknowledgements
The contribution of all the site investigators and clinical co-ordinators is acknowledged.