Relationship between axial and bending behaviors of the human thoracolumbar vertebra

Two-year Skeletal Effects of Sleeve Gastrectomy in Adolescents with Obesity Assessed with Quantitative CT and MR Spectroscopy

Background Sleeve gastrectomy (SG) is effective in the treatment of cardiometabolic complications of obesity but is associated with bone loss. Purpose To determine the long-term effects of SG on vertebral bone strength, density, and bone marrow adipose tissue (BMAT) in adolescents and young adults with obesity. Materials and Methods This 2-year prospective nonrandomized longitudinal study enrolled adolescents and young adults with obesity who underwent either SG (SG group) or dietary and exercise counseling without surgery (control group) at an academic medical center from 2015 to 2020. Participants underwent quantitative CT of the lumbar spine (L1 and L2 levels) to assess bone density and strength, proton MR spectroscopy to assess BMAT (L1 and L2 levels), and MRI of the abdomen and thigh to assess body composition. Student t and Wilcoxon signed-rank tests were used to compare 24-month changes between and within groups. Regression analysis was performed to evaluate associations between body composition, vertebral bone density, strength, and BMAT. Results A total of 25 participants underwent SG (mean age, 18 years ± 2 [SD], 20 female), and 29 underwent dietary and exercise counseling without surgery (mean age, 18 years ± 3, 21 female). Body mass index (BMI) decreased by a mean of 11.9 kg/m2 ± 5.21 [SD] after 24 months in the SG group (P < .001), while it increased in the control group (mean increase, 1.49 kg/m2 ± 3.10; P = .02). Mean bone strength of the lumbar spine decreased after surgery compared with that in control subjects (mean decrease, -728 N ± 691 vs -7.24 N ± 775; P < .001). BMAT of the lumbar spine increased after SG (mean lipid-to-water ratio increase, 0.10 ± 0.13; P = .001). Changes in vertebral density and strength correlated positively with changes in BMI and body composition (R = 0.34 to R = 0.65, P = .02 to P < .001) and inversely with vertebral BMAT (R = -0.33 to R = -0.47, P = .03 to P = .001). Conclusion SG in adolescents and young adults reduced vertebral bone strength and density and increased BMAT compared with those in control participants. Clinical trial registration no. NCT02557438 © RSNA, 2023 See also the editorial by Link and Schafer in this issue.

Mechanical Properties and Bioactivity of Poly(Lactic Acid) Composites Containing Poly(Glycolic Acid) Fiber and Hydroxyapatite Particles

To enhance the mechanical strength and bioactivity of poly(lactic acid) (PLA) to the level that can be used as a material for spinal implants, poly(glycolic acid) (PGA) fibers and hydroxyapatite (HA) were introduced as fillers to PLA composites. To improve the poor interface between HA and PLA, HA was grafted by PLA to form HA-g-PLA through coupling reactions, and mixed with PLA. The size of the HA particles in the PLA matrix was observed to be reduced from several micrometers to sub-micrometer by grafting PLA onto HA. The tensile and flexural strength of PLA/HA-g-PLA composites were increased compared with those of PLA/HA, apparently due to the better dispersion of HA and stronger interfacial adhesion between the HA and PLA matrix. We also examined the effects of the length and frequency of grafted PLA chains on the tensile strength of the composites. By the addition of unidirectionally aligned PGA fibers, the flexural strength of the composites was greatly improved to a level comparable with human compact bone. In the bioactivity tests, the growth of apatite on the surface was fastest and most uniform in the PLA/PGA fiber/HA-g-PLA composite.

Synthesis and Modification of Hydroxyapatite Nanofiber for Poly(Lactic Acid) Composites with Enhanced Mechanical Strength and Bioactivity

To enhance the bioactivity of poly(lactic acid) (PLA), a potential bone repair material, without the lowering of mechanical strength, hydroxyapatite (HA) was introduced in the form of nanofibers as the filler for application in spinal implant materials. HA nanofibers (HANF) with aspect ratio as high as ~100 were synthesized by controlling the starting pH of the reaction. While the tensile and flexural strength of PLA/HANF composites were enhanced compared with those of PLA resin, and were higher for the composites with HANF of higher aspect ratio. To further strengthen the composites, HANF was grafted with PLA chain to form HANF-g-PLA, which could improve the interface between the HANF and matrix PLA. PLA/HANF-g-PLA composites showed even higher tensile and flexural strength than PLA/HANF composites, apparently due to the better dispersion and interfacial adhesion. The composite containing 10 wt% HANF-g-PLA showed the flexural strength of 124 MPa, which was 25% higher than that of PLA resin. In the bioactivity test using a simulated body fluid solution, the rate and uniformity of the apatite growth were observed to be higher for the composites with HANF, and were even higher for those with HANF-g-PLA. This study suggested the possibility of using the PLA/HANF-g-PLA composite in the field of spinal implant materials.

Temporal Trends in Vertebral Dimensions - a case study from Finland

Vertebral fractures and other back problems represent a major, increasing worldwide health problem. This has increased the need to better understand the reasons behind this phenomenon. In addition to a reduction in bone mineral density and overall size of the vertebral body, research has indicated a possible association between the shape of the endplate and spinal disorders. As one previous study has shown changes in vertebral body dimensions between contemporary people and their medieval counterparts, we wanted to examine the potential temporal trends in vertebral size and dimensions in Finnish samples of archaeological and contemporary individuals. To conduct this study, we utilized three archaeological populations from the 16^th–19^th century and clinical materials from two population-based Finnish birth cohorts. As the average height of people has increased greatly since the first time period, we also height-adjusted the dimensions to provide a clearer picture of the dimensional changes that have occurred in the later temporal group. Our results were in agreement with those of the earlier study. The archaeological samples had a larger vertebral size than the contemporary population when height was adjusted for. Vertebral mediolateral width in particular had decreased, and the shape of the vertebral body had changed.

Assessment of vertebral wedge strength using cancellous textural properties derived from digital tomosynthesis and density properties from dual energy X-ray absorptiometry and high resolution computed tomography

The purpose of this study was to examine the potential of digital tomosynthesis (DTS) derived cancellous bone textural measures to predict vertebral strength under conditions simulating a wedge fracture. 40 vertebral bodies (T6, T8, T11, and L3 levels) from 5 male and 5 female cadaveric donors were utilized. The specimens were scanned using dual energy X-ray absorptiometry (DXA) and high resolution computed tomography (HRCT) to obtain measures of bone mineral density (BMD) and content (BMC), and DTS to obtain measures of bone texture. Using a custom loading apparatus designed to deliver a nonuniform displacement resulting in a wedge deformity similar to those observed clinically, the specimens were loaded to fracture and their fracture strength was recorded. Mixed model regressions were used to determine the associations between wedge strength and DTS derived textural variables, alone and in the presence of BMD or BMC information. DTS derived fractal, lacunarity and mean intercept length variables correlated with wedge strength, and individually explained up to 53% variability. DTS derived textural variables, notably fractal dimension and lacunarity, contributed to multiple regression models of wedge strength independently from BMC and BMD. The model from a scan orientation transverse to the spine axis and in the anterior-posterior view resulted in highest explanatory capability (R²_adj = 0.91), with a scan orientation parallel to the spine axis and in the lateral view offering an alternative (R²_adj = 0.88). In conclusion, DTS can be used to examine cancellous texture relevant to vertebral wedge strength, and potentially complement BMD in assessment of vertebral fracture risk.

Quantitative, 3D Visualization of the Initiation and Progression of Vertebral Fractures Under Compression and Anterior Flexion

The biomechanical mechanisms leading to vertebral fractures are not well understood. Clinical and laboratory evidence suggests that the vertebral endplate plays a key role in failure of the vertebra as a whole, but how this role differs for different types of vertebral loading is not known. Mechanical testing of human thoracic spine segments, in conjunction with time-lapsed micro-computed tomography, enabled quantitative assessment of deformations occurring throughout the entire vertebral body under axial compression combined with anterior flexion ("combined loading") and under axial compression only ("compression loading"). The resulting deformation maps indicated that endplate deflection was a principal feature of vertebral failure for both loading modes. Specifically, the onset of endplate deflection was temporally coincident with a pronounced drop in the vertebra's ability to support loads. The location of endplate deflection, and also vertebral strength, were associated with the porosity of the endplate and the microstructure of the underlying trabecular bone. However, the location of endplate deflection and the involvement of the cortex differed between the two types of loading. Under the combined loading, deflection initiated, and remained the largest, at the anterior central endplate or the anterior ring apophysis, depending in part on health of the adjacent intervertebral disc. This deflection was accompanied by outward bulging of the anterior cortex. In contrast, the location of endplate deflection was more varied in compression loading. For both loading types, the earliest progression to a mild fracture according to a quantitative morphometric criterion occurred only after much of the failure process had occurred. The outcomes of this work indicate that for two physiological loading modes, the vertebral endplate and underlying trabecular bone are critically involved in vertebral fracture. These outcomes provide a strong biomechanical rationale for clinical methods, such as algorithm-based qualitative (ABQ) assessment, that diagnose vertebral fracture on the basis of endplate depression. © 2015 American Society for Bone and Mineral Research.

Vertebral Strength and Estimated Fracture Risk Across the BMI Spectrum in Women

Somewhat paradoxically, fracture risk, which depends on applied loads and bone strength, is elevated in both anorexia nervosa and obesity at certain skeletal sites. Factor-of-risk (Φ), the ratio of applied load to bone strength, is a biomechanically based method to estimate fracture risk; theoretically, higher Φ reflects increased fracture risk. We estimated vertebral strength (linear combination of integral volumetric bone mineral density [Int.vBMD] and cross-sectional area from quantitative computed tomography [QCT]), vertebral compressive loads, and Φ at L4 in 176 women (65 anorexia nervosa, 45 lean controls, and 66 obese). Using biomechanical models, applied loads were estimated for: 1) standing; 2) arms flexed 90°, holding 5 kg in each hand (holding); 3) 45° trunk flexion, 5 kg in each hand (lifting); 4) 20° trunk right lateral bend, 10 kg in right hand (bending). We also investigated associations of Int.vBMD and vertebral strength with lean mass (from dual-energy X-ray absorptiometry [DXA]) and visceral adipose tissue (VAT, from QCT). Women with anorexia nervosa had lower, whereas obese women had similar, Int.vBMD and estimated vertebral strength compared with controls. Vertebral loads were highest in obesity and lowest in anorexia nervosa for standing, holding, and lifting (p < 0.0001) but were highest in anorexia nervosa for bending (p < 0.02). Obese women had highest Φ for standing and lifting, whereas women with anorexia nervosa had highest Φ for bending (p < 0.0001). Obese and anorexia nervosa subjects had higher Φ for holding than controls (p < 0.03). Int.vBMD and estimated vertebral strength were associated positively with lean mass (R = 0.28 to 0.45, p ≤ 0.0001) in all groups combined and negatively with VAT (R = -[0.36 to 0.38], p < 0.003) within the obese group. Therefore, women with anorexia nervosa had higher estimated vertebral fracture risk (Φ) for holding and bending because of inferior vertebral strength. Despite similar vertebral strength as controls, obese women had higher vertebral fracture risk for standing, holding, and lifting because of higher applied loads from higher body weight. Examining the load-to-strength ratio helps explain increased fracture risk in both low-weight and obese women.

Effect of the metastatic defect on the structural response and failure process of human vertebrae: an experimental study

BACKGROUND

Pathologic vertebral fractures are associated with intractable pain, loss of function and high morbidity in patients with metastatic spine disease. However, the failure mechanisms of vertebrae with lytic defects and the failed vertebrae's ability to retain load carrying capacity remain unclear.

METHODS

Eighteen human thoracic and lumbar vertebrae with simulated uncontained bone defects were tested under compression-bending loads to failure. Failure was defined as 50% reduction in vertebral body height. The vertebrae were allowed to recover under load and re-tested to failure using the initial criteria. Repeated measure ANOVA was used to test for changes in strength and stiffness parameters.

FINDINGS

Vertebral failure occurred via buckling and fracture of the cortex around the defect, followed by collapse of the defect region. Compared to the intact vertebrae, the failed vertebrae exhibited a significant loss in compressive strength (59%, p<0.001), stiffness (53%, p<0.05) and flexion (70%, p<0.01) strength. Significant reduction in anterior-posterior shear (strength (63%, p<0.01) and stiffness (67%, p<0.01)) and lateral bending strength (134%, p<0.05) were similarly recorded. In the intact vertebrae, apart from flexion strength (r(2)=0.63), both compressive and anterior-posterior shear strengths were weakly correlated with their stiffness parameters (r(2)=0.24 and r(2)=0.31). By contrast, in the failed vertebrae, these parameters were strongly correlated, (r(2)=0.91, r(2)=0.86, and r(2)=0.92, p<0.001 respectively).

INTERPRETATION

Failure of the vertebral cortex at the defect site dominated the initiation and progression of vertebral failure with the vertebrae failing via a consolidation process of the vertebral bone. Once failed, the vertebrae showed remarkable loss of load carrying capacity.

A patient-specific computer tomography-based finite element methodology to calculate the six dimensional stiffness matrix of human vertebral bodies

In most finite element (FE) studies of vertebral bodies, axial compression is the loading mode of choice to investigate structural properties, but this might not adequately reflect the various loads to which the spine is subjected during daily activities or the increased fracture risk associated with shearing or bending loads. This work aims at proposing a patient-specific computer tomography (CT)-based methodology, using the currently most advanced, clinically applicable finite element approach to perform a structural investigation of the vertebral body by calculation of its full six dimensional (6D) stiffness matrix. FE models were created from voxel images after smoothing of the peripheral voxels and extrusion of a cortical shell, with material laws describing heterogeneous, anisotropic elasticity for trabecular bone, isotropic elasticity for the cortex based on experimental data. Validated against experimental axial stiffness, these models were loaded in the six canonical modes and their 6D stiffness matrix calculated. Results show that, on average, the major vertebral rigidities correlated well or excellently with the axial rigidity but that weaker correlations were observed for the minor coupling rigidities and for the image-based density measurements. This suggests that axial rigidity is representative of the overall stiffness of the vertebral body and that finite element analysis brings more insight in vertebral fragility than densitometric approaches. Finally, this extended patient-specific FE methodology provides a more complete quantification of structural properties for clinical studies at the spine.

Failure strength of human vertebrae: prediction using bone mineral density measured by DXA and bone volume by micro-CT

Significant relationships exist between areal bone mineral density (BMD) derived from dual energy X-ray absorptiometry (DXA) and bone strength. However, the predictive validity of BMD for osteoporotic vertebral fractures remains suboptimal. The diagnostic sensitivity of DXA in the lumbar spine may be improved by assessing BMD from lateral-projection scans, as these might better approximate the objective of measuring the trabecular-rich bone in the vertebral body, compared to the commonly-used posterior-anterior (PA) projections. Nowadays, X-ray micro-computed tomography (μCT) allows non-destructive three-dimensional structural characterization of entire bone segments at high resolution. In this study, human lumbar cadaver spines were examined ex situ by DXA in lateral and PA projections, as well as by μCT, with the aims (1) to investigate the ability of bone quantity measurements obtained by DXA in the lateral projection and in the PA projection, to predict variations in bone quantity measurements obtained by μCT, and (2) to assess their respective capabilities to predict whole vertebral body strength, determined experimentally. Human cadaver spines were scanned by DXA in PA projections and lateral projections. Bone mineral content (BMC) and BMD for L2 and L3 vertebrae were determined. The L2 and L3 vertebrae were then dissected and entirely scanned by μCT. Total bone volume (BV(tot)=cortical+trabecular), trabecular bone volume (BV), and trabecular bone volume fraction (BV/TV) were calculated over the entire vertebrae. The vertebral bodies were then mechanically tested to failure in compression, to determine ultimate load. The variables BV(tot), BV, and BV/TV measured by μCT were better predicted by BMC and BMD measured by lateral-projection DXA, with higher R(2) values and smaller standard errors of the estimate (R(2)=0.65-0.90, SEE=11%-18%), compared to PA-projection DXA (R(2)=0.33-0.53, SEE=22%-34%). The best predictors of ultimate load were BV(tot) and BV assessed by μCT (R(2)=0.88 and R(2)=0.81, respectively), and BMC and BMD from lateral-projection DXA (R(2)=0.82 and R(2)=0.70, respectively). Conversely, BMC and BMD from PA-projection DXA were lower predictors of ultimate load (R(2)=0.49 and R(2)=0.37, respectively). This ex vivo study highlights greater capabilities of lateral-projection DXA to predict variations in vertebral body bone quantity as measured by μCT, and to predict vertebral strength as assessed experimentally, compared to PA-projection DXA. This provides basis for further exploring the clinical application of lateral-projection DXA analysis.

Prediction of new clinical vertebral fractures in elderly men using finite element analysis of CT scans

Vertebral strength, as estimated by finite element analysis of computed tomography (CT) scans, has not yet been compared against areal bone mineral density (BMD) by dual-energy X-ray absorptiometry (DXA) for prospectively assessing the risk of new clinical vertebral fractures. To do so, we conducted a case-cohort analysis of 306 men aged 65 years and older, which included 63 men who developed new clinically-identified vertebral fractures and 243 men who did not, all observed over an average of 6.5 years. Nonlinear finite element analysis was performed on the baseline CT scans, blinded to fracture status, to estimate L1 vertebral compressive strength and a load-to-strength ratio. Volumetric BMD by quantitative CT and areal BMD by DXA were also evaluated. We found that, for the risk of new clinical vertebral fracture, the age-adjusted hazard ratio per standard deviation change for areal BMD (3.2; 95% confidence interval [CI], 2.0-5.2) was significantly lower (p < 0.005) than for strength (7.2; 95% CI, 3.6-14.1), numerically lower than for volumetric BMD (5.7; 95% CI, 3.1-10.3), and similar for the load-to-strength ratio (3.0; 95% CI, 2.1-4.3). After also adjusting for race, body mass index (BMI), clinical center, and areal BMD, all these hazard ratios remained highly statistically significant, particularly those for strength (8.5; 95% CI, 3.6-20.1) and volumetric BMD (9.4; 95% CI, 4.1-21.6). The area-under-the-curve for areal BMD (AUC = 0.76) was significantly lower than for strength (AUC = 0.83, p = 0.02), volumetric BMD (AUC = 0.82, p = 0.05), and the load-to-strength ratio (AUC = 0.82, p = 0.05). We conclude that, compared to areal BMD by DXA, vertebral compressive strength and volumetric BMD consistently improved vertebral fracture risk assessment in this cohort of elderly men.

Predicting fracture using 2D finite element modelling

A decrease in bone density at the hip or spine has been shown to increase the risk of fracture. A limitation of the bone mineral density (BMD) measurement is that it provides only a measure of a bone sample's average density when projected onto a 2D surface. Effectively, what determines bone fracture is whether an applied load exceeds ultimate strength, with both bone tissue material properties (can be approximated through bone density), and geometry playing a role. The goal of this project was to use bone geometry and BMD obtained from radiographs and DXA measurements respectively to estimate fracture risk, using a two-dimensional finite element model (FEM) of the sagittal plane of lumbar vertebrae. The Canadian Multicentre Osteoporosis Study (CaMos) data was used for this study. There were 4194 men and women over the age of 50 years, with 786 having fractures. Each subject had BMD testing and radiographs of their lumbar vertebrae. A single two dimensional FEM of the first to fourth lumbar vertebra was automatically generated for each subject. Bone tissue stiffness was assigned based on the BMD of the individual vertebrae, and adjusted for patient age. Axial compression boundary conditions were applied with a force proportional to body mass. The resulting overall strain from the applied force was found. Men and women were analyzed separately. At baseline, the sensitivity of BMD to predict fragility fractures in women and men was 3.77% and 0.86%, while the sensitivity of FEM to predict fragility fractures for women and men was 10.8% and 11.3%. The FEM ROC curve demonstrated better performance compared to BMD. The relative risk of being considered at high fracture risk using FEM at baseline, was a better predictor of 5 year incident fragility fracture risk compared to BMD.

Mechanical contributions of the cortical and trabecular compartments contribute to differences in age-related changes in vertebral body strength in men and women assessed by QCT-based finite element analysis

The biomechanical mechanisms underlying sex-specific differences in age-related vertebral fracture rates are ill defined. To gain insight into this issue, we used finite element analysis of clinical computed tomography (CT) scans of the vertebral bodies of L3 and T10 of young and old men and women to assess age- and sex-related differences in the strength of the whole vertebra, the trabecular compartment, and the peripheral compartment (the outer 2 mm of vertebral bone, including the thin cortical shell). We sought to determine whether structural and geometric changes with age differ in men and women, making women more susceptible to vertebral fractures. As expected, we found that vertebral strength decreased with age 2-fold more in women than in men. The strength of the trabecular compartment declined significantly with age for both sexes, whereas the strength of the peripheral compartment decreased with age in women but was largely maintained in men. The proportion of mechanical strength attributable to the peripheral compartment increased with age in both sexes and at both vertebral levels. Taken together, these results indicate that men and women lose vertebral bone differently with age, particularly in the peripheral (cortical) compartment. This differential bone loss explains, in part, a greater decline in bone strength in women and may contribute to the higher incidence of vertebral fractures among women than men.

Influence of vertical trabeculae on the compressive strength of the human vertebra

Vertebral strength, a key etiologic factor of osteoporotic fracture, may be affected by the relative amount of vertically oriented trabeculae. To better understand this issue, we performed experimental compression testing, high-resolution micro-computed tomography (µCT), and micro-finite-element analysis on 16 elderly human thoracic ninth (T(9)) whole vertebral bodies (ages 77.5 ± 10.1 years). Individual trabeculae segmentation of the µCT images was used to classify the trabeculae by their orientation. We found that the bone volume fraction (BV/TV) of just the vertical trabeculae accounted for substantially more of the observed variation in measured vertebral strength than did the bone volume fraction of all trabeculae (r(2) = 0.83 versus 0.59, p < .005). The bone volume fraction of the oblique or horizontal trabeculae was not associated with vertebral strength. Finite-element analysis indicated that removal of the cortical shell did not appreciably alter these trends; it also revealed that the major load paths occur through parallel columns of vertically oriented bone. Taken together, these findings suggest that variation in vertebral strength across individuals is due primarily to variations in the bone volume fraction of vertical trabeculae. The vertical tissue fraction, a new bone quality parameter that we introduced to reflect these findings, was both a significant predictor of vertebral strength alone (r(2) = 0.81) and after accounting for variations in total bone volume fraction in multiple regression (total R(2) = 0.93). We conclude that the vertical tissue fraction is a potentially powerful microarchitectural determinant of vertebral strength.

Relation of vertebral deformities to bone density, structure, and strength

Because they are not reliably discriminated by areal bone mineral density (aBMD) measurements, it is unclear whether minimal vertebral deformities represent early osteoporotic fractures. To address this, we compared 90 postmenopausal women with no deformity (controls) with 142 women with one or more semiquantitative grade 1 (mild) deformities and 51 women with any grade 2-3 (moderate/severe) deformities. aBMD was measured by dual-energy X-ray absorptiometry (DXA), lumbar spine volumetric bone mineral density (vBMD) and geometry by quantitative computed tomography (QCT), bone microstructure by high-resolution peripheral QCT at the radius (HRpQCT), and vertebral compressive strength and load-to-strength ratio by finite-element analysis (FEA) of lumbar spine QCT images. Compared with controls, women with grade 1 deformities had significantly worse values for many bone density, structure, and strength parameters, although deficits all were much worse for the women with grade 2-3 deformities. Likewise, these skeletal parameters were more strongly associated with moderate to severe than with mild deformities by age-adjusted logistic regression. Nonetheless, grade 1 vertebral deformities were significantly associated with four of the five main variable categories assessed: bone density (lumbar spine vBMD), bone geometry (vertebral apparent cortical thickness), bone strength (overall vertebral compressive strength by FEA), and load-to-strength ratio (45-degree forward bending ÷ vertebral compressive strength). Thus significantly impaired bone density, structure, and strength compared with controls indicate that many grade 1 deformities do represent early osteoporotic fractures, with corresponding implications for clinical decision making.

Variation of trabecular microarchitectural parameters in cranial, caudal and mid-vertebral regions of the ovine L3 vertebra

The lumbar vertebrae are major load-bearing structures within the spinal column. The current understanding of the microstructure of these bodies and their full role in load-bearing is incomplete. There is a need to develop our understanding of these issues to improve fracture prediction in musculoskeletal diseases such as osteoporosis. The lumbar vertebrae consist primarily of trabecular bone enclosed in a thin cortical shell, but little is known about how microstructural parameters vary within these structures, particularly in relation to the trabecular compartment. The specific aim of this study was to use micro-computed tomography to characterize the trabecular microarchitecture of the ovine L3 vertebra in cranial, mid-vertebra and caudal regions. The L3 vertebra was obtained from skeletally mature ewes (n = 18) more than 4 years old. Three-dimensional reconstructions of three pre-defined regions were obtained and microarchitectural parameters were calculated. Whereas there was no difference in bone volume fraction or structural model index between regions, trabecular number, thickness, spacing, connectivity density, degree of anisotropy and bone mineral density all displayed significant regional variations. The observed differences were consistent with the biomechanical hypothesis that in vivo loads are distributed differently at the endplates compared with the mid-vertebra. Thus, a more integrative approach combining biomechanical theory and anatomical features may improve fracture risk assessment in the future.

Contribution of trabecular and cortical components to biomechanical behavior of human vertebrae: an ex vivo study

Whereas there is clear evidence for a strong influence of bone quantity (i.e., bone mass or bone mineral density) on vertebral mechanical behavior, there are fewer data addressing the relative influence of cortical and trabecular bone microarchitecture. The aim of this study was to determine the relative contributions of bone mass, trabecular microarchitecture, and cortical thickness and curvature to the mechanical behavior of human lumbar vertebrae. Thirty-one L3 vertebrae (16 men, 15 women, aged 75 +/- 10 years and 76 +/- 10 years, respectively) were obtained. Bone mineral density (BMD) of the vertebral body was assessed by lateral dual energy X-ray absorptiometry (DXA), and 3D trabecular microarchitecture and anterior cortical thickness and curvature was assessed by micro-computed tomography (microCT). Then compressive stiffness, work to failure, and failure load were measured on the whole vertebral body. BMD was correlated with compressive stiffness (r = 0.60), failure load (r = 0.70), and work to failure (r = 0.55). Except for the degree of anisotropy, all trabecular and cortical parameters were correlated with mechanical behavior (r = 0.36 to 0.58, p = .05 to .001, and r = 0.36 to 0.61, p = .05 to .0001, respectively). Stepwise and multiple regression analyses indicated that the best predictor of (1) failure load was the combination of BMD, structural model index (SMI), and trabecular thickness (Tb.Th) (R = 0.80), (2) stiffness was the combination of BMD, Tb.Th, and curvature of the anterior cortex (R = 0.82), and (3) work to failure was the combination of anterior cortical thickness and BMD (R = 0.68). Our data imply that measurements of cortical thickness and curvature may enhance prediction of vertebral fragility and that therapies that improve both vertebral cortical and trabecular bone properties may provide a greater reduction in fracture risk.

Biomechanical effects of teriparatide in women with osteoporosis treated previously with alendronate and risedronate: results from quantitative computed tomography-based finite element analysis of the vertebral body

Previous antiresorptive treatment may influence the anabolic response to teriparatide. The OPTAMISE (Open-label Study to Determine How Prior Therapy with Alendronate or Risedronate in Postmenopausal Women with Osteoporosis Influences the Clinical Effectiveness of Teriparatide) study reported greater increases in biochemical markers of bone turnover and volumetric bone mineral density (BMD) when 12 months of teriparatide treatment was preceded by 2 years or more of risedronate versus alendronate treatment. The objective of this study was to use quantitative computed tomography (CT)-based nonlinear finite element modeling to evaluate how prior therapy with alendronate or risedronate in postmenopausal women with osteoporosis influences the biomechanical effectiveness of teriparatide. Finite element models of the L1 vertebra were created from quantitative CT scans, acquired before and after 12 months of therapy with teriparatide, from 171 patients from the OPTAMISE study. These models were subjected to uniaxial compression. Total BMD-derived bone volume fraction (BV/TV(d), i.e., bone volume [BV]/total volume [TV]), estimated from quantitative CT-based volumetric BMD, vertebral stiffness, and failure load (strength) were calculated for each time measurement point. The results of this study demonstrated that 12 months of treatment with teriparatide following prior treatment with either risedronate or alendronate increased BMD-derived BV/TV(d), the predicted vertebral stiffness, and failure load. However, the effects of teriparatide were more pronounced in patients treated previously with risedronate, which is consistent with the findings of the OPTAMISE study. The mean (+/-standard error) increase in stiffness was greater in the prior risedronate group than the prior alendronate group (24.6+/-3.2% versus 14.4+/-2.8%, respectively; p=0.0073). Similarly, vertebral failure load increased by 27.2+/-3.5% in the prior risedronate group versus 15.3+/-3.1% in the prior alendronate group (p=0.0042). The mechanical variables increased in greater proportion than BV/TV(d), which increased by 6.9+/-0.9% versus 4.6+/-0.8% in the prior-risedronate and prior-alendronate groups, respectively (p=0.0290). Our finding indicated that while teriparatide can be used with success on patients who have previously undergone treatment with risedronate and alendronate, it demonstrated greater anabolic effect on biomechanical properties in prior-risedronate patients in the first year of teriparatide treatment.

Role of trabecular microarchitecture in whole-vertebral body biomechanical behavior

The role of trabecular microarchitecture in whole-vertebral biomechanical behavior remains unclear, and its influence may be obscured by such factors as overall bone mass, bone geometry, and the presence of the cortical shell. To address this issue, 22 human T(9) vertebral bodies (11 female; 11 male; age range: 53-97 yr, 81.5 +/- 9.6 yr) were scanned with microCT and analyzed for measures of trabecular microarchitecture, BMC, cross-sectional area, and cortical thickness. Sixteen of the vertebrae were biomechanically tested to measure compressive strength. To estimate vertebral compressive stiffness with and without the cortical shell for all 22 vertebrae, two high-resolution finite element models per specimen-one intact model and one with the shell removed-were created from the microCT scans and virtually compressed. Results indicated that BMC and the structural model index (SMI) were the individual parameters most highly associated with strength (R(2) = 0.57 each). Adding microarchitecture variables to BMC in a stepwise multiple regression model improved this association (R(2) = 0.85). However, the microarchitecture variables in that regression model (degree of anisotropy, bone volume fraction) differed from those when BMC was not included in the model (SMI, mean trabecular thickness), and the association was slightly weaker for the latter (R(2) = 0.76). The finite element results indicated that the physical presence of the cortical shell did not alter the relationships between microarchitecture and vertebral stiffness. We conclude that trabecular microarchitecture is associated with whole-vertebral biomechanical behavior and that the role of microarchitecture is mediated by BMC but not by the cortical shell.

Development and validation of a predictive bone fracture risk model for astronauts

There are still many unknowns in the physiological response of human beings to space, but compelling evidence indicates that accelerated bone loss will be a consequence of long-duration spaceflight. Lacking phenomenological data on fracture risk in space, we have developed a predictive tool based on biomechanical and bone loading models at any gravitational level of interest. The tool is a statistical model that forecasts fracture risk, bounds the associated uncertainties, and performs sensitivity analysis. In this paper, we focused on events that represent severe consequences for an exploration mission, specifically that of spinal fracture resulting from a routine task (lifting a heavy object up to 60 kg), or a spinal, femoral or wrist fracture due to an accidental fall or an intentional jump from 1 to 2 m. We validated the biomechanical and bone fracture models against terrestrial studies of ground reaction forces, skeletal loading, fracture risk, and fracture incidence. Finally, we predicted fracture risk associated with reference missions to the moon and Mars that represented crew activities on the surface. Fracture was much more likely on Mars due to compromised bone integrity. No statistically significant gender-dependent differences emerged. Wrist fracture was the most likely type of fracture, followed by spinal and hip fracture.

Prediction of vertebral strength under loading conditions occurring in activities of daily living using a computed tomography-based nonlinear finite element method

STUDY DESIGN

A clinical study on osteoporotic vertebral strength in daily living using a computed tomography (CT)-based nonlinear finite element (FE) model.

OBJECTIVE

To evaluate the differences in predicted fracture strength of osteoporotic vertebral bodies among the different loading conditions that are occurring in the activities of daily living.

SUMMARY OF BACKGROUND DATA

FE model has been reported to predict vertebral strength in uniaxial loading, but forward bending load plays an important role in osteoporotic vertebral fractures.

METHODS

Strengths of the second lumbar vertebra in 41 female patients with postmenopausal osteoporosis were analyzed using a nonlinear CT-based FE method. Three different loading conditions were adopted uniaxial compression, forward bending, and erect standing. The same boundary condition was used for all loading conditions. Predicted strengths under forward bending and erect standing were compared with that under uniaxial compression and differences in strength were statistically analyzed.

RESULTS

The regression equation relating strength under uniaxial compression to that under erect standing was expressed as y = 0.8912x + 19.332 (R = 0.9522), whereas the equation relating uniaxial compression to forward bending was y = 0.7033x + 55.071 (R = 0.8342). Both relationships were significant, but the correlation between forward bending and uniaxial compression was not strong, while strength was lower under forward bending than under uniaxial compression according to the Friedman multiple comparison test (P = 0.00017).

CONCLUSION

Strength under forward bending correlated significantly to that under uniaxial compression, but the correlation was not strong. Therefore, in osteoporotic patients, both uniaxial compression and forward bending should be assessed to evaluate fracture risk in daily living using a CT-based FE method.

Relative strength of thoracic vertebrae in axial compression versus flexion

BACKGROUND CONTEXT

Noninvasive strength assessment techniques are the clinical standard in the diagnosis and treatment of osteoporotic vertebral fractures, and the efficacy of these protocols depends on their ability to predict vertebral strength at all at-risk spinal levels under multiple physiological loading conditions.

PURPOSE

To assess differences in vertebral strength between loading modes and across spinal levels.

STUDY DESIGN/SETTING

This study examined the relative strength of isolated vertebral bodies in compression versus flexion.

METHODS

Destructive biomechanical tests were conducted on 30 pairs of donor-matched, isolated thoracic vertebral bodies (T9 and T10; F=19, M=11; 87+5 years old, max=97 years old, min=80 years old) in both uniform axial compression and flexion using previously described protocols. Quantitative computed tomography (QCT) scans were taken before mechanical testing and used to obtain bone mineral density (BMD) and "mechanics of solids" (MOS) measures, such as axial and bending rigidities.

RESULTS

Compressive strength was higher than flexion strength for each donor by 940+152N (p<.001, paired t test), and vertebral strengths in the two loading modes were moderately correlated (adjusted R(2)=0.50, p<.001). For both compression and flexion loading modes, adjacent-level BMD and MOS metrics had approximately half the predictive capacity as same-level measurements, and BMD and MOS values were only moderately correlated across spinal levels.

CONCLUSIONS

The results of this study are important in designing clinical test protocols for assessing vertebral fracture risk. Because vertebral body flexion and compressive strength are not strongly correlated and flexion strength is significantly less than compressive strength, it is imperative to investigate a patient's spinal structural capacity under bending loading conditions. Furthermore, our work suggests that clinicians using QCT-based measures should perform site-specific strength assessments on each at-risk spinal level. Future work should focus on improving the accuracy of densitometric measures in predicting vertebral strength in flexion and also on examining same- versus adjacent-level strength assessment for radiographic techniques with lower X-ray dosage, such as dual-energy X-ray absorptiometry.

Vertebral strength changes in rheumatoid arthritis patients treated with alendronate, as assessed by finite element analysis of clinical computed tomography scans: a prospective randomized clinical trial

OBJECTIVE

Finite element analysis of clinical computed tomography (CT) scans provides a noninvasive means of assessing vertebral strength that is superior to dual x-ray absorptiometry (DXA)-measured areal bone mineral density. The present study was undertaken to compare strength changes, measured using this newer method, in rheumatoid arthritis (RA) patients who were treated with alendronate (ALN) versus those who were not.

METHODS

Thirty female RA patients without radiologic signs of L3 compression fractures or a history of osteoporosis medication were enrolled in a prospective randomized clinical trial. Patients were randomly assigned to the ALN group (5 mg orally, once daily) or the control group not receiving antiresorptive treatment. All patients were evaluated by DXA and quantitative CT at baseline and reevaluated after a mean of 12.2 months. Nonlinear finite element analysis was performed on the CT scans (n = 29 available for analysis) to compute an estimate of vertebral compressive strength and to assess strength changes associated with changes in the trabecular compartment and the outer 2 mm of bone (peripheral compartment).

RESULTS

On average, vertebral strength was significantly decreased from baseline in the control group (n = 15) (median change -10.6%; P = 0.008) but was maintained in the ALN group (n = 14) (median change +0.4%; P = 0.55), with a significant difference between the 2 groups (P < 0.01). Strength decreased more rapidly within the trabecular bone, and ALN treatment was much more effective in the peripheral than the trabecular compartment.

CONCLUSION

Our results indicate that patients with RA can lose a substantial amount of vertebral strength over a relatively short period of time, and this loss can be prevented by ALN, primarily via its positive effect on the outer 2 mm of vertebral bone.

Structural determinants of vertebral fracture risk

Vertebral fractures are more strongly associated with specific bone density, structure, and strength parameters than with areal BMD, but all of these variables are correlated.

INTRODUCTION

It is unclear whether the association of areal BMD (aBMD) with vertebral fracture risk depends on bone density per se, bone macro- or microstructure, overall bone strength, or spine load/bone strength ratios.

MATERIALS AND METHODS

From an age-stratified sample of Rochester, MN, women, we identified 40 with a clinically diagnosed vertebral fracture (confirmed semiquantitatively) caused by moderate trauma (cases; mean age, 78.6 +/- 9.0 yr) and compared them with 40 controls with no osteoporotic fracture (mean age, 70.9 +/- 6.8 yr). Lumbar spine volumetric BMD (vBMD) and geometry were assessed by central QCT, whereas microstructure was evaluated by high-resolution pQCT at the ultradistal radius. Vertebral failure load ( approximately strength) was estimated from voxel-based finite element models, and the factor-of-risk (phi) was determined as the ratio of applied spine loads to failure load.

RESULTS

Spine loading (axial compressive force on L3) was similar in vertebral fracture cases and controls (e.g., for 90 degrees forward flexion, 2639 versus 2706 N; age-adjusted p = 0.173). However, fracture cases had inferior values for most bone density and structure variables. Bone strength measures were also reduced, and the factor-of-risk (phi) was 35-37% greater (worse) among women with a vertebral fracture. By age-adjusted logistic regression, relative risks for the strongest fracture predictor in each of the five main variable categories were bone density (total lumbar spine vBMD: OR per SD change, 2.2; 95% CI, 1.1-4.3), bone geometry (vertebral apparent cortical thickness: OR, 2.1; 95% CI, 1.1-4.1), bone microstructure (none significant); bone strength ("cortical" [outer 2 mm] compressive strength: OR, 2.5; 95% CI, 1.3-4.8), and factor-of-risk (phi for 90 degrees forward flexion/overall vertebral compressive strength: OR, 3.2; 95% CI, 1.4-7.5). These variables were correlated with spine aBMD (partial r, -0.32 to 0.75), but each was a stronger predictor of fracture in the logistic regression analyses.

CONCLUSIONS

The association of aBMD with vertebral fracture risk is explained by its correlation with more specific bone density, structure, and strength parameters. These may allow deeper insights into fracture pathogenesis.

Locations of bone tissue at high risk of initial failure during compressive loading of the human vertebral body

Knowledge of the location of initial regions of failure within the vertebra - cortical shell, cortical endplates vs. trabecular bone, as well as anatomic location--may lead to improved understanding of the mechanisms of aging, disease and treatment. The overall objective of this study was to identify the location of the bone tissue at highest risk of initial failure within the vertebral body when subjected to compressive loading. Toward this end, micro-CT-based 60-micron voxel-sized, linearly elastic, finite element models of a cohort of thirteen elderly (age range: 54-87 years, 75+/-9 years) female whole vertebrae without posterior elements were virtually loaded in compression through a simulated disc. All bone tissues within each vertebra having either the maximum or minimum principal strain beyond its 90th percentile were defined as the tissue at highest risk of initial failure within that particular vertebral body. Our results showed that such high-risk tissue first occurred in the trabecular bone and that the largest proportion of the high-risk tissue also occurred in the trabecular bone. The amount of high-risk tissue was significantly greater in and adjacent to the cortical endplates than in the mid-transverse region. The amount of high-risk tissue in the cortical endplates was comparable to or greater than that in the cortical shell regardless of the assumed Poisson's ratio of the simulated disc. Our results provide new insight into the micromechanics of failure of trabecular and cortical bone within the human vertebra, and taken together, suggest that, during strenuous compressive loading of the vertebra, the tissue near and including the endplates is at the highest risk of initial failure.

Quantitative computed tomography-based predictions of vertebral strength in anterior bending

STUDY DESIGN

This study examined the ability of QCT-based structural assessment techniques to predict vertebral strength in anterior bending.

OBJECTIVE

The purpose of this study was to compare the abilities of QCT-based bone mineral density (BMD), mechanics of solids models (MOS), e.g., bending rigidity, and finite element analyses (FE) to predict the strength of isolated vertebral bodies under anterior bending boundary conditions.

SUMMARY OF BACKGROUND DATA

Although the relative performance of QCT-based structural measures is well established for uniform compression, the ability of these techniques to predict vertebral strength under nonuniform loading conditions has not yet been established.

METHODS

Thirty human thoracic vertebrae from 30 donors (T9-T10, 20 female, 10 male; 87 +/- 5 years of age) were QCT scanned and destructively tested in anterior bending using an industrial robot arm. The QCT scans were processed to generate specimen-specific FE models as well as trabecular bone mineral density (tBMD), integral bone mineral density (iBMD), and MOS measures, such as axial and bending rigidities.

RESULTS

Vertebral strength in anterior bending was poorly to moderately predicted by QCT-based BMD and MOS measures (R2 = 0.14-0.22). QCT-based FE models were better strength predictors (R2 = 0.34-0.40); however, their predictive performance was not statistically different from MOS bending rigidity (P > 0.05).

CONCLUSIONS

Our results suggest that the poor clinical performance of noninvasive structural measures may be due to their inability to predict vertebral strength under bending loads. While their performance was not statistically better than MOS bending rigidities, QCT-based FE models were moderate predictors of both compressive and bending loads at failure, suggesting that this technique has the potential for strength prediction under nonuniform loads. The current FE modeling strategy is insufficient, however, and significant modifications must be made to better mimic whole bone elastic and inelastic material behavior.

Comparison of quantitative computed tomography-based measures in predicting vertebral compressive strength

Patient-specific measures derived from quantitative computed tomography (QCT) scans are currently being developed as a clinical tool for vertebral strength prediction. QCT-based measurement techniques vary greatly in structural complexity and generally fall into one of three categories: (1) bone mineral density (BMD), (2) "mechanics of solids" (MOS) models, such as minimum axial rigidity (the product of axial stiffness and vertebral height), or (3) three-dimensional finite element (FE) models. There is no clear consensus as to the relative performance of these measures due to differences in experimental protocols, sample sizes and demographics, and outcome metrics. The goal of this study was to directly compare the performance of QCT-based assessment techniques of varying degrees of structural sophistication in predicting experimental vertebral compressive strength. Eighty-one human thoracic vertebrae (T6-T10) from 44 donors cadavers (F=32, M=12; 85+/-8 years old, max=97 years old, min=54 years old) were QCT scanned and destructively tested in uniaxial compression. The QCT scans were processed to generate FE models and various BMD and MOS measures, including trabecular bone mineral density (tBMD), integral bone mineral density (iBMD), and axial rigidity. Bone mineral density was weakly to moderately predictive of compressive strength (R(2)=0.16 and 0.62 for tBMD and iBMD, respectively). In vitro vertebral strength was strongly correlated with both axial rigidity (R(2)=0.81) and FE strength measurements (R(2)=0.80), and the predictive capabilities of these two metrics were statistically equivalent (p>0.05 for differences between FE and axial rigidity). The results of this study indicate that non-invasive predictive measures of vertebral strength should include some level of structural sophistication, specifically, gross geometric and material property distribution information. For uniaxial compression of isolated vertebrae, which is the current biomechanical testing paradigm for new non-invasive strength assessment techniques, QCT-based FE and axial rigidity measures are equivalent predictors of experimental strength. However, before abandoning the FE method in favor of more simplistic techniques, future work should investigate the performance of the FE method versus MOS measures for more physiologically representative loading conditions, e.g., anterior bending or in situ loading with intervertebral discs intact.

Regional trabecular morphology assessed by micro-CT is correlated with failure of aged thoracic vertebrae under a posteroanterior load and may determine the site of fracture

INTRODUCTION

Spinal mobilization is commonly used in the treatment of patients with back pain, including individuals with osteoporosis. Previous data indicated that traditional predictors of skeletal failure-lateral or anteroposterior bone mineral density (BMD) by dual energy X-ray absorptiometry (DXA) or geometry of the spinous process or vertebral body-do not predict failure load during posteroanterior spinal mobilization. Morphological differences and inhomogeneities in BMD may have important effects on vertebral strength but integral BMD values by DXA cannot reflect these potentially important differences. We investigated the determinants of spinal fracture using muCT.

MATERIALS AND METHODS

We measured failure load and failure site in 11 T5-8 cadaveric specimens (mean age 78 years) when a posteroanterior load was applied at the spinous process of T6 using a servohydraulic material testing machine. Radiography and CT scan were used to verify failure site. We observed no damage to the adjacent T7 vertebrae following the T6 posteroanterior failure test. The T7 vertebrae were sectioned to produce regional samples of the spinous process, the lamina and a vertebral body core. Each sample was scanned with muCT to measure bone microarchitectural parameters. We segmented and analysed four trabecular regions (spinous process base and middle, central lamina and central vertebral body). We used one-way repeated measures ANOVA to compare regions and computed Pearson correlations to assess the relation between PA failure load of T6 and the morphological parameters of T7.

RESULTS

The BV/TV at the base or middle of the T7 spinous process (fracture sites), Tb.N and Tb.Th at the base were significantly correlated with posteroanterior failure load of T6 (BV/TV base: r=0.74, p=0.01; BV/TV middle: r=0.73, p=0.01; Tb.N base: r=0.64, p=0.03; Tb.Th base: r=0.65, p=0.03). The Tb.Th of the lamina was significantly greater than Tb.Th of the spinous process base (p=0.002).

CONCLUSIONS

Whereas previous data indicated that BMD by DXA was not a good predictor of posteroanterior failure load, regional BV/TV of the spinous process base and middle regions, the sites of fracture, are correlated with posteroanterior failure load. Trabecular thickness differed significantly between the base of the spinous process and the lamina, and may have influenced the site of fracture.

Effects of teriparatide and alendronate on vertebral strength as assessed by finite element modeling of QCT scans in women with osteoporosis

FE modeling was used to estimate the biomechanical effects of teriparatide and alendronate on lumbar vertebrae. Both treatments enhanced predicted vertebral strength by increasing average density. This effect was more pronounced for teriparatide, which further increased predicted vertebral strength by altering the distribution of density within the vertebra, preferentially increasing the strength of the trabecular compartment.

INTRODUCTION

Teriparatide 20 microg/day (TPTD) and alendronate 10 mg/day (ALN) increase areal, measured by DXA, and volumetric, measured by QCT, lumbar spine BMD through opposite effects on bone remodeling. Using finite element (FE) modeling of QCT scans, we sought to compare the vertebral strength characteristics in TPTD- and ALN-treated patients.

MATERIALS AND METHODS

A subset of patients (N = 28 TPTD; N = 25 ALN) from the Forteo Alendronate Comparator Trial who had QCT scans of the spine at baseline and postbaseline were analyzed. The QCT scans were analyzed for compressive strength of the L(3) vertebra using FE modeling. In addition, using controlled parameter studies of the FE models, the effects of changes in density, density distribution, and geometry on strength were calculated, a strength:density ratio was determined, and a response to bending was also quantified.

RESULTS

Both treatments had positive effects on predicted vertebral strength characteristics. At least 75% of the patients in each treatment group had increased strength of the vertebra at 6 months compared with baseline. Patients in both treatment groups had increased average volumetric density and increased strength in the trabecular bone, but the median percentage increases for these parameters were 5- to 12-fold greater for TPTD. Larger increases in the strength:density ratio were also observed for TPTD, and these were primarily attributed to preferential increases in trabecular strength.

CONCLUSIONS

These results provide new insight into the effects of these treatments on estimated biomechanical properties of the vertebra. Both treatments positively affected predicted vertebral strength through their effects on average BMD, but the magnitudes of the effects were quite different. Teriparatide also affected vertebral strength by altering the distribution of density within the vertebra, so that overall, teriparatide had a 5-fold greater percentage increase in the strength:density ratio.

Trabecular bone density of male human cervical and lumbar vertebrae

The objective of this study was to determine the bone mineral density (BMD) of cervical vertebrae and correlate with the lumbar spine. Fifty-seven young adult healthy male volunteers, ranging from 18 to 41 years of age, underwent quantitative computed tomography (QCT) scanning of C2-T1 and L2-L4 vertebrae. To account for correlations, repeated measures techniques were used to compare data as a function of spinal level and region. Linear regression methods were used (+/-95% CI) to compare data as a function of spinal level and region. The mean age and body height were 25.0 +/- 5.8 years and 181.0 +/- 7.6 cm. BMD decreased from the rostral to caudal direction along the spinal column. Grouped data indicated that the neck is the densest followed by the first thoracic vertebra and low back with mean BMD of 256.0 +/- 48.1, 194.3 +/- 44.2, and 172.2 +/- 28.4 mg/cm(3), respectively; differences were statistically significant. While BMD did not vary significantly between the three lumbar bodies, neck vertebrae demonstrated significant trends. The matrix of correlation coefficients between BMD and spinal level indicated that the relationship is strong in the lumbar (r = 0.92-0.96) and cervical (r = 0.73-0.92) spines. Data from the present study show that the trabecular bony architecture of the neck is significantly different from the low back. These quantitative BMD data from a controlled young adult healthy human male volunteer population may be valuable in establishing normative data specifically for the neck. From a trabecular bone density perspective, these results indicate that lumbar vertebrae cannot act as the best surrogates for neck vertebrae. Significant variations in densities among neck vertebrae, unlike the low back counterpart, may underscore the need to treat these bones as different structures.

Sensitivity of Vertebral Compressive Strength to Endplate Loading Distribution

The sensitivity of vertebral body strength to the distribution of axial forces along the endplate has not been comprehensively evaluated. Using quantitative computed tomography-based finite element models of 13 vertebral bodies, an optimization analysis was performed to determine the endplate force distributions that minimized (lower bound) and maximized (upper bound) vertebral strength for a given set of externally applied axial compressive loads. Vertebral strength was also evaluated for three generic boundary conditions: uniform displacement, uniform force, and a nonuniform force distribution in which the interior of the endplate was loaded with a force that was 1.5 times greater than the periphery. Our results showed that the relative difference between the upper and lower bounds on vertebral strength was 14.2±7.0%(mean±SD). While there was a weak trend for the magnitude of the strength bounds to be inversely proportional to bone mineral density (R2=0.32, p=0.02), both upper and lower bound vertebral strength measures were well predicted by the strength response under uniform displacement loading conditions (R2=0.91 and R2=0.99, respectively). All three generic boundary conditions resulted in vertebral strength values that were statistically indistinguishable from the loading condition that resulted in an upper bound on strength. The results of this study indicate that the uncertainty in strength arising from the unknown condition of the disc is dependent on the condition of the bone (whether it is osteoporotic or normal). Although bone mineral density is not a good predictor of strength sensitivity, vertebral strength under generic boundary conditions, i.e., uniform displacement or force, was strongly correlated with the relative magnitude of the strength bounds. Thus, explicit disc modeling may not be necessary.

Metastatic Burst Fracture Risk Assessment Based on Complex Loading of the Thoracic Spine

The mechanical integrity of vertebral bone is compromised when metastatic cancer cells migrate to the spine, rendering it susceptible to burst fracture under physiologic loading. Risk of burst fracture has been shown to be dependent on the magnitude of the applied load, however limited work has been conducted to determine the effect of load type on the stability of the metastatic spine. The objective of this study was to use biphasic finite element modeling to evaluate the effect of multiple loading conditions on a metastatically-involved thoracic spinal motion segment. Fifteen loading scenarios were analyzed, including axial compression, flexion, extension, lateral bending, torsion, and combined loads. Additional analyses were conducted to assess the impact of the ribcage on the stability of the thoracic spine. Results demonstrate that axial loading is the predominant load type leading to increased risk of burst fracture initiation, while rotational loading led to only moderate increases in risk. Inclusion of the ribcage was found to reduce the potential for burst fracture by 27%. These findings are important in developing a more comprehensive understanding of burst fracture mechanics and in directing future modeling efforts. The results in this study may also be useful in advising less harmful activities for patients affected by lytic spinal metastases.

Cortical and trabecular load sharing in the human vertebral body

The biomechanical role of the vertebral cortical shell remains poorly understood. Using high-resolution finite element modeling of a cohort of elderly vertebrae, we found that the biomechanical role of the shell can be substantial and that the load sharing between the cortical and trabecular bone is complex. As a result, a more integrative measure of the trabecular and cortical bone should improve noninvasive assessment of fracture risk and treatments.

INTRODUCTION

A fundamental but poorly understood issue in the assessment of both osteoporotic vertebral fracture risk and effects of treatment is the role of the trabecular bone and cortical shell in the load-carrying capacity of the vertebral body.

MATERIALS AND METHODS

High-resolution microCT-based finite element models were developed for 13 elderly human vertebrae (age range: 54-87 years; 74.6 +/- 9.4 years), and parameter studies-with and without endplates-were performed to determine the role of the shell versus trabecular bone and the effect of model assumptions.

RESULTS

Across vertebrae, whereas the average thickness of the cortical shell was only 0.38 +/- 0.06 mm, the shell mass fraction (shell mass/total bone mass)-not including the endplates-ranged from 0.21 to 0.39. The maximum load fraction taken by the shell varied from 0.38 to 0.54 across vertebrae and occurred at the narrowest section. The maximum load fraction taken by the trabecular bone varied from 0.76 to 0.89 across vertebrae and occurred near the endplates. Neither the maximum shell load fraction nor the maximum trabecular load fraction depended on any of the densitometric or morphologic properties of the vertebra, indicating the complex nature of the load sharing mechanism. The variation of the shell load-carrying capacity across vertebrae was significantly altered by the removal of endplates, although these models captured the overall trend within a vertebra.

CONCLUSIONS

The biomechanical role of the thin cortical shell in the vertebral body can be substantial, being about 45% at the midtransverse section but as low as 15% close to the endplates. As a result of the complexity of load sharing, sampling of only midsection trabecular bone as a strength surrogate misses important biomechanical information. A more integrative approach that combines the structural role of both cortical and trabecular bone should improve noninvasive assessment of vertebral bone strength in vivo.