Introduction

Based on former large-scale studies, cardiovascular (CV) mortality and morbidity are many times higher in chronic kidney disease (CKD) than in the general population, causing a public health issue worldwide [1, 2]. For patients of young ages (25- to 34-year-old) with end-stage kidney disease, annual mortality is increased 500- to 1000-fold and corresponds to that of the patients above the 80-year-old general population [3, 4]. In the case of GFR decline, the total mortality hazard ratio is 5.9, and the CV events ratio is 3.4 higher in ESKD than in GFR > 60 ml/min/1.73 m2 [5, 6]. All of these increased CV risks are caused by more common traditional and nontraditional risk factors in patients with CKD [7]. The remodeling of the myocardium and blood vessels could be one of these risk factors leading to CV events, heart failure, and progression to end-stage kidney disease (ESKD). Therefore, identifying these risk factors and high-risk patients is very important for interventional strategies and managing patients with CKD.

Echocardiography is a widely used and valuable noninvasive method for the determination of the left ventricular systolic, measuring systolic ejection fraction (LVEF), and diastolic function (LVDD), which has prognostic significance in ischemic heart disease, heart failure, and end-stage renal failure [8,9,10,11]. Tissue Doppler Imaging (TDI) echocardiography is another way to measure the rate of myocardial contractility and help refine diastolic dysfunction. In community-based epidemiological studies, the ratio of E (transmitral E wave velocity) to Ea (early diastolic mitral velocity) has been reported to be significantly associated with LV diastolic function and filling pressure [12].

CKD patients have cardiovascular disease more than twice as often as non-CKD patients. According to the USRDS Annual Data Report, nearly 40 percent of patients with stage 4–5 CKD carried a diagnosis of HF in 2015, according to the latest study showing 41–50 percent [13, 14]. Data regarding the prevalence of heart failure with preserved ejection fraction (HFpEF) in dialysis are scarce and are usually underdiagnosed. Contrary to heart failure with reduced ejection fraction (HFrEF), HFpEF is more prevalent in women and older patients. Over-65 patients account for more than half of all HF cases in the community [15]. In the highest age decile (≥ 90 years old), nearly all patients with HF have preserved EF. Studies focusing on defining the percentage of patients with HFpEF vs. HFrEF incidents or prevalent HF cohorts have demonstrated that approximately half of the patients with HF have HFpEF [16]. Due to a variety of factors, CKD and HFpEF are becoming more prevalent. Whether due to a common etiology or arising independently, CKD and HFpEF are often coincident in patients. Furthermore, the population with both problems is expanding. A study published a decade ago found that renal dysfunction is associated with worse outcomes and higher mortality in HFpEF patients. Despite the association between CKD and adverse outcomes, the interaction between CKD, clinical features, and cardiac structural and functional abnormalities in HFpEF has not been fully understood [17].

Objective

The literature contains limited data about the prognostic effect of TDI parameters on cardiovascular and renal outcomes in CKD. The study's goal was to determine the predictive value of TDI parameters in a homogeneous group of IgAN patients with CKD stages G–3.

Materials and methods

Patients

We monitored 79 patients with IgAN at the University of Pécs' Clinical Center's 2nd Department of Internal Medicine, Nephrology, and Diabetes for an average of 70 months between 2009 and 2018. The diagnosis of IgAN was confirmed by renal biopsy in all patients. The local ethics committee approved the study protocol, and all participants gave their written consent to its completion. (Approval number: 3170/2008).

At the start of patient enrollment, echocardiography was performed, and classic CV risk factors (hypertension, carbohydrate metabolism disorder, obesity, lipid abnormalities, smoking) and patient medication were also recorded. The criteria of the ATP III (Adult Treatment Panel III) were used to identify the metabolic syndrome. The obesity criteria were a BMI over 30 kg/m2. Renal function was calculated using the CKD-EPI formula (eGFR, ml/min, 1.73 m2). Patients with severe comorbidities (active cancer treatment, fever, and kidney transplant recipients) were excluded. Renal replacement therapy, end-stage renal disease (CKD-5), and a history of kidney transplantation were also exclusion criteria. A 24-h blood pressure monitor was used by Meditech ABPM devices to determine the patient’s 24-h average systolic and diastolic blood pressure, pulse pressure, and diurnal index. Additional CV examinations (ergometry, coronagraphy, etc.) were also performed based on the patient's complaints.

Patients were observed regularly, and follow-up exams were performed every 3–6 months (more frequently if necessary). During these visits, medical events that had occurred since the previous visit were discussed, as well as physical status and detailed laboratory tests. Blood pressure values were determined from the average of three measurements taken after 10 min of rest.

The study's primary composite endpoints included cardiovascular outcomes, overall mortality, coronary intervention (due to an acute coronary event or acute myocardial infarction), stroke, and renal outcomes, such as the development of ESKD (eGFR decreased below 15 ml/min/1.73 m2 or renal replacement therapy was started). Then, the CV (overall mortality, coronary intervention, stroke) and renal endpoints (eGFR decreased below 15 ml/min/1.73 m2 or renal replacement therapy was started) were individually examined as secondary endpoints.

Echocardiographic measurement

The Aloka SSD 1400 echocardiography equipment was used. Using 2D images of the length of the apical left ventricular segment and the area of the left ventricular short-axis muscle, the left ventricular mass (LVM) was computed (LVM = (5/6 area * length)). The Cornell criterion was used to determine LVMI, which was then indexed for height (in meters). The left ventricular ejection fraction (LVEF) was calculated by adding the diastolic and systolic volumes of the left ventricle using the unidirectional (Simpson method: EF = (Dvol-Svol) / Dvol * 100)). Based on traditional spectral Doppler measurements, mitral inflow and pulmonary venous flow were used to evaluate diastolic function. We also calculated the isovolumetric relaxation time (IVRT), the E wave deceleration time, and the E wave to A wave ratio (E/A ratio). RWT and/or LVMI abnormalities were used to define LVH. TDI was used to measure the early and late displacements of the lateral and septal basal wall fragments closest to the left ventricle (Ea and Aa) and calculate the average (Fig. 1). Then determine the E/Ea and Ea/Aa ratios. To exclude interindividual differences, two investigators (cardiologist specialists) examined all patients.

Fig. 1
figure 1

High and low early (Ea) and late diastolic mitral annular velocity (Aa) at the interventricular septal annulus (A, B) and lateral side of the mitral annulus (C, D) measured by tissue Doppler echocardiography

Statistical analysis

We divided our patients into two groups according to the Ea cutoff (13 cm/s) for statistical analysis. Unless otherwise specified, all values are the mean and standard deviation. Differences between the two groups were compared by the Student’s t test, the Mann–Whitney U test for continuous variables, and the χ2 test for categorical variables. A bivariate correlation method (Pearson’s correlation) assessed the relationship between two continuous variables. The factors influencing Ea were investigated using univariate and multivariate linear regression analysis. We also performed univariate and multivariate logistic regression analyses to examine the relationships between Ea and other covariates. Survival was assessed by the Mantel-Cox log-rank test. The effect of factors influencing survival was analyzed by Cox regression analysis. SPSS version 22.0 (Statistical Program for Social Sciences for Windows, SPSS Inc., Chicago, IL) was used to analyze the data, and a significance level of 0.05 was used for statistical analysis. Values of p < 0.05 were considered statistically significant.

Results

Based on our results, to determine the best cutoff values of Ea, we made a logistic regression analysis (ROC curve). 13 cm/s was the best cutoff to predict primary and secondary endpoints. The patients were divided into two groups based on this cutoff. The baseline characteristics showed significant differences in age, blood pressure, metabolic parameters (hypertension, dyslipidemia, BMI, carbohydrate metabolism disorder), eGFR, antihypertensive drugs (ACE/ARB, BB, CCB), and statin usage. The two groups had no significant difference in LV ejection fraction, LV end-diastolic diameter, hemoglobin, or lipid levels. However, the lower Ea group had significantly higher diastolic dysfunction, proteinuria, and uric acid levels (Table 1).

Table 1 Baseline characteristics of IgAN patients

There was a significant correlation between eGFR and diastolic function parameters: mitral inflow E and A wave ratio (E/A) (r = 0.466; p < 0.01), E wave deceleration time (EDT) (r = − 0.270; p < 0.01), and tissue Doppler image parameters: early diastolic mitral annular velocity (Ea) (r = 0.544; p < 0.01), late diastolic mitral annular velocity ratio (Ea/Aa) (r = 0.455; p < 0.01) (Fig. 2).

Fig. 2
figure 2

Correlations between diastolic function parameters, tissue Doppler parameters, and eGFR

Patients with decreased Ea (< 13 cm/s) had significantly more endpoint events (20/42 patients versus 3/37 patients, Chi-square: 11.449; p = 0.001 by the Mantel-Cox log-rank test) than patients with higher Ea (≥ 13 cm/s). Analyzing the endpoints separately (cardiovascular or renal), the renal endpoint was significant (Chi-square: 8.441; p = 0.004), but the cardiovascular endpoint was not significant (Chi-square: 3.506; p = 0.061) (Fig. 3). The independent factor that influenced Ea was only EDT by univariate and multivariate regression (Table 2), and the independent factors that influenced Ea were age and albuminuria by logistic regression (Table 3). Using the Cox regression model, the primary endpoint independent predictors of survival were gender, eGFR, diabetes, dyslipidemia, Ea, and E/Ea (Table 4). The secondary renal endpoint predictors were gender, eGFR, dyslipidemia, urine albuminuria, and Ea. The secondary CV endpoint predictors were diabetes, eGFR, E/A, E/Ea, Ea, and Aa.

Fig. 3
figure 3

Kaplan–Meier curves based on Ea in IgAN in primary combined (A), secondary renal (B), and secondary cardiovascular endpoints (C)

Table 2 Uni- and multivariate regression analysis influencing factors for Ea
Table 3 Uni- and multivariate logistic regression analysis influencing factors for Ea
Table 4 Cox regression analysis of primary combined and secondary renal and cardiovascular endpoints influencing parameters

Discussion

Our study examined the relationship of Doppler and tissue Doppler echocardiography parameters with renal function in a homogenous immunocomplex-mediated CKD population of IgAN patients, and we found a correlation between E/A, EDT Ea, Ea/A, and eGFR.

Based on our results, in IgAN, the increase in diastolic dysfunction is observed with decreasing renal function, which is best described by EDT and Ea as defined by TDI, similarly in CKD.

Cardiovascular alterations develop in early CKD stages (echocardiographic parameters). There may be a difference in the dynamics of this, which may stem from the etiology of the disease [18], but there is a small amount of data in the literature on this. Patients at high CV risk should be screened in the early stages of CKD. TDI parameters are sensitive early markers, supported by a close correlation with eGFR, even in relatively few cases.

Previous studies have shown that E/Ea, an estimate of LV filling pressure by Doppler echocardiography, is a predictor of all-cause mortality in patients with LV systolic dysfunction and after acute myocardial infarction [19, 20]. Another study in patients with ESRD also reported that an E/Ea ≥ 15 could predict an increase in LV filling pressure with a sensitivity of 82% and a specificity of 88% and was associated with an increased risk of mortality [21]. In addition to predicting all-cause mortality, a high E/Ea has been reported to provide additional prognostic value in patients with ESRD beyond traditional echocardiographic parameters [21, 22]. Chen et al. also found that a high E/Ea was associated with an increased risk of CV events in patients with CKD [23].

A previous study evaluated the association between E/Ea and the progression of renal dysfunction in patients with and without CKD and found a higher E/Ea in the patients with a more rapid decline in renal function [24]. Chen et al. also reported an independent association between a high E/Ea and an increased risk of starting dialysis in patients with CKD stages 3–4 [25]. This implies that a high E/Ea ratio may lead to high volume status, increase renal efferent pressure, and decrease renal blood flow, subsequently leading to a progressive decline in renal function [25, 26]. A higher preload status may also contribute to a more rapid progression to dialysis.

In the traditional pathophysiological model, pressure overload leads to concentric LV hypertrophy, fibrotic remodeling, and LVDD. Eventually, LVDD leads to left atrial (LA) hypertension and remodeling, pulmonary venous hypertension, and right ventricular and atrial remodeling.

In the other model, proinflammatory cardiovascular and noncardiovascular coexisting conditions lead to systemic microvascular endothelial inflammation, global cardiac and skeletal muscle inflammation, and subsequent fibrosis. Thus, systemic microvascular endothelial inflammation has been proposed as an additional mechanism leading to myocardial inflammation and fibrosis, increasing oxidative stress, and promoting alterations in cardiomyocyte signaling pathways. These changes promote cardiomyocyte remodeling as well as microvascular dysfunction in cardiac and skeletal muscle [27]. The main pathophysiological alteration leading to HFpEF remains incompletely defined.

In CKD, several conditions contribute to the pathogenesis of HFpEF, such as arterial hypertension [28]. One of the most significant cardiac alterations in CKD is left ventricular hypertrophy (LVH), and CKD contributes significantly to its development. It develops early in the progression of kidney dysfunction, is frequently accompanied by myocardial fibrosis and LVDD, and is an independent risk factor for mortality in this population. The role of CKD is well documented, but its basis is not fully understood. The effect of uremia on the myocardium includes structural changes such as cardiomyocyte hypertrophy, myocardial fibrosis, and thickening of the intramural arteries. Together, these structural changes predispose to LVDD in response to the cumulative action of traditional and CKD-related risk factors [29,30,31]. There is good evidence that interstitial fibrosis is related to changes in collagen myocardial metabolism. On the other hand, cardiomyocyte hypertrophy and vascular remodeling may be adaptive responses to pressure and volume overload [29]. Others, such as hyperphosphatemia, hyperparathyroidism, and hypovitaminosis D, play a greater role in more advanced stages of CKD and dialysis [32, 33]. Another important factor is renin–angiotensin–aldosterone system (RAAS) activation, potentially inducing myocardial fibrosis and hypertrophy. Activation of the intracardiac RAAS seems to be critically involved in the overload status observed in dialysis, but angiotensin II and aldosterone can also be involved in myocardial cell hypertrophy and fibrosis independent of afterload [34].

In Escoli et al.’s review, they suggest patients with CKD and ESRD should be monitored regularly (perhaps every 1–2 years) for the development and assessment of the severity of LVH and cardiac fibrosis, most likely with serial echocardiography [35]. Nevertheless, our data suggest that it may also be useful for CKD 1–3 stages.

It is well known that hypertension is a major risk factor for the development of LVDD in chronic hemodialyzed patients [36]. However, in the early stages of CKD, when the blood pressure elevation is not so significant, the mechanism of the relationship between CKD and LVDD is not fully characterized. Therefore, the value of the echocardiography examination should be important in CKD.

Nerpin et al. [37] identified a significant inverse relationship between eGFR and LVDD in a community-based elderly population. LVDD is common in CKD patients, even in the early stages [38]. Several echocardiographic parameters, including LVH, large left atrial volume, and decreased LV ejection fraction (LVEF), have been revealed to be associated with cardiac events among patients with CKD [39, 40]. Furthermore, CKD at this age was reported to be strongly correlated with adverse cardiovascular (CV) outcomes [40]. However, LVDD was not associated with incident heart failure (HF) or all-cause death in the CRIC study [41]. These discrepancies have arisen because the studies used different parameters to evaluate cardiac function, particularly LV diastolic dysfunction.

In Shu’s KNOW-CKD study, LVDD was independently associated with adverse CV outcomes and all-cause mortality in patients with predialysis CKD [42].

Liang et al. proved that systolic dysfunction and LVDD demonstrated mutually augmentative effects on CV mortality and suggested that, together with conventional nephroprotection, early cardioprotection should be emphasized for patients with CKD in the early stage. Therefore, cardioprotective management should be initiated as early as possible after CKD diagnosis [43].

Known CV risk factors such as baseline eGFR, proteinuria, and hypertension are also risk factors for CKD progression, contributing to the acceleration of renal function loss and progression to ESKD. However, the progression of CKD, which is a complex process, cannot be explained in all cases by these traditional risk factors.

In our study, we analyzed the prognostic role of LVDD for CV and renal endpoints. Based on the correlation between Ea and GFR in IgAN, decreased Ea and the development of LVDD may predict the progression of renal disease and CV events before reaching ESRD. In our IgAN patients, we found that the decreased Ea had a significant effect on both combined and individual renal and CV outcomes.

We hypothesize that a decreased Ea due to LVDD development and a decreased eGFR may synergistically affect the poor prognosis. Worsening renal function, a higher incidence of LVDD, and CV complications result in a worse prognosis and impaired LVDD (lower Ea). As a result of all this, a greater deterioration in renal function is expected. However, the clear connection between the deterioration of kidney function and LVDD is not elucidated.

LVDD is closely linked to CKD because the uremic milieu predisposes patients with CKD to systemic arterial stiffness and myocardial interstitial fibrosis, ultimately leading to LVH and impaired LV relaxation and compliance [44, 45].

More than half of the patients with stages 1 and 2 CKD had a normal diastolic function, whereas only 13% of those with stage 5 CKD had a normal diastolic function. This relationship between CKD progression and LVDD burden is consistent with that observed in a population study [46].

The association of LVDD with all-cause mortality remains controversial in CKD. The CRIC investigators revealed that LVDD was not associated with incident HF or all-cause mortality [42]. However, the CRIC study did not use tissue Doppler in the assessment of LVDD. Therefore, it is not surprising that no association was found between LVDD and outcomes in the CRIC analysis. In contrast, the Mayo Clinic group demonstrated a significant association between LVDD and all-cause mortality in patients with CKD; nonetheless, no interaction between LVDD and CKD stage has been identified [46, 47].

Hypertension is a very common complication in IgAN [47, 48], affecting 50–70% of patients, and this was confirmed by our data. Increased RAAS activity and hypertension in CKD also increase the incidence of vascular events; thus, RAAS blockade is the standard treatment (recommended in all guidelines) in these patients in general and for the patients who have IgAN [49,50,51]. Based on our former results and those of others, we thought that RAAS also plays a key role in the development of arterial stiffness and LVDD in renal disease, as in IgAN [49,50,51]. However, there is no data on whether ACEI and/or ARB treatment could afford an LVDD-lowering effect in patients with IgAN. In our study, at the start, more patients received ACEI and/or ARB therapy in the lower Ea group than in the higher Ea group. However, we were unable to distinguish between RAAS inhibitor users and non-users. At the end of the observation, almost all patients received RAAS blockers.

In our study, there was a significant difference in the use of a RAAS inhibitor between lower and higher Ea patients. Based on this observation, RAAS may be important in the evaluation of LVDD. However, it should be noted that the blood pressure of the study population was well controlled. Patients with IgAN exhibiting lower Ea had deteriorated renal function and an increased incidence of ESKD and CV complications in both sexes compared to those with higher Ea values. However, this may be particularly important because the further progression of CKD may be accelerated in older age and with impaired renal function, and more complications may develop with worse CV status, which may also worsen the prognosis by launching a "vicious circle".

Metabolic parameters play an important role in the progression of IgAN, as has been previously demonstrated by our workgroup and also by others [52, 53]. According to our findings, the appearance of metabolic syndrome in CKD is known to worsen the prognosis, and if the TDI parameters worsen, LVDD develops; these alterations cause a worsening in the metabolic parameters (and vice versa). It follows from all this that these parameters should be treated as soon as possible, which confirms the importance of complex metabolic risk reduction in these patients [52, 53].

Limitations of the study

Our results indicated that the Ea value obtained from echocardiography has prognostic significance; however, difficulties may occur during echocardiographic measurement. In some populations, specifically the elderly, a lack of cooperation can be a problem. Renal function was determined by estimating GFR, which is widely accepted in the literature, but the limitations of the formula are also well known. The extent and change of proteinuria were not examined in the present study. The evaluation of the results may also be weakened by the low number of cases, especially the low number of female patients. Our study follow-up time was not long enough to get significant results for the CV endpoints. We did not examine left or right atrial volume, myocardial strain, or strain rate, and we did not use this parameter for LVDD determination. In our study, there was no control group to compare our results with a non-CKD population. Despite these limitations, the results of this study highlight that the onset of target organ damage in CKD is predicted by decreased Ea.

Conclusion

Our results suggest that Ea as assessed by tissue Doppler echocardiography appears suitable for estimating prognosis in IgAN patients and that the lower value is an independent prognostic factor for ESKD and CV events, but there was no control group without IgA nephropathy, so it cannot be concluded that this was specific for this entity. Lower Ea should call attention to those CKD patients who have higher renal and CV risk at an earlier stage of CKD (G2-3) and need to be monitored more closely, referred for further CV tests, and given maximal nephroprotection.

Our findings support the role of echocardiography in the high-CV-risk population of CKD patients, which also helps to understand the relationship between heart abnormalities and renal impairment.

In conclusion, impaired renal function gradually correlates with LVDD and tissue Doppler parameters in patients with IgAN. Decreased renal function is associated with decreased Ea and LVDD, which are responsible for poor prognosis due to worse CV and renal outcomes. In the background, the role of common vascular and myocardial pathological remodeling, which is exacerbated by metabolic changes, could be hypothesized. To confirm our results, further large-scale, multicenter prospective studies are warranted to evaluate the role of CV risk factors in mediating the changes in the TDI parameter, Ea, as well as the complex relationships between CV disease and CKD.