Abstract

Purpose. Acute kidney injury (AKI) is a major and severe complication following donation-after-circulatory-death (DCD) liver transplantation (LT) and is associated with increased postoperative morbidity and mortality. However, the risk factors and the prognosis factors of AKI still need to be further explored, and the relativity of intraoperative hepatic blood inflow (HBI) and AKI following LT has not been discussed yet. The purpose of this study was to investigate the correlation between HBI and AKI and to construct a prediction model of early acute kidney injury (EAKI) following DCD LT with the combination of HBI and other clinical parameters. Methods. Clinical data of 132 patients who underwent DCD liver transplantation at the first hospital of China Medical University from April 2005 to March 2017 were analyzed. Data of 105 patients (the first ten years of patients) were used to develop the prediction model. Then we assessed the clinical usefulness of the prediction models in the validation cohort (27 patients). EAKI according to Kidney Disease Improving Global Outcomes (KDIGO) criteria based on serum creatinine increase during 7-day of postoperative follow-up. Results. After Least Absolute Shrinkage and Selection Operator (LASSO) regression and simplification, a simplified prediction model consisting of the Child-Turcotte-Pugh (CTP) score (p=0.033), anhepatic phase (p=0.014), packed red blood cell (pRBC) transfusion (p=0.027), and the HBI indexed by height (HBI/h) (p=0.002) was established. The C-indexes of the model in the development and validation cohort were 0.823 [95% CI, 0.738-0.908] and 0.921 [95% CI, 0.816-1.000], respectively. Conclusions. In this study, we demonstrated the utility of HBI/h as a predictor for EAKI following DCD LT, as well as the clinical usefulness of the prediction model through the combination of the CTP score, anhepatic phase, pRBC transfusion and HBI/h.

1. Introduction

Acute kidney injury (AKI) is a major and severe complication following orthotopic liver transplantation and is associated with poor graft survival and increased mortality [1, 2]. The incidence of AKI after liver transplantation (LT) is high and ranges from 20 to 64% [1, 38].

Currently, there is no effective therapy or preventive strategy available for AKI after LT [9], although several promising strategies have been investigated [1012]. Therefore, to predict the occurrence of AKI after LT and treat AKI as early as possible will significantly improve the outcome of liver transplantation.

Although a number of studies have evaluated AKI after LT, the clinical risk factors, especially intraoperative risk factors, still need to be further explored. Most of the studies reported risk factors of AKI after LT including longer anhepatic phase [13], larger red blood cell (RBC) transfusion amount [5, 14], and higher Child-Turcotte-Pugh (CTP) score [15]. In addition, some hemodynamic factors such as intraoperative mean arterial blood pressure [16], the use of intraoperative venovenous bypass technology [17], central venous pressure, right ventricular end-diastolic volume, and mixed venous oxygen saturation [6] have been reported to significantly influence the occurrence of postoperative AKI.

Twenty-five percent of the cardiac output flows through the liver, and hepatic blood inflow (HBI) changes will influence systemic hemodynamics [18]. However, whether intraoperative HBI including hepatic artery flow (HAF) and portal vein flow (PVF) following reconstruction has a correlation with AKI following LT is unknown.

Therefore, the aim of this study was to investigate the relationship between HBI and early acute kidney injury (EAKI) after DCD LT as well as develop and validate a prediction model for EAKI after DCD LT using preoperative and intraoperative factors including HBI.

2. Materials and Methods

2.1. Design and Patients

This retrospective observational study was approved by the institutional review board of our institution (the First Hospital of China Medical University). We retrospectively reviewed the electronic medical records of 156 consecutive adult patients who underwent DCD LT at our institution between April 2005 and March 2017. The need for informed consent was waived given the study’s retrospective design. Patients without complete intraoperative hepatic blood flow data or postoperative serum creatinine monitoring data were excluded (n = 15). Patients with renal replacement therapy (RRT) were excluded (n = 1). Patients who underwent combined liver-kidney transplantation (n = 2) or died within 48 hours postoperatively (n = 5) were excluded. Patients with perirenal hematoma after LT (n=1) were excluded. The remaining 132 patients were analyzed. Then, 105 patients were selected from April 2005 to March 2015 to form the development cohort of this study. The remaining 27 patients from April 2015 to March 2017 were evaluated to form a validation cohort.

2.2. Definition of EAKI

AKI was defined according to the KDIGO criteria [19]: an increase in serum creatinine by 0.3 mg/dL ( ≥26.5 μmol/L) within 48 hours or an increase in creatinine to ≥ 1.5 times baseline within the first 7 postoperative days, which have been validated in patients undergoing LT [6, 7, 20]. We determined postoperative early acute kidney injury (EAKI) according to the KDIGO criteria mentioned above during 7 days of postoperative follow-up. A urine output criterion was not used. All patients who met the KDIGO criteria within 7 days after LT were classified as having EAKI.

2.3. Intraoperative Hepatic Blood Flow Monitoring

Our center used Transonic HT313 Flowmeter (Ithaca, New York, USA) to measure the reconstructed HAF and PVF before abdominal closure. The HAF is the sum of the mean blood flow of reconstructed hepatic arteries. The PVF represents intraoperative mean portal vein flow. The HBI is the sum of HAF and PVF.

2.4. Statistical Analysis
2.4.1. Patient Characteristics

The Shapiro test was used to determine the normality of the data. Either Fisher’s exact test or χ2 test was used to compare the categorical variables between EAKI and non-EAKI patients. Comparisons of continuous variables between EAKI and non-EAKI patients were performed with Student’s T tests or Mann-Whitney U tests.

2.4.2. Prediction Model Building

The Least Absolute Shrinkage and Selection Operator (LASSO) regression [21] was used to select predictors and eliminate multicollinearity. Logistic regression was used to identify univariate and multivariate predictors for EAKI. Model A was developed in the development cohort by multivariable logistic regression whose variables were selected by LASSO regression based on the one standard error criteria [22]. Model B is a simplified model of model A.

2.4.3. Models Evaluation

Calibration of the models was assessed using Hosmer-Lemeshow goodness-of-fit statistics. To evaluate the prediction performance of the models, the Harrell’s concordance index (C-index) was measured in both the development and validation cohort. A calibration curve was plotted to compare the agreement between observed outcomes (Y-axis) and the predictions of the model (X-axis). Finally, to facilitate calculation of individualized risk in clinical practice, model B was converted to an easy-to-use nomogram.

2.4.4. Software

R software (Version 3.5.1, http://www.r-project.org) was used for the statistical analysis and plot generation, with rms, foreign, pROC, psych, glmnet, and ResourceSelection packages. P < 0.05 was considered statistically significant.

3. Results

3.1. Patient Characteristics

Patient demographics and perioperative variables according to the diagnosis of EAKI are presented in Table 1. EAKI occurred in 31 patients (29.5%) in the development cohort. There were significant differences in liver tumor (p=0.027), preoperative total bilirubin level (p=0.003), preoperative MELD score (p=0.015), preoperative CTP score (p=0.003), cold ischemic time (CIT) (p=0.028), operation time (OT) (p=0.015), anhepatic phase (AP) (p=0.002), packed red blood cell (pRBC) transfusion (p=0.005), cryoprecipitate transfusion (p=0.035), PVF (p=0.015), and HBI (p=0.015) between EAKI and non-EAKI patients in the development cohort. However, for HAF, there was no significant difference (p=0.905) between EAKI and non-EAKI patients.

3.2. Prediction Model Building

As shown in Figures 1(a) and 1(b), all continuous and categorical variables in the development cohort were performed by LASSO regression analysis using 10-fold cross-validation and the one standard error criterion. When parameter λ corresponds to the minimum-deviance within one standard error, 6 variables including CTP score, CIT, AP, HBI/h, OT, and pRBC transfusion were obtained. Then, the 6 variables were selected to build model A using multivariable logistic regression. Model B, a simplified model of model A, was built via eliminating CIT (p=0.211) and OT (p=0.212) from model A (Table 2). After the Hosmer-Lemeshow test, model A and model B were a good fit for the data (Table 3). The two prediction models can be shown in the following equations.

Model A

Model B: probability of postoperative EAKI occurrence.

3.3. Prediction Model Performance Assessment

The receiver operating characteristic (ROC) curves of model A and model B in the development and validation cohort were shown in Figure 2. In the development cohort, the C-indexes of model A and B were 0.847 [95% confidence interval (CI), 0.765-0.928] and 0.823 [95% CI, 0.738-0.908], respectively, and there was no significant difference (p = 0.167, DeLong’s test) in the C-indexes between the two models. In the validation cohort, the C-indexes of model A and B were 0.921 [95% CI, 0.809-1.000] and 0.921 [95% CI, 0.816-1.000], respectively, and there was no significant difference (p = 1.000, DeLong’s test) in the C-indexes between the two models (Table 4). The calibration curves of models were shown in Figure 3, and slopes of mode A in the development cohort and validation cohort were 1.000 and 1.830, respectively. Slopes of mode B in the development cohort and validation cohort were 1.000 and 2.056. Finally, the nomograms of model B were shown in Figure 4.

4. Discussion

In this retrospective observational study, we made two important clinical findings. First, the HBI/h was a key risk factor in EAKI occurrence after DCD LT. Second, the prediction models combining with the HBI/h and other clinical parameters could accurately predict EAKI occurrence after DCD LT and had clinical value.

A previous study has reported that predictors for postoperative late acute renal failure (ARF) differ from early ARF and correspond to postoperative parameters [23]. Therefore, we built prediction models using preoperative and intraoperative parameters to predict postoperative early AKI.

AKI after LT is multifactorial in origin. Higher CTP score [15], longer cold ischemic time[24], longer anhepatic phase [13], longer operation time [1], and larger pRBC transfusion amount [5, 14] have been reported as high-risk factors or predictors of AKI after living-donor or cadaveric LT.

Higher MELD scores or CTP scores suggest worse pretransplant liver reservation function. And a higher MELD score may suggest worse pretransplant kidney function, either. Many studies have shown an association between preoperative MELD score and postoperative AKI [1, 4, 13]. In our study, MELD score was eliminated and CTP score was selected by LASSO regression; this might because there was no significant difference in baseline renal function between EAKI and non-EAKI patients.

Hepatic ischemia-reperfusion injury (HIRI) plays a critical role in the pathogenesis of AKI after LT [25]. HIRI is associated with a systemic inflammatory response, which may cause AKI through hemodynamic mechanisms and direct tubular cell death [26]. Thus, long warm ischemic time (WIT) and CIT will increase the probability of postoperative AKI theoretically. WIT had been considered as a key predictor of AKI following LT [7, 8]. But there was no significant difference (p=0.623) in WIT between EAKI and non-EAKI patients because WIT was strictly controlled in our center. CIT as a risk factor of AKI after LT remains controversial [7, 24]; this may be because different donor criteria or different cold storage conditions were used in different centers. In our study, there was no significant difference (p=0.883) in CIT between EAKI and non-EAKI patients in the validation cohort. So CIT as a predictor of EAKI needs to be reevaluated using a larger sample size in the future.

Systemic hemodynamic changes after LT are the cause of postoperative AKI [27]. Longer anhepatic phase, longer operation time, and larger pRBC transfusion all influence systemic hemodynamics and then influence renal blood perfusion. HBI accounts for about 25% of cardiac output [18] and HBI will change more or less after blood vessel reconstruction in LT. Our study suggested EAKI was relevant to HBI and PVF. Interestingly, after eliminating multicollinearity in LASSO regression, HBI/h was selected into prediction models, but how HBI influence renal blood flow needs to be investigated in future research.

Studies have shown that the incidence of ARF after classic orthotopic LT is significantly higher than that of ameliorative piggyback liver transplantation (APBLT) [27]. In this study, APBLT was performed in DCD LT.

This study has limitations. First, although our prediction model satisfies the minimum sample size requirement [28], the larger sample size is needed to verify the accuracy of the prediction model. Second, the models should be validated prospectively at other centers to demonstrate its applicability, however, there is no HBI data from other centers could be available in present. Third, during LT, we did not directly monitor cardiac output and kidney flow which would reflect perfusion of kidney and predict postoperative EAKI better. In addition, severe stages of EAKI were not analyzed because of the small sample size.

5. Conclusion

In this study, we demonstrated that HBI can be a reliable predictor of EAKI after DCD LT for the first time. In addition, we developed an easy and accurate prediction model including CTP score, AP, HBI/h and pRBC transfusion for EAKI after DCD LT.

Data Availability

The patients’ data used to support the findings of this study are included within the supplementary information file.

Conflicts of Interest

No conflicts of interest exist.

Authors’ Contributions

Qingpeng Liu, Feng Li, Rui Guo, Bowen Wang, and Xianliang Lu are responsible for data collection and collation. Qingpeng Liu and Ao Jiao are responsible for investigation. Ao Jiao contributed to methodology. Jialin Zhang and Xiaohang Li contributed to project administration. Jialin Zhang is responsible for resources. Jialin Zhang and Xiaohang Li contributed to supervision. Ao Jiao contributed to visualization. Ao Jiao contributed to writing of the paper. Chengshuo Zhang and Ning Sun contributed to proofreading.

Acknowledgments

This work was supported by Liaoning Province Key R&D Steering Program (no. 2017225031).

Supplementary Materials

The csv file is the anonymous data set of patients. (Supplementary Materials)