Cortical bone loss at the tibia in postmenopausal women with osteoporosis is associated with incident non-vertebral fractures: Results of a randomized controlled ancillary study of HORIZON

BACKGROUND In postmenopausal women, yearly intravenous zoledronate (ZOL) compared to placebo (PLB) significantly increased bone mineral density (BMD) at lumbar spine (LS), femoral neck (FN), and total hip (TH) and decreased fracture risk. The effects of ZOL on BMD at the tibial epiphysis (T-EPI) and diaphysis (T-DIA) are unknown. METHODS A randomized controlled ancillary study of the HORIZON trial was conducted at the Department of Osteoporosis of the University Hospital of Berne, Switzerland. Women with ≥1 follow-up DXA measurement who had received ≥1 dose of either ZOL (n=55) or PLB (n=55) were included. BMD was measured at LS, FN, TH, T-EPI, and T-DIA at baseline, 6, 12, 24, and 36 months. Morphometric vertebral fractures were assessed. Incident clinical fractures were recorded as adverse events. RESULTS Baseline characteristics were comparable with those in HORIZON and between groups. After 36 months, BMD was significantly higher in women treated with ZOL vs. PLB at LS, FN, TH, and T-EPI (+7.6%, +3.7%, +5.6%, and +5.5%, respectively, p<0.01 for all) but not T-DIA (+1.1%). The number of patients with ≥1 incident non-vertebral or morphometric fracture did not differ between groups (9 ZOL/11 PLB). Mean changes in BMD did not differ between groups with and without incident fracture, except that women with an incident non-vertebral fracture had significantly higher bone loss at predominantly cortical T-DIA (p=0.005). CONCLUSION ZOL was significantly superior to PLB at T-EPI but not at T-DIA. Women with an incident non-vertebral fracture experienced bone loss at T-DIA.

Finite element analyses of human vertebral bodies embedded in polymethylmethalcrylate or loaded via the hyperelastic intervertebral disc models provide equivalent predictions of experimental strength

Quantitative computer tomography (QCT)-based finite element (FE) models of vertebral body provide better prediction of vertebral strength than dual energy X-ray absorptiometry. However, most models were validated against compression of vertebral bodies with endplates embedded in polymethylmethalcrylate (PMMA). Yet, loading being as important as bone density, the absence of intervertebral disc (IVD) affects the strength. Accordingly, the aim was to assess the strength predictions of the classic FE models (vertebral body embedded) against the in vitro and in silico strengths of vertebral bodies loaded via IVDs. High resolution peripheral QCT (HR-pQCT) were performed on 13 segments (T11/T12/L1). T11 and L1 were augmented with PMMA and the samples were tested under a 4° wedge compression until failure of T12. Specimen-specific model was generated for each T12 from the HR-pQCT data. Two FE sets were created: FE-PMMA refers to the classical vertebral body embedded model under axial compression; FE-IVD to their loading via hyperelastic IVD model under the wedge compression as conducted experimentally. Results showed that FE-PMMA models overestimated the experimental strength and their strength prediction was satisfactory considering the different experimental set-up. On the other hand, the FE-IVD models did not prove significantly better (Exp/FE-PMMA: R²=0.68; Exp/FE-IVD: R²=0.71, p=0.84). In conclusion, FE-PMMA correlates well with in vitro strength of human vertebral bodies loaded via real IVDs and FE-IVD with hyperelastic IVDs do not significantly improve this correlation. Therefore, it seems not worth adding the IVDs to vertebral body models until fully validated patient-specific IVD models become available.

Experimental Validation of Finite Element Analysis of Human Vertebral Collapse Under Large Compressive Strains

Osteoporosis-related vertebral fractures represent a major health problem in elderly populations. Such fractures can often only be diagnosed after a substantial deformation history of the vertebral body. Therefore, it remains a challenge for clinicians to distinguish between stable and progressive potentially harmful fractures. Accordingly, novel criteria for selection of the appropriate conservative or surgical treatment are urgently needed. Computer tomography-based finite element analysis is an increasingly accepted method to predict the quasi-static vertebral strength and to follow up this small strain property longitudinally in time. A recent development in constitutive modeling allows us to simulate strain localization and densification in trabecular bone under large compressive strains without mesh dependence. The aim of this work was to validate this recently developed constitutive model of trabecular bone for the prediction of strain localization and densification in the human vertebral body subjected to large compressive deformation. A custom-made stepwise loading device mounted in a high resolution peripheral computer tomography system was used to describe the progressive collapse of 13 human vertebrae under axial compression. Continuum finite element analyses of the 13 compression tests were realized and the zones of high volumetric strain were compared with the experiments. A fair qualitative correspondence of the strain localization zone between the experiment and finite element analysis was achieved in 9 out of 13 tests and significant correlations of the volumetric strains were obtained throughout the range of applied axial compression. Interestingly, the stepwise propagating localization zones in trabecular bone converged to the buckling locations in the cortical shell. While the adopted continuum finite element approach still suffers from several limitations, these encouraging preliminary results towardsthe prediction of extended vertebral collapse may help in assessing fracture stability in future work.

Clinical versus pre-clinical FE models for vertebral body strength predictions

The finite element analysis is an accepted method to predict vertebral body compressive strength. This study compares measurements obtained from in vitro tests with the ones from two different simulation models: clinical quantitative computer tomography (QCT) based homogenized finite element (hFE) models and pre-clinical high-resolution peripheral QCT-based (HR-pQCT) hFE models. About 37 vertebral body sections were prepared by removing end-plates and posterior elements, scanned with QCT (390/450μm voxel size) as well as HR-pQCT (82μm voxel size), and tested in compression up to failure. Non-linear viscous damage hFE models were created from QCT/HT-pQCT images and compared to experimental results based on stiffness and ultimate load. As expected, the predictability of QCT/HR-pQCT-based hFE models for both apparent stiffness (r2=0.685/0.801r2=0.685/0.801) and strength (r2=0.774/0.924r2=0.774/0.924) increased if a better image resolution was used. An analysis of the damage distribution showed similar damage locations for all cases. In conclusion, HR-pQCT-based hFE models increased the predictability considerably and do not need any tuning of input parameters. In contrast, QCT-based hFE models usually need some tuning but are clinically the only possible choice at the moment.

Image-based biomechanical assessment of vertebral body and intervertebral disc in the human lumbar spine

Life expectancy continuously increases but our society faces age-related conditions. Among musculoskeletal diseases, osteoporosis associated with risk of vertebral fracture and degenerative intervertebral disc (IVD) are painful pathologies responsible for tremendous healthcare costs. Hence, reliable diagnostic tools are necessary to plan a treatment or follow up its efficacy. Yet, radiographic and MRI techniques, respectively clinical standards for evaluation of bone strength and IVD degeneration, are unspecific and not objective. Increasingly used in biomedical engineering, CT-based finite element (FE) models constitute the state-of-art for vertebral strength prediction. However, as non-invasive biomechanical evaluation and personalised FE models of the IVD are not available, rigid boundary conditions (BCs) are applied on the FE models to avoid uncertainties of disc degeneration that might bias the predictions. Moreover, considering the impact of low back pain, the biomechanical status of the IVD is needed as a criterion for early disc degeneration. Thus, the first FE study focuses on two rigid BCs applied on the vertebral bodies during compression test of cadaver vertebral bodies, vertebral sections and PMMA embedding. The second FE study highlights the large influence of the intervertebral disc’s compliance on the vertebral strength, damage distribution and its initiation. The third study introduces a new protocol for normalisation of the IVD stiffness in compression, torsion and bending using MRI-based data to account for its morphology. In the last study, a new criterion (Otsu threshold) for disc degeneration based on quantitative MRI data (axial T2 map) is proposed. The results show that vertebral strength and damage distribution computed with rigid BCs are identical. Yet, large discrepancies in strength and damage localisation were observed when the vertebral bodies were loaded via IVDs. The normalisation protocol attenuated the effect of geometry on the IVD stiffnesses without complete suppression. Finally, the Otsu threshold computed in the posterior part of annulus fibrosus was related to the disc biomechanics and meet objectivity and simplicity required for a clinical application. In conclusion, the stiffness normalisation protocol necessary for consistent IVD comparisons and the relation found between degeneration, mechanical response of the IVD and Otsu threshold lead the way for non-invasive evaluation biomechanical status of the IVD. As the FE prediction of vertebral strength is largely influenced by the IVD conditions, this data could also improve the future FE models of osteoporotic vertebra.

Incidence and risk factors for early adjacent vertebral fractures after balloon kyphoplasty for osteoporotic fractures: analysis of the SWISSspine registry

PURPOSE The SWISSspine registry (SSR) was launched in 2005 to assess the safety and effectiveness of balloon kyphoplasty (BKP). In the meantime, repeated reports on high rates of adjacent vertebral fractures (ASF) after BKP of vertebral insufficiency fractures were published. The causes for ASF and their risk factors are still under debate. The purpose of this study was to report the incidence and potential risk factors of ASF within the SSR dataset. METHODS The SSR data points are collected perioperatively and during follow-ups, with surgeon- and patient-based information. All patients documented with a monosegmental osteoporotic vertebral insufficiency fracture between March 2005 and May 2012 were included in the study. The incidence of ASF, significant associations with co-variates (patient age, gender, fracture location, cement volume, preoperative segmental kyphosis, extent of kyphosis correction, and individual co-morbidities) and influence on quality of life (EQ-5D) and back pain (VAS) were analyzed. RESULTS A total of 375 patients with a mean follow-up of 3.6 months was included. ASF were found in 9.9 % (n = 37) and occurred on average 2.8 months postoperatively. Preoperative segmental kyphosis >30° (p = 0.026), and rheumatoid arthritis (p = 0.038) and cardiovascular disease (p = 0.047) were significantly associated with ASF. Furthermore, patients with ASF had significantly higher back pain at the final follow-up (p = 0.001). No further significant associations between the studied co-variates and ASF were seen in the adjusted analysis. CONCLUSIONS The findings suggest that patients with a preoperative segmental kyphosis >30° or patients with co-morbidities like rheumatoid arthritis and a cardiovascular disease are at high risk of ASF within 6 months after the index surgery. In case of an ASF event, back pain levels are significantly increased. LEVEL OF EVIDENCE IV.

Reconstruction of 3D Vertebral Models from a Single 2D Lateral Fluoroscopic Image

Accurate three-dimensional (3D) models of lumbar vertebrae are required for image-based 3D kinematics analysis. MRI or CT datasets are frequently used to derive 3D models but have the disadvantages that they are expensive, time-consuming or involving ionizing radiation (e.g., CT acquisition). In this chapter, we present an alternative technique that can reconstruct a scaled 3D lumbar vertebral model from a single two-dimensional (2D) lateral fluoroscopic image and a statistical shape model. Cadaveric studies are conducted to verify the reconstruction accuracy by comparing the surface models reconstructed from a single lateral fluoroscopic image to the ground truth data from 3D CT segmentation. A mean reconstruction error between 0.7 and 1.4 mm was found.

Finite element analysis predicts experimental failure patterns in vertebral bodies loaded via intervertebral discs up to large deformation

Vertebral compression fracture is a common medical problem in osteoporotic individuals. The quantitative computed tomography (QCT)-based finite element (FE) method may be used to predict vertebral strength in vivo, but needs to be validated with experimental tests. The aim of this study was to validate a nonlinear anatomy specific QCT-based FE model by using a novel testing setup. Thirty-seven human thoracolumbar vertebral bone slices were prepared by removing cortical endplates and posterior elements. The slices were scanned with QCT and the volumetric bone mineral density (vBMD) was computed with the standard clinical approach. A novel experimental setup was designed to induce a realistic failure in the vertebral slices in vitro. Rotation of the loading plate was allowed by means of a ball joint. To minimize device compliance, the specimen deformation was measured directly on the loading plate with three sensors. A nonlinear FE model was generated from the calibrated QCT images and computed vertebral stiffness and strength were compared to those measured during the experiments. In agreement with clinical observations, most of the vertebrae underwent an anterior wedge-shape fracture. As expected, the FE method predicted both stiffness and strength better than vBMD (R2 improved from 0.27 to 0.49 and from 0.34 to 0.79, respectively). Despite the lack of fitting parameters, the linear regression of the FE prediction for strength was close to the 1:1 relation (slope and intercept close to one (0.86 kN) and to zero (0.72 kN), respectively). In conclusion, a nonlinear FE model was successfully validated through a novel experimental technique for generating wedge-shape fractures in human thoracolumbar vertebrae.

Vertebral artery stenosis in the Basilar Artery International Cooperation Study (BASICS): prevalence and outcome

We assessed the prevalence of vertebral artery (VA) stenosis or occlusion and its influence on outcome in patients with acute basilar artery occlusion (BAO). We studied 141 patients with acute BAO enrolled in the Basilar Artery International Cooperation Study (BASICS) registry of whom baseline CT angiography (CTA) of the intracranial VAs was available. In 72 patients an additional CTA of the extracranial VAs was available. Adjusted risk ratios (aRRs) for death and poor outcome, defined as a modified Rankin Scale score ≥4, were calculated with Poisson regression in relation to VA occlusion, VA occlusion or stenosis ≥50 %, and bilateral VA occlusion. Sixty-six of 141 (47 %) patients had uni- or bilateral intracranial VA occlusion or stenosis ≥50 %. Of the 72 patients with intra- and extracranial CTA, 46 (64 %) had uni- or bilateral VA occlusion or stenosis ≥50 % and 9 (12 %) had bilateral VA occlusion. Overall, VA occlusion or stenosis ≥50 % was not associated with the risk of poor outcome. Patients with intra- and extracranial CTA and bilateral VA occlusion had a higher risk of poor outcome than patients without bilateral VA occlusion (aRR, 1.23; 95 % CI 1.02-1.50). The risk of death did not depend on the presence of unilateral or bilateral VA occlusion or stenosis ≥50 %. In conclusion, in patients with acute BAO, unilateral VA occlusion or stenosis ≥50 % is frequent, but not associated with an increased risk of poor outcome or death. Patients with BAO and bilateral VA occlusion have a slightly increased risk of poor outcome.

The Initial Slope of the Variogram, Foundation of the Trabecular Bone Score, Is Not or Poorly Associated With Vertebral Strength

Trabecular bone score (TBS) rests on the textural analysis of DXA to reflect the decay in trabecular structure characterising osteoporosis. Yet, its discriminative power in fracture studies remains incomprehensible as prior biomechanical tests found no correlation with vertebral strength. To verify this result possibly due to an unrealistic set-up and to cover a wide range of loading scenarios, the data from three previous biomechanical studies using different experimental settings was used. They involved the compressive failure of 62 human lumbar vertebrae loaded 1) via intervertebral discs to mimic the in vivo situation (“full vertebra”), 2) via the classical endplate embedding (“vertebral body”) or 3) via a ball joint to induce anterior wedge failure (“vertebral section”). HR-pQCT scans acquired prior testing were used to simulate anterior-posterior DXA from which areal bone mineral density (aBMD) and the initial slope of the variogram (ISV), the early definition of TBS, were evaluated. Finally, the relation of aBMD and ISV with failure load (Fexp) and apparent failure stress (σexp) was assessed and their relative contribution to a multi-linear model was quantified via ANOVA. We found that, unlike aBMD, ISV did not significantly correlate with Fexp and σexp, except for the “vertebral body” case (r2 = 0.396, p = 0.028). Aside from the “vertebra section” set-up where it explained only 6.4% of σexp (p = 0.037), it brought no significant improvement to aBMD. These results indicate that ISV, a replica of TBS, is a poor surrogate for vertebral strength no matter the testing set-up, which supports the prior observations and raises a fortiori the question of the deterministic factors underlying the statistical relationship between TBS and vertebral fracture risk.

Determination of vertebral range of motion using inertial measurement units in 27 Franches-Montagnes stallions and comparison between conditions and with a mixed population

REASONS FOR PERFORMING STUDY: The diagnosis of equine back disorders is challenging. Objectively determining movement of the vertebral column may therefore be of value in a clinical setting. OBJECTIVES: To establish whether surface-mounted inertial measurement units (IMUs) can be used to establish normal values for range of motion (ROM) of the vertebral column in a uniform population of horses trotting under different conditions. STUDY DESIGN: Vertebral ROM was established in Franches-Montagnes stallions and a general population of horses and the variability in measurements compared between the two groups. Repeatability and the influence of specific exercise condition (on ROM) were assessed. Finally, attempts were made to explain the findings of the study through the evaluation of factors that might influence ROM. METHODS: Dorsoventral (DV) and mediolateral (ML) vertebral ROM was measured at a trot under different exercise conditions in 27 Franches-Montagnes stallions and six general population horses using IMUs distributed over the vertebral column. RESULTS: Variability in the ROM measurements was significantly higher for general population horses than for Franches-Montagnes stallions (both DV and ML ROM). Repeatability was strong to very strong for DV measurements and moderate for ML measurements. Trotting under saddle significantly reduced the ROM, with sitting trot resulting in a significantly lower ROM than rising trot. Age is unlikely to explain the low variability in vertebral ROM recorded in the Franches-Montagnes horses, while this may be associated with conformational factors. CONCLUSIONS: It was possible to establish a normal vertebral ROM for a group of Franches-Montagnes stallions. While within-breed variation was low in this population, further studies are necessary to determine variation in vertebral ROM for other breeds and to assess their utility for diagnosis of equine back disorders.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p ≤ 0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.