Abstract

Background. Patients with congestive heart failure (CHF) are vulnerable to contrast-induced kidney injury (CI-AKI), but few prediction models are currently available. Therefore, we aimed to establish a simple nomogram for CI-AKI risk assessment for patients with CHF undergoing coronary angiography. Methods. A total of 1876 consecutive patients with CHF (defined as New York Heart Association functional class II-IV or Killip class II-IV) were enrolled and randomly (2:1) assigned to a development cohort and a validation cohort. The endpoint was CI-AKI defined as serum creatinine elevation of ≥0.3 mg/dL or 50% from baseline within the first 48–72 hours following the procedure. Predictors for the simple nomogram were selected by multivariable logistic regression with a stepwise approach. The discriminative power was assessed using the area under the receiver operating characteristic (ROC) curve and was compared with the classic Mehran score in the validation cohort. Calibration was assessed using the Hosmer–Lemeshow test and 1000 bootstrap samples. Results. The incidence of CI-AKI was 9.06% (170) in the total sample, 8.64% (n = 109) in the development cohort, and 9.92% (n = 61) in the validation cohort (). The simple nomogram including four predictors (age, intra-aortic balloon pump, acute myocardial infarction, and chronic kidney disease) demonstrated a similar predictive power as the Mehran score (area under the curve: 0.80 vs. 0.75, ), as well as a well-fitted calibration curve. Conclusions. The present simple nomogram including four predictors is a simple and reliable tool to identify CHF patients at risk of CI-AKI, whereas further external validations are needed.

1. Introduction

Contrast-induced kidney injury (CI-AKI) is a common complication following coronary angiography (CAG) or percutaneous coronary intervention (PCI), and the occurrence of CI-AKI has been demonstrated to be related to poor outcomes [13]. Published guidelines recommend hydration therapy during the perioperative period of coronary catherization as one of the important strategies to prevent CI-AKI [4, 5]. However, for patients with congestive heart failure (CHF), the extra fluid load might lead to further adverse events of acute heart failure [6, 7], even though hydration volume is halved in such patients. Therefore, precise stratification and identification of patients at high risk of CI-AKI are required to perform preventive hydration accurately to reduce unnecessary overload complications. Previous studies have reported various prediction models for CI-AKI, but few of them were developed for CHF patients specifically, even though CHF patients are vulnerable to CI-AKI [8]. In addition, the classic Mehran score, which includes 8 variables, might be too complicated for clinical application [9]. Therefore, in this study, we aimed to establish a simple nomogram for CI-AKI risk assessment among patients with CHF undergoing CAG/PCI.

2. Materials and Methods

2.1. Patients

The current study population was based on a prospective observation cohort (PREdictive Value of COntrast voluMe to creatinINe Clearance Ratio, PRECOMIN, NCT01400295) which enrolled consecutive patients undergoing CAG/PCI in Guangdong Provincial People’s Hospital between January 2010 and October 2012. In this study, patients with CHF (defined as New York Heart Association (NYHA) functional class II-IV or Killip class II-IV) undergoing CAG/PCI were enrolled. [10] The exclusion criteria were pregnancy, lactation, contrast exposure within the 7 days before or 3 days after the procedure, cardiovascular surgery, no use of low-osmolarity contrast agents, undergoing hemodialysis, missing preoperative or postoperative creatinine, malignancy, and no use of isotonic saline for hydration). The Ethics Committee of the Guangdong Provincial People’s Hospital approved this study. All the patients involved provided written informed consent and were randomly assigned to a development cohort and a validation cohort in a 2:1 ratio.

2.2. Coronary Angiography and Laboratory Examination

The procedures were conducted by interventional cardiologists following published guidelines, institutional policy, and clinical routine. [11] Patients undergoing a non-emergency procedure received hydration therapy (0.5–1 mL/kg/h for at least 2–12 hours before and 6–24 hours after the procedure). Patients undergoing emergency procedures received unspecified hydration therapy before the procedure. Serum creatinine (SCr) was measured for all patients at admission and at 1, 2, and 3 days after the procedure using the Jaffe method.

2.3. Endpoint and Definitions

The endpoint of this study was CI-AKI defined as a SCr elevation of ≥0.3 mg/dL or 50% from baseline within the first 48 to 72 hours following contrast exposure. [12] Baseline characteristics, angiographic data, and medications were prospectively defined and have been reported in a previous study [13]. The definitions of hypotension, diabetes, anemia, chronic kidney disease (CKD), and intra-aortic balloon pump (IABP) were the same as those used for the Mehran score [9]. Acute myocardial infarction (AMI) was defined according to the third universal definition of myocardial infarction [14]. Hypoalbuminemia was defined as serum albumin <35 g/L [15]. Age, weight, heart rate (HR), and contrast volume were defined as continuous variables. Angiotensin-converting enzyme inhibitors (ACEIs)/angiotensin-receptor blockers (ARBs), diuretics, and beta blockers were defined as the prescription of these medicines during the perioperative period. To calculate the estimated glomerular filtration rate (eGFR), the Modification of Diet in Renal Disease (MDRD) equation was used: eGFR (mL/min/1.73 m2) = 186 (SCr ^  − 1.154) (Age ^  − 0.203) 0.742 (if female). The unit of SCr in this formula is mg/dl. CKD was defined as eGFR <60 mL/min/1.73 m2. Patients were further divided into 5 CKD stages according to the guidelines (CKD G1: eGFR ≥ 90 mL/min/1.73 m2; CKD G2: eGFR: 60–90 mL/min/1.73 m2; CKD G3: eGFR: 30–60 mL/min/1.73 m2; CKD G4: eGFR: 15–30 mL/min/1.73 m2; and CKD G5: eGFR: <15 mL/min/1.73 m2). [16] And they were also divided into 3 categories according to their heart function based on NYHA or Killip class. All eligible patients enrolled were followed up at 1 month, 6 months, and every 1 year after enrollment until April 2019.

2.4. Statistical Analysis

All the patients involved were randomly assigned to a development cohort and a validation cohort in a 2:1 ratio. Continuous variables were compared with an unpaired, 2-tailed t-test and are expressed as the mean ± SD or were compared through the Wilcoxon rank-sum test and are expressed as the median ± interquartile. Categorical variables were compared using the χ2 test or Fisher’s exact test and are expressed as percentages. Kaplan–Meier curve and multivariable Cox proportional hazard regression adjusted for known risk factors in the context of long-term prognosis were performed to explore the association between CI-AKI and long-term mortality. [17, 18].

To build the nomogram, candidate variables that were imbalanced between groups in the development cohort or that are clinically important, such as risk factors included in the traditional Mehran score, were included in the univariable logistic analysis. Variables with >15% missing values were not considered candidates, i.e., low-density lipoprotein-C and HbA1c. Significant variables from the univariable logistic analysis were then included in the multivariable logistic analysis. A backward stepwise approach was performed to screen the variables by successively removal of nonsignificant () covariates until all the remaining variables were statistically significant. Then, we manually investigated the contribution of the remaining variables to determine the final predictors. Collinearity between variables was also evaluated. A nomogram was then formulated based on the results and by using the rms package of R. To form the nomogram, each regression coefficient in the multivariate logistic regression was proportionally converted into a 0- to 100-point scale. The variable with the highest β coefficient (absolute value) was assigned 100 points. The points are added across each variable to calculate the total points, which are finally converted to predicted probabilities. The performance of the nomogram was assessed using the area under the receiver operating characteristic (ROC) curve and concordance C-statistic for discriminative ability. For validation, a bootstrap method (1000 times) was performed in both the development and validation cohorts to evaluate the stability of the C-statistic. Calibration was assessed using the Hosmer–Lemeshow test and the 1000 bootstrap samples to decrease the overfit bias [19, 20]. Area under curve (AUC) comparison between the nomogram and the Mehran score was performed using DeLong’s test in the validation cohort. Missing data were not imputed. In all analyses, was considered statistically significant. All analyses were conducted with R software (version 3.6.2; R Foundation for Statistical Computing, Vienna, Austria) and SPSS (version 26.0).

3. Results

3.1. Baseline Characteristics

The details and missing data of the included patients are listed in Supplementary Table 1. Among the included 1876 CHF patients, 1261 and 615 patients were divided into the development and validation cohorts, respectively. In the total cohort, approximately one-fourth were female (25.11%). The mean age was 64.77 ± 10.70 years, and the mean SCr was 96.86 ± 50.73 μmol/L. No significant differences were identified between the development and validation cohorts, except for the diastolic blood pressure, contrast volume, and the proportion of percutaneous coronary intervention (PCI) during the procedure.

The incidence of CI-AKI was 9.06% (n = 170) in the total sample, 8.64% (109 patients) in the development cohort, and 9.92% (61 patients) in the validation cohort (). The distribution of the SCr value within the time frame (24–72 hours after the procedure) is shown in Supplementary Figure 1. For the development cohort, patients complicated with CI-AKI following CAG tended to be older and had lower left ventricular ejection fraction (LVEF), haemoglobin, and eGFR than those without CI-AKI. Patients with CI-AKI were also more likely to have hypotension, diabetes, CKD, AMI, and intra-aortic balloon pump (IABP), and they were less likely to be prescribed ACEI/ARB and beta blockers. No significant difference between groups was identified in contrast volume (Table 1).

3.2. Development and Validation of the CI-AKI-Predicting Nomogram

The results of univariate logistic analysis are detailed in Table 2. Through multivariate logistic analysis and a backward stepwise approach, age (OR: 1.04 95%, CI: 1.02–1.06), IABP (OR: 3.97 95%, CI: 2.14–7.37), AMI (OR: 3.27 95%, CI: 2.09–5.10), and CKD (OR: 2.83 95%, CI: 1.80–4.44) were selected as predictors of CI-AKI (Table 3). The Hosmer–Lemeshow statistic of multivariable analysis suggested a good fit (χ2 = 10.78, ).

A simple nomogram based on the selected predictors was formed (Figure 1). The nomogram was internally validated with the bootstrap validation method (1000 times). In the development cohort, the nomogram demonstrated good discriminative power for estimating the risk of CI-AKI, with an unadjusted C-statistic of 0.79 (95% CI, 0.75–0.84) and a bootstrap-corrected C statistic of 0.79. In addition, calibration plots graphically showed good agreement on the presence of CI-AKI between the risk estimation and the observed frequency (Figure 2(a)).

In the validation cohort, the simple nomogram demonstrated a C-statistic of 0.80 (95% CI 0.75–0.86), which was similar to that of the Mehran score (AUC: 0.75, 95% CI 0.68–0.81) among patients with CHF (, Figure 3) and a bootstrap-corrected C statistic of 0.80. Moreover, there was also a good calibration curve for the risk estimation (Figure 2(b)).

The discriminative ability and calibration were also tested among patients with different stages of CHF and stages of CKD, respectively. No significant difference in discriminative ability was identified, except for the comparison between “Heart function II” and “Heart function IV.” Hosmer–Lemeshow statistic suggested a good fit in all stages of CHF and stages of CKD (Supplementary Table 2).

3.3. Risk of CI-AKI Based on the Nomogram Scores

Based on the predicted incidence of CI-AKI in relation to different total nomogram scores, we further divided the patients into 5 score categories: scores < 50 (risk = 1.96%), 50 ≤ scores < 100 (risk = 4.71%), 100 ≤ scores < 150 (risk = 13.90%), 150 ≤ scores < 200 (risk = 39.09%), and scores ≥ 200 (risk = 65.78%). The rates of CI-AKI in the validation cohort were close to those in the development cohort inside each of the 5 score categories (Supplementary Figure 2).

The optimal cut-off value of the total nomogram scores to identify patients at risk was determined to be 100. The sensitivity and specificity used to differentiate the presence from the absence of CI-AKI were 69.7% and 75.6% in the development cohort and 70.0% and 74.5% in the validation cohort, respectively. The total nomogram scores of 100 to identify patients at risk also demonstrated a good discriminative power in both the development cohort (AUC: 0.73, 95% CI: 0.68–0.78, and ) and the validation cohort (AUC: 0.72, 95% CI: 0.65–0.79, and ).

4. Discussion

The present study might be the first to develop a simple nomogram for the prediction of CI-AKI among CHF patients undergoing CAG/PCI. Through multivariable logistic regression analysis and a stepwise approach, a nomogram based on four factors (age, IABP use, AMI, and CKD) was constructed. In particular, the simple nomogram demonstrated a similar discriminative power as the classic Mehran score in the validation cohort, as well as good stability and calibration.

To define CHF patients in this study, two classification methods (NYHA and Killip) were used interchangeably. Previous studies had already taken this method to identify CHF patients [21, 22]. In our clinical routine, we evaluate the heart function of AMI and non-AMI patients with Killip and NYHA system, respectively. Using these two systems interchangeably may help identifying patients with impaired heart function more concisely.

In this study of the CHF population, the overall incidence of CI-AKI was 9.06%. The previously reported incidence of CI-AKI has not been consistent among different studies. Among the patients with a low-risk of ordinary coronary artery disease, the CI-AKI incidence was as low as 2-3%, while in the high-risk population, the CI-AKI incidence could be as high as 19% or 50% [2, 3]. Since heart failure patients might be complicated with haemodynamic disturbances, which might further cause kidney dysfunction [23, 24], the relatively high incidence of CI-AKI among our CHF cohort was mostly rational. Moreover, we noticed that some patients were missing SCr data on day 3, which may lead to the underestimation of the incidence of CI-AKI. However, many previous studies had limited the diagnosis time frame of CI-AKI to 48 hours, or even 24 hours [9, 25]. In the study conducted by Kim et al., CI-AKI was defined as an increase in SCr of ≥50% or 0.3 mg/dL, or a decrease in eGFR of ≥25% within 24 hours after PCI [25].

The present simple nomogram was constructed based on age, AMI, IABP, and CKD, which are common risk factors associated with CI-AKI, as reported by numerous previous studies [12, 26, 27]. Thomas T. Tsai et al. conducted a study including 985,737 consecutive patients who underwent PCI at 1,253 sites in America and found that ST-segment elevation myocardial infarction (STEMI) (OR: 2.60; 95% CI: 2.53–2.67) as well as non-ST-segment elevation myocardial infarction (NSTEMI) (OR: 1.81; 95% CI: 1.61–2.04) were independently associated with CI-AKI. Additionally, each 10-year increase in age (OR: 1.15; 95% CI: 1.14–1.16), IABP before the procedure (OR: 2.13; 95% CI: 1.92–2.35), and CKD stage 3–5 were all risk factors for CI-AKI [2]. Pierre Aubry et al. performed a retrospective cross-sectional population-based study involving 1,047,329 cases of contrast exposure in Europe to quantify the effect of risk factors for CI-AKI. By multivariate analysis, age >80 years (OR: 2.7; 95% CI: 2.6–2.8) and CKD (OR: 2.3; 95% CI: 2.2–2.3) were identified as independent risk factors for CI-AKI [28]. For patients in Asia, Pei-Chun Fan et al. developed and validated a risk prediction model (ADVANCIS score) for incident CI-AKI based on 82,186 patients admitted for acute coronary syndrome (ACS) receiving PCI. Through multivariable logistic regression analysis, age (OR: 1.02; 95% CI: 1.01–1.02), CKD (OR: 11.38; 95% CI: 9.86–13.13), and IABP use (OR: 1.80; 95% CI: 1.62–2.00) were strongly associated with the incidence of CI-AKI. [29].

The nomogram developed in this study is a simple but efficient prediction model that includes only four predictors. We noticed that some important risk factors reported by previous results were not included in the current simple nomogram, such as contrast volume and diabetes [30]. The nonsignificant difference between the CI-AKI and non-CI-AKI subjects in contrast volume may be one of the reasons. In addition, some recent studies have found that the dose of contrast agent is not necessarily related to the incidence of CI-AKI [31, 32]. Also, with the development of coronary interventional treatment techniques, the volume of contrast agents is gradually reduced [9]. Moreover, in the current study, all patients included were restricted to use low-osmolarity contrast agents during the procedure, and the risk of kidney impairment was reduced by the usage of low-osmolarity and isotonic contrast agents [33]. The risk of CI-AKI might be due to haemodynamic instability rather than direct injury by contrast agents. Regarding diabetes mellitus, a long course of hyperglycaemia would result in renal microcirculation disorder and renal dysfunction, which might affect the incidence of CI-AKI [34, 35]. However, there is no evidence suggesting that short-term hyperglycaemia associated with diabetes mellitus without renal dysfunction would increase the risk of CI-AKI. Therefore, the novel model we developed directly included CKD as one of the predictors rather than diabetes mellitus. Moreover, adding contrast volume and diabetes to the nomogram did not significantly increase the C-statistic (>10%).

The Mehran score is one of the most classic CI-AKI risk estimating system, which was developed based on 8,357 patients undergoing PCI [9]. Comparing with the Mehran score, the current nomogram has several strengths. First, our previous study showed that the CI-AKI definition of this study has higher population attributable risk (PAR) than that of the CI-AKI definition of the Mehran score, which means the current CI-AKI definition has better prognostic value [36]. Second, the current nomogram with only 4 variables was established for CHF patients specifically. Less variables lead to easier clinical implication. However, the nomogram also has some shortness. First, the nomogram was developed based on a relatively small cohort. Second, the nomogram lacks external validation.

We also noticed that several researches have provided nomograms to predict acute kidney injury (AKI) lately. However, these nomograms were developed for different patients specifically [3739]. By analysing the Medical Information Mart for Intensive Care- (MIMIC-) III v. 1.4 database, Deng et al. established a nomogram to predict the risk of septic AKI within the first 24 h after admission of intensive care unit. Among the included 2,917 sepsis patients without CKD, the nomogram with 7 variables demonstrated well-fitted calibration curves and good C-indexes in both the training and validation cohorts [37]. In another study conducted by Xu et al., they developed a preprocedural nomogram with 5 variables to predict the risk of AKI among patients undergoing nephrectomy [38]. Comparing with these nomograms, our nomogram was developed specifically for the CHF population undergoing CAG/PCI. Moreover, with fewer and readily accessed variables, our nomogram was easier for clinical implication in the department of cardiology.

Overall, in this study, we established a novel CI-AKI risk prediction nomogram specified for the CHF population. In this simple nomogram, patients with more than 100 total points are identified as at-risk patients, accounting for nearly 30% of the entire CHF cohort. The threshold of 100 points (approximately 10% risk) was also similar to the risk threshold of previous models [9, 29].

4.1. Limitations

Our study had some limitations. First, this study was an analysis of single-centre data and the nomogram we established is not a general one for CI-AKI risk evaluation in any radiological study, but specifically for coronary study/intervention. However, no relevant study had so far elucidated the specified risk factors of CI-AKI in CHF patients. Our study provided a new tool for CI-AKI risk assessment among CHF patients, prompting more attention and further studies on risk assessment of patients with comorbidities undergoing coronary catheterization. Second, the present nomogram was not externally validated, although there was a calibration with 1000 bootstrap samples to decrease the overfit bias. Third, our nomogram did not demonstrate better discriminative power than the Mehran score. However, our model has less variables, which makes it easier for clinical implement. Fourth, some other potential risk factors, especially procedural context (i.e., emergency vs. elective procedure), were not included due to the insufficient data. Finally, some patients were discharged within 72 hours after the CAG/PCI, so creatinine levels were not measured on day 3 in these patients, which might lead to the underestimation of the incidence of CI-AKI.

5. Conclusions

The presented nomogram with four predictors (age, IABP, AMI, and CKD) is a simple and reliable tool for CI-AKI risk stratification among CHF patients, which enables physicians to identify patients at risk and implement precise prevention strategies in time. However, further external validations are needed before clinical generalization.

Data Availability

Data relevant to this study are available from the corresponding authors upon reasonable request.

Disclosure

This article has been presented at the European Society of Cardiology (ESC) Congress 2020 as a conference abstract. The funding body did not participate in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Li Lei and Yibo He contributed equally to this work. LL, YH, JL, ZN, SC, JC, and YL substantially contributed to the conception and design of the study; LL, YH, ZG, BL, JL, GC, LL, ML, and WY involved in data collection; LL, YH, JL, ZN, SC, and YL performed data analysis and/or interpretation of data for the work; and LL, YH, YL, and JC drafted the work or revised it critically for important intellectual content; all the authors approved the final version to be published.

Acknowledgments

The authors thank ESC Congress 2020 for giving them the chance to share their work. This work was supported by the Beijing Lisheng Cardiovascular Pilot Foundation (Grant no. LHJJ201612127), the Guangdong Provincial People’s Hospital Dengfeng Project Fund (Grant no. DFJH201919), the Progress in Science and Technology Project of Guangzhou (Grant no. 201904010470), the Guangdong Provincial People’s Hospital Dengfeng Project Fund (Grant no. DFJH2020026), and the China Youth Clinical Research Fund (Grant no. 2017-CCA-VG-020).

Supplementary Materials

Supplementary Table 1. Baseline characteristics. Supplementary Table 2. C-statistic and Hosmer–Lemeshow statistic of the nomogram in patients with different stages of congestive heart failure and stages of chronic kidney disease. Supplementary Figure 1. The distribution of the serum creatinine value within 24–72 hours after the procedure. Supplementary Figure 2. Risk of contrast-induced acute kidney injury according to the total nomogram score categories. (Supplementary Materials)