Authors: Yossi Smorgick, Ruth Pelleg‐Kallevag, Dror Lindner, Yoram Anekstein, Sergey Goldstein, Hila May
Categories: Communication, aging, bone mineral density, fracture progression, osteoporosis, vertebra compression fracture, vertebral body collapse, vertebral body density
Source: Clinical Anatomy (New York, N.y.)
Doi: 10.1002/ca.24219
Vertebral osteoporotic fractures (VOF) are among the most frequent fractures in the elderly, often leading to an impaired lifestyle and a high economic burden. Although a reduced bone mass density is considered one of the main risk factors for VOF, its role in determining the fracture type, using the AO spine‐DGOU classification for osteoporotic thoracolumbar fractures, as well as its progression, is unknown. The current study aimed (1) reveal whether the bone density of the vertebral bodies of fractured and non‐fractured vertebrae predicts the type of fracture, (2) examine whether bone density is associated with the initial and progressive collapse of the vertebral body, and (3) provide predictive measures for fracture progression. The study sample included 124 patients (40 males and 84 females) with an acute osteoporotic vertebral fracture who underwent a computerized tomography scan at the time of diagnosis and an x‐ray at least 3 months later. The bone density of the fractured and adjacent (non‐fractured) vertebrae was measured at diagnosis. The magnitude of the collapse and the progression of the fracture over time were calculated from height measurements of the vertebral bodies at diagnosis and follow‐up. Age was a significant factor in predicting the fracture type and magnitude of collapse, whereas sex and bone density were not. The severity of the fracture was involved in predicting its progression, demonstrating that severe‐type fractures tended to continue to collapse after diagnosis. However, when each type was examined independently, the density of the fractured vertebra had a protective effect on fracture progression. To conclude, identifying the type of fracture is beneficial in determining patient prognosis. Furthermore, the density of the fractured vertebra, the magnitude of collapse, and patient age are valuable predictors of fracture progression.
Keywords: aging, bone mineral density, fracture progression, osteoporosis, vertebra compression fracture, vertebral body collapse, vertebral body density
Vertebral osteoporotic fractures (VOF) are among the major health concerns nowadays, especially with the increase in life expectancy (Ballane et al., 2017; Holroyd et al., 2008). These fractures commonly lead to an impaired quality of life with an increased mortality risk and a high economic burden on the health system and the individual.
Osteoporosis is characterized by a reduction in bone mass and alteration in bone architecture, resulting in increased bone fragility and fracture risk (Ensrud, 2017). The World Health Organization defines osteoporosis as a bone mineral density (BMD) value of more than 2.5 standard deviations below the mean for a young white female reference group. Dual‐energy x‐ray absorptiometry (DEXA) is the gold standard for measuring BMD, and low BMD is considered a strong risk factor for sustaining a fracture (Blake & Fogelman, 2007). Consequently, drug development is directed at improving BMD to reduce the risk of new osteoporosis‐related fractures (Nogués & Martinez‐Laguna, 2018).
Among the most common osteoporotic fractures are those of the spine. These fractures include several types, which were classified by the German Society for Orthopedics and Trauma (DGOU) in the AO Spine‐DGOU classification for osteoporotic Thoracolumbar Fractures (Schnake et al., 2018). Nevertheless, it is unclear whether this classification method has a prognostic value. Moreover, BMD contributes to predicting the risk of osteoporotic vertebral fracture (Cann et al., 1985; Wang et al., 2012), whereby improved BMD reduces the risk of vertebral fractures in osteoporotic patients (Bouxsein et al., 2019; Cummings et al., 2002). However, the effect of BMD on vertebral fracture type and progression (i.e., deterioration) over time is unclear.
Previous studies have shown that bone density in Hounsfield units (HU) obtained from computed tomography (CT) can predict BMD (Aydin Ozturk et al., 2021; Lee et al., 2013; Pinto et al., 2022; Schreiber et al., 2011; Schreiber et al., 2014). Furthermore, HU measurements may be superior to DEXA in evaluating spinal bone quality since they are less influenced by degenerative changes, vascular calcifications, and spinal deformity (Pinto et al., 2022). Recently, several studies have evaluated HU measurements to predict proximal junctional failure (Meredith et al., 2013), pedicle screw loosening (Bredow et al., 2016; Schwaiger et al., 2014), osteoporotic vertebral fractures (Schwaiger et al., 2014), cage subsidence (Mi et al., 2017; Xi et al., 2020), and vertebral collapse fractures in multiple myeloma patients (Zijlstra et al., 2020). Studying the correlation between BMD, as reflected in HU measurements in patients with osteoporotic compression spine fractures, may elucidate whether BMD influences the fracture type and corresponds to the definition in the AO Spine‐DGOU Classification for Osteoporotic Thoracolumbar Fractures (Schnake et al., 2018).
The major aims of this study (1) To determine whether significant differences exist in the bone density (measured in HU) of fractured and non‐fractured vertebral bodies between the various types of osteoporotic vertebral fractures according to the definition of the AO Spine‐DGOU Classification for Osteoporotic Thoracolumbar Fractures (Schnake et al., 2018). (2) To examine whether vertebral body bone density (in HU) or fracture classification is associated with initial and progressive collapse of the vertebral body following acute osteoporotic compression spine fractures. (3) To provide predictive measures for fracture progression.
We hypothesized that the BMD of non‐fractured vertebral bodies would play a significant role in determining the fracture type; namely, a lower BMD would be associated with a more severe fracture. However, we hypothesized the opposite for the fractured vertebra. That is, a higher BMD will be associated with a more severe fracture, and that BMD will be negatively associated with the magnitude of the collapse. Therefore, we hypothesized that BMD and the magnitude of vertebral collapse would serve as predictive variables for fracture progression.
The study included 124 patients (40 males and 84 females between 50 and 93 years) with acute osteoporotic compression spine fractures who underwent a computerized tomography (CT) scan between 2008 and 2013. The study was approved by the Institutional Review Board.
Patients were included in the study if they were 50 years or older and had an acute osteoporotic compression spine fracture. The diagnosis of acute fracture was based on the clinical presentation of acute pain, usually after minimal trauma, typical appearance on the CT scan with new onset fracture lines, and comparison to old radiographs when available (Smorgick et al., 2020). All patients included in the study underwent a CT scan at the time of diagnosis and plain radiography at least 3 months later. The mean length of follow‐up was 17.5 months, with a range of 3 to 70 months. All patients were treated with a spinal orthosis for 3 months. Patients were excluded if they presented with pathological fractures, transverse or spinous process fractures, and fractures due to high‐energy trauma. Patients with old fractures and those with a radiological follow‐up of less than 3 months were also excluded from the study. Consequently, of 305 patients, 34 were excluded due to having a transverse process or a spinous process fracture, 34 due to high‐energy trauma, 44 did not meet the age criterion (i.e., were <50 years old), 21 had old fractures, and five were excluded due to having pathological fractures. Of the 167 remaining patients, 43 did not have the required three‐month follow‐up radiograph.
Vertebral body fractures were classified by a spine surgeon according to the AO spine‐DGOU classification for osteoporotic thoracolumbar fractures (Schnake et al., 2018). Based on the AO spine‐DGOU classification for osteoporotic thoracolumbar fractures, we divided the fractures into three major types, which were also the most frequent (95% of the fractures) (Schnake et al., 2018):
OF1 was excluded from our fracture criteria because it can be diagnosed only via MRI images and OF5 was excluded because it is surgically treated.
The vertebral body height (mm) was measured using the Sectra PACS (v. 23.2.6, 2022, Sweden). Vertebral body density (HU) and homogeneity of the density (HU) were measured using the IntelliSpace Portal, Philips (v. 11.2) (Table 1). All measurements were obtained for the fractured vertebra and adjacent (above and below) non‐fractured vertebrae. Vertebral body height was measured at diagnosis and at the follow‐up examination. Additional variables were calculated from the height measurements, including the collapse magnitude at diagnosis and follow‐up and the progress (deterioration) of collapse over time (i.e., the delta between the values at diagnosis and follow‐up examination). The density and homogeneity measurements were obtained only from CT scans performed at diagnosis (Table 1).
Interobserver variation was performed for all variables measured. Two independent researchers blindly performed the measurements. Intraclass correlation coefficient (ICC) analysis revealed excellent results (0.942 ≤ ICC ≤ 0.968).
Differences in age and vertebral characteristics between males and females were examined using independent sample t‐tests. A chi‐square test was performed to determine the association between sex and fracture type. Spearman or Pearson correlations (according to variable distribution, SI Table 1) were carried out between age and other linear measurements. Kruskal‐Wallis tests were performed to examine the differences in age, follow‐up duration, and vertebral variables between fracture types. Paired‐sample t‐tests were performed to examine differences in variables between the fractured vertebra and the average calculated from the adjacent vertebrae. Multinomial logistic regression and logistic regression (with the Forward LR method using dummy variables) were performed to determine the predictive variables for each fracture type. Only measurements measured at the time of fracture diagnosis were included in these analyses. A linear regression (Forward method) was performed to examine the predictive variables for the deterioration in collapse between diagnosis and follow‐up, for the combined fracture types as well as for each condition separately. In these analyses, variables measured at the time of diagnosis and the duration of the follow‐up were included.
The type of fracture was found to be sex‐independent (χ^2^ = 2.960, df = 2, p = 0.228). Similar frequencies of OF2 were found for males and females (35.0% and 33.3%, respectively). OF3 was slightly more frequent among females (46.4% vs. 32.5%), whereas OF4 was more frequent among males (32.5% vs. 20.2%). Nevertheless, fracture type was age‐dependent (p = 0.005). Accordingly, the mean age of patients with fracture type OF3 and OF4 was higher (71.3 ± 9.41 and 70.1 ± 10.83 years, respectively) than those with a fracture type OF2 (65.0 ± 8.73 years). Therefore, further analyses were performed while controlling for age. No significant differences were found in the duration of the follow‐up between different fracture types (p = 0.230).
A significant difference in the average height of the non‐fractured vertebrae was found between males and females, yet not in any other characteristic of the vertebrae nor in the deterioration of collapse over time (Table 2). However, significant associations were found between the different vertebral variables and age. At diagnosis, associations with age were found only for vertebral density, especially with the average density of non‐fractured vertebrae (r = −0.514, p < 0.001). At follow‐up, all variables measured were significantly associated with age, with the strongest correlation between age and the height of the fractured vertebra (r = −0.427, p < 0.001), followed by the average height of the adjacent vertebrae and the magnitude of the collapse (r = −0.323 and r = 0.349, p < 0.001, respectively). The deterioration in the magnitude of the vertebral collapse was also significantly correlated with age (r = 0.414, p < 0.001) (SI Table 2).
The mean height of the fractured vertebrae of the group with an OF4 fracture type was significantly lower than that of the other two groups, and the collapse was significantly larger (Table 3). In the follow‐up examination, significant differences in the fractured vertebral body height were found between all fracture types. The magnitude of collapse increased over time and differed significantly between all fracture types (Table 3 and Figure 1). Furthermore, significant differences were found in vertebral body collapse between the diagnosis and the follow‐up, whereby the deterioration over time was the largest in the OF4 type, followed by the OF3 type. The smallest progression was evident in the OF2 type, reaching significance when compared with the other two groups (Table 3 and Figure 1).
FIGURE 1 (a) Collapse magnitude (%) of fractured vertebrae at diagnosis and follow‐up by the fracture type. (b) Deterioration in the magnitude of the vertebral collapse (%) between the time of diagnosis and follow‐up by fracture type. Only individuals over 60 years old were included (N = 98).
The fractured vertebra was significantly denser than the average density calculated from the adjacent non‐fractured vertebrae, independent of fracture type (p < 0.01). However, a tendency of increment in the density of the fractured vertebra was found between OF2 and OF4 types, with a significant reduction in the density homogeneity (SD) in the OF4 type (Table 3). Nevertheless, in the adjacent non‐fractured vertebrae, no significant differences in bone density or homogeneity were found between fracture types (Table 3).
Following the multinomial logistic regression, age and vertebrae height (fractured and the average of the adjacent ones) differentiated between OF 2 and OF3. Yet, only age differentiated between OF2 and OF4. Nevertheless, sex, density of the fractured and average of the adjacent vertebrae, density homogeneity of the fractured vertebra, and rate of collapse at diagnosis differentiated between OF3 and OF4 (Table 4). In the binary logistic regression, different variables were included to predict each type of fracture (Table 4). The probability of manifesting the OF2 type decreased with increasing age (Exp(B) = 0.922), as well as with the increase in the magnitude of collapse at diagnosis (Exp(B) = 0.902). The probability of manifesting the OF3 type increased slightly with increasing age (Exp(B) = 1.050), whereas the OF4 type manifestation increased with increasing magnitude of collapse at diagnosis (Exp(B) = 1.112), regardless of the age of the individual (excluded from the discriminant function).
The main factor that predicted the progression of compression fracture over time (i.e., deterioration) was the fracture type; the more severe the fracture, the more it tended to progress. However, the predictive factors for each fracture type vary. The progression of collapse among patients with an OF2 fracture was larger with increased age; however, the collapse magnitude at diagnosis served as a protective factor (the larger the collapse, the smaller the progression). Among patients with OF3, the density of the fractured vertebra served as a protective factor against the progression of collapse over time. In patients with OF4, increased age and reduced density homogeneity advanced the progression of the collapse, whereas increased vertebral body density was a protective factor (Table 5).
Osteoporotic vertebral fractures are multifactorial and vary in definition and diagnosis (Ballane et al., 2017). These characteristics hinder their prediction and prevention. This study is the first to examine the predictive value of the AO Spine‐DGOU Classification for Osteoporotic Thoracolumbar Fractures (Schnake et al., 2018) for progressive collapse. It is also the first study to investigate the role of vertebral body density (in HU) in determining the type of osteoporotic fracture defined by the AO spine‐DGOU classification for osteoporotic thoracolumbar fractures (Schnake et al., 2018), as well as in predicting the progression of vertebral body collapse. Interestingly, our findings refuted our hypothesis that the density of non‐fractured vertebrae will predict the severity of the fracture and will help predict its progression. However, our hypotheses that the density of the fractured vertebra and its magnitude of collapse predict the type of fracture, as well as its progression, were validated.
Our study supports previous findings indicating that the mean BMD of fractured vertebrae, including the old and recent fracture subgroups, was significantly higher than that of the non‐fractured vertebrae (Hsu et al., 2021). In addition, we found that when the fracture was more severe, the density of the fractured vertebra was higher, although it did not reach significance. A possible explanation for the high density of fractured vertebrae in general, and with a more severe manifestation, might be the degree of compression. A more compressed vertebra involves less air space between the trabeculae. Moreover, vertebral body density homogeneity was significantly lower in the OF4 fracture type. The reduced homogeneity in OF4 could be due to the various forms of fractures that are included in this category, as well as their comminution nature (Schnake et al., 2018).
Similar to previous studies (Pinto et al., 2022), and probably due to the higher incidence of osteoporosis in females compared to males (Holroyd et al., 2008), the number of females with an osteoporotic vertebral fracture in our study cohort was more than twice of males. However, there was no significant difference between males and females concerning the fracture type, especially in stable impression fractures (OF2). As opposed to OF2, both OF3 and OF4 are considered unstable fractures (Schnake et al., 2018). These differences in the fracture type stability were evident in our results. Significant differences in the deterioration in collapse magnitude (between the time of diagnosis and follow‐up) were found between OF2 and the other two types, with no significant differences between the OF3 and OF4 types. However, in each fracture type, a larger compression had a protective effect on the progression of the fracture. It is noteworthy that males and females exhibited a similar progression of vertebral collapse following acute osteoporotic compression spine fracture.
Advanced age significantly contributed to a more severe manifestation of fracture type and its progression over time. This is probably due to the reduction in the number, thickness, and interconnectivity of trabeculae in the bones of older patients with osteoporosis (Rao & Singrakhia, 2003). The deterioration in bone quality with age is supported by our study. We found a strong and significant negative correlation between age and bone density of the non‐fractured vertebrae. However, the density of the non‐fractured vertebra did not correlate significantly with fracture progression over time.
Although osteoporosis manifestation and its severity increase the risk of vertebral fracture (Holroyd et al., 2008), according to our results, it did not increase the risk of a specific fracture type since we did not find any significant difference in the density of the non‐fractured vertebrae between the fracture types. Hence, other explanatory factors that affect the magnitude of the collapse of fractured vertebrae should be considered. As opposed to hip fractures, which usually occur due to a fall, vertebral fractures can result from various activities, including lifting, bending, and reaching (Holroyd et al., 2008; Myers & Wilson, 1997). Therefore, the immediate cause of the fracture (e.g., trauma or specific activity) and vertebral characteristics, such as the cross‐sectional area, might have a larger effect on determining the type of vertebral fracture (Gilsanz et al., 1995; Holroyd et al., 2008; Myers & Wilson, 1997; Silman et al., 1997).
This study has several limitations, including the age distribution of the patients; therefore, for some of the analyses, we had to control for age by studying a subsample. However, for most analyses, we examined the effect of age on the fracture type or severity of deterioration over time. Another limitation was the inability to measure bone density at the follow‐up examination, as only X‐rays were available.
In conclusion, identifying the type of fracture using the AO spine‐DGOU classification for osteoporotic thoracolumbar fractures classification is beneficial in determining the prognosis of patients with fractures. Furthermore, the density of the fractured vertebra, its magnitude of collapse, and the age of the individual are beneficial for predicting fracture progression.
The authors declare no financial support.
Yossi Smorgick, Ruth Pelleg‐Kallevag, Yoram Anekstein, Sergey Goldstein, and Hila May, declare that they have no conflict of interest.
We thank Mrs. Ariana Dan for editing this article.
Smorgick, Y. , Pelleg‐Kallevag, R. , Lindner, D. , Anekstein, Y. , Goldstein, S. , & May, H. (2025). Vertebral body density role in determining vertebral osteoporotic fracture type and its progression. Clinical Anatomy, 38(1), 97–104. 10.1002/ca.24219
Yossi Smorgick, Email: ysmorgick@gmail.com.
Hila May, Email: mayhila@tauex.tau.ac.il.
The datasets generated during and analyzed during the current study are available from the corresponding authors.
The datasets generated during and analyzed during the current study are available from the corresponding authors.