Skip to main content
  • Research Article
  • Open access
  • Published:

Establishment of a risk prediction model for residual low back pain in thoracolumbar osteoporotic vertebral compression fractures after percutaneous kyphoplasty



This study aims to identify potential independent risk factors for residual low back pain (LBP) in patients with thoracolumbar osteoporotic vertebral compression fractures (OVCFs) following percutaneous kyphoplasty (PKP) treatment. Additionally, we aim to develop a nomogram that can accurately predict the occurrence of residual LBP.


We conducted a retrospective review of the medical records of thoracolumbar OVCFs patients who underwent PKP treatment at our hospital between July 2021 and December 2022. Residual LBP was defined as the presence of moderate or greater pain (VAS score ≥ 4) in the low back one day after surgery, and patients were divided into two groups: the LBP group and the non-LBP group. These patients were then randomly allocated to either a training or a validation set in the ratio of 7:3. To identify potential risk factors for residual LBP, we employed lasso regression for multivariate analysis, and from this, we constructed a nomogram. Subsequently, the predictive accuracy and practical clinical application of the nomogram were evaluated through a receiver operating characteristic (ROC) curve, a calibration curve, and a decision curve analysis (DCA).


Our predictive model revealed that five variables—posterior fascial oedema, intravertebral vacuum cleft, time from fracture to surgery, sarcopenia, and interspinous ligament degeneration—were correlated with the presence of residual LBP. In the training set, the area under the ROC was 0.844 (95% CI 0.772–0.917), and in the validation set, it was 0.842 (95% CI 0.744–0.940), indicating that the model demonstrated strong discriminative performance. Furthermore, the predictions closely matched actual observations in both the training and validation sets. The decision curve analysis (DCA) curve suggested that the model provides a substantial net clinical benefit.


We have created a novel numerical model capable of accurately predicting the potential risk factors associated with the occurrence of residual LBP following PKP in thoracolumbar OVCFs patients. This model serves as a valuable tool for guiding specific clinical decisions for patients with OVCFs.


Osteoporosis is a systemic metabolic bone disorder characterized by decreased bone mass and the destruction of trabecular bone structure, which often leads to fragility fractures [1]. Osteoporotic vertebral compression fractures (OVCFs), a prevalent complication of osteoporosis-related fragility fractures [2, 3], are frequently associated with chronic and persistent low back pain (LBP), progressive spinal deformity, and other issues, significantly impacting the patients' daily quality of life. In clinical practice, OVCFs can be managed through both conservative and surgical approaches. While conservative treatment of OVCFs has shown relatively high mortality rates, PKP has emerged as an effective treatment [4, 5], providing substantial relief from back pain, spinal stability restoration, early patient mobilization, and prevention of complications such as hypostatic pneumonia and pressure ulcers [6, 7]. However, our clinical observations have revealed that some patients continue to experience residual LBP after undergoing PKP. Previous research studies have reported various potential risk factors associated with the persistence of LBP following PKP [8], including factors such as non-healing bone cement interfaces, thoracolumbar fascia injury (TLFI), cement volume and distribution, and IVC. Some studies have developed risk prediction models for residual LBP after PKP. However, previous models had a limited set of potential independent risk factors. In our study, we conducted a comprehensive analysis of various baseline factors related to postoperative residual LBP in patients with OVCFs. We aimed to develop a predictive model that can be utilized by physicians to inform clinical decision-making for OVCFs patients.



Clinical data of patients with single-segment thoracolumbar OVCFs who underwent a bilateral pedicle approach PKP from July 2021 to December 2022 in our hospital were retrospectively analyzed. The flowchart of our study design is shown in Fig. 1. The specific criteria for inclusion and exclusion were as outlined below. Inclusion criteria: (1) Low back pain due to low-energy injury (e.g., fall, sprain, or lifting a heavy object); (2) X-ray showing a single-segment thoracolumbar vertebral fracture (T10-L2) and prominent high signal intensity on magnetic resonance imaging (MRI) fat-suppressed images of the fractured vertebral body; (3) Dual-energy X-ray absorptiometry (DXA) T-scores of ≤ − 2.5 in the spine/hip; and (4) No symptoms of nerve root and/or spinal cord compression. And exclusion criteria were as follows: Visual analogue scale (VAS) pain score < 5. (2) OVCFs caused by tumor, infection, or tuberculosis, etc. (3) Fracture of the vertebral column in the posterior region, with bone fragments displaced into the spinal canal; and (4) Incomplete follow-up data.

Fig. 1
figure 1

Flowchart of study design

Surgical technique

After the patient was successfully placed under general anesthesia, the patient was positioned in a prone posture. Using the C-arm machine for fluoroscopy, the bilateral arches of the fractured vertebrae were precisely located, and the skin incision sites were marked with a marking pen. Following standard iodine disinfection procedures and the placement of a sterile towel, an incision of approximately 0.5 cm was made at each marked site. Subsequently, a guide needle was carefully inserted into the appropriate position along the outer edge of the vertebral pedicle, allowing the insertion of an expansion cannula to create a working channel. A balloon filled with contrast agent was then positioned to provide support for the vertebral body. Upon reaching the desired height for the vertebral body, the expansion process was halted, and the balloon was carefully removed. At this stage, the bone cement (Mendec, Italy) was introduced slowly into the fractured vertebra. The injection was ceased once the cement had sufficiently diffused to the posterior wall of the vertebral body. Following the solidification of the bone cement, the needle used for the procedure was removed, and the skin incision was covered with sterile gauze, marking the conclusion of the operation. Additionally, lumbar protection through the use of a lumbar brace was routinely employed for one month post-surgery.

Postoperative management

After the surgery, all patients received postoperative anti-osteoporotic treatment [9, 10], which included oral calcium and vitamin D supplementation [11], as well as an annual intravenous infusion of zoledronic acid [12] (Aclasta, 100 ml/5 mg) for a period of three years. And all patients also received biochemical markers of bone turnover (BTMs) testing [13, 14], such as the bone alkaline phosphatase (bALP), procollagen type I N propeptide (PINP). Additionally, patients underwent routine orthopedic radiographic evaluations the day following the surgery to assess the distribution and diffusion of cement within the fractured vertebral body. Subsequent evaluations were conducted at one day, one month, and three months post-surgery to measure the mean severity of LBP and daily dysfunction. To evaluate pain severity and functional status in patients with OVCFs, VAS and Oswestry disability index (ODI) were employed. Residual LBP, as defined in this study, was characterized by a VAS score of ≥ 4 on the first day following surgery. Based on their VAS scores, patients were categorized into a LBP group and a non-LBP group.

Outcome measures

Pre- and postoperative potential risk factors associated with the procedure were extracted from the medical record system, surgical records, radiology image management system, and questionnaires, including (1) demographic characteristics (age, BMI, gender, BMD) and preexisting medical conditions (history of diabetes, history of hypertension, etc.) (2) VAS, ODI, time from fracture to surgery, and preoperative radiological parameters (AVH, Cobb's angle, IVC, posterior fascial oedema, ISLD, and sarcopenia) (3) Postoperative radiographic parameters (cement volume, leakage, cement distribution, and changes in AVH and Cobb angle). All measurements were performed independently by three spine surgeons with 5 years of experience.

The anterior vertebral height (AVH) and Cobb angle of the fractured vertebral body were measured prior to the surgical procedure and one day after surgery. The AVH of the fractured vertebra and the posterior border height of the adjacent normal vertebral body were measured on lateral radiographs using the method described by Teng et al. [15]. Due to the different presentation sizes of the radiographs, the vertebral body height was expressed as relative magnitude: (AVH of the fractured vertebral body /mean of the heights of adjacent posterior vertebral body) × 100% (Fig. 2). To assess the severity of the localized posterior convexity deformity, we used the local kyphotic Cobb’s angle (LKA) (Fig. 2). The definition of posterior fascial oedema was based on MRI (Fig. 3). The type of ISLD on MRI according to the classification of Keorochana et al. was shown in Fig. 4. Vertebral vacuum cleft (IVC) on CT was shown in Fig. 5. For the measurement of sarcopenia, the total psoas area (TPA) was measured based on the lumbar spine CT cross-sections. In brief, the CT image cross-section of the third lumbar vertebral transverse process was reviewed, and the outlines of bilateral psoas major were manually sketched (Fig. 6). The cross-sectional area was measured using the imaging system's software and averaged over three measurements. TPA was determined by dividing the cross-section area by the square of height, expressed as mm2/m2. A diagnosis of sarcopenia was made if the TPA was < 385 mm2/m2 in females or < 545 mm2/m2 in males.

Fig. 2
figure 2

Radiographic evaluation of compressed vertebrae. a Posterior height of the upper vertebrae; b anterior height of the vertebrae; c posterior height of the lower vertebrae; LKA: local kyphotic Cobb’s angle

Fig. 3
figure 3

Posterior fascia oedema was visible on T1-WI (A, orange arrow), T2-WI (B, orange arrow) or STIR images (C, orange arrow) 0f MRI

Fig. 4
figure 4

Example of each grade of interspinous ligament degeneration (orange arrow). A mild (high signal intensity on T1- and T2-WI). B moderate (low signal intensity on T1-WI and high signal intensity on T2-WI). C severe (low or iso-signal intensity on T1- and T2-WI with marked narrowing of the interspinous interval)

Fig. 5
figure 5

Intravertebral vacuum cleft was visible on coronal (A, yellow arrow) and sagittal (B, yellow arrow) views of CT

Fig. 6
figure 6

The TPA at level L3 was measured on CT by tracing the bilateral psoas major muscle outline. A sagittal views of Lumbar CT (L3, yellow line); B cross-sectional views of Lumbar CT (TPA, red fill)

Statistical analysis

SPSS version 22.0 statistical software and R software (version 4.1.2) were used for data processing and generating relevant graphics. If continuous variables followed a normal distribution, they were expressed as  ± S and t-test was used for comparison between groups; otherwise, they were described as median and interquartile range (IQR) and used rank sum test to compare groups. Categorical variables were expressed as n (%), with group comparisons being performed through the Chi-squared (X2) test. The "glmnet" package was employed for the least absolute shrinkage and selection operator (LASSO) regression model, which aided in the selection of the most optimal predictive features [16]. Subsequently, we employed the "glm" package to create a numerical model through multivariate logistic regression analysis. Then, we made use of the "rms" package within the R software to construct the nomogram. Internal validation of the model was carried out through 500 iterations of bootstrap sampling using the "caret" package. ROC curves and calibration curves were plotted to evaluate the model's predictive performance and accuracy by using the "fbroc", "rms" packages. The "rmda" package was utilized to create a DCA plot, assessing the clinical utility of the risk prediction model. Differences were defined as statistically significant at P < 0.05.


In this study, we conducted postoperative follow-up for a total of 267 patients with OVCFs based on our inclusion and exclusion criteria. Among these patients, 46 (17.22%) experienced postoperative residual LBP, while the remaining 221 patients did not. As indicated in Table 1, the variances in preoperative VAS scores and Oswestry Disability Index (ODI) scores between the two groups did not demonstrate statistical significance (P > 0.05). In the postoperative assessments conducted at 1 day, 1 month, and 3 months, a significant difference in both VAS scores and ODI scores between the two groups was observed (P < 0.05).

Table 1 Summary of VAS and ODI after surgery between LBP group and non-LBP group

Baseline patient characteristics in both the training and validation sets

We randomly divided the patients into a training set (186 patients) and a validation set (81 patients) in the ratio of 7:3. The fundamental characteristics of the OVCFs patients in both sets were shown in Table 2. Importantly, no statistically significant differences were observed in any of the variables between these sets (P > 0.05).

Table 2 Comparison of baseline characteristics between training set and validation set

Univariate analysis and variables screening using Lasso regression

The outcomes of the univariate analysis for baseline demographic features, clinical and radiographic factors at baseline, as well as intraoperative factors, were presented in Table 3. When comparing the two groups, significant differences were observed in several variables: posterior fascial oedema (P = 0.022), IVC (P = 0.042), time from fracture to surgery (P = 0.016), sarcopenia (P = 0.001), and ISLD (P = 0.004). However, there were no statistically significant distinctions between the two groups for the remaining variables. To mitigate the influence of multicollinearity, confounding factors, and other issues among these variables, we conducted LASSO regression analysis. Ultimately, we selected five nonzero coefficient variables, as depicted in Fig. 7. These nonzero coefficient variables included posterior fascial oedema, IVC, time from fracture to surgery, sarcopenia, and ISLD.

Table 3 Results of Univariate analysis of variables associated with postoperative residual LBP in patients with OVCFs
Fig. 7
figure 7

Variable selection by the LASSO binary logistic regression model. A LASSO coefficient profiles of the clinical features.  B The optimal penalization coefficient lambda was generated in the LASSO via tenfold cross-validation. We plotted the partial likelihood deviance (binomial deviance) curve versus log(lambda) and drew dotted vertical lines based on 1 standard error criteria. LASSO Least absolute shrinkage and selection operator

Multivariate analysis

A multivariate analysis was conducted using potential predictive features identified through LASSO regression analysis. This analysis revealed that posterior fascial oedema (OR 9.10; 95% CI 2.96–30.73; P = 0.00), IVC (OR 8.67; 95% CI 2.47–32.16; P = 0.001), time from fracture to surgery (OR 6.55; 95% CI 2.52–18.49; P = 0.00), sarcopenia (OR 9.73; 95% CI 3.53–29.67; P = 0.00), and ISLD (OR 5.06, 95% CI 1.92–14.02; P = 0.001) were considered as potential independent risk factors for residual LBP. These findings were summarized in Table 4.

Table 4 Results of multivariate analysis of variables associated with postoperative residual LBP in patients with OVCFs

Predictive model construction

Following the identification of predictors through multivariate analysis, we constructed a nomogram (Fig. 8). Additionally, Receiver Operating Characteristic (ROC) curves were plotted. In the training set, the model's discrimination, measured by the AUC, was 0.844 (95% CI 0.772–0.917), while in the validation set, it was 0.842 (95% CI 0.744–0.940). This indicates that the model exhibited strong predictive accuracy and discrimination. (Fig. 9). To ensure the predictive model's reliability, we performed calibration using calibration plots and Hosmer–Lemeshow tests. The calibration curves exhibited an excellent fit between the predictive model and the validation set. Furthermore, the Hosmer–Lemeshow test demonstrated a high level of agreement between predicted and actual probabilities (training set, P = 0.841; validation set, P = 0.345) (Fig. 10). The Decision Curve Analysis (DCA) results for both the training and validation sets demonstrated a significant net clinical benefit of the numerical model (Fig. 11). For a more comprehensive assessment, we compared the ROC curve of the prediction nomogram with that of the model utilizing single predictors, as depicted in Fig. 12. Remarkably, the AUCs of the single predictors were consistently smaller than those of the predictive model, underscoring the model's robust performance.

Fig. 8
figure 8

Predictive nomogram for residual LBP after percutaneous vertebroplasty (1-Yes, 0-No)

Fig. 9
figure 9

Receiver operating characteristic curve (ROC) validation of the residual LBP risk nomogram prediction. A ROC curve for the prediction model. B ROC curve of the internal validation

Fig. 10
figure 10

Calibration curves of the predictive residual LBP risk nomogram. A calibration curve for the prediction model. B calibration curve of the internal validation

Fig. 11
figure 11

Decision curve analysis for the residual LBP risk nomogram.  A decision curve analysis of the predictive model. B decision curve analysis of the internal validation

Fig. 12
figure 12

Rationality curve analysis for the residual LBP risk nomogram


PKP is a commonly used and effective treatment in clinical practice, surpassing conservative methods in significantly reducing long-term mortality. However, PKP treatment carries certain complications [ [17]], including postoperative adjacent vertebral fractures, pulmonary embolism, cement leakage, and postoperative residual LBP. Among these, the persistence of LBP post-surgery substantially diminishes patients' satisfaction and quality of life. The specific causes of this residual LBP after PKP treatment remain unclear. Our clinical observations have revealed a multitude of factors influencing postoperative LBP. In this study, we delved into the examination of factors impacting postoperative residual LBP, specifically in a total of 267 patients with thoracolumbar OVCFs. Through the development of a predictive model for residual LBP following PKP treatment in patients with thoracolumbar OVCFs, we identified five potential independent risk factors: posterior fascial oedema, IVC, time from fracture to surgery, sarcopenia, and degeneration of the interspinous ligament. Significantly, the ROC analysis demonstrated an AUC of 0.844 (95% CI 0.772–0.917), signifying the model's robust discrimination and calibration. Moreover, the DCA results in both the training and validation sets affirmed the model's substantial net clinical benefit. Moreover, a plausibility analysis of the predictive model again showed that the model was reasonable.

IVC is a condition typically associated with older individuals who have severe osteoporosis, and it is widely believed to be attributed to ischemic necrosis. Furthermore, as noted by Kim [18], the IVC sign was predominantly located at the thoracolumbar junction (81%). Considering the heightened activity and increased mechanical stress experienced by this region, biomechanics was also believed to have a notable impact on the formation of the IVC sign. Therefore, IVC is considered indicative of spinal instability. In the treatment of OVCFs accompanied by IVC, PVP or PKP is often the primary choice [19]. However, several studies have uncovered a strong association between IVC and long-term postoperative issues, such as residual low back pain, vertebral body collapse, and instability [20]. Li et al. [21] demonstrated that patients with OVCFs and IVC typically experience substantially higher VAS scores than those without IVCs. This could be attributed to the cement primarily filling the IVC, potentially leaving the rest of the vertebral body inadequately supported and treated. In alignment with these findings, our results also indicated that the preoperative presence of IVC was a significant risk factor for postoperative residual LBP following PKP.

Furthermore, despite posterior fascial oedema being categorized as a soft tissue injury, it remains a significant consideration in clinical practice. The majority of patients with OVCFs are caused by low-energy injury such as falls. Nevertheless, our study revealed that the presence of posterior fascial oedema in OVCFs patients is not uncommon (19.2%), consistent with previous research findings. Yang [22] and colleagues demonstrated that PKP effectively alleviated LBP in OVCFs patients, but the same relief was not observed in patients with posterior fascial oedema. Their conclusion pointed to a strong connection between posterior fascial oedema and residual LBP. The posterior spinal nerve roots' branches traverse through the thoracolumbar fascia (TLF). When fascial injury occurs, nerve compression, inflammatory factor stimulation, and soft tissue oedema emerge as potential culprits behind the associated pain. Additionally, research on models related to the induction of posterior fascial oedema had consistently shown that reducing oedema and inflammatory factors tended to correlate with gradual pain reduction [23]. Consistent with prior studies, multifactorial analysis underscores the status of posterior fascial oedema as an independent risk factor for post-PKP residual LBP.

Sarcopenia, an age-related condition, is both progressive and widespread [24]. It involves a reduction in muscle mass and function and is often closely linked to osteoporosis. Multiple studies have indicated that nutritional disorders can be the root cause of both sarcopenia and osteoporosis, leading to functional impairment. It is not uncommon for patients to grapple with both conditions [25]. Muscles play a pivotal role in maintaining spinal stability. Sho [26] and colleagues have illustrated how sarcopenia impacts LBP and daily life for osteoporosis patients. Additionally, Yu [27] and team have demonstrated that sarcopenia diminishes the clinical efficacy of PKP. Hence, it is reasonable to hypothesize that sarcopenia may influence the alleviation of LBP following PKP in patients with OVCFs. In this study, we assessed muscle mass using lumbar spine CT. The measurement of the total psoas major muscle area (TPA) at the level of the third lumbar vertebra's transverse process offers a simple and rapid method to evaluate sarcopenia [28, 29]. Our findings indicated a notably higher occurrence of sarcopenia in the residual LBP group when compared to the control group. Multifactorial regression analysis further confirmed that sarcopenia stands as an independent risk factor for post-PKP residual LBP. Consequently, for OVCFs patients with sarcopenia, the focus should extend beyond merely enhancing vertebral stiffness to also encompass muscle mass improvement. In clinical practice, the treatment of sarcopenia recommends a combination of exercise and multiple nutritional supplements, including protein and vitamin D supplementation [30, 31].

We have observed that some patients with OVCFs continue to experience persistent LBP even after a period of conservative treatment, significantly impacting their daily lives. In clinical practice, surgical intervention is often recommended. While numerous previous studies have compared the effectiveness of conservative treatment and vertebral augmentation, little attention has been given to the timing of vertebral augmentation. Akihito [32] and colleagues discovered that early surgical treatment in OVCFs patients reduced the incidence of postoperative LBP and re-fractures. However, for those OVCFs patients who still experience residual LBP after surgery, some studies have concluded that the timing of vertebral augmentation is not an independent risk factor [33]. In our study, we utilized a 4-week period from injury to surgery as a reference point and found that the proportion of patients with an injury-to-surgery interval exceeding 4 weeks was higher in the group with residual LBP compared to the control group. This led us to propose that the timing of vertebral augmentation may indeed be a significant factor influencing the development of residual LBP following surgery in OVCFs patients (OR 6.554, P = 0.00). This could be attributed to various factors: (1) It might be associated with the delayed formation of osseous junctions or pseudoarthrosis, often stemming from the delayed treatment of compression fractures. Such delays are known to contribute to chronic pain that ensues after an injury [34,35,36]. (2) The formation of primitive callus typically takes place around four weeks post-fracture and undergoing vertebral kyphoplasty at a later date (more than 4 weeks) is not conducive to alignment improvement.

Degenerative changes in the lumbar spine rank among the primary causes of LBP and disability in the elderly [37]. Typically, we place our focus on alterations in intervertebral disks while overlooking the significance of non-disk structures within the spine, including the facet joints, paraspinal muscles, and spinal ligaments. Among these structures, the interspinous ligament holds a pivotal position within the posterior spinal ligament complex. It plays a crucial role in spinal flexion. When the body undergoes excessive forward bending or experiences an energetic impact, the posterior column structures of the lumbar spine come under increased tension and load, with the interspinous ligament often bearing the brunt of the injury. Multiple studies have assessed the interspinous ligament's role through biomechanical tests [38] and through the observation of anatomical, biochemical, and pathological changes in degenerative spines. Maes [39] and colleagues have concluded that the degeneration of posterior column ligaments is a significant contributor to LBP. In the investigation of changes within the interspinous ligament, MRI has established itself as the primary diagnostic method of choice. Keorochana et al. categorized ISLD into four types [40, 41]. However, our study encompassed cases beyond mild degeneration. In line with prior research, multifactorial analysis demonstrated that ISLD stands as a significant risk factor for persistent LBP following PKP (OR 5.061, P = 0.001). Consequently, we created a nomogram through the analysis of potential risk factors to predict the occurrence of residual LBP. This model proves to be a viable tool for risk assessment, enabling early identification of OVCFs patients with a high risk of residual LBP. With this knowledge, individualized strategies can be crafted for these patients during the perioperative period.

Our study has several limitations. First, it is important to acknowledge that our study is retrospective and has a restricted sample size. Second, our postoperative follow-up period is relatively short, and we acknowledge the importance of conducting studies with extended follow-up durations and multiple assessment time points. Third, it is essential to highlight that this study is conducted at a single center, and the validation of our findings would benefit from future multicenter studies conducted in different hospitals.


We have created a novel numerical model that can precisely forecast the potential risk factors associated with the onset of LBP following PKP in patients with thoracolumbar OVCFs. This model is poised to aid in the early identification of OVCFs patients at a heightened risk of enduring residual LBP, offering valuable insights to guide clinical decision-making.

Availability of data and materials

The authors will provide the raw data supporting the conclusions of this paper without reservation.



Osteoporotic vertebral compression fractures


Percutaneous kyphoplasty


Low back pain


Receiver operating characteristic


Decision curve analysis


Thoracolumbar fascia injury


Intravertebral vacuum cleft


Interspinous ligament degeneration


Biochemical markers of bone turnover


Bone alkaline phosphatase


Procollagen type I N propeptide


Oswestry disability index


Visual analogue scale


Anterior vertebral height


Local kyphotic angle


Magnetic resonance imaging


Total psoas area


  1. Migliorini F, Giorgino R, Hildebrand F, et al. Fragility fractures: risk factors and management in the elderly. Med Kaunas Lith. 2021;57(10):1119.

    Article  Google Scholar 

  2. Tang B, Zeng H, Hu S, Liu K, Wu L, Shi X. Percutaneous vertebroplasty combined with zoledronic acid in treatment and prevention of osteoporotic vertebral compression fractures: a systematic review and meta-analysis of comparative studies. World Neurosurg. 2022;157:75–87.

    Article  PubMed  Google Scholar 

  3. Osterhoff G, Asatryan G, Spiegl UJA, Pfeifle C, Jarvers JS, Heyde CE. Impact of multifidus muscle atrophy on the occurrence of secondary symptomatic adjacent osteoporotic vertebral compression fractures. Calcif Tissue Int. 2022;110(4):421–7.

    Article  PubMed  CAS  Google Scholar 

  4. Vendeuvre T, Brossard P, Pic JB, et al. Vertebral balloon kyphoplasty versus vertebral body stenting in non-osteoporotic vertebral compression fractures at the thoracolumbar junction: a comparative radiological study and finite element analysis (BONEXP study). Eur Spine J Off Publ Eur Spine Soc Eur Spinal Deform Soc Eur Sect Cerv Spine Res Soc. 2021;30(10):3089–98.

    Article  Google Scholar 

  5. Ning L, Zhu J, Tian S, et al. Correlation analysis between basic diseases and subsequent vertebral fractures after percutaneous kyphoplasty (PKP) for osteoporotic vertebral compression fractures. Pain Physician. 2021;24(6):E803–10.

    Article  PubMed  Google Scholar 

  6. Rabei R, Patel K, Ginsburg M, et al. Percutaneous vertebral augmentation for vertebral compression fractures: national trends in the medicare population (2005–2015). Spine. 2019;44(2):123–33.

    Article  PubMed  Google Scholar 

  7. Balkarli H, Kilic M, Balkarli A, Erdogan M. An evaluation of the functional and radiological results of percutaneous vertebroplasty versus conservative treatment for acute symptomatic osteoporotic spinal fractures. Injury. 2016;47(4):865–71.

    Article  PubMed  Google Scholar 

  8. Li Y, Yue J, Huang M, et al. Risk factors for postoperative residual back pain after percutaneous kyphoplasty for osteoporotic vertebral compression fractures. Eur Spine J Off Publ Eur Spine Soc Eur Spinal Deform Soc Eur Sect Cerv Spine Res Soc. 2020;29(10):2568–75.

    Article  Google Scholar 

  9. Migliorini F, Maffulli N, Colarossi G, Eschweiler J, Tingart M, Betsch M. Effect of drugs on bone mineral density in postmenopausal osteoporosis: a Bayesian network meta-analysis. J Orthop Surg. 2021;16(1):533.

    Article  Google Scholar 

  10. Migliorini F, Colarossi G, Eschweiler J, Oliva F, Driessen A, Maffulli N. Antiresorptive treatments for corticosteroid-induced osteoporosis: a Bayesian network meta-analysis. Br Med Bull. 2022;143(1):46–56.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Conti V, Russomanno G, Corbi G, et al. A polymorphism at the translation start site of the vitamin D receptor gene is associated with the response to anti-osteoporotic therapy in postmenopausal women from southern Italy. Int J Mol Sci. 2015;16(3):5452–66.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Migliorini F, Colarossi G, Baroncini A, Eschweiler J, Tingart M, Maffulli N. Pharmacological management of postmenopausal osteoporosis: a level I evidence based - expert opinion. Expert Rev Clin Pharmacol. 2021;14(1):105–19.

    Article  PubMed  CAS  Google Scholar 

  13. Migliorini F, Maffulli N, Spiezia F, Peretti GM, Tingart M, Giorgino R. Potential of biomarkers during pharmacological therapy setting for postmenopausal osteoporosis: a systematic review. J Orthop Surg. 2021;16(1):351.

    Article  Google Scholar 

  14. Migliorini F, Maffulli N, Spiezia F, Tingart M, Maria PG, Riccardo G. Biomarkers as therapy monitoring for postmenopausal osteoporosis: a systematic review. J Orthop Surg. 2021;16(1):318.

    Article  Google Scholar 

  15. Teng MMH, Wei CJ, Wei LC, et al. Kyphosis correction and height restoration effects of percutaneous vertebroplasty. AJNR Am J Neuroradiol. 2003;24(9):1893–900.

    PubMed  PubMed Central  Google Scholar 

  16. Hu X, Shen F, Zhao Z, Qu X, Ye J. An individualized gait pattern prediction model based on the least absolute shrinkage and selection operator regression. J Biomech. 2020;112:110052.

    Article  PubMed  Google Scholar 

  17. Ma Y, Wu X, Xiao X, et al. Effects of teriparatide versus percutaneous vertebroplasty on pain relief, quality of life and cost-effectiveness in postmenopausal females with acute osteoporotic vertebral compression fracture: a prospective cohort study. Bone. 2020;131:115154.

    Article  PubMed  CAS  Google Scholar 

  18. Kim DY, Lee SH, Jang JS, Chung SK, Lee HY. Intravertebral vacuum phenomenon in osteoporotic compression fracture: report of 67 cases with quantitative evaluation of intravertebral instability. J Neurosurg. 2004;100(1 Suppl Spine):24–31.

    Article  PubMed  Google Scholar 

  19. Wei P, Yao Q, Xu Y, Zhang H, Gu Y, Wang L. Percutaneous kyphoplasty assisted with/without mixed reality technology in treatment of OVCF with IVC: a prospective study. J Orthop Surg. 2019;14(1):255.

    Article  Google Scholar 

  20. Ito Y, Hasegawa Y, Toda K, Nakahara S. Pathogenesis and diagnosis of delayed vertebral collapse resulting from osteoporotic spinal fracture. Spine J Off J North Am Spine Soc. 2002;2(2):101–6.

    Article  Google Scholar 

  21. Li Z, Liu T, Yin P, et al. The therapeutic effects of percutaneous kyphoplasty on osteoporotic vertebral compression fractures with or without intravertebral cleft. Int Orthop. 2019;43(2):359–65.

    Article  PubMed  Google Scholar 

  22. Luo Y, Jiang T, Guo H, Lv F, Hu Y, Zhang L. Osteoporotic vertebral compression fracture accompanied with thoracolumbar fascial injury: risk factors and the association with residual pain after percutaneous vertebroplasty. BMC Musculoskelet Disord. 2022;23(1):343.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Schilder A, Hoheisel U, Magerl W, Benrath J, Klein T, Treede RD. Sensory findings after stimulation of the thoracolumbar fascia with hypertonic saline suggest its contribution to low back pain. Pain. 2014;155(2):222–31.

    Article  PubMed  Google Scholar 

  24. Cruz-Jentoft AJ, Baeyens JP, Bauer JM, et al. Sarcopenia: European consensus on definition and diagnosis: report of the European Working Group on Sarcopenia in Older People. Age Ageing. 2010;39(4):412–23.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Balagopal P, Ljungqvist O, Nair KS. Skeletal muscle myosin heavy-chain synthesis rate in healthy humans. Am J Physiol. 1997;272(1 Pt 1):E45-50.

    Article  PubMed  CAS  Google Scholar 

  26. Iwahashi S, Hashida R, Matsuse H, et al. The impact of sarcopenia on low back pain and quality of life in patients with osteoporosis. BMC Musculoskelet Disord. 2022;23(1):142.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Peng Y, Wu X, Ma X, Xu D, Wang Y, Xia D. Comparison between the clinical effect of percutaneous kyphoplasty for osteoporosis vertebral compression fracture patient with or without sarcopenia: a retrospective cohort study. Int J Gen Med. 2023;16:3095–103.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Jones KI, Doleman B, Scott S, Lund JN, Williams JP. Simple psoas cross-sectional area measurement is a quick and easy method to assess sarcopenia and predicts major surgical complications. Colorectal Dis Off J Assoc Coloproctol G B Irel. 2015;17(1):O20-26.

    Article  CAS  Google Scholar 

  29. Yamazaki Y, Kanaji S, Takiguchi G, et al. Skeletal muscle loss after laparoscopic gastrectomy assessed by measuring the total psoas area. Surg Today. 2020;50(7):693–702.

    Article  PubMed  Google Scholar 

  30. Kwon YN, Yoon SS. Sarcopenia: neurological point of view. J Bone Metab. 2017;24(2):83–9.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Tintignac LA, Brenner HR, Rüegg MA. Mechanisms regulating neuromuscular junction development and function and causes of muscle wasting. Physiol Rev. 2015;95(3):809–52.

    Article  PubMed  CAS  Google Scholar 

  32. Minamide A, Maeda T, Yamada H, et al. Early versus delayed kyphoplasty for thoracolumbar osteoporotic vertebral fractures: the effect of timing on clinical and radiographic outcomes and subsequent compression fractures. Clin Neurol Neurosurg. 2018;173:176–81.

    Article  PubMed  Google Scholar 

  33. Lin M, Wen X, Huang Z, et al. A nomogram for predicting residual low back pain after percutaneous kyphoplasty in osteoporotic vertebral compression fractures. Osteoporos Int J Establ Result Coop Eur Found Osteoporos Natl Osteoporos Found USA. 2023;34(4):749–62.

    Article  Google Scholar 

  34. Kanchiku T, Imajo Y, Suzuki H, Yoshida Y, Taguchi T. Usefulness of an early MRI-based classification system for predicting vertebral collapse and pseudoarthrosis after osteoporotic vertebral fractures. J Spinal Disord Tech. 2014;27(2):E61-65.

    Article  PubMed  Google Scholar 

  35. Wu CT, Lee SC, Lee ST, Chen JF. Classification of symptomatic osteoporotic compression fractures of the thoracic and lumbar spine. J Clin Neurosci Off J Neurosurg Soc Australas. 2006;13(1):31–8.

    Article  CAS  Google Scholar 

  36. Tsujio T, Nakamura H, Terai H, et al. Characteristic radiographic or magnetic resonance images of fresh osteoporotic vertebral fractures predicting potential risk for nonunion: a prospective multicenter study. Spine. 2011;36(15):1229–35.

    Article  PubMed  Google Scholar 

  37. Kelsey JL, White AA. Epidemiology and impact of low-back pain. Spine. 1980;5(2):133–42.

    Article  PubMed  CAS  Google Scholar 

  38. Iida T, Abumi K, Kotani Y, Kaneda K. Effects of aging and spinal degeneration on mechanical properties of lumbar supraspinous and interspinous ligaments. Spine J Off J North Am Spine Soc. 2002;2(2):95–100.

    Article  Google Scholar 

  39. Maes R, Morrison WB, Parker L, Schweitzer ME, Carrino JA. Lumbar interspinous bursitis (Baastrup disease) in a symptomatic population: prevalence on magnetic resonance imaging. Spine. 2008;33(7):E211-215.

    Article  PubMed  Google Scholar 

  40. Keorochana G, Taghavi CE, Tzeng ST, et al. MRI classification of interspinous ligament degeneration of the lumbar spine: intraobserver and interobserver reliability and the frequency of disagreement. Eur Spine J Off Publ Eur Spine Soc Eur Spinal Deform Soc Eur Sect Cerv Spine Res Soc. 2010;19(10):1740–5.

    Article  Google Scholar 

  41. Keorochana G, Taghavi CE, Tzeng ST, et al. Magnetic resonance imaging grading of interspinous ligament degeneration of the lumbar spine and its relation to aging, spinal degeneration, and segmental motion. J Neurosurg Spine. 2010;13(4):494–9.

    Article  PubMed  Google Scholar 

Download references


Not applicable.


There was no direct funding source aligned to this study.

Author information

Authors and Affiliations



WT and YN conceived and designed the study. Data collection and manuscript writing were undertaken by WT and DL. YN, PS and FL analyzed and interpreted the data. YS critically reviewed and revised the manuscript. All authors played substantial roles in contributing to and approving the final version of the manuscript.

Corresponding author

Correspondence to Fanguo Lin.

Ethics declarations

Ethics approval and consent to participate

We obtained ethical approval for this retrospective study from the Medical Ethics Committee of the Second Affiliated Hospital of Soochow University (No. JD-HG-2023-94). All patients' medical information was anonymized, thus waiving the requirement for informed consent. We conducted our study in compliance with the principles outlined in the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that the study was conducted in the absence of any business or financial relationships that could be interpreted as potential conflicts of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tu, W., Niu, Y., Su, P. et al. Establishment of a risk prediction model for residual low back pain in thoracolumbar osteoporotic vertebral compression fractures after percutaneous kyphoplasty. J Orthop Surg Res 19, 41 (2024).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: