- 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 (1–3). 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 (10–13).
Postoperative hungry bone syndrome (HBS), a complication affecting 25%-75% of PTX cases (14–16), 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 (21–23). 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.
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).
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. 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.
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. 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.
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.
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. 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 (27–29). 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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
18. Jain N and Reilly RF. Hungry bone syndrome. Curr Opin Nephrol hypertension. (2017) 26:250–5. doi: 10.1097/mnh.0000000000000327
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
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
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
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
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
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
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
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
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
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
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
30. Vimalraj S. Alkaline phosphatase: Structure, expression and its function in bone mineralization. Gene. (2020) 754:144855. doi: 10.1016/j.gene.2020.144855
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
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
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
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
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
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
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
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, IndiaReviewed by:
Anna Eremkina, Endocrinology Research Center, RussiaShuai 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
Yalin Chai