- 1Department of Endocrinology, The First Medical Center of People's Liberation Army General Hospital, Beijing, China
- 2Department of Endocrinology, Central Hospital of Dalian University of Technology, Dalian, Liaoning, China
- 3Department of Neurology, Johns Hopkins Hospital, Baltimore, MD, United States
- 4Yidu Cloud Technology Inc., Beijing, China
- 5School of Software Technology, Dalian University of Technology, Dalian, Liaoning, China
Background: Accurately predicting intraoperative hemodynamic instability (HI) in patients with pheochromocytoma is essential for improving prognosis; however, clinically applicable large-sample, high-precision predictive models remain limited. This study develops and validates an ensemble learning (EL) model to predict HI risk.
Methods: This cohort study included a derivation cohort (n = 353) and an external validation cohort (n = 51), from January 2011 to February 2023. General clinical and intraoperative hemodynamic data were collected. Ensemble feature selection was used to identify key predictors. 5-fold cross-validation was repeated 1000 times to develop the EL model. Shapley Additive Explanations was used to analyze feature contributions, and the model was implemented as a web calculator. The primary outcome was the occurrence of intraoperative HI, evaluated by area under the curve (AUC), sensitivity, specificity, and calibration.
Results: Of 45 variables, tumor size, preoperative systolic blood pressure, age, fasting plasma glucose, and body mass index were top predictors. The developed EL model achieved AUC, sensitivity, and specificity values of 0.886, 0.776, and 0.836 and 0.744, 0.733, and 0.667 in Training set and external validations, respectively. Higher SBP (≥ 125 mmHg), larger tumor size (≥ 60 mm), older age (≥ 55 years), higher FPG (≥ 6 mmol/L), and BMI <22 or >30 kg/m² increased HI risk. The model demonstrated strong calibration and is accessible at http://60.205.91.235/.
Conclusions: This study identified five key predictors of intraoperative HI in patients with pheochromocytoma. The developed EL model provides an accurate, clinically applicable HI risk estimation tool, potentially improving clinical management.
1 Introduction
Pheochromocytoma is a rare neuroendocrine tumor characterized by excessive catecholamine production, storage, and secretion (1). Surgery currently remains the only effective treatment; however, several factors, such as general anesthesia during surgery, ventilator-mediated abdominal pressure fluctuations, and direct treatment of the tumor, can trigger catecholamine release, leading to intraoperative hypertensive crises. Severe cases may result in stroke, myocardial infarction, heart failure, multiple organ failure, or mortality (2–4). Advancements in preoperative preparation, anesthesia, and surgical techniques for pheochromocytoma have significantly reduced perioperative mortality rates from 20%–45% to 0%–2.9% (5–12). Despite these improvements, hemodynamic instability (HI) still affects 39.1–100% of patients (13–18).
Accurately predicting HI before surgery and implementing timely interventions can significantly reduce its incidence, mitigate complications, and improve surgical outcomes. While current clinical research primarily investigates risk factors for HI, few studies have focused on predictive models for perioperative HI risk or complications in pheochromocytoma cases (12, 14, 19). Existing studies are notably limited by small sample sizes, challenges in acquiring model variables, and the absence of external validation or advanced analytical approaches such as machine learning (ML). This underscores the need for data-driven research utilizing larger sample sizes and more sophisticated analytical techniques.
Compared to conventional statistical models, ML offers superior classification and regression capabilities, effectively handling complex data relationships while simulating the comprehensive diagnostic reasoning of clinicians. Ensemble learning (EL) further enhances predictive performance by integrating multiple ML models. It constructs an optimal predictive framework by leveraging the strengths of diverse algorithms, thereby improving stability and prediction accuracy.
This study employed an EL prediction model to analyze the risk factors for intraoperative HI in patients with pheochromocytoma and validated model performance on an external dataset. The findings provide a foundation for optimizing preoperative preparation, reducing intraoperative HI incidence, and ultimately improving patient outcomes.
2 Methods
2.1 Patients
This retrospective study included data from 439 patients with pheochromocytoma who underwent unilateral laparoscopic adrenalectomy at the First Medical Centre of the Chinese PLA General Hospital (Beijing, China) between January 2011 and November 2021. Additionally, for external validation, data from 61 patients with pheochromocytoma who underwent the same procedure at the affiliated Central Hospital of the Dalian University of Technology (Dalian, China) between January 2011 and February 2023 were included.
Ethical approval was obtained from the clinical research ethics committees of both hospitals (S2022-085–01 and YN2021-023-01). Due to the retrospective nature of the study, the requirement for informed consent was waived. This study adhered to the principles outlined in the Helsinki Declaration. As inclusion criteria, patients should have been diagnosed with unilateral pheochromocytoma of the adrenal gland, underwent laparoscopic adrenalectomy performed for an adrenal tumor, and received postoperative histopathological confirmation of pheochromocytoma. Exclusion criteria included patients with paraganglioma, malignant pheochromocytoma, a family history of pheochromocytoma-related diseases, or concurrent tumors. Preoperative management included the administration of alpha-blockers and, for some patients, a combination of beta-blockers or even calcium channel blockers for blood pressure control, along with fluid infusion for volume expansion.
Ultimately, 353 patients from the First Medical Centre of the Chinese PLA General Hospital were included in the modelling cohort, while 51 patients from the affiliated Central Hospital of the Dalian University of Technology comprised the external validation set.
2.2 Characteristics and outcomes
2.2.1 Case data collection
Clinical data, including symptoms, preoperative medication use, medical history, surgical history, and family history were retrospectively retrieved from electronic medical records of the hospital upon patient admission. Catecholamine and metabolite levels were measured using high-performance liquid chromatography for 24 h urinary catecholamines, including norepinephrine, epinephrine, and dopamine. High-performance liquid chromatography-tandem mass spectrometry was employed to measure catecholamine metabolites, including normetanephrine (NMN), metanephrine (MN), and 3-methoxytyramine. Details regarding testing times and equipment are in Supplementary Table 1. Imaging results from computed tomography (CT), magnetic resonance imaging (MRI), and positron emission tomography-CT scans were accessed via the hospital electronic records, documenting tumor location, morphology, and maximum diameter. The tumor size was defined based on the largest diameter identified in the imaging results.
2.2.2 Intraoperative hemodynamic data collection
Perioperative nursing records, temperature charts, and anesthesia records were reviewed to document intraoperative and perioperative complications, surgical techniques, surgery duration, intraoperative hemodynamics (blood pressure and heart rate [HR]), and fluid input and output (crystalloids, colloids, blood products, drainage, and urine output) during and after surgery.
Preoperative blood pressure and HR were measured on the morning of surgery. Intraoperative monitoring included electrocardiography and pulse oximetry, with blood pressure and HR automatically and continuously recorded every 5 minutes via a radial artery catheter. Based on prior studies (11, 20–22), intraoperative hemodynamic characteristics were evaluated using the cumulative frequency and duration of the following events: systolic blood pressure (SBP) ≥ 200 mmHg, SBP ≥ 180 mmHg, SBP ≥ 160 mmHg, SBP < 90 mmHg, diastolic blood pressure (DBP) < 50 mmHg, mean arterial pressure < 60 mmHg, tachycardia (HR ≥ 110 bpm), and bradycardia (HR ≤ 50 bpm).
2.3 Model development and evaluation
2.3.1 Model overview
In medical research, EL methodologies provide a refined approach to decision-making processes. Within this framework, the base-learner outcomes resemble the collective assessments of diverse medical practitioners regarding the health status of the patient, while the meta-learner integrates these varied assessments into a cohesive final evaluation. The EL predictive model in this study comprises three core components: feature engineering, base learning, and meta-learning.
2.3.2 Data processing
Data on 73 characteristics from each patient in the derivation cohort were collected. Categorical variables such as diabetes mellitus, hypertension, and the primary outcome variable (HI) were encoded as 0 (not occurring, negative) and 1 (occurring, positive). Owing to significant changes in catecholamine detection methods over the past decade, variables related to catecholamines could not be treated as continuous variables, which is consistent with previous pheochromocytoma-related studies (14). Instead, these variables were encoded as 0 (not elevated) and 1 (elevated) and subsequently combined into composite variables.
The missing proportions of features in the original dataset were analyzed. Features with missing proportions >10% were excluded, while those with <10% were imputed using multiple imputation methods.
Based on the complete dataset, the derivation cohort was randomly divided into a training set and an internal validation set in an 8:2 ratio. This process was repeated 1000 times for robust model evaluation.
2.4 Feature selection
To enhance feature screening applicability and identify relevant input features, we carried out ensemble feature selection (EFS) using random forest (RF), adaptive boosting (AdaBoost), gradient boosting decision tree (GBDT), and extremely randomized trees (Extra-Trees) algorithms (23). To minimize the impact of dataset partitioning on feature screening, a 5-fold cross-validation method was applied. The top five features were selected as inputs for the EL prediction model.
2.5 Selection of base and meta-learners
The optimal parameters of the RF, AdaBoost, GBDT, Extra-Trees, and Lasso models were determined via grid search. For learner selection, the top five features were used as inputs for the RF, AdaBoost, GBDT, Extra-Trees, and Lasso models, and all the features were used for the Lasso model. Model performance was evaluated using 5-fold cross-validation on the training set. Standard performance metrics, including the area under the curve (AUC), F1 score, sensitivity, specificity, accuracy, precision, positive predictive value, and negative predictive value, were employed.
2.6 Model evaluation
The top three models based on performance were selected as candidate-based learners, with the best-performing model serving as the meta-learner to construct an integrated prediction framework. Base learner outputs served as inputs to the meta-learner. The EL prediction models were evaluated on the training and internal validation sets, with the optimal one selected as the final EL prediction model. An external validation set was used to further validate the model. The Shapley Additive Explanations (SHAP) dependence method was used to analyze the relationships between features and their contributions to the HI risk (24). Additionally, a calibration curve assessed the relationship between predicted and actual values. Finally, the EL prediction model was presented as a web calculator.
2.7 Model construction and statistical analyzes
Continuous variables with normal distributions were compared using unpaired two-tailed t-tests, while non-normally distributed variables were analyzed using Wilcoxon rank-sum tests. Categorical variables were compared using Chi-squared tests. Model construction and evaluation metrics calculations were performed using Python 3.8 software and libraries, including “sklearn.ensemble,” “sklearn.model_selection,” and “sklearn.metrics.” One thousand random dataset partitions were generated using the train_test_split function in the sklearn.model_selection package, with the current time as the random seed. The statistical significance level was set to P < 0.05. Statistical analyzes were performed using R version 4.0.1 software (R Foundation for Statistical Computing, Vienna, Austria) and the “randomForest,” “ggplot2,” and “rms,” packages.
3 Results
3.1 Dataset
The derivation cohort included 353 patients (183, 51.8% positive). The dataset was randomly divided into training and internal validation sets of 282 and 71 patients, respectively. The external validation set included 51 patients (30, 58.8% positive). Table 1 presents the clinical and biological characteristics of pheochromocytoma patients with and without HI.
Table 1. Clinical and biological characteristics of patients with pheochromocytoma with and without hemodynamic instability (HI).
3.2 Missing data imputation
The dataset included 73 features, with missing data percentages ranging from 1.13% to 79.06%. Features with missing data rates exceeding 10% (28 features) were excluded, while those with missing rates below 10% (three features) were imputed using multiple imputation techniques. Supplementary Figure 1 shows the specific missing feature rates.
After imputation, no significant outliers were observed in the feature set. The collinearity correlation matrix heatmap in Supplementary Figure 2 shows the Pearson correlation coefficients between each pair of variables, with colors representing correlation strength.
3.3 EFS
The feature importance scores were summed and ranked (Figure 1), yielding importance scores for 45 features of each model. The top five features, preoperative SBP, tumor size, baseline fasting plasma glucose (FPG), age, and body mass index (BMI), were selected for training the learners.
Figure 1. Variables ranking based on ensemble feature selection (EFS) scores. The y-axis shows the top ten variables ranked by importance, while the x-axis shows their cumulative importance values calculated using ensemble feature selection methods such as Random Forest, AdaBoost, Gradient Boosting, and Extremely randomized trees (Extra-Trees).
3.4 Base and meta-learners
Among the evaluated models, Extra-Trees yielded the best overall performance, followed by the RF model, which employs a similar generation process. The five-feature Lasso method outperformed the all-feature Lasso approach and ranked third. These top three models served as candidate base learners and were combined into meta-learners with various combinations of base learners. The following EL prediction models were implemented: ET-RF, which employs Extra-Trees and RF as base learners, with Extra-Trees as the meta-learner; ET-LA, which uses Extra-Trees and Lasso as base learners, with Extra-Trees as the meta-learner; and ET-LA-RF, which integrates Extra-Trees, Lasso, and RF as base learners, with Extra-Trees as the meta-learner. All model performances were compared using the metrics listed in Supplementary Table 2.
3.5 Selection of prediction models
Three integrated learning prediction models were compared using 1000 randomly divided training and internal validation sets (Table 2). The ET-RF model had the highest AUC (0.886) and the best sensitivity, positive predictive value, and all evaluation indices. On the external validation set, ET-RF outperformed the other two EL prediction models in seven out of eight evaluation metrics (Table 3). Figure 2 shows the SHAP dependence plot for the ET-RF model, where each dot represents a patient. The plot depicts how the attributed importance of each variable changes with its value. SHAP values exceeding zero represent increasing HI risk. Key findings from the SHAP analysis reveal that Higher SBP (≥ 125 mmHg), larger tumor size (≥ 60 mm), older age (≥ 55 years), higher FPG (≥ 6 mmol/L), and BMI <22 or >30 kg/m² increased HI risk.
Table 2. Performance of ensemble learning prediction models on training set and internal validation set.
Figure 2. Shapley additive explanations (SHAP) dependence plots for clinical variables predicting hemodynamic instability in patients with pheochromocytoma. The plots show SHAP values for (A) Age, (B) Body Mass Index (BMI), (C) Tumor Size (maximum CT/MRI measurement), (D) Systolic Blood Pressure (SBP) before anesthesia, and (E) Fasting Plasma Glucose (FPG) before drug administration. Positive SHAP values indicate that higher variable values are associated with an increased risk of hemodynamic instability (HI), while negative values suggest a protective or reduced risk.
3.6 Web calculator
To facilitate clinical application, we developed an online web calculator tool for physicians. Physicians need only to input the tumor size, preoperative SBP, age, baseline FPG level, and BMI of the patient. The calculator estimates HI risk based on a selected threshold probability: high-risk estimate probability ≥0.5, low-risk estimate probability <0.5. The calculator is accessible at http://60.205.91.235/.
3.7 Evaluation of predictive models
Supplementary Figure 3 shows the calibration curves of the three models for the internal (left) and external validation sets (right). The ET-RF prediction model, which uses Extra-trees and RF as base learners and Extra-trees as a meta-learner, closely matched the actual outcomes.
4 Discussion
In this study, we selected a research center with high surgical volume for pheochromocytoma, utilizing a large sample size and rigorous design that included only patients undergoing laparoscopic surgery. Our resulting predictive model demonstrated strong discriminative capability, enabling clinicians to more effectively implement targeted interventions and enhance patient outcomes.
Previous studies indicate that intraoperative HI can lead to perioperative complications and mortality (8, 11, 19). Although our study aligns with prior findings of zero mortality (5, 8, 9), HI resulted in significant adverse events within our cohort, including cardiovascular complications in eight patients, massive blood loss in five, and cerebrovascular complications in three. Additionally, 16 patients required intensive care due to critical HI episodes, further underscoring the clinical significance of this outcome.
While previous studies have investigated predictive models for intraoperative HI in pheochromocytoma patients (12, 14, 19), many have limited generalizability due to small sample sizes, insufficient external validation, difficulties in acquiring comprehensive model features, and neglecting variability in surgical approaches.
Our study addresses these gaps by integrating a robust sample size, comprehensive patient characteristics, and rigorous methodological standards. Our model has several advantages. First, EFS integrates multiple feature selection algorithms, reducing bias from single-algorithm use (25). Second, previous ML models, like Lasso or linear regression, inadequately capture complex feature interactions (19, 26). T We addressed these limitations by employing ensemble learning (EL), involving feature engineering, base learners, and meta-learners to generate predictions. This integrated approach outperforms single models, providing more accurate clinical predictions. Furthermore, given the established influence of anesthesia and surgical team experience on HI risk (27), we conducted external validation with data from multiple institutions, confirming robust predictive performance and broader clinical applicability (25).
This study reveals that age, BMI, tumor size, preoperative SBP, and baseline FPG are independent risk factors for intraoperative HI in patients with pheochromocytoma. Unlike previous models (12, 14, 19), our analysis using SHAP dependence plots revealed significant associations: a higher preoperative SBP (≥125 mmHg) notably increased HI risk, reinforcing prior findings that optimal preoperative blood pressure management is critical for minimizing HI, related complications, and mortality (16, 28, 29). Current guidelines consistently advise maintaining blood pressure within safe ranges to reduce HI risk (3).
Tumor size is a common risk factor for HI (4, 5, 7, 8, 15). Larger pheochromocytomas typically produce and metabolize more catecholamines, resulting in severe hypertension during adrenalectomy (30). Tumor size also emerged as a critical factor; risk transitioned from negative to positive at a threshold of 55–60 mm. Patients with tumors ≥60 mm are known to experience increased intraoperative blood loss compared to those with smaller tumors (31). Larger pheochromocytomas correlate positively with HI risk, partly due to greater intraoperative blood loss from the surgical procedure and rapid reductions in serum catecholamine levels and circulating blood volume (32, 33).
Four prior studies have utilized preoperative data to predict intraoperative HI, identifying age, tumor size, and BMI as significant independent risk factors (12, 14, 19, 29). Additionally, Previous research on risk factors also suggests a correlation between age and HI (16). Our SHAP dependence analysis corroborates these findings, demonstrating that age above 55 years is associated with increased HI risk, consistent with previous reports identifying age >45 years as an independent risk factor for intraoperative hemodynamic instability (14, 34). Furthermore, BMI is both a risk factor for HI and an independent risk factor for severe postoperative and cardiovascular complications (35, 36). The SHAP dependence plot showed that a BMI outside the specific range 22–30 Kg/m2, whether higher or lower, was associated with an increased HI risk. Specifically, lower BMI is often associated with a reduced circulating capacity, making patients more prone to HI.
Glucose regulation plays a pivotal role as well. Patients with concurrent diabetes and pheochromocytoma exhibit elevated perioperative HI risk (17, 18). Our findings indicate increased HI risk when baseline FPG elevated. Insulin resistance in these patients, driven by excessive catecholamine secretion, significantly contributes to HI (18, 37, 38). Additionally, elevated blood glucose levels may reflect underlying cardiovascular dysfunction, such as myocardial impairment and arteriosclerosis, thereby exacerbating intraoperative hemodynamic instability (39, 40).
This study has several limitations. First, it was retrospective, and prospective studies are necessary to further evaluate the applicability of the model. Second, because the retrospective period spanned many years, newer detection methods for catecholamines have since emerged, and measurement standards vary across institutions. Variables related to catecholamines and their metabolites were initially converted into categorical variables but were eventually excluded from feature selection due to substantial missing data. Nevertheless, our model still demonstrated excellent predictive performance for hemodynamic instability. Moreover, since catecholamine and their metabolites measurements are complex and not routinely performed in all centers, the ability of our model to achieve high accuracy without relying on these variables enhances its generalizability and practical utility in clinical settings.
In conclusion, this study analyzed a large sample of patients with varied characteristics and identified five key risk factors for intraoperative HI in pheochromocytoma patients. Using an EL-based prediction model, we developed and validated a risk model for intraoperative HI utilizing internal and external datasets. Additionally, we created a user-friendly and accurate web calculator to assess individualized HI risk prior to surgery. This tool enables clinicians to proactively predict HI risk in patients undergoing laparoscopic unilateral pheochromocytoma resection, enabling timely interventions to reduce HI incidence and improve patient outcomes.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
Ethical approval was obtained from the clinical research ethics committees of both hospitals (S2022-085-01 and YN2021-023-01). The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because Due to the retrospective nature of the study, the requirement for informed consent was waived.
Author contributions
YL: Conceptualization, Investigation, Validation, Writing – original draft, Writing – review & editing. CL: Writing – review & editing, Investigation, Methodology. SL: Investigation, Writing – review & editing, Resources, Validation. LZh: Data curation, Formal analysis, Methodology, Software, Visualization, Writing – original draft. RL: Data curation, Formal analysis, Methodology, Software, Visualization, Writing – original draft, Validation. KC: Conceptualization, Investigation, Writing – review & editing. LZa: Data curation, Investigation, Writing – review & editing, Project administration. WG: Writing – review & editing, Conceptualization, Resources. YM: Conceptualization, Resources, Writing – review & editing. ZL: Project administration, Supervision, Writing – review & editing. ZG: Conceptualization, Investigation, Project administration, Supervision, Writing – review & editing. JD: Conceptualization, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This study was supported by the Outstanding Young Science and Technology Talent in Dalian (2024RY001), the Peak Climbing Plan Project of Dalian Municipal Central Hospital (2024ZZ035, 2023ZZ029, 2021ZZ006) and the Dalian Science and Technology Innovation Fund (2023JJ13SN060).
Acknowledgments
We thank the surgical and anesthesia teams of the First Medical Center, Chinese PLA General Hospital and Central Hospital of Dalian University of Technology.
Conflict of interest
Author CL was employed by the company Yidu Cloud Technology Inc.
The remaining 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 author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.1670909/full#supplementary-material
References
1. Lenders JW, Eisenhofer G, Mannelli M, and Pacak K. Phaeochromocytoma. Lancet. (2005) 366:665–75. doi: 10.1016/S0140-6736(05)67139-5
2. Lentschener C, Gaujoux S, Tesniere A, and Dousset B. Point of controversy: perioperative care of patients undergoing pheochromocytoma removal-time for a reappraisal? Eur J Endocrinol. (2011) 165:365–73. doi: 10.1530/EJE-11-0162
3. Lenders JW, Duh QY, Eisenhofer G, Gimenez-Roqueplo AP, Grebe SK, Murad MH, et al. Pheochromocytoma and paraganglioma: an endocrine society clinical practice guideline. J Clin Endocrinol Metab. (2014) 99:1915–42. doi: 10.1210/jc.2014-1498
4. Bruynzeel H, Feelders RA, Groenland TH, van den Meiracker AH, van Eijck CH, Lange JF, et al. Risk Factors for Hemodynamic Instability during Surgery for Pheochromocytoma. J Clin Endocrinol Metab. (2010) 95:678–85. doi: 10.1210/jc.2009-1051
5. Kiernan CM, Du L, Chen X, Broome JT, Shi C, Peters MF, et al. Predictors of hemodynamic instability during surgery for pheochromocytoma. Ann Surg Oncol. (2014) 21:3865–71. doi: 10.1245/s10434-014-3847-7
6. Livingstone M, Duttchen K, Thompson J, Sunderani Z, Hawboldt G, Sarah Rose M, et al. Hemodynamic stability during pheochromocytoma resection: lessons learned over the last two decades. Ann Surg Oncol. (2015) 22:4175–80. doi: 10.1245/s10434-015-4519-y
7. Buitenwerf E, Osinga TE, Timmers H, Lenders JWM, Feelders RA, Eekhoff EMW, et al. Efficacy of α-blockers on hemodynamic control during pheochromocytoma resection: a randomized controlled trial. J Clin Endocrinol Metab. (2020) 105:2381–91. doi: 10.1210/clinem/dgz188
8. Pacak K. Preoperative management of the pheochromocytoma patient. J Clin Endocrinol Metab. (2) 92:4069–79. doi: 10.1210/jc.2007-1720
9. Wachtel H, Kennedy EH, Zaheer S, Bartlett EK, Fishbein L, Roses RE, et al. Preoperative Metyrosine Improves Cardiovascular Outcomes for Patients Undergoing Surgery for Pheochromocytoma and Paraganglioma. Ann Surg Oncol. (2015) 22 Suppl 3:S646–54. doi: 10.1245/s10434-015-4862-z
10. Lenders JWM, Kerstens MN, Amar L, Prejbisz A, Robledo M, Taieb D, et al. Genetics, diagnosis, management and future directions of research of phaeochromocytoma and paraganglioma: a position statement and consensus of the Working Group on Endocrine Hypertension of the European Society of Hypertension. J Hypertens. (2020) 38:1443–56. doi: 10.1097/HJH.0000000000002438
11. Brunaud L, Nguyen-Thi PL, Mirallie E, Raffaelli M, and Bresler L. Predictive factors for postoperative morbidity after laparoscopic adrenalectomy for pheochromocytoma: a multicenter retrospective analysis in 225 patients. Surg Endoscopy. (2016) 30:1051–9.
12. Zhao Y, Fang L, Cui L, and Bai S. Application of data mining for predicting hemodynamics instability during pheochromocytoma surgery. BMC Med Inf decision making. (2020) 20:1–13.
13. Kim WW, Kim DH, Cho JW, Rah CS, Lee YM, Chung KW, et al. The association between the type of anesthesia and hemodynamic instability during pheochromocytoma surgery: a retrospective cohort study. Surg Endosc. (2022) 36:5491–500. doi: 10.1007/s00464-021-08910-3
14. Zhang Z, Ye Y, Yu J, Liao S, Pan W, Guo Y, et al. A Nomogram for Predicting Intraoperative Hemodynamic Instability in Patients With Pheochromocytoma. Front Endocrinol (Lausanne). (2021) 12:787786. doi: 10.3389/fendo.2021.787786
15. Urabe F, Kimura S, Iwatani K, Takahashi K, Ito K, Tashiro K, et al. Risk Factors for Perioperative Hemodynamic Instability in Pheochromocytoma: A Systematic Review and Meta-Analysis. J Clin Med. (2021) 10. doi: 10.3390/jcm10194531
16. Tariel F, Dourmap C, Prudhomme T, Hascoet J, Soulie M, Moreau B, et al. Adrenalectomy for pheochromocytoma: complications and predictive factors of intraoperative hemodynamic instability. Am Surg. (2022) 27:31348221135774. doi: 10.1177/00031348221135774
17. Takeda T, Hakozaki K, Yanai Y, Masuda T, Yasumizu Y, Tanaka N, et al. Risk factors for haemodynamic instability and its prolongation during laparoscopic adrenalectomy for pheochromocytoma. Clin Endocrinol (Oxf). (2021) 95:716–26. doi: 10.1111/cen.14557
18. Pisarska-Adamczyk M, Zawadzka K, Wieckowski K, Przeczek K, Major P, Wysocki M, et al. Risk factors for hemodynamic instability during laparoscopic pheochromocytoma resection: a retrospective cohort study. Gland Surg. (2021) 10:892–900. doi: 10.21037/gs-20-783
19. Bai S, Wu B, Yao Z, Zhu X, Jiang Y, and Wang H. Development and validation of a clinical model to predict intraoperative hemodynamic instability in patients with pheochromocytomas surgery. Endocr J. (2020) 67:81–9. doi: 10.1507/endocrj.EJ19-0278
20. Lafont M, Fagour C, Haissaguerre M, Darancette G, Wagner T, Corcuff JB, et al. Per-operative hemodynamic instability in normotensive patients with incidentally discovered pheochromocytomas. J Clin Endocrinol Metab. (2015) 100:417–21. doi: 10.1210/jc.2014-2998
21. Brunaud L, Boutami M, Nguyen-Thi PL, Finnerty B, Germain A, Weryha G, et al. Both preoperative alpha and calcium channel blockade impact intraoperative hemodynamic stability similarly in the management of pheochromocytoma. Surgery. (2014) 156:1410–7. doi: 10.1016/j.surg.2014.08.022
22. Gaujoux S, Bonnet S, Lentschener C, Thillois JM, Duboc D, Bertherat J, et al. Preoperative risk factors of hemodynamic instability during laparoscopic adrenalectomy for pheochromocytoma. Surg Endosc. (2016) 30:2984–93. doi: 10.1007/s00464-015-4587-x
23. Scholten A, Vriens MR, Cromheecke GJ, Borel Rinkes IH, and Valk GD. Hemodynamic instability during resection of pheochromocytoma in MEN versus non-MEN patients. Eur J Endocrinol. (2011) 165:91–6. doi: 10.1530/EJE-11-0148
24. Lundberg SM and Lee S-I. A unified approach to interpreting model predictions. Adv Neural Inf Process Syst. (2017) 30.
25. Ramspek CL, Jager KJ, Dekker FW, Zoccali C, and van Diepen M. External validation of prognostic models: what, why, how, when and where? Clin Kidney J. (2021) 14:49–58. doi: 10.1093/ckj/sfaa188
26. Li M, Sun K, Dai W, Xiang W, Zhang Z, Zhang R, et al. Preoperative prediction of peritoneal metastasis in colorectal cancer using a clinical-radiomics model. Eur J Radiol. (2020) 132:109326. doi: 10.1016/j.ejrad.2020.109326
27. Vorselaars W, Postma EL, Mirallie E, Thiery J, Lustgarten M, Pasternak JD, et al. Hemodynamic instability during surgery for pheochromocytoma: comparing the transperitoneal and retroperitoneal approach in a multicenter analysis of 341 patients. Surgery. (2018) 163:176–82. doi: 10.1016/j.surg.2017.05.029
28. Araujo-Castro M, Garcia Centeno R, López-García MC, Lamas C, Álvarez-Escolá C, Calatayud Gutiérrez M, et al. Risk factors for intraoperative complications in pheochromocytomas. Endocr Relat Cancer. (2021) 28:695–703. doi: 10.1530/ERC-21-0230
29. Liu Y, Liu C, Wang Y, Li S, Li X, Liu X, et al. Nomogram for predicting intraoperative hemodynamic instability in patients with normotensive pheochromocytoma. J Clin Endocrinol Metab. (2023) 108:1657–65. doi: 10.1210/clinem/dgad024
30. Guerrero MA, Schreinemakers JM, Vriens MR, Suh I, Hwang J, Shen WT, et al. Clinical spectrum of pheochromocytoma. J Am Coll Surgeons. (2009) 209:727–32.
31. Natkaniec M, Pedziwiatr M, Wierdak M, Major P, Migaczewski M, Matlok M, et al. Laparoscopic transperitoneal lateral adrenalectomy for large adrenal tumors. Urol Int. (2016) 97:165–72. doi: 10.1159/000444146
32. Natkaniec M, Pedziwiatr M, Wierdak M, Bialas M, Major P, Matlok M, et al. Laparoscopic adrenalectomy for pheochromocytoma is more difficult compared to other adrenal tumors. Wideochir Inne Tech Maloinwazyjne. (2015) 10:466–71. doi: 10.5114/wiitm.2015.52869
33. Bozkurt IH, Arslan M, Yonguc T, Degirmenci T, Koras O, Gunlusoy B, et al. Laparoscopic adrenalectomy for large adrenal masses: Is it really more complicated? Kaohsiung J Med Sci. (2015) 31:644–8. doi: 10.1016/j.kjms.2015.09.005
34. Jiang M, Ding H, Liang Y, Tang J, Lin Y, Xiang K, et al. Preoperative risk factors for haemodynamic instability during pheochromocytoma surgery in Chinese patients. Clin Endocrinology. (2018) 88:498–505. doi: 10.1111/cen.13544
35. Bai S, Yao Z, Zhu X, Li Z, Jiang Y, Wang R, et al. Risk factors for postoperative severe morbidity after pheochromocytoma surgery: A single center retrospective analysis of 262 patients. Int J Surg. (2018) 60:188–93. doi: 10.1016/j.ijsu.2018.11.019
36. Stolk RF, Bakx C, Mulder J, Timmers HJ, and Lenders JW. Is the excess cardiovascular morbidity in pheochromocytoma related to blood pressure or to catecholamines? J Clin Endocrinol Metab. (2013) 98:1100–6. doi: 10.1210/jc.2012-3669
37. Petrák O, Klímová J, Mráz M, Haluzíková D, Doležalová RP, Kratochvílová H, et al. Pheochromocytoma with adrenergic biochemical phenotype shows decreased GLP-1 secretion and impaired glucose tolerance. J Clin Endocrinol Metab. (2020) 105. doi: 10.1210/clinem/dgaa154
38. Komada H, Hirota Y, So A, Nakamura T, Okuno Y, Fukuoka H, et al. Insulin secretion and insulin sensitivity before and after surgical treatment of pheochromocytoma or paraganglioma. J Clin Endocrinol Metab. (2017) 102:3400–5. doi: 10.1210/jc.2017-00357
39. Punthakee Z, Iglesias PP, Alonso-Coello P, Gich I, India I, Malaga G, et al. Association of preoperative glucose concentration with myocardial injury and death after non-cardiac surgery (GlucoVISION): a prospective cohort study. Lancet Diabetes Endocrinol. (2018) 6:790–7. doi: 10.1016/S2213-8587(18)30205-5
Keywords: ensemble learning, prediction model, hemodynamic instability, intraoperative risk prediction, pheochromocytoma, web-based calculator
Citation: Liu Y, Liu C, Li S, Zhao L, Lin R, Chen K, Zang L, Gu W, Mu Y, Lyu Z, Gao Z and Dou J (2025) Ensemble learning prediction model for intraoperative hemodynamic instability in patients with pheochromocytoma. Front. Endocrinol. 16:1670909. doi: 10.3389/fendo.2025.1670909
Received: 22 July 2025; Accepted: 31 October 2025;
Published: 12 December 2025.
Edited by:
Sergei Tevosian, University of Florida, United StatesReviewed by:
Livia Palmieri, Sapienza University of Rome, ItalyMaria Eliana Rapicano, Federico II University Hospital, Italy
Copyright © 2025 Liu, Liu, Li, Zhao, Lin, Chen, Zang, Gu, Mu, Lyu, Gao and Dou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jingtao Dou, amluZ3Rhb2RvdUAxNjMuY29t; Zhengnan Gao, Z2FvMjAwOEAxNjMuY29t
†These authors have contributed equally to this work
‡ORCID: Zhengnan Gao, orcid.org/0000-0002-0077-8758
Jingtao Dou, orcid.org/0000-0003-0375-6527
Chao Liu4