Vertebral Failure Research Articles

CT-based finite element simulating spatial bone damage accumulation predicts metastatic human vertebrae strength and stiffness.

Introduction: Pathologic vertebral fractures are devastating for patients with spinal metastases. However, the mechanical process underlying these fractures is poorly understood, limiting physician's ability to predict which vertebral bodies will fail. Method: Here, we show the development of a damage-based finite element framework producing highly reliable pathologic vertebral strength and stiffness predictions from X-Ray computed tomography (CT) data. We evaluated the performance of specimen-specific material calibration vs. global material calibration across osteosclerotic, osteolytic, and mixed lesion vertebrae that we derived using a machine learning approach. Results: The FE framework using global calibration strongly predicted the pathologic vertebrae stiffness (R 2 = 0.90, p < 0.0001) and strength (R 2 = 0.83, p = 0.0002) despite the remarkable variance in the pathologic bone structure and density. Specimen-specific calibration produced a near-perfect prediction of both stiffness and strength (R 2 = 0.99, p < 0.0001, for both), validating the FE approach. The FE damage-based simulations highlighted the differences in the pattern of spatial damage evolution between osteosclerotic and osteolytic vertebral bodies. Discussion: With failure, the FE simulation suggested a common damage evolution pathway progressing largely localized to the low bone modulus regions within the vertebral volume. Applying this FE approach may allow us to predict the onset and anatomical location of vertebral failure, which is critical for developing image-based diagnostics of impending pathologic vertebral fractures.

Modified post-installed adhesive anchorage test: Investigation of in-place compressive strength of concrete and influential factors

This paper presents the Modified Post-installed Adhesive Anchorage test (MPAA) as a method for evaluating the in-place compressive strength of concrete. It considers various influential factors, including test surface characteristics, concrete strength, anchorage depth, bolt diameter, and rebar location. Additionally, a formula for calculating concrete compressive strength, applicable to various influencing factors, was suggested. The findings demonstrate the effectiveness of MPAA in assessing in-place compressive strength of concrete with a high degree of precision. Following pull-out test, the prevalent failure modes observed were vertebral failure and vertebral-bond composite failure. Moreover, diverse concrete strengths, anchorage depths, and rebar locations demonstrated differing impacts on failure modes, while test surface characteristics and bolt diameter exerted less noticeable effects. The pull-out force-displacement curve during pull-out displayed three discernible phases: the elastic phase, crack propagation phase, and descent phase. When conducting on-site assessments, it is advisable to situate test points on the side surfaces of the members. Furthermore, factors like cover thickness and rebar location should be thoroughly evaluated when positioning test points. This encompasses the prudent choice of anchorage depth and bolt diameter to improve test precision. The suggested formula for calculating compressive strength exhibited minimal error and is applicable in determining concrete compressive strength for a wide range of practical engineering contexts. This study not only establishes a theoretical underpinning and offers support for standard revisions but also furnishes valuable guidance for conducting on-site testing procedures.

Vertebral Fracture Risk Thresholds from Phantom-Less Quantitative Computed Tomography-Based Finite Element Modeling Correlate to Phantom-Based Outcomes

Introduction: Osteoporosis indicates weakened bones and heightened fracture susceptibility due to diminished bone quality. Dual-energy x-ray absorptiometry is unable to assess bone strength. Volumetric bone mineral density (vBMD) from quantitative computed tomography (QCT) has been used to establish guidelines as equivalent measurements for osteoporosis. QCT-based finite element analysis (FEA) has been implemented using calibration phantoms to establish bone strength thresholds based on the established vBMD. The primary aim was to validate vertebral failure load thresholds using a phantom-less approach with previously established thresholds, advancing a phantom-free approach for fracture risk prediction.Methodology: A controlled cohort of 108 subjects (68 females) was used to validate sex-specific vertebral fracture load thresholds for normal, osteopenic, and osteoporotic subjects, obtained using a QCT/FEA-based phantom-less calibration approach and two material equations.Results: There were strong prediction correlations between the phantom-less and phantom-based methods (R2: 0.95 and 0.97 for males, and R2: 0.96 and 0.98 for females) based on the two equations. Bland Altman plots and paired t-tests showed no significant differences between methods. Predictions for bone strengths and thresholds using the phantom-less method matched those obtained using the phantom calibration and those previously established, with ≤4500 N (fragile) and ≥6000 N (normal) bone strength in females, and ≤6500 N (fragile) and ≥8500 N (normal) bone strength in males.Conclusion: Phantom-less QCT-based FEA can allow for prospective and retrospective studies evaluating incidental vertebral fracture risk along the spine and their association with spine curvature and/or fracture etiology. The findings of this study further supported the application of phantom-less QCT-based FEA modeling to predict vertebral strength, aiding in identifying individuals prone to fractures. This reinforces the rationale for adopting this method as a comprehensive approach in predicting and managing fracture risk.

A Comprehensive Flexural Analysis for Sustainable Concrete Structure Reinforced by Embedded Parts

In this research, the plate embedded parts and grooved embedded parts reinforced concrete structures were investigated. Two types of plate embedded parts and three types of grooved embedded parts experienced coating treatment to enable sustainable function. Later, the ultimate failure capacity by bending experiments was conducted and compared with the theoretically calculated results. Moreover, three grooved embedded parts were simulated by ABAQUS to compare the results with the experimental exploration results, which was in close agreement with the theoretically calculated results and finite element analysis results. The result indicated that the failure modes of the embedded specimens under the five working conditions are all concrete vertebral failure. The plate-type embedded components were proved to exhibit higher ultimate bearing capacity than the grooved embedded parts. Moreover, the flexural and shear capacity of these five types of embedded parts has not been fully developed. The ultimate flexural and shear capacity of these five types of embedded parts could be further explored by adjusting the higher concrete grade.

Multi-detector computed tomography (MDCT) imaging: association of bone texture parameters with finite element analysis (FEA)-based failure load of single vertebrae and functional spinal units.

Osteoporosis is a systemic skeletal disease that is characterized by low bone mass and microarchitectural deterioration, predisposing affected individuals to fragility fractures. Yet, standard measurement of areal bone mineral density (BMD) in dual-energy X-ray absorptiometry (DXA) as the current reference standard has limitations for correctly detecting osteoporosis and fracture risk, with opportunistic osteoporosis screening using computed tomography (CT) showing increasing importance. This study's objective is to compare finite element analysis (FEA)-based vertebral failure load with parameters of texture analysis (TA) derived from multi-detector CT (MDCT). MDCT data of seven subjects (mean age: 71.9±7.4 years) were included for FEA and TA. Manual segmentation was performed for the vertebral bodies T11, T12, L1, and L2 and the intervertebral discs (IVDs) T11/12, T12/L1, L1/2, and L2/3. Correlation analyses between FEA-derived failure loads and parameters of TA for the single vertebrae and two functional spinal units (FSUs) were calculated, defining FSU-1 as T11-IVD-T12-IVD-L1 and FSU-2 as T12-IVD-L1-IVD-L2. Furthermore, multivariate regressions were performed to identify the texture parameters that predicted the failure load best. For single vertebrae, the strongest correlations were observed for skewnessglobal, kurtosisglobal, and gray level variance (rho =-0.7668 to -0.7362; P<0.001), while for FSUs, SumAverage, long-run emphasis, long-run low gray-level emphasis, homogeneity, and energy showed the strongest correlations (rho =-0.8187 to 0.8407; P<0.05) to failure loads. SumAverage best predicted the failure load for single vertebrae (R2 adj =0.523, P<0.001). For the two FSUs, kurtosisglobal (FSU-1: R2 adj =0.611, P=0.001) and skewnessglobal (FSU-2: R2 adj =0.579, P=0.002) were the best predictors. TA using MDCT data of the spine was significantly associated with FEA-derived failure loads of both, single vertebrae and FSUs. Texture parameters predicted failure loads of FSUs as a more realistic in-vivo scenario equally well as compared to single vertebrae analyses. TA may reflect a less complex and time-consuming approach to accurately and non-invasively evaluate vertebral bone strength.

Bisphosphonate treatment changes regional distribution of trabecular microstructure in human lumbar vertebrae.

In osteoporosis patients, antiresorptive treatments such as alendronate reduce the resorption of trabecular bone and thus minimize vertebral fracture risk. However, fracture risk reduction efficacy of antiresorptive drugs varies between skeletal sites and is highest for vertebral bone. In human vertebrae, cancellous bone is distributed heterogeneously between regions. This microstructural heterogeneity is changing with patient age and is likely to play a major role in vertebral failure mechanisms and fracture susceptibility. Whether antiresorptive treatment affects the heterogeneity of vertebral microstructure in osteoporosis has not been unraveled. Our aim was to assess whether antiresorptive treatment would have a region-dependent influence on vertebral trabecular bone. Therefore, we used high-resolution peripheral quantitative computed tomography (HR-pQCT), microcomputed tomography (microCT) and uniaxial compression testing to determine the structure and mechanical properties of trabecular bone cores from anterior and posterior regions of 22 lumbar vertebrae from elderly osteoporotic women. We analyzed age-matched ex vivo bone samples from bisphosphonate-treated female osteoporosis patients (age: 82 ± 7y, bisphosphonate treatment period: 4 ± 2 years) along treatment-naïve female controls (82 ± 7y). MicroCT analysis showed a significantly lower bone volume fraction (p = 0.006) and lower trabecular number (p = 0.003) for the anterior bone cores compared to posterior bone cores in the treatment-naïve group. The bisphosphonate-treated group had a more homogeneous bone volume distribution and did not show significant regional differences in bone volume, it however also displayed significantly different trabecular numbers (p = 0.016). In bone cores of the bisphosphonate-treated group, trabeculae were thicker in comparison to treatment-naïve controls (p = 0.011). Differences in bone volume further resulted in different maximum forces during compression testing between the samples. In addition, the percental difference between BV/TVμCT in anterior and posterior bone cores was lower in bisphosphonate-treated vertebrae when vertebrae with directly adjacent fractures (n = 3) were excluded. In conclusion, regional trabecular bone microstructure in lumbar vertebrae of bisphosphonate-treated women was more homogeneous compared to treatment-naïve controls. Bisphosphonate treatment, which specifically targets resorption surfaces common in anterior vertebral bone, might have resulted in a region-specific preservation of vertebral microstructure and loading capacity. This could have positive implications for the reduction of wedge fracture risk and add to the explanation of the higher efficacy of fracture risk reduction in vertebrae in comparison to other fracture regions.

Evaluation of the capability of the simulated dual energy X-ray absorptiometry-based two-dimensional finite element models for predicting vertebral failure loads

Prediction of the vertebral failure load is of great importance for the prevention and early treatment of bone fracture. However, an efficient and effective method for accurately predicting the failure load of vertebral bones is still lacking. The aim of the present study was to evaluate the capability of the simulated dual energy X-ray absorptiometry (DXA)-based finite element (FE) model for predicting vertebral failure loads.Thirteen dissected spinal segments (T11/T12/L1) were scanned using a HR-pQCT scanner and then were mechanically tested until failure. The subject-specific three-dimensional (3D) and two-dimensional (2D) FE models of T12 were generated from the HR-pQCT scanner and the simulated DXA images, respectively. Additionally, the areal bone mineral density (aBMD) and areal bone mineral content (aBMC) of T12 were calculated. The failure loads predicted by the simulated DXA-based 2D FE models were more moderately correlated with the experimental failure loads (R2 = 0.66) than the aBMC (R2 = 0.61) and aBMD (R2 = 0.56). The 2D FE models were slightly outperformed by the HR-pQCT-based 3D FE models (R2 = 0.71). The present study demonstrated that the simulated DXA-based 2D FE model has better capability for predicting the vertebral failure loads than the densitometric measurements but is outperformed by the 3D FE model. The 2D FE model is more suitable for clinical use due to the low radiation dose and low cost, but it remains to be validated by further in vitro and in vivo studies.

A novel technique with reduced computed tomography exposure to predict vertebral compression fracture: a finite element study based on rat vertebrae.

Vertebral compression fractures are a significant clinical issue with an annual incidence of approximately 750,000 cases in the USA alone. Mechanical properties of vertebrae are successfully evaluated through finite element (FE) models based on vertebrae CT. However, clinical drawbacks associated to radiation transmission encouraged to explore the possibility to use selected or reduced portions of the vertebra. The objective of our study was to develop a new procedure to predict vertebral compression fracture from sub-volumes. We reconstructed rat vertebras from micro-CT of thoracic and lumbar groups. Each vertebra was partitioned into three sub-volumes of different axial thickness. FE simulating compression tests were performed on each model to evaluate their failure load and stiffness. Using a power function, a high correlation was found for stiffness and strength. The sub-volume with three fifths thickness had a failure load of 180.7 ± 19.2N for thoracic and of 209.5 ± 27.4N for the lumbar vertebra. These values were not significantly different from the values found for the entire vertebra (p > 0.05). Based on our findings, failure loads and stiffnesses obtained with reduced CT scans can be successfully used to predict full vertebral failure. This sub-region analysis and power relationship suggests that one can limit radiation exposure to patients when bone characterization is needed. Graphical abstract Estimated mechanical properties in relation to the extent of the computed tomography reconstruction.

MDCT-based Finite Element Analysis of Vertebral Fracture Risk: What Dose is Needed?

The aim of this study was to compare vertebral failure loads, predicted from finite element (FE) analysis of patients with and without osteoporotic vertebral fractures (OVF) at virtually reduced dose levels, compared to standard-dose exposure from multidetector computed tomography (MDCT) imaging and evaluate whether ultra-low dose derived FE analysis can still differentiate patient groups. An institutional review board (IRB) approval was obtained for this retrospective study. A total of 16patients were evaluated at standard-dose MDCT; eight with and eight without OVF. Images were reconstructed at virtually reduced dose levels (i. e. half, quarter and tenth of the standard dose). Failure load was determined at L1-3 from FE analysis and compared between standard, half, quarter, and tenth doses and used to differentiate between fracture and control groups. Failure load derived at standard dose (3254 ± 909 N and 3794 ± 984 N) did not significantly differ from half (3390 ± 890 N and 3860 ± 1063 N) and quarter dose (3375 ± 915 N and 3925 ± 990 N) but was significantly higher for one tenth dose (4513 ± 1762 N and 4766 ± 1628 N) for fracture and control groups, respectively. Failure load differed significantly between the two groups at standard, half and quarter doses, but not at tenth dose. Receiver operating characteristic (ROC) curve analysis also demonstrated that standard, half, and quarter doses can significantly differentiate the fracture from the control group. The use of MDCT enables adose reduction of at least 75% compared to standard-dose for an adequate prediction of vertebral failure load based on non-invasive FE analysis.

The trabecular effect: A population-based longitudinal study on age and sex differences in bone mineral density and vertebral load bearing capacity

BackgroundApproximately 16–24% of postmenopausal women are affected by vertebral fractures, negatively affecting their quality of life. Trabecular and cortical bones in vertebrae decline differently with age, thus having a distinct impact on vertebral failure loads. The purpose of this study was to investigate the effect of trabecular and cortical volumetric bone mineral density loss over time on estimated failure loads; and to evaluate the effect of sex and age. MethodFracture properties from a cohort of 82 patients were evaluated for L1–L3 vertebrae at baseline and 6th year using an image-based method that implements axial rigidity analysis. Cortical and trabecular volumetric bone mineral density were obtained, as well as their individual contribution to total failure load. Regression analyses were performed to determine the effect of age and sex on volumetric bone mineral density and failure loads. FindingsDecline in trabecular and cortical volumetric bone mineral density, and failure load was sex-dependent (p ≤ 0.0095). Cortical and trabecular volumetric bone mineral density reduced 2.08 (g/cm3)/year and 2.02 (g/cm3)/year, respectively. A 1012 N difference in failure load, ~70% attributed to trabecular bone, was found between men and women of similar age. Over 6 years, this difference increased by 287 N. Areal bone mineral density measured by dual X-ray absorptiometry explained ~60% of the vertebral failure load. InterpretationTrabecular bone has a significantly greater effect than cortical bone on the structural integrity and load bearing capacity of vertebrae. This might lead to a higher incidence of fragility fractures in osteoporotic women. Our non-invasive, quantitative computed tomography image-based approach may improve prevention, monitoring, and management of fractures.

Stress and strain analyses of single and segmental lumbar spines based on an accurate finite element model for vertebrae

Two types of accurate 3D finite element models were developed for lumbar spine L5 and L3-L5, and a series of computational analyses were carried out to identify the mechanical behavior of lumbar spine under compressive loads. In order to establish a precise FE model for single and segmental lumbar vertebrae, a group of computed tomography images was converted into a 3D geometric surface, and a tetrahedral mesh using the image processing software. For the FE analysis, two types of bones, namely, a fully filled type with cortical bone, and a cortical bone with a cancellous bone core, was considered for the investigation of the effects of the existence of central core of bone. The computational results for L5 indicate that a larger von Mises stress distribution was found on the pedicles and on the vertebral body with respect to the remaining vertebral parts. In particular, the highest equivalent stress was occurred on the surface of cortical bone in the case of a mixture of cortical and cancellous bone. On the other hand, the numerical results of L3-L5 demonstrate that the highest tensile and compressive principal strains were generated on the posterior vertebral rims, and on the pedicles and pars interarticularis. Through the comparative study, it was confirmed that the calculation results related to L5 and L3-L5 were in close agreement with the experiments. Based on the introduced method, the vertebral failure as well as the stress-strain behavior can be specifically identified.

What is the future of patient-specific vertebral fracture prediction?

This paper aims to introduce a few alternative methodologies for prediction of vertebral fractures, the most common type being fragility fracture in the elderly. Current methods, such as DXA, for diagnosing osteoporosis and predicting the risk of vertebral failure, are often not accurate thereby preventing those patients at risk from receiving adequate treatment. Robust fracture prediction models for vertebral fracture risk should not only include BMD, as measured by DXA, but should incorporate a wide range of factors including bone geometry, bone mineral distribution within the vertebral body, daily living activities, and spine musculature. One promising technique is finite element modeling, which has been developed over the past several decades and implements clinical imaging, such as quantitative computed tomography (QCT), and engineering fundamentals to more accurately predict the risk of fracture. Other imaging tools that assess bone mineral distribution and structure at the microscopic level include micro-CT or high-resolution peripheral QCT (HR-pQCT). These newer techniques hold the promise of more accurate diagnosis of osteoporosis and those at risk for vertebral insufficiency fractures before they occur.

Obesity and Obesity Shape Markedly Influence Spine Biomechanics: A Subject-Specific Risk Assessment Model.

Underlying mechanisms of obesity-related back pain remain unexplored. Thus, we aim to determine the effect of obesity and its shapes on the spinal loads and the associated risks of injury. Obesity shapes were initially constructed by principal component analysis based on datasets on 5852 obese individuals. Spinal loads, cycles to vertebral failure and trunk stability margin were estimated in a subject-specific trunk model taking account of personalized musculature, passive ligamentous spine, obesity shapes, segmental weights, spine kinematics and bone mineral density. Three obesity shapes (mean and extreme abdominal circumferences) at three body weights (BWs) of 86, 98 and 109kg were analyzed. Additional BW (12kg) increased spinal loads by ~11.8%. Higher waist circumferences at identical BW increased spinal forces to the tune of ~20kg additional BW and the risk of vertebral fatigue compression fracture by 3-7 times when compared with smaller waist circumferences. Forward flexion, greater BW and load in hands increased the trunk stability margin. Spinal loads markedly increased with BW, especially at greater waist circumferences. The risk of vertebral fatigue fracture also substantially increased at greater waist circumferences though not at smaller ones. Obesity and its shape should be considered in spine biomechanics.

Mechanical assessment of the effects of metastatic lytic defect on the structural response of human thoracolumbar spine.

To investigate the effects of a clinical lytic defect on the structural response of human thoracolumbar functional spinal unit. A novel CT-compatible mechanical test system was used to image the deformation of a T12-L1 motion segment and measure the change in strain response under compressive loads ranging from 50 to 750 N. A lytic lesion (LM) with cortex involvement (33% by volume) was introduced to the upper vertebral body and the CT experiments were repeated. Finite element models, established from the CT volumes, were used to investigate the defect's effects on the structural response and the state of principal and shear stresses within the affected and adjacent vertebrae. The lytic lesion resulted in severe loss of the vertebral structural competence, resulting in significant, non-linear, and asymmetric increase in the experimentally measured strains and computed stresses within both vertebrae (p < 0.01). At the cortex, the tensile strains were significantly increased, while compressive strains significantly decreased, (p < 0.05). Both the vertebral bone and cortex regions adjacent to the defect showed significant increase in computed compressive, tensile, and shear stresses (p < 0.01). Changes in stress and strain distribution within the affected and adjacent vertebral bone and the experimentally observed bulging and buckling of the vertebral cortices suggested that initiation of catastrophic vertebral failure may occur under load magnitudes encountered in daily living. Although the effect of LM on the global deformation of the spine was well-predicted, our results show that FE predictions of local strain changes must be carefully assessed for clinical relevance. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1808-1819, 2016.

Accuracy of finite element analyses of CT scans in predictions of vertebral failure patterns under axial compression and anterior flexion

Finite element (FE) models built from quantitative computed tomography (QCT) scans can provide patient-specific estimates of bone strength and fracture risk in the spine. While prior studies demonstrate accurate QCT-based FE predictions of vertebral stiffness and strength, the accuracy of the predicted failure patterns, i.e., the locations where failure occurs within the vertebra and the way in which the vertebra deforms as failure progresses, is less clear. This study used digital volume correlation (DVC) analyses of time-lapse micro-computed tomography (μCT) images acquired during mechanical testing (compression and anterior flexion) of thoracic spine segments (T7–T9, n=28) to measure displacements occurring throughout the T8 vertebral body at the ultimate point. These displacements were compared to those simulated by QCT-based FE analyses of T8. We hypothesized that the FE predictions would be more accurate when the boundary conditions are based on measurements of pressure distributions within intervertebral discs of similar level of disc degeneration vs. boundary conditions representing rigid platens. The FE simulations captured some of the general, qualitative features of the failure patterns; however, displacement errors ranged 12–279%. Contrary to our hypothesis, no differences in displacement errors were found when using boundary conditions representing measurements of disc pressure vs. rigid platens. The smallest displacement errors were obtained using boundary conditions that were measured directly by DVC at the T8 endplates. These findings indicate that further work is needed to develop methods of identifying physiological loading conditions for the vertebral body, for the purpose of achieving robust, patient-specific FE analyses of failure mechanisms.

A BIOMECHANICAL STUDY OF SHEAR LOAD ON BONE–SCREW INTERFACE OF THORACOLUMBAR VERTEBRAE

Vertebral screw failure by cutting through or pulling out of the vertebral body is common in intraoperative and postoperative stages of anterior thoracolumbar surgery, especially in osteoporotic patients. This biomechanical study was conducted to investigate the maximum shear force thoracolumbar vertebrae can withstand and the corresponding displacement of screws and analyze their correlation with vertebral bone mineral density. Forty individual vertebra specimens (T11-L3) were obtained from eight fresh adult cadaveric thoracolumbar spine specimens and randomly divided into an experimental group and a control group. Screws were placed in the center point of vertebrae and penetrated the contralateral cortex. Shear loading and axial pullout tests were conducted on the experimental group while only axial pullout test was conducted on the control group. The maximum shear forces and maximum axial pullout forces were recorded. The conditions of vertebral body destruction and screw channel were observed and the maximum axial pullout forces were recorded and analyzed. A large amount of thread bone debris was observed in the control group. In the experimental group, however, only a small amount of thread bone debris was observed; the widths of screw paths were larger than those in the control group and gradually increased from screw tips in the direction of loading. The vertebral bodies had an average shear strength of [Formula: see text][Formula: see text]N, and the corresponding average screw displacement was [Formula: see text][Formula: see text]mm. Linear regression analysis showed that the shear strength had a significant positive correlation with vertebral bone mineral density (BMD) (r[Formula: see text][Formula: see text][Formula: see text]0.958, P[Formula: see text][Formula: see text][Formula: see text]0.01), while the screw displacement had a significant negative correlation with vertebral BMD (r[Formula: see text][Formula: see text][Formula: see text]−0.933, P[Formula: see text][Formula: see text][Formula: see text]0.01). No significant difference in bone density was found between the destruction and the control groups (P[Formula: see text][Formula: see text][Formula: see text]0.05); the difference in maximum axial pullout strength between the destruction and the control groups was significant (P[Formula: see text][Formula: see text][Formula: see text]0.01). These results indicated that vertebral BMD positively correlated with the maximum shear force and negatively with the screw displacement. It is important that the corrective strength for spinal deformity may affect bone–screw interface in anterior thoracolumbar surgery.

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.

Prediction of vertebral failure load by using x-ray vector radiographic imaging.

To examine whether x-ray vector radiographic (XVR) parameters could predict the biomechanically determined vertebral failure load. Local institutional review boards approved the study and donors provided written informed consent before death. Twelve thoracic vertebral bodies were removed from three human cadavers and embedded in resin. XVR measurements were performed by using a Talbot-Lau grating interferometer with the beam direction in anterior-posterior and lateral direction. The mean anisotropy and the mean local average scattering power were calculated for a region of interest within each vertebra. Trabecular bone mineral density (BMD) was determined in each vertebra by using a clinical multidetector computed tomographic scanner. Failure load of the vertebral bodies was determined from destructive biomechanical tests. Statistical analyses were performed with statistical software with a two-sided Pvalue of .05 to calculate Pearson correlation coefficients and multiple regression model. Statistically significant correlations (P < .05) for failure load with XVR parameters in the lateral direction (r = -0.84 and 0.68 for anisotropy and local average scattering power, respectively) and for failure load and anisotropy in anteroposterior direction (r = -0.65) were found. A multiple regression model showed that the combination of the local average scattering power in lateral direction and BMD predicted failure load significantly better than BMD alone (adjusted R = 0.88 compared with 0.78, respectively; P < .001). The study results imply that XVR can improve the prediction of osteoporosis.

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

BackgroundPathologic vertebral fractures are associated with intractable pain, loss of function and high morbidity in patients with metastatic spine disease. However, the failure mechanisms of vertebrae with lytic defects and the failed vertebrae's ability to retain load carrying capacity remain unclear. MethodsEighteen human thoracic and lumbar vertebrae with simulated uncontained bone defects were tested under compression-bending loads to failure. Failure was defined as 50% reduction in vertebral body height. The vertebrae were allowed to recover under load and re-tested to failure using the initial criteria. Repeated measure ANOVA was used to test for changes in strength and stiffness parameters. FindingsVertebral failure occurred via buckling and fracture of the cortex around the defect, followed by collapse of the defect region. Compared to the intact vertebrae, the failed vertebrae exhibited a significant loss in compressive strength (59%, p<0.001), stiffness (53%, p<0.05) and flexion (70%, p<0.01) strength. Significant reduction in anterior-posterior shear (strength (63%, p<0.01) and stiffness (67%, p<0.01)) and lateral bending strength (134%, p<0.05) were similarly recorded. In the intact vertebrae, apart from flexion strength (r2=0.63), both compressive and anterior-posterior shear strengths were weakly correlated with their stiffness parameters (r2=0.24 and r2=0.31). By contrast, in the failed vertebrae, these parameters were strongly correlated, (r2=0.91, r2=0.86, and r2=0.92, p<0.001 respectively). InterpretationFailure of the vertebral cortex at the defect site dominated the initiation and progression of vertebral failure with the vertebrae failing via a consolidation process of the vertebral bone. Once failed, the vertebrae showed remarkable loss of load carrying capacity.

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. 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. 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. 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. 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.