Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Endocrinol., 05 September 2025

Sec. Thyroid Endocrinology

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

Machine learning-based predictive model for hungry bone syndrome following parathyroidectomy in secondary hyperparathyroidism

Yalin ChaiYalin ChaiNan YuanNan YuanJiaming YinJiaming YinBing ShenBing ShenLijie SunLijie SunLin ZhangLin ZhangLe YinLe YinXuan WangXuan WangFeng Luo*Feng Luo*Congjuan Luo*&#x;Congjuan Luo*†
  • The Affiliated Hospital of Qingdao University, Qingdao, China

Objective: To develop an interpretable machine learning model for predicting hungry bone syndrome (HBS) risk following parathyroidectomy in secondary hyperparathyroidism (SHPT) patients.

Methods: This retrospective study analyzed 181 SHPT patients who underwent parathyroidectomy at the Affiliated Hospital of Qingdao University (2015 - 2025). Participants were randomly divided into a training group (70%) and a validation group (30%). From 46 candidate variables, five key predictors were selected through logistic regression and Boruta algorithm. Seven machine learning models were trained, evaluated by ROC curves, calibration curves, and decision curve analysis (DCA). Model interpretability was quantified via SHapley Additive exPlanations (SHAP).

Results: The XGBoost algorithm demonstrated excellent predictive performance, with an AUC of 0.878 (95% CI: 0.779 - 0.973) and an F1 score of 0.871 for the validation cohort. The key predictors included preoperative parathyroid hormone (Pre-PTH), the percentage decay between Pre-PTH and PTH at skin closure (%PTH), alkaline phosphatase, serum calcium, and age. Additionally, we designed a web application to estimate HBS risk.

Conclusions: This interpretable machine-learning model is effective in predicting the risk of HBS in SHPT patients after parathyroidectomy, thereby providing guidance for postoperative surveillance strategies.

1 Introduction

Secondary hyperparathyroidism (SHPT), a prevalent complication of chronic kidney disease (CKD), contributes to bone lesions, vascular calcification, and elevated risks of cardiovascular events and mortality (13). Epidemiological studies indicate that SHPT affects 20%-80% of CKD patients, with prevalence rates correlating with disease severity and dialysis duration (4). Current therapeutic strategies encompass vitamin D analogs (e.g., calcitriol) (5), calcimimetics (e.g., cinacalcet), and phosphate binders (6). For refractory cases, parathyroidectomy (PTX) remains the definitive intervention per KDIGO guidelines, particularly in CKD patients with persistent hyperparathyroidism (intact parathyroid hormone [iPTH] >800 pg/mL) (7, 8). PTX has been shown to improve survival rates by 15%-57% in dialysis-dependent patients (9) and alleviate symptoms such as pruritus, bone pain, and fracture risk (1013).

Postoperative hungry bone syndrome (HBS), a complication affecting 25%-75% of PTX cases (1416), is clinically characterized by prolonged hypocalcemia (corrected serum calcium <2.1 mmol/L for >4 days) (17, 18). The pathophysiology involves accelerated bone remodeling under chronic PTH stimulation, followed by abrupt mineralization and calcium influx into osteoid tissue after rapid postoperative PTH decline, resulting in severe hypocalcemia (19, 20). Elevated preoperative PTH levels are strongly associated with HBS development (2123). Notably, recent studies propose that %PTH (the percentage decay between Pre-PTH and PTH at skin closure) predicts hypocalcemia following thyroidectomy (24). However, this metric has not yet been evaluated as a predictor of HBS in SHPT patients undergoing PTX. We investigated the relationship between %PTH and postoperative HBS and used machine learning to develop a predictive model for the occurrence of HBS after parathyroidectomy in SHPT patients.

2 Methods

2.1 Patients and designs

This retrospective cohort study enrolled patients diagnosed with secondary hyperparathyroidism who underwent parathyroidectomy at the Affiliated Hospital of Qingdao University between January 2015 and May 2025. Inclusion criteria required comprehensive preoperative clinical evaluation and documented serum calcium measurements within 72 hours postoperatively. The study protocol was reviewed and approved by the Ethics Committee of the Affiliated Hospital of Qingdao University (Approval No. QYFY WZLL 29980). Sample size determination adhered to the 10 events per variable (EPV) criterion, ensuring that the minimum number of outcomes for binary classifications exceeded 10-fold the number of independent variables in the predictive model.

Surgical indications for parathyroidectomy were defined as meeting at least one of the following criteria:(1) Severe SHPT (persistent intact parathyroid hormone [iPTH] >800 pg/mL) refractory to pharmacological therapy with calcitriol or vitamin D analogs;(2)Severe SHPT accompanied by hyperphosphatemia (serum phosphate >2.0 mmol/L);(3)Severe SHPT with clinically significant symptoms (intractable pruritus and/or bone pain);(4)Radiologically confirmed parathyroid hyperplasia (maximum gland diameter >1.0 cm).Exclusion criteria comprised:(1)Concurrent hepatobiliary or pancreatic disorders;(2)Cognitive impairment impairing informed consent or follow-up compliance;(3)Recurrent SHPT with prior parathyroidectomy history. The flow chart of the study is shown in Figure 1.

Figure 1
Flowchart showing the patient selection process. Initially, 204 patients with secondary hyperparathyroidism who underwent parathyroidectomy from January 2015 to May 2025 were considered. Exclusions included one patient with cognitive impairment, four with liver, biliary, or pancreatic diseases, two with a history of PTX surgery, and sixteen with extensive missing data, leaving 181 enrolled patients. These were divided into two groups: HBS group (n=105) and Non-HBS group (n=76).

Figure 1. The patient flow chart in our study.

2.2 Data collection and definitions

Demographic and clinical characteristics were systematically collected, including: General characteristics: Sex, age, height, body mass index (BMI), dialysis modality (hemodialysis/peritoneal dialysis), dialysis vintage, and clinical manifestations (bone pain, height reduction, pruritus); Preoperative management: Use of calcimimetics (Cinacalcet), vitamin D analogs (Calcitriol, Paricalcitol, Alfacalcitol), phosphate binders(Lanthanum carbonate), and medication for osteoporosis (Bisphosphonates); Laboratory parameters: Pre-PTH, PTH at skin closure, calcium, potassium, phosphate, albumin, alkaline phosphatase (ALP), cystatin C, creatinine, urea, uric acid, triglycerides, total cholesterol, low-density lipoprotein (LDL), high-density lipoprotein (HDL), and hemoglobin; Surgery-related factors: Procedure type, number of excised parathyroid glands; Comorbidities:hypertension, diabetes mellitus, coronary artery disease, cerebrovascular disease, osteoporosis, fracture.

Blood samples were obtained within 72 hours postoperatively. Patients were stratified into two groups based on corrected serum calcium levels:HBS group: Corrected calcium <2.1 mmol/L; Non-HBS group: Corrected calcium ≥2.1 mmol/L. Corrected calcium formula: Corrected Ca (mmol/L) = ionized Ca + (40 − serum albumin [g/L]) × 0.02.The percentage decay between Pre-PTH and PTH at skin closure(%PTH): [(Preoperative PTH – PTH at skin closure)/Preoperative PTH] × 100%. Blood was collected uniformly 20 minutes after total thyroidectomy (before skin closure) to control half-life effects.

2.3 Calcium management

Heparin-free dialysis was implemented during the final preoperative session and the initial postoperative session. Postoperative calcium monitoring followed a tiered protocol. The average medication dosage for managing hypocalcemia after parathyroidectomy varies among patients, but typically involves a combination of intravenous and oral calcium supplements. The recommended dose of intravenous calcium gluconate is 2 mg/kg/h, with the median total dose reported in studies being approximately 8.2g (Interquartile Range 6.1 - 10.3g). Oral calcium supplements usually consist of 3 - 6g elemental calcium daily (7.5 - 15g calcium carbonate), which should be combined with active vitamin D (e.g., calcitriol 0.5 - 2.0 μg/day). Discharge criteria included: Stable corrected calcium within normal range; Absence of hypocalcemia-related symptoms; No surgery-associated complications.

2.4 Statistical analysis and model development

All statistical analyses were conducted using R software (version 4.3.1, R Foundation for Statistical Computing, Austria), with a significance level set at P < 0.05.The normality test and the chi-square test were used for the measurement data, and the normally distributed data were expressed as mean ± standard deviation, and the t-test was used for comparison between groups; the non-normally distributed data were expressed as median (interquartile spacing), and the Mann-Whitney U-test was used for comparison between groups; the counting data were expressed as percentage (%), and the χ2-test or Fisher’s exact probability method was used for comparison between groups. ‘s exact probability method. MissForest (random forest interpolation technique) was used to interpolate the data, and data with more than 20% missing were discarded. The training set and test set are separated and interpolated respectively to avoid information leakage.

Use LASSO regression, Boruta method and logistic regression to select the final feature variables used in the model. In logistic analysis, variables with p < 0.05 are considered as potential risk factors. The LASSO method selects features and reduces dimension by narrowing the coefficients, retaining the features with large contributions and eliminating redundant features. Boruta is a feature selection method that determines the importance of each variable by comparing its Z-score. Three methods were selected to select the common characteristic variables as the final variables of the model. This method improves the accuracy of the model, reduces the risk of overfitting and excludes irrelevant predictors. Machine learning is more advantageous for classification and prediction and can handle features more efficiently than traditional statistical methods (22, 25, 26). Data were trained on the following seven ML models: Logistic Regression (LR), Adaptive Boosting (AdaBoost), Support Vector Machine (SVM), eXtreme Gradient Boosting (XGBoost), Categorical Boosting (CatBoost), K-Nearest Neighbors (KNN), and Neural Network (NN). We used grid search to optimize the parameters. The performance of the prediction model was evaluated by the AUC of the ROC curve, calibration curve, specificity, accuracy, recall, F1 score. Decision curve analysis (DCA) was plotted to assess clinical utility. For the most performance-driven diagnostic models, we revalidated their generalization capabilities using 10-fold cross-validation to prevent overfitting. The SHAP feature importance ranking and SHAP beeplot visually demonstrate each feature’s contribution to the final prediction, while the SHAP plot intuitively visualizes how different features influence individual predictions.

3 Results

3.1 Characteristics of the participants

From January 2015 to May 2025, 181 eligible secondary hyperparathyroidism (SHPT) patients undergoing parathyroidectomy (PTX) were included in the final analysis. These patients were randomly allocated into a training cohort (n=128) and a validation cohort (n=53) at a 7:3 ratio. Among the 181 participants, 105 (58%) developed postoperative hungry bone syndrome (HBS). The overall missing data situation is shown in Supplementary Table 1. Compared with the non-HBS group, patients in the HBS group had a significantly lower median age (P <0.001). Regarding comorbidities, the HBS group showed higher prevalence of osteoporosis (P = 0.017). In terms of preoperative medication use, the HBS group demonstrated higher usage rates of paricalciferol (P = 0.037) and bisphosphonates (P = 0.040). For preoperative symptoms, the HBS group exhibited higher proportions of osteodynia, eight reduction, and pruritus. These findings are detailed in Table 1. In terms of surgical methods, the proportion of SPTX was lower in HBS group (P <0.05), and the proportion of TPTX and TPTX+AT was higher. The proportion of HBS group with ≥4 parathyroid glands removed was higher, and the proportion of <4 parathyroid glands removed was lower. The results are shown in Table 2. Laboratory data comparison (Table 3) showed that compared with the non-HBS group, the HBS group had a higher percentage of Ca, Bun, Pre-PTH, PTH at skin closure, %PTH than the non-HBS group (P <0.05).

Table 1
www.frontiersin.org

Table 1. Comparison of General characteristics in the Non-HBS and HBS groups.

Table 2
www.frontiersin.org

Table 2. Comparison of Surgery-related factors in the Non-HBS and HBS groups.

Table 3
www.frontiersin.org

Table 3. Comparison of Biochemical indicators in the Non-HBS and HBS groups.

3.2 Feature selection

The selection of variables was based on the overlapping results of logistic regression, Lasso regression and Boruta’s algorithm. Univariate analysis showed (Table 4) that %PTH had a p-value <0.001, indicating that %PTH was significantly associated with postoperative HBS. In addition, Pre-PTH (p<0.001)), age (p<0.001), Ca (p=0.014), ALP (p<0.001), TPTX (p=0.045), SPTX (p=0.029) and the number of parathyroid glands resected (p=0.004) were significantly different between the HBS and non-HBS groups. All the above significant parameters were included in the multifactorial logistic regression analysis. In the LASSO regression analysis, variable coefficients are shown in Figure 2A, while the relationship between regularization parameter (λ) and mean cross-validation error (CVM) is illustrated in Figure 2B. The five variables identified through LASSO regression as strongly associated with HBS include Pre-PTH, ALP, %PTH, Age, and Ca. The regression coefficients for these variables are detailed in Supplementary Table 2. The Boruta method identified variables including Pre-PTH, ALP,%PTH, Age, PTH at skin closure, number of parathyroid glands removed, and Ca (Figure 2C). Ultimately, these five variables—Pre-PTH, ALP, %PTH, Age, and Ca—were selected for subsequent analyses.

Figure 2
Panel A shows a plot of binomial deviance against log(lambda), with a red curve and error bars illustrating the relationship. Panel B displays a coefficient path plot with multiple colored lines against log(lambda). Panel C is a boxplot comparing different variables, with red, blue, yellow, and green boxes representing different categorical data points on the x-axis.

Figure 2. Feature selection. (A) The relationship between Lambda (regularization parameter) and CVM (mean cross validation error) in Lasso regression; (B) Lasso regression Lambda and Coefficients Plot. (C) Feature selection based on Boruta principle. Green boxes indicate important variables, red boxes indicate unimportant variables.

Table 4
www.frontiersin.org

Table 4. Univariate and multivariate logistic regression analysis.

3.3 Performance comparison of different ML methods

Following feature selection, seven machine learning algorithms—Logistic Regression (LR), Adaptive Boosting (AdaBoost), Support Vector Machine (SVM), eXtreme Gradient Boosting (XGBoost), Categorical Boosting (CatBoost), K-Nearest Neighbors (KNN), and Neural Network (NN)—were trained to predict postoperative HBS.

To compare the need to include %PTH, we compared the difference in AUC values between whether or not to exclude %PTH. The results showed that the AUC value of the model including %PTH was higher (Supplementary Table 3), so we included %PTH for further analysis.

Among machine learning models, the XGBoost model (AUC = 0.878) showed the best performance, followed by LR (AUC = 0.876), CatBoost (AUC = 0.874), KNN (AUC = 0.869), SVM (AUC = 0.865), NN (AUC = 0.833) and AdaBoost (AUC = 0.821). (Figure 3A). Classification in the validation cohort outperformed all models with an F1 score of 0.871 and an accuracy of 0.849. Detailed performance metrics for all models, including specificity, precision, and accuracy, are summarized in Table 5. The calibration curves show a very good agreement between predicted probabilities and observations (Figure 3B). In addition, decision curve analysis (DCA) confirmed its optimal net benefit across clinically relevant probability thresholds (Figure 3C).

Figure 3
Three-panel data visualization comparing machine learning models. Panel A: ROC curves show sensitivity versus 1-specificity for models like Logistic, SVM, Neural Network, Xgboost, KNN, Adaboost, and CatBoost, with AUC values indicated. Panel B: Calibration plots depicting observed event percentage against bin midpoint with models following a diagonal reference line. Panel C: Decision curves showing standardized net benefit against the high-risk threshold for the models. Each line is color-coded per model.

Figure 3. Predictive performance of the model in validation cohorts. (A) Receiver Operating Characteristic (ROC) curves for the validation cohorts. The curves illustrate the discriminatory ability of different predictive models, with the area under the ROC curve (AUC) values displayed for each model: Logistic Regression (LR), eXtreme Gradient Boosting (XGBoost), Adaptive Boosting (AdaBoost), Support Vector Machine (SVM), K-Nearest Neighbor (KNN), Gradient Boosting Machine (GBM), Categorical Boosting (CatBoost). (B) Calibration Curve. Exhibited excellent alignment between predicted probabilities and observed outcomes. (C) Decision Curve Analysis (DCA). The curves show the net benefit of each model across various threshold probabilities.

Table 5
www.frontiersin.org

Table 5. Detailed performance metrics for all models in validation cohorts.

To evaluate the generalization capability of the top-performing diagnostic models, we conducted cross-validation with 10-fold resampling to prevent overfitting. The ROC curve results (Figure 4) confirm that our constructed model demonstrates excellent generalization performance without signs of overfitting or underfitting. Based on these findings, the XGBoost model was identified as the most effective model in this study.

Figure 4
ROC Curve Plot showing sensitivity versus 1-specificity for ten XGB models. The curves are close to each other, with various area under the curve (AUC) values ranging from 0.863 to 0.884. Key models and their AUCs are listed: XGB_BR10TEST, 0.881; XGB_BR1TEST, 0.874; XGB_BR2TEST, 0.873; and XGB_BR3TEST, 0.884. A diagonal reference line indicates random performance.

Figure 4. Ten fold cross validation evaluation. ROC curve plot of XGBoost.

3.4 Clinical decision threshold analysis

ROC curve analysis with the AUC-ROC index identified a statistically optimal threshold of 0.489 (95% CI: 0.483 - 0.509), achieving 87.1% sensitivity and 81.8% specificity at this level (Supplementary Table 3). The decision curve demonstrated a net benefit of 0.3614 at the 35% threshold, confirming its clinical applicability. We established a three-tier risk stratification system: Low-risk (<20%): Standard calcium supplementation regimen; Moderate-risk (20 - 35%): Close monitoring; High-risk (≥35%): Intensive calcium supplementation regimen. The XGBoost prediction model developed in this study reduced HBS undiagnosed rate from 19.4% to 0% at the 35% decision threshold, providing a robust tool for developing personalized calcium supplementation strategies in clinical practice.

3.5 Interpretability analysis

The SHAP interpretability method system reveals the decision-making mechanism of XGBoost model in predicting HBS risk. Figure 5A (feature importance bar chart) shows that preoperative PTH level (mean SHAP = 0.32) is the most significant predictive factor, followed by Age,%PTH, Ca, and ALP. Figure 5B (bee swarm plot) further clarifies feature directionality: Preoperative high calcium levels significantly reduce HBS risk (SHAP value negative skew), while younger patients and elevated ALP drive risk increase (SHAP value positive clustering). The individualized decision mechanism is illustrated in Figure 5C (force plot): Taking a 47-year-old patient (preoperative PTH 1974 pg/mL, %PTH 84, calcium 2.35 mmol/L, ALP 292 U/L) as an example, the model predicts a 98.0% probability of HBS.

Figure 5
Panel A shows a bar chart of mean SHAP values for an XGBoost model, indicating the importance of features: Pre_PTH, Age, %PTH, Ca, and ALP. Panel B displays a scatter plot of SHAP values colored by feature value, illustrating the impact of variables on predictions. Panel C presents a force plot, explaining prediction contributions with features ALP, Pre_PTH, Ca, Age, and %PTH, highlighting a base value of 0.52 and a final prediction of 1.

Figure 5. Interpretability of machine learning predictions using SHAP. (A) SHAP importance plot; (B) SHAP beeswarm plot; (C) SHAP force plot. Pre-PTH, preoperative parathyroid hormone; %PTH, [(Preoperative PTH – PTH at skin closure)/Preoperative PTH] × 100%; ALP, Alkaline phosphatase (U/L); Ca, Serum calcium (mmol/L).

3.6 Application of model

To bridge the gap between predictive analytics and clinical utility, we operationalized the XGBoost model through an intuitive web-based platform (https://chaiyalin.shinyapps.io/make_web/). This tool enables clinicians to input patient-specific parameters—including preoperative PTH, ALP, serum calcium, age, %PTH—to generate real-time HBS risk stratification. The platform’s device-agnostic architecture ensures seamless access across desktop and mobile environments without requiring specialized software installation.

4 Discussion

Postoperative hungry bone syndrome (HBS) is a frequent complication following parathyroidectomy (PTX) in secondary hyperparathyroidism (SHPT) patients. We found that preoperative PTH was the most important risk factor for postoperative HBS. The possible reason is the rapid reduction in PTH levels shifts the bone turnover equilibrium toward mineralization, resulting in massive calcium influx into osteoid tissue and subsequent hypocalcemia (19, 20). This study identifies %PTH as independent predictors of HBS, leveraging PTH’s short half-life (3 – 5 minutes) to reflect real-time parathyroid functional dynamics.

Alkaline phosphatase (ALP), a biomarker of osteoblast activity, emerged as a critical predictor of HBS, consistent with prior studies (2729). Elevated preoperative ALP levels signify heightened bone turnover, which exacerbates postoperative calcium sequestration (30). Notably, ALP demonstrates stronger correlations with bone histomorphometric changes and dialysis patient mortality than PTH (31), underscoring its clinical utility in risk stratification. While bone-specific ALP (B-ALP) offers superior specificity for bone metabolism (21), its exclusion from routine clinical assays in primary care settings limits its practical applicability. Our study corroborates prior evidences (23, 32) demonstrating that preoperative hypocalcemia independently predicts postoperative HBS, likely reflecting an elevated baseline bone remodeling state in SHPT patients (33). Mechanistically, abrupt postoperative PTH suppression disrupts the bone resorption-formation equilibrium, triggering accelerated osteoblast-driven mineralization and subsequent calcium efflux from circulation to bone tissue. Notably, younger age emerged as an independent HBS risk factor, aligning with findings by Kritmetapak et al. (34) and He et al. (23). This association may stem from heightened osteoblast activity and skeletal calcium utilization efficiency in younger individuals. Conversely, older patients exhibited reduced susceptibility, potentially attributable to age-related declines in 1α-hydroxylase activity, chronic malnutrition, and impaired osteoblast calcium uptake capacity. However, conflicting data from Gong et al. (35) associate advanced age with HBS risk, possibly due to prevalent vitamin D deficiency and protein-energy wasting in elderly populations. These discrepancies likely arise from cohort heterogeneity and surgical technique variations, underscoring the need for multicenter studies to clarify age-specific risk profiles. Additionally, we observed a reduced HBS rate among patients using bisphosphonates preoperatively, which aligns with other studies. Some authors have documented that bisphosphonates can improve HBS (19, 36, 37). Davenport et al. (19)described a study showing that pamidronate could decrease postoperative HBS rates.

In recent years, machine learning has shown great potential in disease diagnosis and prognosis. However, few studies have been conducted to predict postoperative HBS using ML models. This study introduces an interpretable machine learning model that individualizes HBS risk prediction using five clinically accessible variables. The model’s operationalization as a web-based calculator enables real-time risk stratification, aligning with KDIGO guidelines for proactive postoperative management. Internal validation via cohort partitioning demonstrated robust discriminatory performance, supporting its reliability for clinical deployment. Moreover, traditional evaluation methods often exhibit high rates of missed diagnoses. Our XGBoost model effectively reduces diagnostic gaps by integrating multidimensional features such as%PTH levels while maintaining appropriate specificity. We have further developed a risk stratification framework (low, medium, and high-risk) based on this threshold, establishing corresponding preventive measures and monitoring strategies. This model is designed for intraoperative decision support, enabling surgeons to adjust calcium prophylaxis strategies based on real-time %PTH dynamics. High-risk patients trigger immediate IV calcium infusion protocols.

However, several limitations warrant consideration. This study has the following limitations: 1. The single-center retrospective design limits model generalizability, requiring multicenter prospective validation; 2. The HBS definition is based solely on biochemical indicators without incorporating symptom assessment, which may lead to minor misclassification; 3. Key bone metabolism markers such as bone-specific alkaline phosphatase and osteocalcin are absent (38); 4. Model calibration performance may drift over time, necessitating periodic recalibration; 5. Differences in surgical techniques across centers and variations in surgeons’ expertise may affect model applicability.

5 Conclusion

Our findings demonstrate that %PTH exhibits a strong independent association with postoperative hungry bone syndrome (HBS). Leveraging machine learning algorithms, we developed and validated a predictive model incorporating five clinically accessible preoperative parameters: 1) preoperative PTH, 2) %PTH, 3) age, 4) serum calcium, and 5) alkaline phosphatase (ALP). This model serves as a practical clinical tool for early identification of high-risk patients, enabling targeted prophylactic interventions such as preoperative calcium optimization and intensified postoperative monitoring. By utilizing routinely available preoperative biomarkers, our approach bridges the gap between predictive analytics and actionable clinical decision-making in CKD-associated SHPT management.

Data availability statement

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

Ethics statement

The studies involving humans were approved by the Ethics Committee of the Affiliated Hospital of Qingdao University (Approval No. QYFY WZLL 29980). 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 due to the use of anonymized retrospective data and the absence of additional risks.

Author contributions

YC: Formal analysis, Visualization, Data curation, Validation, Project administration, Methodology, Writing – review & editing, Conceptualization, Writing – original draft. NY: Writing – review & editing, Methodology, Project administration, Visualization, Data curation. JY: Writing – review & editing, Formal analysis, Data curation, Methodology. BS: Data curation, Writing – review & editing, Formal analysis, Methodology. LS: Supervision, Writing – review & editing, Validation, Software. LZ: Writing – review & editing, Project administration, Data curation, Formal analysis. LY: Formal analysis, Writing – review & editing, Data curation. XW: Formal analysis, Data curation, Writing – review & editing. FL: Supervision, Conceptualization, Writing – review & editing, Methodology, Funding acquisition, Writing – original draft, Formal analysis. CL: Writing – original draft, Funding acquisition, Writing – review & editing, Formal analysis, Methodology, Supervision, Conceptualization.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by the National Natural Science Foundation of China (grant numbers 81700585); the Qingdao Key Health Discipline Development Fund; and the Qingdao Key Clinical Specialty Elite Discipline.

Acknowledgments

We thank the participants in the study.

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.

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.1635451/full#supplementary-material

Abbreviations

AdaBoost, Adaptive Boosting; ALP, Alkaline Phosphatase; ROC, Receiver operating characteristic curve; CatBoost, Categorical Boosting; CKD, Chronic kidney disease; DCA, Decision curve analysis; HBS, Hungry bone syndrome; KNN:K-nearest neighbors; LightGBM, Light gradient boosting machine; LR, Logistic regression; NN, Neural Network; PTH, Parathyroid hormone; PTX, Parathyroidectomy; Scr, Serum creatinine; SHAP, Shapley additive explanations; SHPT, Secondary hyperparathyroidism; SVM, Support vector machine; XGBoost, Extreme gradient boosting.

References

1. Kalantar-Zadeh K, Kuwae N, Regidor DL, Kovesdy CP, Kilpatrick RD, Shinaberger CS, et al. Survival predictability of time-varying indicators of bone disease in maintenance hemodialysis patients. Kidney Int. (2006) 70:771–80. doi: 10.1038/sj.ki.5001514

PubMed Abstract | Crossref Full Text | Google Scholar

2. Fligor SC, Li C, Hamaguchi R, William J, and James BC. Decreasing surgical management of secondary hyperparathyroidism in the United States. J Surg Res. (2021) 264:444–53. doi: 10.1016/j.jss.2021.03.013

PubMed Abstract | Crossref Full Text | Google Scholar

3. Ding Y, Wang H, Zou Q, Jin Y, Zhang Z, and Huang J. Factors associated with calcium requirements after parathyroidectomy in chronic kidney disease patients. Int Urol nephrol. (2018) 50:535–40. doi: 10.1007/s11255-017-1759-8

PubMed Abstract | Crossref Full Text | Google Scholar

4. Xu Y, Evans M, Soro M, Barany P, and Carrero JJ. Secondary hyperparathyroidism and adverse health outcomes in adults with chronic kidney disease. Clin Kidney J. (2021) 14:2213–20. doi: 10.1093/ckj/sfab006

PubMed Abstract | Crossref Full Text | Google Scholar

5. Thadhani RI, Rosen S, Ofsthun NJ, Usvyat LA, Dalrymple LS, Maddux FW, et al. Conversion from intravenous vitamin D analogs to oral calcitriol in patients receiving maintenance hemodialysis. Clin J Am Soc Nephrol. (2020) 15:384–91. doi: 10.2215/cjn.07960719

PubMed Abstract | Crossref Full Text | Google Scholar

6. Alexander ST, Hunter T, Walter S, Dong J, Maclean D, Baruch A, et al. Critical cysteine residues in both the calcium-sensing receptor and the allosteric activator AMG 416 underlie the mechanism of action. Mol Pharmacol. (2015) 88:853–65. doi: 10.1124/mol.115.098392

PubMed Abstract | Crossref Full Text | Google Scholar

7. Hiramitsu T, Hasegawa Y, Futamura K, Okada M, Goto N, Narumi S, et al. Treatment for secondary hyperparathyroidism focusing on parathyroidectomy. Front endocrinol. (2023) 14:1169793. doi: 10.3389/fendo.2023.1169793

PubMed Abstract | Crossref Full Text | Google Scholar

8. Ketteler M, Block GA, Evenepoel P, Fukagawa M, Herzog CA, McCann L, et al. Diagnosis, evaluation, prevention, and treatment of chronic kidney disease-mineral and bone disorder: synopsis of the kidney disease: improving global outcomes 2017 clinical practice guideline update. Ann Internal Med. (2018) 168:422–30. doi: 10.7326/m17-2640

PubMed Abstract | Crossref Full Text | Google Scholar

9. Lau WL, Obi Y, and Kalantar-Zadeh K. Parathyroidectomy in the management of secondary hyperparathyroidism. Clin J Am Soc Nephrol: CJASN. (2018) 13:952–61. doi: 10.2215/cjn.10390917

PubMed Abstract | Crossref Full Text | Google Scholar

10. Komaba H, Nakamura M, and Fukagawa M. Resurgence of parathyroidectomy: evidence and outcomes. Curr Opin Nephrol hypertension. (2017) 26:243–9. doi: 10.1097/mnh.0000000000000326

PubMed Abstract | Crossref Full Text | Google Scholar

11. Cheng SP, Lee JJ, Liu TP, Yang T-L, Chen H-H, Wu C-J, et al. Parathyroidectomy improves symptomatology and quality of life in patients with secondary hyperparathyroidism. Surgery. (2014) 155:320–8. doi: 10.1016/j.surg.2013.08.013

PubMed Abstract | Crossref Full Text | Google Scholar

12. Filho WA, van der Plas WY, Brescia MDG, Nascimento CP, Goldenstein PT, Neto LMM, et al. Quality of life after surgery in secondary hyperparathyroidism, comparing subtotal parathyroidectomy with total parathyroidectomy with immediate parathyroid autograft: Prospective randomized trial. Surgery. (2018) 164:978–85. doi: 10.1016/j.surg.2018.06.032

PubMed Abstract | Crossref Full Text | Google Scholar

13. van der Plas W, Kruijff S, Sidhu SB, Delbridge LW, Sywak MS, and Engelsman AF. Parathyroidectomy for patients with secondary hyperparathyroidism in a changing landscape for the management of end-stage renal disease. Surgery. (2021) 169:275–81. doi: 10.1016/j.surg.2020.08.014

PubMed Abstract | Crossref Full Text | Google Scholar

14. Carsote M and Nistor C. Forestalling hungry bone syndrome after parathyroidectomy in patients with primary and renal hyperparathyroidism. Diagnostics (Basel Switzerland). (2023) 13. doi: 10.3390/diagnostics13111953

PubMed Abstract | Crossref Full Text | Google Scholar

15. Tai YL, Shen HY, Nai WH, Fu J-F, Wang I-K, Huang C-C, et al. Hungry bone syndrome after parathyroid surgery. Hemodialysis Int Int Symposium Home Hemodialysis. (2023) 27:134–45. doi: 10.1111/hdi.13067

PubMed Abstract | Crossref Full Text | Google Scholar

16. Amjad W, Ginzberg SP, Passman JE, Heintz J, Kelz RR, and Wachtel H. Predictive risk score for postparathyroidectomy hungry bone syndrome in patients with secondary hyperparathyroidism. J Clin Endocrinol Metab. (2024) 109:603–10. doi: 10.1210/clinem/dgad636

PubMed Abstract | Crossref Full Text | Google Scholar

17. Witteveen JE, van Thiel S, Romijn JA, and Hamdy NA. Hungry bone syndrome: still a challenge in the post-operative management of primary hyperparathyroidism: a systematic review of the literature. Eur J endocrinol. (2013) 168:R45–53. doi: 10.1530/eje-12-0528

PubMed Abstract | Crossref Full Text | Google Scholar

18. Jain N and Reilly RF. Hungry bone syndrome. Curr Opin Nephrol hypertension. (2017) 26:250–5. doi: 10.1097/mnh.0000000000000327

PubMed Abstract | Crossref Full Text | Google Scholar

19. Schnyder MA, Stolzmann P, Huber GF, and Schmid C. A patient with a history of breast cancer and multiple bone lesions: a case report. J Med Case Rep. (2017) 11:127. doi: 10.1186/s13256-017-1296-1

PubMed Abstract | Crossref Full Text | Google Scholar

20. Yajima A, Ogawa Y, Takahashi HE, Tominaga Y, Inou T, and Otsubo O. Changes of bone remodeling immediately after parathyroidectomy for secondary hyperparathyroidism. Am J Kidney Dis. (2003) 42:729–38. doi: 10.1016/s0272-6386(03)00909-0

PubMed Abstract | Crossref Full Text | Google Scholar

21. Wang M, Chen B, Zou X, Wei T, Gong R, Zhu J, et al. A nomogram to predict hungry bone syndrome after parathyroidectomy in patients with secondary hyperparathyroidism. J Surg Res. (2020) 255:33–41. doi: 10.1016/j.jss.2020.05.036

PubMed Abstract | Crossref Full Text | Google Scholar

22. Cao R, Jiang H, Liang G, and Zhang W. Dynamic nomogram for predicting hungry bone syndrome before parathyroidectomy. Endocrine. (2024) 83:196–204. doi: 10.1007/s12020-023-03493-6

PubMed Abstract | Crossref Full Text | Google Scholar

23. He C, Zhang Y, Li L, Cheng G, Zhang W, Tang Y, et al. Risk Factor Analysis and Prediction of Severe Hypocalcemia after Total Parathyroidectomy without Auto-Transplantation in Patients with Secondary Hyperparathyroidism. Int J Endocrinol. (2023) 2023:1901697. doi: 10.1155/2023/1901697

PubMed Abstract | Crossref Full Text | Google Scholar

24. Muller O, Bauvin P, Bacoeur O, Michailos T, Bertoni M-V, Demory C, et al. Machine learning-based algorithm for the early prediction of postoperative hypocalcemia risk after thyroidectomy. Ann surgery. (2024) 280:835–41. doi: 10.1097/sla.0000000000006480

PubMed Abstract | Crossref Full Text | Google Scholar

25. Weissler EH, Naumann T, Andersson T, Ranganath R, Elemento O, Luo Y, et al. The role of machine learning in clinical research: transforming the future of evidence generation. Trials. (2021) 22:537. doi: 10.1186/s13063-021-05489-x

PubMed Abstract | Crossref Full Text | Google Scholar

26. An Q, Rahman S, and Zhou J. Kang JJ. A comprehensive review on machine learning in healthcare industry: classification, restrictions, opportunities and challenges. Sensors (Basel Switzerland). (2023) 23. doi: 10.3390/s23094178

PubMed Abstract | Crossref Full Text | Google Scholar

27. Ge Y, Yang G, Wang N, Zha X, Yu X, Mao H, et al. Bone metabolism markers and hungry bone syndrome after parathyroidectomy in dialysis patients with secondary hyperparathyroidism. Int Urol nephrol. (2019) 51:1443–9. doi: 10.1007/s11255-019-02217-y

PubMed Abstract | Crossref Full Text | Google Scholar

28. Gao D, Liu Y, Cui W, Lu X, and Lou Y. A nomogram prediction model for hungry bone syndrome in dialysis patients with secondary hyperparathyroidism after total parathyroidectomy. Eur J Med Res. (2024) 29:208. doi: 10.1186/s40001-024-01801-y

PubMed Abstract | Crossref Full Text | Google Scholar

29. Yang G, Zha X, Mao H, Yu X, Wang N, and Xing C. Hypocalcemia-based prediction of hungry bone syndrome after parathyroidectomy in hemodialysis patients with refractory secondary hyperparathyroidism. J Int Med Res. (2018) 46:4985–94. doi: 10.1177/0300060518788744

PubMed Abstract | Crossref Full Text | Google Scholar

30. Vimalraj S. Alkaline phosphatase: Structure, expression and its function in bone mineralization. Gene. (2020) 754:144855. doi: 10.1016/j.gene.2020.144855

PubMed Abstract | Crossref Full Text | Google Scholar

31. Lau WL, Kalantar-Zadeh K, Kovesdy CP, and Mehrotra R. Alkaline phosphatase: better than PTH as a marker of cardiovascular and bone disease? Hemodialysis Int Int Symposium Home Hemodialysis. (2014) 18:720–4. doi: 10.1111/hdi.12190

PubMed Abstract | Crossref Full Text | Google Scholar

32. Cheng J, Lv Y, Zhang L, and Liu Y. Construction and validation of a predictive model for hypocalcemia after parathyroidectomy in patients with secondary hyperparathyroidism. Front Endocrinol (Lausanne). (2022) 13:1040264. doi: 10.3389/fendo.2022.1040264

PubMed Abstract | Crossref Full Text | Google Scholar

33. Habas EC, Eledrisi M, Khan F, and Elzouki AY. Secondary hyperparathyroidism in chronic kidney disease: pathophysiology and management. Cureus. (2021) 13:e16388. doi: 10.7759/cureus.16388

PubMed Abstract | Crossref Full Text | Google Scholar

34. Kritmetapak K, Kongpetch S, Chotmongkol W, Raruenrom Y, Sangkhamanon S, and Pongchaiyakul C. Incidence of and risk factors for post-parathyroidectomy hungry bone syndrome in patients with secondary hyperparathyroidism. Ren Fail. (2020) 42:1118–26. doi: 10.1080/0886022x.2020.1841655

PubMed Abstract | Crossref Full Text | Google Scholar

35. Gong W, Lin Y, Xie Y, Meng Z, and Wang Y. Predictors of early postoperative hypocalcemia in patients with secondary hyperparathyroidism undergoing total parathyroidectomy. J Int Med Res. (2021) 49:3000605211015018. doi: 10.1177/03000605211015018

PubMed Abstract | Crossref Full Text | Google Scholar

36. Fackelmayer OJ, Wu JX, and Yeh MW. Endocrine surgery: management of postoperative complications following endocrine surgery of the neck. Surg Clinics North America. (2021) 101:767–84. doi: 10.1016/j.suc.2021.05.019

PubMed Abstract | Crossref Full Text | Google Scholar

37. Kaderli RM, Riss P, Dunkler D, Pietschmann P, Selberherr A, Scheuba C, et al. The impact of vitamin D status on hungry bone syndrome after surgery for primary hyperparathyroidism. Eur J endocrinol. (2018) 178:1–9. doi: 10.1530/eje-17-0416

PubMed Abstract | Crossref Full Text | Google Scholar

38. Ko WC, Liu CL, Lee JJ, Liu TP, Wu CJ, and Cheng SP. Osteocalcin is an independent predictor for hungry bone syndrome after parathyroidectomy. World J surgery. (2020) 44:795–802. doi: 10.1007/s00268-019-05251-0

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: hungry bone syndrome, secondary hyperparathyroidism, parathyroidectomy, risk factors, machine learning

Citation: Chai Y, Yuan N, Yin J, Shen B, Sun L, Zhang L, Yin L, Wang X, Luo F and Luo C (2025) Machine learning-based predictive model for hungry bone syndrome following parathyroidectomy in secondary hyperparathyroidism. Front. Endocrinol. 16:1635451. doi: 10.3389/fendo.2025.1635451

Received: 27 May 2025; Accepted: 18 August 2025;
Published: 05 September 2025.

Edited by:

Sree Bhushan Raju, Nizam’s Institute of Medical Sciences, India

Reviewed by:

Anna Eremkina, Endocrinology Research Center, Russia
Shuai Lu, Beijing Jishuitan Hospital, China
Petru Adrian Radu, Nephrology Clinical Hospital “Dr. Carol Davila”, Romania

Copyright © 2025 Chai, Yuan, Yin, Shen, Sun, Zhang, Yin, Wang, Luo and Luo. 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: Feng Luo, bHVvZmVuZzE5ODIxMTI0QDEyNi5jb20=; Congjuan Luo, bHVvY29uZ2p1YW4yMDE4QDE2My5jb20=

ORCID: Congjuan Luo, orcid.org/0000-0003-1952-8668

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.