# Development an Inflammation-Related Factor-Based Model for Predicting Organ Failure in Acute Pancreatitis: A Retrospective Cohort Study

Authors: Yunpeng Peng; Xiaole Zhu; Chaoqun Hou; Chenyuan Shi; Dongya Huang; Zipeng Lu; Yi Miao; Qiang Li
Journal: Mediators of Inflammation
 (2021)

---

## Abstract

Several inflammation-related factors (IRFs) have been reported to predict organ failure of acute pancreatitis (AP) in previous clinical studies. However, there are a few shortcomings in these models. The aim of this study was to develop a new prediction model based on IRFs that could accurately identify the risk for organ failure in AP.Methods. 100 patients with their clinical information and IRF data (levels of 10 cytokines, percentages of different immune cells, and data obtained from white blood cell count) were retrospectively enrolled in this study, and 94 patients were finally selected for further analysis. Univariate and multivariate analysis were applied to evaluate the potential risk factors for the organ failure of AP. The area under the ROC curve (AUCs), sensitivity, and specificity of the relevant model were assessed to evaluate the prediction ability of IRFs. A new scoring system to predict the organ failure of AP was created based on the regression coefficient of a multivariate logistic regression model. Results. The incidence of OF in AP patients was nearly 16% (15/94) in our derivation cohort. Univariate analytic data revealed that IL6, IL8, IL10, MCP1, CD3+ CD4+ T lymphocytes, CD19+ B lymphocytes, PCT, APACHE II score, and RANSON score were potential predictors for AP organ failure, and IL6 (P=0.038), IL8 (P=0.043), and CD19+B lymphocytes (P=0.045) were independent predictors according to further multivariate analysis. In addition, a preoperative scoring system (0-11 points) was constructed to predict the organ failure of AP using these three factors. The AUC of the new score system was 0.86. The optimal cut-off value of the new scoring system was 6 points. Conclusions. Our prediction model (based on IL6, IL8, and CD19+ B Lymphocyte) has satisfactory working efficiency to identify AP patients with high risk of organ failure.


---

## Body

## 1. Introduction

Generally, acute pancreatitis (AP) commonly presents as a mild and self-limiting course with slighter clinical symptoms. However, nearly 20% AP patients present with severe disease, which is correlated with higher rates of organ failure and mortality [1, 2]. According to previous studies, organ failure typically not only develops early in the course of acute pancreatitis but also may develop later due to infected pancreatic necrosis-induced sepsis. Organ failure is the most important determinant of outcome in acute pancreatitis [3]. Therefore, early identification of potential risk factors for OF occurrence in AP patients is helpful for physicians to select accurate and effective treatment interventions.In previous studies, many factors have been reported as independent risk factors for the development of OF in AP, such as lipase, albumin, BUN, and pleural effusion. Moreover, the role of inflammation-related factors (IRFs) in secondary OF is getting more and more attention [4]. The severity of AP is significantly associated with the occurrence and development of systemic inflammatory response syndrome (SIRS). In the early stage of acute pancreatitis (AP), abundant cytokines induced by local pancreatic inflammation enter the bloodstream and further cause SIRS by “trigger effect,” which eventually leads to organ failure [5]. Different immune cells and cytokines are the key members, and abnormal alteration of some members has been confirmed as potential indicators for the progression of AP. For example, lower level of CD3+ CD4+ T lymphocytes and higher level of CD19+ B lymphocytes could predict organ failure in the early phase of AP [6], and escalating levels of cytokines like IL6, IL8, MCP-1, and TNF-a were frequently observed in AP patients with organ failure [7]. Furthermore, the most commonly used inflammation biomarkers in clinical practice, C-reactive protein (CRP) and PCT, play similar roles in the early prediction of AP-related OF [8, 9]. It has been reported that patients with early elevated levels of CRP and PCT tended to suffer with organ failure in several studies [10].In general, various types of IRFs were potential predictive indicators of OF in patients with AP. However, the forecast performance of each indicator was insufficient in clinical practice. Therefore, the purpose of this work was to develop and validate a simple and effective risk model based on IRFs for the early recognition of OF.


## 2. Patients and Methods

Consecutive patients with AP treated in the First Affiliated Hospital of Nanjing Medical University (Nanjing, Jiangsu province, China) from 2017.06 to 2018.06 were retrospectively enrolled in our study, and patients referred from other hospitals were excluded in our study. Three criteria were applied to diagnose AP, including persistent abdominal pain, serum amylase or lipase elevation of more than 3 times upper limit of normal, and characteristic ultrasound and/or CT findings. Patients were diagnosed as AP when two or more signs mentioned above were observed. Furthermore, patients with recent surgery (less than 1 week), chronic, traumatic, endoscopic, or recurrent acute pancreatitis, immunosuppression, or immune deficiency were excluded. The principles of the Helsinki Declaration were applied in the performance of our study, and this study was approved by Ethics Committee of the First Affiliated Hospital of Nanjing Medical University.According to the modified Marshall scoring system, serum creatinine was applied to define renal failure (more than 1.9 mg/dL), systolic blood pressure was used to define cardiovascular failure (less than 90 mmHg), and ratio of PaO2/FiO2 was used to confirm respiratory failure (less than 300 mmHg). The duration was used to distinguish transient OF (≤48 h) from persistent OF (>48 h) [11]. Patients who have persistent OF (>48 h) was finally enrolled into this study.


### 2.1. Data Collection

The peripheral blood samples of all patients were obtained within 24 h after admission to hospital. Peripheral blood lymphocyte subset assay was applied to detect the percentage of CD3+T lymphocytes, CD3+ CD4+T lymphocytes, CD3+ CD8+ cytotoxic T lymphocytes, CD16+ CD56+ natural killer cells, and CD19+ B lymphocytes. Blood routine examination was used to obtain white blood cell count and percentage. Blood samples were routinely collected on AP patients admitted to our hospital and then stored and subsequently analyzed for the concentration of 10 cytokines (IL1a, IL1b, IL4, IL6, IL8, IL10, IL13, MCP-1, IFNg, and TNFa) by using the Human Inflammation Array Q1 (QAH-INF-1-1) (RayBiotech, Norcross, America). Moreover, other IRF (CRP and PCT) data, several biochemical parameters, and basic clinical information were also collected.


### 2.2. Statistical Analysis

All statistical analyses were conducted by Stata/SE version 10.0 for Windows. Descriptive data were presented asmean±standarddeviation for quantitative variables or absolute and relative frequencies for qualitative variables. Student’s t-test or nonparametric tests (Chi-square and Mann–Whitney) were used to compare OF and NOF group. Firstly, the analysis of univariate regression was conducted to evaluate the association between the predictive factors and OF. Second, factors with a P<0.1 in the univariate regression analysis were included in multivariate logistic regression analysis. A P value of 0.05 were deemed to be statistically significant. The logistic regression coefficients as well as the allocation of scoring points for each predictive factor, based on the regression coefficients, are given in Table 1. Our calculation of IRF model is based on the disease risk prediction model of the Framingham Heart Study [12]: (1) estimation of the regression coefficients of the multiple logistic regression; (2) organize the risk factors into meaningful categories and determine a reference value (midpoint) for each category, continuous variables were divided into groups, and the group median value was used as the reference value. Then, we need to determine the basic risk reference value of all risk factors. The basic risk reference value refers to that if a certain risk factor of the patient takes this value, the risk score is 0. The higher the value is, the higher the score is, and the higher the risk is. (3) Calculate the distance (D) between the classification of each risk factor and the basic risk. (4) We set the unit distance (B) for 1 point. (5) Each risk factor was scored, point=D/B.Table 1
Predictive scoring system for OF.
Preoperative factorβ coefficientPoints contributedIL 6 (pg/mL)<7007-1500.0152571411150-2503>2504IL 8 (pg/mL)<8.1008.1-21.30.2482667042>21.34CD19+ B lymphocyte s(%)<9009-140.1229277621>143


## 2.1. Data Collection

The peripheral blood samples of all patients were obtained within 24 h after admission to hospital. Peripheral blood lymphocyte subset assay was applied to detect the percentage of CD3+T lymphocytes, CD3+ CD4+T lymphocytes, CD3+ CD8+ cytotoxic T lymphocytes, CD16+ CD56+ natural killer cells, and CD19+ B lymphocytes. Blood routine examination was used to obtain white blood cell count and percentage. Blood samples were routinely collected on AP patients admitted to our hospital and then stored and subsequently analyzed for the concentration of 10 cytokines (IL1a, IL1b, IL4, IL6, IL8, IL10, IL13, MCP-1, IFNg, and TNFa) by using the Human Inflammation Array Q1 (QAH-INF-1-1) (RayBiotech, Norcross, America). Moreover, other IRF (CRP and PCT) data, several biochemical parameters, and basic clinical information were also collected.


## 2.2. Statistical Analysis

All statistical analyses were conducted by Stata/SE version 10.0 for Windows. Descriptive data were presented asmean±standarddeviation for quantitative variables or absolute and relative frequencies for qualitative variables. Student’s t-test or nonparametric tests (Chi-square and Mann–Whitney) were used to compare OF and NOF group. Firstly, the analysis of univariate regression was conducted to evaluate the association between the predictive factors and OF. Second, factors with a P<0.1 in the univariate regression analysis were included in multivariate logistic regression analysis. A P value of 0.05 were deemed to be statistically significant. The logistic regression coefficients as well as the allocation of scoring points for each predictive factor, based on the regression coefficients, are given in Table 1. Our calculation of IRF model is based on the disease risk prediction model of the Framingham Heart Study [12]: (1) estimation of the regression coefficients of the multiple logistic regression; (2) organize the risk factors into meaningful categories and determine a reference value (midpoint) for each category, continuous variables were divided into groups, and the group median value was used as the reference value. Then, we need to determine the basic risk reference value of all risk factors. The basic risk reference value refers to that if a certain risk factor of the patient takes this value, the risk score is 0. The higher the value is, the higher the score is, and the higher the risk is. (3) Calculate the distance (D) between the classification of each risk factor and the basic risk. (4) We set the unit distance (B) for 1 point. (5) Each risk factor was scored, point=D/B.Table 1
Predictive scoring system for OF.
Preoperative factorβ coefficientPoints contributedIL 6 (pg/mL)<7007-1500.0152571411150-2503>2504IL 8 (pg/mL)<8.1008.1-21.30.2482667042>21.34CD19+ B lymphocyte s(%)<9009-140.1229277621>143


## 3. Results

### 3.1. Characteristics of Enrolled Patients

A total of 100 patients were enrolled in this study, 94 patients were finally enrolled in the statistical analysis, and 6 patients were excluded due to the lack of IRF data (3 patients for PCT, 1 patient for CRP, 2 patients for PCT and CRP). Related information of clinical parameters was shown in Table2, including age, gender, etiology, number of OF, length of hospital stay, hospital cost, APACHE II score, RANSON score, and IRFs levels. As described in Table 2, biliary pancreatitis was the most common type of AP in this cohort (54/94); the incidence of persistent OF was nearly 16% (15/94). Another 5 patients had transient organ failure and were placed in the control group.Table 2
Characteristics of enrolled patients.
VariableNOF (n=79)OF (n=15)ZPAge (years) (mean±SD)50.7±16.347.1±18.20.774Sex, male/female33/467/80.821Time of onset (hours)0.7050-24 h28724-48 h27448-72 h244Etiology0.417Biliary,N (%)4311Alcohol,N (%)82Hypertriglyceridemia,N (%)131Other,N (%)151APPACHE II score3.8±2.55.7±3.20.011RANSON score1.1±0.92.3±9.8<0.001Length of hospital stay (days)9 (7, 12)16 (15, 36)-4.06<0.001Hospital cost (CNY∗)32435 (22203, 41689)75075 (62198, 174771)<0.001Inflammatory markers (pg/mL)Median (P25, P75)Median (P25, P75)IL1a2.35 (1.55, 4.25)1.55 (1.28, 2.56)-1.760.474IL1b3.32 (1.73, 8.18)3.16 (1.69, 11.38)-0.030.895IL42.33 (1.42, 4.75)2.08 (1.57, 3.67)-0.490.431IL624.43 (10.18, 42.3)106.9 (78.97, 188.13)-4.82<0.001IL84.52 (2.99, 7.86)14.14 (5.68, 18.03)-3.760.002IL101.51 (0.83, 2.62)8.25 (4.33, 14.98)-4.670.003IL130.58 (0.32, 0.76)0.31 (0.2, 0.65)-1.360.348MCP1152.64 (125.18, 192.46)234.11 (165.02, 280.94)-1.010.001IFNg1.56 (0.93, 2.49)1.21 (0.78, 2.24)-2.710.505TNFa5.78 (2.34, 11.16)4.4 (2.2, 11.55)-0.450.624PCT (ng/mL)0.42 (0.12, 1.18)2.07 (0.33, 23.38)-2.65<0.001CRP (mg/mL)90 (40.1, 90)90 (84.9, 90)-1.020.84Routine blood testWBC, ×109/L (mean±SD)10.7±4.113.3±4.90.244LY, ×109/L median (P25, P75)1.13 (0.85, 1.57)0.99 (0.85, 1.25)-1.1720.241MO, ×109/L median (P25, P75)0.53 (0.44, 0.78)0.45 (0.34, 0.95)-0.7490.454NE, ×109/L (mean±SD)8.7±3.911.4±4.40.161Immunity markers (%)CD3 T lymphocytes67.8±9.558.6±14.10.002CD3CD4 T lymphocytes39.8±10.031.3±8.30.007CD3CD8 cytotoxic T lymphocytes23.06 (17.56, 28.53)21.56 (13.21, 26.84)-1.10.272CD16+ CD56+ natural killer cells10.9 (6.7, 18.39)14.67 (7.63, 17.39)-0.490.624CD19+ B lymphocytes14.5 (10.2, 19.8)25.57 (17.9, 29.97)-3.929<0.001CD4CD8 cytotoxic T lymphocytes1.76 (1.29, 2.24)1.55 (1.09, 2.87)-0.0830.934∗CNY (Chinese Yuan), white blood cell counts (WBC), lymphocyte (LY), neutrophils (NE), and monocyte (MO).


### 3.2. Univariate and Multivariate Analysis Based on Clinical Information and IRFs

All patients enrolled in final analysis were divided into two groups according to the occurrence of OF (NOF vs. OF). Univariate analysis was selected to evaluate the association between OF occurrence and IRF levels (Table3). The analytic data demonstrated that the levels of cytokines IL6, IL8, IL10, and MCP1 were abnormally increased in OF group, and the variation tendency of CD19+ B lymphocytes, PCT, APACHE II score, and RANSON score were similar to these cytokines. However, the percentages of CD3+ CD4+ T lymphocytes were inversely decreased in OF group than that in NOF group. The other IRF levels were not markedly correlated with the presence of OF. These results indicated that IL6 (OR 1.01, 95% CI 1.004-1.024, P=0.003), IL8 (OR 1.20, 95% CI 1.08-1.33, P=0.001), IL10 (OR 1.26, 95% CI 1.08-1.47, P=0.003), MCP1 (OR 1.01, 95% CI 1.003-1.026, P=0.004), CD3+ CD4+ T lymphocytes (OR 0.92, 95% CI 0.87-0.98, P=0.01), CD19+ B lymphocytes (OR 1.11, 95% CI 1.04-1.18, P=0.001), PCT (OR 1.14, 95% CI 1.04-1.26, P=0.006), APACHE II score (OR 1.28, 95% CI 1.05-1.55, P=0.015), and RANSON score (OR 3.57, 95% CI 1.86-6.85, P<0.001) were potential predictors for OF in AP.Table 3
Univariate analyses of factors predicting OF.
Odds ratioStd. err.zP>∣z∣95% conf. intervalAge0.98688810.0171669-0.760.4480.9538086, 1.021115Gender0.81987580.4637188-0.350.7250.2705899, 2.484189APACHE II score1.2752590.12806482.420.0151.047415, 1.552667RANSON score3.5702971.1868673.83<0.0011.860978, 6.849637IL1a0.92045420.1064778-0.720.4740.7337286, 1.154699IL1b0.99755490.0183145-0.130.8940.9622973, 1.034104IL40.94865960.0635995-0.790.4320.8318493, 1.081873IL61.0143110.00492152.930.0031.004711, 1.024003IL81.2003820.06389833.430.0011.081455, 1.332386IL101.2636060.09836193.010.0031.084806, 1.471875IL130.50397050.3617765-0.950.340.123415, 2.057986MCP11.0118750.00415372.880.0041.003767, 1.020049IFNg0.87522060.1774049-0.660.5110.5882761, 1.302129TNFa0.98448690.0314253-0.490.6240.9247815, 1.048047PCT1.1445020.05633092.740.0061.039253, 1.260409CRP1.011830.01115441.070.2860.990202, 1.03393Cause0.62760130.1894851-1.540.1230.3472878, 1.13417CD3 T lymphocytes0.92985390.0238692-2.830.0050.8842285, 0.9778336CD3CD4 T lymphocytes0.92415350.0283805-2.570.010.8701696, 0.9814864CD3CD8 cytotoxic T lymphocytes0.95752430.0332156-1.250.2110.8945867, 1.02489CD16+ CD56+ natural killer cells1.0213150.0320970.670.5020.9603039, 1.086201CD19+ B lymphocytes1.1106880.03644593.20.0011.041504, 1.184468CD4CD8 cytotoxic T lymphocytes1.0500040.28064810.180.8550.6218407, 1.772975WBC, ×109/L1.0784870.06991621.170.2440.9498018, 1.224606LY, ×109/L1.5165590.99359760.640.5250.4199339, 5.476934MO, ×109/L0.39133180.253425-1.450.1470.1099791, 1.392451NE, ×109/L1.1012040.07619251.390.1640.961552, 1.261138Ca0.04738020.0504324-2.860.0040.0058825, 0.3816217To assess whether the potential risk factors (P<0.01) were independently associated with the progression of OF, multivariate logistic regression analysis was subsequently conducted. The related results suggested that IL6 (OR 1.10, 95% CI 1.001-1.03, P=0.038), IL8 (OR 1.28, 95% CI 1.01-1.63, P=0.043), and CD19+ B lymphocytes (OR 1.13, 95% CI 1.003-1.275, P=0.045) were independent predictors for OF development, as well as RANSON scores (OR 9.38, 95% CI 1.90-46.39, P=0.006) (Table 4).Table 4
Multivariate analyses of factors predicting OF.
Odds ratioStd. err.zP>∣z∣95% conf. intervalRANSON score9.3806737.6506682.740.0061.896762, 46.39328PCT1.0993120.06126631.70.0890.9855586, 1.226196IL61.0153740.00746282.080.0381.000852, 1.030107IL81.2818020.15748172.020.0431.007494, 1.630794IL101.0564820.14494620.40.6890.8073833, 1.382434MCP10.9821850.0102192-1.730.0840.9623586, 1.00242CD19+ B lymphocytes1.1308030.06937720.0451.002684, 1.275292ROC curve analysis were applied to evaluate the predictive value of IRF model for OF, as well as APACHE II score and RANSON score. As shown in Table5 and Figure 1, the AUC of IRFs model was 0.905 while the AUCs of APACHE II, RANSON, IL6, IL8, and CD19+ B lymphocytes, respectively, were 0.684, 0.821, 0.895, 0.808, and 0.821. In total, the AUC of IRFs model was higher compared to APACHE II and RANSON score, as well as individual factor enrolled in IRF model.Table 5
predictive value of IRFs model and other clinical factors for OF.
Preoperative factorAUCStd. err.95% conf. intervalLower limitUpper limitAPACHEII0.6840.0720.5410.826RANSON0.8210.0550.7120.929IL60.8950.0440.8080.981IL80.8080.0560.6990.916CD19+ B lymphocytes0.8210.0470.7290.913IFRs0.9050.0360.8350.976Figure 1
ROC curve analysis of the IRFs and relative clinical factors in prediction of OF.Therefore, these three inflammation-related variables (IL6, IL8, and CD19+ B lymphocytes) were finally selected to develop a prediction model (IRFs model). All details of IRF model are shown in Table1, including classification method, regression coefficient, and assigned score for each predictive factor.


### 3.3. Point Allocation for Predictors of Organ Failure Based on Regression Coefficients

Results of multivariate logistic regression analysis suggested that IL6, IL8, and CD19+ B lymphocytes are independent predictors for organ failure. Therefore, these three factors were selected into the final model for the construction of the prediction score (IRF score). The prediction scores based on this model range from 1 to 11, and AP patients with higher score had increasing risk to suffer with OF. The risk of AP patients with OF progressively increased as the score increased. ROC curve analysis was applied to evaluate the predictive value of new score system for OF, as well as APACHE II score and RANSON score. The AUC 0.86 (95% CI: 0.78-0.94) of ROC analysis demonstrated a satisfactory discrimination of new score system and higher than other two criteria (APACHE II and RANSON score) (Figure2). A score of 6 points was suggested to be the optimal cut-off value (Youdenindex=0.511) to divide the risk strata.Figure 2
Comparisons of the AUCs between scoring systems in prediction of OF.The mean ROC curve of tenfold crossvalidation of the new prediction score (IRF score) is shown in Figure3, which gave an AUC of 0.80±0.10 (mean±1std.dev.) indicating satisfactory separating capacity for our model.Figure 3
Mean receiver-operating characteristic (ROC) curve of the new prediction score for prediction of OF in tenfold crossvalidation.


## 3.1. Characteristics of Enrolled Patients

A total of 100 patients were enrolled in this study, 94 patients were finally enrolled in the statistical analysis, and 6 patients were excluded due to the lack of IRF data (3 patients for PCT, 1 patient for CRP, 2 patients for PCT and CRP). Related information of clinical parameters was shown in Table2, including age, gender, etiology, number of OF, length of hospital stay, hospital cost, APACHE II score, RANSON score, and IRFs levels. As described in Table 2, biliary pancreatitis was the most common type of AP in this cohort (54/94); the incidence of persistent OF was nearly 16% (15/94). Another 5 patients had transient organ failure and were placed in the control group.Table 2
Characteristics of enrolled patients.
VariableNOF (n=79)OF (n=15)ZPAge (years) (mean±SD)50.7±16.347.1±18.20.774Sex, male/female33/467/80.821Time of onset (hours)0.7050-24 h28724-48 h27448-72 h244Etiology0.417Biliary,N (%)4311Alcohol,N (%)82Hypertriglyceridemia,N (%)131Other,N (%)151APPACHE II score3.8±2.55.7±3.20.011RANSON score1.1±0.92.3±9.8<0.001Length of hospital stay (days)9 (7, 12)16 (15, 36)-4.06<0.001Hospital cost (CNY∗)32435 (22203, 41689)75075 (62198, 174771)<0.001Inflammatory markers (pg/mL)Median (P25, P75)Median (P25, P75)IL1a2.35 (1.55, 4.25)1.55 (1.28, 2.56)-1.760.474IL1b3.32 (1.73, 8.18)3.16 (1.69, 11.38)-0.030.895IL42.33 (1.42, 4.75)2.08 (1.57, 3.67)-0.490.431IL624.43 (10.18, 42.3)106.9 (78.97, 188.13)-4.82<0.001IL84.52 (2.99, 7.86)14.14 (5.68, 18.03)-3.760.002IL101.51 (0.83, 2.62)8.25 (4.33, 14.98)-4.670.003IL130.58 (0.32, 0.76)0.31 (0.2, 0.65)-1.360.348MCP1152.64 (125.18, 192.46)234.11 (165.02, 280.94)-1.010.001IFNg1.56 (0.93, 2.49)1.21 (0.78, 2.24)-2.710.505TNFa5.78 (2.34, 11.16)4.4 (2.2, 11.55)-0.450.624PCT (ng/mL)0.42 (0.12, 1.18)2.07 (0.33, 23.38)-2.65<0.001CRP (mg/mL)90 (40.1, 90)90 (84.9, 90)-1.020.84Routine blood testWBC, ×109/L (mean±SD)10.7±4.113.3±4.90.244LY, ×109/L median (P25, P75)1.13 (0.85, 1.57)0.99 (0.85, 1.25)-1.1720.241MO, ×109/L median (P25, P75)0.53 (0.44, 0.78)0.45 (0.34, 0.95)-0.7490.454NE, ×109/L (mean±SD)8.7±3.911.4±4.40.161Immunity markers (%)CD3 T lymphocytes67.8±9.558.6±14.10.002CD3CD4 T lymphocytes39.8±10.031.3±8.30.007CD3CD8 cytotoxic T lymphocytes23.06 (17.56, 28.53)21.56 (13.21, 26.84)-1.10.272CD16+ CD56+ natural killer cells10.9 (6.7, 18.39)14.67 (7.63, 17.39)-0.490.624CD19+ B lymphocytes14.5 (10.2, 19.8)25.57 (17.9, 29.97)-3.929<0.001CD4CD8 cytotoxic T lymphocytes1.76 (1.29, 2.24)1.55 (1.09, 2.87)-0.0830.934∗CNY (Chinese Yuan), white blood cell counts (WBC), lymphocyte (LY), neutrophils (NE), and monocyte (MO).


## 3.2. Univariate and Multivariate Analysis Based on Clinical Information and IRFs

All patients enrolled in final analysis were divided into two groups according to the occurrence of OF (NOF vs. OF). Univariate analysis was selected to evaluate the association between OF occurrence and IRF levels (Table3). The analytic data demonstrated that the levels of cytokines IL6, IL8, IL10, and MCP1 were abnormally increased in OF group, and the variation tendency of CD19+ B lymphocytes, PCT, APACHE II score, and RANSON score were similar to these cytokines. However, the percentages of CD3+ CD4+ T lymphocytes were inversely decreased in OF group than that in NOF group. The other IRF levels were not markedly correlated with the presence of OF. These results indicated that IL6 (OR 1.01, 95% CI 1.004-1.024, P=0.003), IL8 (OR 1.20, 95% CI 1.08-1.33, P=0.001), IL10 (OR 1.26, 95% CI 1.08-1.47, P=0.003), MCP1 (OR 1.01, 95% CI 1.003-1.026, P=0.004), CD3+ CD4+ T lymphocytes (OR 0.92, 95% CI 0.87-0.98, P=0.01), CD19+ B lymphocytes (OR 1.11, 95% CI 1.04-1.18, P=0.001), PCT (OR 1.14, 95% CI 1.04-1.26, P=0.006), APACHE II score (OR 1.28, 95% CI 1.05-1.55, P=0.015), and RANSON score (OR 3.57, 95% CI 1.86-6.85, P<0.001) were potential predictors for OF in AP.Table 3
Univariate analyses of factors predicting OF.
Odds ratioStd. err.zP>∣z∣95% conf. intervalAge0.98688810.0171669-0.760.4480.9538086, 1.021115Gender0.81987580.4637188-0.350.7250.2705899, 2.484189APACHE II score1.2752590.12806482.420.0151.047415, 1.552667RANSON score3.5702971.1868673.83<0.0011.860978, 6.849637IL1a0.92045420.1064778-0.720.4740.7337286, 1.154699IL1b0.99755490.0183145-0.130.8940.9622973, 1.034104IL40.94865960.0635995-0.790.4320.8318493, 1.081873IL61.0143110.00492152.930.0031.004711, 1.024003IL81.2003820.06389833.430.0011.081455, 1.332386IL101.2636060.09836193.010.0031.084806, 1.471875IL130.50397050.3617765-0.950.340.123415, 2.057986MCP11.0118750.00415372.880.0041.003767, 1.020049IFNg0.87522060.1774049-0.660.5110.5882761, 1.302129TNFa0.98448690.0314253-0.490.6240.9247815, 1.048047PCT1.1445020.05633092.740.0061.039253, 1.260409CRP1.011830.01115441.070.2860.990202, 1.03393Cause0.62760130.1894851-1.540.1230.3472878, 1.13417CD3 T lymphocytes0.92985390.0238692-2.830.0050.8842285, 0.9778336CD3CD4 T lymphocytes0.92415350.0283805-2.570.010.8701696, 0.9814864CD3CD8 cytotoxic T lymphocytes0.95752430.0332156-1.250.2110.8945867, 1.02489CD16+ CD56+ natural killer cells1.0213150.0320970.670.5020.9603039, 1.086201CD19+ B lymphocytes1.1106880.03644593.20.0011.041504, 1.184468CD4CD8 cytotoxic T lymphocytes1.0500040.28064810.180.8550.6218407, 1.772975WBC, ×109/L1.0784870.06991621.170.2440.9498018, 1.224606LY, ×109/L1.5165590.99359760.640.5250.4199339, 5.476934MO, ×109/L0.39133180.253425-1.450.1470.1099791, 1.392451NE, ×109/L1.1012040.07619251.390.1640.961552, 1.261138Ca0.04738020.0504324-2.860.0040.0058825, 0.3816217To assess whether the potential risk factors (P<0.01) were independently associated with the progression of OF, multivariate logistic regression analysis was subsequently conducted. The related results suggested that IL6 (OR 1.10, 95% CI 1.001-1.03, P=0.038), IL8 (OR 1.28, 95% CI 1.01-1.63, P=0.043), and CD19+ B lymphocytes (OR 1.13, 95% CI 1.003-1.275, P=0.045) were independent predictors for OF development, as well as RANSON scores (OR 9.38, 95% CI 1.90-46.39, P=0.006) (Table 4).Table 4
Multivariate analyses of factors predicting OF.
Odds ratioStd. err.zP>∣z∣95% conf. intervalRANSON score9.3806737.6506682.740.0061.896762, 46.39328PCT1.0993120.06126631.70.0890.9855586, 1.226196IL61.0153740.00746282.080.0381.000852, 1.030107IL81.2818020.15748172.020.0431.007494, 1.630794IL101.0564820.14494620.40.6890.8073833, 1.382434MCP10.9821850.0102192-1.730.0840.9623586, 1.00242CD19+ B lymphocytes1.1308030.06937720.0451.002684, 1.275292ROC curve analysis were applied to evaluate the predictive value of IRF model for OF, as well as APACHE II score and RANSON score. As shown in Table5 and Figure 1, the AUC of IRFs model was 0.905 while the AUCs of APACHE II, RANSON, IL6, IL8, and CD19+ B lymphocytes, respectively, were 0.684, 0.821, 0.895, 0.808, and 0.821. In total, the AUC of IRFs model was higher compared to APACHE II and RANSON score, as well as individual factor enrolled in IRF model.Table 5
predictive value of IRFs model and other clinical factors for OF.
Preoperative factorAUCStd. err.95% conf. intervalLower limitUpper limitAPACHEII0.6840.0720.5410.826RANSON0.8210.0550.7120.929IL60.8950.0440.8080.981IL80.8080.0560.6990.916CD19+ B lymphocytes0.8210.0470.7290.913IFRs0.9050.0360.8350.976Figure 1
ROC curve analysis of the IRFs and relative clinical factors in prediction of OF.Therefore, these three inflammation-related variables (IL6, IL8, and CD19+ B lymphocytes) were finally selected to develop a prediction model (IRFs model). All details of IRF model are shown in Table1, including classification method, regression coefficient, and assigned score for each predictive factor.


## 3.3. Point Allocation for Predictors of Organ Failure Based on Regression Coefficients

Results of multivariate logistic regression analysis suggested that IL6, IL8, and CD19+ B lymphocytes are independent predictors for organ failure. Therefore, these three factors were selected into the final model for the construction of the prediction score (IRF score). The prediction scores based on this model range from 1 to 11, and AP patients with higher score had increasing risk to suffer with OF. The risk of AP patients with OF progressively increased as the score increased. ROC curve analysis was applied to evaluate the predictive value of new score system for OF, as well as APACHE II score and RANSON score. The AUC 0.86 (95% CI: 0.78-0.94) of ROC analysis demonstrated a satisfactory discrimination of new score system and higher than other two criteria (APACHE II and RANSON score) (Figure2). A score of 6 points was suggested to be the optimal cut-off value (Youdenindex=0.511) to divide the risk strata.Figure 2
Comparisons of the AUCs between scoring systems in prediction of OF.The mean ROC curve of tenfold crossvalidation of the new prediction score (IRF score) is shown in Figure3, which gave an AUC of 0.80±0.10 (mean±1std.dev.) indicating satisfactory separating capacity for our model.Figure 3
Mean receiver-operating characteristic (ROC) curve of the new prediction score for prediction of OF in tenfold crossvalidation.


## 4. Discussion

As reported in previous studies, organ failure was a common but dangerous clinical manifestations in AP, and organ failure was also regarded as the main reason for the death of AP patients. Therefore, the identification of individuals with organ failure in the early stage was the key step for the treatment of these patients. Different kinds of IRF-induced SIRS were reported to be responsible for the occurrence and development of organ failure in AP, and it has been confirmed that some of these IRFs could function as independent indicators for such pathogenetic process. Several researches suggested that abnormal alteration of some immune cells could play the predictive role in the organ failure process, such as CD4+ lymphocyte, CD19+ B lymphocyte, neutrophil-lymphocyte ratio, and CD14hiCD16-monocytes [6, 13, 14]. Some studies focused on cytokines demonstrated that various cytokine levels were significantly associated with the incidence of organ failure, including IL6, IL8, IL17, MCP-1, and TNF-a [15]. Moreover, the inflammatory factors PCT and CRP were also reported to work as predictors for the organ failure of AP [10, 16]. However, nearly all studies evaluated the predictive efficiency of potential factors, respectively, and the combined value of these factors for organ failure prediction was still to be investigated.In our study, we firstly assess the predictive value of several IRFs for the development of organ failure in AP, including the percentages and/or cell counts of immune cells, the levels of 10 cytokines, and 2 common inflammation biomarkers. Similar to previous studies, the results based on this aspect suggested that the increased levels of IL6, IL8, and CD19+ B lymphocyte were, respectively, correlated with the increased incidence rate of organ failure [6]. IL6 and IL8 were frequently used cytokines which could closely participate in the differentiation of certain immune cells. For example, IL6 was the key accelerator in the final differentiation of B-cells into Ig-secreting cells, and IL8 had the function on neutrophil activation [17]. IL6 and IL8 could take part in a wide variety of inflammation-associated disease by inducing the inflammatory response, such as rheumatoid arthritis [18], ulcerative colitis [19], and acute pancreatitis [20]. In acute pancreatitis, the levels of these cytokines in peripheral blood were abnormally increased and significantly associated with the disease development. And some basic studies focused on AP pathogenesis also found that IL6 and IL8 could function as the important medium for different incentives induced and aggravated inflammatory process [21]. In addition, activated CD19+ B lymphocyte could also influence the process of various inflammatory diseases by promoting antigen-mediated humoral immunity [22]. By comprehensively considering the important roles of these three factors in inflammatory process and our analytic data, we selected them to further conduct prediction model.In clinical practice, two AP scoring systems, APACHE II and RANSON, were widely used to evaluate to the severity of such disease. In fact, these two models were relatively cumbersome due to the excessive parameters enrolled. Therefore, researchers specialized in AP have been trying to build simple and feasible models to replace them. In recent years, several simple models have been developed to evaluate the severity of acute pancreatitis in previous investigations. For example, Hong et al. published a prediction model for severe acute pancreatitis by using SIRS, albumin, BUN, and pleural effusion [23]; Qiu et al. conducted a model for the prediction of multiple organ failure in moderately severe and severe acute pancreatitis by using HCT, K-time, IL-6, and creatinine [24]; intra-abdominal pressure and body mass index were applied to establish prediction model to evaluate the AP severity [25]. However, there was no similar model based on IRFs (promote the inflammatory process in AP) ever published. In our study, we initiatively conducted an effective model for the prediction of AP-related organ failure by using three IRFs, including IL6, IL8, and CD19+ B lymphocyte. This model was more simple and convenient than the widely accepted AP points-scoring systems (APACHE II and RANSON), and the AUCs, sensitivity, and specificity of this model were similar to others, sometimes even better. This shows that the prediction efficiency of our model was partially better than existing points-scoring model. On the other hand, we can use this simple model to identify AP patients with high risk of organ failure rapidly and effectively by only detecting IL6, IL8, and CD19+ B lymphocyte.However, our study also has some limitations. Firstly, there were data missing in this study, especially the information of some other IRFs related to AP which were not enrolled, such as IL17, IL12, and M2 monocyte, which might be responsible for the selection bias. Secondly, all data applied in our study were collected retrospectively, which might reduce the reliability of results based on these data. Lastly, because of the funding problem, the sample size of our study was insufficient, not only the total number of all AP patients but also the number of AP patients with organ failure, and the validation cohort for IRF model was also absent. In the future research, we will expand the sample size and set up a validation cohort to further confirm the working efficiency of our prediction model, as well as enroll more potential IRFs.


## 5. Conclusions

In conclusions, IL6, IL8, and CD19+ B lymphocyte were reliable predictors for the organ failure of AP, and the prediction model has satisfactory working efficiency to identify AP patients with high risk of organ failure. This study may improve clinical prevention and treatment strategies of AP progression.



---

*Source: 4906768-2021-09-11.xml*