Effects of suppression of bone turnover on cortical and trabecular load sharing in the canine vertebral body

The relative biomechanical effects of antiresorptive treatment on cortical thickness vs. trabecular bone microarchitecture in the spine are not well understood. To address this, T-10 vertebral bodies were analyzed from skeletally mature female beagle dogs that had been treated with oral saline (n=8 control) or a high dose of oral risedronate (0.5 mg/kg/day, n=9 RIS-suppressed) for 1 year. Two linearly elastic finite element models (36-μm voxel size) were generated for each vertebral body—a whole-vertebra model and a trabecular-compartment model—and subjected to uniform compressive loading. Tissue-level material properties were kept constant to isolate the effects of changes in microstructure alone. Suppression of bone turnover resulted in increased stiffness of the whole vertebra (20.9%, p=0.02) and the trabecular compartment (26.0%, p=0.01), while the computed stiffness of the cortical shell (difference between whole-vertebra and trabecular-compartment stiffnesses, 11.7%, p=0.15) was statistically unaltered. Regression analyses indicated subtle but significant changes in the relative structural roles of the cortical shell and the trabecular compartment. Despite higher average cortical shell thickness in RIS-suppressed vertebrae (23.1%, p=0.002), the maximum load taken by the shell for a given value of shell mass fraction was lower (p=0.005) for the RIS-suppressed group. Taken together, our results suggest that—in this canine model—the overall changes in the compressive stiffness of the vertebral body due to suppression of bone turnover were attributable more to the changes in the trabecular compartment than in the cortical shell. Such biomechanical studies provide an unique insight into higher-scale effects such as the biomechanical responses of the whole vertebra.     

Topography of the mineralization front on different surfaces of the body of a vertebra

Peculiarities of the relief of the mineralization front have been investigated on the periostal surface of the human vertebral body in several zones. The material has been obtained from male persons at the age of 20-84 years and studied by means of the light and scanning electron microscopy. The size of lateral surfaces of the vertebral body does not differ essentially from the relief of the periostal surfaces of other bones. In people of middle age certain changes in structure of the mineralized cartilage plate in the area adjoining the nucleus pulposus and the fibrous ring of the intervertebral disc are demonstrated. In persons of elderly and old age a definite decrease in thickness of the cortical layer of the vertebra is noted. At the same time, the plate of the mineralized cartilage adjoining the cortical layer grows thicker and collagene fibers in the spinal column ligaments undergo mineralization. Sometimes, microfractures of the cortical plate of the vertebral body and Schmorl noduli are revealed.     

Mid-sagittal dimensions of cervical vertebral bodies.

A series of lateral radiographs of the cervical spinal column was evaluated in order to determine vertebral body dimensions. The sample included males (N=30) and females (N=31) 18 to 24 years old, comprising three stature percentile ranges (1-20; 40-60; 80-99) of the U.S. adult population. A two-dimensional analysis of vertebral body height (average distance between superior-inferior surgaces), depth (average distance between anteriorposterior surfaces), and area (average height X average depth) revealed minimal effects due to stature. In all subjects, average depth exceeded average height for vertebral bodies C3 through C7. Upon combining stature groups, both sexes revealed maximum average values for these dimensions at the seventh cervical vertebral body. Minimum average height occurred at C5 whereas minimum average depth was found at C3. Significant correlation (alpha greater than 0.05) was found for males between ponderal index and height and depth of the C7 vertebra. Male head weight correlated significantly with C3, C4, C5 and C6 vertebral body height and with C3, C5 and C6 vertebral body depth. For females, C7 height and C6 depth correlated significantly with ponderal index and head weight respectively. Probable biomechanical relationships of specific cervical vertebral bodies are noted     

Vertebral morphology of Nacholapithecus kerioi based on KNM-BG 35250

This paper describes the morphology of the vertebral remains of the KNM-BG 35250 Nacholapithecus kerioi individual from the Middle Miocene of Kenya. Cervical vertebrae are generally large relative to presumed body mass, suggesting a heavy head with large jaws and well-developed neck muscles. The atlas retains the lateral and posterior bridges over the vertebral artery. The axis has a robust dens and a large angle formed by superior articular surfaces. The thoracic vertebral specimens include the diaphragmatic vertebra and one post-diaphragmatic vertebra. The thoracic vertebral bodies are much smaller that those of male Papio cynocephalus, whereas many of the dorsal elements are large and robust, exceeding those of male P. cynocephalus. Lumbar vertebral bodies are small relative to body mass, craniocaudally moderately long, and have a median ventral keel. The transverse process is craniocaudally long and arises from the widest part of the body cranially and the pedicle above the inferior vertebral notch caudally. Anapophyses are present in one of the preserved lumbar vertebrae. The postzygapophyses are thick dorsoventrally. These lumbar features are broadly shared with Proconsul. However, the base of the spinous process is longer and more caudally positioned in N. kerioi compared to Proconsul, and is more similar to the condition in Pongo. They are not dorsally (or moderately caudally) directed as is seen in P. nyanzae, Pan, and most other extant primates. A caudally directed spinous process does not permit a broad range of spinal dorsiflexion. The presumed stiff back in N. kerioi suggests a different locomotor repertoire than in Proconsul. Morotopithecus bishopi, although not possessing the same features, exhibits another morphological suite of characters for lumbar stiffness. Diverse functional adaptations of the lumbar spine were present in African hominoids during the Early to Middle Miocene.     

Medieval trabecular bone architecture: the influence of age, sex, and lifestyle

Osteoporosis has become a growing health concern in developed countries and an extensive area of research in skeletal biology. Despite numerous paleopathological studies of bone mass, few studies have measured bone quality in past populations. In order to examine age- and sex-related changes in one aspect of bone quality in the past, a study was made of trabecular bone architecture in a British medieval skeletal sample. X-ray images of 5-mm-thick coronal lumbar vertebral bone sections were taken from a total of 54 adult individuals divided into three age categories (18-29, 30-49, and 50+ years), and examined using image analysis to evaluate parameters related to trabecular bone structure and connectivity. Significant age-related changes in trabecular bone structure (trabecular bone volume (BV/TV), trabecular number (Tb.N), trabecular separation (Tb.Sp), and anisotropic ratio (Tb.An)) were observed to occur primarily by middle age with significant differences between the youngest and two older age groups. Neither sex showed continuing change in trabecular structure between the middle and old age groups. Age-related changes in bone connectivity (number of nodes (N.Nd) and node-to-node strut length (Nd.Nd)) similarly indicated a change in bone connectivity only between the youngest and two older age groups. However, females showed no statistical differences among the age groups in bone connectivity. These patterns of trabecular bone loss and fragility contrast with those generally found in modern populations that typically report continuing loss of bone structure and connectivity between middle and old age, and suggest greater loss in females. The patterns of bone loss in the archaeological samples must be interpreted cautiously. We speculate that while nutritional factors may have initiated some bone loss in both sexes, physical activity could have conserved bone architecture in old age in both sexes, and reproductive factors such as high parity and extended periods of lactation could have played a key role in female bone maintenance in this historic population. The study of qualitative elements (such as trabecular architecture) is vital if we are to understand bone maintenance and fragility in the past.     

Parametric finite element analysis of vertebral bodies affected by tumors

The vertebral column is the most frequent site of metastatic involvement of the skeleton. Due to the proximity to the spinal cord, from 5% to 10% of all cancer patients develop neurologic manifestations. As a consequence, fracture risk prediction has significant clinical importance. In this study, we model the metastatically involved vertebra so as to parametrically investigate the effects of tumor size, material properties and compressive loading rate on vertebral strength. A two-dimensional axisymmetric finite element model of a spinal motion segment consisting of the first lumbar vertebral body (no posterior elements) and adjacent intervertebral disc was developed to allow the inclusion of a centrally located tumor in the vertebral body. After evaluating elastic, mixed, and poroelastic formulations, we concluded that the poroelastic representation was most suitable for modeling the metastatically involved vertebra's response to compressive load. Maximum principal strains were used to localize regions of potential vertebral trabecular bone failure. Radial and axial vertebral body displacements were used as relative indicators of spinal canal encroachment and endplate failure. Increased tumor size and loading rate, and reduced trabecular bone density all elevated axial and radial displacements and maximum tensile strains. The results of this parametric study suggest that vertebral tumor size and bone density contribute significantly to a patients risk for vertebral fracture and should be incorporated in clinical assessment paradigms.     

Mechanical loading of mouse caudal vertebrae increases trabecular and cortical bone mass-dependence on dose and genotype

Most in vivo studies addressing the skeletal responses of mice to mechanical loading have targeted cortical bone. To investigate trabecular bone responses also we have developed a caudal vertebral axial compression device (CVAD) that transmits mechanical loads to compress the fifth caudal vertebra via stainless steel pins inserted into the forth and sixth caudal vertebral bodies. Here, we used the CVAD in C57BL/6 (B6) and C3H/Hej (C3H) female mice (15 weeks of age) to investigate whether the effect of regular bouts of mechanical stimulation on bone modeling and bone mass was dependent on dose and genotype. A combined micro-computed tomographic and dynamic histomorphometric analysis was carried out at the end of a 4-week loading regimen (3,000 cycles, 10 Hz, 3 x week) for load amplitudes of 0N, 2N, 4N and 8N. Significant increases in trabecular bone mass of 9 and 21% for loads of 4N and 8N, respectively, were observed in B6 mice. A significant increase of 10% in trabecular bone mass occurred for a load of 8N in the C3H strain. For other loads, no significant increases were detected. Both mouse strains exhibited substantial increases in trabecular bone formation rates for all loads, B6: 111% (2N), 86% (4N), 164% (8N), C3H: 41% (2N), 38% (4N), 141% (8N). Significant decreases in osteoclast number of 146 and 93% for a load of 8N were detected in B6 and C3H mice, respectively. These findings demonstrate that the effect of loading on the structural and functional parameters of bone is dose and genotype dependent. The caudal vertebral loading model established here is proposed for further studies addressing the molecular processes involved in the skeletal responses to mechanical stimuli.     

Examination of continuum and micro-structural properties of human vertebral cancellous bone using combined cellular solid models

Two- and three-dimensional structural models of the vertebral body have been used to estimate the mechanical importance of parameters that are difficult to quantify experimentally such as lattice disorder, trabecular thickness, trabecular spacing, connectivity, and fabric. Many of the models that investigate structure–function relationships of the vertebral body focus only on the trabecular architecture and neglect solid–fluid interactions. We developed a cellular solid model composed of two idealized unit cell geometries to investigate the continuum and micro-structural properties of human vertebral cancellous bone in a mathematically tractable model. Using existing histomorphological data we developed structure–function relationships for the mechanical properties of the solid phase, estimated the micro-structural strains, and predicted the fluid flow characteristics. We found that the micro-structural strains may be 1.7 to 2.2 times higher than the continuum level strains between the ages of 40 and 80. In addition, the predicted permeability agrees well with the experimental data.     

Comparison of trabecular bone behavior in core and whole bone samples using high-resolution modeling of a vertebral body

Computational analysis of trabecular bone normally involves the modeling of (experimental tests of) cored samples. However, the lack of constraint on the sides of the extracted trabecular bone samples limits the information that can be inferred regarding true in situ behavior. Here, the element-by-element voxel-based finite element method was applied via, a custom-written software suite (FEEBE), to a 72 μm resolution model of an ovine vertebra. The difference between the apparent modulus of eight concentric core cylinders when modeled as part of the whole bone (containing 84 × 10⁶ degrees of freedom) and independent of the whole bone was investigated. The results showed that cored trabecular bone apparent modulus depended significantly on the core diameter when modeled as an extracted core (r ² = 0.975) and as part of a whole bone (r ² = 0.986). The cause of this result was separated into the side-artifact effect and bone volume fraction (BV/TV) effect. For the independently modeled cores, the apparent modulus of an inner core region of interest varied with increasing thickness of the outer annulus. This was attributed to the side-artifact effect, given that the BV/TV of the core region was constant. Within the whole trabecular structure, the side artifact was eliminated as the entire bone structure was modeled. However, a BV/TV effect influenced the apparent modulus depending on the size of the core selected for determining apparent modulus. Changing the size of the core varied the overall BV/TV of the core, and this significantly (r ² = 0.999) influences the apparent modulus. Therefore, determining a ‘true’ apparent modulus for trabecular bone was not achievable. The independently modeled cores consistently under-predict the in vivo apparent modulus. It is recommended that if a ‘true’ apparent modulus is required, the BV/TV at which it is required needs to be first determined. Apparent modeling of entire bones at microscale resolution allowed regions of low and high tissue strains to be identified, consistent with patterns of trabecular bone remodeling and resorption reported in literature. The basivertebral vein cavity underwent the highest strains within the entire vertebral body, suggesting that failure might initiate here, despite containing visibly thicker struts and plate trabeculae. Although computationally expensive, analysis of the entire vertebral body provided a full picture of in situ trabecular bone deformation.     

Structural and mechanical adaptations of immature trabecular bone to strenuous exercise

Effects of strenuous exercise on immature bone were examined in two clinically important regions, femoral neck (FN) and lumbar vertebra (L6). Female Sprague-Dawley rats (n = 20, 8 wk of age, 150-170 g) were exercised progressively 5 days/wk for approximately 1 h/day for 10 wk at 75-80% of maximum oxygen capacity on a motor-driven treadmill. Caged age-matched rats served as controls (n = 20). Rat FNs were tested in cantilever bending, and vertebral bodies were compressed to 50% of their initial height at a fast strain rate. In response to the strenuous exercise, the relative area of the FN trabecular core increased significantly at the expense of the cortical shell. With that structural change, the exercised FN had significantly less energy to proportional limit than controls. The FN material properties (normal stresses at proportional limit and maximum) were significantly diminished after 10 wk of strenuous exercise. At the same time, no differences were found in vertebral geometry or structural and material properties. In the immature rate, the differential responses of the FN vs. L6 may relate to load history rather than a general systemic response to the strenuous exercise.     

Mechanisms of initial endplate failure in the human vertebral body

Endplate failure occurs frequently in osteoporotic vertebral fractures and may be related to the development of high tensile strain. To determine whether the highest tensile strains in the vertebra occur in the endplates, and whether such high tensile strains are associated with the material behavior of the intervertebral disc, we used micro-CT-based finite element analysis to assess tissue-level strains in 22 elderly human vertebrae (81.5±9.6 years) that were compressed through simulated intervertebral discs. In each vertebra, we compared the highest tensile and compressive strains across the different compartments: endplates, cortical shell, and trabecular bone. The influence of Poisson-type expansion of the disc on the results was determined by compressing the vertebrae a second time in which we suppressed the Poisson expansion. We found that the highest tensile strains occurred within the endplates whereas the highest compressive strains occurred within the trabecular bone. The ratio of strain to assumed tissue-level yield strain was the highest for the endplates, indicating that the endplates had the greatest risk of initial failure. Suppressing the Poisson expansion of the disc decreased the amount of highly tensile-strained tissue in the endplates by 79.4±11.3%. These results indicate that the endplates are at the greatest risk of initial failure due to the development of high tensile strains, and that such high tensile strains are associated with the Poisson expansion of the disc. We conclude that initial failure of the vertebra is associated with high tensile strains in the endplates, which in turn are influenced by the material behavior of the disc.     

Micromechanics of the human vertebral body for forward flexion

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

Contribution of vertebral [corrected] bodies, endplates, and intervertebral discs to the compression creep of spinal motion segments

Spinal segments show non-linear behavior under axial compression. It is unclear to what extent this behavior is attributable to the different components of the segment. In this study, we quantified the separate contributions of vertebral bodies and intervertebral discs to creep of a segment. Secondly, we investigated the contribution of bone and osteochondral endplate (endplates including cartilage) to the deformation of the vertebral body. From eight porcine spines a motion segment, a disc and a vertebral body were dissected and subjected to mechanical testing. In an additional test, cylindrical samples, machined from the lowest thoracic vertebrae of 11 porcine spines, were used to compare the deformation of vertebral bone and endplate. All specimens were subjected to three loading cycles, each comprising a loading phase (2.0 MPa, 15 min) and a recovery phase (0.001 MPa, 30 min). All specimens displayed substantial time-dependent height changes. Average creep was the largest in motion segments and smallest in vertebral bodies. Bone samples with endplates displayed substantially more creep than samples without. In the early phase, behavior of the vertebra was similar to that of the disc. Visco-elastic deformation of the endplate therefore appeared dominant. In the late creep phase, behavior of the segment was similar to that of isolated discs, suggesting that in this phase the disc dominated creep behavior, possibly by fluid flow from the nucleus. We conclude that creep deformation of vertebral bodies contributes substantially to creep of motion segments and that within a vertebral body endplates play a major role.     

Morphological variation of the thoracolumbar vertebrae in Macropodidae and its functional relevance

The present study was designed to investigate how the form of the marsupial thoracolumbar vertebrae varies to cope with the particular demands of diverse loading and locomotor behaviors. The vertebral columns of 10 species of Macropodidae, with various body masses and modes of locomotion, together with two other arboreal marsupials, koala and cuscus, were selected. Seventy-four three-dimensional landmark coordinates were acquired on each of the 10 last presacral vertebrae of the 70 vertebral columns. The interspecific variations of the third lumbar vertebra (L3, which approximates the mean) and the transitional patterns of the thoracolumbar segments were examined using the combined approach of generalized Procrustes analysis (GPA) and principal components analysis (PCA). The results of analyses of an individual vertebra (L3) and of the transitional patterns indicate significant interspecific differences. In the L3 study the first PC shows allometric shape variation, while the second PC seems to relate to adaptation for terrestrial versus arboreal locomotion. When the L3 vertebrae of the common spotted cuscus and koala are included for comparison, the vertebra of the tree kangaroo occupies an intermediate position between the hopping kangaroo and these arboreal marsupials. The L3 vertebrae in the arboreal marsupials possess a distinct dorsoventrally expanded vertebral body, and perpendicularly orientated spinous and transverse processes. The results of the present study suggest that vertebral shape in the kangaroo and wallaroos provides a structural adaptation to hopping through a relatively enlarged loading area and powerful lever system. In contrast, the small-sized bettongs (or rat kangaroos) have a relatively flexible column and elongated levers for the action of back muscles that extend and laterally flex the spine. The complex pattern of vertebral shape transition in the last 10 presacral vertebrae was examined using PCAs that compare between species information about vertebral shape variation along the thoracolumbar column. The results reinforce and emphasize important aspects of the patterns of variation seen in the detailed analysis of the third lumbar vertebra. The results also imply that size, spinal loading pattern, and locomotor behavior exert an influence on shaping the vertebra. Further, the morphological adaptations are consistent within these marsupials and this opens up the possibility that this kind of analysis may be useful in making functional inferences from fossil material.     

Biomechanical alterations in intact osteoporotic spine due to synthetic augmentation: finite element investigation

A three-dimensional nonlinear finite element model (FEM) was developed for a parametric study that examined the effect of synthetic augmentation on nonfractured vertebrae. The objective was to isolate those parameters primarily responsible for the effectiveness of the procedure; bone cement volume and bone density were expected to be highly important. Injection of bone cement was simulated in the FEM of a vertebral body that included a cellular model for the trabecular core. The addition of 10% and 20% cement by volume resulted in an increase in failure load, and the larger volume resulted in an increase in stiffness for the vertebral body. Placement of cement within the vertebral body was not as critical a parameter as cement amount. Simulated models of very poor bone quality saw the best therapeutic benefits.     

太行山猕猴第Ⅶ颈椎变量的异速生长分析

以肱骨最大长为参照,对成年太行山猕猴第Ⅶ颈椎变量进行了异速生长分析。结果表明,椎体后高、椎体上矢径、全宽呈正异速生长;椎体下横径接近等速生长;椎体前高、椎体下矢径、椎体上横径、椎孔矢径、椎孔横径、矢径呈负异速生长。     

Local distribution and dosimetry of 226Ra in the trabecular skeleton of the beagle.

Young adult beagle dogs received a single injection of 38.1 kBq/kg body wt 226Ra and were serially sacrificed at 4 to 2955 days postinjection. Samples of sites of trabecular bone in the lumbar vertebral body, proximal ulna, and distal femoral metaphysis and epiphysis were analyzed autoradiographically. The time-dependent changes in the average 226Ra concentrations in the four regions were analyzed in terms of a compartmental model. The clearance rate from the lumbar vertebral body was about four times more rapid than for the proximal ulna and distal femoral epiphysis. Ratios of hotspot to diffuse label concentrations varied from about 10 to 23. The dose rate to the endosteum ranged between 8.7 and 39.5 mGy/day initially and 4 and 10.5 mGy/day toward the end of the observation period. Mean marrow dose rates were lower by a factor of 3 to 9.5. During their residence time the nuclei of bone lining cells receive a maximum dose of 8 Gy in the proximal ulna (2955 days after injection) and a minimum dose of 0.63 Gy in the lumbar vertebra (2955 days after injection). This corresponds on the average to 17 and 1.4 alpha-particle hits to the cell nuclei, respectively.     

Rib fractures under anterior–posterior dynamic loads: Experimental and finite-element study

The purpose of this study was to investigate whether using a finite-element (FE) mesh composed entirely of hexahedral elements to model cortical and trabecular bone (all-hex model) would provide more accurate simulations than those with variable thickness shell elements for cortical bone and hexahedral elements for trabecular bone (hex–shell model) in the modeling human ribs. First, quasi-static non-injurious and dynamic injurious experiments were performed using the second, fourth, and tenth human thoracic ribs to record the structural behavior and fracture tolerance of individual ribs under anterior–posterior bending loads. Then, all-hex and hex–shell FE models for the three ribs were developed using an octree-based and multi-block hex meshing approach, respectively. Material properties of cortical bone were optimized using dynamic experimental data and the hex–shell model of the fourth rib and trabecular bone properties were taken from the literature. Overall, the reaction force–displacement relationship predicted by both all-hex and hex–shell models with nodes in the offset middle-cortical surfaces compared well with those measured experimentally for all the three ribs. With the exception of fracture locations, the predictions from all-hex and offset hex–shell models of the second and fourth ribs agreed better with experimental data than those from the tenth rib models in terms of reaction force at fracture (difference <15.4%), ultimate failure displacement and time (difference <7.3%), and cortical bone strains. The hex–shell models with shell nodes in outer cortical surfaces increased static reaction forces up to 16.6%, compared to offset hex–shell models. These results indicated that both all-hex and hex–shell modeling strategies were applicable for simulating rib responses and bone fractures for the loading conditions considered, but coarse hex–shell models with constant or variable shell thickness were more computationally efficient and therefore preferred.     

Cortical somatosensory potentials evoked by magnetic stimulation: effect of body height,age, and stimulus intensity

Cortical somatosensory evoked potentials elicited by non-invasive magnetic stimulation at Th10, Th12, and L5 vertebral levels and at mid-gluteus and ankle were studied in 45 normal controls varying in body height and age. Peak latencies of P2 evoked by ankle, mid-gluteus and L5 stimulation were correlated with body height. However, P2 latencies evoked by Th10 and Th12 stimulation were independent of body height. In contrast, all N2 latencies evoked by Th10, Th12, L5, mid-gluteus and ankle were correlated with body height. There were no correlations between latencies of P3 and N3 and body height. No peak latencies of cortical components were affected by age within the age group studied (15–67 years). In the intensity study, N2 showed double peaks when low to moderate stimulus intensities were applied to Th12 vertebra. Furthermore, in 2 patients with mild somatosensory dysfunction of lumbar roots or cord, N2 evoked by L5 stimulation revealed double peaks with prolonged interpeak latencies between the first peak of N2 and P3, which might be related to the temporal dispersion of the afferent fibers. These data suggest that the N2 component is induced by both fast and slow afferent fibers.