Participants properties
A total of 577 patients in the PNDRFAP cohort were enrolled from 667 subjects (Fig. 1A). There were no significant differences in gender, BMI, drinking history, smoking history, hypertension, diabetes, MMSE score, duration of anesthesia, duration of surgery, estimated volume of infusion, and intraoperative blood loss between the postoperative delirium and non-postoperative delirium groups (Table 1). For the PNDABLE cohort, we screened 667 participants, of whom 621 finally met the study criteria and were included (Fig. 1B). The two groups, postoperative delirium and non-postoperative delirium, were compared using basic clinical data such as gender, BMI, smoking history, MMSE score, duration of anesthesia, estimated volume of infusion, intraoperative blood loss and duration of surgery. The disparity was not statistically significant (Table 2).
Table 1
Characteristics of participants in PNDRFAP
| PNDRFAP |
Characteristic | NPOD (n = 474) | POD (n = 103) | P |
Age [year, M(Q)] | 62(11) | 75(8) | < 0.001 |
Male [n(%)] | 277 (58.40) | 70 (68.00) | 0.074 |
BMI [kg.m− 2, M(Q)] | 25.00(5.00) | 26.00 (5.00) | 0.182 |
Education [year, M(Q)] | 12(3) | 6(4) | < 0.001 |
Smoking history [n(%)] | 131(27.60) | 22 (21.40) | 0.191 |
Drinking history [n(%)] | 161(34.00) | 26 (25.20) | 0.086 |
Hypertension [n(%)] | 228(48.10) | 57(55.30) | 0.183 |
Diabetes [n(%)] | 95(20.00) | 26(25.20) | 0.240 |
MDAS [scores, M(Q)] | 1(6) | 15(8) | < 0.001 |
MMSE [scores, M(Q)] | 28(3) | 28(3) | 0.826 |
Total bilirubin [µmol/L, M(Q)] | 8.5(3.30) | 13.70(3.60) | < 0.001 |
Direct bilirubin [µmol/L, M(Q)] | 5.60(3.53) | 3.40(2.10) | < 0.001 |
Indirect bilirubin [µmol/L, M(Q)] | 2.50(1.33) | 10.20(3.80) | < 0.001 |
Duration of anesthesia (min) | 145.20(79.80) | 139.80(70.20) | 0.615 |
Duration of surgery (min) | 94.80(64.80) | 90.00(57.60) | 0.896 |
Estimated volume of infusion (ml) | 1100.00(500.00) | 1100.00(600.00) | 0.058 |
Intraoperative blood loss (ml) | 20.00(195.00) | 50.00(190.00) | 0.648 |
Continuous variable use Student's t test or Mann-Whitney U, Categorical variable use Chi-square test. |
The numerical variables of normal distribution are statistically described by Average±Standard deviation [‾x±s]. |
Non-normally distributed numerical variables are statistically described by Median (Interquartile spacing) [M(Q)]. |
Categorical variables are statistically described by Sample size (Percent) [n(%)]. |
*P < 0.05 |
Table 2
Characteristics of participants in PNDABLE
| PNDABLE |
Characteristic | NPOD (n = 504) | POD (n = 117) | P |
Age [year, M(Q)] | 61.50(11) | 74(9) | < 0.001 |
Male [n(%)] | 291 (57.74) | 76 (64.96) | 0.152 |
BMI [kg.m− 2, M(Q)] | 25.32 (4.77) | 25.71(4.95) | 0.111 |
Education [year, M(Q)] | 12(3) | 7(4) | < 0.001 |
Smoking history [n(%)] | 138(27.38) | 26 (22.22) | 0.254 |
Drinking history [n(%)] | 177 (35.12) | 26 (22.22) | 0.007 |
Hypertension [n(%)] | 169 (33.53) | 54(46.15) | 0.010 |
Diabetes [n(%)] | 75 (14.88) | 31(26.50) | 0.003 |
MDAS [scores, M(Q)] | 1(6) | 15(7) | < 0.001 |
MMSE [scores, M(Q)] | 28(3) | 29(3.5) | 0.826 |
Total bilirubin [µmol/L, M(Q)] | 11.90(5.25) | 12.40(5.10) | < 0.001 |
Direct bilirubin [µmol/L, M(Q)] | 2.40(1.40) | 2.50(1.40) | < 0.001 |
Indirect bilirubin [µmol/L, M(Q)] | 9.40(4.30) | 9.70(4.50) | < 0.001 |
Aβ42 [pg/ml, M(Q)] | 357.46(238.64) | 288.73(153.00) | < 0.001 |
P-tau [pg/ml, M(Q)] | 39.98(15.79) | 45.18 (20.84) | < 0.001 |
T-tau [pg/ml, M(Q)] | 179.48(82.17) | 221.55(165.71) | < 0.001 |
Aβ42/P-tau [M(Q)] | 8.85(5.60) | 6.34(4.31) | < 0.001 |
Aβ42/T-tau [M(Q)] | 1.97(1.26) | 1.31(1.01) | < 0.001 |
Duration of anesthesia (min) | 150.00(79.80) | 145.20(65.70) | 0.522 |
Duration of surgery (min) | 100.20(70.20) | 90.00(62.70) | 0.973 |
Estimated volume of infusion (ml) | 1100.00(600.00) | 1100.00(500.00) | 0.288 |
Intraoperative blood loss (ml) | 20.00(195.00) | 50.00(190.00) | 0.438 |
PNDABLE, the Perioperative Neurocognitive Disorder and Biomarker Lifestyle study. |
POD, postoperative delirium; NPOD, non-postoperative delirium; |
Aβ42, amyloid beta42; t-Tau, total Tau; p-Tau, phosphorylated Tau; |
Continuous variable use Student's t test or Mann-Whitney U, Categorical variable use Chi-square test. |
The numerical variables of normal distribution are statistically described by Average±Standard deviation [‾x±s]. |
Non-normally distributed numerical variables are statistically described by Median (Interquartile spacing) [M(Q)]. |
Categorical variables are statistically described by Sample size (Percent) [n(%)]. |
*P < 0.05 |
In the PNDRFAP study, it was found that there were significant differences in the levels of total bilirubin, direct bilirubin, and indirect bilirubin between POD patients and NPOD patients (Total bilirubin, OR (95% CI): 1.839(1.645–2.056) P < 0.001; Direct bilirubin, OR (95% CI): 1.725(1.448–2.054) P < 0.001; Indirect bilirubin, OR (95% CI): 1.558(1.424–1.705) P < 0.001) (Table 3). In the PNDABLE cohort, the levels of total bilirubin, direct bilirubin, indirect bilirubin, Aβ42, tau, and P-tau were also found to be higher in the postoperative delirium group compared to those in the non-postoperative delirium group. Other variations in the data are illustrated in Fig. 2. We compared bilirubin markers and biomarkers in cerebrospinal fluid between patients with postoperative delirium (POD) and those without postoperative delirium (NPOD). The PNDABLE cohort study found that preoperative levels of total bilirubin, indirect bilirubin, direct bilirubin, P-tau, and T-tau were higher in patients with postoperative delirium (POD) compared to those without postoperative delirium (NPOD). [Total bilirubin, OR (95% CI): 1.203 (1.137–1.274), P < 0.001; Direct bilirubin, OR (95% CI): 1.724 (1.465–2.030), P < 0.001; Indirect bilirubin, OR (95% CI): 1.182 (1.109–1.260), P < 0.001; P-tau, OR (95% CI): 1.042 (1.025–1.060), P < 0.001; T-tau, OR (95% CI): 1.008 (1.006–1.011), P < 0.001] (Table 4).
Table 3
Logistic regression on analysis and sensitivity analysis in PNDRFAP study
| Model 1a | Model 2b | Model 3c | Model 4d |
OR (95% CI) | P | OR (95% CI) | P | OR (95% CI) | P | OR (95% CI) | P |
Total bilirubin, µmol/L | 1.839(1.645–2.056) | < 0.001 | 1.673(1.434–1.952) | < 0.001 | 1.680(1.430–1.973) | < 0.001 | 1.686(1.424–1.995) | < 0.001 |
Direct bilirubin, µmol/L | 1.725(1.448–2.054) | < 0.001 | 1.773(1.261–2.494) | < 0.001 | 1.824(1.281–2.595) | < 0.001 | 1.856(1.264–2.725) | 0.002 |
Indirect bilirubin, µmol/L | 1.558(1.424–1.705) | < 0.001 | 1.494(1.312–1.702) | < 0.001 | 1.504(1.310–1.727) | < 0.001 | 1.520(1.314–1.758) | < 0.001 |
Note: OR = odds ratio; 95%CI = 95% confidence interval; Aβ42 = β-amyloid42; T-tau = total tau; P-tau = phosphorylated tau |
aModel 1: Unadjusted |
bModel 2: Adjusted for age(50–90), sex, years of education and MMSE scores |
cModel 3: Adjusted for age(50–90), sex, years of education, MMSE scores, smoking history, drinking history, hypertension and diabetes |
dModel 4: Adjusted for age ≥ 65, sex, years of education, MMSE scores, smoking history, drinking history, hypertension and diabetes |
*P < 0.05 |
Table 4
Logistic regression on analysis and sensitivity analysis in PNDABLE study
| Model 1a | Model 2b | Model 3c | Model 4d |
OR (95% CI) | P | OR (95% CI) | P | OR (95% CI) | P | OR (95% CI) | P |
Total bilirubin, µmol/L | 1.203(1.137–1.274) | < 0.001 | 1.208(1.111–1.312) | < 0.001 | 1.210(1.110–1.318) | < 0.001 | 1.168(1.064–1.284) | < 0.001 |
Direct bilirubin, µmol/L | 1.724(1.465–2.030) | < 0.001 | 1.602(1.284–1.999) | < 0.001 | 1.772(1.395–2.251) | < 0.001 | 1.641(1.267–2.125) | < 0.001 |
Indirect bilirubin, µmol/L | 1.182(1.109–1.260) | < 0.001 | 1.2180(1.075–1.294) | < 0.001 | 1.170(1.064–1.286) | 0.001 | 1.127(1.014–1.253) | 0.026 |
Aβ42, pg/ml | 0.996(0.995–0.998) | < 0.001 | 0.995(0.993–0.997) | < 0.001 | 0.995(0.993–0.997) | < 0.001 | 0.994(0.992–0.997) | < 0.001 |
P-tau, pg/ml | 1.042(1.025–1.060) | < 0.001 | 1.030(1.008–1.053) | 0.008 | 1.031(1.007–1.056) | 0.010 | 1.027(1.0001–1.1055) | 0.046 |
T-tau, pg/ml | 1.008(1.006–1.011) | < 0.001 | 1.006(1.003–1.010) | < 0.001 | 1.007(1.003–1.010) | < 0.001 | 1.005(1.001–1.009) | 0.011 |
Aβ42/P-tau | 0.834(0.783–0.889) | < 0.001 | 0.827(0.762–0.898) | < 0.001 | 0.824(0.756–0.899) | < 0.001 | 0.809(0.731–0.896) | < 0.001 |
Aβ42/T-tau | 0.390(0.293–0.519) | < 0.001 | 0.434(0.303–0.621) | < 0.001 | 0.418(0.288–0.606) | < 0.001 | 0.442(0.293–0.666) | < 0.001 |
Note: OR = odds ratio; 95%CI = 95% confidence interval; Aβ42 = β-amyloid42; T-tau = total tau; P-tau = phosphorylated tau |
aModel 1: Unadjusted |
bModel 2: Adjusted for age(50–90), sex, years of education and MMSE scores |
cModel 3: Adjusted for age(50–90), sex, years of education, MMSE scores, smoking history, drinking history, hypertension and diabetes |
dModel 4: Adjusted for age ≥ 65, sex, years of education, MMSE scores, smoking history, drinking history, hypertension and diabetes |
*P < 0.05 |
Sensitivity analysis
The PNDRFAP (Table 3) and PNDABLE (Table 4) studies were analyzed in a sensitive manner with three models, namely Model 2, Model 3, and Model 4, in order to verify their reliability.
Adjustment for age, gender, educational level, and MMSE (model 2), there was significant correlation between total bilirubin, indirect bilirubin, direct bilirubin, P-tau, T-tau and the occurrence of POD based on a binary logistic regression model [Total bilirubin, OR (95% CI) 1.208 (1.111–1.312), P < 0.001; Direct bilirubin, OR (95% CI) 1.602 (1.284–1.999), P < 0.001; Indirect bilirubin, OR (95% CI) 1.218 (1.075–1.294), P < 0.001]. The results remained consistent in Model 3 after adjusting for hypertension, diabetes, smoking history, and alcohol consumption history. [Total bilirubin, OR (95% CI) 1.210 (1.110–1.318), P < 0.001; Direct bilirubin, OR (95% CI) 1.772 (1.395–2.251), P < 0.001; Indirect bilirubin, OR (95% CI) 1.170 (1.064–1.286), P < 0.001]. Based on Model 3, the age of the patients was controlled, and the results were consistent. [Total bilirubin, OR (95%CI)1.168(1.064–1.284) P < 0.001; Direct bilirubin, OR (95%CI)1.641(1.267–2.125) P < 0.001; Indirect bilirubin, OR (95%CI)1.127(1.014–1.253) P < 0.026] (Table 3). In PNDRFAP studies, the levels of total bilirubin, direct bilirubin, and indirect bilirubin were compared between the POD and NPOD groups. The results showed that the levels of total bilirubin, direct bilirubin, and indirect bilirubin were significantly higher in the POD group compared to the NPOD group [Total bilirubin, OR (95% CI) 1.839 (1.645–2.056), P < 0.001; Direct bilirubin, OR (95% CI) 1.725 (1.448–2.054), P < 0.001; Indirect bilirubin, OR (95% CI) 1.558 (1.424–1.705), P < 0.001]. These findings were consistent across Model 2, Model 3, and Model 4.
Calibration curve and nomogram
A calibration curve of the incidence of postoperative delirium (POD) in the PNDABLE study demonstrated strong agreement between predicted and observed values. Elevated total bilirubin, direct bilirubin, and indirect bilirubin may be the risk factors of postoperative delirium (POD). (Fig. 3A) There is an increased risk for postoperative delirium (POD) within certain ranges along with increases in total bilirubin, direct bilirubin and indirect bilirubin. Aβ42 may be the protective factor of POD after operation. As Aβ42 content decreases, there is an increase in the likelihood of POD development. Furthermore, regression coefficients were used to assign a value level to each factor in the model. The total score was then obtained by adding up each score. Based on the data from this study, the patient's score ranged from 159.58 to 301.59. In the end, the functional transformation relationship between the total score and the probability of postoperative delirium could be calculated. (Fig. 3B)
Receiver operating characteristic curve(ROC)
The ROC curve indicated that total bilirubin, direct bilirubin, and indirect bilirubin all exhibited strong predictive properties for postoperative delirium (total bilirubin AUC = 0.681; direct bilirubin AUC = 0.706; indirect bilirubin AUC = 0.633). Additionally, the combination of total bilirubin, direct bilirubin, and indirect bilirubin was a more accurate predictor of postoperative delirium than a single indicator (AUC = 0.713). If bilirubin markers are combined with biomarkers in cerebrospinal fluid, the prediction of postoperative delirium is more accurate (AUC = 0.812) (Fig. 4).
Decision curve analysis
The net benefit (y-axis) of biological indicators is plotted against various threshold probabilities (x-axis) on the decision curve (Fig. 5). It can be observed in Fig. 5 that total bilirubin, direct bilirubin, indirect bilirubin, and the combination of the three bilirubin indicators show differences in net benefit. The combination of multiple indicators provides a greater net benefit compared to total bilirubin, direct bilirubin, and indirect bilirubin, as well as in comparison to the clinical default strategies of "treat all" or "treat none". Therefore, compared to a single index, the combination of multiple indicators can better discriminate and guide clinical diagnosis and treatment.
Causal mediation analyses
Binary logistic regression in PNDABLE studies revealed that total bilirubin, direct bilirubin, indirect bilirubin, T-tau, P-tau, Aβ42/P-tau, and Aβ42/tau were positively correlated with POD, while Aβ42 was negatively correlated with POD. Therefore, total bilirubin, direct bilirubin, and indirect bilirubin, as well as T-tau, P-tau, Aβ42/P-tau, and Aβ42/tau, are all risk factors for POD, while Aβ42 is a protective factor for POD. We are aware that Aβ42, T-tau, P-tau, Aβ42/P-tau, and Aβ42/tau have a specific impact on POD. Therefore, we are further investigating whether total bilirubin, direct bilirubin, and indirect bilirubin can mediate the effects of T-tau, P-tau, Aβ42/P-tau, and Aβ42/tau on POD. Mediation analysis revealed that the association between total bilirubin and POD was mediated by T-tau, P-tau, Aβ42/P-tau, and Aβ42/tau. Similarly, the connection between direct bilirubin, indirect bilirubin, and POD was also mediated by T-tau, P-tau, Aβ42/P-tau, and Aβ42/tau (Fig. 6). The study revealed that the association between total bilirubin, direct bilirubin, indirect bilirubin, and postoperative delirium (POD) was influenced by CSF t-tau, with respective proportions of 9.08%, 10.26%, and 8.70%.
Prediction model development
In order to enhance clinical diagnosis and treatment, a dynamic nomogram is created. We developed a dynamic nomogram model to forecast the likelihood of postoperative delirium (POD)occurrence using eight independent predictors, which include Aβ42, T-tau, P-tau, Aβ42/tau, total bilirubin, direct bilirubin, and indirect bilirubin. Users can access the web page (https://bylqdsslyy.shinyapps.io/dynnomapp/), input variables on the page, and then click "Predict" in the lower left corner of the page to obtain the probability of POD occurring. In addition, the dynamic nomogram also provides a graphical summary, numerical summary, and model summary. (Fig. 7)