Abstract

Background. The ratio of early mitral inflow velocity to the global diastolic strain rate (E/E’sr) and global longitudinal systolic strain (GLS) of the left ventricle (LV) are emerging indices of diastolic and systolic functions, respectively, for the LV. Their prognostic significance in the prediction of mortality and cardiovascular (CV) outcomes remains underexplored in hemodialysis (HD) patients. Methods. This prospective study included 190 maintenance HD patients. The E/E’sr ratio and GLS were assessed using two-dimensional speckle tracking echocardiography. The clinical outcomes included overall mortality, CV mortality, and major adverse cardiovascular events (MACE). The associations between the E/E’sr ratio, GLS, and clinical outcomes were evaluated using multivariate Cox regression analysis. The incremental values of the E/E’sr ratio and GLS in outcome prediction were assessed by changes in Cox models. Results. Over a median follow-up period of 3.7 years, there were 35 overall deaths, 16 CV deaths, and 45 MACE. Impaired diastolic function with a higher E/E’sr ratio was associated with overall mortality (HR, 1.484; 95% CI, 1.201−1.834; ), CV mortality (HR, 1.584; 95% CI, 1.058–2.371; ), and MACE (HR, 1.205; 95% CI, 1.040−1.397; ) in multivariate adjusted Cox analysis. Worsening GLS was associated with overall mortality (HR, 1.276; 95% CI, 1.101−1.480; ), CV mortality (HR, 1.513; 95% CI, 1.088−2.104; ), and MACE (HR, 1.214; 95% CI, 1.103−1.337; ). The E/E’sr ratio and GLS had better outcome prediction than the E to early diastolic mitral annular velocity (E/E’) ratio and left ventricular ejection fraction (LVEF). Moreover, adding the E/E’sr ratio and GLS to Cox models containing relevant clinical and conventional echocardiographic parameters improved the prediction of overall mortality (), CV mortality (), and MACE (). Conclusion. The E/E’sr ratio and GLS, as emerging indices of LV diastolic and systolic functions, significantly predict mortality and CV outcomes and outperform conventional echocardiographic parameters in outcome prediction in HD patients.

1. Introduction

Cardiovascular (CV) disease is the leading cause of mortality in patients undergoing hemodialysis (HD) [1]. Higher prevalence of traditional risk factors and functional abnormalities of the heart may contribute to this high CV risk in end-stage renal disease (ESRD) [2, 3]. Pressure and volume overload could cause such cardiac abnormalities [4, 5]. Two-dimensional (2D) speckle tracking echocardiography (STE) allows for angle-independent quantification of myocardial deformation to more accurately reflect systolic and diastolic performances of all myocardial segments [6, 7].

2D STE can assess the left ventricular (LV) early global diastolic strain rate (E’sr). Furthermore, the early mitral inflow velocity (E) to E’sr ratio has been reported to be an emerging index of LV diastolic function [8, 9] and is strongly correlated with invasively measured LV filling pressure [1012]. The E/E’sr ratio is associated with unfavorable outcomes among patients with acute myocardial infarction [13] and systolic heart failure [14]. Global LV longitudinal systolic strain (GLS) has been recognized as a proper indicator of LV systolic function [7, 15]. Less negative GLS is associated with an increased risk of death in patients undergoing dialysis [16, 17]. However, the associations between the E/E’sr ratio and the risk of mortality and CV outcomes have never been investigated in chronic HD patients. Therefore, this study is aimed at examining the prognostic significance of the E/E’sr ratio and GLS in the prediction of overall mortality, CV mortality, and major adverse cardiovascular events (MACE) in maintenance HD patients. We further explored whether the emerging indices of LV diastolic and systolic functions outperformed the conventional echocardiographic parameters in the prediction of mortality and CV outcomes.

2. Materials and Methods

2.1. Study Patients

The inclusion criterion of the present study was patients with maintenance hemodialysis months at the outpatient HD unit. This study enrolled 219 maintenance HD patients at a regional hospital in Taiwan from March to October 2014. Patients with refusal of examinations (), lack of STE measurements (), and atrial fibrillation () were excluded. No patient was excluded because of poor echogenicity. Overall, 190 study patients were included (Figure 1). The study adhered to the Declaration of Helsinki and was approved by the Institutional Review Board of Kaohsiung Medical University Hospital, and all participants provided their written informed consent.

2.2. Echocardiographic Measurements

Patients received echocardiographic measurements in the left decubitus position by one well-experienced cardiologist, using a Vivid 7 system (GE Vingmed Ultrasound AS, Horten, Norway). The cardiologist was blind to patients’ clinical information. Early diastolic velocities (E’) of lateral and septal mitral annuli were averaged using Doppler tissue imaging to calculate the E/E’ ratio. The LV ejection fraction (LVEF) and LV mass were calculated using the biplane Simpson’s and Devereux’s methods, respectively [6]. The LV mass index (LVMI) was calculated as LV mass divided by the body surface area. Left atrial volume was calculated using the biplane area-length method. The left atrial volume index (LAVI) was defined as left atrial volume divided by the body surface area. Relative wall thickness was calculated by . All volumetric measurements and analyses were performed in accordance with EAE/ASE recommendations [6].

LV apical two- and four-chamber and long-axis views were obtained. The endocardial border was defined manually, and epicardial surface tracing was automatically performed by the system to create a region of interest [18]. The LV chamber was divided into six segments, with their strain and strain rate curves being analyzed. The peak segmental longitudinal systolic strain and early diastolic strain rates were determined from these curves (Figure 2). The E’sr and GLS were assessed and averaged in 18 LV segments from the three standard apical views (four-, two-, and three-chamber views). All the 18 speckle tracking segments were kept into the analysis in all patients. LV dimensions, LVEF, LAVI, LVMI, E’sr, and GLS were measured from the index beat [9, 19, 20]. A single beat was analyzed each time, and the values from three cardiac cycles were average to obtain each index. All STE data were recorded and analyzed offline using EchoPAC version 08.

2.3. Assessment of the Ankle-Brachial Index (ABI) and Brachial-Ankle Pulse Wave Velocity (baPWV)

ABI and baPWV were measured 10−30 minutes before the HD session using an ABI-form analyzer (Colin VP1000, Komaki, Japan) which simultaneously measured blood pressure in both arms and ankles. ABI was calculated as systolic blood pressure of the ankles divided by systolic blood pressure of the arms. The baPWV was automatically calculated as the transmission distance divided by the transmission time. The ABI is a simple and noninvasive test for establishing the diagnosis of peripheral artery disease (PAD) and a marker of generalized atherosclerosis, which are prevalent among HD patients and associated with worse clinical outcomes [21]. PAD was defined as an [22, 23].

2.4. Evaluation of Aortic Arch Calcification (AoAC) and the Cardiothoracic Ratio by Chest X-Rays

An experienced radiologist blind to the patients’ clinical information reviewed their chest X-rays and assessed AoAC using the scale proposed by Ogawa et al. [24]. The aortic arch was divided into 16 sections on the chest X-rays, and the number of sections with calcification was counted. The cardiothoracic ratio, assessed by the radiologist, was defined as the ratio of the transverse diameter of the cardiac shadow to the transverse diameter of the chest on the chest X-rays.

2.5. Demographic, Medical, and Laboratory Data

Demographic and medical data including age, gender, and comorbidities were obtained from patients’ medical records and interviews. Laboratory tests were conducted using overnight fasting blood samples obtained within 1 month of enrollment. Information of the use of medications, including angiotensin-converting enzyme (ACE) inhibitors, angiotensin II receptor blockers (ARBs), β-blockers, and statins, was obtained from medical records.

2.6. Definition of Overall Mortality, CV Mortality, and MACE

All study patients’ medical records and hospital courses were reviewed by two cardiologists to define the cause of death and MACE. CV mortality was defined as sudden cardiac death, fatal myocardial infarction, ventricular arrhythmia, fatal stroke, and heart failure. MACE was defined as follows: hospitalization for unstable angina, nonfatal myocardial infarction, sustained ventricular arrhythmia, hospitalization for congestive heart failure, transient ischemia attack or stroke, hospitalization for peripheral artery occlusive disease, and CV death [25]. If the patients experienced more than one MACE, only the first was analyzed. The model for MACE was censored at the development of MACE or the end of follow-up. All patients were followed until December 31, 2017, or the study endpoint (overall or CV mortality).

2.7. Statistical Analysis

All statistical analyses were carried out using SPSS version 19.0 (SPSS Inc., Chicago, IL, USA) for Windows. Data are expressed as percentages, deviation, or median (25th–75th percentile) for the dialysis vintage, triglycerides, high-sensitivity C-reactive protein (hs-CRP), and AoAC. The study patients were stratified into three groups according to tertiles of the E/E’sr ratio, with the 1st tertile as the reference category. Multiple comparisons among the study groups were performed using one-way analysis of variance followed by the post hoc Bonferroni test. These study patients were also classified into two groups according to a preserved GLS (≤−16%) or impaired GLS (>−16%) [26, 27]. Differences between groups were analyzed using the chi-square test for categorical variables and the independent -test for continuous variables. Cox proportional hazards analysis was performed to evaluate the associations between the E/E’sr ratio and GLS and development of overall mortality, CV mortality, or MACE. The adjusted covariates included age, sex, dialysis vintage, diabetes mellitus (DM), coronary artery disease, cerebrovascular disease, body mass index, systolic blood pressure, , baPWV, the cardiothoracic ratio, AoAC, albumin, triglycerides, total cholesterol, hemoglobin, creatinine, calcium-phosphorous product, hs-CRP, and the use of ACE inhibitors or ARBs, β-blockers, and statins. Survival curves for cumulative probability of overall mortality, CV mortality and MACE were illustrated using the Kaplan-Meier method and compared among mentioned groups of patients by the log-rank test. Direct comparisons between the E/E’ ratio and the E/E’sr ratio, as well as GLS and LVEF, were performed using multivariate models and assessed by changes in . Incremental model performance of the E/E’sr ratio and GLS over clinical and echocardiographic variables was assessed by comparing the model at each step. A value<0.05 was considered statistically significant.

3. Results

A total of 190 study patients were included. The mean age was years, and there were 98 males and 92 females. The mean values of the E/E’sr ratio and GLS were and , respectively. The patients were stratified into three groups according to tertiles of the E/E’sr ratio (<50.4 cm, 50.4–67.9 cm, and >67.9 cm). Patients belonging to the 3rd tertile of the E/E’sr ratio had higher prevalence of DM and coronary artery disease, higher cardiothoracic ratio, higher AoAC, lower serum creatinine, higher prevalence of ACE inhibitors or ARB and β-blocker use, higher LAVI, higher LVMI, lower LVEF, higher E/E’ ratio, higher E/A ratio, lower E’ wave, lower E’sr, and less negative GLS compared to patients in the 1st tertile of the E/E’sr ratio (Table 1).

Table 2 shows the comparison of baseline characteristics between patients with preserved GLS (≤−16%) or impaired GLS (>−16%). Compared to patients with preserved GLS, those with impaired GLS had higher prevalence of DM, coronary artery disease, and cerebrovascular disease; higher systolic blood pressure; higher prevalence of ; higher AoAC; higher hemoglobin; higher calcium-phosphorous product; higher LVMI; lower LVEF; lower E’ wave; lower E’sr; and higher E/E’sr ratio.

As shown in Figure 3(a), the cumulative incidence rates of overall mortality, CV mortality, and MACE were highest among patients in the 3rd tertile of the E/E’sr ratio, intermediate among those in the 2nd tertile and lowest among those in the 1st tertile ( for the trend). The incidence rates of overall mortality (), CV mortality (), and MACE () were higher among patients with impaired GLS in comparison with those with preserved GLS (Figure 3(b)).

3.1. Risk of Overall Mortality

Over a median follow-up of 3.7 years (interquartile range: 3.3-3.8 years), there were 35 (18.4%) deaths, including fatal CV events (), sepsis or septic shock (), gastrointestinal bleeding (), malignancy (), and liver failure () among study patients. The Kaplan-Meier curves for the cumulative probability of overall mortality according to tertiles of the E/E’sr ratio (Figure 4(a)) show substantially higher probability of overall mortality among patients in the 3rd tertile of the E/E’sr ratio compared to patients in the 1st or 2nd tertile of the E/E’sr ratio ( by the log-rank test). As shown in Figure 5(a), patients with impaired GLS had higher probability of overall mortality compared to those with preserved GLS ( by the log-rank test).

Table 3 displays the hazard ratios (HR) of the E/E’sr ratio and GLS for overall mortality with and without adjustment for demographic, clinical, biochemical factors. A high E/E’sr ratio (per 10 cm) was significantly associated with overall mortality in the unadjusted model (HR, 1.191; 95% confidence interval (CI), 1.090−1.301; ), in the age- and sex-adjusted model (HR, 1.185; 95% CI, 1.074−1.308; ), and in the multivariable model adjusted for age, sex, dialysis vintage, DM, coronary artery disease, cerebrovascular disease, BMI, systolic blood pressure, , baPWV, cardiothoracic radio, AoAC, albumin, triglycerides, total cholesterol, hemoglobin, creatinine, calcium-phosphorous product, and hs-CRP. This association holds significant (HR, 1.484; 95% CI, 1.201−1.834; ) after being further adjusted for medication use of ACE inhibitors or ARBs, β-blockers, and statins.

Furthermore, GLS (per 1%) was significantly associated with overall mortality in the unadjusted model (HR, 1.127; 95% CI, 1.043−1.219; ), in the age- and sex-adjusted model (HR, 1.139; 95% CI, 1.046−1.239; ), and in the full multivariable adjusted model (HR, 1.276; 95% CI, 1.101−1.480; ).

3.2. Risk of CV Mortality

Sixteen (8.4%) CV deaths were recorded during the follow-up period, including sudden cardiac death (), myocardial infarction (), ventricular arrhythmia (), fatal stroke (), and heart failure (). The Kaplan-Meier curves (Figure 4(b)) show higher cumulative probability of CV mortality among patients in the 3rd tertile of the E/E’sr ratio compared to patients in the 1st or 2nd tertile of the E/E’sr ratio ( by the log-rank test). Figure 5(b) shows that compared to patients with preserved GLS, those with impaired GLS had higher cumulative probability of CV mortality ( by the log-rank test).

As shown in Table 3, a high E/E’sr ratio (per 10 cm) was associated with CV mortality in the unadjusted model (HR, 1.271; 95% CI, 1.135−1.423; ), in the age- and sex-adjusted model (HR, 1.261; 95% CI, 1.112−1.429; ), and in the full multivariable adjusted model (HR, 1.584; 95% CI, 1.058−2.371; ). GLS (per 1%) was associated with CV mortality in the unadjusted model (HR, 1.202; 95% CI, 1.075−1.343; ), in the age- and sex-adjusted model (HR, 1.208; 95% CI, 1.072−1.360; ), and in the full multivariable adjusted model (HR, 1.513; 95% CI, 1.088−2.104; ).

3.3. Risk of MACE

Forty-five (23.7%) MACE were documented during the follow-up period, including hospitalization for heart failure (), coronary artery disease (), ventricular arrhythmia (), stroke (), peripheral artery disease (), and CV deaths (). As shown in Figures 4(c) and 5(c), patients in the 3rd tertile of the E/E’sr ratio had higher cumulative probability of MACE compared to patients in the 1st or 2nd tertile of the E/E’sr ratio ( by the log-rank test) and higher cumulative probability of MACE among patients with impaired GLS compared to those with preserved GLS over the follow-up period ( by the log-rank test).

Table 3 shows that the high E/E’sr ratio (per 10 cm) was associated with MACE in the unadjusted model (HR, 1.188; 95% CI, 1.101−1.282; ) and in the full multivariable adjusted model (HR, 1.205; 95% CI, 1.040−1.397; ). GLS (per 1%) was associated with MACE in the unadjusted model (HR, 1.174; 95% CI, 1.100−1.254; ) and in the full multivariable adjusted model (HR, 1.214; 95% CI, 1.103−1.337; ).

3.4. Comparison of the E/E’ Ratio and the E/E’sr Ratio to Overall Mortality, CV Mortality, and MACE

As shown in Table 4, the addition of the E/E’ ratio to the basic model (comprises age, sex, dialysis vintage, DM, coronary artery disease, cerebrovascular disease, BMI, systolic blood pressure, , baPWV, cardiothoracic ratio, AoAC, album, triglycerides, total cholesterol, hemoglobin, creatinine, calcium-phosphorous product, hs-CRP, and the use of ACE inhibitors or ARBs, β-blockers, and statins) did not significantly improve the prediction for overall mortality, CV mortality, and MACE. In contrast, the addition of the E/E’sr ratio to the basic model showed significant improvement of the prediction for overall mortality (, ), CV mortality (, ), and MACE (, ).

3.5. Comparison of LVEF and GLS to Overall Mortality, CV Mortality, and MACE

As shown by a direct comparison in Table 4, the addition of LVEF to the basic model did not significantly improve the outcome prediction. However, the addition of GLS to the basic model showed significant improvement of the prediction for overall mortality (, ), CV mortality (, ), and MACE (, ).

3.6. Incremental Value of the E/E’sr Ratio and GLS in Relation to Overall Mortality, CV Mortality, and MACE

The incremental values of the E/E’sr ratio and GLS in the prediction of overall mortality, CV mortality, and MACE are shown in Figures 6(a)6(c), respectively. The addition of the echo model (comprises LAVI, LVMI, LVEF, and the E/E’ ratio) to the basic model did not result in a significant improvement in the prediction of adverse outcomes. Moreover, the addition of the E/E’sr ratio and GLS to the basic model plus the echo model resulted in a further significant improvement in the prediction of overall mortality (), CV mortality (), and MACE ().

4. Discussion

In this study, we found that the higher E/E’sr ratio and GLS were independently associated with increased risk of overall mortality, CV mortality, and MACE in HD patients. The E/E’sr ratio was better than the E/E’ ratio, and GLS was better than LVEF in predicting adverse outcomes. Furthermore, the E/E’sr ratio and GLS had significant incremental prognostic values beyond clinical and conventional echocardiographic parameters.

An important finding of this study highlights that the E/E’sr ratio is a novel risk factor for overall mortality, CV mortality, and MACE in patients undergoing HD. At present, the E/E’ ratio is a recommended modality to assess LV diastolic function [8, 28] and associated with mortality in HD patients [29] but it still has some drawbacks such as angle dependency and risk of errors with . The E’sr obtained by STE from the whole left ventricle could overcome these limitations and more accurately represent global LV relaxation. Thus, the E/E’sr ratio correlates the LV filling pressures better than the E/E’ ratio does [10, 11]. The deformation-based E/E’sr ratio provides more important information with regard to global myocardial relaxation than the velocity-based E/E’ ratio, and the E/E’sr ratio was independently associated with adverse outcomes in several disease states and the general population [14, 3032]. Furthermore, we found the superiority of the E/E’sr ratio over the E/E’ ratio in predicting mortality and CV outcomes and the incremental prognostic value of the E/E’sr ratio and GLS over the conventional echocardiographic parameters in HD patients. The E’sr angle independently detects subtle myocardial motion, and it more precisely reflects LV global diastolic function compared with the E/E’ ratio [32].

Evaluating LV systolic function is fundamental on echocardiography [6], and LVEF remains the most widely utilized indicator. Technical limitations in the measurement of LVEF include suboptimal endocardial definition and the formulas that make assumptions with regard to the geometry of the left ventricle [33]. Furthermore, LVEF as a measure of contractility is affected by load dependency and LVEF is insensitive to identify the subtle degree of systolic dysfunction in patients with [34]. As a result, the association between LVEF and mortality was inconsistent in certain studies [34, 35]. Our study indicates the independent effect of GLS on overall mortality, CV mortality, and MACE, and GLS outperformed LVEF in predicting unfavorable outcomes in chronic HD patients. These findings are in line with previous studies on nondialyzed CKD and on ESRD patients [16, 17, 27, 36]. GLS can assess the function of longitudinally orientated myofibers, which are most vulnerable because of their subendocardial location. Although GLS is load dependent, GLS is sensitive to detect early subendocardial changes with better reproducibility than LVEF and reflect the extent of myocardial fibrosis and uremic cardiomyopathy, even in those with preserved LVEF [16, 17, 37].

Another important finding is that the addition of the E/E’sr ratio and GLS to models containing markers of atherosclerosis and conventional indicators for LV systolic and diastolic functions as risk factors of mortality in HD patients [2, 21, 38, 39] offered incremental value in the prediction of adverse outcomes. The speckle tracking imaging is based on frame-by-frame tracking of the displacement of speckles within the myocardium during the cardiac cycle and subsequent measurement of LV deformations. This technique makes it independent on imaging factors including reverberation artifacts and attenuation. Thus, the E/E’sr ratio and GLS may be more representative of global LV function. Therefore, we suggest that the E/E’sr ratio and GLS should be measured during echocardiographic examinations to provide important prognostic information for chronic HD patients.

There are several limitations in the present study. First, the number of study patients is relatively small and the observation period may be not long enough. Second, the echocardiographic parameters were measured from the index beat. This method has been proved to be as accurate as the time-consuming method of averaging echocardiographic parameters from multiple cardiac cycles [9, 19]. Third, 2D STE generates longitudinal, radial, and circumferential deformation measurements and LV twist [40]. However, only E’sr and GLS were measured and analyzed in this study. The comparisons between these parameters with radial and circumferential strains and LV twist in predicting outcomes are warranted in the future study.

5. Conclusion

The E/E’sr ratio and GLS, as emerging indices of LV diastolic and systolic functions obtained from 2D STE, are useful parameters and are superior to the E/E’ ratio and LVEF in the prediction of mortality and CV outcomes in maintenance HD patients and may offer an incremental value of prognostic significance over relevant clinical and conventional echocardiographic parameters.

Data Availability

The data supporting the findings of this study are available within the article or are available from the corresponding author upon reasonable request.

Disclosure

A part of this work has been presented as poster in the 2017 Annual General Meeting and Academic Lecture of the Taiwan Society of Internal Medicine.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

The research idea and study design were from J.C.H., H.M.S., and S.C.C.; data acquisition was performed by J.C.H., H.M.S., P.Y.W., J.J.L., W.H.L., S.C.C., and Y.W.C.; data analysis/interpretation was performed by J.C.H., H.M.S., J.J.L., and Y.W.C.; statistical analysis was performed by J.C.H., P.Y.W., and S.C.C.; supervision or mentorship was done by S.C.C., Y.L.H., J.M.C., and H.C.C. All authors contributed important intellectual content during manuscript drafting or revision and approved the final version of the manuscript.

Acknowledgments

The research presented in this article is supported by the grant from Kaohsiung Municipal Hsiao-Kang Hospital (grant number: kmhk-107-004), Kaohsiung Medical University, Kaohsiung, Taiwan.