- 1Department of Orthopaedics, The First Affiliated Hospital of Soochow University, Suzhou, China
- 2Department of Otorhinolaryngology, The First Affiliated Hospital of Soochow University, Suzhou, China
- 3Department of Orthopaedics, The Second Affiliated Hospital of Soochow University, Suzhou, China
Background: Osteoporosis in postmenopausal women is characterized by significant bone mass loss due to reduced oestrogen, leading to an increased risk of osteoporotic vertebral compression fractures (OVCF). Comprehensive risk prediction models for diagnosing and predicting fracture risk in this population are still lacking.
Objective: This study aims to identify key risk factors for OVCF in postmenopausal osteoporotic women and develop a machine learning model to predict OVCF risk by integrating clinical, biological, and musculoskeletal data.
Methods: This retrospective case-control study included 486 postmenopausal women diagnosed with osteoporosis between 2015 and 2018. The patients were divided into a non-fracture group (Group A) and a vertebral fracture group (Group B) based on whether they developed OVCF during the subsequent 5 years of treatment and follow-up. Univariate and multivariate logistic regression analyses were performed to identify independent risk factors for OVCF. Furthermore, a comprehensive risk prediction model was constructed using multiple machine learning algorithms.
Results: Among the 486 postmenopausal women, 269 (55.35%) experienced OVCF. Low bone mineral density (BMD), chronic inflammation, and sarcopenia were identified as independent risk factors, while regular anti-osteoporotic treatment was associated with a reduced fracture incidence. The Balanced Bagging machine learning model demonstrated an accuracy of 98.98%, a sensitivity of 98.24%, a specificity of 100.00%, and the model’s F1-score was 0.99. The deployed model outputs calibrated, patient-specific probabilities with case-level explanations and supports dynamic re-scoring as new BMD/CTx/NLR results become available, enabling personalized risk management in routine care.
Conclusion: The development of OVCF in postmenopausal osteoporotic women is influenced by a combination of bone metabolism, inflammatory processes, and muscle health. The machine learning model developed in this study provides a reliable and accurate tool for personalized OVCF risk prediction, allowing clinicians to optimize prevention and treatment strategies. Future large-scale prospective studies are required to validate these findings and enhance the model’s predictive capabilities.
Highlights
• Low bone mineral density (BMD), chronic inflammation (elevated NLR), and sarcopenia are identified as key independent risk factors for OVCF, while regular anti-osteoporotic treatment was associated with a reduced fracture incidence.
• The Balanced Bagging machine learning model integrates clinical, biological, and musculoskeletal data to predict osteoporotic vertebral compression fracture (OVCF) risk in postmenopausal women with 98.98% accuracy, 98.24% sensitivity, 100.00% specificity and the model’s F1-score was 0.99.
• Comprehensive multidimensional anti-osteoporosis therapy is expected to significantly reduce the incidence of OVCF, emphasizing the importance of comprehensive and multidimensional treatment of high-risk populations.
1 Introduction
Osteoporosis is characterized by reduced bone mass and disruption of bone microarchitecture. In postmenopausal women, oestrogen decline accelerates bone resorption, lowering bone mineral density (BMD) and increasing fragility fracture risk (1, 2). Among them, OVCF is one of the most frequent osteoporotic fractures (3). These fractures often occur in the thoracolumbar spine, leading to significant somatic pain and activity limitation, increasing the risk of disability and affecting patients’ quality of life (4, 5). Therefore, identifying the risk causes for OVCF in postmenopausal osteoporotic women and reducing the probability of fracture through effective preventive and therapeutic means are of great humanistic and socioeconomic importance. Despite numerous studies examining the relationship between osteoporosis and fragility fractures (4–8), the factors are not comprehensive enough to accurately assess the way they interact. Although BMD is widely used to diagnose osteoporosis and estimate fracture risk, reliance on BMD alone is insufficient because risk is also influenced by bone quality, systemic inflammation and bone turnover (9, 10). Moreover, the specific determinants of OVCF in postmenopausal women remain incompletely defined, and few studies have integrated diverse dimensions of risk in this population (11–14).
Bone turnover markers, including type I procollagen N-terminal propeptide (P1NP; formation) and carboxy-terminal cross-linked telopeptide of type I collagen (CTx; resorption), reflect skeletal metabolism, while osteocalcin (BGP) indicates osteoblastic activity (14–17). Calcium (Ca), phosphorus (P) and vitamin D3 (Vit D3) contribute to bone health maintenance, and clinical factors such as nutrition, hypertension and diabetes mellitus may modify OVCF risk (18, 19). Low-grade inflammation, captured by the neutrophil-to-lymphocyte ratio (NLR) and the systemic immune-inflammation index (SII), may potentiate bone resorption and fracture susceptibility (20). Emerging evidence links sarcopenia to higher vertebral fracture risk and poorer recovery. Postmenopausal women show greater age- and sex-related losses in muscle mass and strength, with more fatty infiltration of paraspinal muscles, supporting a role for muscle status in OVCF susceptibility. We therefore assessed the psoas major and multifidus muscle indices (PMI/MMI) as imaging surrogates (21, 22). Future work incorporating muscle quality and functional measures may clarify the incremental contribution of sarcopenia in postmenopausal osteoporosis. Against this background, we developed a machine-learning model that integrates clinical variables, densitometry, bone-turnover markers, inflammatory indices, and muscle metrics to identify the risk factors for OVCF and explore the correlations among these factors, thereby enabling personalized risk prediction and informing targeted prevention and treatment strategies (23, 24).
2 Methods
2.1 Management of the sample
The study protocol was approved by the Ethics Committee of the First Affiliated Hospital of Soochow University (2024 Lun Research Grant No. 728). This study investigated factors associated with osteoporotic vertebral compression fractures (OVCF) in postmenopausal women with osteoporosis. We analyzed 486 postmenopausal women diagnosed with primary osteoporosis between January 2015 and December 2018, all of whom received up to 5 years of outpatient follow-up. Patients were categorized into two groups based on the occurrence of OVCF. Inclusion criteria were: postmenopausal women with a menopause duration of 1–20 years, a diagnosis of primary osteoporosis, and regular follow-up for at least 5 years. Exclusion criteria included severe congenital or acquired spinal deformities, conditions significantly affecting bone metabolism, and incomplete medical or imaging records. Among the diseases that significantly affect bone metabolism are (1) chronic kidney disease (eGFR <60 mL−1 min−1). (2) Hyper- or hypoparathyroidism and thyroid dysfunction. (3) Chronic liver disease (Child–Pugh B/C). (4) Malabsorption syndrome. (5) Active malignancy or recent chemotherapy/radiotherapy. (6) Long-term systemic glucocorticoid therapy (≥5 mg prednisolone >3 months), aromatase inhibitors, antiepileptic drugs, heparin/warfarin. (7) Severe rheumatoid arthritis or other inflammatory arthritis. (8) Chronic alcohol abuse.
2.2 Observation indicators
We retrospectively collected demographic data, laboratory results, and imaging information, including age, BMI, fracture history, histories of hypertension and diabetes, serum calcium, vitamin D3, osteocalcin (BGP), type I collagen amino-terminal propeptide (P1NP), carboxy-terminal cross-linked telopeptide (CTx), bone mineral density (BMD), neutrophil-to-lymphocyte ratio (NLR), psoas major index (PMI), and multifidus muscle index (MMI). All patients underwent serological tests and dual-energy X-ray absorptiometry (DEXA) at diagnosis and during follow-up to assess BMD and bone metabolism markers. Time 0 was defined as the date of the index DXA. A standardised baseline assessment pathway was implemented, whereby bone turnover markers and lumbar MRI (for PMI/MMI) were obtained routinely on the same day as DXA whenever feasible, and otherwise within four weeks. Osteoporosis was defined according to the World Health Organization (WHO) criteria using dual-energy X-ray absorptiometry (DXA) at the lumbar spine (L1–L4) or the total hip. A T-score of ≤−2.5 indicates osteoporosis. In addition, a T-score between −1.0 and −2.5, in conjunction with fragility fractures of the pelvis and lumbar spine, can also serve as diagnostic criteria for osteoporosis. Given that this study focuses on vertebral fragility fractures caused by decreased bone density, and based on prior research and clinical experience, we set a T-score of <−1.5 as the criterion for sample inclusion.
Anti-osteoporosis medications at baseline and during follow-up (bisphosphonates, denosumab, teriparatide, or SERMs) were extracted from the electronic record and summarised in Table 1. Regimens followed product labelling (denosumab 60 mg SC every six months; teriparatide 20 μg SC daily) with vitamin D and calcium supplementation. BMD was reassessed at least annually (six-monthly in selected cases). Incident OVCF events were identified from the electronic record and confirmed radiographically using predefined criteria. The case records documented that participants underwent scheduled annual lateral thoracolumbar radiographs and symptom-triggered imaging to identify incident vertebral fractures. An incident vertebral fracture was defined according to Genant semi-quantitative criteria, or as a ≥20% and ≥4 mm reduction in vertebral body height (anterior, middle, or posterior) on follow-up compared with the most recent negative radiograph. Cases involving high-energy trauma were excluded, and each patient underwent lumbar magnetic resonance imaging (MRI) within 4 weeks of a radiograph showing a lumbar fracture to confirm the diagnosis.

Table 1. General data comparison between No vertebral fracture group (Group A) and Osteoporotic vertebral compression fracture group (Group B) patients.
Additionally, cross-sectional areas of the psoas and multifidus muscles at the L3 level were measured on T2-weighted MRI, with PMI and MMI calculated by normalizing muscle area to height squared (cm2/m2). Sarcopenia was defined as a PMI below 3.9 cm2/m2. All images were independently reviewed by two fellowship-trained musculoskeletal radiologists who were blinded to the clinical predictors, with disagreements resolved by consensus.
2.3 Statistical analysis
Statistical analysis was performed using SPSS software (version 27.0, IBM, United States), presenting continuous data as mean ± standard deviation analyzed with independent t-tests. Categorical data were evaluated using chi-square tests or Fisher’s exact tests. A p-value <0.05 was considered significant. Variables that demonstrated significant differences in the univariate analysis were subsequently included in a multivariate logistic regression to identify independent risk factors for OVCF among postmenopausal women with osteoporosis. Prior to multivariable modelling, we assessed pairwise dependencies and applied a VIF threshold of <5.0 to evaluate collinearity. For clusters of biologically related predictors (e.g., SII with NLR; CTx with P1NP/BGP), we retained the clinically more parsimonious and statistically more robust variable (NLR and CTx, respectively). For transparency, BGP, P1NP, SII, and PMI were included in descriptive statistics and exploratory visualisations but were not carried forward into the comparative modelling process.
2.4 Model analysis
In this study, various machine learning models—including Logistic Regression, Naive Bayes, SVM, Decision Tree, AdaBoost, Gradient Boosting, and Balanced Bagging—were developed using Python 3.11.7 with the scikit-learn library. Since the retrospectively collected laboratory dataset contained only a small fraction of missing values, we employed the missForest algorithm, thereby minimizing information loss while preserving the underlying distribution structure. Model training employed 10-fold cross-validation and grid search for hyperparameter optimization: the former partitions the training data into 10 subsets, iteratively using nine for training and one for validation to reduce partition bias, while grid search systematically explores a predefined parameter space to identify the optimal parameter combination. Model performance was evaluated using metrics such as AUC, sensitivity, specificity, PPV, NPV, accuracy, F1-score, and Brier Score, ensuring a comprehensive assessment of classification ability, stability, and calibration. To control overfitting in high-variance learners, we imposed a priori constraints (for tree/boosting models: max_depth ≤3, min_samples_leaf ≥20, learning_rate ≤0.05) and used class-balanced resampling within a Balanced Bagging framework. Learning curve analyses showed that performance plateaued once the training sample exceeded approximately 200 cases; permutation tests (1,000 shuffles) yielded an AUC of approximately 0.50, indicating that the observed discrimination was not due to chance. Additionally, SHAP (SHapley Additive exPlanations) was applied to quantify each feature’s contribution to the model’s predictions, enhancing interpretability and providing insights into the decision-making process. The same individual can be re-scored at follow-up when new densitometry or laboratory results are available, allowing trajectory monitoring after treatment initiation or adjustment. A five-feature configuration (BMD T-score, CTx, NLR, age, BMI) preserved >90% of full-model discrimination, supporting deployment where selected assays are unavailable.
3 Result
3.1 Patient characteristics
A total of 486 postmenopausal osteoporotic women (mean age 63.18 ± 4.57 years) were included, with 217 in the non-fracture group (Group A) and 269 in the OVCF group (Group B). Significant differences between groups were observed in BMI (p = 0.045), fracture history (p = 0.001), BMD (p = 0.001), calcium (p = 0.027), Vit D3 (p = 0.001), and bone turnover markers (CTx and P1NP, p = 0.001) (Table 1).
3.2 Visualization analysis
Box plots showed that Group B had significantly higher age and CTx levels, and lower BMD and Vit D3 levels (Figure 1). Bar charts revealed that Group B had a higher fracture history and lower osteoporosis treatment rates compared to Group A (Figure 2), suggesting that age, BMD, fracture history, and treatment status are key predictors of OVCF.

Figure 1. Box plot comparison of continuous variables by groups. Using box plots and scatter points, it displays the median, quartiles, and outliers of each variable in different groups, providing a visual comparison of the differences between these variables.

Figure 2. Categorical variable distribution between Group A and Group B. The graph shows the distribution of several categorical variables: treatment for osteoporosis, history of fracture, hypertension, and diabetes. The bar graph compares the counts for each category in these variables, providing a clear visual comparison of the frequency distribution between the two groups.
3.3 Univariate and multivariate analysis
Univariate analysis demonstrated significant associations (p < 0.05) between OVCF and variables such as age, BMI, fracture history, diabetes, BMD, osteoporosis treatment, calcium, Vit D3, CTx, P1NP, NLR, SII, and PMI. Multivariate logistic regression (Hosmer–Lemeshow p = 0.565) identified age, BMI, fracture history, BMD, osteoporosis treatment, calcium, Vit D3, CTx, and NLR as independent risk factors for OVCF (Table 2).

Table 2. Results of multi-factor logistic regression analysis of risk factors associated with osteoporotic vertebral compression fractures (OVCF) in postmenopausal women.
3.4 Multi-model predictive analysis
Internal cross-validation (Figures 3, 4) showed that model performance stabilized with over 200 training samples. On the external dataset (Figure 5), the Balanced Bagging model outperformed others, achieving an AUC of 0.9956 and an accuracy of 98.98%, with high sensitivity and specificity. SHAP analysis (Figures 6, 7) indicated that CTx, BMD, BMI, NLR, age, and Vit D3 were the most influential features—particularly CTx, BMD, and NLR—in discriminating between classes. Overall, the model effectively captured multidimensional patterns, supporting robust OVCF risk prediction.

Figure 6. SHAP analysis (negative)—feature importance for negative class. (a) SHAP summary dot plot. (b) SHAP waterfall plot for the 11th sample. (c) SHAP decision plot.

Figure 7. SHAP analysis (positive)—feature importance for positive class. (a) SHAP summary dot plot. (b) SHAP waterfall plot for the 11th sample. (c) SHAP decision plot.
4 Discussion
In this single-centre cohort of 486 postmenopausal women, we identified age, BMI, prior fracture, lower BMD, absence of anti-osteoporotic treatment, lower calcium, lower Vit D3, higher CTx and higher NLR as independent correlates of OVCF, reflecting the interplay of skeletal fragility, bone turnover and systemic inflammation (25). Our machine-learning analysis, using these nine features, achieved excellent discrimination and calibration; Balanced Bagging performed best (temporal hold-out AUC 0.9956; accuracy 98.98%), and learning curves plateaued beyond approximately 200 cases, supporting model stability. Brier score and calibration plots indicated well-calibrated probabilities. SHAP analyses consistently ranked CTx, BMD and NLR as the most influential contributors, with age, BMI and Vit D3 providing additional signal.
OVCF is the most prevalent complication of osteoporosis in postmenopausal women, typically affecting the thoracolumbar spine and causing pain, functional limitation and healthcare burden (26–29). Lower BMD remained a dominant risk factor in our data, consistent with its status as the diagnostic standard and a core determinant of vertebral fragility. Reduced mineral content is accompanied by deterioration of microarchitecture (trabecular thinning, loss of connectivity and endplate microcracks), decreasing load-sharing capacity and raising the likelihood of wedge or biconcave deformities under routine activities (30–33). CTx, a marker of type I collagen degradation, reflects systemic osteoclastic activity. Elevated CTx signifies uncoupled remodelling with increased resorption cavities, transient porosity and weaker bone packets before secondary mineralisation is complete (34, 35). In the vertebral body—rich in trabecular bone—this manifests as reduced stiffness and earlier yield under axial load. CTx outperformed formation markers (P1NP, BGP) for discrimination, implying that resorption acceleration, rather than formation dynamics, is the proximate driver of vertebral failure in this setting. SHAP showed monotonic risk escalation with rising CTx, supporting a dose-response relation. Chronic low-grade inflammation likely amplifies this risk: oestrogen decline is associated with increased pro-inflammatory cytokines (e.g., TNF-α, IL-6), which elevate NLR, suppress osteoblastic activity and enhance osteoclastic activity, thereby weakening bone structural strength (36–40). BMI showed a positive association with OVCF, a finding compatible with sarcopenic obesity and visceral adiposity, which promote systemic inflammation despite higher absolute mass; this may offset any mechanical protection conferred by weight alone (41). Although PMI was lower in the OVCF group at the univariate level, it did not remain independent after adjustment, suggesting that global anthropometry and biochemical drivers captured the incremental risk more parsimoniously in this dataset.
Among tested learners, the Balanced Bagging model yielded the lowest cross-validated deviance with favorable calibration, likely because bagging reduces variance under class imbalance and captures non-linear interactions without overfitting when depth and leaf size are constrained (30). The Balanced Bagging framework can be retrained with reduced feature sets, allowing deployment when certain biomarkers (e.g., CTx) are temporarily unavailable. SHAP analyses show that a core subset—BMD, CTx, NLR, age, BMI—retains >90% of the model’s discriminative capacity, supporting a “tiered” prediction approach adaptable to real-world data completeness. These findings support a multidomain approach to risk stratification that integrates densitometry, resorption activity and inflammatory tone. Embedded within electronic health records, the model can deliver patient-specific risk with SHAP-based explanations, guiding optimisation of anti-resorptive therapy, calcium-vitamin D repletion and resistance exercise, and prioritising high-risk patients for early intervention (36, 42). Our findings suggest that the combination of bone-specific translational markers (especially CTx), systemic inflammation index (NLR), and lumbar paravertebral muscle health index (PMI) can lead to a multifactorial and comprehensive assessment model of skeletal vulnerability.
In clinical applications, machine learning-based risk prediction tools can help physicians identify at-risk individuals at an early stage and implement personalized interventions for treatment. By comprehensively evaluating individual patients, the model can be embedded into electronic health record systems, providing clinicians with a more intuitive risk assessment to identify women with osteoporosis whose biochemical or imaging parameters indicate fracture vulnerability. Helps develop more effective prevention and treatment strategies. Comprehensive factor analysis helps identify patient deficiencies and provide multidimensional interventions. Early intensive anti-osteoporosis treatment, optimization of calcium-vitamin D status, and targeted resistance training can effectively stop the transition from osteoporosis to bone loss, thereby preventing the first OVCF. For future clinical deployments, our strategy is to reduce the need for bespoke data collection by leveraging data generated from routine care and providing a calibrated, patient-specific probability at the point of care. The production model uses nine routine variables, and a five-feature configuration (BMD T-score, CTx, NLR, age, BMI) retains >90% of the full model’s discriminative power. The tool remains deployable in resource-limited settings, while federated updating and periodic recalibration support sustainable, privacy-preserving model improvement across hospitals. Unlike single-modality calculators (e.g., DEXA thresholds or turnover markers alone) or physics-only morphometric models, our framework integrates densitometry, biochemical resorption activity, inflammatory tone, and clinical context to produce calibrated, patient-specific probabilities with case-level explanations. This multimodal, EHR-embeddable approach supports targeted, modifiable interventions and remains deployable in resource-limited settings via a five-feature configuration that preserves >90% of full-model discrimination.
Although this study explored the risk factors for OVCF in postmenopausal women with osteoporosis through multifactorial analysis and machine learning modelling, there are still some limitations. First, this study was a retrospective study with a single sample source. Future large-scale prospective studies are needed to validate the broad applicability of these risk factors and predictive models. All participants in the study were recruited from China; therefore, genetic backgrounds, lifestyle factors, and medical care may differ from other ethnic groups, which may not make the model generalisable worldwide. Future prospective multicentre validation in European and North American populations will be sought, as well as recalibration of model thresholds for different ethnic backgrounds. In addition, given the low prevalence and incomplete ascertainment of some secondary causes of altered bone metabolism, we treated these factors as exclusion criteria rather than including them as covariates in the model. As a result, our model may not be generalisable to these subgroups and could potentially underestimate the risk in such contexts. Thirdly, although the baseline calibration was favorable, the time-varying treatment exposure (osteoporosis treatment regimen) and adherence were not explicitly modeled. Future studies should explore the influence of more potential factors on secondary OVCF in postmenopausal osteoporotic women to further improve the prediction model.
5 Conclusion
This study identifies low BMD, chronic inflammation, and bone breaking indicator as key risk causes for osteoporotic vertebral compression fractures (OVCF) in postmenopausal women, while regular anti-osteoporotic treatment reduces fracture risk. Multiple machine learning developed provides a reliable tool for personalized OVCF risk prediction, integrating clinical, biological, and musculoskeletal data to enhance prevention and treatment strategies. Future studies are needed to validate the model and explore additional factors, improving its accuracy and clinical application.
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
The studies involving humans were approved by the Ethics Committee of the First Affiliated Hospital of Soochow University (2024 Lun Research Grant No. 728). 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 ethics committee exemption from informed consent.
Author contributions
XS: Resources, Software, Supervision, Validation, Visualization, Writing – review & editing, Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Writing – original draft. PJ: Conceptualization, Formal analysis, Investigation, Software, Supervision, Validation, Visualization, Writing – original draft. YY: Data curation, Formal analysis, Methodology, Software, Validation, Visualization, Writing – review & editing. HS: Conceptualization, Formal analysis, Investigation, Methodology, Software, Validation, Writing – original draft, Writing – review & editing. WT: Data curation, Investigation, Software, Supervision, Visualization, Writing – review & editing, Validation. JM: Formal analysis, Project administration, Validation, Writing – original draft. PZ: Methodology, Supervision, Validation, Visualization, Writing – original draft. QY: Data curation, Investigation, Resources, Software, Supervision, Visualization, Writing – original draft. HY: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing. YQ: Funding acquisition, Resources, Software, Supervision, Validation, Visualization, Writing – review & editing, Writing – original draft.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This project was funded by the Science and Education Strong Health Project of Suzhou in China (No. MSXM2024011).
Acknowledgments
Thank you for all the support from the First Affiliated Hospital of Soochow University.
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 authors declare that no Gen 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.
References
1. Rachner, TD, Khosla, S, and Hofbauer, LC. Osteoporosis: now and the future. Lancet. (2011) 377:1276–87. doi: 10.1016/S0140-6736(10)62349-5
2. Eastell, R, O’Neill, TW, Hofbauer, LC, Langdahl, B, Reid, IR, Gold, DT, et al. Postmenopausal osteoporosis. Nat Rev Dis Primers. (2016) 2:16069. doi: 10.1038/nrdp.2016.69
3. Geusens, P, Feldman, R, Oates, M, Thomas, T, Makras, P, Jakob, F, et al. Romosozumab reduces incidence of new vertebral fractures across severity grades among postmenopausal women with osteoporosis. Bone. (2022) 154:116209. doi: 10.1016/j.bone.2021.116209
4. Hansen, DG, Tutaworn, T, and Lane, JM. What’s new in osteoporosis and fragility fractures. J Bone Joint Surg Am. (2022) 104:1509–15. doi: 10.2106/JBJS.22.00530
5. Franic, D, and Verdenik, I. Risk factors for osteoporosis in postmenopausal women—from the point of view of primary care gynecologist. Zdr Varst. (2018) 57:33–8. doi: 10.2478/sjph-2018-0005
6. Wells, GA, Hsieh, SC, Zheng, C, Peterson, J, Tugwell, P, and Liu, W. Risedronate for the primary and secondary prevention of osteoporotic fractures in postmenopausal women. Cochrane Database Syst Rev. (2022) 2022:CD004523. doi: 10.1002/14651858.CD004523.pub4
7. Lin, YH, and Teng, MMH. Association of possible sarcopenic obesity with osteoporosis and fragility fractures in postmenopausal women. Arch Osteoporos. (2022) 17:65. doi: 10.1007/s11657-022-01107-8
8. Lane, JM, and Witayakom, W. What’s new in osteoporosis and fragility fractures. J Bone Joint Surg Am. (2023) 105:1303–8. doi: 10.2106/JBJS.23.00475
9. Richards, JB, Rivadeneira, F, Inouye, M, Pastinen, TM, Soranzo, N, Wilson, SG, et al. Bone mineral density, osteoporosis, and osteoporotic fractures: a genome-wide association study. Lancet. (2008) 371:1505–12. doi: 10.1016/S0140-6736(08)60599-1
10. Haseltine, KN, Chukir, T, Smith, PJ, Jacob, JT, Bilezikian, JP, and Farooki, A. Bone mineral density: clinical relevance and quantitative assessment. J Nucl Med. (2021) 62:446–54. doi: 10.2967/jnumed.120.256180
11. Qu, L, Zuo, X, Yu, J, Duan, R, and Zhao, B. Association of inflammatory markers with all-cause mortality and cardiovascular mortality in postmenopausal women with osteoporosis or osteopenia. BMC Womens Health. (2023) 23:487. doi: 10.1186/s12905-023-02631-6
12. Tang, Y, Peng, B, Liu, J, Liu, Z, Xia, Y, and Geng, B. Systemic immune-inflammation index and bone mineral density in postmenopausal women: a cross-sectional study of the national health and nutrition examination survey (NHANES) 2007–2018. Front Immunol. (2022) 13:975400. doi: 10.3389/fimmu.2022.975400
13. Chen, S, Sun, X, Jin, J, Zhou, G, and Li, Z. Association between inflammatory markers and bone mineral density: a cross-sectional study from NHANES 2007–2010. J Orthop Surg Res. (2023) 18:305. doi: 10.1186/s13018-023-03795-5
14. Pan, C, Liu, X, Li, T, Wang, G, and Sun, J. Kinetic of bone turnover markers after osteoporotic vertebral compression fractures in postmenopausal female. J Orthop Surg Res. (2018) 13:314. doi: 10.1186/s13018-018-1025-5
15. Kuo, TR, and Chen, CH. Bone biomarker for the clinical assessment of osteoporosis: recent developments and future perspectives. Biomark Res. (2017) 5:18. doi: 10.1186/s40364-017-0097-4
16. Tian, A, Ma, J, Feng, K, Liu, Z, Chen, L, Jia, H, et al. Reference markers of bone turnover for prediction of fracture: a meta-analysis. J Orthop Surg Res. (2019) 14:68. doi: 10.1186/s13018-019-1100-6
17. Ruiz-Gaspà, S, Blanch-Rubió, J, Ciria-Recasens, M, Monfort, J, Tío, L, Garcia-Giralt, N, et al. Reduced proliferation and osteocalcin expression in osteoblasts of male idiopathic osteoporosis. Calcif Tissue Int. (2010) 86:220–6. doi: 10.1007/s00223-010-9336-5
18. Martiniakova, M, Babikova, M, Mondockova, V, Blahova, J, Kovacova, V, and Omelka, R. The role of macronutrients, micronutrients and flavonoid polyphenols in the prevention and treatment of osteoporosis. Nutrients. (2022) 14:523. doi: 10.3390/nu14030523
19. Wang, Z, Wang, H, Zhuang, C, Chen, W, Hoang, TM, Li, J, et al. The effect of type 2 diabetes mellitus on the prognosis of osteoporotic vertebral compression fracture with osteoporotic fracture classification after vertebroplasty. J Orthop Surg Res. (2023) 18:342. doi: 10.1186/s13018-023-03792-8
20. Golsorkhtabaramiri, M, McKenzie, J, and Potter, J. Predictability of neutrophil to lymphocyte ratio in preoperative elderly hip fracture patients for post-operative short-term complications: a retrospective study. BMC Musculoskelet Disord. (2023) 24:227. doi: 10.1186/s12891-023-06211-5
21. Wang, M, Liu, J, Chen, X, Wang, X, and Chen, X. Muscle quality and spine fractures. J Cachexia Sarcopenia Muscle. (2022) 13:1426–8. doi: 10.1002/jcsm.12915
22. Sun, H, Tang, W, You, X, Deng, L, Chen, L, Qian, Z, et al. The role of the lumbar paravertebral muscles in the development of short-term residual pain after lumbar fusion surgery. Spine. (2025) 50:537–47. doi: 10.1097/BRS.0000000000005303
23. Oumer, KS, Liu, Y, Yu, Q, Wu, F, and Yang, S. Awareness of osteoporosis among 368 residents in China: a cross-sectional study. BMC Musculoskelet Disord. (2020) 21:197. doi: 10.1186/s12891-020-03217-1
24. Wang, L, Li, Y, Gu, J, Xiao, L, and Wang, J. Knowledge, awareness and perception towards osteoporosis risk in China: a systematic review. Iran J Public Health. (2024) 53:1009–20. doi: 10.18502/ijph.v53i5.15581
25. Wang, J, Shu, B, Tang, DZ, Li, CG, Xie, XW, Jiang, LJ, et al. The prevalence of osteoporosis in China, a community based cohort study of osteoporosis. Front Public Health. (2023) 11:1084005. doi: 10.3389/fpubh.2023.1084005
26. Hiligsmann, M, Cornelissen, D, Vrijens, B, Abrahamsen, B, al-Daghri, N, Biver, E, et al. Determinants, consequences and potential solutions to poor adherence to anti-osteoporosis treatment: results of an expert group meeting organized by the European Society for Clinical and Economic Aspects of Osteoporosis, Osteoarthritis and Musculoskeletal Diseases (ESCEO) and the International Osteoporosis Foundation (IOF). Osteoporos Int. (2019) 30:2155–65. doi: 10.1007/s00198-019-05104-5
27. Naik-Panvelkar, P, Norman, S, Elgebaly, Z, Elliott, J, Pollack, A, Thistlethwaite, J, et al. Osteoporosis management in Australian general practice: an analysis of current osteoporosis treatment patterns and gaps in practice. BMC Fam Pract. (2020) 21:32. doi: 10.1186/s12875-020-01103-2
28. Lee, BG, Choi, JH, Kim, DY, Choi, WR, Lee, SG, and Kang, CN. Risk factors for newly developed osteoporotic vertebral compression fractures following treatment for osteoporotic vertebral compression fractures. Spine J. (2019) 19:301–5. doi: 10.1016/j.spinee.2018.06.347
29. Mao, W, Dong, F, Huang, G, He, P, Chen, H, Qin, S, et al. Risk factors for secondary fractures to percutaneous vertebroplasty for osteoporotic vertebral compression fractures: a systematic review. J Orthop Surg Res. (2021) 16:644. doi: 10.1186/s13018-021-02722-w
30. Verma, I, and Prasad, SK. Exploring ensemble learning techniques for infant mortality prediction: a technical analysis of XGBoost stacking AdaBoost and bagging models. Birth Defects Res. (2025) 117:e2443. doi: 10.1002/bdr2.2443
31. Looker, AC, Sarafrazi Isfahani, N, Fan, B, and Shepherd, JA. Trends in osteoporosis and low bone mass in older US adults, 2005–2006 through 2013–2014. Osteoporos Int. (2017) 28:1979–88. doi: 10.1007/s00198-017-3996-1
32. Shevroja, E, Cafarelli, FP, Guglielmi, G, and Hans, D. DXA parameters, trabecular bone score (TBS) and bone mineral density (BMD), in fracture risk prediction in endocrine-mediated secondary osteoporosis. Endocrine. (2021) 74:20–8. doi: 10.1007/s12020-021-02806-x
33. Smets, J, Shevroja, E, Hügle, T, Leslie, WD, and Hans, D. Machine learning solutions for osteoporosis—a review. J Bone Miner Res. (2021) 36:833–51. doi: 10.1002/jbmr.4292
34. Gurban, CV, Balaş, MO, Vlad, MM, Caraba, AE, Jianu, AM, Bernad, ES, et al. Bone turnover markers in postmenopausal osteoporosis and their correlation with bone mineral density and menopause duration. Romanian J Morphol Embryol. (2019) 60:1127–35.
35. Vasikaran, S, Eastell, R, Bruyère, O, Foldes, AJ, Garnero, P, Griesmacher, A, et al. Markers of bone turnover for the prediction of fracture risk and monitoring of osteoporosis treatment: a need for international reference standards. Osteoporos Int. (2011) 22:391–420. doi: 10.1007/s00198-010-1501-1
36. Weaver, CM, Gordon, CM, Janz, KF, Kalkwarf, HJ, Lappe, JM, Lewis, R, et al. The National Osteoporosis Foundation’s position statement on peak bone mass development and lifestyle factors: a systematic review and implementation recommendations. Osteoporos Int. (2016) 27:1281–386. doi: 10.1007/s00198-015-3440-3
37. Zhao, Y, Zhang, S, Yi, Y, Qu, T, Gao, S, Lin, Y, et al. Neutrophil-to-lymphocyte ratio as a predictor for cardiovascular diseases: a cohort study in Tianjin. China J Hum Hypertens. (2023) 37:576–83. doi: 10.1038/s41371-022-00724-7
38. Öztürk, ZA, Yesil, Y, Kuyumcu, ME, Bilici, M, Öztürk, N, Yeşil, NK, et al. Inverse relationship between neutrophil lymphocyte ratio (NLR) and bone mineral density (BMD) in elderly people. Arch Gerontol Geriatr. (2013) 57:81–5. doi: 10.1016/j.archger.2013.02.005
39. Vargiu, V, Amar, ID, Rosati, A, Dinoi, G, Turco, LC, Capozzi, VA, et al. Hormone replacement therapy and cervical cancer: a systematic review of the literature. Climacteric. (2021) 24:120–7. doi: 10.1080/13697137.2020.1826426
40. Sipilä, S, Törmäkangas, T, Sillanpää, E, Aukee, P, Kujala, UM, Kovanen, V, et al. Muscle and bone mass in middle-aged women: role of menopausal status and physical activity. J Cachexia Sarcopenia Muscle. (2020) 11:698–709. doi: 10.1002/jcsm.12547
41. Phillips, LK, and Prins, JB. The link between abdominal obesity and the metabolic syndrome. Curr Hypertens Rep. (2008) 10:156–64. doi: 10.1007/s11906-008-0029-7
Keywords: osteoporotic vertebral compression fracture, percutaneous kyphoplasty, postoperative refracture, osteoporosis, postmenopausal women, machine learning
Citation: Sun X, Jing P, Yang Y, Sun H, Tang W, Mi J, Zong P, Yan Q, Yang H and Qiao Y (2025) Risk prediction of osteoporotic vertebral compression fractures in postmenopausal osteoporotic women by machine learning modelling. Front. Med. 12:1664219. doi: 10.3389/fmed.2025.1664219
Edited by:
Q. Wang, Guangzhou University of Chinese Medicine, ChinaReviewed by:
Tawika Kaewchur, Chiang Mai University, ThailandCesar Libanati, Independent Researcher, Fort Lauderdale, FL, United States
Copyright © 2025 Sun, Jing, Yang, Sun, Tang, Mi, Zong, Yan, Yang and Qiao. 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: Qi Yan, eWFucWlzcGluZUAxMjYuY29t; Huilin Yang, aGx5YW5nQHN1ZGEuZWR1LmNu; Yusen Qiao, cWlhb3l1c2VuODYxMkBzdWRhLmVkdS5jbg==
†These authors have contributed equally to this work