ORIGINAL RESEARCH article

Front. Endocrinol., 01 December 2025

Sec. Clinical Diabetes

Volume 16 - 2025 | https://doi.org/10.3389/fendo.2025.1501530

Risk stratification and predictive value of glucose variability for the development of post-acute pancreatitis diabetes mellitus

  • 1. Department of Critical Care Medicine, Nanjing Jinling Hospital, Affiliated Hospital of Medical School, Nanjing University, Nanjing, China

  • 2. Department of Critical Care Medicine, Nanjing Drum Tower Hospital, Affiliated Hospital of Medical School, Nanjing University, Nanjing, China

  • 3. Department of Critical Care Medicine, Affiliated Jinling Hospital, School of Medicine, Southeast University, Nanjing, China

  • 4. Department of Endocrinology, Nanjing Jinling Hospital, Affiliated Hospital of Medical School, Nanjing University, Nanjing, China

  • 5. Health Medicine Department, Jinling Hospital, Affiliated Hospital of Medical School, Nanjing University, Nanjing, China

  • 6. Department of Critical Care Medicine, Jinling Hospital, Nanjing Medical University, Nanjing, Jiangsu, China

Abstract

Background:

To assess whether glucose variability (GV) during acute pancreatitis (AP) predicts post-acute pancreatitis diabetes mellitus (PPDM-A), significantly affecting patient life quality.

Methods:

This study was performed during 2016–2020 at Jinling Hospital, with a 3-year follow-up for each patient. Cox proportional hazards model was used to evaluate the association of GV with the possibility of developing PPDM-A. Dose-response relationships of GV with the three-year probability of PPDM-A were characterized based on a restricted cubic splines (RCS) model. GV was analyzed to predict the ability for PPDM-A by calculating area under the receiver operating characteristic curves (AUCs).

Results:

PPDM-A rates rose from 16% at one year to 27.3% at three years post-AP. Multivariate Cox analysis indicated that the largest amplitude of glycemic excursions (LAGE) exhibited independent association with an increased PPDM-A risk within 3 years (HR = 1.21, 95% CI: 1.05–1.38, P <0.01). RCS results identified optimum LAGE threshold as 5.1, with significantly higher 3-year PPDM-A rates of abnormal LAGE group (LAGE ≥5.1 mmol/L) when compared with normal LAGE group (LAGE <5.1 mmol/L, P <0.001). AUCs for LAGE in predicting PPDM-A incidence in 12, 24, and 36 months were 0.883 (95% CI: 0.862–0.930), 0.916 (95% CI: 0.887-0.945), and 0.926 (95% CI: 0.895-0.948), respectively.

Conclusions:

LAGE in hospital stay accurately predicts PPDM-A. Further investigation plays an essential role in determining whether GV-targeting interventions can confer favorable clinical outcomes.

Introduction

Post-acute pancreatitis diabetes mellitus (PPDM-A) is a specifical diabetes subtype that develops following acute pancreatitis (AP) and is characterized by impaired insulin secretion and glucose metabolism (). It is associated with higher rates of morbidity and mortality, as well as a decreased quality of life (). The reported incidence of PPDM-A varies widely in different studies, ranging from 5% to 70% (). This variability can be caused by differences in study populations, follow-up durations, and diagnostic criteria. In the Chinese population, recent studies have provided specific epidemiological data: A multicenter retrospective cohort study involving 1,804 patients reported a 6.2% PPDM-A prevalence over a median follow-up of 3.04 years (); Additionally, Ma et al. constructed a nomogram for predicting PPDM-A in a Chinese cohort of 616 patients, and identified that 20.0% of participants developed diabetes within 3 months of AP onset (). While traditional risk factors for PPDM-A, such as etiology, obesity, and severity of AP, have been widely investigated, there is still conflicting evidence among a variety of studies (, , ). Therefore, it is essential to identify individuals at high risk for PPDM-A during the recovery phase from an AP episode, which is a vital component of the ‘holistic prevention of pancreatitis’ approach ().

Several studies showed that in‐hospital hyperglycemia might be a potential indicator for the development of PPDM-A and its associated complications (, , ). However, these findings have primarily depended on single time point assessments of admission blood glucose or fasting blood glucose levels. This approach, focusing on static blood glucose measurements, may not fully capture the association of dynamic glucose fluctuations with patient prognosis. In recent years, there has been an increasing interest in exploring the impact of glucose variability (GV) on the clinical prognosis of hospitalized individuals. Studies have exhibited significant correlations between increased GV during hospitalization and poorer short- and long-term outcomes. Mendez et al. investigated non-critically ill hospitalized patients (medical/surgical units), finding GV independently predicted prolonged hospitalization (), while Akirov et al. studied general surgery ward populations (non-cardiac/non-trauma), reporting high GV increased 30-day mortality risk by 78% (). A study specifically highlighted the predictive value of early-stage GV in plasma glucose levels for the development of diabetes (). Despite the obtained findings, the specific role of GV in the context of risk stratification and prediction of PPDM-A has not been thoroughly investigated. In addition, the lack of an optimum threshold for GV indices in accurately predicting PPDM-A and effectively stratifying risk has limited the clinical application.

Therefore, a study was performed on Chinese adults following AP by collecting their medical records and conducting telephone follow-ups. This cohort study was used to determine PPDM-A cumulative rates within a 3-year period and to investigate potential risk factors in association with its development. Therefore, this study specifically aimed to (1) propose the novel GV dynamic index and determine optimum threshold of GV indices for enhancing risk stratification, and (2) explore the association of GV indices with the development PPDM-A.

Methods

Study design and ethics statement

This cohort study was performed to explore all inpatients who experienced a first attack of AP between January 1, 2016, and December 31, 2020, at the Therapy Center for Severe Acute Pancreatitis (SAP), Jinling Hospital. Data were acquired based on Electronic Medical Record system (AP-Database), with an approval from the Institutional Review Board (2019NZKY-003-03). Analysis was conducted following the specific regulations. Moreover, no informed consents were needed in this retrospective study. To protect privacy, patient information was anonymized and deidentified. Our study was carried out strictly following Declaration of Helsinki and corresponding ethical guidelines.

Study participants

Subjects enrolled in this study were adults (>18 years) experiencing the initial AP attack. Exclusion criteria were presented as follows: I. Individuals with impaired fasting glucose (IFG, 5.6-6.9 mmol/L) or HbA1c levels ≥5.7% (indicating prediabetes) (). It should be noted that baseline BMI and HbA1c were not included as risk factors. Most participants had SAP with hypercatabolic states (), making acute-phase BMI unreliable due to fluid shifts and weight loss (). HbA1c accuracy may be compromised by anemia, acute kidney injury, or hypoalbuminemia (). II. Those with other pancreatic injuries, including pancreatectomy, trauma, neoplasms, hemochromatosis, cystic fibrosis, rare genetic diseases, and fibrocalculous pancreatopathy; III. Individuals undergoing severe cardiac/hepatic/renal insufficiency, or those diagnosed with malignancies; IV. Pregnant or postpartum individuals; and V. Individuals with over10% missing data. Defined as missing 3 or more items among 26 core clinical indicators; these indicators include demographic characteristics, AP-related factors, glucose-related metrics, biochemical/inflammatory markers, and complication/prognosis data, detailed in Data Extraction section.

Data extraction

Demographics (such as sex and age), duration of hospitalization, AP cause and severity, levels of pancreatic enzymes (amylase and lipase), blood glucose levels, serum calcium, interleukin-6 (IL-6), C-reactive protein (CRP), hepatic and renal function markers, lipid profiles, infection status, white blood cell count (WBC), hematocrit (HCT), and platelet count (PLT) were extracted from the AP-Database. Trained professionals were responsible for follow-up visits of subjects after hospital discharge. Follow-up visits were performed through telephone for those who were unable to visit our center. The above indicators constitute 26 core items, with glucose-related metrics (LAGE, SD, CV) merged as one to avoid overcounting. An individual is deemed to have ‘over 10% data missingness’ if missing 3+ of these 26 items and thus excluded.

Definitions and classifications

Based on Atlanta classification after revision, AP was diagnosed according to two among these three criteria (): (1) sudden severe and persistent upper abdominal pain usually radiating to the back. (2) serum lipase and/or amylase levels ≥ 3-fold upper limit of normal, and (3) representative imaging results consistent with AP on abdominal ultrasound, contrast-enhanced CT, and MRI. AP severity was classified into mild, moderate and severe following revised Atlanta classification (). Subjects of study were categorized into mild AP when they did not have systemic or local complications or organ failure. Moderate AP was defined in patients developing systemic or local complications without persistent organ failure (< 48 h). Severe AP was defined in patients with persistent single/multiple organ failure (> 48h). Consistent with 2021 American Diabetes Association (ADA) standards, PPDM-A was defined as new-onset diabetes following AP (>3 months) without underlying diabetes prior to AP attack ().

Definition of GV indices and abnormal GV

The GV indices, including Standard Deviation of Blood Glucose (SDBG), Largest Amplitude of Glycemic Excursions (LAGE), and Coefficient of Variation (CV) are employed to assess blood glucose level variability among people with diabetes. Firstly, serum glucose levels measured in the course of disease were detected in all patients. Subsequently, SDBG was defined as arithmetic SDBG values. CV was defined as the SDBG to average blood glucose level ratio. LAGE was calculated as the difference between the highest and lowest blood glucose values (mmol/L) recorded during a 24-hour observation period, serving as a key metric of glycemic variability independent of average glucose concentrations (). In healthy adults with normal glucose tolerance, a LAGE value of under 5.7 mmol/L is considered within the normal physiological range ().In our study, frequent capillary blood glucose (CBG) measurements (finger-prick) offer a clinically feasible method for LAGE calculation. CBG was measured every 2 hours and LAGE was calculated as the difference between the maximum and minimum glucose value across a 24-hour period (requiring 12 measurements). LAGE values were averaged across three monitoring days to enhance reliability. To identify the GV index most strongly associated with the occurrence of PPDM-A within three years, least absolute shrinkage and selection operator (LASSO) regression were used. This model allowed us to screen the most relevant GV index from the above-mentioned indices. To further analyze the optimum thresholds of these GV indices, RCS regression was used for evaluating dose-response associations of indices of GV with the 3-year probability of PPDM-A. Regression involved possible confounding factor adjustment at 5, 25, 75, and 95th percentiles of GV indices. The optimum threshold was GVs value that corresponded to lower limit of HR 95% CI > 1. Therefore, GV value ≥ threshold was deemed as having aberrant GV; otherwise, they were classified as having normal GV.

Statistical analyses

A Shapiro–Wilk test was used to evaluate continuous variables were normally distributed. Based on the findings, continuous variables were indicated by mean and standard deviation or median and interquartile range. Normally-distributed data were compared by t-test, whereas non-normally-distributed counterparts were compared by adopting Mann-Whitney test.

Hazard ratio (HR) having 95% confidence intervals (CI) was calculated by using univariate and multivariate Cox proportional hazards regression models. Variables satisfying P < 0.05 through univariate regression and clinically relevant characteristics were selected as candidate risk factors. LASSO regression satisfying P < 0.05 was used to screen best-fit covariates in multivariate regression (Supplementary Figure 1). Correlation heatmaps were plotted to test multi-collinearity among variables in the multivariate COX regression models (Supplementary Figure 2).

In addition, 3-year PPDM-A incidence was evaluated by drawing Kaplan–Meier survival curves, which was compared with Logrank test in abnormal and normal LAGE groups. Area under receiver operating characteristic (ROC) curves (AUC) was calculated to assess the PPDM prediction ability. Delong test was utilized to compare ROC curves.

Subgroup analyses stratified by different factors, including age (≤50, >50 years), gender, individuals with or without walled-off necrosis(WON), and those with recurring AP, were performed. Moreover, the association of abnormal GV with the above subgroup variables was analyzed, and the correlation of abnormal GV with 3-year PPDM-A incidence was analyzed in every subgroup.

P < 0.05 (two-tailed) represented statistical significance. R Programming Language (version 4.0.3) and STATA (version.15.0) were used for statistical analyses. R packages “ggrcs” [0.2.4 version] and “rms” [6.2–0 version] were utilized to visualize and analyze the RCS model; “survival” [3.2–11 version] was used in survival analysis; “timeROC” [0.4 version] was applied to estimate time-dependent ROC curve and AUC); whereas “forestplot” [2.0.1 version] was used to visualize forest plots during sensitivity analysis.

Results

Finally, totally 388 subjects satisfying the eligibility criteria were enrolled for analysis. DM was observed in 27.3% (106/388) of patients after the onset of AP(Figure 1). PPDM-A developed in 62, 91, and 106 patients at 1, 2, and 3 years after the onset of AP, with cumulative rates of 16%, 23.5%, and 27.3%, respectively.

Figure 1

Patient features

Table 1 displays patient demographic and clinical data, with 388 cases being included in this study. Based on diabetes, patients were classified into two groups. There were no significant differences in terms of gender and age between the two groups. However, individuals with PPDM-A showed longer hospital stays (12 (6–39) vs 7 (4–18), p<0.001), elevated levels of uric acid (UA), triglycerides (TG), procalcitonin (PCT), amylase, IL-6, CRP, admission blood glucose (ABG), LAGE, coefficient of variation (CV), and HCT (all P < 0.05), as well as a higher severity of AP (P = 0.006) in relative to non-diabetic patients. Individuals with PPDM-A had a higher incidence of hyperlipidemic AP (p<0.001) than those who had normal glucose levels, and they also exhibited a higher AP recurrence rate, as well as a higher occurrence of Wall-off necrosis (WON), acute kidney injury (AKI), acute respiratory distress syndrome (ARDS), and infected pancreatic necrosis (IPN) when compared with individuals without diabetes (all P < 0.05).

Table 1

CharacteristicAP without DM (n=282)PPDM-A (n=106)P value
Age [median (IQR)], years45 (36, 54)44(35, 52)0.62
Male, n(%)178 (63)72 (68)0.378
Hospital stay[median (IQR)], days7(4, 18)12(6, 39)<0.001
Severity, n(%)0.006
Mild52 (18)6 (6)
Moderately severe71 (25)27 (25)
Severe159 (56)73 (69)
Etiology, n(%)<0.001
Biliary161 (57)36 (34)
Hyperlipidemic94 (33)62 (58)
Alcoholic6 (2)4 (4)
Idiopathic21 (7)4 (4)
ABG [median (IQR)], mmol/L6.9 (5.7, 8.5)10.5 (7.9, 13.9)<0.001
LAGE [median (IQR)], mmol/L3.9 (2.3, 5.8)10.4 (8.2, 13.9)<0.001
SDBG1.94 (1.13, 2.74)5.18 (4.11, 6.25)<0.001
CV%28 (17, 36)47 (38, 58)<0.001
IL-6 [median (IQR)], ng/L55 (25, 151)136 (83, 270)<0.001
PCT [median (IQR)], ug/L0.4 (0.2, 1.3)0.5 (0.2, 1.9)0.013
Amylase [median (IQR)], U/L87 (41, 228)57 (35, 140)0.046
Lipase [median (IQR)], U/L302(64, 674)243(72,631)0.351
LDH [median (IQR)], U/L671 (460, 1,162)712 (426, 1,297)0.336
TG [median (IQR)], mmol/L1.72 (1.14, 2.92)2.58 (1.56, 5.41)<0.001
TC [median (IQR)], mmol/L3.26 (2.44, 4.16)3.38 (2.55, 4.89)0.145
Albumin [median (IQR)], g/L30.2 (26.5, 34.6)31.2 (28.0, 34.8)0.175
Cr [median (IQR)], umol/l53 (41, 71)55 (42, 77)0.179
BUN [median (IQR)], mmol/l4.3 (3.0, 6.5)4.9 (3.2, 7.7)0.189
UA [median (IQR)], umol/L204 (136, 288)241 (162, 327)0.010
WBC, ×109/L10.5 (7.6, 14.1)11.4 (8.1, 15.0)0.171
PLT, ×106/L193 (143, 255)214 (151, 299)0.032
HCT
CRP [median (IQR)], mg/L114 (52, 185)210 (135, 263)<0.001
WON, n(%)154 (55)78 (74)<0.001
IPN, n(%)85 (30)46 (43)0.014
ARDS, n(%)65 (23)45 (42)<0.001
AKI, n(%)54 (19)32 (30)0.02
Recurrence, n(%)30 (11)22 (21)0.009

Patient demographics and clinical characteristics.

Values are presented as means ± standard deviation (SD), median (IQR) or frequencies (percentages). *Statistically significant (P < 0.05) for patients with PPDM-A vs. normoglycemia.

CRP, C-reactive protein; LDH, lactate dehydrogenase; Cr, creatinine; BUN, blood urea nitrogen; UA, uric acid; TG, triglyceride; TC, total cholesterol; ABG, admission blood glucose; LAGE, Largest amplitude of glycemic excursions; CV, coefficient of variation; PCT, procalcitonin; WON, Wall-off necrosis; IPN, infected pancreatic necrosis; ARDS, acute respiratory distress syndrome; AKI, acute kidney injury.

The meaning of the bold values is statistical significance (P < 0.05) for patients with PPDM-A vs AP without DM.

Relationship of GV indices and clinical features to PPDM-A incidence

Univariate Cox regression analysis indicated that elevated levels of ABG, LAGE, SDBG, and CV, as well as a higher recurrence rate of AP, occurrence of WON, IPN, ARDS, and AKI, longer hospital stays, elevated IL-6, TG, Cr, BUA, HCT, PLT, and CRP, as well as the severity and etiology of AP, were all associated with an increased risk of 3-year probability of PPDM-A (Table 2).

Table 2

VariablesHR[95% CI]P1aHR[95% CI]P2
ABG1.161.13-1.19p<0.001--
LAGE1.251.21-1.29p<0.0011.271.22-1.32p<0.001
SD1.471.40-1.55p<0.001--
CV%1.061.05-1.07p<0.001--
Severityreference----
3.051.26-7.380.014--
3.451.50-7.920.004--
Etiologyreference----
2.541.68-3.83p<0.001--
2.720.97-7.65p=0.058--
0.840.30-2.37p=0.748--
WON2.081.35-3.21p=0.001--
ARDS2.091.42-3.07p<0.001--
AKI1.681.11-2.54p=0.014--
Age1.000.98-1.01p=0.516--
LOS, days1.011.01-1.02<0.001--
recurrence1.901.19-3.040.007--
Gender, male1.190.79-1.780.411--
IL-6, per 10 ng/L1.011.01-1.02p<0.0011.011.00-1.010.004
PCT, ug/L1.000.99-1.010.553--
Amylase, per 10 ng/L0.990.99-1.000.253--
LDH, per 10 U/L1.001.00-1.010.0531.011.00-1.01<0.001
TG, mmo/L1.081.05-1.11<0.001--
TC, mmo/L1.000.98-1.030.926--
Albumin, g/L1.020.99-1.050.242--
Cr, per 101.021.01-1.040.002--
BUN, mmol/L1.041.02-1.06<0.001--
UA, per 10 ng/L1.021.01-1.04<0.001--
WBC, ×109/L1.020.99-1.060.129--
HCT0.060.00-0.920.044--
PLT, per 10*109/L1.021.01-1.040.0031.021.01-1.040.003
CRP, per 10 mg/L1.081.05-1.10<0.0011.041.02-1.07<0.001

Univariate and multivariate association of LAGE and 3-year incidence of PPDM-A in COX regression model.

HR Hazard ratio, aHR adjustive Hazard ratio, P1 p-value for univariate analysis, P2 p-value for multivariate analysis, CRP, C-reactive protein; LDH, lactate dehydrogenase; Cr, creatinine; BUN, blood urea nitrogen; UA, uric acid; TG, triglyceride; TC, total cholesterol; ABG, admission blood glucose; LAGE, Largest amplitude of glycemic excursions; CV, coefficient of variation; WON, Wall-off necrosis; IPN, infected pancreatic necrosis; ARDS, acute respiratory distress syndrome; AKI,acute kidney injury. - denoted variables that were screened out by the least absolute shrinkage and selection operator (LASSO) regression.

Consistent with the -2log-likelihood and binomial family type measure, LASSO regression was conducted out with R software through 10-fold K cross-validation to centralize and normalize those involved variables and to select the optimum lambda value. “Lambda.lse” provides the model with least independent variables yet high performance. Finally, those 5 potential risk factors from the candidate variables were selected by LASSO regression, including LAGE, IL-6, PLT, CRP, and LDH. In addition, the most significant features selected by LASSO were adopted for multivariate Cox proportional hazards analyses.

The multivariate Cox regression analysis identified the higher levels LAGE (aHR = 1.21, 95% CI: 1.05–1.38, P <0.001), CRP (aHR = 1.04, 95% CI: 1.02–1.07, P = 0.001), IL-6 (aHR = 1.01, 95% CI: 1.00–1.01, P = 0.004), LDH (aHR=1.01, 95% CI: 1.00–1.01, P <0.001) and PLT (aHR = 1.02, 95% CI: 1.01–1.04, P = 0.003) as independent prognostic factors for PPDM-A.

Multivariate RCS analyses adjusted for IL-6, CRP, LDH, and platelet count suggested that LAGE was correlated with 3-year probability of PPDM-A in AP cases (Figure 2). The dose–response curve morphology exhibited monotonic increase. The lower limits of 95% CI associated with HR for 3-year probability of PPDM-A were > 1 at LAGE ≥ 5.1 mmol/L. Therefore, the optimum threshold was LAGE = 5.1 mmol/L. LAGE < 5.1 mmol/L and LAGE ≥ 5.1mmol/L were deemed as normal and abnormal LAGE, respectively.

Figure 2

In the normal LAGE group, increasing LAGE levels were not significantly correlated with the 3-year likelihood of developing PPDM-A. However, for abnormal LAGE group, the increasing LAGE levels were related to a higher PPDM-A risk over 3 years.

Clinical outcomes for both LAGE groups

Participants were categorized into normal (205, 52.8%) or abnormal LAGE group (183, 47.2%). The baseline clinical characteristics showed significant variances between the two groups. Compared to the normal LAGE group, patients with abnormal LAGE (≥5.1 mmol/L) had significantly longer hospital stays, higher rates of severe acute pancreatitis, more frequent hyperlipidemic etiology, and increased complications such as ARDS, AKI, and walled-off necrosis (Table 3). Moreover, patients with abnormal LAGE exhibited a higher incidence of PPDM-A within 3 years compared with those with normal LAGE. (Supplementary Figure 3, Log-rank P < 0.001).

Table 3

CharacteristicNormal LAGE (205)Abnormal LAGE (183)P-value
Age, (mean ± SD), years46 ± 1445 ± 120.535
LOS, (mean ± SD), years12 ± 1627 ± 33<0.001
IL-6, (mean ± SD), ng/L118 ± 171200 ± 3260.002
PCT, (mean ± SD), ug/L3.2 ± 18.63.7 ± 11.20.757
Amylase, (mean ± SD), U/L193 ± 260177 ± 2630.567
LDH, (mean ± SD), U/L787 ± 5801,259 ± 1,284<0.001
TG, (mean ± SD), mmo/l2.41 ± 2.584.09 ± 5.59<0.001
TC, (mean ± SD), mmo/l3.54 ± 1.494.48 ± 9.730.200
Albumin, (mean ± SD), g/l31.3 ± 6.230.9 ± 6.20.518
Cr, (mean ± SD), umol/l66 ± 6991 ± 1010.006
BUN, (mean ± SD), mmol/l5.2 ± 3.86.8 ± 7.50.010
UA, (mean ± SD), umol/l225 ± 108249 ± 1470.077
WBC (×109), mean ± SD11.5 ± 5.411.9 ± 5.80.524
HCT, mean ± SD0.34 ± 0.070.33 ± 0.070.153
Plt(×109), mean ± SD218 ± 111226 ± 1170.523
CRP, (mean ± SD), mg/L121 ± 83166 ± 92<0.001
Severity<0.001
Mild49 (24)9 (5)
Moderately severe52 (25)46 (25)
Severe104 (51)128 (70)
Etiology<0.001
Biliary125 (61)72 (39)
Hyperlipidemic62 (30)94 (51)
Alcoholic6 (3)4 (2)
Idiopathic12 (6)13 (7)
WON, n(%)106 (52)126 (69)<0.001
IPN, n(%)59 (29)72 (39)0.028
ARDS, n(%)34 (17)76 (42)<0.001
AKI, n(%)29 (14)57 (31)<0.001
Recurrence, n(%)20 (10)32 (17)0.026
Male, n(%)129 (63)121 (66)0.512
PPDM-A, n(%)4 (2)102 (56)<0.001

Baseline clinical characteristics of patients stratified by the optimal cutoff point of LAGE.

Data are presented as n (%) or mean ± SD. CRP, C-reactive protein; LDH, lactate dehydrogenase; Cr,creatinine; BUN, blood urea nitrogen; UA, uric acid; TG, triglyceride; TC, total cholesterol; ABG, admission blood glucose; LAGE, Largest amplitude of glycemic excursions; CV, coefficient of variation; WON, Wall-off necrosis; IPN, infected pancreatic necrosis; ARDS, acute respiratory distress syndrome; AKI,acute kidney injury; LOS length of hospital stay.

Predicting ability of LAGE for PPDM-A incidence

Differences in predicting ability of LAGE, and the etiology of AP and IL-6 for PPDM-A at different time points were compared by using time-dependent ROC curve analysis. The AUCs of the LAGE in determining 1-, 2-, 3-year PPDM-A risk were 0.883 (95% CI: 0.862–0.930), 0.916 (95% CI: 0.887-0.945), 0.926 (95% CI: 0.895-0.948), separately (Figure 3a). The CRP discrimination was 0.66 (95% CI: 0.575–0.735), 0.70 (95% CI: 0.633–0.765), and 0.69 (95% CI: 0.636–0.752) in determining 1-year, 2-year, 3-year probability of PPDM-A respectively. (Figure 3b). In previous studies, the etiology of AP has been regarded as the risk factor related to DM occurrence after AP, and AUC of the etiology of AP for 3-year incidence of PPDM-A was 0.60 (95% CI: 0.544-0.657) (Figure 3c).

Figure 3

Sensitivity analysis on the correlation of abnormal LAGE with 3-year PPDM-A incidence

Participants of abnormal LAGE group exhibited the significantly increased PPDM-A risk over 3 years compared with those in normal LAGE group, even when some confounders were adjusted (aHR =41.75, 95% CI:15.36–113.49, P < 0.001). Moreover, abnormal LAGE was consistently associated with an elevated possibility of developing PPDM-A within 3 years across various subgroups, including patients aged >50 and ≤50 years, men and women, participants with/without WON, and those with recurring AP. It is vital to emphasize that even in patients with non-severe AP, the abnormal LAGE group still exhibited a relatively higher incidence of PPDM-A over 3 years (HR = 41.66, 95% CI:9.95–174.47) (Figure 4).

Figure 4

Discussion

Individuals with PPDM-A confront a significant risk to their health and lives due to frequent hypoglycemia episodes and complications (, ). Our study was aimed at providing a more comprehensive understanding of PPDM-A, building on our 2016 research that had limitations in scope and definition (, ). It was found that the prevalence of PPDM-A increased from 16% at one year to 27.3% at three years after AP onset, which is higher than the 23% cumulative incidence reported in meta-analyses (, ). This discrepancy may be caused by methodological differences and regional variations (, ).

In this study, we innovatively demonstrated the predictive value of GV for PPDM-A, identifying LAGE = 5.1 mmol/L as optimum threshold in risk stratification. Our findings were strong in different subgroups, even for non-severe AP patients. Elevated glucagon and pancreatic polypeptide levels, alongside decreased ghrelin and peptide YY responses, are characteristics of PPDM-A and could serve as potential biomarkers (). Our research builds on the known association of GV with mortality and diabetes, demonstrating it as the significant risk factor related to PPDM-A. Different from previous studies concentrating on early glucose fluctuations, our analysis considered the entire disease course, from admission to discharge, providing a more comprehensive assessment of GV’s predictive value for PPDM-A development (, , , ).

Abnormal blood glucose levels in hospitalization are related to the higher PPDM-A risk (). Univariate analysis suggested that the glucose level at admission served as the risk factor related to PPDM-A development within 3 years. However, the glucose level at admission lacked significant correlation with PPDM-A in the multivariate model. Lasso regression was employed to eliminate the interference of multiple covariates between admission glucose and glycemic variability, enhancing the credibility of our results. Nevertheless, single-point measurements could not sufficiently indicate the real association of blood glucose with PPDM-A, while LAGE was the dynamic factor overcoming static data limitations. Sensitivity analysis suggested that our results were robust.

This study aimed to focus on clarifying the candidate mechanisms between blood glucose fluctuations and new-onset DM development in AP patients. A plausible explanation for our findings refers to that exposure to fluctuating glucose levels heightens oxidative stress, causing pronounced endothelial dysfunction (, ). This, in turn, exacerbates the systemic inflammatory response and insulin resistance, therefore elevating an individual’s risk of developing hyperglycemia in the future (, ). Clinically, elevated LAGE (>5.1 mmol/L) predicts hyperglycemic excursions (>4.4 mmol/L) and hypoglycemic risk (>3.48 mmol/L), while inversely associating with time-in-range in type 2 diabetes ().

The association between PPDM-A and elevated serum LDH levels can be attributed to several factors. AP results in extensive tissue damage and an inflammatory response, causing the release of LDH from damaged cells into the blood. Hypoxia and metabolic disturbances during pancreatic inflammation further exacerbate tissue damage, which can increase LDH levels. Elevated LDH levels also reflect the severity of pancreatic necrosis and inflammation, which are related to pancreatic beta-cell dysfunction and increased insulin resistance, leading to diabetes ().

In addition to LAGE and LDH, increased IL-6 and CRP were also significantly associated with the PPDM-A development. This aligns with research findings indicating that increased IL-6 levels after AP are associated with a higher risk of hyperglycemia (). IL-6 and CRP increase during inflammation, impacting glucose metabolism through insulin resistance (, ). This can lead to impaired insulin receptor phosphorylation, inhibiting glycogen synthesis and promoting glucose release (). Further investigations are necessary to understand the link of abnormal glucose metabolism with adipocytokines for better preventive and therapeutic strategies for PPDM-A.

The higher basic platelet level is the candidate indicator for a higher thrombosis risk, enhanced inflammatory response, while decreased vascular endothelial function (). The increase in platelet count has been shown as a potential predictor for PPDM-A in this study. This is likely resulted from the inflammatory response triggered by AP which can cause insulin resistance and impaired glucose metabolism. The elevated platelet count may reflect the ongoing inflammation and oxidative stress in the body, which are known risk factors associated with diabetes occurrence (). It is important to monitor patients with elevated platelet counts following AP closely for developing diabetes, as early intervention and management can contribute to preventing or delaying the onset of the disease.

Furthermore, the relationship between AP severity and the risk of PPDM-A requires careful interpretation. Evidence regarding AP severity as a risk factor is mixed, with some studies supporting the link and others finding no significant association(6,7,10,11). In our study, AP severity was significantly associated with PPDM-A in univariate analysis (Table 2), and patients in the high LAGE group had a higher incidence of severe AP (Table 3). However, AP severity was not retained as an independent predictor in the multivariate model after LASSO regression. We propose a plausible explanation wherein the impact of AP severity on PPDM-A is largely mediated through its effect on LAGE. More intense inflammatory insults and metabolic stress during AP likely predispose patients to greater glycemic excursions. The resulting glucose variability, quantified by LAGE, then acts as a more direct and measurable driver of future diabetes risk. Therefore, these findings suggest that while the initial injury and severity set the stage, the subsequent metabolic and inflammatory dysregulation, particularly glucose variability, may be more proximal determinants of PPDM-A risk.

However, this study had several limitations that need to be acknowledged. First, as a retrospective analysis conducted at a single institution using electronic medical records, this study may be subject to recall and selection biases commonly associated with historically documented patient information. Therefore, the generalizability of the conclusions should be confirmed in future multi-center, prospective studies with larger cohorts. Second, the inability to obtain blood specimens during the follow-up period precluded the assessment of islet function (e.g., by insulin or C-peptide release tests), although such functional evaluations would have substantially strengthened the investigation of PPDM-A’s pathological mechanisms. Third, while in-hospital glucose monitoring enabled the assessment of GV during the acute phase of AP, the absence of glucose monitoring data in the chronic phase limited a comprehensive evaluation of its association with PPDM-A and restricted insights into the temporal dynamics of GV’s impact. In addition, conventional metabolic indices such as BMI and HOMA-IR were not included in our analysis due to concerns regarding their reliability in the acute phase of severe pancreatitis. This decision was based on two main considerations: (1) the impracticality of obtaining accurate anthropometric measurements in critically ill patients, and (2) the substantial confounding effect of acute pancreatitis–related β-cell dysfunction on insulin resistance indices ().

Conclusion

To conclude, this study indicated that secondary diabetes is prevalent among AP patients. GV exerts a vital role in risk prediction for PPDM-A. Abnormal GV is the efficient way to stratify and manage risk, which can be used to predict the possibility of PPDM-A clinically. More prospective research is needed to clarify the mechanisms of GV promoting the incidence of PPDM-A and to determine whether GV interventions positively influences the improvement of clinical outcome.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

Data were acquired based on Electronic Medical Record system (AP-Database), with an approval from the Institutional Review Board (2019NZKY-003-03). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

SQ: Conceptualization, Formal Analysis, Investigation, Writing – original draft. YL: Data curation, Methodology, Writing – review & editing. WW: Investigation, Writing – review & editing. WN: Investigation, Writing – review & editing. SZ: Investigation, Writing – review & editing. YX: Supervision, Writing – review & editing. PG: Funding acquisition, Investigation, Supervision, Writing – review & editing. WY: Supervision, Writing – review & editing. ZT: Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by funding from the National Natural Science Foundation of China (No. 82270678and 82102275).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer LL declared a shared affiliation, with no collaboration, with one of the authors, YL, to the handling editor at the time of the review.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2025.1501530/full#supplementary-material

References

Summary

Keywords

acute pancreatitis, post-acute pancreatitis diabetes mellitus, glucose variability, risk stratification, largest amplitude of glycemic excursions

Citation

Qian S, Liu Y, Wang W, Ni W, Zhao S, Xu Y, Gu P, Yu W and Tong Z (2025) Risk stratification and predictive value of glucose variability for the development of post-acute pancreatitis diabetes mellitus. Front. Endocrinol. 16:1501530. doi: 10.3389/fendo.2025.1501530

Received

25 September 2024

Revised

07 November 2025

Accepted

13 November 2025

Published

01 December 2025

Volume

16 - 2025

Edited by

Prof. Dr. med. Peter C Ambe, Witten/Herdecke University, Germany

Reviewed by

Brittany Bruggeman, University of Florida, United States

Ling Li, Southeast University, China

Ahmed E. Sherif, University of Edinburgh, United Kingdom

Updates

Copyright

*Correspondence: Zhihui Tong, ; Wenkui Yu, ; Ping Gu,

† These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics