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

Finite element analysis of biomechanical effects of percutaneous cement discoplasty in scoliosis

OBJECTIVE

To investigate the effect of bone cement on the vertebral body and biomechanical properties in percutaneous cement discoplasty (PCD) for degenerative lumbar disc disease.

METHODS

Three-dimensional reconstruction of L2 ~ L3 vertebral bodies was performed in a healthy volunteer, and the corresponding finite element model of the spine was established. Biomechanical analysis was performed on the changes in stress distribution in different groups of models by applying quantitative loads.

RESULTS

Models with percutaneous discoplasty (PCD) showed improved stability under various stress conditions, and intervertebral foraminal heights were superior to models without discoplasty.

CONCLUSION

Cement discoplasty can improve the stability of the vertebral body to a certain extent and restore a certain height of the intervertebral foramen, which has a good development prospect and potential.

Spinal stability analysis of lumbar interbody fusion according to pelvic type and cage angle based on simplified spinal model with various pelvic indices

Recently, the objectives of lumbar interbody fusion (LIF) have been extended to include the correction of broader/relative indications in addition to spinal fixation. Accordingly, LIF must be optimized for sagittal alignment while simultaneously achieving decompression. Therefore, a representative model classified into three pelvic types, i.e., neutral pelvis (NP), anterior pelvis (AP), and retroverted pelvis (RP), was selected according to the pelvic index, and LIF was performed on each representative model to analyze Lumbar lordosis (LL) and the corresponding equivalent stress. The finite element (FE) model was based on a sagittal 2D X-ray image. The calculation efficiency and convergence were improved by simplifying the modeling of the vertebral body in general and its posterior portion in particular. Based on the position of the pelvis, according to the pelvic shape, images of patients were classified into three types: AP, RP, and NP. Subsequently, representative images were selected for each type. The fixation device used in the fusion model was a pedicle screw and a spinal rod of a general type. PEEK was used as the cage material, and the cage shape was varied by using three different cage angles: 0°, 4°, and 8°. Spinal mobility: The pelvic type with the highest range of motion (ROM) for the spine was the NP type; the AP type had the highest LL. Under a combination load, the NP type exhibited the highest lumbar flexibility (LF), which was 2.46° lower on average compared to the case where a pure moment was applied. Equivalent stress on the spinal fixation device: The equivalent stress acting on the vertebrae was lowest when cage 0 was used for the NP and AP type. For the RP type, the lowest equivalent stress on the vertebrae was observed when cage 4 was used. Finally, for the L5 upper endplate, the stress did not vary significantly for a given type of cage. In conclusion, there was no significant difference in ROM according to cage angle, and the highest ROM, LL and LF were shown in the pelvic shape of NP type. However, when comparing the results with other pelvic types, it was not possible to confirm that LF is completely dependent on LL and ROM.

Romosozumab Enhances Vertebral Bone Structure in Women With Low Bone Density

Romosozumab monoclonal antibody treatment works by binding sclerostin and causing rapid stimulation of bone formation while decreasing bone resorption. The location and local magnitude of vertebral bone accrual by romosozumab and how it compares to teriparatide remains to be investigated. Here we analyzed the data from a study collecting lumbar computed tomography (CT) spine scans at enrollment and 12 months post-treatment with romosozumab (210 mg sc monthly, n = 17), open-label daily teriparatide (20 μg sc, n = 19), or placebo (sc monthly, n = 20). For each of the 56 women, cortical thickness (Ct.Th), endocortical thickness (Ec.Th), cortical bone mineral density (Ct.bone mineral density (BMD)), cancellous BMD (Cn.BMD), and cortical mass surface density (CMSD) were measured across the first lumbar vertebral surface. In addition, color maps of the changes in the lumbar vertebrae structure were statistically analyzed and then visualized on the bone surface. At 12 months, romosozumab improved all parameters significantly over placebo and resulted in a mean vertebral Ct.Th increase of 10.3% versus 4.3% for teriparatide, an Ec.Th increase of 137.6% versus 47.5% for teriparatide, a Ct.BMD increase of 2.1% versus a -0.1% decrease for teriparatide, and a CMSD increase of 12.4% versus 3.8% for teriparatide. For all these measurements, the differences between romosozumab and teriparatide were statistically significant (p < 0.05). There was no significant difference between the romosozumab-associated Cn.BMD gains of 22.2% versus 18.1% for teriparatide, but both were significantly greater compared with the change in the placebo group (-4.6%, p < 0.05). Cortical maps showed the topographical locations of the increase in bone in fracture-prone areas of the vertebral shell, walls, and endplates. This study confirms widespread vertebral bone accrual with romosozumab or teriparatide treatment and provides new insights into how the rapid prevention of vertebral fractures is achieved in women with osteoporosis using these anabolic agents. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Stochastic parametric skeletal dosimetry model for humans: General approach and application to active marrow exposure from bone-seeking beta-particle emitters

The objective of this study is to develop a skeleton model for assessing active marrow dose from bone-seeking beta-emitting radionuclides. This article explains the modeling methodology which accounts for individual variability of the macro- and microstructure of bone tissue. Bone sites with active hematopoiesis are assessed by dividing them into small segments described by simple geometric shapes. Spongiosa, which fills the segments, is modeled as an isotropic three-dimensional grid (framework) of rod-like trabeculae that “run through” the bone marrow. Randomized multiple framework deformations are simulated by changing the positions of the grid nodes and the thickness of the rods. Model grid parameters are selected in accordance with the parameters of spongiosa microstructures taken from the published papers. Stochastic modeling of radiation transport in heterogeneous media simulating the distribution of bone tissue and marrow in each of the segments is performed by Monte Carlo methods. Model output for the human femur at different ages is provided as an example. The uncertainty of dosimetric characteristics associated with individual variability of bone structure was evaluated. An advantage of this methodology for the calculation of doses absorbed in the marrow from bone-seeking radionuclides is that it does not require additional studies of autopsy material. The biokinetic model results will be used in the future to calculate individual doses to members of a cohort exposed to ^89,90Sr from liquid radioactive waste discharged to the Techa River by the Mayak Production Association in 1949–1956. Further study of these unique cohorts provides an opportunity to gain more in-depth knowledge about the effects of chronic radiation on the hematopoietic system. In addition, the proposed model can be used to assess the doses to active marrow under any other scenarios of ⁹⁰Sr and ⁸⁹Sr intake to humans.

The micro-structural analysis of lumbar vertebrae in alcoholic liver cirrhosis

Although vertebral fracture is more common among alcoholic liver cirrhosis patients when compared to general population, current data on three-dimensional micro-architecture are scarce. Our study showed significant trabecular deterioration in lumbar vertebrae obtained from alcoholic liver cirrhosis donors, suggesting that they should be advised to undergo early-stage screening for osteoporosis.

PURPOSE

Recent studies showed an increased incidence of vertebral fractures in alcoholic liver cirrhosis (ALC) patients, while data about vertebral micro-structure are still limited. The aim of this study was to compare trabecular and cortical micro-architecture of lumbar vertebrae between ALC patients and healthy age- and sex-matched controls.

METHODS

Our study included lumbar vertebral samples of male cadaveric donors, divided into ALC (n = 20, age: 59 ± 6 years) and control group (n = 20, age: 59 ± 8 years). Following pathohistological verification of liver cirrhosis, trabecular and cortical bone micro-architecture was analyzed by micro-computed tomography (micro-CT).

RESULTS

Micro-CT evaluation of the trabecular bone in lumbar vertebrae showed a significant decrease in bone volume fraction, trabecular thickness, trabecular number, and connectivity (p < 0.01). In contrast to trabecular deterioration, prominent alteration in cortical parameters was not observed in lumbar vertebrae of ALC patients (p > 0.05).

CONCLUSIONS

Our data indicate that susceptibility to non-traumatic fractures in ALC patients could be explained by alterations in trabecular bone micro-architecture. Thus, we genuinely recommend osteological screening of the lumbar spine for all ALC patients in order to evaluate individual fracture risk. Graphical abstract.

Influence of different postures under vertical impact load on thoracolumbar burst fracture

Clinical studies have extensively shown that burst fractures can cause severe and long-term neurological deficits. However, the mechanism of burst fracture is not clear, and the influence of different spinal postures on burst fracture is still unidentified. The study aimed at investigating the influence of different postures under vertical impact load on thoracolumbar burst fracture. A detailed nonlinear finite element model of T12-L2 segment was developed to investigate these problems. In this work, a rigid ball was used to vertically impact the finite element spinal segment, which emulated the process of burst fracture as in experimental condition. During the process, amounting to 9 different postures (normal, flexion, extension, right/left lateral bending of 8°, right/left axial rotation of 4° and 8°) were studied. Totally five failure modes were observed. Six different parameters, including vertebral height, vertebral bulge, interpedicular widening, vertebral kyphotic angle, posterior vertebral body angle, and joint facet contact force, were observed to evaluate the corresponding severity of burst fracture. Burst fracture in extension was the severest, and the loss of vertebral height in flexion was the most. The different postures in these simulations changed the morphology of intervertebral disc and facet joints force, resulting in different types of fracture.

Differences in Trabecular Bone, Cortical Shell, and Endplate Microstructure Across the Lumbar Spine

Background

The quality of the vertebral body structures such as endplate, cortex, and trabecular bone is important for understanding the performance of implants, particularly at the bone-implant interface. Although vertebral body structures have been analyzed separately in the literature, there is no comprehensive study to assess these anatomical measurements along with their interrelationships in the lumbar spine. Therefore, the purpose of this study was to assess variations in trabecular microstructure, vertebral endplate thickness and concavity, and vertebral body cortex thickness within the lumbar spine.

Methods

A total of 80 lumbar vertebrae (L1-L5) were dissected from 16 human cadaver specimens and imaged with microcomputed tomography to determine trabecular microstructure, vertebral cortex thickness, endplate thickness, and maximum endplate concavity depth. A paired t test and regression analysis were used to determine significant differences (P < .05) between different vertebral levels and correlations between the analyzed anatomical parameters.

Results

L1 vertebra had significantly better (P < .02) trabecular bone microstructure (eg, trabecular bone volume fraction) than all other lumbar vertebrae. However, L1 vertebra also had significantly thinner (P ≤ .02) anterior, left, and right cortices compared to all other vertebral levels. Within L3-L5 intervertebral disc spaces, cranial endplates had significantly greater (P ≤ .03) thickness and maximum concavity depth compared to their respective caudal endplates. No strong correlations were observed between trabecular bone microstructure, maximum endplate concavity depth, vertebral cortex, and endplate thickness parameters.

Conclusions

Detailed reference data of these anatomical parameters for each lumbar vertebral body can aid in improved understanding of bone quality, particularly when assessing different implant designs and fixation approaches. Moreover, such anatomical knowledge may help clinicians with optimal implant design selection and surgical placement of these devices into their respective locations.

A biomechanical investigation of thoracolumbar burst fracture under vertical impact loads using finite element method

BACKGROUND

A sudden vertical impact load on spine can cause spinal burst fracture, especially in the thoracolumbar junction region. This study aimed at investigating the mechanism of spinal burst fracture under different energy vertical impact loads, producing the failure risk region to understand burst fracture, reducing nervous system damage and guiding clinical treatment.

METHODS

A nonlinear finite element model of T12-L1 motion segment was created to analyze the response of the vertical impact load. A rigid ball was used to impact the segment vertically to simulate the vertical impact load in practice. There were three different mass balls to represent the different loads: low energy, intermediate energy and high energy (respectively 13 J, 30 J and 56 J). The results of impact force, vertical displacement, stress, intradiscal pressure and contact force were obtained during the process.

FINDINGS

At low energy condition, the rigid ball rebounded rapidly. At intermediate energy condition, fractures were initiated in vertebral foramen and left rear regions on the superior cortical bone near the superior endplate of L1. At high energy condition, burst fracture occurred and a part of L1 was isolated from the model.

INTERPRETATION

The fracture occurred on the L1 segment only at the intermediate energy and high energy. The strength of vertebral body under low and intermediate energy was enough to support the impact. The burst fracture pattern at high energy was also observed in clinical practice. The findings may explain the mechanism of burst fracture.

Improved accuracy in the assessment of vertebral cortical thickness by quantitative computed tomography using the Iterative Convolution OptimizatioN (ICON) method

Vertebral whole bone strength is substantially affected by cortical bone properties. Disease and therapy may affect cancellous and cortical bone differently. Unlike Dual X-ray Absorptiometry (DXA), Quantitative Computed Tomography (QCT) permits selective assessment of cortical and cancellous bone, but image quality limits the accuracy. We present an image processing method specifically adopted to thin cortices that substantially improves accuracy. Ten human vertebrae embedded in epoxy resin were imaged using clinical QCT and High-Resolution QCT (HR-QCT) protocols, both acquired on a clinical whole body CT scanner, whereas high resolution peripheral QCT (HR-pQCT) was used as gold standard. Microstructural variables and BMD were calculated using in-house software StructuralInsight for QCT and HR-QCT and the manufacturer's μCT evaluation software for HR-pQCT. An adjusted measure, a deconvolved cortical thickness (dcCt.Th), corrected for partial volume effects, was derived applying the new Iterative Convolution OptimizatioN (ICON) method. Direct measurements of cortical thickness (Ct.Th) showed substantial overestimation with mean ± standard deviation of 1.8 ± 0.5 mm for QCT and 1.5 ± 0.3 mm for HR-QCT compared to 0.37 ± 0.07 mm using HR-pQCT. Correlations of both QCT (r² = 0.05, p > 0.5.) and HR-QCT (r² = 0.38, p = 0.060) with the gold standard HR-pQCT were not significant. Also QCT-based BMD and BMC as well as HR-QCT-based BMD did not show a significant correlation with the gold standard approach. Only HR-QCT-based BMC showed a modest correlation (r² = 0.59, p = 0.01) After applying ICON corrections, dcCt.Th resulted in 0.52 ± 0.09 mm for QCT and 0.43 ± 0.07 mm for HR-QCT, both significantly correlated to HR-pQCT (r² = 0.75, p = 0.0012 and r² = 0.93, p < 0.0001, respectively). The average overestimation bias of Ct.Th was reduced from (402 ± 157)% to (45 ± 17)% for QCT and from (330 ± 69)% to (19 ± 8)% for HR-QCT. Due to inaccurate segmentation uncorrected QCT-based Ct.Th measures as well as BMD and BMC showed no correlation to HR-pQCT and thus such bias cortical data can be misleading. The application of ICON reduced random overestimation bias to about 50 μm and 20 μm for QCT and HR-QCT, respectively, leading to a recovery of a significant correlation with the reference data of HR-pQCT. This reveals the potential for fairly accurate assessment of cortical thickness, allowing to better characterize cortical mechanical competence. These results warrant testing of the performance in vivo.

Fundamental biomechanics of the spine--What we have learned in the past 25 years and future directions

Since the publication of the 2nd edition of White and Panjabi׳s textbook, Clinical Biomechanics of the Spine in 1990, there has been considerable research on the biomechanics of the spine. The focus of this manuscript will be to review what we have learned in regards to the fundamentals of spine biomechanics. Topics addressed include the whole spine, the functional spinal unit, and the individual components of the spine (e.g. vertebra, intervertebral disc, spinal ligaments). In these broad categories, our understanding in 1990 is reviewed and the important knowledge or understanding gained through the subsequent 25 years of research is highlighted. Areas where our knowledge is lacking helps to identify promising topics for future research. In this manuscript, as in the White and Panjabi textbook, the emphasis is on experimental research using human material, either in vivo or in vitro. The insights gained from mathematical models and animal experimentation are included where other data are not available. This review is intended to celebrate the substantial gains that have been made in the field over these past 25 years and also to identify future research directions.

Comparative role of disc degeneration and ligament failure on functional mechanics of the lumbar spine

Understanding spinal kinematics is essential for distinguishing between pathological conditions of spine disorders, which ultimately lead to low back pain. It is of high importance to understand how changes in mechanical properties affect the response of the lumbar spine, specifically in an effort to differentiate those associated with disc degeneration from ligamentous changes, allowing for more precise treatment strategies. To do this, the goals of this study were twofold: (1) develop and validate a finite element (FE) model of the lumbar spine and (2) systematically alter the properties of the intervertebral disc and ligaments to define respective roles in functional mechanics. A three-dimensional non-linear FE model of the lumbar spine (L3-sacrum) was developed and validated for pure moment bending. Disc degeneration and sequential ligament failure were modelled. Intersegmental range of motion (ROM) and bending stiffness were measured. The prediction of the FE model to moment loading in all three planes of bending showed very good agreement, where global and intersegmental ROM and bending stiffness of the model fell within one standard deviation of the in vitro results. Degeneration decreased ROM for all directions. Stiffness increased for all directions except axial rotation, where it initially increased then decreased for moderate and severe degeneration, respectively. Incremental ligament failure produced increased ROM and decreased stiffness. This effect was much more pronounced for all directions except lateral bending, which is minimally impacted by ligaments. These results indicate that lateral bending may be more apt to detect the subtle changes associated with degeneration, without being masked by associated changes of surrounding stabilizing structures.

Theoretical effects of fully ductile versus fully brittle behaviors of bone tissue on the strength of the human proximal femur and vertebral body

The influence of the ductility of bone tissue on whole-bone strength represents a fundamental issue of multi-scale biomechanics. To gain insight, we performed a computational study of 16 human proximal femurs and 12 T9 vertebral bodies, comparing the whole-bone strength for the two hypothetical bounding cases of fully brittle versus fully ductile tissue-level failure behaviors, all other factors, including tissue-level elastic modulus and yield stress, held fixed. For each bone, a finite element model was generated (60-82 μm element size; up to 120 million elements) and was virtually loaded in habitual (stance for femur, compression for vertebra) and non-habitual (sideways fall, only for femur) loading modes. Using a geometrically and materially non-linear model, the tissue was assumed to be either fully brittle or fully ductile. We found that, under habitual loading, changing the tissue behavior from fully ductile to fully brittle reduced whole-bone strength by 38.3±2.4% (mean±SD) and 39.4±1.9% for the femur and vertebra, respectively (p=0.39 for site difference). These reductions were remarkably uniform across bones, but (for the femur) were greater for non-habitual (57.1±4.7%) than habitual loading (p<0.001). At overall structural failure, there was 5-10-fold less failed tissue for the fully brittle than fully ductile cases. These theoretical results suggest that the whole-bone strength of the proximal femur and vertebra can vary substantially between fully brittle and fully ductile tissue-level behaviors, an effect that is relatively insensitive to bone morphology but greater for non-habitual loading.

Of mice, rats and men: trabecular bone architecture in mammals scales to body mass with negative allometry

Body mass (BM) in mammal species spans over six orders of magnitude. Although trabecular bone contributes to the mechanical properties of bones, we know much less about how trabecular bone scales with BM than about how cortical bone scales with BM. We therefore conducted a meta-analysis of the existing literature to test in rodents, humans and other mammals, predicted scaling properties between BM and several trabecular parameters: bone volume fraction (BV/TV), trabecular number (Tb.N), trabecular thickness (Tb.Th), trabecular separation (Tb.Sp), connectivity density (ConnD) and degree of anisotropy (DA). Our results show that BV/TV and DA are independent of BM and that Tb.N, Tb.Th and Tb.Sp scale with negative allometry relative to BM. Rodents appear to have relatively thicker and fewer trabeculae than humans, and we propose it that is due to a minimum thickness threshold "imposed" on mechanically functional trabeculae. Consequently, rodents (mice and rats) and humans demonstrate two distinct mechanisms to achieve variations in BV/TV. Although Tb.Th variation is the main contributing factor for differences in BV/TV in humans, Tb.N variation is the main contributing factor for differences in BV/TV in rodents. Our results also demonstrate no correlation between Tb.N and Tb.Th within each taxon (mice, rats and humans). Since rodents are a common animal model for research on bone biomechanics, the evidence that trabecular bone parameters scale and correlate differently in rodents than in humans suggests that care should be applied when extrapolating bone biomechanical results from small animals to large-bodied humans.

Calibration of hyperelastic material properties of the human lumbar intervertebral disc under fast dynamic compressive loads

Under fast dynamic loading conditions (e.g. high-energy impact), the load rate dependency of the intervertebral disc (IVD) material properties may play a crucial role in the biomechanics of spinal trauma. However, most finite element models (FEM) of dynamic spinal trauma uses material properties derived from quasi-static experiments, thus neglecting this load rate dependency. The aim of this study was to identify hyperelastic material properties that ensure a more biofidelic simulation of the IVD under a fast dynamic compressive load. A hyperelastic material law based on a first-order Mooney-Rivlin formulation was implemented in a detailed FEM of a L2-L3 functional spinal unit (FSU) to represent the mechanical behavior of the IVD. Bony structures were modeled using an elasto-plastic Johnson-Cook material law that simulates bone fracture while ligaments were governed by a viscoelastic material law. To mimic experimental studies performed in fast dynamic compression, a compressive loading velocity of 1 m/s was applied to the superior half of L2, while the inferior half of L3 was fixed. An exploratory technique was used to simulate dynamic compression of the FSU using 34 sets of hyperelastic material constants randomly selected using an optimal Latin hypercube algorithm and a set of material constants derived from quasi-static experiments. Selection or rejection of the sets of material constants was based on compressive stiffness and failure parameters criteria measured experimentally. The two simulations performed with calibrated hyperelastic constants resulted in nonlinear load-displacement curves with compressive stiffness (7335 and 7079 N/mm), load (12,488 and 12,473 N), displacement (1.95 and 2.09 mm) and energy at failure (13.5 and 14.7 J) in agreement with experimental results (6551 ± 2017 N/mm, 12,411 ± 829 N, 2.1 ± 0.2 mm and 13.0 ± 1.5 J respectively). The fracture pattern and location also agreed with experimental results. The simulation performed with constants derived from quasi-static experiments showed a failure energy (13.2 J) and a fracture pattern and location in agreement with experimental results, but a compressive stiffness (1580 N/mm), a failure load (5976 N) and a displacement to failure (4.8 mm) outside the experimental corridors. The proposed method offers an innovative way to calibrate the hyperelastic material properties of the IVD and to offer a more realistic simulation of the FSU in fast dynamic compression.

Foam phantom development for artificial vertebrae used for surgical training

Currently the surgical training of kyphoplasty and vertebroplasty is performed on patients or specimens. To improve patient safety, a project was initiated to develop an Augmented Reality simulator for the surgical training of these interventions. Artificial vertebral segments should be integrated to provide realistic haptic feedback. To reach this, resulting forces during needle insertions (trans- and extrapedicular) into formalin-fixed vertebral specimens were measured. The same insertion procedure was also performed on six customized polyurethane blocks with varying mechanical parameters. Based on the results of these measurements, a specific foam phantom was generated and the insertion force measured. Additionally a parametric model for the needle insertion into bone was designed calculating three characteristic parameters for all insertion measurements. The resulting insertion force for the foam phantom was comparable to the specimen measurements and the parametric model provided comprehensible characteristic parameters. Based on the resulting force during needle insertion into human vertebrae, a possible foam recipe for manufacturing artificial segments was found. Furthermore, the parametric model provides characteristic parameters for the assessment of phantoms as well as the development of its production process.

Biomechanics of vertebral compression fractures and clinical application

Local biomechanical factors in the etiology of vertebral compression fractures are reviewed. The vertebral body is particularly vulnerable to compression fracture when its bone mineral density (BMD) falls with age. However, the risk of fracture, and the type of fracture produced, does not depend simply on BMD. Equally important is the state of degeneration of the adjacent intervertebral discs, which largely determines how compressive forces are distributed over the vertebral body. Disc height also influences load-sharing between the vertebral body and neural arch, and hence by Wolff's Law can influence regional variations in trabecular density within the vertebral body. Vertebral deformity is not entirely attributable to trauma: it can result from the gradual accumulation of fatigue damage, and can progress by a quasi-continuous process of "creep". Cement injection techniques such as vertebroplasty and kyphoplasty are valuable in the treatment of these fractures. Both techniques can stiffen a fractured vertebral body, and kyphoplasty may contribute towards restoring its height. The presence of cement can limit endplate deformation, and thereby partially reverse the adverse changes in load-sharing which follow vertebral fracture. Cement also reduces time-dependent "creep" deformation of damaged vertebrae.

Voluntary exercise has long-term in vivo protective effects on osteocyte viability and bone strength following ovariectomy

Osteocytes are recognized as having a pivotal role in bone tissue homeostasis, and stimuli that increase osteocyte death result in decreased bone tissue quality. Previous in vitro studies have shown that mechanical stimulation prevents osteocyte death; however, in vivo evidence of this protective effect is limited. The aim of this study was to investigate if mechanical stimulation provided by voluntary exercise reduces osteocyte death caused by estrogen deficiency. Thirty-two female Wistar rats (5 months old) were either sacrificed as baseline controls (BSL, n = 7), ovariectomized or sham-operated and housed in cages with a voluntary running wheel (OVXEX, n = 7; SHAMEX, n = 6), or ovariectomized or sham-operated and housed in standard cages of equivalent size (OVXSED, n = 6; SHAMSED, n = 6) and sacrificed at age 14 months. Histomorphometric analysis of femur mid-diaphysis cortical bone revealed a significantly higher osteocyte number (N.Ot) and lower empty lacunae number (N.Lc) in both the OVXEX and SHAMEX groups compared to their SED counterparts. Intracortical porosity (Po.Ar) was also lower in both EX groups compared to their SED counterparts and significantly correlated with N.Lc (r = 0.616; P < 0.001). Three-point bending testing showed a significantly higher Young's modulus and ultimate stress in OVXEX compared to OVXSED and significant correlations between N.Lc and both yield stress (r = -0.376, P < 0.05) and ultimate stress (r = -0.369, P < 0.05) and between intracortical porosity and bone ultimate stress (r = -0.451, P < 0.05). Our results show that voluntary exercise prevented osteocyte death and that this protective effect was associated with increases in femur ultimate stress, which could be partially explained by decreases in Po.Ar.

Biomechanical Influence of Disk Properties on the Load Transfer of Healthy and Degenerated Disks Using a Poroelastic Finite Element Model

Spine degeneration is a pathology that will affect 80% of the population. Since the intervertebral disks play an important role in transmitting loads through the spine, the aim of this study was to evaluate the biomechanical impact of disk properties on the load carried by healthy (Thompson grade I) and degenerated (Thompson grades III and IV) disks. A three-dimensional parametric poroelastic finite element model of the L4/L5 motion segment was developed. Grade I, grade II, and grade IV disks were modeled by altering the biomechanical properties of both the annulus and nucleus. Models were validated using published creep experiments, in which a constant compressive axial stress of 0.35 MPa was applied for 4 h. Pore pressure (PP) and effective stress (SE) were analyzed as a function of time following loading application (1 min, 5 min, 45 min, 125 min, and 245 min) and discal region along the midsagittal profile for each disk grade. A design of experiments was further implemented to analyze the influence of six disk parameters (disk height (H), fiber proportion (%F), drained Young's modulus (Ea,En), and initial permeability (ka,kn) of both the annulus and nucleus) on load-sharing for disk grades I and IV. Simulations of grade I, grade III, and grade IV disks agreed well with the available published experimental data. Disk height (H) had a significant influence (p<0.05) on the PP and SE during the entire loading history for both healthy and degenerated disk models. Young’s modulus of the annulus (Ea) significantly affected not only SE in the annular region for both disk grades in the initial creep response but also SE in the nucleus zone for degenerated disks with further creep response. The nucleus and annulus permeabilities had a significant influence on the PP distribution for both disk grades, but this effect occurred at earlier stages of loading for degenerated than for healthy disk models. This is the first study that investigates the biomechanical influence of both geometrical and material disk properties on the load transfer of healthy and degenerated disks. Disk height is a significant parameter for both healthy and degenerated disks during the entire loading. Changes in the annulus stiffness, as well as in the annulus and nucleus permeability, control load-sharing in different ways for healthy and degenerated disks.

The contribution of trabecular bone to the stiffness and strength of rat lumbar vertebrae

STUDY DESIGN

In vitro compressive load-displacement experiments on intact rat lumbar vertebrae and on the same vertebrae after part of their trabecular bone was removed.

OBJECTIVE

To determine the contribution of the trabecular bone component to the stiffness and strength of rat lumbar vertebrae.

SUMMARY OF BACKGROUND DATA

Vertebral fractures are common in the aging population, possibly resulting from the deterioration of the mechanical properties of vertebral bone. Studies of the contribution of trabecular bone to the mechanical behavior of whole vertebra were published, but yielded mixed results. Here, we propose a novel optical metrology approach to address this important question.

METHODS

The bodies of intact rat lumbar vertebrae and the bodies of the same vertebrae after part of their trabecular bone was removed were loaded within their elastic region in a wet environment. The amount of trabecular bone removed was determined by micro-computer tomography scanning. Deformation maps of the dorsal vertebral surface of the intact and manipulated vertebrae were obtained using an optical metrology method, and compared. Intact and manipulated vertebrae were also loaded to failure in compression and their strengths and stiffness were compared.

RESULTS

The preferred trabecular orientation was found to be along the anterior-posterior axis, which is similar to humans. Removal of up to 42% of the trabecular tissue in the intact vertebrae did not significantly affect lumbar vertebral stiffness. However, removal of even smaller amounts of the intact trabecular tissue significantly reduced vertebral strength.

CONCLUSION

Trabeculae in rat lumbar vertebrae fulfill an important role in failure resistance (strength), but have little or no effect on the deformational behavior (stiffness) of the bone. These results differ from previous results we reported for rat femora, where removal of trabecular bone surprisingly increased the stiffness of the whole bone, and suggest that trabecular tissue may have different functions depending on anatomic location, bone function and morphology, and mode of loading.

Increase in vertebral body size in postmenopausal women with osteoporosis

Bone geometry plays a prominent role in bone strength. Cross-sectional studies have shown that advancing age is associated with increasing diameter of long bones, related to both periostal apposition and endosteal resorption. However, there are few data provided by prospective studies, especially concerning the changes in vertebral body dimensions. The objective of this prospective study was to measure the changes occurring in the vertebral body size of women with postmenopausal osteoporosis. Three-year data from placebo groups of the SOTI and TROPOS trials, performed in women with postmenopausal osteoporosis, were used for this study. In these trials, patients underwent lateral radiographs of the thoracic and lumbar spine at baseline and annually over 3 years, according to standardized procedures. Six-point digitization method was used: the four corner points of the vertebral body from T4 to L4 are marked, as well as an additional point in the middle of the upper and lower endplates. From these 6 points, the vertebral body perimeter, area and depth were measured at baseline and at 3 years. The analysis excluded all vertebrae with prevalent or incident fracture. A total of 2017 postmenopausal women (mean age 73.4+/-6.1 years) with a mean lumbar spine T score of -3.1+/-1.5, and a mean femoral neck T score of -3.0+/-0.7 are included in the analysis. Vertebral body dimensions increased over 3 years, by 2.1+/-5.5% (mean depth+/-SD), by 1.7+/-8.3% (mean area+/-SD) and by 1.5+/-4.9% (mean perimeter+/-SD) at the thoracic level (T4 to T12). At the lumbar level (L1 to L4), these dimensions increased as well: 1.4+/-3.6% (mean depth+/-SD), 1.4+/-5.7% (mean area+/-SD), 0.7+/-2.9% (mean perimeter+/-SD). A significant increase in vertebral body size was observed for each vertebral level from T5 to L4 for each of these parameters (p<0.01). These prospective results demonstrate that vertebral body dimensions increase over 3 years in women with postmenopausal osteoporosis.

Translational aspects of bone quality--vertebral fractures, cortical shell, microdamage and glycation: a tribute to Pierre D. Delmas

Among vertebral deformities, the prevalence of wedge fractures is about twice that of endplate (biconcave) deformities, both of which are greater than that of crush deformities. The anterior cortex is, therefore, a site of interest for understanding mechanisms of vertebral fracture. Despite its importance to vertebral mechanics, there are limited data describing the role of cortical shell, microdamage, and bone matrix parameters in vertebral fragility. This review of literature emphasizes the translational aspects of bone quality and demonstrates that a greater understanding of bone fractures will be gained through bone quality parameters related to both cortical and cancellous compartments as well as from microdamage and bone matrix parameters. In the context of vertebral fractures, measures of cortical shell and bone matrix parameters related to the organic matrix (advanced glycation products and alpha/beta CTX ratio) are independent of BMD measurements and can therefore provide an additional estimate of fracture risk in older patients.

A comparison of fatigue failure responses of old versus middle-aged lumbar motion segments in simulated flexed lifting

STUDY DESIGN

Survival analysis techniques were used to compare the fatigue failure responses of elderly motion segments to a middle-aged sample.

OBJECTIVES

To compare fatigue life of a middle-aged sample of lumbosacral motion segments to a previously tested elderly cohort. An additional objective was to evaluate the influence of bone mineral content on cycles to failure.

SUMMARY OF BACKGROUND DATA

A previous investigation evaluated fatigue failure responses of 36 elderly lumbosacral motion segments (average age, 81 +/- 8 years) subjected to spinal loads estimated when lifting a 9-kg load in 3 torso flexion angles (0 degrees, 22.5 degrees, and 45 degrees). Results demonstrated rapid fatigue failure with increased torso flexion; however, a key limitation of this study was the old age of the specimens.

METHODS

Each lumbosacral spine was dissected into 3 motion segments (L1-L2, L3-L4, and L5-S1). Motion segments within each spine were randomly assigned to a spinal loading condition corresponding to lifting 9 kg in 3 torso flexion angles (0 degrees, 22.5 degrees, or 45 degrees). Motion segments were statically loaded and allowed to creep for 15 minutes, then cyclically loaded at 0.33 Hz. Fatigue life was taken as the number of cycles to failure (10 mm displacement after creep loading).

RESULTS

Compared with the older sample of spines, the middle-aged sample exhibited increased fatigue life (cycles to failure) in all the torso flexion conditions. Increased fatigue life of the middle-aged specimens was associated with the increased bone mineral content (BMC) in younger motion segments (mean +/- SD, 30.7 +/- 11.1 g per motion segment vs. 27.8 +/- 9.4 g). Increasing bone mineral content had a protective influence with each additional gram increasing survival times by approximately 12%.

CONCLUSION

Younger motion segments survive considerably longer when exposed to similar spine loading conditions that simulate repetitive lifting in neutral and flexed torso postures, primarily associated with the increased bone mineral content possessed by younger motion segments. Cycles to failure of young specimens at 22.5 degrees flexion were similar to that of older specimens at 0 degrees flexion, and survivorship of young specimens at 45 degrees flexion was similar to the older cohort at 22.5 degrees.