Data-Driven Predictive Modeling for Massive Intraoperative Blood Loss during Liver Transplantation: Integrating Machine Learning Techniques
1Division of Gastroenterological, General and Transplant Surgery, Department of Surgery, Jichi Medical University, Shimotsuke, Tochigi, Japan
2Department of Medical Informatics, Hirosaki University Hospital, Hirosaki, Aomori, Japan
*Correspondence author: Taiichi Wakiya, M.D., Ph.D. Division of Gastroenterological, General and Transplant Surgery, Department of Surgery, Jichi Medical University 3311-1 Yakushiji, Shimotsuke city, Tochigi, 329-0498, Japan Telephone: +81-285-58-7069, Fax: +81-285-58-7069 E-mail: wakiya.taiichi@jichi.ac.jpABSTRACT
Background
Massive intraoperative bleeding (IBL) in liver transplantation (LT) poses serious risks and strains healthcare resources necessitating better predictive models for risk stratification. As traditional models often fail to capture the complex, non-linear patterns underlying bleeding risk, this study aimed to develop data-driven machine learning models for predicting massive IBL during LT using preoperative factors.
Methods
Two hundred ninety consecutive LT cases from a prospective database were analyzed. Logistic regression models were built using 73 preoperative demographic and laboratory variables to predict massive IBL (≥ 80 mL/kg). The dataset was randomly split (70% training, 30% testing). The model was trained and validated through three-fold cross-validation, with backward stepwise feature selection iterated 100 times across unique random splits. The final model, based on a high stability index, was evaluated using the area under the curve (AUC).
Results
Massive IBL was observed in 141 patients (48.6%). In standard logistic regression, significant differences were found in 42 of 73 factors between groups stratified by massive IBL, however, substantial multicollinearity limited interpretability. In the feature selection across 100 iterations, the data-driven model achieved an average AUC of 0.840 in the validation and 0.738 in the test datasets. The final model, based on 11 selected features with a high stability index, achieved an AUC of 0.844. An easy-to-use online risk calculator for massive IBL was developed and is available at: https://tai1wakiya.shinyapps.io/ldlt_bleeding_ml/.
Conclusions
Our findings highlight the potential of machine learning in capturing complex risk factor interactions for predicting massive IBL in LT.
Article notes
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
The author(s) received no specific funding for this work.
Introduction
Liver transplantation (LT) is a high-risk procedure frequently accompanied by substantial intraoperative bleeding (IBL), which can lead to adverse outcomes.(1–3) Additionally, perioperative bleeding and transfusions consume significant human and financial resources.(4, 5) Accordingly, effective strategies to anticipate and reduce bleeding are essential for optimizing outcomes and maintaining the sustainability of transplant programs.
Despite ongoing advances, predicting critical IBL during LT remains a major clinical challenge. Risk stratification tools using conventional statistical methods, such as logistic regression, have been developed to identify patients at increased risk.(6–9) However, these models often assume linearity and independence among variables, limiting their ability to capture the complex, multifactorial nature of LT-associated bleeding. As a result, their accuracy and clinical utility remain suboptimal.
Machine learning offers a promising alternative by allowing the integration of numerous interrelated factors without requiring pre-specified assumptions about variable relationships. ML models can capture nonlinear associations and detect complex patterns that traditional approaches or clinical intuition may overlook. This capability also extends to handling high-dimensional datasets and identifying subtle but clinically relevant trends, potentially offering novel insights into perioperative risk.(10, 11) This capability may enhance the predictive performance of bleeding models and facilitate more personalized perioperative management.
Although ML has shown value in various surgical contexts,(12–14) its application in predicting bleeding outcomes in LT is still limited. Therefore, we aimed to develop a data-driven, ML-based model using prospectively collected preoperative variables to predict massive IBL in LT. Our objective was to establish a clinically applicable framework that outperforms conventional models and supports individualized surgical planning.
Materials and methods
Patients
We conducted a retrospective monocentric observational study using a prospectively maintained database. This study was approved by the Ethics Committee of Jichi Medical University (Approval No. 20-008). This study was designed and conducted in accordance with the principles of the Declaration of Helsinki and Istanbul. The need for written informed consent for the present study was waived by the Institutional Review Board of Jichi Medical University in view of its retrospective design, in accordance with national and local guidelines, considering the fact that all clinical/laboratory measurements and procedures were part of routine care. This study included 290 consecutive patients who underwent living donor LT at our facility between 2008 and 2024.
Surgical procedures and operative management
Donor hepatectomy was selected based on the recipient’s standard liver volume, weight, and graft volume determined using preoperative computed tomographic volumetry. For the recipient’s operation, inverted T-shaped or transverse incisions were made, and a total hepatectomy was performed. The graft hepatic vein was anastomosed to the stump of the recipient’s hepatic veins, which formed a single orifice, in an end-to-end manner. Hepatic artery reconstruction was routinely performed using microsurgical techniques. Choledocho-choledochostomy was the first choice for biliary reconstruction, except in cases in which the bile duct could not be used, such as with biliary atresia (BA) or primary sclerosing cholangitis.
Definition of massive intraoperative bleeding
IBL was calculated on the basis of the in/out balance of the operative field. At our institution, fluid loss from the abdominal cavity, including that from bile and lymphatics, is considered IBL. IBL was divided by body weight and referred to as adjusted intraoperative bleeding (aIBL, mL/kg). Massive aIBL was defined as aIBL ≥ 80 mL/kg. In clinical practice, the estimated circulating blood volume is approximately 70–80 mL/kg.(15) Therefore, aIBL exceeding 80 mL/kg is effectively equivalent to the loss of the entire circulating blood volume, indicating a clinically significant hemorrhagic event. Accordingly, 80 mL/kg serves as a physiologically meaningful and practical cutoff.
Data collection
For each patient, we collected demographic and laboratory data immediately before LT. The data were accessed for research purposes on April 22, 2025. Investigators had access to identifiable participant information during data collection; however, all data were de-identified prior to analysis. Seventy-three perioperative variables were extracted from the prospectively collected database. Patient demographic data included age, sex, body height, body weight, and etiology. Etiology was categorized as follows: acute liver failure (ALF), BA, graft failure, or other. Additionally, information regarding the history of previous abdominal surgery, ABO incompatibility, rituximab desensitization, and graft type was collected. Graft type was categorized as left lateral segment (LLS), left liver (LL), right liver (RL), right posterior segment (RPS), reduced LLS, or monosegment graft. To prevent statistical bias and enhance the generalizability of the model, incidence factors with a prevalence of less than 5%, such as spontaneous bacterial/fungal peritonitis and portal vein thrombosis, were excluded from the predictive variables. The 55 laboratory variables and their abbreviations are detailed in S1 Table. Prior to binary logistic regression analysis, the variables were standardized to have a mean of 0 and a standard deviation of 1, thereby ensuring that the logistic coefficients reflected the magnitude of influence on prediction. Given the study’s focus on identifying modifiable risk factors, particular emphasis was placed on preoperative blood test parameters, as these variables can be influenced by conditions such as infection and thrombosis and can be actively managed and optimized before surgery.
Statistical analyses
Continuous variables were expressed as medians (ranges) and analyzed using nonparametric methods for non-normally distributed data (Mann–Whitney U test). Categorical variables were reported as numbers (percentages) and analyzed using the chi-square test or Fisher’s exact test, as appropriate. Variables with a significant relationship to massive IBL in univariate analysis were used in a binary logistic regression model. The correlation between the two parameters was analyzed using the Spearman rank-order method. Differences were considered significant at P < 0.05. Statistical analyses were performed using GraphPad Prism (v10.2.3; GraphPad Software, San Diego, CA; USA, https://www.graphpad.com).
Data preparation and splitting
The overall workflow is illustrated using a block diagram in Fig 1. The dataset comprised 73 features along with corresponding binary labels. Data preprocessing involved a stratified random split, allocating 70% of the data for training and 30% for testing. To ensure a robust evaluation across different random splits, this process was repeated 100 times using unique random seeds (random_state = rr, where rr = 0, 1, …, 99).
Backward feature selection and model optimization
To identify the most informative predictors of massive intraoperative bleeding, we applied backward stepwise feature selection within a logistic regression framework. Model optimization and evaluation were conducted using stratified 3-fold cross-validation, repeated across 100 random data splits to ensure generalizability. At each iteration, variables were progressively removed based on their contribution to model performance, assessed using the area under the receiver operating characteristic curve (AUC). The feature set yielding the highest validation AUC was selected as the final model. Details of the feature elimination procedure, coefficient ranking, and model tuning are provided in the Supplementary Methods. All analyses were conducted using Python libraries including Scikit-learn and Pandas.(16)
Development of an online calculator for estimating the risk of massive aIBL
We developed an online calculator based on the final logistic regression model to estimate the probability of massive aIBL. To ensure model robustness, feature stability was evaluated across multiple random data splits using occurrence rate and absolute mean weight. Features consistently selected with strong predictive weights were incorporated into the final model. Prior to application, selected features were standardized and weighted using their regression coefficients to calculate the predicted probability of massive aIBL. The calculator is available at: https://tai1wakiya.shinyapps.io/ldlt_bleeding_ml/. Details of the feature selection metrics, weighting procedures, and standardization process are provided in the Supplementary Methods.
Results
Clinical significance of massive aIBL
Among the 290 patients, those with aIBL greater than 80 mL/kg exhibited significantly higher rates of graft loss compared to those with lower aIBL (Fig 2). These findings support the clinical validity of using 80 mL/kg as a reasonable and actionable threshold for defining massive aIBL in LT.
Comparison of patient characteristics in the massive aIBL and non-massive aIBL groups
The demographic data of the 290 enrolled patients are presented in Table 1, which also includes information on the IBL. Among all cases, massive aIBL was observed in 141 patients (48.6%). Except for indication and graft type, no significant differences were observed in the demographic data between the groups.
Comparison of laboratory data in the massive aIBL and non-massive aIBL groups
Laboratory data are presented in Table 2. Significant differences (P < 0.05) were observed in 40 of the 55 laboratory test items between the two groups. To predict the occurrence of massive aIBL, we conducted a binary logistic regression analysis with the occurrence of massive aIBL as the dependent variable. Forty significant predictor variables associated with massive aIBL, identified using univariate analysis (P < 0.05), were included in the binary logistic regression analysis. The analysis identified total bile acids (P = 0.013, odds ratio [OR] = 1.010, 95% confidence interval [CI]: 1.000–1.010), albumin (P = 0.014, OR = 0.228, 95% CI: 0.071–0.736), and ammonia (P = 0.043, OR = 1.020, 95% CI: 1.000–1.050) as significant predictors of incidences of massive aIBL. However, as noted in S2 Table, there was significant multicollinearity among many factors in this analysis. The presence of multicollinearity can lead to instability in estimates, wider confidence intervals, and increased errors in coefficient estimates, complicating interpretation and potentially leading to incorrect conclusions. These results indicate that to construct a predictive model that exhibits multicollinearity with our dataset, it is crucial to select variables objectively, excluding subjectivity, and to utilize alternative analytical methods.
Correlations between preoperative laboratory tests and intraoperative bleeding volume
Next, we examined the relationship between the preoperative laboratory test results and aIBL. Of the 55 parameters tested, 42 showed significant correlations with the bleeding volume. Positive correlations were found for 43.6% of the parameters, whereas 56.4% exhibited negative correlations. The top five positively correlated biomarkers were D-dimer (r = 0.588), prothrombin time international normalized ratio (r = 0.495), total bilirubin (r = 0.461), direct bilirubin (r = 0.438), and total bile acids (r = 0.363). Conversely, antithrombin III (r = -0.560), alpha 2-plasmin inhibitor (r = -0.555), protein C (r = -0.527), prothrombin time percent (r = -0.495), and plasminogen (r = -0.458) were highly negatively correlated (Fig 3). These findings suggest that these markers can potentially guide preoperative assessment and intervention.
Feature selection and model performance in logistic regression for massive aIBL prediction
We further explored a more robust prediction model for massive aIBL. We used the backward stepwise method for feature selection to construct our binary logistic regression model. During this process, we calculated the absolute magnitudes of the coefficients and removed the least influential covariates. This analysis incorporated laboratory variables as well as demographic data, including age, categorical age (under 6 years, 6-18 years, or 18 years and older), sex, body height, body weight, etiologies (ALF, BA, graft failure, or other), history of previous abdominal surgery, ABO incompatibility, rituximab desensitization, and graft type (LLS, LL, RL, RPS, reduced LLS, or monosegment). We started with all 73 features and systematically removed the least significant ones from the set. One of the representative feature selection processes is presented in S3 Table and S1 Fig. The process for narrowing down the variables is detailed in S3 Table. Antithrombin III was identified as the final covariate in the dataset. S1 Fig illustrates the performance of the model at each stage. The highest validation AUC of 0.867 was observed at the stage with 12 covariates. The AUC obtained in the test dataset was 0.728. To ensure a robust evaluation, this process was repeated 100 times across different random splits (S2 Fig). The average number of selected features across these 100 trials was 15.1. The average AUC in the validation and test datasets was 0.840 (standard deviation [SD] = 0.022) and 0.738 (SD = 0.046), respectively. These findings support a systematic approach for preoperatively predicting massive IBL in LT, suggesting the robustness and reliability of feature selection and model performance.
Identifying key predictors and constructing an online calculator for massive aIBL
To evaluate the consistency of selected features across 100 random data splits, index stability was assessed using occurrence rate and absolute mean weight (Fig 4). Based on the results from 100 iterations, features meeting the predefined criteria of an occurrence rate of at least 50% and an absolute mean weight of 0.5 or greater were selected. A total of 11 features were identified as robust predictors. Among them, six features exhibited positive mean weights, indicating a direct association with massive aIBL occurrence: D-dimer, ferritin, age, total bile acids, gamma-glutamyl transferase, and ammonia. Conversely, five features demonstrated negative mean weights, suggesting a potential protective effect: fibrinogen, platelet count, albumin, LLS graft, and antithrombin III. To enhance clinical applicability and ensure streamlined implementation, we developed a final prediction model for massive aIBL using these 11 selected features. In the test across all 290 cases, the model achieved an AUC of 0.844, with a sensitivity of 74.5%, specificity of 78.5%, positive predictive value of 76.6%, and negative predictive value of 76.5%. To further facilitate clinical adoption, we developed a user-friendly online calculator based on this risk model and made available online at: https://tai1wakiya.shinyapps.io/ldlt_bleeding_ml/.
Discussion
In this study, we developed and validated ML models to predict massive IBL during LT. Using preoperative demographic and laboratory data, we employed logistic regression models with backward feature selection to achieve high predictive accuracy. Furthermore, we developed an online calculator that provides physicians with an accessible tool to estimate an individualized risk of massive IBL in the preoperative setting. To our knowledge, this is the first study to apply a ML model for IBL prediction in LT and to translate these findings into a clinically applicable prediction tool. The true innovation of this study lies not only in the excellent accuracy of the model but also in showcasing a groundbreaking artificial intelligence-based approach with significant future potential for management in LT.
Previous studies have attempted to predict bleeding and transfusion requirements in LT by using various clinical parameters and traditional standard statistical models.(6–9) While these studies have provided valuable insights, their models often exhibit suboptimal predictive performance and limited clinical applicability. Standard methods in clinical studies have several limitations, including the selection of variables, confounding factors, and multicollinearity, as demonstrated in this study. Even in prospective studies, the issue of multicollinearity remains a challenge to resolve. Furthermore, standard methods often fail to capture the complex, non-linear relationships between variables that significantly influence bleeding risk. In contrast, our study goes beyond traditional analyses by employing a data-driven approach, allowing us to systematically identify not only well-established risk factors but also previously unrecognized predictors.
Our study leveraged ML techniques, which are used to identify intricate patterns and interactions among multiple factors, demonstrating superior accuracy and robustness in predicting massive IBL. ML offers several key benefits, including increased accuracy through the simultaneous analysis of multiple variables and their interrelationships, and the ability to identify complex patterns and correlations that traditional methods might overlook. By incorporating a comprehensive set of preoperative variables and employing advanced feature selection methods, our ML models offer a more reliable and practical tool for preoperative risk stratification in LT. These models effectively manage complex, high-dimensional datasets, uncover hidden correlations, and explore non-linear medical relationships, enhancing predictive capabilities and potentially leading to new scientific insights.
Crucially, we went beyond model development by translating these findings into a tangible clinical tool. To ensure practical utility in real-world clinical settings, we developed an online calculator based on our predictive model, providing an accessible platform for physicians to rapidly assess massive IBL risk. This helps bridge the gap between advanced computational modeling and clinical decision-making, while integrating ML into clinical practice, which may contribute to improving patient outcomes and optimizing resource management.
Antithrombin III and ferritin emerged as significant predictors in our machine learning-based model, providing new insights into bleeding risk in LT. In cardiovascular surgery, low preoperative antithrombin III levels have been consistently associated with increased intraoperative blood loss.(17) This has led to the consideration of antithrombin III supplementation as a potential intervention. However, the outcomes of such supplementation remain inconclusive, with several studies and meta-analyses reporting inconsistent results.(18–21) Similarly, although the available studies are both dated and limited in number, evidence supporting a positive effect of antithrombin III supplementation in LT remains limited.(22, 23) Regarding ferritin, no direct evidence has previously established its role in predicting bleeding in liver surgery, including LT. Ferritin reflects both iron metabolism and systemic inflammation, and its elevation may indicate a preoperative inflammatory state that predisposes patients to coagulopathy.(24) Inflammatory responses are known to disrupt hemostasis through mechanisms such as endothelial dysfunction, hyperfibrinolysis, and altered coagulation cascades.(25) Supporting this, our analysis identified a strong positive correlation between IBL and C-reactive protein, another established inflammatory marker; however, ferritin demonstrated greater stability as a predictive index than C-reactive protein. Together, these findings illustrate how machine learning can reveal biologically and clinically relevant predictors that are often overlooked by conventional approaches. Antithrombin III and ferritin may serve as valuable biomarkers for risk stratification and perioperative management in LT.
The findings of this current study should be interpreted in light of several limitations. First, it was a single-institution cohort study with a relatively small patient population. Additional training data could potentially enhance the prediction accuracy. Furthermore, the lack of external validation using an independent dataset is also a limitation. One reason for the absence of external validation is the argument made by some biostatistics experts in predictive research stating that independent verification can be misleading and should be omitted as a model evaluation step.(26, 27) These experts report that simulations confirming at least 100 events and 100 non-events are required for a reliable assessment of predictive performance. They suggest using all available data for model development, with some form of cross-validation or bootstrap validation to assess the statistical optimism in average predictive performance.(27, 28) Based on these biostatistical perspectives, we chose to build our ML model using all data with cross-validation, performing 100 iterations with randomly split datasets to ensure robustness and generalizability. However, we acknowledge that the results of predictive research with small sample sizes are exploratory in nature.(27) Nevertheless, external validation in various clinical settings, covering the heterogeneity among cases, is essential for clinical application. Given the high prediction accuracy of our method, further development using large databases, such as national or regional datasets, is expected and necessary.
Conclusions
Our study demonstrates the potential of ML to predict massive IBL during LT. These data-driven predictive models could revolutionize preoperative planning and intraoperative management, ultimately improving patient care and resource efficiency. Future research should focus on external validation using large datasets and the integration of these models into clinical practice.
Data Availability
Data cannot be shared publicly because of institutional policy and patient confidentiality. Data are available from the Jichi Medical University Institutional Ethics Committee (contact via https://www.jichi.ac.jp/kenkyushien/clinical/clinical_human/) for researchers who meet the criteria for access to confidential data.
Financial Disclosure
The authors received no specific funding for this work.
Competing Interests
The authors have declared that no competing interests exist.
Ethics statement
This study was approved by the Ethics Committee of Jichi Medical University (Ethics Committee Approval Case Number 20-008). Informed consent was obtained in the form of opt-out on our website (https://www.jichi.ac.jp/transplant/contents/disclosure.html), with the approval of the Ethics Committee of Jichi Medical University.
Data availability statement
The data used in this study are not publicly available due to institutional policy and patient confidentiality. As participant consent for public sharing was not obtained, the data cannot be shared openly. However, anonymized data may be available from the corresponding author (TW) upon reasonable request.
Figure Legends
Supplemental Figure 1: Representative performance of the binary logistic regression model at each stage of the backward stepwise feature selection process. This figure illustrates a single iteration out of 100 repetitions. AUC, area under the receiver operating characteristic curve.
Supplemental Figure 2: Model performance evaluation based on a feature selection process repeated 100 times. This figure presents the most accurate outcome observed among the selection processes conducted in each random split. AUC, area under the receiver operating characteristic curve.