# Nomograms Predicting the Occurrence of Sepsis in Patients following Major Hepatobiliary and Pancreatic Surgery

Authors: Haoyun Zhang; Fanyu Meng; Shichun Lu
Journal: Gastroenterology Research and Practice
 (2020)

---

## Abstract

Purpose. Sepsis is a severe complication in patients following major hepatobiliary and pancreatic surgery. The purpose of this study was to develop and validate a nomogram based on inflammation biomarkers and clinical characteristics. Methods. Patients who underwent major hepatobiliary and pancreatic surgery between June 2015 and April 2017 were retrospectively collected. Multivariate logistic regression was used to identify the independent risk factors associated with postoperative sepsis. A training cohort of 522 patients in an earlier period was used to develop the prediction models, and a validation cohort of 136 patients thereafter was used to validate the nomograms. Results. Sepsis developed in 55 of 522 patients of the training cohort and 19 of 136 patients in the validation cohort, respectively. In the training cohort, one nomogram based on clinical characteristics was developed. The clinical independent risk factors for postoperative sepsis include perioperative blood transfusion, diabetes, operative time, direct bilirubin, and BMI. Another nomogram was based on both clinical characteristics and inflammation biomarkers. Multivariate regression analyses showed that previous clinical risk factors, PCT, and CRP were independent risk factors for postoperative sepsis. The last nomogram showed a good C-index of 0.844 (95% CI, 0.787-0.900) compared with the previous one of 0.777 (95% CI, 0.713-0.840). Patients with a total score more than 109 in the second model are at high risk. The positive predictive value and negative predictive value of the second nomogram were 27% and 97%, respectively. Conclusion. The nomogram achieved good performances for predicting postoperative sepsis in patients by combining clinical and inflammation risk factors. This model can provide the early risk estimation of sepsis for patients following major hepatobiliary and pancreatic surgery.


---

## Body

## 1. Introduction

Severe sepsis and organ failure are major causes of morbidity and mortality after major hepatobiliary and pancreatic surgery [1]. Aggressive hepatobiliary pancreatic surgery has been associated with high complication rates of 40%–50% [1]. In a study of 583 patients who underwent surgical resection for intrahepatic cholangiocarcinoma, a postoperative complication occurred in 259 patients (44%) [2]. Moreover, sepsis is not uncommon. In one study of 1933 patients who underwent PD (pancreaticoduodenectomy), sepsis was found in 482 patients (24.9%) [3]. With advances in monitoring and prompt initiation of therapy, the mortality of severe sepsis remains higher than 25%~30%, even reaching 40-50% in the presence of shock [4]. These patients who survive to hospital discharge after sepsis are still at risk for death in the following time [5]. Those who survive often suffer from impaired physical or neurocognitive functions, emotional disorders, and a poor quality of life [6].Sepsis was once defined as systemic inflammation caused by infection [7], and then, international consensus guidelines provide a series of warning signs of early sepsis [8]. Recently, it was defined as life-threatening organ dysfunction caused by dysregulated host response to infection [9]. Improper or delayed antibiotic treatment for sepsis patients can increase mortality [10]. Therefore, intravenous antibiotic therapy should be started as early as possible and should cover all likely pathogens [5]. This situation makes it meaningful to early predict the occurrence of sepsis.In this study, we aim to develop and validate nomograms to predict sepsis of patients following major hepatobiliary and pancreatic surgery. Both serum biomarkers and clinical characteristics were analyzed by logistic regression analysis to identify risk factors. Then, one nomogram based on simple clinical characteristics and another one combining serum biomarkers and clinical characteristics were developed in a training cohort that comprises 522 patients. Finally, the nomograms were validated by a validation cohort containing 136 patients.


## 2. Materials and Methods

### 2.1. Patients

Patients who underwent major hepatobiliary and pancreatic surgery from June 2015 to April 2017 in the Department of Hepatobiliary Surgery of the First Medical Center of Chinese People’s Liberation Army General Hospital were retrospectively studied. This study was approved by the Institutional Ethics Committee of Chinese People’s Liberation Army General Hospital, and patients’ informed consent was obtained. Patients in training and validation groups come from a 4 : 1 split of all patients by time and then were divided into the training and validation groups, respectively. Blood samples were collected 1 day after surgery and tested for procalcitonin (PCT) and C-reactive protein (CRP), cytokines, and peripheral neutrophil to lymphocyte ratio (NLR).The main outcome was whether postoperative sepsis occurred in the included patients. The secondary outcomes included the time when sepsis occurred, various infection events such as pneumonia and biliary tract infection, and microbiological tests like body fluids and blood culture.Patients with signs of infection and aSOFAscore≥2 points were identified as sepsis [9]. In clinical use, patients with suspected and culture-proven infection are both in consideration. The SOFA comprises a score relevant to the major organ systems and graded from 0 to 4 according to the degree of dysfunction or failure [11].Major surgery which was defined as organ removal or normal surgical anatomy which has changed after peritoneal access as previously reported, including operations creating any gastrointestinal anastomosis or involving parenchymal resection of the liver, bile duct, spleen, or pancreas, were included [12, 13]. In this study, we mainly included hepatectomy, pancreaticoduodenectomy, hilar cholangiocarcinoma radical surgery, and other hepatobiliary operations, and less extensive surgery such as cholecystectomy and liver tumor radiofrequency ablation was excluded. But surgery including both hepatectomy and cholecystectomy was included.The exclusion criteria were patients withage<18, with liver radiofrequency ablation only, with laparoscopic cholecystectomy only, and with incomplete clinical and inflammation biomarker data.


### 2.2. Methods

In this study, we collected patients’ clinical characteristic data on gender, age, body mass index (BMI), and diabetes. Perioperative data included durations of anesthesia, operative time, ASA physical status score, hemorrhage, and blood transfusion volumes. Perioperative blood transfusion in this study was based on data from 3 days after surgery and before.Laboratory tests include direct bilirubin (DBil), aspartate aminotransferase (AST), alanine aminotransferase (ALT), lymphocyte counts, neutrophil counts, white blood cell count (WBC), C-reactive protein (CRP), procalcitonin (PCT), results of blood bacterial culture, and identified bacteria. Serum biomarkers (IL-1, IL-2, IL-6, IL-8, IL-10, and TNF-α) and peripheral neutrophil to lymphocyte ratio (NLR) were detected at 1 day after surgery.Two surgeons separately collected clinical and inflammatory biomarker data of all patients who met the inclusion criteria. For inconsistencies in data, final results were confirmed by carefully reviewing the electronic medical records.Serum IL-1, IL-2, IL-6, IL-8, IL-10, and TNF-α were measured by chemiluminescent immunoassay technology on the IMMULITE 1000 Immunoassay System as previously reported [14]. IL-6 was measured by the IMMULITE 2000 Immunoassay System.


### 2.3. Statistical Analysis

Normally distributed data of continuous variables were expressed as themean±standarddeviation and compared by an unpaired, 2-tailed t-test. Otherwise, the data were represented as the median (P25 and P75) and compared by the Mann-Whitney test. Categorical variables were compared by the chi-squared test or Fisher exact test. Univariate logistic regression was performed in the training cohort to identify variables that were associated with postoperative sepsis, and then, multivariate logistic regression analysis was used to determine whether the variables included were independent.Nomograms were built on multivariate logistic regression results using the “rms” package of R. One nomogram based on clinical characteristics and another one based on both clinical characteristics and inflammation biomarkers were developed. The predictive performance of the 2 nomograms in the training cohort was measured byC-index and with 1000 bootstrap samples to decrease the overfit bias. Then, the predicted probabilities of the 2 models in the validation cohort were shown as ROC curves.For clinical use, the total scores of each patient were calculated based on the nomograms. The cutoff value for total points of postoperative sepsis was calculated using the “OptimalCutpoints” package of R based on maximizing the Youden index. The net reclassification improvement (NRI) was calculated using the “PredictABEL” package of R.P<0.05 was considered to indicate statistical significance. All analyses were performed in SPSS and R, version 3.6.1.


## 2.1. Patients

Patients who underwent major hepatobiliary and pancreatic surgery from June 2015 to April 2017 in the Department of Hepatobiliary Surgery of the First Medical Center of Chinese People’s Liberation Army General Hospital were retrospectively studied. This study was approved by the Institutional Ethics Committee of Chinese People’s Liberation Army General Hospital, and patients’ informed consent was obtained. Patients in training and validation groups come from a 4 : 1 split of all patients by time and then were divided into the training and validation groups, respectively. Blood samples were collected 1 day after surgery and tested for procalcitonin (PCT) and C-reactive protein (CRP), cytokines, and peripheral neutrophil to lymphocyte ratio (NLR).The main outcome was whether postoperative sepsis occurred in the included patients. The secondary outcomes included the time when sepsis occurred, various infection events such as pneumonia and biliary tract infection, and microbiological tests like body fluids and blood culture.Patients with signs of infection and aSOFAscore≥2 points were identified as sepsis [9]. In clinical use, patients with suspected and culture-proven infection are both in consideration. The SOFA comprises a score relevant to the major organ systems and graded from 0 to 4 according to the degree of dysfunction or failure [11].Major surgery which was defined as organ removal or normal surgical anatomy which has changed after peritoneal access as previously reported, including operations creating any gastrointestinal anastomosis or involving parenchymal resection of the liver, bile duct, spleen, or pancreas, were included [12, 13]. In this study, we mainly included hepatectomy, pancreaticoduodenectomy, hilar cholangiocarcinoma radical surgery, and other hepatobiliary operations, and less extensive surgery such as cholecystectomy and liver tumor radiofrequency ablation was excluded. But surgery including both hepatectomy and cholecystectomy was included.The exclusion criteria were patients withage<18, with liver radiofrequency ablation only, with laparoscopic cholecystectomy only, and with incomplete clinical and inflammation biomarker data.


## 2.2. Methods

In this study, we collected patients’ clinical characteristic data on gender, age, body mass index (BMI), and diabetes. Perioperative data included durations of anesthesia, operative time, ASA physical status score, hemorrhage, and blood transfusion volumes. Perioperative blood transfusion in this study was based on data from 3 days after surgery and before.Laboratory tests include direct bilirubin (DBil), aspartate aminotransferase (AST), alanine aminotransferase (ALT), lymphocyte counts, neutrophil counts, white blood cell count (WBC), C-reactive protein (CRP), procalcitonin (PCT), results of blood bacterial culture, and identified bacteria. Serum biomarkers (IL-1, IL-2, IL-6, IL-8, IL-10, and TNF-α) and peripheral neutrophil to lymphocyte ratio (NLR) were detected at 1 day after surgery.Two surgeons separately collected clinical and inflammatory biomarker data of all patients who met the inclusion criteria. For inconsistencies in data, final results were confirmed by carefully reviewing the electronic medical records.Serum IL-1, IL-2, IL-6, IL-8, IL-10, and TNF-α were measured by chemiluminescent immunoassay technology on the IMMULITE 1000 Immunoassay System as previously reported [14]. IL-6 was measured by the IMMULITE 2000 Immunoassay System.


## 2.3. Statistical Analysis

Normally distributed data of continuous variables were expressed as themean±standarddeviation and compared by an unpaired, 2-tailed t-test. Otherwise, the data were represented as the median (P25 and P75) and compared by the Mann-Whitney test. Categorical variables were compared by the chi-squared test or Fisher exact test. Univariate logistic regression was performed in the training cohort to identify variables that were associated with postoperative sepsis, and then, multivariate logistic regression analysis was used to determine whether the variables included were independent.Nomograms were built on multivariate logistic regression results using the “rms” package of R. One nomogram based on clinical characteristics and another one based on both clinical characteristics and inflammation biomarkers were developed. The predictive performance of the 2 nomograms in the training cohort was measured byC-index and with 1000 bootstrap samples to decrease the overfit bias. Then, the predicted probabilities of the 2 models in the validation cohort were shown as ROC curves.For clinical use, the total scores of each patient were calculated based on the nomograms. The cutoff value for total points of postoperative sepsis was calculated using the “OptimalCutpoints” package of R based on maximizing the Youden index. The net reclassification improvement (NRI) was calculated using the “PredictABEL” package of R.P<0.05 was considered to indicate statistical significance. All analyses were performed in SPSS and R, version 3.6.1.


## 3. Results

### 3.1. Clinical Characteristics and Inflammation Biomarkers

A total of 658 patients who underwent major hepatobiliary and pancreatic surgery and met the inclusion criteria were included. 522 patients and 136 patients were divided into the training and validation groups, respectively.The inflammation biomarkers were detected at 1 day after surgery. Patients’ clinical characteristics in the training and validation cohorts were given in Table1. There were no significant differences between the 2 cohorts in sepsis prevalence (P=0.26). The finally diagnosed sepsis was found in 55 and 19 patients in the training and validation cohorts, respectively. Tables 2 and 3 show the characteristics of sepsis and nonseptic patients in the training and validation cohorts, respectively.Table 1
Characteristics of patients in the training and validation cohorts.
VariablesCohortP valueTraining (n=522)Validation (n=136)Age55.82±12.9156.60±11.110.52GenderMale339970.16Female18339DiabetesYes95200.34No427116BMI23.62±3.2223.55±3.730.83Perioperative blood transfusionYes203450.21No31991Preoperative ALT (U/l)45.05 (21.52-121)40.5 (20.2-117.2)0.47Preoperative AST (U/l)37.55 (21.2-96.75)36.7 (19.57-79.45)0.48Preoperative direct bilirubin (μmol/l)5.4 (3.3-38.8)5.25 (3.2-29.43)0.83Intraoperative blood loss (ml)300 (200-600)300 (200-500)0.51Operative time (min)285 (200-345)285 (200-365)0.75ASA score11140.14243410737522422501Table 2
Clinical characteristics and inflammation biomarkers in the training cohort. Data are shown as themean±SD or median (P25 and P75).
VariableTraining cohortSepsis (-),n=467Sepsis (+),n=55Age55.41±13.0559.29±11.240.02∗SexMale305340.61Female16221Perioperative blood transfusionYes16538≤0.001∗No30217ASA scores1101≤0.001∗23924236411411DiabetesYes78170.01∗No38938BMI23.5±3.2124.69±3.120.01∗Operative time (min)280 (210-340)345 (262.5-431)0.03∗Intraoperative blood loss300 (200-600)400 (200-800)Preoperative ALT44 (20.45-122.5)62 (24.5-104.9)0.53Preoperative AST35.6 (20.85-98.95)42.6 (26.4-81.9)0.47Preoperative direct bilirubin5.1 (3.2-31)35.5 (5.1-132.05)≤0.001∗Inflammation biomarkersIL-1 (pg/ml)5 (5-5)5 (5-5)0.46IL-2 (U/ml)767 (600-1023)1034 (699.5-1378.5)≤0.001∗IL-6 (pg/ml)82.2 (46.15-164.5)171 (102-274)≤0.001∗IL-8 (pg/ml)89 (40.2-190)105 (55.55-242)0.18IL-10 (pg/ml)5.78 (5-10)9.56 (5.86-14.35)≤0.001∗TNF-α (pg/ml)16.3 (11-30.55)16.5 (10.5-24.95)0.50PCT (ng/ml)0.55 (0.28-1.01)1.53 (0.809-3.29)≤0.001∗CRP (mg/dl)3.73 (1.96-6.29)6.62 (4.62-9.35)≤0.001∗NLR11.1 (7.73-15.73)12.17 (8.10-16.94)0.39Table 3
Clinical characteristics and inflammation biomarkers of patients in the validation cohort. Data are shown as themean±SD or median (P25 and P75).
VariableValidation cohortSepsis (-),n=117Sepsis (+),n=19Age55.8±11.0761.53±10.300.04∗SexMale83140.81Female345Perioperative blood transfusionYes3312≤0.001∗No847ASA scores1400.04∗295123175411501DiabetesYes1550.16No10214BMI23.47±3.7724.02±3.540.56Operative time265 (190-345)395 (330-485)≤0.001∗Intraoperative blood loss300 (150-500)300 (300-500)0.10Preoperative ALT39.5 (19.6-116.6)57 (26.9-117.5)0.52Preoperative AST35 (20-79.1)47.6 (17.9-104.85)0.77Preoperative direct bilirubin4.6 (3.2-14.9)14.8 (5.05-144.75)0.03∗Inflammation biomarkersIL-1 (pg/ml)5 (5-5)5 (5-5)0.31IL-2 (U/ml)805 (591-1053)1226 (978-1963)≤0.001∗IL-6 (pg/ml)64.6 (37.5-154)151 (97.75-350.5)0.01∗IL-8 (pg/ml)96 (36.9-187)112 (81.95-231.5)0.27IL-10 (pg/ml)5 (5-5)6.57 (5-11.5)0.25TNF-α (pg/ml)14 (10.2-25.8)22.2 (17.8-35.25)0.02∗PCT (ng/ml)0.544 (0.26-0.911)1.23 (0.89-2.43)≤0.001∗CRP (mg/dl)3.28 (1.75-5.18)5.93 (3.88-13.1)≤0.001∗NLR10.96 (7.87-16.96)9.51 (4.72-12.49)0.18∗P<0.05.All patients in this study underwent postoperative (one day postsurgery) blood test for PCT, CRP, and cytokines (IL-1, IL-2, IL-6, IL-8, IL-10, and TNF-α). NLR was also collected. The median levels and interquartile ranges of inflammation biomarkers are listed in Table 4. There were no significant differences among the 2 cohorts.Table 4
The median levels and interquartile ranges of inflammation biomarkers in the training and validation cohorts.
VariablesCohortP valueTraining (n=522)Validation (n=136)IL-1 (pg/ml)5 (5-5)5 (5-5)0.75IL-2 (U/ml)788 (607.5-1074.8)851 (625-1106)0.33IL-6 (pg/ml)90.55 (48.12-178.5)70.7 (39.52-191.25)0.20IL-8 (pg/ml)91.15 (41.95-194.5)99 (38.08-188)0.90IL-10 (pg/ml)6.02 (5-10.38)6.33 (5-10.4)0.66TNF-α (pg/ml)16.35 (11-29.9)15.75 (10.3-28.45)0.41PCT (ng/ml)0.60 (0.32-1.14)0.57 (0.31-1.10)0.90CRP (mg/dl)4.02 (2.14-6.71)3.65 (2.03-5.41)0.13NLR11.22 (7.76-15.85)10.88 (7.62-16.66)0.63


### 3.2. Development and Validation of Nomograms

The results of univariate logistic regression analysis are presented in Table5. The results were reported as odds ratio (95% CI). Two different multivariate logistic regression analyses were performed. The first one only included clinical characteristics as shown in Table 6, and the second one included both clinical characteristics and inflammation biomarkers as shown in Table 7. In the first regression results, perioperative blood transfusion (OR: 3.638 (1.921-6.891)), diabetes (OR: 2.378 (1.221-4.633)), operative time (OR: 1.003 (1.001-1.005)), direct bilirubin (OR: 1.004 (1.001-1.007)), and BMI (OR: 1.164 (1.057-1.281)) were independently associated with sepsis. In the second regression results, these clinical characteristics plus PCT (OR: 1.143 (1.071-1.221)) and CRP (OR: 1.175 (1.081-1.277)) were independently associated with sepsis.Table 5
Univariate regression analysis of sepsis based on clinical characteristics and inflammation biomarkers in patients following major surgery.
VariablesORPIL-1 (pg/ml)0.997 (0.963-1.032)0.862IL-2 (U/ml)1.001 (1.000-1.001)0.001∗IL-6 (pg/ml)1.003 (1.001-1.004)≤0.001∗IL-8 (pg/ml)1.000 (1.000-1.001)0.125IL-10 (pg/ml)1.005 (1.001-1.009)0.016∗TNF-α (pg/ml)1.001 (0.994-1.008)0.697PCT (ng/ml)1.134 (1.066-1.205)≤0.001∗CRP (mg/dl)1.209 (1.125-1.300)≤0.001∗NLR1.030 (1.002-1.057)0.033∗Perioperative blood transfusion (yes vs. no)4.091 (2.240-7.474)≤0.001∗ASA scores1.715 (0.920-3.197)0.090Diabetes (yes vs. no)2.231 (1.198-4.153)0.011∗BMI1.12 (1.027-1.221)0.01∗Age1.026 (1.002-1.050)0.036∗Gender (male vs. female)0.86 (0.483-1.53)0.608Operative time1.003 (1.002-1.005)≤0.001∗Intraoperative blood loss1.000 (1.000-1.001)0.013∗ALT (U/l)1.001 (0.999-1.002)0.278AST (U/l)1.001 (1.000-1.002)0.17Direct bilirubin (μmol/l)1.005 (1.002-1.008)0.001∗ASA: American Society of Anesthesiologists.∗P<0.05.Table 6
Multivariate regression analysis of sepsis based on clinical characteristics in patients following major surgery.
VariableBPOR (95% CI)Blood transfusion (yes vs. no)1.291≤0.0013.638 (1.921-6.891)Diabetes (yes vs. no)0.8660.0112.378 (1.221-4.633)Operative time (min)0.0030.0121.003 (1.001-1.005)Preoperative direct bilirubin (μmol/l)0.0040.0101.004 (1.001-1.007)BMI0.1510.0021.164 (1.057-1.281)Table 7
Multivariate regression analysis of sepsis based on clinical characteristics and inflammation biomarkers in patients following major surgery.
VariableBPOR (95% CI)Blood transfusion (yes vs. no)1.276≤0.0013.583 (1.815-7.072)Diabetes (yes vs. no)0.9100.0122.485 (1.224-5.046)Operative time (min)0.0020.0231.002 (1.000-1.005)Preoperative direct bilirubin (μmol/l)0.0040.0191.004 (1.001-1.007)BMI0.1770.0011.193 (1.075-1.325)PCT (ng/ml)0.134≤0.0011.143 (1.071-1.221)CRP (mg/dl)0.162≤0.0011.175 (1.081-1.277)Sepsis risk estimation nomograms based on the 2 multivariate regression analysis results were developed as shown in Figure1. Both nomograms showed good accuracy in predicting sepsis. The 2 nomograms were validated using the bootstrap method in the training cohort. The first one has a C-index of 0.777 (95% CI, 0.713-0.840) and a bootstrap-corrected C-index of 0.761, and the second one has a C-index of 0.844 (95% CI, 0.787-0.900) and a bootstrap-corrected C-index of 0.831.Figure 1
(a) Nomogram based on clinical variables for predicting sepsis in patients following major hepatobiliary and pancreatic surgery. (b) Nomogram based on clinical and inflammation variables for predicting sepsis in patients following major hepatobiliary and pancreatic surgery. (c) Bootstrap validation of the predictive performance of the first nomogram in estimating the risk of sepsis in the training cohort. (d) Bootstrap validation of the predictive performance of the second nomogram in estimating the risk of sepsis in the training cohort.
(a)(b)(c)(d)


### 3.3. Compare the Performance of Nomograms for Predicting Sepsis in the Validation Cohort

The predicted sepsis probabilities of the 2 different nomograms for patients in the validation cohort were calculated. The diagnostic performances of the 2 different models for the validation cohort were evaluated by ROC curves and shown in Figure2. In the validation cohort, the first nomogram showed an AUC of 0.756 (95% CI, 0.647-0.864); the second nomogram showed an AUC of 0.839 (95% CI, 0.745-0.932). There is a significant difference between these two ROC curves, with a P value of 0.048. Besides, we could see that the AUC of model 2 was larger.Figure 2
“Model 1” showed the diagnostic performance of predicted probabilities in the validation cohort based on clinical variables; “model 2” showed the diagnostic performance of predicted probabilities in the validation cohort based on clinical variables and inflammation biomarkers.


### 3.4. Cutoff Value of the Predicting Risk for Estimating Sepsis

To better use these nomograms for patients, the best cutoff total scores for risk estimation of the 2 nomograms were calculated based on maximizing the Youden index. We found that at cutoff total scores of 106 for the first model and 109 for the second model, respectively, these two prediction models have the largest Youden index as shown in Table8. Using 106 for the first model and 109 for the second model as our cutoff values for high risk, the positive predictive value was 33% (validation cohort) and the negative predictive value was 91% (validation cohort) for model 1; the positive predictive value was 35% (validation cohort) and the negative predictive value was 95% (validation cohort) for model 2. The positive likelihood ratio and negative likelihood ratio in the training or validation cohort for both models are also shown in Table 8.Table 8
Diagnostic performances of the 2 nomograms for postoperative sepsis estimation.
VariablesTraining cohortValidation cohortModel 1 (clinical variable)Cutoff score106106Sensitivity0.640.47Specificity0.820.85Positive predictive value0.290.33Negative predictive value0.950.91Positive likelihood ratio3.503.13Negative likelihood ratio0.440.62Model 2 (clinical and inflammation variables)Cutoff score109109Sensitivity0.840.74Specificity0.730.78Positive predictive value0.270.35Negative predictive value0.970.95Positive likelihood ratio3.153.36Negative likelihood ratio0.220.33To demonstrate the superiority of PCT and CRP in model 2, we then compare the performance of the 2 models for the validation cohort, and the net reclassification improvement (NRI) was 0.1125 (P=0.43).


## 3.1. Clinical Characteristics and Inflammation Biomarkers

A total of 658 patients who underwent major hepatobiliary and pancreatic surgery and met the inclusion criteria were included. 522 patients and 136 patients were divided into the training and validation groups, respectively.The inflammation biomarkers were detected at 1 day after surgery. Patients’ clinical characteristics in the training and validation cohorts were given in Table1. There were no significant differences between the 2 cohorts in sepsis prevalence (P=0.26). The finally diagnosed sepsis was found in 55 and 19 patients in the training and validation cohorts, respectively. Tables 2 and 3 show the characteristics of sepsis and nonseptic patients in the training and validation cohorts, respectively.Table 1
Characteristics of patients in the training and validation cohorts.
VariablesCohortP valueTraining (n=522)Validation (n=136)Age55.82±12.9156.60±11.110.52GenderMale339970.16Female18339DiabetesYes95200.34No427116BMI23.62±3.2223.55±3.730.83Perioperative blood transfusionYes203450.21No31991Preoperative ALT (U/l)45.05 (21.52-121)40.5 (20.2-117.2)0.47Preoperative AST (U/l)37.55 (21.2-96.75)36.7 (19.57-79.45)0.48Preoperative direct bilirubin (μmol/l)5.4 (3.3-38.8)5.25 (3.2-29.43)0.83Intraoperative blood loss (ml)300 (200-600)300 (200-500)0.51Operative time (min)285 (200-345)285 (200-365)0.75ASA score11140.14243410737522422501Table 2
Clinical characteristics and inflammation biomarkers in the training cohort. Data are shown as themean±SD or median (P25 and P75).
VariableTraining cohortSepsis (-),n=467Sepsis (+),n=55Age55.41±13.0559.29±11.240.02∗SexMale305340.61Female16221Perioperative blood transfusionYes16538≤0.001∗No30217ASA scores1101≤0.001∗23924236411411DiabetesYes78170.01∗No38938BMI23.5±3.2124.69±3.120.01∗Operative time (min)280 (210-340)345 (262.5-431)0.03∗Intraoperative blood loss300 (200-600)400 (200-800)Preoperative ALT44 (20.45-122.5)62 (24.5-104.9)0.53Preoperative AST35.6 (20.85-98.95)42.6 (26.4-81.9)0.47Preoperative direct bilirubin5.1 (3.2-31)35.5 (5.1-132.05)≤0.001∗Inflammation biomarkersIL-1 (pg/ml)5 (5-5)5 (5-5)0.46IL-2 (U/ml)767 (600-1023)1034 (699.5-1378.5)≤0.001∗IL-6 (pg/ml)82.2 (46.15-164.5)171 (102-274)≤0.001∗IL-8 (pg/ml)89 (40.2-190)105 (55.55-242)0.18IL-10 (pg/ml)5.78 (5-10)9.56 (5.86-14.35)≤0.001∗TNF-α (pg/ml)16.3 (11-30.55)16.5 (10.5-24.95)0.50PCT (ng/ml)0.55 (0.28-1.01)1.53 (0.809-3.29)≤0.001∗CRP (mg/dl)3.73 (1.96-6.29)6.62 (4.62-9.35)≤0.001∗NLR11.1 (7.73-15.73)12.17 (8.10-16.94)0.39Table 3
Clinical characteristics and inflammation biomarkers of patients in the validation cohort. Data are shown as themean±SD or median (P25 and P75).
VariableValidation cohortSepsis (-),n=117Sepsis (+),n=19Age55.8±11.0761.53±10.300.04∗SexMale83140.81Female345Perioperative blood transfusionYes3312≤0.001∗No847ASA scores1400.04∗295123175411501DiabetesYes1550.16No10214BMI23.47±3.7724.02±3.540.56Operative time265 (190-345)395 (330-485)≤0.001∗Intraoperative blood loss300 (150-500)300 (300-500)0.10Preoperative ALT39.5 (19.6-116.6)57 (26.9-117.5)0.52Preoperative AST35 (20-79.1)47.6 (17.9-104.85)0.77Preoperative direct bilirubin4.6 (3.2-14.9)14.8 (5.05-144.75)0.03∗Inflammation biomarkersIL-1 (pg/ml)5 (5-5)5 (5-5)0.31IL-2 (U/ml)805 (591-1053)1226 (978-1963)≤0.001∗IL-6 (pg/ml)64.6 (37.5-154)151 (97.75-350.5)0.01∗IL-8 (pg/ml)96 (36.9-187)112 (81.95-231.5)0.27IL-10 (pg/ml)5 (5-5)6.57 (5-11.5)0.25TNF-α (pg/ml)14 (10.2-25.8)22.2 (17.8-35.25)0.02∗PCT (ng/ml)0.544 (0.26-0.911)1.23 (0.89-2.43)≤0.001∗CRP (mg/dl)3.28 (1.75-5.18)5.93 (3.88-13.1)≤0.001∗NLR10.96 (7.87-16.96)9.51 (4.72-12.49)0.18∗P<0.05.All patients in this study underwent postoperative (one day postsurgery) blood test for PCT, CRP, and cytokines (IL-1, IL-2, IL-6, IL-8, IL-10, and TNF-α). NLR was also collected. The median levels and interquartile ranges of inflammation biomarkers are listed in Table 4. There were no significant differences among the 2 cohorts.Table 4
The median levels and interquartile ranges of inflammation biomarkers in the training and validation cohorts.
VariablesCohortP valueTraining (n=522)Validation (n=136)IL-1 (pg/ml)5 (5-5)5 (5-5)0.75IL-2 (U/ml)788 (607.5-1074.8)851 (625-1106)0.33IL-6 (pg/ml)90.55 (48.12-178.5)70.7 (39.52-191.25)0.20IL-8 (pg/ml)91.15 (41.95-194.5)99 (38.08-188)0.90IL-10 (pg/ml)6.02 (5-10.38)6.33 (5-10.4)0.66TNF-α (pg/ml)16.35 (11-29.9)15.75 (10.3-28.45)0.41PCT (ng/ml)0.60 (0.32-1.14)0.57 (0.31-1.10)0.90CRP (mg/dl)4.02 (2.14-6.71)3.65 (2.03-5.41)0.13NLR11.22 (7.76-15.85)10.88 (7.62-16.66)0.63


## 3.2. Development and Validation of Nomograms

The results of univariate logistic regression analysis are presented in Table5. The results were reported as odds ratio (95% CI). Two different multivariate logistic regression analyses were performed. The first one only included clinical characteristics as shown in Table 6, and the second one included both clinical characteristics and inflammation biomarkers as shown in Table 7. In the first regression results, perioperative blood transfusion (OR: 3.638 (1.921-6.891)), diabetes (OR: 2.378 (1.221-4.633)), operative time (OR: 1.003 (1.001-1.005)), direct bilirubin (OR: 1.004 (1.001-1.007)), and BMI (OR: 1.164 (1.057-1.281)) were independently associated with sepsis. In the second regression results, these clinical characteristics plus PCT (OR: 1.143 (1.071-1.221)) and CRP (OR: 1.175 (1.081-1.277)) were independently associated with sepsis.Table 5
Univariate regression analysis of sepsis based on clinical characteristics and inflammation biomarkers in patients following major surgery.
VariablesORPIL-1 (pg/ml)0.997 (0.963-1.032)0.862IL-2 (U/ml)1.001 (1.000-1.001)0.001∗IL-6 (pg/ml)1.003 (1.001-1.004)≤0.001∗IL-8 (pg/ml)1.000 (1.000-1.001)0.125IL-10 (pg/ml)1.005 (1.001-1.009)0.016∗TNF-α (pg/ml)1.001 (0.994-1.008)0.697PCT (ng/ml)1.134 (1.066-1.205)≤0.001∗CRP (mg/dl)1.209 (1.125-1.300)≤0.001∗NLR1.030 (1.002-1.057)0.033∗Perioperative blood transfusion (yes vs. no)4.091 (2.240-7.474)≤0.001∗ASA scores1.715 (0.920-3.197)0.090Diabetes (yes vs. no)2.231 (1.198-4.153)0.011∗BMI1.12 (1.027-1.221)0.01∗Age1.026 (1.002-1.050)0.036∗Gender (male vs. female)0.86 (0.483-1.53)0.608Operative time1.003 (1.002-1.005)≤0.001∗Intraoperative blood loss1.000 (1.000-1.001)0.013∗ALT (U/l)1.001 (0.999-1.002)0.278AST (U/l)1.001 (1.000-1.002)0.17Direct bilirubin (μmol/l)1.005 (1.002-1.008)0.001∗ASA: American Society of Anesthesiologists.∗P<0.05.Table 6
Multivariate regression analysis of sepsis based on clinical characteristics in patients following major surgery.
VariableBPOR (95% CI)Blood transfusion (yes vs. no)1.291≤0.0013.638 (1.921-6.891)Diabetes (yes vs. no)0.8660.0112.378 (1.221-4.633)Operative time (min)0.0030.0121.003 (1.001-1.005)Preoperative direct bilirubin (μmol/l)0.0040.0101.004 (1.001-1.007)BMI0.1510.0021.164 (1.057-1.281)Table 7
Multivariate regression analysis of sepsis based on clinical characteristics and inflammation biomarkers in patients following major surgery.
VariableBPOR (95% CI)Blood transfusion (yes vs. no)1.276≤0.0013.583 (1.815-7.072)Diabetes (yes vs. no)0.9100.0122.485 (1.224-5.046)Operative time (min)0.0020.0231.002 (1.000-1.005)Preoperative direct bilirubin (μmol/l)0.0040.0191.004 (1.001-1.007)BMI0.1770.0011.193 (1.075-1.325)PCT (ng/ml)0.134≤0.0011.143 (1.071-1.221)CRP (mg/dl)0.162≤0.0011.175 (1.081-1.277)Sepsis risk estimation nomograms based on the 2 multivariate regression analysis results were developed as shown in Figure1. Both nomograms showed good accuracy in predicting sepsis. The 2 nomograms were validated using the bootstrap method in the training cohort. The first one has a C-index of 0.777 (95% CI, 0.713-0.840) and a bootstrap-corrected C-index of 0.761, and the second one has a C-index of 0.844 (95% CI, 0.787-0.900) and a bootstrap-corrected C-index of 0.831.Figure 1
(a) Nomogram based on clinical variables for predicting sepsis in patients following major hepatobiliary and pancreatic surgery. (b) Nomogram based on clinical and inflammation variables for predicting sepsis in patients following major hepatobiliary and pancreatic surgery. (c) Bootstrap validation of the predictive performance of the first nomogram in estimating the risk of sepsis in the training cohort. (d) Bootstrap validation of the predictive performance of the second nomogram in estimating the risk of sepsis in the training cohort.
(a)(b)(c)(d)


## 3.3. Compare the Performance of Nomograms for Predicting Sepsis in the Validation Cohort

The predicted sepsis probabilities of the 2 different nomograms for patients in the validation cohort were calculated. The diagnostic performances of the 2 different models for the validation cohort were evaluated by ROC curves and shown in Figure2. In the validation cohort, the first nomogram showed an AUC of 0.756 (95% CI, 0.647-0.864); the second nomogram showed an AUC of 0.839 (95% CI, 0.745-0.932). There is a significant difference between these two ROC curves, with a P value of 0.048. Besides, we could see that the AUC of model 2 was larger.Figure 2
“Model 1” showed the diagnostic performance of predicted probabilities in the validation cohort based on clinical variables; “model 2” showed the diagnostic performance of predicted probabilities in the validation cohort based on clinical variables and inflammation biomarkers.


## 3.4. Cutoff Value of the Predicting Risk for Estimating Sepsis

To better use these nomograms for patients, the best cutoff total scores for risk estimation of the 2 nomograms were calculated based on maximizing the Youden index. We found that at cutoff total scores of 106 for the first model and 109 for the second model, respectively, these two prediction models have the largest Youden index as shown in Table8. Using 106 for the first model and 109 for the second model as our cutoff values for high risk, the positive predictive value was 33% (validation cohort) and the negative predictive value was 91% (validation cohort) for model 1; the positive predictive value was 35% (validation cohort) and the negative predictive value was 95% (validation cohort) for model 2. The positive likelihood ratio and negative likelihood ratio in the training or validation cohort for both models are also shown in Table 8.Table 8
Diagnostic performances of the 2 nomograms for postoperative sepsis estimation.
VariablesTraining cohortValidation cohortModel 1 (clinical variable)Cutoff score106106Sensitivity0.640.47Specificity0.820.85Positive predictive value0.290.33Negative predictive value0.950.91Positive likelihood ratio3.503.13Negative likelihood ratio0.440.62Model 2 (clinical and inflammation variables)Cutoff score109109Sensitivity0.840.74Specificity0.730.78Positive predictive value0.270.35Negative predictive value0.970.95Positive likelihood ratio3.153.36Negative likelihood ratio0.220.33To demonstrate the superiority of PCT and CRP in model 2, we then compare the performance of the 2 models for the validation cohort, and the net reclassification improvement (NRI) was 0.1125 (P=0.43).


## 4. Discussion

In this study, we developed and validated 2 different nomograms for early prediction of sepsis in patients following major hepatobiliary and pancreatic surgery. The last nomogram incorporates 2 inflammation biomarkers including the serum level of PCT and CRP. Incorporating the inflammation biomarkers and clinical risk factors into nomograms facilitates the early prediction of sepsis.Similar to our results, a meta-analysis of over 30 million patients identified that the risk factors associated with postoperative sepsis also include perioperative blood transfusion (OR: 1.90) and diabetes (OR: 1.41) [15]. Other risk factors reported include male gender and emergency surgery [16]. Most of these studies were based on patient and surgery-related risk factors for postoperative sepsis. In our study, inflammation biomarkers like cytokines (IL-1, IL-2, IL-6, IL-8, IL-10, and TNF-α), PCT, CRP, and NLR were studied. IL-2, IL-6, IL-10, and NLR were significant after univariate regression analysis, but they were not independent risk factors for sepsis in our study.PCT and CRP are both widely used in diagnosis of sepsis [17]. PCT is a peptide released in response to proinflammatory stimuli, particularly bacterial-related inflammatory mediators [18], although there is disagreement on the accuracy of PCT for differentiating sepsis from other noninfectious causes of SIRS (systemic inflammatory response syndrome) [19]. In this study, we found that PCT (P≤0.001, 1.143 (1.071-1.221)) was an independent risk factor for the occurrence of sepsis among patients who underwent major hepatobiliary and pancreatic surgery. No better markers are available that outperform PCT in diagnosis of sepsis [20, 21]. CRP showed significant differences in both training (P≤0.001) and validation (P≤0.001) cohorts. In multivariate analysis, CRP (P≤0.001, OR: 1.175 (1.081-1.277)) served as an independent risk factor. It was reported that an elevated serum CRP level is correlated with increased risk of organ failure and death [22].Regarding clinical risk factors, preoperative direct bilirubin levels are independent risk factors for sepsis. It was reported that jaundice (totalbilirubin>2.5mg/dl) was associated with 30-day mortality in patients with bacteraemic cholangitis [23]. Patients with biliary obstruction are common in hepatobiliary surgery, and much attention should be paid to these patients. BMI was also independent risk factors for sepsis in our study; BMI was reported to be associated with increased risk of infections [24]. Obese patients with BMI>30 kg/m2 should be noted. In our study, operative time was also an independent risk factor for sepsis. Duration of operation was reported to correlate with complications, and long duration procedures had greater risk of sepsis/sepsis shock [25]. Blood transfusion was also reported to be independently associated with higher odds of sepsis and septic shock and increased overall 30-day mortality [3]. In our study, blood transfusion was an independent risk factor for the occurrence of sepsis.As shown in Figure2, the last nomogram showed a bigger AUC for the validation group. And the 2 nomograms showed significant differences in predicting the sepsis risk of the validation cohort. The addition of inflammation biomarkers into the predictive model obtained satisfactory improvement for sepsis prediction (the C-index of model 1 was 0.777 in the training cohort and 0.756 in the validation cohort; the C-index of model 2 was 0.844 in the training cohort and 0.839 in the validation cohort). For clinical use of the 2 models, we evaluate the sensitivity, specificity, positive predictive value, and negative predictive value in risk estimation of sepsis using the method of Youden index. Patients with total scores more than 109 in model 2 are at a high-risk subgroup of sepsis. A high negative predictive value (97% in the training cohort) was yielded, but the positive predictive value (27% in the training cohort) was less impressive. Considering that the cost for early interventions such as the timely use of antibiotics and the improvement of bacterial culture tests is low and the patients will get the maximum benefit, we believe that the last model will allow surgeons to estimate the sepsis risk of patients who underwent major hepatobiliary and pancreatic surgery soon after surgery, and corrective treatment can be applied in time.A limitation of our study mainly includes its retrospective design and a relatively smaller sample size. Patients in training and validation groups come from a split of all patients by time in our institution rather than an independent group of patients. Therefore, more strict external validation and prospective study are needed.


## 5. Conclusion

By combining clinical and inflammatory risk factors, we constructed a nomogram for postoperative sepsis for patients following major hepatobiliary and pancreatic surgery. This model can provide the early risk estimation of sepsis for patients after major hepatobiliary and pancreatic surgery.



---

*Source: 9761878-2020-08-01.xml*