Three-dimensional photographic study of cancellous bone in human fourth lumbar vertebral bodies

Determination of anisotropic elastic parameters from morphological parameters of cancellous bone for osteoporotic lumbar spine

In biomechanics, large finite element models with macroscopic representation of several bones or joints are necessary to analyze implant failure mechanisms. In order to handle large simulation models of human bone, it is crucial to homogenize the trabecular structure regarding the mechanical behavior without losing information about the realistic material properties. Accordingly, morphology and fabric measurements of 60 vertebral cancellous bone samples from three osteoporotic lumbar spines were performed on the basis of X-ray microtomography (μCT) images to determine anisotropic elastic parameters as a function of bone density in the area of pedicle screw anchorage. The fabric tensor was mapped in cubic bone volumes by a 3D mean-intercept-length method. Fabric measurements resulted in a high degree of anisotropy (DA = 0.554). For the Young’s and shear moduli as a function of bone volume fraction (BV/TV, bone volume/total volume), an individually fit function was determined and high correlations were found (97.3 ≤ R² ≤ 99.1,p < 0.005). The results suggest that the mathematical formulation for the relationship between anisotropic elastic constants and BV/TV is applicable to current μCT data of cancellous bone in the osteoporotic lumbar spine. In combination with the obtained results and findings, the developed routine allows determination of elastic constants of osteoporotic lumbar spine. Based on this, the elastic constants determined using homogenization theory can enable efficient investigation of human bone using finite element analysis (FEA).

Graphical Abstract

The plate-to-rod transition in trabecular bone loss is elusive

Changes in trabecular micro-architecture are key to our understanding of osteoporosis. Previous work focusing on structure model index (SMI) measurements have concluded that disease progression entails a shift from plates to rods in trabecular bone, but SMI is heavily biased by bone volume fraction. As an alternative to SMI, we proposed the ellipsoid factor (EF) as a continuous measure of local trabecular shape between plate-like and rod-like extremes. We investigated the relationship between EF distributions, SMI and bone volume fraction of the trabecular geometry in a murine model of disuse osteoporosis as well as from human vertebrae of differing bone volume fraction. We observed a moderate shift in EF median (at later disease stages in mouse tibia) and EF mode (in the vertebral samples with low bone volume fraction) towards a more rod-like geometry, but not in EF maximum and minimum. These results support the notion that the plate to rod transition does not coincide with the onset of bone loss and is considerably more moderate, when it does occur, than SMI suggests. A variety of local shapes not straightforward to categorize as rod or plate exist in all our trabecular bone samples.

Maintaining Bone Health in the Lumbar Spine: Routine Activities Alone Are Not Enough

Public health organisations typically recommend a minimum amount of moderate intensity activities such as walking or cycling for two and a half hours a week, combined with some more demanding physical activity on at least 2 days a week to maintain a healthy musculoskeletal condition. For populations at risk of bone loss in the lumbar spine, these guidelines are particularly relevant. However, an understanding of how these different activities are influential in maintaining vertebral bone health is lacking. A predictive structural finite element modelling approach using a strain-driven algorithm was developed to study mechanical stimulus and bone adaptation in the lumbar spine under various physiological loading conditions. These loading conditions were obtained with a previously developed full-body musculoskeletal model for a range of daily living activities representative of a healthy lifestyle. Activities of interest for the simulations include moderate intensity activities involving limited spine movements in all directions such as, walking, stair ascent and descent, sitting down and standing up, and more demanding activities with large spine movements during reaching and lifting tasks. For a combination of moderate and more demanding activities, the finite element model predicted a trabecular and cortical bone architecture representative of a healthy vertebra. When more demanding activities were removed from the simulations, areas at risk of bone degradation were observed at all lumbar levels in the anterior part of the vertebral body, the transverse processes and the spinous process. Moderate intensity activities alone were found to be insufficient in providing a mechanical stimulus to prevent bone degradation. More demanding physical activities are essential to maintain bone health in the lumbar spine.

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.

Influences of endplate removal and bone mineral density on the biomechanical properties of lumbar spine

Purpose

To investigate (1) effects of endplate removal and bone mineral density (BMD) on biomechanical properties of lumbar vertebrae (2) whether the distributions of mechanical strength and stiffness of endplate are affected by BMD.

Methods

A total of thirty-one lumbar spines (L1-L5) collected from fresh cadavers were used in this study. Bone density was measured using lateral DEXA scans and parts of samples were performed with partial or entire endplate removal. All the specimens were divided into three BMD groups. According to endplate integrity of the lumbar vertebrae, each BMD group was then divided into three subgroups: subgroup A: intact endplate; subgroup B: central region of endplate removal; subgroup C: entire endplate removal. The axial compression test was conducted with material testing system at a speed of 2mm/min. The experimental results were statistically analyzed using SPSS 17.0.

Results

(1) Significant differences of biomechanical properties occurred among normal BMD, osteoporotic and serious osteoporotic group (P<0.05). (2) Spearman analysis showed that BMD was positively correlated with the failure load and stiffness of lumbar vertebrae. (3) For each BMD group, significant differences of biomechanical properties were found between subgroup A and C, and between subgroup B and C (P<0.05). (4) For each BMD group, there was no statistical difference of biomechanical properties between subgroup A and B (P>0.05).

Conclusions

Entire endplate removal can significantly decrease the structural properties of lumbar vertebrae with little change in biomechanical properties by preservation of peripheral region of the endplate. BMD is positively correlated to the structural properties of the lumbar vertebrae.

Age-related changes of vertical and horizontal lumbar vertebral trabecular 3D bone microstructure is different in women and men

The study presents a 3D method for subdividing a trabecular network into horizontal and vertical oriented bone. This method was used to investigate the age related changes of the bone volume fraction and thickness of horizontal and vertical trabeculae in human lumbar vertebral bone estimated with unbiased 3D methods in women and men over a large age-range. The study comprised second lumbar vertebral body bone samples from 40 women (aged 21.7-96.4years, median 56.6years) and 39 men (aged 22.6-94.6years, median 55.6years). The bone samples were μCT scanned and the 3D microstructure was quantified. A voxel based algorithm inspecting the local neighborhood is presented and used to segment the trabecular network into horizontal and vertical oriented bone. For both women and men BV/TV decreased significantly with age, Tb.Th* was independent of age, while SMI increased significantly with age. Vertical (BV.vert/TV) and horizontal (BV.horz/TV) bone volume fraction decreased significantly with age for both sexes. BV.vert/TV decreased significantly faster with age for women than for men. Vertical (Tb.Th*.vert) and horizontal (Tb.Th*.horz) trabecular thickness were independent of age, while Tb.Th*.horz/Tb.Th*.vert decreased significantly with age for both sexes. Additionally, the 95th percentile of the trabecular thickness distribution increased significantly with age for vertical trabeculae in women, whereas it was independent of age in men. In conclusion, we have shown that vertical and horizontal oriented bone density decreases with age in both women and men, and that vertical oriented bone is lost more quickly in women than in men. Furthermore, vertical and horizontal trabecular thickness were independent of age, whereas the horizontal to vertical trabecular thickness ratio decreased significantly with age indicating a relatively more pronounced thinning of horizontal trabeculae. Finally, the age-related loss of trabecular elements appeared to result in a compensatory hypertrophy of vertical trabeculae in women, but not in men.

Subsidence of polyetheretherketone intervertebral cages in minimally invasive lateral retroperitoneal transpsoas lumbar interbody fusion

STUDY DESIGN

A retrospective review.

OBJECTIVE

The objective is to evaluate subsidence related to minimally invasive lateral retroperitoneal lumbar interbody fusion by reviewing our experience with this procedure.

SUMMARY OF BACKGROUND DATA

Polyetheretherketone intervertebral cages of different lengths, widths, and heights filled with various allograft types are commonly used as spacers in lumbar fusions. Subsidence is a potential complication. To date, there are no published reports specifically addressing subsidence, because it relates to a series of patients undergoing minimally invasive lateral retroperitoneal transpsoas lumbar interbody fusion.

METHODS

An institutional review board-approved, retrospective review of a prospectively collected database was conducted. One hundred forty consecutive patients who underwent this procedure between L1 and L5 during a 2-year period were included. All patients had T scores of -2.5 or more. Postoperative radiographs during routine follow-ups were reviewed for subsidence, defined as any violation of the vertebral end plate.

RESULTS

Radiographical subsidence occurred in 14.3% (20 of 140), whereas clinical subsidence occurred in 2.1%. Subsidence occurred in 8.8% (21 of 238) of levels fused. Construct length had a significant positive correlation with increasing subsidence rates. Subsidence rates decreased progressively with lower levels in the lumbar spine, but had a higher than expected rate at L4-L5. Subsidence rates of 14.1% (19 of 135) and 1.9% (2 of 103) were associated with 18-and 22-mm-wide cages, respectively. No significant trends were observed with cage lengths. Supplemental lateral plates had a higher rate of subsidence than bilateral pedicle screws. Subsidence occurred at the superior end plate 70% of the time.

CONCLUSION

The use of wider intervertebral cages leads to a significantly lower rate of subsidence, but a longer cage does not necessarily offer a similar advantage. Wide cages are protective against subsidence, and the widest cages should be used whenever feasible for interbody fusion in the lumbar spine to protect indirect compression and promote arthrodesis.

Towards quantitative 3D imaging of the osteocyte lacuno-canalicular network

Osteocytes are the most abundant cells in bone and the only cells embedded in the bone mineral matrix. They form an extended, three-dimensional (3D) network, whose processes interconnecting the cell bodies reside in thin canals, the canaliculi. Together with the osteocyte lacunae, the canaliculi form the lacuno-canalicular network (LCN). As the negative imprint of the cellular network within bone tissue, the LCN morphology is considered to play a central role for bone mechanosensation and mechanotransduction. However, the LCN has neither been visualized nor quantified in an adequate way up to now. On this account, this article summarizes the current state of knowledge of the LCN morphology and then reviews different imaging methods regarding the quantitative 3D assessment of bone tissue in general and of the LCN in particular. These imaging methods will provide new insights in the field of bone mechanosensation and mechanotransduction and thus, into processes of strain sensation and transduction, which are tightly associated with osteocyte viability and bone quality.

Simulation of vertebral trabecular bone loss using voxel finite element analysis

Trabecular bone loss in human vertebral bone is characterised by thinning and eventual perforation of the horizontal trabeculae. Concurrently, vertical trabeculae are completely lost with no histological evidence of significant thinning. Such bone loss results in deterioration in apparent modulus and strength of the trabecular core. In this study, a voxel-based finite element program was used to model bone loss in three specimens of human vertebral trabecular bone. Three sets of analyses were completed. In Set 1, strain adaptive resorption was modelled, whereby elements which were subject to the lowest mechanical stimulus (principal strain) were removed. In Set 2, both strain adaptive and microdamage mechanisms of bone resorption were included. Perforation of vertical trabeculae occurred due to microdamage resorption of elements with strains that exceeded a damage threshold. This resulted in collapse of the trabecular network under compression loading for two of the specimens tested. In Set 3, the damage threshold strain was gradually increased as bone loss progressed, resulting in reduced levels of microdamage resorption. This mechanism resulted in trabecular architectures in which vertical trabeculae had been perforated and which exhibited similar apparent modulus properties compared to experimental values reported in the literature. Our results indicate that strain adaptive remodelling alone does not explain the deterioration in mechanical properties that have been observed experimentally. Our results also support the hypothesis that horizontal trabeculae are lost principally by strain adaptive resorption, while vertical trabeculae may be lost due to perforation from microdamage resorption followed by rapid strain adaptive resorption of the remaining unloaded trabeculae.

A critical damping approach for assessing the role of marrow fat on the mechanical strength of trabecular bone

Several clinical findings revealed that post-menopausal osteoporosis and age-related osteopenia are accompanied by trabecular bone marrow fat (BMF) increase. To help understand this phenomenon, a vibrating string model is proposed, based on the hypothesis that, when bone marrow properties change, the trabecular bone structure remodels itself to preserve its critical damping state. It is found that an inverse relationship holds between trabecular average length and marrow damping coefficient. Such a result leads us to hypothesize the following bone-weakening mechanism. Since fat-rich bone marrow is a worse damper, a BMF increment causes an increase of trabecular average length, which is accomplished by the absorption of horizontal trabeculae (structurally less important than vertical trabeculae). The resulting bone patterns are in excellent agreement with clinical observations of osteoporotic bone. A definitive confirmation of the proposed mechanism will support a therapeutical approach to widespread osteopenic diseases aimed at avoiding, or limiting, BMF increase.

Trabecular Bone Structure Obtained From Multislice Spiral Computed Tomography of the Calcaneus Predicts Osteoporotic Vertebral Deformities

PURPOSE

To compare multislice computed tomography (MSCT)-derived parameters of the trabecular bone structure of the calcaneus with bone mineral density (BMD) in their ability to differentiate between donors with and without osteoporotic fractures of the spine and to optimize CT scan protocols.

METHODS

Forty-two postmortem calcanei (81.2 +/- 10 years) were imaged with a 16-detector row MSCT system using 4 different scan protocols varying spatial resolution (12-24 lp/cm) and radiation dose. Structural parameters of trabecular bone were derived from these images, and BMDs of the calcanei were determined using dual x-ray absorptiometry. Vertebral deformities of the spine were radiographically classified using the Spinal Fracture Index. Diagnostic performance in differentiation between donors with and without vertebral fractures was assessed using receiver operating characteristic (ROC) analysis.

RESULTS

There were significant case-control differences for many of the structural parameters measured (P < 0.05). The highest ROC values were found for apparent trabecular thickness using the high-resolution and high-dose protocols. Statistically significant correlations were found between most structure parameters and BMD (up to r = 0.85, P < 0.01).

CONCLUSION

Structural parameters of trabecular bone as obtained from high-resolution MSCT images of the calcaneus can be used to differentiate between donors with and without osteoporotic vertebral fractures, using a high-resolution and high-dose CT protocol.

First inventory of resorption lacunae on rods and plates of trabecular bone as observed by scanning electron microscopy

In the present study a novel systematic distribution scheme of resorption lacunae (RL) was applied using scanning electron microscopy. RL, classified as either reticulate patch resorption lacunae (RPR) or as longitudinally extended resorption lacunae (LER) [11, were analyzed and quantified according to their localizations on rods (middle, nodes or both) and plates (central or peripheral) in standardized segments from the femoral head of 24 Caucasian subjects without bone disease. Age and gender variations were explored. No clear gender-related distribution pattern could be detected on plates. On rods of males, however, the distribution of RL tended to be higher at the nodes, but seemed to be more prevalent in the middle or extended from the middle to the nodes of rods in females. Certain other non-conclusive tendencies in relation to age, gender, type of RL and localization were observed.

An orientation distribution function for trabecular bone

We describe a new method for quantifying the orientation of trabecular bone from three-dimensional images. Trabecular lattices from five human vertebrae were decomposed into individual trabecular elements, and the orientation, mass, and thickness of each element were recorded. Continuous functions that described the total mass (M(phi,theta)) and mean thickness (tau(phi,theta)) of all trabeculae as a function of orientation were derived. The results were compared with experimental measurements of the elastic modulus in three principal anatomic directions. A power law scaling relationship between the anisotropies in mass and elastic modulus was observed; the scaling exponent was 1.41 (R2=0.88). As expected, the preponderance of trabecular mass was oriented along the cranial-caudal direction; on average, there was 3.4 times more mass oriented vertically than horizontally. Moreover, the vertical trabeculae were 30% thicker, on average, than the horizontal trabeculae. The vertical trabecular thickness was inversely related to connectivity (R2=0.70; P=0.07), suggesting a possible organization into either few, thick trabeculae or many thin trabeculae. The method, which accounts for the mechanical connectedness of the lattice, provides a rapid way to both visualize and quantify the three-dimensional organization of trabecular bone.

Trabecular rod thickness by direct measurement from 3D SEM anaglyphs

This study presents a methodology for measuring the thickness of trabecular rods directly from anaglyphs. Macerated sagittal slices of T12 vertebral bodies from 15 subjects were examined by scanning electron microscopy (SEM). Two digital images (the second image tilted 5 degrees ) were recorded, and a 3D anaglyph was created. The thickness of the trabecular rods (Tb.Th((rods))), and the anatomical orientation of the trabecular rods were measured using an image analyser. Conventional 2D histomorphometry was performed on adjacent bone slices. A total of 1559 rod measurements were made from the 15 vertebral bone slices, with a mean Tb.Th((rods)) of 123 +/- 36 microm. The rod thickness in males (128 +/- 34 microm) was significantly greater than that in females (119 +/- 37 microm, P < 0.001). Tb.Th((rods)) changed significantly with age in the males: the thicker rods in the younger men reduced with age to a thickness similar to that in women. 3D measurements were significantly larger than the 2D estimates, and there was no correlation between the two methods of measurement. An inverse correlation was found between the number of rods and the bone volume fraction (BV/TV), indicating that decreased BV/TV is associated with an increased number of rods. The vertical rods (132 +/- 39 microm) were significantly thicker than the horizontal rods (116 +/- 33 microm, P < 0.001). The determination of rod numbers, and their orientation and individual thicknesses enables a greater understanding of cancellous bone architecture in both individuals and populations, and will allow more reliable finite element modelling. Direct measurements from 3D anaglyphs of intact specimens provide new data that show previously unrecognised age- and sex-related changes.

The effects of bone density and disc degeneration on the structural property distributions in the lower lumbar vertebral endplates

In this study, we hypothesized that vertebral bone density and disc degeneration would affect the structural property distributions of the lower lumbar vertebral endplates (L3-L5). The results may have implications for improving interbody implant designs to better resist subsidence. A 3 mm diameter hemispherical indenter was used to perform indentation tests at 0.2 mm/s to a depth of 3 mm at 27 standardized locations in 55 bony endplates of intact human lumbar vertebrae (L3-L5). The resulting load-displacement curves were used to extract the failure load and stiffness of each test site. Bone density was measured using lateral DEXA scans. Disc condition was determined using a four-point grading scale. Three-way analyses of variance were used to analyze the relationships between the data. The overall failure load decreased with bone mineral density (BMD) in the superior (p < 0.0001) and inferior (p = 0.011) lumbar endplates. In both endplates, the posterolateral regions were significantly stronger than more central regions. With increasing BMD, this difference became more pronounced in the superior endplates only (p = 0.005). Increased disc degeneration was associated with an overall failure load decrease in the inferior lumbar endplates (p = 0.002). The strength in the central regions of the superior endplates was reduced with increasing degeneration, but this was not observed peripherally (p = 0.001). Stiffness magnitude or distribution was not significantly affected by BMD or disc degeneration. The locations of the strongest regions of the endplate did not change with either bone density or disc degeneration. This implies that implant shapes designed using the basic structural property maps for the L3-L5 endplates are appropriate for use in patients with a wide range of pathologies, even though overall failure loads are generally lower in patients with reduced bone density and greater degrees of disc degeneration.

Age-related differences between thinning of horizontal and vertical trabeculae in human lumbar bone as assessed by a new computerized method

To investigate whether vertical trabeculae undergo compensatory thickening with age in the human vertebral body, a new computerized method was developed that is able to distinguish between horizontal and vertical trabeculae on normal histological sections. Study subjects included 48 individuals (24 women aged 19-97 years, and 24 men aged 23-95 years). From each L-2, thick frontal sections of half of the vertebra were embedded undecalcified in methylmetacrylate and cut into 10-microm-thick sections. A simple method able to classify image pixels as belonging to either vertical or horizontal trabeculae was developed and implemented in a computer program. The parallel-plate model was modified so that it was able to determine trabecular thickness, number, and separation (Tb.Th, Tb.N, and Tb.Sp) for horizontal and vertical trabeculae separately. The histomorphometric parameters were measured in three different regions: whole section, mid-third, and sub-endplate, for both horizontal and vertical trabeculae. It was found that the trabecular thickness of vertical trabeculae was independent of age in all investigated regions. The thickness of the horizontal trabeculae, in contrast, decreased significantly with age in all these regions. Tb.N decreased significantly with age for both horizontal and vertical trabeculae in all regions. However, the relative loss of trabeculae per unit length was largest for the horizontal trabeculae, whereas the absolute loss of trabeculae per unit length was largest for the vertical trabeculae. Tb.Sp was found to increase significantly with age for both vertical and horizontal trabeculae in all regions. No significant gender-related differences were found. In conclusion, in this study we describe a new, simple method for separation of horizontal and vertical bone tissue. This method has been applied on frontal vertebral sections. Trabecular bone thickness measured with the parallel-plate model showed that the thickness of horizontal trabeculae decreases significantly with age, whereas the thickness of vertical trabeculae did not decrease significantly with age. Furthermore, although the relative loss of trabeculae was larger for horizontal trabeculae than for vertical trabeculae, the absolute loss of vertical trabeculae was higher than that of horizontal trabeculae.

An investigation of thoracic and lumbar cancellous vertebral architecture using power-spectral analysis of plain radiographs

The internal architecture of the vertebral bodies spanning the levels T1 to L5 in seven male columns was studied using mammographic-resolution radiographs of 2.5-mm-thick planar parasagittal slices. The overlapping radiographic shadows of vertebral trabeculae combined in the image to form a series of 'elements', broadly representative of the cancellous structure. The orientations and sizes of these elements were analysed by applying the Fast Fourier transform (FFT) to the digitized radiographic images. Elements aligned in the 'vertical' orientation, along the long axis of the column, were the most prominent for all vertebral levels. The relative prominence of horizontal to vertical elements was generally constant along the column below T5. In contrast, the relative prominence of oblique to vertical elements declined in the cranio-caudal direction, particularly in individuals aged > or = 60 years. The ratio of 'large' (x > 0.3 mm) to 'small' (0.15 mm < or = x < or = 0.3 mm) elements was unchanged cranio-caudally in specimens < 60 years. However, in individuals > or = 60 years, large elements increased in relative prominence in the caudal direction. These results suggest that a basic orthogonal pattern of trabeculae is found along the male human spine, regardless of differences in vertebral body size. Power-spectral analysis is shown to yield information summarizing the predominant orientations and sizes of radiographically rendered architectural elements of vertebral cancellous bone, to define the effects of ageing on architecture, and to identify broad structural differences between vertebral levels in the adult male spine.

Zone-dependent changes in human vertebral trabecular bone: clinical implications

We have previously shown that there are pronounced age-related changes in human vertebral cancellous bone density and microarchitecture. However, the magnitude of these changes seemed to be dependent on zone location in the vertebral body-the central third vs. the areas adjacent to the endplates. The aim of the present study was, therefore, to investigate whether such zone-specific differences could be identified by static histomorphometric measures. The material comprised 48 individuals (24 women aged 19-97 years, and 24 men aged 23-95 years). Three of the women had a known fracture of the L-2. From each L-2, thick frontal sections of half of the vertebra were embedded undecalcified in methylmethacrylate, cut into 10-microm-thick sections, and stained with aniline blue. The sections were scanned into a computer, and classic static histomorphometry was performed on the images. The histomorphometry was performed on both the whole section and on the separate zones (central and sub-endplate zone). The results showed that trabecular bone volume, trabecular number, and connectivity density decreased significantly faster with age, whereas marrow space star volume increased significantly faster with age in the zones adjacent to the endplates than in the central zone. The other histomorphometric measures showed no zone specificity in relation to aging. However, trabecular thickness and trabecular separation were both higher at all ages in the central zone than in the sub-endplate zone, although this was significant only for trabecular separation. The described differences might have significant clinical implications concerning quantitative computed tomography (QCT) scanning, X-ray analyses, and assessment of fracture liability in the human spine, but the underlying pathogenesis is still not known. This study shows that the human vertebral body can be described as two distinct zones with very specific age-related changes in density and microstructure. This zone-specificity is important for the correct interpretation of clinical data.

Morphologic detail of aging bone in human vertebrae

To investigate aging bone structure of humans--here, in the lumbar vertebral bodies-requires methodologies that have sufficiently high resolving power yet still have sufficient width and depth of field. No clinical imaging method can come close to meeting the first requirement, leading to the disadvantage of being limited to postmortem studies. Few microscopic methods meet the second and third requisites. The three-dimensional (3D) images of bone in this article were obtained using deep-field 3D optical imaging, X-ray imaging, and scanning electron microscopy (SEM) of macerated plane parallel slices. The study of bone as a 3D object provides a different perspective from conventional two-dimensional images, and enriches our understanding of how modeling and remodeling processes regulate bone structure and connectivity. The study of ultraflat block surfaces by quantitative back-scattered electron imaging permits acquisition of data on mineral distributions and densities within a very thin layer (a continuous and perfect very thin section) in the block face. With this information, bone can be viewed as a spectrum of tissue types varying in degree of mineralization.

Inhomogeneity of human vertebral cancellous bone: systematic density and structure patterns inside the vertebral body

In the spine, cancellous bone quality is usually assessed for the whole vertebral body in a transverse central slice. Correct identification and assessment of the weakest parts of the cancellous bone may lead to better prediction of fracture risk. The density and structural parameters were systematically investigated inside the thoracic (T-9), thoracolumbar (T12-L1), and lumbar (L-4) vertebral bodies of nine subjects. On both sides of the median sagittal plane, anterior and posterior 8.2 mm vertical cores were harvested in the thoracic vertebra. In the thoracolumbar and lumbar vertebrae, external samples were also cored. Peripheral quantitative computed tomographic (pQCT) density analysis of the 136 cores was performed at four different levels, from the lower to the upper endplate. The relatively thin slice thickness (300 microm) and small pixel size (70 microm x 70 microm) was considered sufficient to investigate the structural parameters on the four transverse slices and in the sagittal and coronal planes (total of 816 images). Using a constant threshold a binary image was generated and the morphometric data were extracted. The binary image was further skeletonized and classical strut analysis was performed. Cancellous bone density was 20% higher in the posterior cores than in the anterior and external cores. Moreover, clear vertical inhomogeneity was noted because the lowest half of the vertebral body presented lower density than the upper half (differences ranging from 25% to 15%). All structural parameters were strongly dependent on the location of the measurement. Structural differences between anterior, posterior, and external areas were mild and followed the density patterns. On the other hand, vertical inhomogeneity of the structural parameters was important. For example, in the thoracolumbar and lumbar vertebrae, the numbers of nodes or node-to-node struts were almost twofold higher in the inferior half than in the superior half (p < 0.01), whereas trabecular thickness and number of free-ends presented a center/close-to-endplate structural pattern, with central trabeculae being 15% thicker (p < 0.05) and presenting 30% fewer free-ends (p < 0.01) than the close-to-endplate ones. Variability of density and structural parameters was high and a substantial part of this variability could be explained by the place inside the vertebral body where the measurement was made. The weak part was not in the center of the body but in its upper half where the lower density did not seem to be compensated by a higher structural architecture. Further clinical investigation could enhance fracture prediction by tracking and focusing on the weakest part of the vertebral body.

Mapping the structural properties of the lumbosacral vertebral endplates

STUDY DESIGN

A biomechanical investigation using indentation tests in a human cadaveric model to seek variation in the structural properties across the lower lumbar and sacral endplates.

OBJECTIVES

To determine 1) if there are regional differences in endplate strength and 2) whether any differences identified are affected by spinal level (lumbar spine vs. sacrum) or endplate (superior vs. inferior).

SUMMARY OF BACKGROUND DATA

It has been postulated that some regions of the vertebral body may be stronger than others. Conclusive data, either supporting or disproving this theory, would be valuable for both spine surgeons and implant designers because one mode of failure of interbody implants is subsidence into one or both adjacent vertebrae.

METHODS

Indentation tests were performed at 27 standardized test sites in 62 bony endplates of intact human vertebrae (L3-S1) using a 3-mm-diameter, hemispherical indenter with a test rate of 0.2 mm/sec to a depth of 3 mm. The failure load and stiffness at each test site were determined using the load-displacement curves. Three-way analyses of variance were used to analyze the resulting data.

RESULTS

Both the failure load and stiffness varied significantly across the endplate surfaces (P < 0.0001), with posterolateral regions being stronger and stiffer than the central regions. Characteristic distributions were identified in the lumbar superior, lumbar inferior, and sacral endplates. The failure load distributions were found to differ in 1) the superior lumbar and sacral endplates (P = 0.0077), 2) the inferior lumbar and sacral endplates (P = 0.0014), and 3) the superior and inferior lumbar endplates (P < 0.0001). The sacral and inferior lumbar endplates were both found to be stronger than the superior lumbar endplates (sacrum, P = 0.054; inferior, P = 0.008) but were not themselves significantly different (P = 0.89).

CONCLUSIONS

Highly significant regional strength and stiffness variations were identified in the lumbar and sacral endplates. The center of the bone, where implants are currently placed, is the weakest part of the lumbar endplates and is not the strongest region of the sacral endplate.

Accuracy limits for the determination of cortical width and density: the influence of object size and CT imaging parameters

In this study we analysed the accuracy of computed tomography (CT) measurements in assessing cortical bone. We determined the dependency of thickness and density measurements on the true width and density of the cortex and on the spatial resolution in the CT images using two optimized segmentation methods. As a secondary goal, we assessed the ability of CT to reflect small changes in cortical thickness. Two different bone-mimicking phantoms with varying cortical thickness were scanned with single-slice CT on a Somatom Plus 4 scanner. Images were reconstructed with both a standard and a high-resolution convolution kernel. Two special operator-independent segmentation methods were used to automatically detect the edges of the cortical shell. We measured cortical thickness and density and compared the phantom measurements with theoretical computations by simulating a cross-sectional shape of the cortical shell. Based on the simulations, we calculated CT's power to detect small changes in cortical thickness. Simulations and phantom measurements were in very good agreement. Cortical thickness could be measured with an error of less than 10% if the true thickness was larger than 0.9 (0.7) mm for the standard (high-resolution) kernel which is close to the full width at half maximum (FWHM) of the point spread functions for these kernels and our scanner. Density measurements yielded errors of less than 10% for true cortical thickness values above two to three times the FWHM corresponding to 2.5 (2) mm in our case. The simulations showed that a 10% change in cortical width would not be detected with satisfying probability in bones with a cortical shell thinner than 1.2 mm. An accurate determination of the cortical thickness is limited to bones with a thickness higher than the FWHM of the scanner's point spread function. Therefore, the use of a high-resolution reconstruction kernel is crucial. Cortical bone mineral density can only be measured accurately in bones two to three times thicker than this number. In thinner bones, the measured density becomes dependent on the thickness. Changes in cortical thickness can only be assessed if the change is rather large or if the measured bone has sufficient thickness. Therefore, assessing density or thickness of the vertebral shell by CT should be treated with caution.

Effect of estrogen suppression on the mineralization density of iliac crest biopsies in young women as assessed by backscattered electron imaging

The effects of estrogen suppression on bone mineralization in young women were studied by quantitative backscattered electron (BSE) imaging of transiliac biopsies taken before and after treatment for endometriosis. Treatment (6 months) was with analogs of gonadotrophin releasing hormone (GnRH) given either alone (six paired biopsies), which resulted in a marked reduction in the levels of circulating estrogen, or in conjunction with tibolone, a synthetic steroid with estrogenic, progestrogenic, and androgenic properties (four paired biopsies). Estrogen withdrawal increased (p < 0.01) and concomitant tibolone treatment decreased (p < 0.05) the overall mean bone density. Estrogen withdrawal increased the fraction of bone with a high mineralization density [pretreatment: 0.236+/-0.007; GnRH: 0.279+/-0.009, mean +/- standard error of the mean (SEM); p < 0.01]. The concomitant addition of tibolone reversed these effects and increased the proportion of bone with a low mineralization density (pretreatment: 0.198 +/- 0.005; tibolone: 0.230 +/-0.008, p < 0.01). Using previously published data, the mean bone density was inversely correlated with mean wall thickness in cancellous bone (p = 0.030) and with the percentage of active osteons (p = 0.023) in cortical bone. Although treatment had similar effects on the mean bone mineralization density of cortical and cancellous bone, there were different distributions of mineralization between the two sites, with cancellous bone having more skewed and kurtotic distributions both before and after estrogen withdrawal. This study indicates that a short-term estrogen suppression results in the accumulation of bone with a higher mineralization density. As bone with a high mineral content has a decreased impact resistance, this might increase fracture risk. Understanding the cellular and biochemical mechanisms responsible for the local distribution of bone mineral when estrogen is withdrawn may allow the development of new strategies for maintaining bone quality after menopause.

The thickness of human vertebral cortical bone and its changes in aging and osteoporosis: a histomorphometric analysis of the complete spinal column from thirty-seven autopsy specimens

The object of this study was to analyze the cortical thickness (Ct.Th) of the ventral and dorsal shell of the vertebral bodies throughout the human spine in aging and in osteoporosis. Therefore, the complete front column of the spine of 26 autopsy cases (aged 17-90, mean 42 years) without diseases affecting the skeleton and of 11 cases (aged 58-92, mean 77 years) with proven osteoporosis were removed. A sagittal segment prepared through the center of all vertebral bodies was undecalcified, embedded in plastic, ground to a 1 mm thick block, and stained using a modification of the von Kossa method. The analysis included the measurement of the mean cortical thickness of both the ventral and dorsal shell, respectively (from the third cervical to the fifth lumbar vertebral body). The qualitative investigation of the structure of the cortical ring completed the analysis. The presented data revealed a biphasic curve for both the ventral and dorsal shell, skeletally intact with high values of the cortical thickness in the cervical spine (285 microm), and a decrease in the thoracic (244 microm) and an increase in the lumbar spine (290 microm). The mean thickness of the ventral shell is in general greater than the thickness of the dorsal shell in both skeletally normal and osteoporotic cases. The cortical thickness of the spine showed no gender-specific differences (p = NS). There was a slight decrease of the cortical thickness with aging; however, this decrease and the correlation of cortical thickness to age was only significant below vertebral body T8 (r = 0.225-0.574; p(r) < 0.05-0.005). Most interestingly, however, osteoporosis presents itself with a highly significant loss of cortical thickness throughout the whole spine. This decrease of cortical thickness was more marked in the dorsal shell (p < 0.05) than in the ventral shell (ventral from C3 to T6 [p < 0.05] below T6 [p = NS]). We therefore conclude that in osteoporosis the loss of spinal bone mass is not only a loss of trabecular structure but also a loss of cortical thickness. Furthermore, these results may explain the development of regions of least resistance within the spine in aging and the clustering of osteoporotic fractures in the lower thoracic and lumbar spine.

Heterogeneity of the skeleton: comparison of the trabecular microarchitecture of the spine, the iliac crest, the femur, and the calcaneus

The objective of this study was to elucidate the structure of cancellous bone and its age-related changes at different skeletal sites. Therefore, the lumbar spine, iliac crest, femur, and calcaneus of 12 age- and sex-matched skeletal healthy autopsy cases (6 females, 6 males, aged 28-84 years, mean 54 years) were removed. The following analysis includes an evaluation of the trabecular bone volume (BV/TV, %) and the trabecular interconnection (TBPf, mm-1) as well as a qualitative investigation of the structure of trabecular bone. BV/TV shows the following mean values: lumbar spine, 8.3% (+/- 0.8%); iliac crest, 11.5% (+/- 1.6%); intertrochanteric, 10.2% (+/- 1.2%); femoral neck, 15.8% (+/- 1.6%); and calcaneus, 15.4% (+/- 2.0%). There are significant differences between the BV/TV of the femoral neck and that of the lumbar spine as well as between that of the calcaneus and the lumbar spine (p < 0.01). However, a positive correlation can be seen between the bone mass of the spine and that of all other investigated sites (r = 0.67 to r = 0.80; pr < 0.1). The trabecular interconnection of the lumbar spine (2.7 mm-1, SEM +/- 0.2 mm-1) and the femoral neck (0.3 mm-1, SEM +/- 0.3 mm-1) differs significantly. Only these two sites show a significant positive correlation of TBPf (r = 0.60; pr < 0.1). Age-dependent alteration of the spine and the femoral neck in bone mass and bone structure is nearly the same. The trabecular microarchitecture of the iliac crest varies systematically. A region 4-10 cm behind and 1-3 cm below the anterior superior iliac spine turns out to be the most suitable biopsy site because of its closest relation to the lumbar bone mass. However, drawing information about the trabecular interconnection within the lumbar spine by measurement of the iliac crest at any site seems to be impossible. The horizontal specimens reveal a vertical running tubular spongiosa pattern that is arranged in concentric rings starting from the dorsal shell like a honeycomb. The comparison of TBPf in horizontal and vertical planes and its age-related changes indicates the age-related bone loss to be predominantly a loss of horizontal trabeculae. Thus, the presented data provide further information about the skeletal distribution and heterogeneity of the trabecular microarchitecture.