Skip to main content

Survival and analysis of predictors of mortality in patients undergoing replacement renal therapy: a 20-year cohort



optimal management of end-stage renal disease (ESRD) in hemodialysis (HD) patients should be more studied because it is a serious risk factor for mortality, being considered an unquestionable global priority.


we performed a retrospective cohort study from the Nephrology Service in Brazil evaluating the survival of patients with ESRD in HD during 20 years. Kaplan-Meier method with the Log-Rank and Cox’s proportional hazards model explored the association between survival time and demographic factors, quality of treatment and laboratory values.


Data from 422 patients were included. The mean survival time was 6.79 ± 0.37. The overall survival rates at first year was 82,3%. The survival time correlated significantly with clinical prognostic factors. Prognostic analyses with the Cox proportional hazards regression model and Kaplan-Meier survival curves further identified that leukocyte count (HR = 2.665, 95% CI: 1.39–5.12), serum iron (HR = 8.396, 95% CI: 2.02–34.96), serum calcium (HR = 4.102, 95% CI: 1.35–12.46) and serum protein (HR = 4.630, 95% CI: 2.07–10.34) as an independent risk factor for the prognosis of survival time, while patients with chronic obstructive pyelonephritis (HR = 0.085, 95% CI: 0.01–0.74), high ferritin values (HR = 0.392, 95% CI: 0.19–0.80), serum phosphorus (HR = 0.290, 95% CI: 0.19–0.61) and serum albumin (HR = 0.230, 95% CI: 0.10–0.54) were less risk to die.


survival remains low in the early years of ESRD treatment. The present study identified that elevated values of ferritin, serum calcium, phosphorus, albumin, leukocyte, serum protein and serum iron values as a useful prognostic factor for the survival time.

Peer Review reports


Chronic Kidney Disease (CKD) and its evolution to end-stage renal disease (ESRD) is increasingly recognized as a serious risk factor for mortality, being considered an unquestionable global public priority [1,2,3,4]. In 2016 CKD was the 13th in the list of causes of death, with projection that in 2040 it will be the 5th leading cause of death worldwide [2, 5].

Worldwide, the total number of individuals with acute kidney injury, CKD and Renal Replacement Therapy (RRT) exceeds 850 million, a figure that is double the estimated number of people with diabetes worldwide [6]. Data from the 2017 Annual Report of the European Renal Association - European Dialysis and Transplant Association (ERA-EDTA) [7] report that 83.311 individuals from all countries of Europe started RRT in 2016. At the end of the same year, the total number of individuals requiring RRT was 564.638, with more than 80% being on hemodialysis (HD). In Brazil, also in 2017, the total number of patients on chronic dialysis was equivalent to 162.583, with 91.8% of patients undergoing HD and 85% of this treatment modality funded by the Unified Health System (SUS) [8, 9].

HD is the major treatment modality for RRT worldwide, and it is also the one with the highest mortality rates, followed by kidney transplantation and peritoneal dialysis [4, 10, 11]. For these reasons, one of the most worrying outcomes of CKD is ESRD, in which there is a need for RRT.

Identifying the prognostic factors that are associated to and worsen ESRD certainly helps to reduce morbidity and mortality, especially in the first year of patients undergoing HD [4, 12]. For this, it is necessary to analyze the dynamics of biochemical parameter and different factors in the survival time of patients at baseline and for a long follow-up period, which few studies focus on doing [13]. Most epidemiological studies evaluate the prognostic value of a single complication to relate it to survival. However, considering them in isolation is inconsistent in practice, since patients have several irregularities concomitantly.

Thus, considering that the relationship of some biochemical parameter and baseline conditions in the first year of HD leads to lower short and long-term survival, the present study sought to investigate the survival time and the prognosis factors associated with this time in patients with ESRD undergoing HD.


Study design and population

This is a retrospective cohort study, which evaluated the data of 422 patients on HD, over a period of 20 years (January 1998 to December 2018).

All 463 patients registered at the Nephrology Service were investigated. All patients were followed retrospectively until death or censorship (end of follow-up, loss of follow-up, transfer or kidney transplant). Of the 463 patients, 22 were censored because received a transplant and 19 were censored because they were transferred or for loss of follow-up. Thereby, 422 patients had all the data available for analysis. The monitoring ended on December 31, 2018.

Data were obtained from medical records of patients undergoing treatment and from the computerized record system of the Nephrology Service in the municipality of Viçosa, Brazil, where the study was conducted and which works with a total of three sessions per week in HD with an average of 4 h each. This Nephrology Service, located in a hospital, is the only one located in the region and covers the entire Health Microregion of Viçosa, which has nine neighboring municipalities, totalizing 138.336 inhabitants.

For this study, we included consecutive patients with ESRD at any time during the study follow-up, with a minimum treatment time of one month, age ≥ 18 years and with all available baseline data. To be diagnosed with CKD, patients would have to show an irreversible decline in kidney function, that is, a Glomerular Filtration Rate (GFR) ≤ 20 ml/min/1.73m2 estimated by the CKD-EPI formula for more than a month. The diagnosis of underlying kidney disease was based on clinical, laboratory and radiological characteristics.

Demographic and clinical data of the study patients

The main baseline data set collected (one month after starting hemodialysis) were: sex, age, race, occupation, marital status, length of treatment, hemodialysis characteristics and the main cause of ESRD.

After the initiation of therapy, information on biochemical investigation was collected, including: pre and post-dialysis urea (mg/dl), serum creatinine (mg/dl), serum calcium (mg/dl), serum phosphorus (meq/lit), glucose (mg/dl), serum potassium (meq/dl), hemoglobin (g/dl), hematocrit (%), erythrocyte (g/dl), leukocyte (μg/l), serum albumin (g/dl), transferrin saturation (%), ferritin (μg/l), iron serum (μg/l), alkaline phosphatase (U/l), serum protein (g/dl), globulin (g/dl), parathyroid hormone (PTH) (pg/ml), aluminum (μg/l), cholesterol (mg/dl), triglycerides (mg/dl), calcium x phosphorus ratio (mg/dl), PRU (rheumatological profile) and nPCR (normalized protein catabolism rate) (g/dl). The analyzes were made at the baseline, that is, when the patients started treatment.

The cutoff points for each biochemical parameter were defined in terciles or according to the recommendations for patients undergoing HD [14,15,16].

Details of hemodialysis

All patients were dialyzed using standard bicarbonate hemodialysis, performed three times a week with an average duration of 4 h. Dialysis machines of individual proportion were used with water treated by reverse osmosis. Volumetric ultrafiltration control was available on all machines. The standard dialysate flow rate was 500 ml/min and blood flow rates were directed according to the patient’s needs. Dialyzer reuse was performed uniformly using automated methods.

The adequacy of dialysis was determined using the unique Kt/V pool (spKt/V) obtained by the Fresenius 4008S HD machine (Fresenius Medical Care AG, Bad Homburg, Germany) equipped with an online clearance monitor, in the same day of blood collection. Based on the guidelines of the Kidney Disease Outcome Quality Initiative (KDOQI) 2015 [16] in HD, a Kt/V ≥ 1.2 is considered as the minimum dialysis dose for three times a week. Therefore, the Kt/V cutoff point of 1.2 was used to study the effect of varying doses of dialysis on survival for this study.

Statistical analysis

The percentage, mean ± standard deviation (SD), or median and interquartile range (IQR) were used to describe and summarize the baseline data. We tabulated characteristics of patients at the time of their first record of the HD session. The primary objective of the analysis was death.

The initial event (E0) was characterized by the date on which the patient was admitted for HD treatment; the second event (E1), the date on which he stopped performing the procedure for any of the outcomes described above. The outcome of all patients in terms of survival and mortality, as well as the prognostic factors that had the greatest impact on patient survival were analyzed. The strength of association was measured by the Hazard Ratio (HR).

The Kaplan-Meier method was used to calculate the cumulative survival rate and the Log-Rank test, to extract factors that impacted the survival rate and examine the data. Log-Rank test is most effective at the beginning and in the middle of following the survival curve. In addition, Cox’s proportional hazards model was used with HD treatment as a time-dependent exposure to calculate risk rate and the influence of several clinical and demographic variables on patient survival.

The selection of variables was based on scientific evidence, clinical importance and the result of univariate analysis, so that variables with a p-value ≤0.200 were determined to be significant and included in the multivariate regression model. To find the optimal survival regression model, we used backward selection procedure by Wald’s test, which starts from a complete model and removing the weakest predictors from the candidate list one by one until only statistically significant (p ≤ 0.05) predictors remain (the selection cut-off value was from default as well as the importance of clinical concern) [17, 18]. The confidence interval (CI) adopted was 95% and a p < 0.05 was used for statistical significance.

The proportional hazards assumptions were tested using time-dependent covariate analysis in the adjusted models in SPSS and Kaplan Meier curves and have not been violated.

To avoid the problem of overfitting owing to the number of outcomes, we performed bootstrapping validation, in order to determine the CIs for estimating β in the Cox proportional hazard regression (1000 bootstraps) [19, 20]. The significance of the differences among different methods was determined with the use of the signed-rank test for bias and the bootstrap method for the interquartile range from the 1000 bootstrap samples. All analyzes were performed using IBM SPSS Statistics 23.0.


Beta coefficients of the independent variables obtained by multivariate Cox proportional hazards analysis were assessed using 1000 bootstrap resamples. No significant difference between original beta and bootstrapped coefficients was observed.

Table 1 shows some characteristics of the patient’s baseline and the laboratory values ​​measured. The total number of patients was 422. The average age at the beginning of HD was 64.02 ± 15.21 years (mean ± standard deviation), with 132 (31.3%) patients started HD before age 60 and 223 (52.7%) are men. Hypertensive renal disease with renal failure was responsible for 33.4% of the cases with varying etiology. The recommended Kt/v values ​​ (≥ 1.2) were found in 237 (56.2%) patients and the average GFR, calculated using the CKD-EPI formula, was 10.13 mL/min/1.73 m2.

Table 1 Baseline characteristics and laboratory values of patients on hemodialysis

Regarding biochemical parameter, the averages of serum calcium, hemoglobin, hematocrit, serum albumin and pre-dialysis urea were below recommendations for these patients, while PTH and triglycerides were above recommended and the other exams, within the normal range [14,15,16].

Figure 1 shows cumulative survival rate for all patients. The number of patients who died was 254 (60.2%); therefore, 39.8% of the survival time was censored. The mean survival time is 6.79 ± 0.37 years (mean ± standard deviation) and 95% CI 6.06–7.51. The survival rate was 82.3% in the first year, 49.1% in 5 years, 22.5% in 10 years and 13.3% at the end of the follow-up. The median survival time was 4.92 (3.914–5.924) (median ± interquartile range); this means that 50% of individuals have survived for at least 4.92 years.

Fig. 1

- Kaplan–Meier survival curves of chronic hemodialysis

Cumulative survival rate calculated by Kaplan Meier method is displayed according to each factor, and the Log-rank tests were used to compare these factors through univariate analysis (Tables 2 and 3).

Table 2 Influence of sociodemographic factors and baseline dialysis on the survival time
Table 3 Influence of biochemical parameter on the survival time of patients on hemodialysis

The univariate analysis showed, between the socio-demographic variables, there is a significant difference (p < 0.20) in age (< 60 years and ≥ 60 years), marital status and primary cause of ESRD. On the other hand, the analysis didn’t show statistical difference between men and women.

In the biochemical parameter, the variables with significant differences were: pre-dialysis urea, post-dialysis urea, serum creatinine, serum phosphorus, serum albumin, ferritin, serum iron, serum protein, PRU, nPCR, calcium x phosphorus ratio, serum calcium and leukocytes (Table 3).

As can be seen in Table 4, in this multivariate analysis it is considered that, among the clinical and demographic variables, those that significantly influenced patient survival were: the primary cause of ESRD, leukocytes, ferritin, serum iron, serum albumin, serum protein, serum calcium and serum phosphorus.

Table 4 Factors affecting survival time

Patients who presented chronic obstructive pyelonephritis as a cause of ESRD resulted in better prognosis of survival time when compared to glomerular disorder in diabetes mellitus (HR = 0.085, 95% CI: 0.010–0.744). Patients with other causes of the disease, such as urinary tract disorders, hypertensive nephrosclerosis and non-specific chronic renal failure, resulted in worst prognosis of survival when compared to glomerular disorder in diabetes mellitus (HR = 3.599, 95% CI 1.577–8.166) (Table 4 and Fig. 2a).

Fig. 2

Cox regression and the risk factors; ESRD, leukocyte, serum protein, serum calcium

Prognostic factors of survival were explored suggesting that patients with adequate values ​​of ferritin (178.4 μg/l ≥ 494.0 μg/l) and serum albumin (≥ 4.01 g/dl), resulted in better prognosis, with respectively, 2.85 (HR = 0.392, 95% CI 0.192–0.799) and 4.34 (HR = 0.230, 95% CI 0.097–0.541) compared to values ​​below the recommendation for HD patients (Table 4 and Fig. 3a and c).

Fig. 3

Cox regression and risk factors; ferritin, serum iron, serum albumin, serum phosphorus

In contrast, patients with elevated serum iron values ​​(≥ 150.0 μg / l), serum protein (≥ 7.01 mg/dl) and serum calcium (≥ 11.01 mg/dl) and normal values ​​of leukocytes (5.81 μg/l to 7.70 μg/l), have, respectively, a low survival prognosis: 8.396 (95% CI 2.016–10.341), 4.630 (95% CI 2.073–10.341), 4.102 (95% CI 1.351–12.457) and 2665 (95% CI 1.387–5.119) compared to patients with values ​​below the recommended (Table 4; Fig. 2b, c and d and Fig. 3b). Patients with serum phosphorus above recommended (≥ 6.01 meg / lit) have an HR of 0.290 (95% CI 0.138–0.609) (Table 4; Fig. 3d).


In both scientific and clinical practice, prediction aims at accurately predicting the risk of an outcome using multiple predictors collectively, where the final prediction model is usually based on statistically significant, but not necessarily causal, associations in the data at hand [21]. Choice in a dataset have impact of the predictimand of patients with ESRD, so, defining the estimand is equally important in prediction research as in causal inference [22].

Our study provided a set of indicators and parameters that showed to be was an important prognostic factor for survival of patients in the first year of treatment in HD, such as: basic cause of ESRD, serum phosphorus, serum calcium, leukocytes, serum protein, serum albumin, serum iron and ferritin.

The median survival time was 6.79 ± 0.37 years (mean ± standard deviation), and 254 (60.2%) patients died. In the study by Nyuyen et al. (2017) [23] the survival time was 5.27 ± 0.31 years (mean ± standard deviation), in the study by Browne et al. (2014) [1] the death rate was 60%, in the study by Ebhahimi et al. (2019) [4] of 52% and in the study by Chandrashekar et al. (2014) [24] 19.8% over 2 years. These periods are relatively short when compared to US rapporteurs [24, 25]. Although our study has similar values ​​to cohorts performed in other regions, there is variability in the survival rate of HD patients in different countries, therefore, these numbers cannot be compared to simple terms.

Regarding survival rate, it was 82.3% in the first year, 49.1% in 5 years, 22.5% in 10 years and 13.3% at the end of the follow-up. Nyuyen (2017) [23] found in her cohort an 85% survival rate at 1 year, 58% at 5 years and 20% at 10 years. Belino et al. (2017) [26] found a rate of 93.4% at the end of the second year of treatment and Teixeira et al., (2015) [27], found a survival rate in the first similar to our study: 84.71%. In the study by Chandrashekar (2014) [24], of the 19.8% death rate, a significant number of these deaths occurred in the first six months of treatment.

Comparing survival with that of normal population, adjusting mortality from all causes, is 6.1 to 7.8 times higher in HD patients than normal individuals, especially in the first year of treatment [28]. The high mortality rate in the first year of treatment can be justified by the lack of vascular access for hemodialysis, the age of patients when starting treatment and the lack of early diagnosis of CKD, when the first risk factors start to appear [26, 28, 29]. In addition, unplanned start of hemodialysis requires the use of catheters as a vascular access, which is proven to be one of the causes of the lowest survival rate in patients due to the urgency that this procedure requires, making it impossible for the patient to adapt adequately [30].

Currently, it is estimated that 425 million people worldwide are diagnosed with diabetes mellitus and it is estimated that more than 20 and up to 40% of diabetic patients develop CKD, with a significant number of people who develop CKD [30,31,32,33]. It is currently the leading cause of ESRD in developed countries [34]. In our study, individuals with glomerular disorder in diabetes mellitus had lower survival (p ≤ 0.001), which calls our attention to the monitoring of this disease as a prognostic factor for CKD development and subsequently ESRD.

In the analysis of the patient’s blood count, only the leukocyte was an independent prognostic factor for patient survival and its normal to high values (5.81–7.70 and ≥ 7.71) had an HR of, respectively 2.665 (95% CI 1.387–5.119) and 1.435 (95% CI 0.687–2.997). The study by Ebrahimi et al. (2019) [4], confirmed our findings by indicating that as the level of leukocyte count increases, the survival time of patients undergoing HD decreases. A study carried out in Taiwan, showed that the total leukocyte count predicts all causes and cardiovascular mortality in just one year [35], while in another study carried out in the United States, the reduction in the lymphocyte count are risk factors for mortality [36].

These findings can be justified by the fact that patients with ESRD on HD exhibit a chronic inflammatory state with a consequent increase in oxidative stress and impaired immune response [37,38,39]. The dialysis procedure itself contributes strongly to these adverse conditions due to exposure to synthetic materials contained in the dialysis filters and tubes, which favors a repetitive contact between circulating leukocytes and dialysis membranes and, clinically, this repetitive activation of cells can contribute to the increase in morbidity and mortality associated with HD [40, 41].

This predominant inflammatory condition of HD patients also explains the results we found in our study in relation to serum protein (95% CI 2.073–10.341). Elevated protein catabolism is common in patients with ESRD and its underlying etiology includes other complications, such as metabolic acidosis, uremia, systemic inflammation, anemia and other factors that further worsen kidney function and provide greater risks of mortality [42, 43].

Still on proteins, serum albumin is an independent and powerful prognostic index for HD patients, with evidence that at low levels, it also predicts low survival in ESRD [4]. We found a significant difference in the values of serum albumin, in which individuals in HD with serum albumin ≥4.0, had an HR of 0.230 (95% CI 0.097–0.541). Many studies have shown that hypoalbuminemia is a good predictor of mortality.

The 10-year cohort of Kato, Castro and Natarajan (2013) [45] indicated a high risk of mortality also in HD patients with serum albumin levels ≤3.8 g/dl; Teixeira et al. (2015) [27] also found that albumin is related to low survival; Msaad et al. (2019) [44] reported that 77.27% of deceased patients had low albumin and Ebhahimi et al. (2019) [4] demonstrated that for each unit (in g/dl) of increase in serum albumin, the survival time for HD patients increased by approximately 23%.

In addition, hypoalbuminemia is considered a marker of malnutrition and a strong predictor of death in these patients. However, serum albumin levels should be used with caution as nutritional markers for these patients, because low levels of this protein in HD patients are also associated with malnutrition, inflammation and, consequently, with other complications such as anemia secondary to ESRD [42, 45, 46].

Patients with ESRD have a series of hematopoietic abnormalities, most commonly, anemia, which is a multifactorial disease that affects almost all patients undergoing HD and is related to the risk of early death [27, 46,47,48]. HD patients are predisposed to iron deficiency due to residual blood loss during treatment [49, 50].

Serum ferritin is also an established marker for detecting absolute serum iron deficiency in patients on dialysis [48]. The 2012 KDIGO [16] guidelines for anemia management recommend that you assess iron stores by measuring ferritin at least every three months, which was done in our study. We found that patients with adequate ferritin values and with high serum iron might be associated with a poor prognosis of survival time in the treatment of HD.

Regarding calcium, individuals with high values (≥ 11.01 mg/dl) had an HR of 4.102 (95% CI 1.351–12.457). In contrast, patients with high serum phosphorus (≥ 6,01 meq/lit) had an HR of 0,290 (95% CI 0.138–0.609). Phosphorus excretion is reduced with the kidney damage that HD users have, which leads to increased levels of phosphorus in the blood. When there is an imbalance of phosphorus, what happened in our study, it leads to a significant loss of calcium and, consequently, to a debilitating bone disease [51]. Although in the advanced stages most patients present hyperphosphatemia, studies show [52] that the reduction in serum phosphorus improves patient survival, which was not found in our study. This can be explained by the fact that as GFR decreases, serum phosphorus levels increase [53].

However, a positive calcium balance arises easily because the intestinal absorption of calcium is greater than the renal excretion capacity, therefore, the excretion of this mineral in the kidneys decreases and does not increase [52]. In addition, because hypocalcemia is common in ESRD, routine calcium supplementation exists. However, high levels are associated with undesirable and harmful effects, such as vascular calcification, cardiovascular disease and mortality [53]. As with our findings, Inaguma et al. (2017) [53] evidenced in their multivariate analysis that high serum calcium levels are an independent risk factor for a poor survival prognosis.

Our study has several strengths, one of which is the complete coverage of patients referred by nephrologists with detailed baseline information on diagnosis, data and laboratory treatments. In addition, we have a long follow-up period (20 years). However, one weakness is missing data for some variables (e.g. cause of mortality, medications used, body weight and patient’s food consumption). Because this is a 20-year retrospective cohort, many medical records were no longer available for access. In this sense, as these data were not found in the medical records, we did not use these data in our analysis. Furthermore, the study was performed in a point-cohort design in Brazil; so, the results on the influence of baseline indicators on patient’s survival time, in particular, may not reflect those from other countries.


In summary, with the presented prediction model of the dynamics of biochemical parameter and clinical parameters at baseline (first month of admission on HD), it is possible to give a reasonably accurate estimation on the survival of patients with HD. The management of those biochemical parameter in the HD population remains challenging, given the proportion of high risk of death in patients with low levels of serum albumin, ferritin, calcium and phosphorus, evidencing the multifactorial and interrelated nature between several mechanisms that contribute to the progression of CKD, such as inflammation and malnutrition.

This prediction models of survival time can aid nephrologists in providing patients with well-founded information on their future prospects. Even this prognosis be only an absolute risk estimate, it can be reassuring for patients to know where they stand. Furthermore, this model can also guide clinical decision-making, identifying patients at high risk and allow optimization of preventive measures. Intensive intervention to reduce death of patients in the first year of RRT treatment is highly recommended, as well as outpatient nephrological care in pre-dialysis. If no known treatment is available, a nephrologist and her team could still change his approach by intensifying follow-up of the high-risk patient.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.



End-stage renal disease


Chronic kidney disease


Replacement renal therapy




Continuous ambulatory peritoneal dialysis


Glomerular Filtration Rate




Rheumatologic profile


Normalized catabolic protein rate


Renal Disease Outcome Quality Initiative


Standard deviation


Interquartile range


Confidence interval


Odds ratio


Red blood cell count


  1. 1.

    Browne OT, Allgar V, Bhandari S. Analysis of factors predicting mortality of new patients commencing renal replacement therapy 10 years of follow-up. BMC Nephrol. 2014;15(1):1–10.

  2. 2.

    Levin A, Tonelli M, Bonventre J, Coresh J, Donner JA, Fogo AB, et al. Global kidney health 2017 and beyond: a roadmap for closing gaps in care, research, and policy. Lancet. 2017;390(10105):1888–917.

  3. 3.

    Crews DC, Bello AK, Saadi G. 2019 world kidney day editorial - burden, access, and disparities in kidney disease. J Bras Nefrol. 2019;41(1):1–9.

    PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Ebrahimi V, Khademian MH, Masoumi SJ, Morvaridi MR, Ezzatzadegan JS. Factors influencing survival time of hemodialysis patients; time to event analysis using parametric models: a cohort study. BMC Nephrol. 2019;20(1):1–9.

    Article  Google Scholar 

  5. 5.

    Foreman KJ, Marquez N, Dolgert A, Fukutaki K, Fullman N, McGaughey M, et al. Forecasting life expectancy, years of life lost, and all-cause and cause-specific mortality for 250 causes of death: reference and alternative scenarios for 2016–40 for 195 countries and territories. Lancet [Internet]. 2018;392(10159):2052–90.

    Article  Google Scholar 

  6. 6.

    Jager KJ, Kovesdy C, Langham R, Rosenberg M, Jha V, Zoccali C. A single number for advocacy and communication-worldwide more than 850 million individuals have kidney diseases. Nephrol Dial Transplant. 2019;34(11):1803–5.

    PubMed  Article  Google Scholar 

  7. 7.

    Kramer A, Pippias M, Noordzij M, Stel VS, Andrusev AM, Aparicio-Madre MI, et al. The European renal association - European Dialysis and transplant association (ERA-EDTA) registry annual report 2016: a summary. Clin Kidney J. 2019;12(5):702–20..

    PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    Thomé FS, Sesso RC, Lopes AA, Lugon JR, Martins CT. Brazilian chronic dialysis survey 2017. Brazilian J Nephrol [Internet]. 2019;41(2):208–14.

    Article  Google Scholar 

  9. 9.

    Alcalde PR, Kirsztajn GM. Expenses of the Brazilian public healthcare system with chronic kidney disease. Brazilian J Nephrol [Internet]. 2018;40(2):122–9.

    Article  Google Scholar 

  10. 10.

    Carrero JJ, Hecking M, Chesnaye NC, Jager KJ. Sex and gender disparities in the epidemiology and outcomes of chronic kidney disease. Nat Rev Nephrol [Internet]. 2018;14(3):151–64.

    Article  Google Scholar 

  11. 11.

    Kainz A, Berner C, Ristl R, Simon A, Stamm T, Zitt E, et al. Sex-specific analysis of haemodialysis prevalence, practices and mortality over time: the Austrian Dialysis registry from 1965 to 2014. Nephrol Dial Transplant. 2019;34(6):1026–35.

    PubMed  Article  Google Scholar 

  12. 12.

    Liabeuf S, Sajjad A, Kramer A, Bieber B, McCullough K, Pisoni R, et al. Guideline attainment and morbidity/mortality rates in a large cohort of European haemodialysis patients (EURODOPPS). Nephrol Dial Transplant. 2019;34(12):2105–10.

    CAS  PubMed  Article  Google Scholar 

  13. 13.

    Villa-Bellosta R, Rodriguez-Osorio L, Mas S, Abadi Y, Rubert M, De La Piedra C, et al. A decrease in intact parathyroid hormone (iPTH) levels is associated with higher mortality in prevalent hemodialysis patients. PLoS One. 2017;12(3):1–12.

    Google Scholar 

  14. 14.

    Riella MC, Martins C. Nutrição e o Rim. Guanabara; 2011.

    Google Scholar 

  15. 15.

    Riella MC, Martins C. Nutrição e o Rim. Guanabara; 2013.

    Google Scholar 

  16. 16.

    Rocco M, Daugirdas JT, Depner TA, Inrig J, Mehrotra R, Rocco MV, et al. KDOQI clinical practice guideline for hemodialysis adequacy: 2015 update. Am J Kidney Dis. 2015;66(5):884–930.

    Article  Google Scholar 

  17. 17.

    Fan Y, Jin X, Jiang M, Fang N. Elevated serum alkaline phosphatase and cardiovascular or all-cause mortality risk in dialysis patients: a meta-analysis. Sci Rep. 2017;7(1):1–8.

    Article  CAS  Google Scholar 

  18. 18.

    Kleinbaum DG, Klein M. Análise de sobrevivência. Nova Iorque: Springer-Verlag; 2010. p. 3.

    Google Scholar 

  19. 19.

    Efron B, Tibshirani RJ. An introduction to the bootstrap. New York: Chapman e Hall; 1993. p. 436.

    Google Scholar 

  20. 20.

    Sauerbrei W, Schumacher M. A bootstrap resampling procedure for model building: application to the cox regression model. Stat Med. 1992;11(16):2093–109.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Diepen MV, Ramspek CL, Jager KJ, Zoccali C, Dekker FW. Prediction versus aetiology: common pitfalls and how to avoid them. Nephrology Dialysis Transplantation. 2017;32(2):1–5.

    Article  Google Scholar 

  22. 22.

    Geloven NV, et al. Prediction meets causal inference: the role of treatment in clinical prediction models. Eur J Epidemiol. 2020;35:619–30.

    PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Nguyen B, Fukuuchi F. Survival rates and causes of death in Vietnamese chronic hemodialysis patients. Ren Replace Ther. 2017;3(1):22.

    Article  Google Scholar 

  24. 24.

    Chandrashekar A, Ramakrishnan S, Rangarajan D. Survival analysis of patients on maintenance hemodialysis. Indian J Nephrol. 2014;24(4):206–2013.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Collins AJ, Foley RN, Herzog C, Chavers B, Gilbertson D, Ishani A, et al. United States renal data system 2008 annual data report abstract. Am J Kidney Dis. 2009;53(1):A6–7.

    Article  Google Scholar 

  26. 26.

    Valavi E, Rekabi F, Mousavi M, Beladi Mousavi S, Hayati F. Comparison of survival in patients with end-stage renal disease receiving hemodialysis versus peritoneal dialysis. Saudi J Kidney Dis Transplant. 2015;26(2):392.

    Article  Google Scholar 

  27. 27.

    Teixeira FIR, Lopes MLH, Silva GA dos S, Santos RF. Survival of hemodialysis patients at a university hospital. J Bras Nefrol [Internet]. 2015;37(1):64–71.

    Google Scholar 

  28. 28.

    Magalhães LP. Análise dos fatores preditores de mortalidade em pacientes incidentes em hemodiálise; 2015.

    Google Scholar 

  29. 29.

    dos Santos VFC, Borges ZN, Lima SO, Reis FP. Percepções, significados e adaptações à hemodiálise como um espaço liminar: a perspectiva do paciente. Interface - Comun Saúde, Educ. 2018;22(66):853–63.

    Article  Google Scholar 

  30. 30.

    Brodowska-Kania D, Rymarz A, Gibinski K, Kielczewska J, Smoszna J, Saracyn M, et al. First year survival of patients on maintenance dialysis treatment in Poland. Nagoya J Med Sci. 2015;77(4):629–35.

    PubMed  PubMed Central  Google Scholar 

  31. 31.

    Alicic RZ, Rooney MT, Tuttle KR. Diabetic kidney disease: challenges, Progress, and possibilities. Clin J Am Soc Nephrol [Internet]. 2017;12(12):2032–45.

    CAS  Article  Google Scholar 

  32. 32.

    Papadopoulou-Marketou N, Kanaka-Gantenbein C, Marketos N, Chrousos GP, Papassotiriou I. Biomarkers of diabetic nephropathy: a 2017 update. Crit Rev Clin Lab Sci. 2017;54(5):326–42.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Amorim RC, Guedes GS, Vasconcelos SML, Santos JCF. Doença Renal do Diabetes: Cross-Linking entre Hiperglicemia, Desequilíbrio Redox e Inflamação. Arq Bras Cardiol. 2019;112(5):577–87.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Pugliese G, Penno G, Natali A, et al. Diabetic kidney disease: new clinical and therapeutic issues. Joint position statement of the Italian Diabetes Society and the Italian Society of Nephrology on "The natural history of diabetic kidney disease and treatment of hyperglycemia in patients with type 2 diabetes and impaired renal function". J Nephrol. 2020;33(1):9–35.

  35. 35.

    Hsu RK, Chai B, Roy JA, Anderson AH, Bansal N, Feldman HI, et al. Abrupt decline in kidney function before initiating hemodialysis and all-cause mortality: the chronic renal insufficiency cohort (CRIC) study. Am J Kidney Dis [Internet]. 2016;68(2):193–202.

    Article  Google Scholar 

  36. 36.

    Reddan DN, Klassen PS, Szczech LA, Coladonato JA, O’Shea S, Owen WF, et al. White blood cells as a novel mortality predictor in haemodialysis patients. Nephrol Dial Transplant. 2003;18(6):1167–73.

    CAS  PubMed  Article  Google Scholar 

  37. 37.

    Hauser AB, Stinghen AEM, Kato S, Bucharles S, Aita C, Yuzawa Y, et al. Characteristics and causes of immune dysfunction related to uremia and dialysis. Perit Dial Int. 2008;28(SUPP. 3):S183–7.

    PubMed  Article  Google Scholar 

  38. 38.

    Davenport A. Differences in dialysis center practices in determining hemodialysis patient postdialysis target weight and patient survival and hospitalizations. Clin J Am Soc Nephrol. 2019;14(3):321–3.

    PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    Akchurin OM, Kaskel F. Update on inflammation in chronic kidney disease. Blood Purif. 2015;39(1–3):84–92.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Kaupke CJ, Zhang J, Cesario T, Yousefi S, Akeel N, Vaziri ND. Effect of hemodialysis on leukocyte adhesion receptor expression. Am J Kidney Dis [Internet]. 1996;27(2):244–52.

    CAS  PubMed  Article  Google Scholar 

  41. 41.

    Liakopoulos V, Jeron A, Shah A, Bruder D, Mertens PR, Gorny X. Hemodialysis-related changes in phenotypical features of monocytes. Sci Rep. 2018;8(1):1–12.

    Article  CAS  Google Scholar 

  42. 42.

    Alfonso A, Castillo R, Gomez Jimenez F, Negrillo A. Evaluation of nutritional biochemical parameters in Haemodialysis patients over a ten-year period. West Indian Med J. 2015;64(3):213–7.

    PubMed  Google Scholar 

  43. 43.

    Zha Y, Qian Q. Protein nutrition and malnutrition in CKD and ESRD. Nutrients. 2017;9(3):1–19.

    Article  CAS  Google Scholar 

  44. 44.

    Msaad R, Essadik R, Mohtadi K, Meftah H, Lebrazi H, Taki H, et al. Predictors of mortality in hemodialysis patients. Pan Afr Med J. 2019;33:1–14.

    Article  Google Scholar 

  45. 45.

    Kato M, Castro NC, Natarakan R. MicroRNAs:potentialmediatorsandbiomarkers of diabeticcomplications. Free Radic Biol Med. 2013;64:85–94.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Kamal Al shafei N, Nour A. CBC, serum proteins, and Immunoglobulins in chronic hemodialysis patients with or without pruritus in Egypt. Biochem Anal Biochem. 2016;05(01):1–7.

    Article  CAS  Google Scholar 

  47. 47.

    Cases A, Egocheaga MI, Tranche S, Pallarés V, Ojeda R, Górriz JL, et al. Anemia of chronic kidney disease: protocol of study, management and referral to nephrology. Nefrologia. 2018;38(1):8–12.

    PubMed  Article  Google Scholar 

  48. 48.

    Diebold M, Kistler AD. Evaluation of iron stores in hemodialysis patients on maintenance ferric Carboxymaltose dosing. BMC Nephrol. 2019;20(1):1–9.

    Article  Google Scholar 

  49. 49.

    Tsukamoto T, Matsubara T, Akashi Y, Kondo M, Yanagita M. Annual Iron loss associated with hemodialysis. Am J Nephrol [Internet]. 2016;43(1):32–8.

    CAS  Article  Google Scholar 

  50. 50.

    Son R, Fujimaru T, Kimura T, Taki F, Futatsuyama M, Nagahama M, et al. Association between serum ferritin levels and clinical outcomes in maintenance hemodialysis patients: a retrospective single-center cohort study. Ren Replace Ther. 2019;5(1):1–8.

    Article  Google Scholar 

  51. 51.

    Sulaiman MK. Diabetic nephropathy: recent advances in pathophysiology and challenges in dietary management. Diabetol Metab Syndr. 2019;11(1):7.

    PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Langman CB, Cannata-Andía JB. Calcium in Chronic Kidney Disease: Myths and Realities. Clin J Am Soc Nephrol. 2010;5(Supplement 1):S1–2.

    PubMed  Article  Google Scholar 

  53. 53.

    Inaguma D, Koide S, Takahashi K, Hayashi H, et al. Relationship between serum calcium level at dialysis initiation and subsequent prognosis. Renal Replacement Therapy. 2017;3:2.

    Article  Google Scholar 

Download references


The authors greatly appreciate all the investigators for their contribution to the trial, the statistical support and all the staff for their efforts in treating of patients and thank all patients involved in this study.


Not applicable.

Author information




E.S.F. and R.M.M.C. contributed to conception, design, conduct of the trial and interpretation of data, drafting and revision of the manuscript for important intellectual content, approval of the final version for submission and agreement to be accountable for all aspects of the work. T.R.M., L.S.S., G.D.C. and R.G.S. contributed to the interpretation of data, revision of the manuscript for important intellectual content, approval of the final version for submission and agreement to be accountable for all aspects of the work. S.C.B.O., B.O.S., H.H.D., L.D.B., J.C.M. contributed to interpretation of data, revision of the manuscript for important intellectual content, approval of the final version for submission and agreement to be accountable for all aspects of the work. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Emily de Souza Ferreira.

Ethics declarations

Ethics approval and consent to participate

This study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving research study participants were approved by the Ethics Committee on Research with Human Beings (CEP) of the Federal University of Viçosa under the report number 2.459.555. Written informed consent was obtained from all subjects/patients.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ferreira, E.d., Moreira, T.R., da Silva, R.G. et al. Survival and analysis of predictors of mortality in patients undergoing replacement renal therapy: a 20-year cohort. BMC Nephrol 21, 502 (2020).

Download citation


  • Hemodialysis
  • Renal replacement therapy
  • End-stage renal disease
  • Chronic kidney disease; biochemical parameter