Assessment of incident spine and hip fractures in women and men using finite element analysis of CT scans

Discrimination of osteoporosis-related vertebral fractures by DXA-derived 3D measurements: a retrospective case-control study

A retrospective case-control study assessing the association of DXA-derived 3D measurements with osteoporosis-related vertebral fractures was performed. Trabecular volumetric bone mineral density was the measurement that best discriminates between fracture and control groups.

INTRODUCTION

The aim of the present study was to evaluate the association of DXA-derived 3D measurements at the lumbar spine with osteoporosis-related vertebral fractures.

METHODS

We retrospectively analyzed a database of 74 postmenopausal women: 37 subjects with incident vertebral fractures and 37 age-matched controls without any type of fracture. DXA scans at the lumbar spine were acquired at baseline (i.e., before the fracture event for subjects in the fracture group), and areal bone mineral density (aBMD) was measured. DXA-derived 3D measurements, such as volumetric BMD (vBMD), were assessed using a DXA-based 3D modeling software (3D-SHAPER). vBMD was computed at the trabecular, cortical, and integral bone. Cortical thickness and cortical surface BMD were also measured. Differences in DXA-derived measurements between fracture and control groups were evaluated using unpaired t test. Odds ratio (OR) and area under the receiver operating curve (AUC) were also computed. Subgroup analyses according to fractured vertebra were performed.

RESULTS

aBMD of fracture group was 9.3% lower compared with control group (p < 0.01); a higher difference was found for trabecular vBMD in the vertebral body (- 16.1%, p < 0.001). Trabecular vBMD was the measurement that best discriminates between fracture and control groups, with an AUC of 0.733, against 0.682 for aBMD. Overall, similar findings were observed within the subgroup analyses. The L1 vertebral fractures subgroup had the highest AUC at trabecular vBMD (0.827), against aBMD (0.758).

CONCLUSION

This study showed the ability of cortical and trabecular measurements from DXA-derived 3D models to discriminate between fracture and control groups. Large cohorts need to be analyzed to determine if these measurements could improve fracture risk prediction in clinical practice.

Longitudinal changes in hip geometry in relation to the final menstrual period: Study of Women's Health Across the Nation (SWAN)

BACKGROUND

In SWAN, we showed that accelerated loss of bone mineral density (BMD) begins 1 year before the final menstrual period (FMP) to 2 years after the FMP and slows thereafter. However, the risk of fracture depends on both BMD and bone geometry. The hip structural analysis (HSA) measures important geometric properties of bone. Changes in HSA parameters across the menopausal transition have not been previously assessed.

METHODS

The current analysis uses data from SWAN, 5 years before to 5 years after FMP (N = 900, Age (mean(SD)) = 46.85(2.60), 44% White). HSA parameters at the femoral narrow neck were obtained from 2D DXA scans and normalized to baseline values. FMP was determined from annual interviews. Changes in HSA were assessed over 3 periods, 5 to 2 years before FMP (pre-transmenopausal), 2 years before to 1 years after FMP (transmenopausal), 1 to 5 years after FMP (postmenopausal). Mixed linear models with random slopes were used to estimate the rate of change in HSA parameters relative to FMP.

RESULTS

Loss of BMD, cross-sectional area (CSA), and section modulus (SM) and increases in outer diameter (OD) were greatest in the transmenopausal period (p for all<0.05). Changes continued in the postmenopausal period but were not statistically significant. The cumulative percentage changes over 10 years in BMD (-10.67%), CSA (-9.01), SM (-7.03) and OD (+1.95) were statistically significant.

CONCLUSION

Changes in hip geometry across the menopause transition parallel changes in BMD and provide insight into mechanisms that may increase risk of fragility fracture.

In vitro injection of osteoporotic cadaveric femurs with a triphasic calcium-based implant confers immediate biomechanical integrity

Current pharmaceutical therapies can reduce hip fractures by up to 50%, but compliance to treatment is low and therapies take up to 18 months to reduce risk. Thus, alternative or complementary approaches to reduce the risk of hip fracture are needed. The AGN1 local osteo-enhancement procedure (LOEP) is one such alternative approach, as it is designed to locally replace bone lost due to osteoporosis and provide immediate biomechanical benefit. This in vitro study evaluated the initial biomechanical impact of this treatment on human cadaveric femurs. We obtained 45 pairs of cadaveric femurs from women aged 77.8 ± 8.8 years. One femur of each pair was treated, while the contralateral femur served as an untreated control. Treatment included debridement, irrigation/suction, and injection of a triphasic calcium-based implant (AGN1). Mechanical testing of the femora was performed in a sideways fall configuration 24 h after treatment. Of the 45 pairs, 4 had normal, 16 osteopenic, and 25 osteoporotic BMD T-scores. Altogether, treatment increased failure load on average by 20.5% (p < 0.0001). In the subset of osteoporotic femurs, treatment increased failure load by 26% and work to failure by 45% (p < 0.01 for both). Treatment did not significantly affect stiffness in any group. These findings provide evidence that local delivery of the triphasic calcium-based implant in the proximal femur is technically feasible and provides immediate biomechanical benefit. Our results provide strong rationale for additional studies investigating the utility of this approach for reducing the risk of hip fracture. © 2019 The Authors. Journal of Orthopaedic Research® Published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society.

Prediction of lumbar vertebral body compressive strength of overweight and obese older adults using morphed subject-specific finite-element models to evaluate the effects of weight loss

BACKGROUND

Diet and exercise can promote weight loss in older adults; however, there is potential to increase fracture risk due to loss of bone mineral density (BMD) known to accompany weight loss. Weight loss effects on measures of bone quality and strength are currently unknown.

AIMS

The purpose of this study is to develop subject-specific finite-element (FE) models of the lumbar spine and study the effect of intentional weight loss on bone strength in a pilot data set.

METHODS

Computed tomography (CT) scans of the lumbar spine of 30 overweight and obese (mean BMI = 29.7 ± 3.9 kg/m²), older adults (mean age = 65.9 ± 4.6 years) undergoing an 18-month intentional weight loss intervention were obtained at baseline and post-intervention. Measures of volumetric BMD (vBMD) and variable cortical thickness were derived from each subject CT scan. Development of the subject-specific FE models of the lumbar spine involved model morphing techniques to accelerate the development of the models. vBMD-derived material properties and cortical thickness measures were directly mapped to baseline and post-intervention models. Bone strength was estimated through simulation of a quasi-static uniaxial compression test.

RESULTS

From baseline to 18-month post-weight loss intervention, there were statistically significant decreases in estimated bone strength (6.5% decrease; p < 0.05). Adjusting for baseline bone measures and gender revealed no statistically significant correlations between weight change and change in vBMD, cortical thickness, or bone strength.

CONCLUSION

Integration of CT-based measures and FE models with conventional areal BMD can improve the understanding of the effects of intentional weight loss on bone health.

Hip load capacity cut-points for Astronaut Skeletal Health NASA Finite Element Strength Task Group Recommendations

Concerns raised at a 2010 Bone Summit held for National Aeronautics and Space Administration Johnson Space Center led experts in finite element (FE) modeling for hip fracture prediction to propose including hip load capacity in the standards for astronaut skeletal health. The current standards for bone are based upon areal bone mineral density (aBMD) measurements by dual X-ray absorptiometry (DXA) and an adaptation of aBMD cut-points for fragility fractures. Task Group members recommended (i) a minimum permissible outcome limit (POL) for post-mission hip bone load capacity, (ii) use of FE hip load capacity to further screen applicants to astronaut corps, (iii) a minimum pre-flight standard for a second long-duration mission, and (iv) a method for assessing which post-mission physical activities might increase an astronaut’s risk for fracture after return. QCT-FE models of eight astronaut were analyzed using nonlinear single-limb stance (NLS) and posterolateral fall (NLF) loading configurations. QCT data from the Age Gene/Environment Susceptibility (AGES) Reykjavik cohort and the Rochester Epidemiology Project were analyzed using identical modeling procedures. The 75^th percentile of NLS hip load capacity for fractured elderly males of the AGES cohort (9537N) was selected as a post-mission POL. The NLF model, in combination with a Probabilistic Risk Assessment tool, was used to assess the likelihood of exceeding the hip load capacity during post-flight activities. There was no recommendation to replace the current DXA-based standards. However, FE estimation of hip load capacity appeared more meaningful for younger, physically active astronauts and was recommended to supplement aBMD cut-points.

Biofidelic finite element models for accurately classifying hip fracture in a retrospective clinical study of elderly women from the AGES Reykjavik cohort

Clinical retrospective studies have only reported limited improvements in hip fracture classification accuracy using finite element (FE) models compared to conventional areal bone mineral density (aBMD) measurements. A possible explanation is that state-of-the-art quasi-static models do not estimate patient-specific loads. A novel FE modeling technique was developed to improve the biofidelity of simulated impact loading from sideways falling. This included surrogate models of the pelvis, lower extremities, and soft tissue that were morphed based on subject anthropometrics. Hip fracture prediction models based on aBMD and FE measurements were compared in a retrospective study of 254 elderly female subjects from the AGES-Reykjavik study. Subject fragility ratio (FR) was defined as the ratio between the ultimate forces of paired biofidelic models, one with linear elastic and the other with non-linear stress-strain relationships in the proximal femur. The expected end-point value (EEV) was defined as the FR weighted by the probability of one sideways fall over five years, based on self-reported fall frequency at baseline. The change in maximum volumetric strain (ΔMVS) on the surface of the femoral neck was calculated between time of ultimate femur force and 90% post-ultimate force in order to assess the extent of tensile tissue damage present in non-linear models. After age-adjusted logistic regression, the area under the receiver-operator curve (AUC) was highest for ΔMVS (0.72), followed by FR (0.71), aBMD (0.70), and EEV (0.67), however the differences between FEA and aBMD based prediction models were not deemed statistically significant. When subjects with no history of falling were excluded from the analysis, thus artificially assuming that falls were known a priori with no uncertainty, a statistically significant difference in AUC was detected between ΔMVS (0.85), and aBMD (0.74). Multivariable linear regression suggested that the variance in maximum elastic femur force was best explained by femoral head radius, pelvis width, and soft tissue thickness (R2 = 0.79; RMSE = 0.46 kN; p < 0.005). Weighting the hip fracture prediction models based on self-reported fall frequency did not improve the models' sensitivity, however excluding non-fallers lead to significant differences between aBMD and FE based models. These findings suggest that an accurate assessment of fall probability is necessary for accurately identifying individuals predisposed to hip fracture.

Prediction of incident vertebral fracture using CT-based finite element analysis

Prior studies show vertebral strength from computed tomography-based finite element analysis may be associated with vertebral fracture risk. We found vertebral strength had a strong association with new vertebral fractures, suggesting that vertebral strength measures identify those at risk for vertebral fracture and may be a useful clinical tool.

INTRODUCTION

We aimed to determine the association between vertebral strength by quantitative computed tomography (CT)-based finite element analysis (FEA) and incident vertebral fracture (VF). In addition, we examined sensitivity and specificity of previously proposed diagnostic thresholds for fragile bone strength and low BMD in predicting VF.

METHODS

In a case-control study, 26 incident VF cases (13 men, 13 women) and 62 age- and sex-matched controls aged 50 to 85 years were selected from the Framingham multi-detector computed tomography cohort. Vertebral compressive strength, integral vBMD, trabecular vBMD, CT-based BMC, and CT-based aBMD were measured from CT scans of the lumbar spine.

RESULTS

Lower vertebral strength at baseline was associated with an increased risk of new or worsening VF after adjusting for age, BMI, and prevalent VF status (odds ratio (OR) = 5.2 per 1 SD decrease, 95% CI 1.3-19.8). Area under receiver operating characteristic (ROC) curve comparisons revealed that vertebral strength better predicted incident VF than CT-based aBMD (AUC = 0.804 vs. 0.715, p = 0.05) but was not better than integral vBMD (AUC = 0.815) or CT-based BMC (AUC = 0.794). Additionally, proposed fragile bone strength thresholds trended toward better sensitivity for identifying VF than that of aBMD-classified osteoporosis (0.46 vs. 0.23, p = 0.09).

CONCLUSION

This study shows an association between vertebral strength measures and incident vertebral fracture in men and women. Though limited by a small sample size, our findings also suggest that bone strength estimates by CT-based FEA provide equivalent or better ability to predict incident vertebral fracture compared to CT-based aBMD. Our study confirms that CT-based estimates of vertebral strength from FEA are useful for identifying patients who are at high risk for vertebral fracture.

Mechanical behavior of metastatic femurs through patient-specific computational models accounting for bone-metastasis interaction

This paper proposes a computational model based on a finite-element formulation for describing the mechanical behavior of femurs affected by metastatic lesions. A novel geometric/constitutive description is introduced by modelling healthy bone and metastases via a linearly poroelastic constitutive approach. A Gaussian-shaped graded transition of material properties between healthy and metastatic tissues is prescribed, in order to account for the bone-metastasis interaction. Loading-induced failure processes are simulated by implementing a progressive damage procedure, formulated via a quasi-static displacement-driven incremental approach, and considering both a stress- and a strain-based failure criterion. By addressing a real clinical case, left and right patient-specific femur models are geometrically reconstructed via an ad-hoc imaging procedure and embedding multiple distributions of metastatic lesions along femurs. Significant differences in fracture loads, fracture mechanisms, and damage patterns, are highlighted by comparing the proposed constitutive description with a purely elastic formulation, where the metastasis is treated as a pseudo-healthy tissue or as a void region. Proposed constitutive description allows to capture stress/strain localization mechanisms within the metastatic tissue, revealing the model capability in describing possible strain-induced mechano-biological stimuli driving onset and evolution of the lesion. The proposed approach opens towards the definition of effective computational strategies for supporting clinical decision and treatments regarding metastatic femurs, contributing also to overcome some limitations of actual standards and procedures.

Effect of Statistically Iterative Image Reconstruction on Vertebral Bone Strength Prediction Using Bone Mineral Density and Finite Element Modeling: A Preliminary Study

Statistical iterative reconstruction (SIR) using multidetector computed tomography (MDCT) is a promising alternative to standard filtered back projection (FBP), because of lower noise generation while maintaining image quality. Hence, we investigated the feasibility of SIR in predicting MDCT-based bone mineral density (BMD) and vertebral bone strength from finite element (FE) analysis. The BMD and FE-predicted bone strength derived from MDCT images reconstructed using standard FBP (FFBP) and SIR with (FSIR) and without regularization (FSIRB0) were validated against experimental failure loads (Fexp). Statistical iterative reconstruction produced the best quality images with regard to noise, signal-to-noise ratio, and contrast-to-noise ratio. Fexp significantly correlated with FFBP, FSIR, and FSIRB0. FFBP had a significant correlation with FSIRB0 and FSIR. The BMD derived from FBP, SIRB0, and SIR were significantly correlated. Effects of regularization should be further investigated with FE and BMD analysis to allow for an optimal iterative reconstruction algorithm to be implemented in an in vivo scenario.

Patient-Specific Phantomless Estimation of Bone Mineral Density and Its Effects on Finite Element Analysis Results: A Feasibility Study

Objectives

This study proposes a regression model for the phantomless Hounsfield units (HU) to bone mineral density (BMD) conversion including patient physical factors and analyzes the accuracy of the estimated BMD values.

Methods

The HU values, BMDs, circumferences of the body, and cross-sectional areas of bone were measured from 39 quantitative computed tomography images of L2 vertebrae and hips. Then, the phantomless HU-to-BMD conversion was derived using a multiple linear regression model. For the statistical analysis, the correlation between the estimated BMD values and the reference BMD values was evaluated using Pearson's correlation test. Voxelwise BMD and finite element analysis (FEA) results were analyzed in terms of root-mean-square error (RMSE) and strain energy density, respectively.

Results

The HU values and circumferences were statistically significant (p < 0.05) for the lumbar spine, whereas only the HU values were statistically significant (p < 0.05) for the proximal femur. The BMD values estimated using the proposed HU-to-BMD conversion were significantly correlated with those measured using the reference phantom: Pearson's correlation coefficients of 0.998 and 0.984 for the lumbar spine and proximal femur, respectively. The RMSEs of the estimated BMD values for the lumbar spine and hip were 4.26 ± 0.60 (mg/cc) and 8.35 ± 0.57 (mg/cc), respectively. The errors of total strain energy were 1.06% and 0.91%, respectively.

Conclusions

The proposed phantomless HU-to-BMD conversion demonstrates the potential of precisely estimating BMD values from CT images without the reference phantom and being utilized as a viable tool for FEA-based quantitative assessment using routine CT images.

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

PURPOSE

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.

MATERIALS AND METHODS

An institutional review board (IRB) approval was obtained for this retrospective study. A total of 16 patients 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.

RESULTS

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.

CONCLUSION

The use of MDCT enables a dose reduction of at least 75% compared to standard-dose for an adequate prediction of vertebral failure load based on non-invasive FE analysis.

Fracture Prediction by Computed Tomography and Finite Element Analysis: Current and Future Perspectives

PURPOSE OF REVIEW

This review critiques the ability of CT-based methods to predict incident hip and vertebral fractures.

RECENT FINDINGS

CT-based techniques with concurrent calibration all show strong associations with incident hip and vertebral fracture, predicting hip and vertebral fractures as well as, and sometimes better than, dual-energy X-ray absorptiometry areal biomass density (DXA aBMD). There is growing evidence for use of routine CT scans for bone health assessment. CT-based techniques provide a robust approach for osteoporosis diagnosis and fracture prediction. It remains to be seen if further technical advances will improve fracture prediction compared to DXA aBMD. Future work should include more standardization in CT analyses, establishment of treatment intervention thresholds, and more studies to determine whether routine CT scans can be efficiently used to expand the number of individuals who undergo evaluation for fracture risk.

A novel sideways fall simulator to study hip fractures ex vivo

Falls to the side are the leading cause of hip fractures in the elderly. The load that a person experiences during a fall cannot be measured with volunteers for ethical reasons. To evaluate injurious loads, while considering relevant energy input and body posture for a sideways fall, a subject-specific cadaveric impact experiment was developed. Full cadaveric femur-pelvis constructs (N = 2) were embedded in surrogate soft tissue material and attached to metallic surrogate lower limbs. The specimens were then subjected to an inverted pendulum motion, simulating a fall to the side with an impact to the greater trochanter. The load at the ground and the deformation of the pelvis were evaluated using a 6-axis force transducer and two high-speed cameras. Post-test, a trauma surgeon (PG) evaluated specimen injuries. Peak ground contact forces were 7132 N and 5641 N for the fractured and non-fractured specimen, respectively. We observed a cervical fracture of the femur in one specimen and no injuries in a second specimen, showing that the developed protocol can be used to differentiate between specimens at high and low fracture risk.

Mapping anisotropy improves QCT-based finite element estimation of hip strength in pooled stance and side-fall load configurations

Hip fractures are one of the most severe consequences of osteoporosis. Compared to the clinical standard of DXA-based aBMD at the femoral neck, QCT-based FEA delivers a better surrogate of femoral strength and gains acceptance for the calculation of hip fracture risk when a CT reconstruction is available. Isotropic, homogenised voxel-based, finite element (hvFE) models are widely used to estimate femoral strength in cross-sectional and longitudinal clinical studies. However, fabric anisotropy is a classical feature of the architecture of the proximal femur and the second determinant of the homogenised mechanical properties of trabecular bone. Due to the limited resolution, fabric anisotropy cannot be derived from clinical CT reconstructions. Alternatively, fabric anisotropy can be extracted from HR-pQCT images of cadaveric femora. In this study, fabric anisotropy from HR-pQCT images was mapped onto QCT-based hvFE models of 71 human proximal femora for which both HR-pQCT and QCT images were available. Stiffness and ultimate load computed from anisotropic hvFE models were compared with previous biomechanical tests in both stance and side-fall configurations. The influence of using the femur-specific versus a mean fabric distribution on the hvFE predictions was assessed. Femur-specific and mean fabric enhance the prediction of experimental ultimate force for the pooled, i.e. stance and side-fall, (isotropic: r²=0.81, femur-specific fabric: r²=0.88, mean fabric: r²=0.86,p<0.001) but not for the individual configurations. Fabric anisotropy significantly improves bone strength prediction for the pooled configurations, and mapped fabric provides a comparable prediction to true fabric. The mapping of fabric anisotropy is therefore expected to help generate more accurate QCT-based hvFE models of the proximal femur for personalised or multiple load configurations.

Osteoporosis and Hip Fracture Risk From Routine Computed Tomography Scans: The Fracture, Osteoporosis, and CT Utilization Study (FOCUS)

Methods now exist for analyzing previously taken clinical computed tomography (CT) scans to measure a dual-energy X-ray absorptiometry (DXA)-equivalent bone mineral density (BMD) at the hip and a finite element analysis-derived femoral strength. We assessed the efficacy of this "biomechanical CT" (BCT) approach for identifying patients at high risk of incident hip fracture in a large clinical setting. Using a case-cohort design sampled from 111,694 women and men aged 65 or older who had a prior hip CT scan, a DXA within 3 years of the CT, and no prior hip fracture, we compared those with subsequent hip fracture (n = 1959) with randomly selected sex-stratified controls (n = 1979) and analyzed their CT scans blinded to all other data. We found that the age-, race-, and body mass index (BMI)-adjusted hazard ratio (HR; per standard deviation) for femoral strength was significant before (women: HR = 2.8, 95% confidence interval [CI] 2.2-3.5; men: 2.8, 2.1-3.7) and after adjusting also for the (lowest) hip BMD T-score by BCT (women: 2.1, 1.4-3.2; men: 2.7, 1.6-4.6). The hazard ratio for the hip BMD T-score was similar between BCT and DXA for both sexes (women: 2.1, 1.8-2.5 BCT versus 2.1, 1.7-2.5 DXA; men: 2.8, 2.1-3.8 BCT versus 2.5, 2.0-3.2 DXA) and was higher than for the (lowest) spine/hip BMD T-score by DXA (women: 1.6, 1.4-1.9; men: 2.1, 1.6-2.7). Compared with the latter as a clinical-practice reference and using both femoral strength and the hip BMD T-score from BCT, sensitivity for predicting hip fracture was higher for BCT (women: 0.66 versus 0.59; men: 0.56 versus 0.48), with comparable respective specificity (women: 0.66 versus 0.67; men: 0.76 versus 0.78). We conclude that BCT analysis of previously acquired routine abdominal or pelvic CT scans is at least as effective as DXA testing for identifying patients at high risk of hip fracture. © 2018 American Society for Bone and Mineral Research.

Can patient-specific finite element models better predict fractures in metastatic bone disease than experienced clinicians?: Towards computational modelling in daily clinical practice

OBJECTIVES

In this prospective cohort study, we investigated whether patient-specific finite element (FE) models can identify patients at risk of a pathological femoral fracture resulting from metastatic bone disease, and compared these FE predictions with clinical assessments by experienced clinicians.

METHODS

A total of 39 patients with non-fractured femoral metastatic lesions who were irradiated for pain were included from three radiotherapy institutes. During follow-up, nine pathological fractures occurred in seven patients. Quantitative CT-based FE models were generated for all patients. Femoral failure load was calculated and compared between the fractured and non-fractured femurs. Due to inter-scanner differences, patients were analyzed separately for the three institutes. In addition, the FE-based predictions were compared with fracture risk assessments by experienced clinicians.

RESULTS

In institute 1, median failure load was significantly lower for patients who sustained a fracture than for patients with no fractures. In institutes 2 and 3, the number of patients with a fracture was too low to make a clear distinction. Fracture locations were well predicted by the FE model when compared with post-fracture radiographs. The FE model was more accurate in identifying patients with a high fracture risk compared with experienced clinicians, with a sensitivity of 89% versus 0% to 33% for clinical assessments. Specificity was 79% for the FE models versus 84% to 95% for clinical assessments.

CONCLUSION

FE models can be a valuable tool to improve clinical fracture risk predictions in metastatic bone disease. Future work in a larger patient population should confirm the higher predictive power of FE models compared with current clinical guidelines.Cite this article: F. Eggermont, L. C. Derikx, N. Verdonschot, I. C. M. van der Geest, M. A. A. de Jong, A. Snyers, Y. M. van der Linden, E. Tanck. Can patient-specific finite element models better predict fractures in metastatic bone disease than experienced clinicians? Towards computational modelling in daily clinical practice. Bone Joint Res 2018;7:430-439. DOI: 10.1302/2046-3758.76.BJR-2017-0325.R2.

Are CT-Based Finite Element Model Predictions of Femoral Bone Strength Clinically Useful?

PURPOSE OF REVIEW

This study reviews the available literature to compare the accuracy of areal bone mineral density derived from dual X-ray absorptiometry (DXA-aBMD) and of subject-specific finite element models derived from quantitative computed tomography (QCT-SSFE) in predicting bone strength measured experimentally on cadaver bones, as well as their clinical accuracy both in terms of discrimination and prediction. Based on this information, some basic cost-effectiveness calculations are performed to explore the use of QCT-SSFE instead of DXA-aBMD in (a) clinical studies with femoral strength as endpoint, (b) predictor of the risk of hip fracture in low bone mass patients.

RECENT FINDINGS

Recent improvements involving the use of smooth-boundary meshes, better anatomical referencing for proximal-only scans, multiple side-fall directions, and refined boundary conditions increase the predictive accuracy of QCT-SSFE. If these improvements are adopted, QCT-SSFE is always preferable over DXA-aBMD in clinical studies with femoral strength as the endpoint, while it is not yet cost-effective as a hip fracture risk predictor, although pathways that combine both QCT-SSFE and DXA-aBMD are promising.

How does the femoral cortex depend on bone shape? A methodology for the joint analysis of surface texture and shape

In humans, there is clear evidence of an association between hip fracture risk and femoral neck bone mineral density, and some evidence of an association between fracture risk and the shape of the proximal femur. Here, we investigate whether the femoral cortex plays a role in these associations: do particular morphologies predispose to weaker cortices? To answer this question, we used cortical bone mapping to measure the distribution of cortical mass surface density (CMSD, mg/cm2) in a cohort of 125 females. Principal component analysis of the femoral surfaces identified three modes of shape variation accounting for 65% of the population variance. We then used statistical parametric mapping (SPM) to locate regions of the cortex where CMSD depends on shape, allowing for age. Our principal findings were increased CMSD with increased gracility over much of the proximal femur; and decreased CMSD at the superior femoral neck, coupled with increased CMSD at the calcar femorale, with increasing neck-shaft angle. In obtaining these results, we studied the role of spatial normalization in SPM, identifying systematic misregistration as a major impediment to the joint analysis of CMSD and shape. Through a series of experiments on synthetic data, we evaluated a number of registration methods for spatial normalization, concluding that only those predicated on an explicit set of homologous landmarks are suitable for this kind of analysis. The emergent methodology amounts to an extension of Geometric Morphometric Image Analysis to the domain of textured surfaces, alongside a protocol for labelling homologous landmarks in clinical CT scans of the human proximal femur.

Association of incident hip fracture with the estimated femoral strength by finite element analysis of DXA scans in the Osteoporotic Fractures in Men (MrOS) study

Finite element model can estimate bone strength better than BMD. This study used such a model to determine its association with hip fracture risk and found that the strength estimate provided limited improvement over the hip BMDs in predicting femoral neck (FN) fracture risk only.

INTRODUCTION

Bone fractures occur only when it is loaded beyond its ultimate strength. The goal of this study was to determine the association of femoral strength, as estimated by finite element (FE) analysis of DXA scans, with incident hip fracture as a single condition or with femoral neck (FN) and trochanter (TR) fractures separately in older men.

METHODS

This prospective case-cohort study included 91 FN and 64 TR fracture cases and a random sample of 500 men (14 had a hip fracture) from the Osteoporotic Fractures in Men study during a mean ± SD follow-up of 7.7 ± 2.2 years. We analysed the baseline DXA scans of the hip using a validated plane-stress, linear-elastic FE model of the proximal femur and estimated the femoral strength during a sideways fall.

RESULTS

The estimated strength was significantly (P < 0.05) associated with hip fracture independent of the TR and total hip (TH) BMDs but not FN BMD, and combining the strength with BMD did not improve the hip fracture prediction. The strength estimate was associated with FN fractures independent of the FN, TR and TH BMDs; the age-BMI-BMD adjusted hazard ratio (95% CI) per SD decrease of the strength was 1.68 (1.07-2.64), 2.38 (1.57, 3.61) and 2.04 (1.34, 3.11), respectively. This association with FN fracture was as strong as FN BMD (Harrell's C index for the strength 0.81 vs. FN BMD 0.81) and stronger than TR and TH BMDs (0.8 vs. 0.78 and 0.81 vs. 0.79). The strength's association with TR fracture was not independent of hip BMD.

CONCLUSIONS

Although the strength estimate provided additional information over the hip BMDs, its improvement in predictive ability over the hip BMDs was confined to FN fracture only and limited.

Biomechanical effects of metastasis in the osteoporotic lumbar spine: A Finite Element Analysis

BACKGROUND

Cancer patients are likely to undergo osteoporosis as consequence of hormone manipulation and/or chemotherapy. Little is known about possible increased risk of fracture in this population. The aim of this study was to describe the biomechanical effect of a metastatic lesion in an osteoporotic lumbar spine model.

METHODS

A finite element model of two spinal motion segments (L3-L5) was extracted from a previously developed L3-Sacrum model and used to analyze the effect of metastasis size and bone mineral density (BMD) on Vertebral bulge (VB) and Vertebral height (VH). VB and VH represent respectively radial and axial displacement and they have been correlated to burst fracture. A total of 6 scenarios were evaluated combining three metastasis sizes (no metastasis, 15% and 30% of the vertebral body) and two BMD conditions (normal BMD and osteoporosis).

RESULTS

15% metastasis increased VB and VH by 178% and 248%, respectively in normal BMD model; while VB and VH increased by 134% and 174% in osteoporotic model. 30% metastasis increased VB and VH by 88% and 109%, respectively, when compared to 15% metastasis in normal BMD model; while VB and VH increased by 59% and 74% in osteoporotic model.

CONCLUSION

A metastasis in the osteoporotic lumbar spine always leads to a higher risk of vertebral fracture. This risk increases with the size of the metastasis. Unexpectedly, an increment in metastasis size in the normal BMD spine produces a greater impact on vertebral stability compared to the osteoporotic spine.

Fully automated segmentation of a hip joint using the patient-specific optimal thresholding and watershed algorithm

BACKGROUND AND OBJECTIVE

Automated segmentation with high accuracy and speed is a prerequisite for FEA-based quantitative assessment with a large population. However, hip joint segmentation has remained challenging due to a narrow articular cartilage and thin cortical bone with a marked interindividual variance. To overcome this challenge, this paper proposes a fully automated segmentation method for a hip joint that uses the complementary characteristics between the thresholding technique and the watershed algorithm.

METHODS

Using the golden section method and load path algorithm, the proposed method first determines the patient-specific optimal threshold value that enables reliably separating a femur from a pelvis while removing cortical and trabecular bone in the femur at the minimum. This provides regional information on the femur. The watershed algorithm is then used to obtain boundary information on the femur. The proximal femur can be extracted by merging the complementary information on a target image.

RESULTS

For eight CT images, compared with the manual segmentation and other segmentation methods, the proposed method offers a high accuracy in terms of the dice overlap coefficient (97.24 ± 0.44%) and average surface distance (0.36 ± 0.07 mm) within a fast timeframe in terms of processing time per slice (1.25 ± 0.27 s). The proposed method also delivers structural behavior which is close to that of the manual segmentation with a small mean of average relative errors of the risk factor (4.99%).

CONCLUSION

The segmentation results show that, without the aid of a prerequisite dataset and users' manual intervention, the proposed method can segment a hip joint as fast as the simplified Kang (SK)-based automated segmentation, while maintaining the segmentation accuracy at a similar level of the snake-based semi-automated segmentation.

Automated DXA-based finite element analysis for hip fracture risk stratification: a cross-sectional study

Fracture risk indices (FRIs) generated from DXA-based finite element analysis were associated with hip fracture independent of FRAX score computed with femoral neck bone mineral density (BMD). Prospective studies are warranted to determine whether FRIs represent an improvement over BMD for predicting incident hip fractures.

INTRODUCTION

The study aims to examine the association between prior hip fracture and FRIs derived from automated finite element analysis (FEA) of DXA hip scans. Femoral neck, intertrochanteric, and subtrochanteric FRIs were calculated as the von Mises stress induced by a sideways fall divided by the bone yield stress over the specified region of interest (ROI).

METHODS

Using the Manitoba Bone Mineral Density Database, we selected women age ≥ 65 years with femoral neck T-scores below - 1 and no osteoporosis treatment. From this population, we identified 324 older women with hip fracture before DXA testing and a random sample of 658 non-fracture controls. FRIs were derived from the anonymized DXA scans. Logistic regression models were used to estimate odds ratios (ORs) and 95% confidence intervals (95% CIs) for the associations between FRIs (per SD increase) and hip fracture.

RESULTS

After adjusting for FRAX score (hip fracture with BMD), femoral neck FRI (OR 1.36, 95% CI 1.13, 1.64), intertrochanteric FRI (OR 1.81, 95% CI 1.44, 2.27), and subtrochanteric FRI (OR 2.09, 95% CI 1.68, 2.60) were associated with hip fracture. Intertrochanteric and subtrochanteric FRIs gave significantly higher c-statistics (all P ≤ 0.05) than femoral neck BMD. Subgroup analyses showed that all FRIs were more strongly associated with hip fracture in women who were younger and had higher body mass index (BMI) or non-osteoporotic BMD (all P _interaction < 0.1).

CONCLUSIONS

FRIs derived from DXA-based FEA were independently associated with prior hip fracture, suggesting that they could potentially improve hip fracture risk assessment.

Vertebral Fracture Risk in Diabetic Elderly Men: The MrOS Study

Type 2 diabetes (T2DM) is associated with a significant increase in risk of nonvertebral fractures, but information on risk of vertebral fractures (VFs) in subjects with T2DM, particularly among men, is lacking. Furthermore, it is not known whether spine bone mineral density (BMD) can predict the risk of VF in T2DM. We sought to examine the effect of diabetes status on prevalent and incident vertebral fracture, and to estimate the effect of lumbar spine BMD (areal and volumetric) as a risk factor for prevalent and incident morphometric vertebral fracture in T2DM (n = 875) and nondiabetic men (n = 4679). We used data from the Osteoporotic Fractures in Men (MrOS) Study, which enrolled men aged ≥65 years. Lumbar spine areal BMD (aBMD) was measured with dual-energy X-ray absorptiometry (DXA), and volumetric BMD (vBMD) by quantitative computed tomography (QCT). Prevalence (7.0% versus 7.7%) and incidence (4.4% versus 4.5%) of VFs were not higher in T2DM versus nondiabetic men. The risk of prevalent (OR, 1.05; 95% CI, 0.78 to 1.40) or incident vertebral-fracture (OR, 1.28; 95% CI, 0.81 to 2.00) was not higher in T2DM versus nondiabetic men in models adjusted for age, clinic site, race, BMI, and aBMD. Higher spine aBMD was associated with lower risk of prevalent VF in T2DM (OR, 0.55; 95% CI, 0.48 to 0.63) and nondiabetic men (OR, 0.66; 95% CI, 0.5 to 0.88) (p for interaction = 0.24) and of incident VF in T2DM (OR, 0.50; 95% CI, 0.41 to 0.60) and nondiabetic men (OR, 0.54; 95% CI, 0.33 to 0.88) (p for interaction = 0.77). Results were similar for vBMD. In conclusion, T2DM was not associated with higher prevalent or incident VF in older men, even after adjustment for BMI and BMD. Higher spine aBMD and vBMD are associated with lower prevalence and incidence of VF in T2DM as well as nondiabetic men. © 2017 American Society for Bone and Mineral Research.

The QUALYOR (QUalité Osseuse LYon Orléans) study: a new cohort for non invasive evaluation of bone quality in postmenopausal osteoporosis. Rationale and study design

The diagnostic performance of densitometry is inadequate. New techniques of non-invasive evaluation of bone quality may improve fracture risk prediction. Testing the value of these techniques is the goal of the QUALYOR cohort.

INTRODUCTION

The bone mineral density (BMD) of postmenopausal women who sustain osteoporotic fracture is generally above the World Health Organization definition for osteoporosis. Therefore, new approaches to improve the detection of women at high risk for fracture are warranted.

METHODS

We have designed and recruited a new cohort to assess the predictive value of several techniques to assess bone quality, including high-resolution peripheral quantitative computerized tomography (HRpQCT), hip QCT, calcaneus texture analysis, and biochemical markers. We have enrolled 1575 postmenopausal women, aged at least 50, with an areal BMD femoral neck or lumbar spine T-score between - 1.0 and - 3.0. Clinical risk factors for fracture have been collected along with serum and blood samples.

RESULTS

We describe the design of the QUALYOR study. Among these 1575 women, 80% were aged at least 60. The mean femoral neck T-score was - 1.6 and the mean lumbar spine T-score was -1.2. This cohort is currently being followed up.

CONCLUSIONS

QUALYOR will provide important information on the relationship between bone quality variables and fracture risk in women with moderately decreased BMD.

Bone disease in monoclonal gammopathy of undetermined significance: results from a screened population-based study

Previous studies have shown that individuals with monoclonal gammopathy of undetermined significance (MGUS) have an increased risk of fractures, although the underlying mechanisms remain unknown. Our aim was to analyze bone mineral density (BMD), bone volume, and risk of fractures among individuals with MGUS. We performed a screening using the Age, Gene/Environment Susceptibility (AGES)-Reykjavik Study cohort, consisting of 5764 elderly individuals, identifying 300 individuals with MGUS, and 275 with light-chain MGUS. Quantitative computerized tomography was performed in the lumbar spine and hip to evaluate BMD and bone geometry. Analysis of variance and the Tukey honest significance test were used to compare the groups. Hospital records were used to record fractures, with a mean follow-up of 6.9 years. Cox proportional hazard was used to compare fracture risk. No difference was found in BMD between subjects with MGUS and others in the spine (P = .34) or in total hip (P = .30). Individuals with MGUS had a significant increase in bone volume compared with others in the spine (P < .001) and total hip (P < .001). Overall, the risk of fractures was not significantly increased in individuals with MGUS (hazard ratio [HR], 1.19; 95% confidence interval [CI], 0.94-1.50). Men with MGUS had a significantly increased fracture risk, compared with other men (HR, 1.46; 95% CI, 1.03-2.08). Our results show that although individuals with MGUS do not have decreased BMD, bone volume is increased, and MGUS men have a 50% increased fracture risk. These results indicate that bone disease and fractures in MGUS differ from processes known from osteoporosis.

Bone mineral density measurements in vertebral specimens and phantoms using dual-layer spectral computed tomography

To assess whether phantomless calcium-hydroxyapatite (HA) specific bone mineral density (BMD) measurements with dual-layer spectral computed tomography are accurate in phantoms and vertebral specimens. Ex-vivo human vertebrae (n = 13) and a phantom containing different known HA concentrations were placed in a semi-anthropomorphic abdomen phantom with different extension rings simulating different degrees of obesity. Phantomless dual-layer spectral CT was performed at different tube current settings (500, 250, 125 and 50 mAs). HA-specific BMD was derived from spectral-based virtual monoenergetic images at 50 keV and 200 keV. Values were compared to the HA concentrations of the phantoms and conventional qCT measurements using a reference phantom, respectively. Above 125 mAs, errors for phantom measurements ranged between -1.3% to 4.8%, based on spectral information. In vertebral specimens, high correlations were found between BMD values assessed with spectral CT and conventional qCT (r ranging between 0.96 and 0.99; p < 0.001 for all) with different extension rings, and a high agreement was found in Bland Altman plots. Different degrees of obesity did not have a significant influence on measurements (P > 0.05 for all). These results suggest a high validity of HA-specific BMD measurements based on dual-layer spectral CT examinations in setups simulating different degrees of obesity without the need for a reference phantom, thus demonstrating their feasibility in clinical routine.

Critical issues and current challenges in osteoporosis and fracture prevention. An overview of unmet needs

Osteoporosis is a silent disease with increasing prevalence due to the global ageing population. Decreased bone strength and bone quality is the hallmark of osteoporosis which leads to an increased risk of fragility fractures in elderly. It has been estimated that approximately ~50% of women will suffer during their lifetime from an osteoporotic fracture. This must be considered as a major health concern, as it has previously been established that fragility fracture has been associated with decreased quality of life due to increased disability, more frequent hospital admission and most importantly, osteoporotic fractures have been related to an augmented mortality risk. Anti-osteoporotic drugs are available for improving bone quality. Although there is access to these therapeutic options, there remain multiple unmet needs in the field of osteoporosis and fracture care, for example, the primary prevention of osteoporosis in young individuals (to reach a high peak bone mass), the optimization of the use of imaging techniques [dual-energy X-ray absorptiometry (DXA), vertebral fracture assessment (VFA) and new techniques measuring bone quality], the use of nonmedical treatment options and surgical techniques of fracture healing. In this review, we will discuss topics that play a role in the occurrence and prevention of fractures, and we give an overview of and insight into the critical issues and challenges around osteoporosis and fracture prevention.

Comparison of non-invasive assessments of strength of the proximal femur

It is not clear which non-invasive method is most effective for predicting strength of the proximal femur in those at highest risk of fracture. The primary aim of this study was to compare the abilities of dual energy X-ray absorptiometry (DXA)-derived aBMD, quantitative computed tomography (QCT)-derived density and volume measures, and finite element analysis (FEA)-estimated strength to predict femoral failure load. We also evaluated the contribution of cortical and trabecular bone measurements to proximal femur strength. We obtained 76 human cadaveric proximal femurs (50 women and 26 men; age 74±8.8years), performed imaging with DXA and QCT, and mechanically tested the femurs to failure in a sideways fall configuration at a high loading rate. Linear regression analysis was used to construct the predictive model between imaging outcomes and experimentally-measured femoral strength for each method. To compare the performance of each method we used 3-fold cross validation repeated 10 times. The bone strength estimated by QCT-based FEA predicted femoral failure load (R²_adj=0.78, 95%CI 0.76-0.80; RMSE=896N, 95%CI 830-961) significantly better than femoral neck aBMD by DXA (R²_adj=0.69, 95%CI 0.66-0.72; RMSE=1011N, 95%CI 952-1069) and the QCT-based model (R²_adj=0.73, 95%CI 0.71-0.75; RMSE=932N, 95%CI 879-985). Both cortical and trabecular bone contribute to femoral strength, the contribution of cortical bone being higher in femurs with lower trabecular bone density. These findings have implications for optimizing clinical approaches to assess hip fracture risk. In addition, our findings provide new insights that will assist in interpretation of the effects of osteoporosis treatments that preferentially impact cortical versus trabecular bone.

Cost-effectiveness of Virtual Bone Strength Testing in Osteoporosis Screening Programs for Postmenopausal Women in the United States

Purpose To investigate whether assessment of bone strength with quantitative computed tomography (CT) in combination with dual-energy x-ray absorptiometry (DXA) is cost-effective as a screening tool for osteoporosis in postmenopausal women. Materials and Methods A state-transition microsimulation model of osteoporosis for postmenopausal women aged 55 years or older was developed with a lifetime horizon and U.S. societal perspective. All model inputs were derived from published literature. Three strategies were compared: no screening, DXA with T score-dependent rescreening intervals, and a combination of DXA and quantitative CT with different intervals (3, 5, and 10 years) at different screening initiation ages (55-65 years). Oral bisphosphonate therapy was started if DXA hip T scores were less than or equal to -2.5, 10-year risk for hip fracture was greater than 3% (World Health Organization Fracture Risk Assessment Tool score, or FRAX), 10-year risk for major osteoporotic fracture was greater than 20% (FRAX), quantitative CT femur bone strength was less than 3000 N, or occurrence of first fracture (eg, hip, vertebral body, wrist). Outcome measures were incremental cost-effectiveness ratios (ICERs) in 2015 U.S. dollars per quality-adjusted life year (QALY) gained and number of fragility fractures. Probabilistic sensitivity analysis was also performed. Results The most cost-effective strategy was combined DXA and quantitative CT screening starting at age 55 with quantitative CT screening every 5 years (ICER, $2000 per QALY). With this strategy, 12.8% of postmenopausal women sustained hip fractures in their remaining life (no screening, 18.7%; DXA screening, 15.8%). The corresponding percentages of vertebral fractures for DXA and quantitative CT with a 5-year interval, was 7.5%; no screening, 11.1%; DXA screening, 9%; for wrist fractures, 14%, 17.8%, and 16.4%, respectively; for other fractures, 22.6%, 30.8%, and 27.3%, respectively. In probabilistic sensitivity analysis, DXA and quantitative CT at age 55 years with quantitative CT screening every 5 years was the best strategy in more than 90% of all 1000 simulations (for thresholds of $50 000 per QALY and $100 000 per QALY). Conclusion Combined assessment of bone strength and bone mineral density is a cost-effective strategy for osteoporosis screening in postmenopausal women and has the potential to prevent a substantial number of fragility fractures. © RSNA, 2017 Online supplemental material is available for this article.

Critical issues and current challenges in osteoporosis and fracture prevention. An overview of unmet needs

Osteoporosis is a silent disease with increasing prevalence due to the global ageing population. Decreased bone strength and bone quality is the hallmark of osteoporosis which leads to an increased risk of fragility fractures in elderly. It has been estimated that approximately ~50% of women will suffer during their lifetime from an osteoporotic fracture. This must be considered as a major health concern, as it has previously been established that fragility fracture has been associated with decreased quality of life due to increased disability, more frequent hospital admission and most importantly, osteoporotic fractures have been related to an augmented mortality risk. Anti-osteoporotic drugs are available for improving bone quality. Although there is access to these therapeutic options, there remain multiple unmet needs in the field of osteoporosis and fracture care, for example, the primary prevention of osteoporosis in young individuals (to reach a high peak bone mass), the optimization of the use of imaging techniques [dual-energy X-ray absorptiometry (DXA), vertebral fracture assessment (VFA) and new techniques measuring bone quality], the use of nonmedical treatment options and surgical techniques of fracture healing. In this review, we will discuss topics that play a role in the occurrence and prevention of fractures, and we give an overview of and insight into the critical issues and challenges around osteoporosis and fracture prevention.

Incorporation of CT-based measurements of trunk anatomy into subject-specific musculoskeletal models of the spine influences vertebral loading predictions

We created subject-specific musculoskeletal models of the thoracolumbar spine by incorporating spine curvature and muscle morphology measurements from computed tomography (CT) scans to determine the degree to which vertebral compressive and shear loading estimates are sensitive to variations in trunk anatomy. We measured spine curvature and trunk muscle morphology using spine CT scans of 125 men, and then created four different thoracolumbar spine models for each person: (i) height and weight adjusted (Ht/Wt models); (ii) height, weight, and spine curvature adjusted (+C models); (iii) height, weight, and muscle morphology adjusted (+M models); and (iv) height, weight, spine curvature, and muscle morphology adjusted (+CM models). We determined vertebral compressive and shear loading at three regions of the spine (T8, T12, and L3) for four different activities. Vertebral compressive loads predicted by the subject-specific CT-based musculoskeletal models were between 54% lower to 45% higher from those estimated using musculoskeletal models adjusted only for subject height and weight. The impact of subject-specific information on vertebral loading estimates varied with the activity and spinal region. Vertebral loading estimates were more sensitive to incorporation of subject-specific spinal curvature than subject-specific muscle morphology. Our results indicate that individual variations in spine curvature and trunk muscle morphology can have a major impact on estimated vertebral compressive and shear loads, and thus should be accounted for when estimating subject-specific vertebral loading. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:2164-2173, 2017.

Phantomless calibration of CT scans for measurement of BMD and bone strength-Inter-operator reanalysis precision

Patient-specific phantomless calibration of computed tomography (CT) scans has the potential to simplify and expand the use of pre-existing clinical CT for quantitative bone densitometry and bone strength analysis for diagnostic and monitoring purposes. In this study, we quantified the inter-operator reanalysis precision errors for a novel implementation of patient-specific phantomless calibration, using air and either aortic blood or hip adipose tissue as internal calibrating reference materials, and sought to confirm the equivalence between phantomless and (traditional) phantom-based measurements. CT scans of the spine and hip for 25 women and 15 men (mean±SD age of 67±9years, range 41-86years), one scan per anatomic site per patient, were analyzed independently by two analysts using the VirtuOst software (O.N. Diagnostics, Berkeley, CA). The scans were acquired at 120kVp, with a slice thickness/increment of 3mm or less, on nine different CT scanner models across 24 different scanners. The main parameters assessed were areal bone mineral density (BMD) at the hip (total hip and femoral neck), trabecular volumetric BMD at the spine, and vertebral and femoral strength by finite element analysis; other volumetric BMD measures were also assessed. We found that the reanalysis precision errors for all phantomless measurements were ≤0.5%, which was as good as for phantom calibration. Regression analysis indicated equivalence of the phantom- versus phantomless-calibrated measurements (slope not different than unity, R²≥0.98). Of the main parameters assessed, non-significant paired mean differences (n=40) between the two measurements ranged from 0.6% for hip areal BMD to 1.1% for mid-vertebral trabecular BMD. These results indicate that phantom-equivalent measurements of both BMD and finite element-derived bone strength can be reliably obtained from CT scans using patient-specific phantomless calibration.

Fracture during oral bisphosphonate therapy is associated with deteriorated bone material strength index

BACKGROUND

Some patients experience fractures while receiving oral bisphosphonates (BPs) treatment. Clinical risk factors, advanced bone density loss, and microarchitecture deterioration have been associated with such fractures but bone tissue properties other than bone mineral density (BMD) have not been assessed.

METHODS

In a cross-sectional study of postmenopausal women on bisphosphonates for at least 4years with good adherence to treatment, 21 patients with incident fractures were compared with 18 treated patients without new fractures. Demographic and clinical variables, BMD, laboratory tests, and bone material strength index (BMSi) assessed by impact microindentation at the tibial diaphysis were recorded for all participants.

RESULTS

Clinical and laboratory results did not differ between patients taking BPs with incident fractures and those without new fractures. However, BMSi was significantly lower (mean±SD) in those who fractured (73.76±6.49) than in no-fracture patients (81.64±6.26; p=0.001). Lumbar spine (LS) BMD was also lower in fractured patients (p=0.03). Adjusted models including age, body mass index, years on BP treatment, and LS-BMD confirmed an increase in fracture risk per BMSi standard deviation decrease: adjusted OR 23.5 [95% CI 2.16 to 255.66], p=0.01. ROC analyses showed an area under the curve of 0.82 (95% CI 0.68 to 0.95) for BMSi, higher than that for BMD at any location, which ranged from 0.64 (95% CI 0.47 to 0.82) for femoral neck (FN) BMD to 0.71 (95% CI 0.55 to 0.87) for LS-BMD.

CONCLUSIONS

Patients who fracture while receiving BPs treatment have worse BMSi scores than BP-treated patients without fractures. The potential for BMSi to provide an additional osteoporosis treatment target should be explored.

Greater Gains in Spine and Hip Strength for Romosozumab Compared With Teriparatide in Postmenopausal Women With Low Bone Mass

Romosozumab is a monoclonal antibody that inhibits sclerostin and has been shown to reduce the risk of fractures within 12 months. In a phase II, randomized, placebo-controlled clinical trial of treatment-naïve postmenopausal women with low bone mass, romosozumab increased bone mineral density (BMD) at the hip and spine by the dual effect of increasing bone formation and decreasing bone resorption. In a substudy of that trial, which included placebo and teriparatide arms, here we investigated whether those observed increases in BMD also resulted in improvements in estimated strength, as assessed by finite element analysis. Participants received blinded romosozumab s.c. (210 mg monthly) or placebo, or open-label teriparatide (20 μg daily) for 12 months. CT scans, obtained at the lumbar spine (n = 82) and proximal femur (n = 46) at baseline and month 12, were analyzed with finite element software (VirtuOst, O.N. Diagnostics) to estimate strength for a simulated compression overload for the spine (L₁ vertebral body) and a sideways fall for the proximal femur, all blinded to treatment assignment. We found that, at month 12, vertebral strength increased more for romosozumab compared with both teriparatide (27.3% versus 18.5%; p = 0.005) and placebo (27.3% versus -3.9%; p < 0.0001); changes in femoral strength for romosozumab showed similar but smaller changes, increasing more with romosozumab versus teriparatide (3.6% versus -0.7%; p = 0.027), and trending higher versus placebo (3.6% versus -0.1%; p = 0.059). Compartmental analysis revealed that the bone-strengthening effects for romosozumab were associated with positive contributions from both the cortical and trabecular bone compartments at both the lumbar spine and hip. Taken together, these findings suggest that romosozumab may offer patients with osteoporosis a new bone-forming therapeutic option that increases both vertebral and femoral strength within 12 months. © 2017 American Society for Bone and Mineral Research.

Vertebral and femoral bone mineral density and bone strength in prostate cancer patients assessed in phantomless PET/CT examinations

PURPOSE

Bone fracture risk assessed ancillary to positron emission tomography with computed tomography co-registration (PET/CT) could provide substantial clinical value to oncology patients with elevated fracture risk without introducing additional radiation dose. The purpose of our study was to investigate the feasibility of obtaining valid measurements of bone mineral density (BMD) and finite element analysis-derived bone strength of the hip and spine using PET/CT examinations of prostate cancer patients by comparing against values obtained using routine multidetector-row computed tomography (MDCT) scans-as validated in previous studies-as a reference standard.

MATERIALS AND METHODS

Men with prostate cancer (n=82, 71.6±8.3 years) underwent Fluorine-18 NaF PET/CT and routine MDCT within three months. Femoral neck and total hip areal BMD, vertebral trabecular BMD and femur and vertebral strength based on finite element analysis were assessed in 63 paired PET/CT and MDCT examinations using phantomless calibration and Biomechanical-CT analysis. Men with osteoporosis or fragile bone strength identified at either the hip or spine (vertebral trabecular BMD ≤80mg/cm³, femoral neck or total hip T-score ≤-2.5, vertebral strength ≤6500N and femoral strength ≤3500N, respectively) were considered to be at high risk of fracture. PET/CT- versus MDCT-based BMD and strength measurements were compared using paired t-tests, linear regression and by generating Bland-Altman plots. Agreement in fracture-risk classification was assessed in a contingency table.

RESULTS

All measurements from PET/CT versus MDCT were strongly correlated (R²=0.93-0.97; P<0.0001 for all). Mean differences for total hip areal BMD (0.001g/cm², 1.1%), femoral strength (-60N, 1.3%), vertebral trabecular BMD (2mg/cm³, 2.6%) and vertebral strength (150N; 1.7%) measurements were not statistically significant (P>0.05 for all), whereas the mean difference in femoral neck areal BMD measurements was small but significant (-0.018g/cm²; -2.5%; P=0.007). The agreement between PET/CT and MDCT for fracture-risk classification was 97% (0.89 kappa for repeatability).

CONCLUSION

Ancillary analyses of BMD, bone strength, and fracture risk agreed well between PET/CT and MDCT, suggesting that PET/CT can be used opportunistically to comprehensively assess bone integrity. In subjects with high fracture risk such as cancer patients this may serve as an additional clinical tool to guide therapy planning and prevention of fractures.

Quantitative Computed Tomography-Current Status and New Developments

This review focuses on new developments and current controversies in the field of quantitative computed tomography. Recent positions of the International Society for Clinical Densitometry acknowledged the clinical value of quantitative computed tomography of the spine and the hip using clinical whole-body computed tomography (CT) scanners. Opportunistic screening summarizes a number of new approaches describing the dual use of clinical CT scans. For example, CT scans may have been taken for tumor diagnosis but may also be used for the prediction of high or low fracture risks as an additional benefit for the patient. The assessment of the cortical parameters is another topic of current research. In CT images of the spine and the hip, a number of techniques have been developed to determine the thickness, mass, and bone density of the cortex. In higher-spatial resolution peripheral CT images of the radius and tibia obtained from special purpose scanners, 1 focus is the measurement of cortical porosity. Two different approaches, one based on the direct segmentation of the pores and one based on cortical density, will be reviewed.

Predicting Future Hip Fractures on Routine Abdominal CT Using Opportunistic Osteoporosis Screening Measures: A Matched Case-Control Study

OBJECTIVE

Hip fracture is a major consequence of low bone mineral density, which is treatable but underdiagnosed. The purpose of this case-control study is to determine whether lumbar vertebral trabecular attenuation, vertebral compression fractures, and femoral neck T scores readily derived from abdominopelvic CT scans obtained for various indications are associated with future hip fragility fracture.

MATERIALS AND METHODS

A cohort of 204 patients with hip fracture (130 women and 74 men; mean age, 74.3 years) who had undergone abdominopelvic CT before fracture occurred (mean interval, 24.8 months) was compared with an age- and sex-matched control cohort without hip fracture. L1 trabecular attenuation, vertebral compression fractures of grades 2 and 3, and femoral neck T scores derived from asynchronous quantitative CT were recorded. The presence of one or more clinical risk factor for fracture was also recorded. Multivariate logistic regression models were used to determine the association of each measurement with the occurrence of hip fracture.

RESULTS

The mean L1 trabecular attenuation value, the presence of one or more vertebral compression fracture, and CT-derived femoral neck T scores were all significantly different in patients with hip fracture versus control subjects (p < 0.01). Logistic regression models showed a significant association of all measurements with hip fracture outcome after adjustments were made for age, sex, and the presence of one or more clinical risk factor. L1 trabecular attenuation and CT-derived femoral neck T scores showed moderate accuracy in differentiating case and control patients (AUC, 0.70 and 0.78, respectively).

CONCLUSION

L1 trabecular attenuation, CT-derived femoral neck T scores, and the presence of at least one vertebral compression fracture on CT are all associated with future hip fragility fracture in adults undergoing routine abdominopelvic CT for a variety of conditions.

Low-trauma fractures without osteoporosis

In clinical practice, areal bone mineral density (aBMD) is usually measured using dual-energy X-ray absorptiometry (DXA) to assess bone status in patients with or without osteoporotic fracture. As BMD has a Gaussian distribution, it is difficult to define a cutoff for osteoporosis diagnosis. Based on epidemiological considerations, WHO defined a DXA-based osteoporosis diagnosis with a T-score <-2.5. However, the majority of individuals who have low-trauma fractures do not have osteoporosis with DXA (i.e., T-score <-2.5), and some of them have no decreased BMD at all. Some medical conditions (spondyloarthropathies, chronic kidney disease and mineral bone disorder, diabetes, obesity) or drugs (glucocorticoids, aromatase inhibitors) are more prone to cause fractures with subnormal BMD. In the situation of fragility fractures with subnormal or normal BMD, clinicians face a difficulty as almost all the pharmacologic treatments have proved their efficacy in patients with low BMD. However, some data are available in post hoc analyses in patients with T score >-2. Overall, in patients with a previous fragility fracture (especially vertebra or hip), treatments appear to be effective. Thus, the authors recommend treating some patients with a major fragility fracture even if areal BMD T score is above -2.5.

Virtual stress testing of fracture stability in soldiers with severely comminuted tibial fractures

Virtual stress testing (VST) provides a non-invasive estimate of the strength of a healing bone through a biomechanical analysis of a patient's computed tomography (CT) scan. We asked whether VST could improve management of patients who had a tibia fracture treated with external fixation. In a retrospective case-control study of 65 soldier-patients who had tibia fractures treated with an external fixator, we performed VST utilizing CT scans acquired prior to fixator removal. The strength of the healing bone and the amount of tissue damage after application of an overload were computed for various virtual loading cases. Logistic regression identified computed outcomes with the strongest association to clinical events related to nonunion within 2 months after fixator removal. Clinical events (n = 9) were associated with a low tibial strength for compression loading (p < 0.05, AUC = 0.74) or a low proportion of failed cortical bone tissue for torsional loading (p < 0.005, AUC = 0.84). Using post-hoc thresholds of a compressive strength of four times body-weight and a proportional of failed cortical bone tissue of 5%, the test identified all nine patients who failed clinically (100% sensitivity; 40.9% positive predictive value) and over three fourths of those (43 of 56) who progressed to successful healing (76.8% specificity; 100% negative predictive value). In this study, VST identified all patients who progressed to full, uneventful union after fixator removal; thus, we conclude that this new test has the potential to provide a quantitative, objective means of identifying tibia-fracture patients who can safely resume weight bearing. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:805-811, 2017.

Advances in imaging approaches to fracture risk evaluation

Fragility fractures are a growing problem worldwide, and current methods for diagnosing osteoporosis do not always identify individuals who require treatment to prevent a fracture and may misidentify those not a risk. Traditionally, fracture risk is assessed using dual-energy X-ray absorptiometry, which provides measurements of areal bone mineral density at sites prone to fracture. Recent advances in imaging show promise in adding new information that could improve the prediction of fracture risk in the clinic. As reviewed herein, advances in quantitative computed tomography (QCT) predict hip and vertebral body strength; high-resolution HR-peripheral QCT (HR-pQCT) and micromagnetic resonance imaging assess the microarchitecture of trabecular bone; quantitative ultrasound measures the modulus or tissue stiffness of cortical bone; and quantitative ultrashort echo-time MRI methods quantify the concentrations of bound water and pore water in cortical bone, which reflect a variety of mechanical properties of bone. Each of these technologies provides unique characteristics of bone and may improve fracture risk diagnoses and reduce prevalence of fractures by helping to guide treatment decisions.

Clinical experience with microindentation in vivo in humans

Densitometry and imaging techniques are currently used in clinical settings to measure bone quantity and spatial structure. Recently, Reference Point Indentation has opened the possibility of directly assessing the mechanical characteristics of cortical bone in living individuals, adding a new dimension to the assessment of bone strength. Impact microindentation was specifically developed for clinical studies and has been tested in several populations where there are discrepancies between bone density and fracture propensity, such as type 2 diabetes, atypical femoral fracture, stress fractures, glucocorticoid treatment, patients with osteopenia and fragility fractures, and individuals infected with HIV, among others. Microindentation will complement, not replace, existing bone analysis methods, particularly where bone mineral density does not fully explain fracture propensity. The available evidence provides solid proof of concept; future studies will fully define the role of microindentation for the assessment of bone health both in clinics and in research.

Risk of vertebral compression fractures in multiple myeloma patients: A finite-element study

The purpose of this study was to develop and validate a finite element (FE) model to predict vertebral bone strength in vitro using multidetector computed tomography (MDCT) images in multiple myeloma (MM) patients, to serve as a complementing tool to assess fracture risk. In addition, it also aims to differentiate MM patients with and without vertebral compression fractures (VCFs) by performing FE analysis on vertebra segments (T1-L5) obtained from in vivo routine MDCT imaging scans. MDCT-based FE models were developed from the in vitro vertebrae samples and were then applied to the in vivo vertebrae segments of MM patients (n = 4) after validation. Predicted fracture load using FE models correlated significantly with experimentally measured failure load (r = 0.85, P < 0.001). Interestingly, an erratic behavior was observed in patients with fractures (n = 2) and a more gradual change in FE-predicted strength values in patients without fractures (n = 2). Severe geometric deformations were also observed in models that have already attained fractures. Since BMD is not a reliable parameter for fracture risk prediction in MM subjects, it is necessary to use advanced tools such as FE analysis to predict individual fracture risk. If peaks are observed between adjacent segments in an MM patient, it can be safe to conclude that the spine is experiencing regions of structural instability. Such an FE visualization may have therapeutic consequences to prevent MM associated vertebral fractures.

Focal osteoporosis defects play a key role in hip fracture

BACKGROUND

Hip fractures are mainly caused by accidental falls and trips, which magnify forces in well-defined areas of the proximal femur. Unfortunately, the same areas are at risk of rapid bone loss with ageing, since they are relatively stress-shielded during walking and sitting. Focal osteoporosis in those areas may contribute to fracture, and targeted 3D measurements might enhance hip fracture prediction. In the FEMCO case-control clinical study, Cortical Bone Mapping (CBM) was applied to clinical computed tomography (CT) scans to define 3D cortical and trabecular bone defects in patients with acute hip fracture compared to controls. Direct measurements of trabecular bone volume were then made in biopsies of target regions removed at operation.

METHODS

The sample consisted of CT scans from 313 female and 40 male volunteers (158 with proximal femoral fracture, 145 age-matched controls and 50 fallers without hip fracture). Detailed Cortical Bone Maps (c.5580 measurement points on the unfractured hip) were created before registering each hip to an average femur shape to facilitate statistical parametric mapping (SPM). Areas where cortical and trabecular bone differed from controls were visualised in 3D for location, magnitude and statistical significance. Measures from the novel regions created by the SPM process were then tested for their ability to classify fracture versus control by comparison with traditional CT measures of areal Bone Mineral Density (aBMD). In women we used the surgical classification of fracture location ('femoral neck' or 'trochanteric') to discover whether focal osteoporosis was specific to fracture type. To explore whether the focal areas were osteoporotic by histological criteria, we used micro CT to measure trabecular bone parameters in targeted biopsies taken from the femoral heads of 14 cases.

RESULTS

Hip fracture patients had distinct patterns of focal osteoporosis that determined fracture type, and CBM measures classified fracture type better than aBMD parameters. CBM measures however improved only minimally on aBMD for predicting any hip fracture and depended on the inclusion of trabecular bone measures alongside cortical regions. Focal osteoporosis was confirmed on biopsy as reduced sub-cortical trabecular bone volume.

CONCLUSION

Using 3D imaging methods and targeted bone biopsy, we discovered focal osteoporosis affecting trabecular and cortical bone of the proximal femur, among men and women with hip fracture.

Effects of dose reduction on bone strength prediction using finite element analysis

This study aimed to evaluate the effect of dose reduction, by means of tube exposure reduction, on bone strength prediction from finite-element (FE) analysis. Fresh thoracic mid-vertebrae specimens (n = 11) were imaged, using multi-detector computed tomography (MDCT), at different intensities of X-ray tube exposures (80, 150, 220 and 500 mAs). Bone mineral density (BMD) was estimated from the mid-slice of each specimen from MDCT images. Differences in image quality and geometry of each specimen were measured. FE analysis was performed on all specimens to predict fracture load. Paired t-tests were used to compare the results obtained, using the highest CT dose (500 mAs) as reference. Dose reduction had no significant impact on FE-predicted fracture loads, with significant correlations obtained with reference to 500 mAs, for 80 mAs (R²  = 0.997, p < 0.001), 150 mAs (R² = 0.998, p < 0.001) and 220 mAs (R² = 0.987, p < 0.001). There were no significant differences in volume quantification between the different doses examined. CT imaging radiation dose could be reduced substantially to 64% with no impact on strength estimates obtained from FE analysis. Reduced CT dose will enable early diagnosis and advanced monitoring of osteoporosis and associated fracture risk.