Determinants of the mechanical behavior of human lumbar vertebrae after simulated mild fracture

Finite element models with automatic computed tomography bone segmentation for failure load computation

Bone segmentation is an important step to perform biomechanical failure load simulations on in-vivo CT data of patients with bone metastasis, as it is a mandatory operation to obtain meshes needed for numerical simulations. Segmentation can be a tedious and time consuming task when done manually, and expert segmentations are subject to intra- and inter-operator variability. Deep learning methods are increasingly employed to automatically carry out image segmentation tasks. These networks usually need to be trained on a large image dataset along with the manual segmentations to maximize generalization to new images, but it is not always possible to have access to a multitude of CT-scans with the associated ground truth. It then becomes necessary to use training techniques to make the best use of the limited available data. In this paper, we propose a dedicated pipeline of preprocessing, deep learning based segmentation method and post-processing for in-vivo human femurs and vertebrae segmentation from CT-scans volumes. We experimented with three U-Net architectures and showed that out-of-the-box models enable automatic and high-quality volume segmentation if carefully trained. We compared the failure load simulation results obtained on femurs and vertebrae using either automatic or manual segmentations and studied the sensitivity of the simulations on small variations of the automatic segmentation. The failure loads obtained using automatic segmentations were comparable to those obtained using manual expert segmentations for all the femurs and vertebrae tested, demonstrating the effectiveness of the automated segmentation approach for failure load simulations.

Computational comparison study of virtual compression and shear test for estimation of apparent elastic moduli under various boundary conditions

Virtual compression tests based on finite element analysis are representative noninvasive methods to evaluate bone strength. However, owing to the characteristic porous structure of bones, the material obtained from micro-computed tomography images in the finite-element model is not uniformly distributed. These characteristics cause differences in the apparent elastic moduli depending on the boundary conditions and affect the accuracy of bone-strength evaluation. Therefore, this study aimed to evaluate and compare the apparent elastic moduli under various, virtual-compression and shear-test boundary conditions. Four, nonuniform models were constructed with increasing model complexity. For representative boundary conditions, two, different, testing directions, and constrained surfaces were applied. As a result, the apparent elastic moduli of the nonuniform model varied up to 55.2% based on where the constrained surface was located in the single-end-cemented condition. Additionally, when connectivity in the test direction was lost, the accuracy of the apparent elastic moduli was low. A graphical comparison showed that the equivalent-stress distribution was more advantageous for analyzing load transferability and physical behavior than the strain-energy distribution. These results clearly show that the prediction accuracy of the apparent elastic moduli can be guaranteed if the boundary condition on the constraint and loading surfaces of the nonuniform model are applied symmetrically and the connectivity of the elements in the testing direction is well maintained. This study will aid in precision improvement of bone-strength-indicator determination for osteoporosis prevention.

Study on the scanning protocols for measuring bone mineral density by gemstone CT spectral imaging based on European spine phantom

BACKGROUND

Gemstone spectral computed tomography (GSCT) has been used to measure bone mineral density (BMD) in human vertebrae and animal models gradually.

PURPOSE

To investigate the effect of scanning protocols for BMD measurements by GSCT using the European spine phantom (ESP) and its accuracy and precision.

MATERIAL AND METHODS

The ESP number 145 containing three hydroxyapatite (HAP) inserts with densities of 50, 100, and 200 mg/cm³ were labeled as L₁, L₂, and L₃, respectively. Quantitative CT (QCT) protocol and 14 groups of scanning protocols configured by GSCT were used to repeatedly scan the ESP 10 times. Their measurements were compared with the true values of ESP and their relative standard deviation and relative error were calculated.

RESULTS

The measured values of the three inserts at different exposure levels were statistically significant (P < 0.05). The measured values in the 0.8 s/r 260 mA group, 0.5 s/r 630 mA group, and 0.6 s/r 640 mA group were not significantly different from the actual ESP values for L₁ and L₂. However, the measured values at all the parameters were significantly different from the actual values for the L₃.

CONCLUSION

CT gemstone spectral imaging can accurately and quantitatively measure the HAP value of ESP, but the results of BMD will be affected by the scanning protocols. The best scanning parameter of ESP measured by GSCT was 0.8 s/r 260 mA, taking dose into consideration, and the measurement accuracy of vertebrae with low BMD was higher than that of QCT under this parameter.

Dual-energy X-ray Absorptiometry Does Not Represent Bone Structure in Patients with Osteoporosis: A Comparison of Lumbar Dual-Energy X-Ray Absorptiometry with Vertebral Biopsies

STUDY DESIGN

Prospective cross-sectional exploratory study.

OBJECTIVE

To evaluate the correlation between in vivo lumbar dual-energy x-ray absorptiometry (DXA) and parameters of bone architecture in micro-computed tomography (micro-CT) in patients with osteoporosis.

SUMMARY OF BACKGROUND DATA

DXA is the current diagnostic standard for evaluating osteoporosis. However, there are various concerns regarding its validity, especially in the spine. No study has so far investigated whether in vivo DXA correlates with the actual lumbar bone architecture.

METHODS

Lumbar DXA scans were compared with micro-CT analysis of vertebral biopsies in patients with osteoporotic vertebral fractures (fracture group) and those without (control group). Preoperatively, all patients underwent a DXA scan (L1-L4). Intraoperative biopsies from nonfractured vertebrae (preferably L3) were analyzed by micro-CT regarding bone quantity and quality. The groups were compared regarding differences in DXA and micro-CT results. In each group, a correlation analysis was performed between DXA and micro-CT.

RESULTS

The study included 66 patients (33 per group). Preoperative DXA results were worse in the fracture group than the control group (areal bone mineral density [aBMD] 0.95 vs. 1.31, T-score -1.97 vs. 0.92, each P < 0.001). Micro-CT analysis confirmed differences regarding quantitative parameters (bone/total volume: 0.09 vs. 0.12, P < 0.001) and qualitative parameters (connectivity index: 15.73 vs. 26.67, P < 0.001; structure model index: 2.66 vs. 2.27, P < 0.001; trabecular number: 2.11 vs. 2.28, P = 0.014) of bone architecture between both groups. The DXA results did not correlate with micro-CT parameters in the fracture group. In the control group, correlations were found for some parameters (bone/total volume vs. aBMD: r = 0.51, P = 0.005; trabecular number vs. aBMD: r = 0.56, P = 0.001).

CONCLUSION

These data constitute the first comparison of DXA measurements with microstructural analysis of vertebral biopsies in patients with osteoporosis. Our results indicate that lumbar DXA neither qualitatively nor quantitatively represents microstructural bone architecture and is therefore not a reliable tool for the evaluation of bone quality in the spine.Level of Evidence: 3.

Relationship between diffuse idiopathic skeletal hyperostosis and fragility vertebral fracture: a prospective study in older men

OBJECTIVE

To analyse the risk of incident vertebral and non-vertebral fracture in men with DISH.

METHODS

In 782 men ages 50-85 years, DISH was diagnosed using Resnick's criteria. In men followed prospectively for 7.5 years, a radiographic incident vertebral fracture was defined by a decrease of ≥20% or ≥4mm in any vertebral height vs baseline. Self-reported incident non-vertebral fractures were confirmed by medical records.

RESULTS

Men with DISH had higher BMD at the lumbar spine (P < 0.05), but not at other skeletal sites. After adjustment for confounders including disc space narrowing (DSN) and endplate irregularity, the risk of vertebral fracture was higher in men with DISH vs men without DISH [10/164 (6.1%) vs 16/597 (2.7%), P < 0.05; odds ratio (OR) 2.89 (95% CI 1.15, 7.28), P < 0.05]. DISH and low spine BMD were each associated with a higher vertebral fracture risk. The vertebral fracture risk was higher in men who had both DISH and severe DSN. DISH and endplate irregularities (EIs) were each associated with higher vertebral fracture risk. DISH, DSN and EIs define the intervertebral space dysfunction, which was associated with higher vertebral fracture risk [OR 3.99 (95% CI 1.45, 10.98), P < 0.01]. Intervertebral space dysfunction improved the vertebral fracture prediction (ΔAUC = +0.111, P < 0.05), mainly in men with higher spine BMD (>0.9 g/cm2; ΔAUC = +0.189, P < 0.001). DISH was not associated with the risk of non-vertebral fracture.

CONCLUSION

DISH is associated with higher vertebral fracture risk, independently of other risk factors. Assessment of the intervertebral space dysfunction components may improve the vertebral fracture prediction in older men.

Local and global microarchitecture is associated with different features of bone biomechanics

Purpose

Beside areal bone mineral density (aBMD), evaluation of fragility fracture risk mostly relies on global microarchitecture. However, microarchitecture is not a uniform network. Therefore, this study aimed to compare local structural weakness to global microarchitecture on whole vertebral bodies and to evaluate how local and global microarchitecture was associated with bone biomechanics.

Methods

From 21 human L3 vertebrae, aBMD was measured using absorptiometry. Parameters of global microarchitecture were measured using HR-pQCT: trabecular bone volume fraction (Tb.BV/TV_global), trabecular number, structure model index and connectivity density (Conn.D). Local minimal values of aBMD and Tb.BV/TV were identified in the total (Tt) or trabecular (Tb) area of each vertebral body. “Two dimensional (2D) local structural weakness” was defined as Tt.BMD_min, Tt.BV/TV_min and Tb.BV/TV_min. Mechanical testing was performed in 3 phases: 1/ initial compression until mild vertebral fracture, 2/ unloaded relaxation, and 3/ second compression until failure.

Results

Initial and post-fracture mechanics were significantly correlated with bone mass, global and local microarchitecture. Tt.BMD_min, Tt.BV/TV_min, Tb.BV/TV_min, and initial and post-fracture mechanics remained significantly correlated after adjustment for aBMD or Tb.BV/TV_global (p < 0.001 to 0.038). The combination of the most relevant parameter of bone mass, global and local microarchitecture associated with failure load and stiffness demonstrated that global microarchitecture explained initial and post-fracture stiffness, while local structural weakness explained initial and post-fracture failure load (p < 0.001).

Conclusion

Local and global microarchitecture was associated with different features of vertebral bone biomechanics, with global microarchitecture controlling stiffness and 2D local structural weakness controlling strength. Therefore, determining both localized low density and impaired global microarchitecture could have major impact on vertebral fracture risk prediction.

•

Global and local microarchitecture were associated with different features of bone biomechanics.

•

Localized low density and/or impaired microarchitecture regions could have major impact on bone mechanical behavior.

•

Global microarchitecture determined initial and post-fracture vertebral stiffness.

•

Local microarchitecture determined initial and post-fracture vertebral failure load.

Bone structure determined by HR-MDCT does not correlate with micro-CT of lumbar vertebral biopsies: a prospective cross-sectional human in vivo study

Background

Osteoporosis is characterized by a deterioration of bone structure and quantity that leads to an increased risk of fractures. The primary diagnostic tool for the assessment of the bone quality is currently the dual-energy X-ray absorptiometry (DXA), which however only measures bone quantity. High-resolution multidetector computed tomography (HR-MDCT) offers an alternative approach to assess bone structure, but still lacks evidence for its validity in vivo. The objective of this study was to assess the validity of HR-MDCT for the evaluation of bone architecture in the lumbar spine.

Methods

We conducted a prospective cross-sectional study to compare the results of preoperative lumbar HR-MDCT scans with those from microcomputed tomography (μCT) analysis of transpedicular vertebral body biopsies. For this purpose, we included patients undergoing spinal surgery in our orthopedic department. Each patient underwent preoperative HR-MDCT scanning (L1-L4). Intraoperatively, transpedicular biopsies were obtained from intact vertebrae. Micro-CT analysis of these biopsies was used as a reference method to assess the actual bone architecture. HR-MDCT results were statistically analyzed regarding the correlation with results from μCT.

Results

Thirty-four patients with a mean age of 69.09 years (± 10.07) were included in the study. There was no significant correlation for any of the parameters (bone volume/total volume, trabecular separation, trabecular thickness) between μCT and HR-MDCT (bone volume/total volume: r = − 0.026 and p = 0.872; trabecular thickness: r = 0.074 and r = 6.42; and trabecular separation: r = − 0.18 and p = 0.254).

Conclusion

To our knowledge, this is the first study comparing in vivo HR-MDCT with μCT analysis of vertebral biopsies in human patients. Our findings suggest that lumbar HR-MDCT is not valid for the in vivo evaluation of bone architecture in the lumbar spine. New diagnostic tools for the evaluation of osteoporosis and preoperative orthopedic planning are urgently needed.

Standardizing compression testing for measuring the stiffness of human bone

Objectives

The ability to determine human bone stiffness is of clinical relevance in many fields, including bone quality assessment and orthopaedic prosthesis design. Stiffness can be measured using compression testing, an experimental technique commonly used to test bone specimens in vitro. This systematic review aims to determine how best to perform compression testing of human bone.

Methods

A keyword search of all English language articles up until December 2017 of compression testing of bone was undertaken in Medline, Embase, PubMed, and Scopus databases. Studies using bulk tissue, animal tissue, whole bone, or testing techniques other than compression testing were excluded.

Results

A total of 4712 abstracts were retrieved, with 177 papers included in the analysis; 20 studies directly analyzed the compression testing technique to improve the accuracy of testing. Several influencing factors should be considered when testing bone samples in compression. These include the method of data analysis, specimen storage, specimen preparation, testing configuration, and loading protocol.

Conclusion

Compression testing is a widely used technique for measuring the stiffness of bone but there is a great deal of inter-study variation in experimental techniques across the literature. Based on best evidence from the literature, suggestions for bone compression testing are made in this review, although further studies are needed to establish standardized bone testing techniques in order to increase the comparability and reliability of bone stiffness studies.

Cite this article: S. Zhao, M. Arnold, S. Ma, R. L. Abel, J. P. Cobb, U. Hansen, O. Boughton. Standardizing compression testing for measuring the stiffness of human bone. Bone Joint Res 2018;7:524–538. DOI: 10.1302/2046-3758.78.BJR-2018-0025.R1.

Low Dose of Bisphosphonate Enhances Sclerostin Antibody-Induced Trabecular Bone Mass Gains in Brtl/+ Osteogenesis Imperfecta Mouse Model

Osteogenesis imperfecta (OI) is a genetic disorder characterized by altered bone quality and imbalanced bone remodeling, leading to skeletal fractures that are most prominent during childhood. Treatments for OI have focused on restoring pediatric bone density and architecture to recover functional strength and consequently reduce fragility. Though antiresorptive agents like bisphosphonates (BPs) are currently the most common intervention for the treatment of OI, a number of studies have shown efficacy of sclerostin antibody (SclAb) in inducing gains in bone mass and reducing fragility in OI mouse models. In this study, the effects of the concurrent use of BP and SclAb were evaluated during bone growth in a mouse harboring an OI-causing Gly→Cys mutation on col1a1. A single dose of antiresorptive BP facilitated the anabolic action of SclAb by increasing availability of surfaces for new bone formation via retention of primary trabeculae that would otherwise be remodeled. Chronic effects of concurrent administration of BP and SclAb revealed that accumulating cycles conferred synergistic gains in trabecular mass and vertebral stiffness, suggesting a distinct advantage of both therapies combined. Cortical gains in mass and strength occurred through SclAb alone, independent of presence of BP. In conclusion, these preclinical results support the scientific hypothesis that minimal antiresorptive treatment can amplify the effects of SclAb during early stages of skeletal growth to further improve bone structure and rigidity, a beneficial outcome for children with OI. © 2018 American Society for Bone and Mineral Research.

Distal skeletal tibia assessed by HR-pQCT is highly correlated with femoral and lumbar vertebra failure loads

Dual energy X-ray absorptiometry (DXA) is the standard for assessing fragility fracture risk using areal bone mineral density (aBMD), but only explains 60-70% of the variation in bone strength. High-resolution peripheral quantitative computed tomography (HR-pQCT) provides 3D-measures of bone microarchitecture and volumetric bone mineral density (vBMD), but only at the wrist and ankle. Finite element (FE) models can estimate bone strength with 86-95% precision. The purpose of this study is to determine how well vBMD and FE bone strength at the wrist and ankle relate to fracture strength at the hip and spine, and to compare these relationships with DXA measured directly at those axial sites. Cadaveric samples (radius, tibia, femur and L4 vertebra) were compared within the same body. The radius and tibia specimens were assessed using HR-pQCT to determine vBMD and FE failure load. aBMD from DXA was measured at the femur and L4 vertebra. The femur and L4 vertebra specimens were biomechanically tested to determine failure load. aBMD measures of the axial skeletal sites strongly correlated with the biomechanical strength for the L4 vertebra (r=0.77) and proximal femur (r=0.89). The radius correlated significantly with biomechanical strength of the L4 vertebra for vBMD (r=0.85) and FE-derived strength (r=0.72), but not with femur strength. vBMD at the tibia correlated significantly with femoral biomechanical strength (r=0.74) and FE-estimated strength (r=0.83), and vertebral biomechanical strength for vBMD (r=0.97) and FE-estimated strength (r=0.91). The higher correlations at the tibia compared to radius are likely due to the tibia's weight-bearing function.

Vertebral body morphology is associated with incident lumbar vertebral fracture in postmenopausal women. The OFELY study

We investigate the predictive role of vertebral anterior cortical curvature and height heterogeneity in the occurrence of vertebral fractures in postmenopausal women. Women who will fracture had shorter vertebral height, greater heterogeneity of height than those who will not fracture, and their anterior vertebral body edge was less concave.

INTRODUCTION

Vertebral morphology has been demonstrated to be associated with further risk of fracture. The aim of this study was to analyze vertebral anterior cortical curvature (Ct.curv) and vertebral height heterogeneity in postmenopausal women before the occurrence of a vertebral fracture.

METHODS

This case-control study included 29 postmenopausal women who have underwent incident lumbar vertebral fractures (mean age 71 ± 9 years, mean time to fractures 9 ± 4 years), age-matched with 57 controls. From lateral X-rays of lumbar spine radiographs (T12 to L4), the following parameters were measured: (1) the posterior, middle, and anterior vertebral heights; (2) the heterogeneity of heights evaluated by the coefficient of variation of these three variables; (3) antero-posterior width, a 2D estimator of cross-sectional area; and (4) Ct.curv.

RESULTS

Mean vertebral heights were significantly lower among women who fractured than in controls (p < 0.05). The anterior and middle heights were significantly lower at L4 and L3 levels in fracture group (p = 0.02). The heterogeneity of vertebral height was significantly greater in the fracture group (p = 0.003). In addition, fractured patients had a significantly higher Ct.curv on L3 (p = 0.04). After adjustment for bone mineral density (BMD), only the heterogeneity of vertebral height remained significant (p = 0.005).

CONCLUSION

The current case-control study confirmed the association between vertebral height and occurrence of future vertebral fracture in postmenopausal women. The vertebrae with the smallest Ct.curv tended to fracture less often, and the heterogeneity of vertebral heights was associated with future fracture independently of BMD. An additional validation in a prospective study would be needed to confirm these initial results.

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

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

Optimal sample volumes of human trabecular bone in μCT analysis within vertebral body and femoral head

Trabecular bones of different skeletal sites have different bone morphologies. How to select an appropriate volume of region of interest (ROI) to reflect the microarchitecture of trabecular bone in different skeletal sites was an interesting problem. Therefore, in this study, the optimal volumes of ROI within vertebral body and femoral head, and if the relationships between volumes of ROI and microarchitectural parameters were affected by trabecular bone morphology were studied. Within vertebral body and femoral head, different cubic volumes of ROI (from (1 mm)(3) to (20 mm)(3)) were set to compare with control groups(whole volume of trabecular bone). Five microarchitectural parameters (BV/TV, Tb.N, Tb.Th, Tb.Sp, and BS/BV) were obtained. Nonlinear curve fitting functions were used to explore the relationships between the microarchitectural parameters and the volumes of ROI. The volumes of ROI could affect the microarchitectural parameters when the volume was smaller than (8 mm)(3) within the vertebral body and smaller than (13 mm)(3) within the femoral head. As the volume increased, the variable tendencies of BV/TV, Tb.N, and Tb.Sp were different between these two skeletal sites. The curve fitting functions between these two sites were also different. The relationships between volumes of ROI and microarchitectural parameters were affected by the different trabecular bone morphologies within lumbar vertebral body and femoral head. When depicting the microarchitecture of human trabecular bone within lumbar vertebral body and femoral head, the volume of ROI would be larger than (8 mm)(3) and (13 mm)(3).

Specimen-specific vertebral fracture modeling: a feasibility study using the extended finite element method

Osteoporotic vertebral body fractures are an increasing clinical problem among the aging population. Specimen-specific finite element models, derived from quantitative computed tomography (QCT), have the potential to more accurately predict failure loads in the vertebra. Additionally, the use of extended finite element modeling (X-FEM) allows for a detailed analysis of crack initiation and propagation in various materials. Our aim was to study the feasibility of QCT/X-FEM analysis to predict fracture properties of vertebral bodies. Three cadaveric specimens were obtained, and the L3 vertebrae were excised. The vertebrae were CT scanned to develop computational models and mechanically tested in compression to measure failure load, stiffness and to observe crack location. One vertebra was used for calibration of the material properties from experimental results and CT gray-scale values. The two additional specimens were used to assess the model prediction. The resulting QCT/X-FEM model of the specimen used for calibration had 2 and 4% errors in stiffness and failure load, respectively, compared with the experiment. The predicted failure loads of the additional two vertebrae were larger by about 41-44% when compared to the measured values, while the stiffness differed by 129 and 40%. The predicted fracture patterns matched fairly well with the visually observed experimental cracks. Our feasibility study indicated that the QCT/X-FEM method used to predict vertebral compression fractures is a promising tool to consider in future applications for improving vertebral fracture risk prediction in the elderly.

The role of patient-mode high-resolution peripheral quantitative computed tomography indices in the prediction of failure strength of the elderly women's thoracic vertebral body

The correlations between the failure load of 20 T12 vertebral bodies, their patient-mode high-resolution peripheral quantitative computed tomography (HR-pQCT) indices, and the L1 areal bone mineral density (aBMD) were investigated. For the prediction of the T12 vertebral failure load, the T12 HR-pQCT microarchitectural parameters added significant information to that of L1 aBMD and to that of cortical BMD, but not to that of T12 vertebral BMD and not to that of T12 trabecular BMD.

INTRODUCTION

HR-pQCT is a new in vivo imaging technique for assessing the three-dimensional microarchitecture of cortical and trabecular bone at the distal radius and tibia. But little is known about this technique in the direct measurement of vertebral body.

METHODS

Twenty female donors with the mean age of 80.1 (7.6) years were included in the study. Dual X-ray absorptiometry of the lumbar spine and femur was performed. The spinal specimens (T11/T12/L1) were dissected, scanned using HR-pQCT scanner, and mechanically tested under 4° wedge compression. The L1 aBMD, T12 patient-mode HR-pQCT indices, and T12 vertebral failure loads were analyzed.

RESULTS

For the prediction of vertebral failure load, the inclusion of BV/TV into L1 aBMD was the best model (R (2) = 0.52), Tb.N and Tb.Sp added significant information to the L1 aBMD and to the cortical BMD, but none of the vertebral microarchitectural parameters yielded additional significant information to the trabecular BMD (or BV/TV) and to the vertebral BMD.

CONCLUSION

Vertebral microarchitectural parameters obtained from the patient-mode HR-pQCT analysis provide significant information on bone strength complementary to that of aBMD and to that of cortical BMD, but not to that of vertebral BMD and not to that of trabecular BMD.

The role of bone intrinsic properties measured by infrared spectroscopy in whole lumbar vertebra mechanics: organic rather than inorganic bone matrix?

Whole bone strength is determined by bone mass, microarchitecture and intrinsic properties of the bone matrix. However, few studies have directly investigated the contribution of bone tissue material properties to whole bone strength in humans. This study assessed the role of bone matrix composition on whole lumbar vertebra mechanics. We obtained 17 fresh frozen human lumbar spines (8 W, 9 M, aged 76±11years). L3 bone mass was measured by DXA and microarchitecture by μ-CT with a 35 μm-isotropic resolution. Microarchitectural parameters were directly measured: Tb.BV/TV, SMI, Tb.Th, DA, Ct.Th, Ct.Po and radius of anterior cortical curvature. Failure load (N), stiffness (N/mm) and work to failure (N.mm) were extracted from quasi-static uniaxial compressive testing performed on L3 vertebral bodies. FTIRM analysis was performed on 2 μm-thick sections from L2 trabecular cores, with a Perkin-Elmer GXII Auto-image Microscope equipped with a wide band detector. Twenty measurements per sample were performed at 30∗100 μm of spatial resolution. Each spectrum was collected at 4 cm(-1) resolution and 50 scans in transmission mode. Mineral and collagen maturity, and mineralization and crystallinity index were measured. There was no association between the bone matrix characteristics and bone mass or microarchitecture. Mineral maturity, mineralization and crystallinity index were not related to whole vertebra mechanics. However, collagen maturity was positively correlated with whole vertebra failure load and stiffness (r=0.64, p=0.005 and r=0.54, p=0.025, respectively). The collagen maturity (3rd step) in combination with bone mass (i.e., BMC, 1st step) and microarchitecture (i.e., Tb.Th, 2nd step) improved the prediction of whole vertebra mechanical properties in forward stepwise multiple regression models, together explaining 71% of the variability in whole vertebra stiffness (p=0.001). In conclusion, we demonstrated a substantial contribution of collagen maturity, but not mineralization parameters, to whole bone strength of human lumbar vertebrae that was independent of bone mass and microarchitecture.

Vertebral fragility and structural redundancy

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

Trabecular architecture and vertebral fragility in osteoporosis

Osteoporosis heightens vertebral fragility owing to the biomechanical effects of diminished bone structure and composition. These biomechanical effects are only partially explained by loss in bone mass, so additional factors that are independent of bone mass are also thought to play an important role in vertebral fragility. Recent advances in imaging equipment, imaging-processing methods, and computational capacity allow researchers to quantify trabecular architecture in the vertebra at the level of the individual trabecular elements and to derive biomechanics-based measures of architecture that are independent of bone mass and density. These advances have shed light on the role of architecture in vertebral fragility. In addition to the adverse biomechanical consequences associated with trabecular thinning and loss of connectivity, a reduction in the number of vertical trabecular plates appears to be particularly harmful to vertebral strength. In the clinic, detailed architecture analysis is primarily applied to peripheral sites such as the distal radius and tibia. Analysis of trabecular architecture at these peripheral sites has shown mixed results for discriminating between patients with and without a vertebral fracture independent of bone mass, but has the potential to provide unique insight into the effects of therapeutic treatments. Overall, it does appear that trabecular architecture has an independent role on vertebral strength. Additional research is required to determine how and where architecture should be measured in vivo and whether assessment of trabecular architecture in a clinical setting improves prospective fracture risk assessment for the vertebra.

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

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

Lathyrism-induced alterations in collagen cross-links influence the mechanical properties of bone material without affecting the mineral

In the present study a rat animal model of lathyrism was employed to decipher whether anatomically confined alterations in collagen cross-links are sufficient to influence the mechanical properties of whole bone. Animal experiments were performed under an ethics committee approved protocol. Sixty-four female (47 day old) rats of equivalent weights were divided into four groups (16 per group): Controls were fed a semi-synthetic diet containing 0.6% calcium and 0.6% phosphorus for 2 or 4 weeks and β-APN treated animals were fed additionally with β-aminopropionitrile (0.1% dry weight). At the end of this period the rats in the four groups were sacrificed, and L2-L6 vertebra were collected. Collagen cross-links were determined by both biochemical and spectroscopic (Fourier transform infrared imaging (FTIRI)) analyses. Mineral content and distribution (BMDD) were determined by quantitative backscattered electron imaging (qBEI), and mineral maturity/crystallinity by FTIRI techniques. Micro-CT was used to describe the architectural properties. Mechanical performance of whole bone as well as of bone matrix material was tested by vertebral compression tests and by nano-indentation, respectively. The data of the present study indicate that β-APN treatment changed whole vertebra properties compared to non-treated rats, including collagen cross-links pattern, trabecular bone volume to tissue ratio and trabecular thickness, which were all decreased (p<0.05). Further, compression tests revealed a significant negative impact of β-APN treatment on maximal force to failure and energy to failure, while stiffness was not influenced. Bone mineral density distribution (BMDD) was not altered either. At the material level, β-APN treated rats exhibited increased Pyd/Divalent cross-link ratios in areas confined to a newly formed bone. Moreover, nano-indentation experiments showed that the E-modulus and hardness were reduced only in newly formed bone areas under the influence of β-APN, despite a similar mineral content. In conclusion the results emphasize the pivotal role of collagen cross-links in the determination of bone quality and mechanical integrity. However, in this rat animal model of lathyrism, the coupled alterations of tissue structural properties make it difficult to weigh the contribution of the anatomically confined material changes to the overall mechanical performance of whole bone. Interestingly, the collagen cross-link ratio in bone forming areas had the same profile as seen in actively bone forming trabecular surfaces in human iliac crest biopsies of osteoporotic patients.