Localized damage in vertebral bone is most detrimental in regions of high strain energy density

An additively manufactured model for preclinical testing of cervical devices

Purpose

Composite models have become commonplace for the assessment of fixation and stability of total joint replacements; however, there are no comparable models for the cervical spine to evaluate fixation. The goal of this study was to create the framework for a tunable non-homogeneous model of cervical vertebral body by identifying the relationships between strength, in-fill density, and lattice structure and creating a final architectural framework for specific strengths to be applied to the model.

Methods

The range of material properties for cervical spine were identified from literature. Using additive manufacturing software, rectangular prints with three lattice structures, gyroid, triangle, zig-zag, and a range of in-fill densities were 3D-printed. The compressive and shear strengths for all combinations were calculated in the axial and coronal planes. Eleven unique vertebral regions were selected to represent the distribution of density. Each bone density was converted to strength and subsequently correlated to the lattice structure and in-fill density with the desired material properties. Finally, a complete cervical vertebra model was 3D-printed to ensure sufficient print quality.

Results

Materials testing identified a relationship between in-fill densities and strength for all lattice structures. The axial compressive strength of the gyroid specimens ranged from 1.5 MPa at 10% infill to 31.3 MPa at 100% infill and the triangle structure ranged from 2.7 MPa at 10% infill to 58.4 MPa at 100% infill. Based on these results, a cervical vertebra model was created utilizing cervical cancellous strength values and the corresponding in-fill density and lattice structure combination. This model was then printed with 11 different in-fill densities ranging from 33% gyroid to 84% triangle to ensure successful integration of the non-homogeneous in-fill densities and lattice structures.

Conclusions

The findings from this study introduced a framework for using additive manufacturing to create a tunable, customizable biomimetic model of a cervical vertebra.

Differences in Trabecular Microarchitecture and Simplified Boundary Conditions Limit the Accuracy of Quantitative Computed Tomography-Based Finite Element Models of Vertebral Failure

Vertebral fractures are common in the elderly, but efforts to reduce their incidence have been hampered by incomplete understanding of the failure processes that are involved. This study's goal was to elucidate failure processes in the lumbar vertebra and to assess the accuracy of quantitative computed tomography (QCT)-based finite element (FE) simulations of these processes. Following QCT scanning, spine segments (n = 27) consisting of L1 with adjacent intervertebral disks and neighboring endplates of T12 and L2 were compressed axially in a stepwise manner. A microcomputed tomography scan was performed at each loading step. The resulting time-lapse series of images was analyzed using digital volume correlation (DVC) to quantify deformations throughout the vertebral body. While some diversity among vertebrae was observed on how these deformations progressed, common features were large strains that developed progressively in the superior third and, concomitantly, in the midtransverse plane, in a manner that was associated with spatial variations in microstructural parameters such as connectivity density. Results of FE simulations corresponded qualitatively to the measured failure patterns when boundary conditions were derived from DVC displacements at the endplate. However, quantitative correspondence was often poor, particularly when boundary conditions were simplified to uniform compressive loading. These findings suggest that variations in trabecular microstructure are one cause of the differences in failure patterns among vertebrae and that both the lack of incorporation of these variations into QCT-based FE models and the oversimplification of boundary conditions limit the accuracy of these models in simulating vertebral failure.

Dependence of anisotropy of human lumbar vertebral trabecular bone on quantitative computed tomography-based apparent density

Most studies investigating human lumbar vertebral trabecular bone (HVTB) mechanical property-density relationships have presented results for the superior-inferior (SI), or "on-axis" direction. Equivalent, directly measured data from mechanical testing in the transverse (TR) direction are sparse and quantitative computed tomography (QCT) density-dependent variations in the anisotropy ratio of HVTB have not been adequately studied. The current study aimed to investigate the dependence of HVTB mechanical anisotropy ratio on QCT density by quantifying the empirical relationships between QCT-based apparent density of HVTB and its apparent compressive mechanical properties--elastic modulus (E(app)), yield strength (σ(y)), and yield strain (ε(y))--in the SI and TR directions for future clinical QCT-based continuum finite element modeling of HVTB. A total of 51 cylindrical cores (33 axial and 18 transverse) were extracted from four L1 human lumbar cadaveric vertebrae. Intact vertebrae were scanned in a clinical resolution computed tomography (CT) scanner prior to specimen extraction to obtain QCT density, ρ(CT). Additionally, physically measured apparent density, computed as ash weight over wet, bulk volume, ρ(app), showed significant correlation with ρ(CT) [ρ(CT) = 1.0568 × ρ(app), r = 0.86]. Specimens were compression tested at room temperature using the Zetos bone loading and bioreactor system. Apparent elastic modulus (E(app)) and yield strength (σ(y)) were linearly related to the ρ(CT) in the axial direction [E(SI) = 1493.8 × (ρ(CT)), r = 0.77, p < 0.01; σ(Y,SI) = 6.9 × (ρ(CT)) − 0.13, r = 0.76, p < 0.01] while a power-law relation provided the best fit in the transverse direction [E(TR) = 3349.1 × (ρ(CT))(1.94), r = 0.89, p < 0.01; σ(Y,TR) = 18.81 × (ρ(CT))(1.83), r = 0.83, p < 0.01]. No significant correlation was found between ε(y) and ρ(CT) in either direction. E(app) and σ(y) in the axial direction were larger compared to the transverse direction by a factor of 3.2 and 2.3, respectively, on average. Furthermore, the degree of anisotropy decreased with increasing density. Comparatively, ε(y) exhibited only a mild, but statistically significant anisotropy: transverse strains were larger than those in the axial direction by 30%, on average. Ability to map apparent mechanical properties in the transverse direction, in addition to the axial direction, from CT-based densitometric measures allows incorporation of transverse properties in finite element models based on clinical CT data, partially offsetting the inability of continuum models to accurately represent trabecular architectural variations.

Locally measured microstructural parameters are better associated with vertebral strength than whole bone density

Whole vertebrae areal and volumetric bone mineral density (BMD) measurements are not ideal predictors of vertebral fractures. We introduce a technique which enables quantification of bone microstructural parameters at precisely defined anatomical locations. Results show that local assessment of bone volume fraction at the optimal location can substantially improve the prediction of vertebral strength.

INTRODUCTION

Whole vertebrae areal and volumetric BMD measurements are not ideal predictors of vertebral osteoporotic fractures. Recent studies have shown that sampling bone microstructural parameters in smaller regions may permit better predictions. In such studies, however, the sampling location is described only in general anatomical terms. Here, we introduce a technique that enables the quantification of bone volume fraction and microstructural parameters at precisely defined anatomical locations. Specific goals of this study were to investigate at what anatomical location within the vertebrae local bone volume fraction best predicts vertebral-body strength, whether this prediction can be improved by adding microstructural parameters and to explore if this approach could better predict vertebral-body strength than whole bone volume fraction and finite element (FE) analyses.

METHODS

Eighteen T12 vertebrae were scanned in a micro-computed tomography (CT) system and FE meshes were made using a mesh-morphing tool. For each element, bone microstructural parameters were measured and correlated with vertebral compressive strength as measured experimentally. Whole bone volume fraction and FE-predicted vertebral strength were also compared to the experimental measurements.

RESULTS

A significant association between local bone volume fraction measured at a specific central region and vertebral-body strength was found that could explain up to 90% of the variation. When including all microstructural parameters in the regression, the predictive value of local measurements could be increased to 98%. Whole bone volume fraction could explain only 64% and FE analyses 76% of the variation in bone strength.

CONCLUSIONS

A local assessment of volume fraction at the optimal location can substantially improve the prediction of bone strength. Local assessment of other microstructural parameters may further improve this prediction but is not clinically feasible using current technology.

Trabecular bone structure in lumbosacral transitional vertebrae: distribution and densities across sagittal vertebral body segments

BACKGROUND CONTEXT

Lumbosacral transitional vertebrae (LSTV) are associated with altered articular morphology at the L5-S1 junction. Studies related to lumbo-sacral trabecular architecture in LSTV are few. Altered lumbosacral load bearing at these anomalous junctions possibly results in changes in the number, density, and trajectory of the trabecular bone in transitional lumbosacral vertebral bodies.

PURPOSE

To investigate the pattern, distribution, and density of trabecular bone in the terminal lumbar vertebrae and the first sacral segments in LSTV-affected spines. Measurements were compared with those obtained from normal lumbosacral specimens.

STUDY DESIGN

Observational and descriptive human cadaveric study of vertebral trabecular architecture.

METHODS

Blocks of tissues were obtained from normal (n=20) and LSTV cadaveric specimens (n=16) by sectioning vertically through the fifth lumbar and the first sacral vertebra on either side of the midsagittal plane. Photographs of the cut surfaces were computationally enlarged and mapped for vertical and transverse trabecular numbers and surface areas using the software Image J. All parameters including the trabecular density were computed for anterior, middle, and posterior segments of each of the vertebral elements.

RESULTS

The anterior and the posterior segments showed greater number of trabeculae across all LSTV subtypes in both the terminal lumbar and first sacral vertebrae in comparison with the middle segment. L5 exhibited greater number of vertical trabeculae, whereas the first sacral segments demonstrated greater number and densities of transverse trabeculae. Transition-associated vertebrae showed overall reduced number of the lumbar trabeculae but relatively compact sacral posterior segments with greater number of horizontal trabeculae.

CONCLUSIONS

Findings suggest that some of these variations have overall reduced number of trabeculae across lumbo-sacral vertebrae in LSTV. Screw placements and subsequent pullouts in LSTV may be reviewed in light of different trabecular patterns as reported in this study.

Vertebral fragility and structural redundancy

The mechanisms of age-related vertebral fragility remain unclear, but may be related to the degree of "structural redundancy" of the vertebra; ie, its ability to safely redistribute stress internally after local trabecular failure from an isolated mechanical overload. To better understand this issue, we performed biomechanical testing and nonlinear micro-CT-based finite element analysis on 12 elderly human thoracic ninth vertebral bodies (age 76.9 ± 10.8 years). After experimentally overloading the vertebrae to measure strength, we used nonlinear finite element analysis to estimate the amount of failed tissue and understand the failure mechanisms. We found that the amount of failed tissue per unit bone mass decreased with decreasing bone volume fraction (r(2)  = 0.66, p < 0.01). Thus, for the weak vertebrae with low bone volume fraction, overall failure of the vertebra occurred after failure of just a tiny proportion of the bone tissue (<5%). This small proportion of failed tissue had two sources: the existence of fewer vertically oriented load paths to which load could be redistributed from failed trabeculae; and the vulnerability of the trabeculae in these few load paths to undergo bending-type failure mechanisms, which further weaken the bone. Taken together, these characteristics suggest that diminished structural redundancy may be an important aspect of age-related vertebral fragility: vertebrae with low bone volume fraction are highly susceptible to collapse because so few trabeculae are available for load redistribution if the external loads cause any trabeculae to fail.

Micromechanics of the human vertebral body for forward flexion

To provide mechanistic insight into the etiology of osteoporotic wedge fractures, we investigated the spatial distribution of tissue at the highest risk of initial failure within the human vertebral body for both forward flexion and uniform compression loading conditions. Micro-CT-based linear elastic finite element analysis was used to virtually load 22 human T9 vertebral bodies in either 5° of forward flexion or uniform compression; we also ran analyses replacing the simulated compliant disc (E=8 MPa) with stiff polymethylmethacrylate (PMMA, E=2500 MPa). As expected, we found that, compared to uniform compression, forward flexion increased the overall endplate axial load on the anterior half of the vertebra and shifted the spatial distribution of high-risk tissue within the vertebra towards the anterior aspect of the vertebral body. However, despite that shift, the high-risk tissue remained primarily within the central regions of the trabecular bone and endplates, and forward flexion only slightly altered the ratio of cortical-to-trabecular load sharing at the mid-vertebral level (mean±SD for n=22: 41.3±7.4% compression; 44.1±8.2% forward flexion). When the compliant disc was replaced with PMMA, the anterior shift of high-risk tissue was much more severe. We conclude that, for a compliant disc, a moderate degree of forward flexion does not appreciably alter the spatial distribution of stress within the vertebral body.

Time-dependent compressive deformation of the ageing spine: relevance to spinal stenosis

STUDY DESIGN

Mechanical testing of cadaveric spines.

OBJECTIVE

To test the hypothesis that, in the ageing spine, vertebrae deform more than discs, and contribute to time-dependent creep.

SUMMARY OF BACKGROUND DATA

Intervertebral discs and vertebrae deform under load, narrowing the intervertebral foramen and increasing the risk of nerve root entrapment. Little is known about compressive deformations when elderly spines are subjected to sustained physiologic loading.

METHODS

A total of 117 thoracolumbar motion segments, aged 19 to 96 yrs (mean, 69), were subjected to 1kN compressive loading for 0.5, 1, or 2 hours. Deformations during the first 7 seconds were designated "elastic" and subsequent deformations as "creep". A 3-parameter model was fitted to experimental data in order to characterize their viscous modulus E1, elastic modulus E2 (initial stiffness), and viscosity eta (resistance to fluid flow). Intradiscal pressure (IDP) was measured using a miniature needle-mounted transducer. In 17 specimens loaded for 0.5 hours, an optical MacReflex system measured compressive deformations separately in the disc and each vertebral body.

RESULTS

On average, the disc contributed 28% of the spine's elastic deformation, 51% of the creep deformation, and 38% of total deformation. Elastic, creep, and total deformations of 84 motion segments in 2-hour tests averaged 0.87, 1.37, and 2.24 mm respectively. Measured deformations were predicted accurately by the model (average r2 = 0.97), but E1, E2, and eta depended on the duration of loading. E1 and eta decreased with advancing age and disc degeneration, in proportion to falling IDP (P < 0.001). Total compressive deformation increased with age, but rarely exceeded 3 mm.

CONCLUSION

When the ageing spine is compressed, vertebral bodies show greater elastic deformations than intervertebral discs, and creep by a similar amount. Responses to axial compression depend largely on IDP, but deformations appear to be limited by impaction of adjacent neural arches. Total compressive deformations are sufficient to cause foraminal stenosis in some individuals.

Predicting trabecular bone microdamage initiation and accumulation using a non-linear perfect damage model

Studies evaluating the mechanical behavior of the trabecular microstructure play an important role in our understanding of pathologies such as osteoporosis, and in increasing our understanding of bone fracture and bone adaptation. Understanding of such behavior in bone is important for predicting and providing early treatment of fractures. The objective of this study is to present a numerical model for studying the initiation and accumulation of trabecular bone microdamage in both the pre- and post-yield regions. A sub-region of human vertebral trabecular bone was analyzed using a uniformly loaded anatomically accurate microstructural three-dimensional finite element model. The evolution of trabecular bone microdamage was governed using a non-linear, modulus reduction, perfect damage approach derived from a generalized plasticity stress-strain law. The model introduced in this paper establishes a history of microdamage evolution in both the pre- and post-yield regions.

The Effect of Regional Variations of the Trabecular Bone Properties on the Compressive Strength of Human Vertebral Bodies

Cancellous centrum is a major component of the vertebral body and significantly contributes to its structural strength and fracture risk. We hypothesized that the variability of cancellous bone properties in the centrum is associated with vertebral strength. Microcomputed tomography (micro-CT)-based gray level density (GLD), bone volume fraction (BV/TV), and finite element modulus (E) were examined for different regions of the trabecular centrum and correlated with vertebral body strength determined experimentally. Two sets of images in the cancellous centrum were digitally prepared from micro-CT images of eight human vertebral bodies (T10-L5). One set included a cubic volume (1 per vertebral centrum, n = 8) in which the largest amount of cancellous material from the centrum was included but all the shell materials were excluded. The other set included cylindrical volumes (6 per vertebral centrum, n = 48) from the anterior (4 regions: front, center, left, and right of the midline of vertebra) and the posterior (2 regions: left and right) regions of the centrum. Significant positive correlations of vertebral strength with GLD (r (2) = 0.57, p = 0.03) and E (r (2) = 0.63, p = 0.02) of the whole centrum and with GLD (r (2) = 0.65, p = 0.02), BV/TV (r (2) = 0.72, p = 0.01) and E (r (2) = 0.85, p = 0.001) of the central region of the vertebral centrum were found. Vertebral strength decreased with increasing coefficient of variation of GLD, BV/TV, and E calculated from subregions of the vertebral centrum. The values of GLD, BV/TV, and E in centrum were significantly smaller for the anterior region than for the posterior region. Overall, these findings supported the significant role of regional variability of centrum properties in determining the whole vertebral strength.

The influence of material property and morphological parameters on specimen-specific finite element models of porcine vertebral bodies

The use of finite element (FE) methods in spinal research is increasing, but there is only limited information available on the influence of different input parameters on the model predictions. The aim of this study was to investigate the role of these parameters in FE models of the vertebral body. Experimental tests were undertaken on porcine lumbar vertebral bodies and scans of the specimens were used to create specimen-specific FE models. Three models were created for each specimen with combinations of generic and specimen-specific parameters. Stiffness and strength predictions were also made directly from the specimen trabecular bone volume fraction (BVF) and cross-sectional area (CSA). The agreement between the experimental results and the FE models with generic morphology was poorer (concordance coefficients = 0.058, 0.125 for stiffness, strength) than those made from the BVF and CSA (concordance coefficients = 0.638, 0.609). The greatest levels of agreement were found with the morphologically specific models including element-specific material properties (concordance coefficients = 0.881, 0.752). This indicates that highly specific models, both in terms of morphology and bone quality, are necessary if the FE tool is to be used effectively for spinal research and clinical practice.

Antero–postero differences in cortical thickness and cortical porosity of T12 to L5 vertebral bodies

OBJECTIVE

This study will investigate interrelationships between the cortical shell and cancellous bone trabecular thickness, in vertebral bodies.

METHODS

One hundred and sixty vertebral bodies from T12 to L5 were obtained at autopsy. The average age of the cohort was 59.3+/-22.1 years (range = 20-94 years). Cortical thickness, cortical porosity and trabecular thickness from the adjacent cancellous bone were measured.

RESULTS

At the mid-vertebral body anterior cortical thickness was significantly greater than posterior cortical thickness (524 +/- 352 vs. 370 +/- 283 microm, respectively, P < 0.0001) and mid-anterior cortical porosity was significantly less than mid-posterior cortical porosity (24 +/- 14% vs. 32 +/- 16%, respectively, P < 0.0001). There were no anterior/posterior differences in trabecular thickness of the cancellous bone adjacent to the cortical walls.

CONCLUSION

This study provides a novel perspective of T12 to L5 vertebral body bone, where measurement of cortical thickness and cortical porosity in a cohort of skeletally normal individuals revealed structural differences between load bearing anterior and posterior cortical walls. The data suggest that modulators of change to vertebral body bone may affect the cortical and trabecular bone differently. The relationships between cortical and cancellous bone suggest that the middle sectors of the vertebral body play a critical role in load bearing.

Regional variations in microstructural properties of vertebral trabeculae with structural groups

STUDY DESIGN

Micro-computed tomography (CT) scanning to investigate three-dimensional microstructural properties of L4 vertebral bodies.

OBJECTIVE

To identify the regional variations in the three-dimensional microstructural properties of vertebral cancellous bones with respect to structural types for the prediction of related regional fracture risks.

SUMMARY OF BACKGROUND DATA

The literature contains no reports on regional variations in morphologic properties of vertebral trabeculae with microstructural types, which may shed light on the patterns of osteoporotic fractures.

METHODS

Ninety cubic cancellous specimens were obtained from 6 normal L4 vertebral bodies of 6 male donors 62 to 70 years of age and were scanned using a high-resolution micro-CT system. These specimens were further divided into two groups according to the average structure model index (SMI) of the 15 trabecular specimens in each vertebral body. Adjustment for age differences was done for the microstructural parameters, i.e.-, bone volume fraction, trabecular number, trabecular thickness, structure model index, degree of architectural anisotropy, and connectivity density, to allow investigation on the regional variations in different transverse layers and vertical columns independent of age.

RESULTS

Trabecular specimens with lower mass were liable to form high-SMI group and the differences in all parameters reached significance level either between columns or between layers from two groups.

CONCLUSIONS

The anterior column in the high-SMI group is more susceptible to vertebral body wedge fracture; and in the low-SMI group, off-axis bone damage is most harmful to the central column of vertebral trabeculae. The data obtained may help to identify the most critical locations of fracture risks at an early stage and provide a microstructural basis for the repair and clinical treatment of vertebral fractures.

Automatic generation of accurate subject-specific bone finite element models to be used in clinical studies

Most of the finite element models of bones used in orthopaedic biomechanics research are based on generic anatomies. However, in many cases it would be useful to generate from CT data a separate finite element model for each subject of a study group. In a recent study a hexahedral mesh generator based on a grid projection algorithm was found very effective in terms of accuracy and automation. However, so far the use of this method has been documented only on data collected in vitro and only for long bones. The present study was aimed at verifying if this method represents a procedure for the generation of finite element models of human bones from data collected in vivo, robust, accurate, automatic and general enough to be used in clinical studies. Robustness, automation and numerical accuracy of the proposed method were assessed on five femoral CT data sets of patients affected by various pathologies. The generality of the method was verified by processing a femur, an ileum, a phalanx, a proximal femur reconstruction, and the micro-CT of a small sample of spongy bone. The method was found robust enough to cope with the variability of the five femurs, producing meshes with a numerical accuracy and a computational weight comparable to those found in vitro. Even when the method was used to process the other bones the levels of mesh conditioning remained within acceptable limits. Thus, it may be concluded that the method presents a generality sufficient to cope with almost any orthopaedic application.

Finite element modeling of trabecular bone damage

This paper presents a finite element-based, computational model for analysis of structural damage to trabecular bone tissues. A modulus reduction method was formulated from elasto-plasticity theory, and was used to account for site-specific trabecular bone tissue damage. Trabecular bone tissue damage is illustrated using a large-scale, anatomically accurate, two-dimensional, microstructural finite element model of a human thoracic vertebral body. Four models with varying specifications for damage accumulation were subjected to compressive loading and unloading cycles. The numerical results and experimental validation demonstrated that the modulus reduction method reproduced the non-linear mechanical behaviour of vertebal trabecular bone. The iterative computational approach presented provides a methodology to study trabecular bone damage, and should provide researchers with a computational approach to study bone fracture and repair and to predict vertebral fragility.

Quantitative computed tomography estimates of the mechanical properties of human vertebral trabecular bone

The objective of this study was to report our quantitative computed tomography (QCT) density-mechanical property regressions for trabecular bone for use in biomechanical modelling of the human spine. Cylindrical specimens of human vertebral trabecular bone (from T10 to L4) were cored from 32 cadavers (mean +/- SD age = 70.1 +/- 16.8; 13 females, 19 males) and scanned using QCT. Mechanical tests were conducted using a protocol that minimized end-artifacts over the apparent density range tested (0.09-0.38 g/cm3). To account for the presence of multiple specimens per donor in this data set, donor was treated as a random effect in the regression model. Mean modulus (319 +/- 189 MPa) was higher and mean yield strain (0.78 +/- 0.06%) was lower than typical values reported previously due to minimization of the end-artifact errors. QCT density showed a strong positive correlation with modulus (n = 76) and yield stress (r2 = 0.90-0.95, n = 53, p < 0.001). There was a weak positive linear correlation with yield strain (r2 = 0.58, n = 53, p = 0.07). Prediction errors, incurred when estimating modulus or strength for specimens from a new donor, were 30-36% of the mean values of these properties. Direct QCT density-mechanical property regressions gave more precise predictions of mechanical properties than if physically measured wet apparent density was used as an intermediate variable to predict mechanical properties from QCT density. Use of these QCT density-mechanical property regressions should improve the fidelity of QCT-based biomechanical models of the human spine for whole bone and bone-implant analyses.

Intervertebral disc disorganization is related to trabecular bone architecture in the lumbar spine

Cancellous bone morphometry was investigated in the sagittal plane of lumbar vertebrae using histoquantitation. The aim of this study was to identify variations in cancellous bone architecture at increasing states of intervertebral disc (IVD) disorganization after age adjustment and to investigate regional variations within the whole vertebral body. Measurements were taken of the ratio of bone volume (BV) to total volume (TV), trabecular thickness (Tb.Th), trabecular separation (Tb.Sp), and trabecular number (Tb.N). Lumbar spines (T12-L5) of 19 men and 8 women were removed at autopsy from an adult sample with no clinical history of bone-related disease or histologically identifiable bone disease. It was found that degeneration of the IVD becomes more common with increasing age. After age-adjustment, significant increases in the proportion of BV/TV were observed in the presence of advancing IVD disorganization. Significant architectural changes were observed in the anterior regions of the vertebral body with increases in Tb.Th and Tb.N and decreases in Tb.Sp. Minimal alterations were found at posterior regions. Bone loss was observed in central regions (most distant from the cortex) as IVD disorganization increased through reduction in both Tb.N and Tb.Th. The BV/TV increase in anterior areas of the centrum may be a response to a redistribution of load to the vertebral body periphery as a result of IVD disorganization. It appears that trabecular morphology is related to the condition of the associated IVD, rather than being the sole consequence of a loss of BV/TV with age. This relationship could influence the occurrence of vertebral body crush fracture.