A novel registration-based methodology for prediction of trabecular bone fabric from clinical QCT: A comprehensive analysis

2D-3D reconstruction of the proximal femur from DXA scans: Evaluation of the 3D-Shaper software

Introduction: Osteoporosis is currently diagnosed based on areal bone mineral density (aBMD) computed from 2D DXA scans. However, aBMD is a limited surrogate for femoral strength since it does not account for 3D bone geometry and density distribution. QCT scans combined with finite element (FE) analysis can deliver improved femoral strength predictions. However, non-negligible radiation dose and high costs prevent a systematic usage of this technique for screening purposes. As an alternative, the 3D-Shaper software (3D-Shaper Medical, Spain) reconstructs the 3D shape and density distribution of the femur from 2D DXA scans. This approach could deliver a more accurate estimation of femoral strength than aBMD by using FE analysis on the reconstructed 3D DXA. Methods: Here we present the first independent evaluation of the software, using a dataset of 77 ex vivo femora. We extend a prior evaluation by including the density distribution differences, the spatial correlation of density values and an FE analysis. Yet, cortical thickness is left out of this evaluation, since the cortex is not resolved in our FE models. Results: We found an average surface distance of 1.16 mm between 3D DXA and QCT images, which shows a good reconstruction of the bone geometry. Although BMD values obtained from 3D DXA and QCT correlated well (r ² = 0.92), the 3D DXA BMD were systematically lower. The average BMD difference amounted to 64 mg/cm³, more than one-third of the 3D DXA BMD. Furthermore, the low correlation (r ² = 0.48) between density values of both images indicates a limited reconstruction of the 3D density distribution. FE results were in good agreement between QCT and 3D DXA images, with a high coefficient of determination (r ² = 0.88). However, this correlation was not statistically different from a direct prediction by aBMD. Moreover, we found differences in the fracture patterns between the two image types. QCT-based FE analysis resulted mostly in femoral neck fractures and 3D DXA-based FE in subcapital or pertrochanteric fractures. Discussion: In conclusion, 3D-Shaper generates an altered BMD distribution compared to QCT but, after careful density calibration, shows an interesting potential for deriving a standardized femoral strength from a DXA scan.

A combined active shape and mean appearance model for the reconstruction of segmental bone loss

This study investigates the novel combination of an active shape and mean appearance model to estimate missing bone geometry and density distribution from sparse inputs simulating segmental bone loss of the femoral diaphysis. An active shape Gaussian Process Morphable model was trained on healthy right femurs of South African males to model shape. The density distribution was approximated based on the mean appearance of computed tomography images from the training set. Estimations of diaphyseal resections were obtained by probabilistic fitting of the active shape model to sparse inputs consisting of proximal and distal femoral data on computed tomography images. The resulting shape estimates of the diaphyseal resections were then used to map the mean appearance model to the patients' missing bone geometry, constructing density estimations. In this way, resected bone surfaces were estimated with an average error of 2.24 (0.5) mm. Density distributions were approximated within 87 (0.7) % of the intensity of the original target images before the simulated segmental bone loss. These results fall within the acceptable tolerances required for surgical planning and reconstruction of long bone defects.

The effect of cement augmentation on pedicle screw fixation under various load cases : results from a combined experimental, micro-CT, and micro-finite element analysis

Aims

Anchorage of pedicle screw rod instrumentation in the elderly spine with poor bone quality remains challenging. Our study aims to evaluate how the screw bone anchorage is affected by screw design, bone quality, loading conditions, and cementing techniques.

Methods

Micro-finite element (µFE) models were created from micro-CT (μCT) scans of vertebrae implanted with two types of pedicle screws (L: Ennovate and R: S⁴). Simulations were conducted for a 10 mm radius region of interest (ROI) around each screw and for a full vertebra (FV) where different cementing scenarios were simulated around the screw tips. Stiffness was calculated in pull-out and anterior bending loads.

Results

Experimental pull-out strengths were excellently correlated to the µFE pull-out stiffness of the ROI (R² > 0.87) and FV (R² > 0.84) models. No significant difference due to screw design was observed. Cement augmentation increased pull-out stiffness by up to 94% and 48% for L and R screws, respectively, but only increased bending stiffness by up to 6.9% and 1.5%, respectively. Cementing involving only one screw tip resulted in lower stiffness increases in all tested screw designs and loading cases. The stiffening effect of cement augmentation on pull-out and bending stiffness was strongly and negatively correlated to local bone density around the screw (correlation coefficient (R) = -0.95).

Conclusion

This combined experimental, µCT and µFE study showed that regional analyses may be sufficient to predict fixation strength in pull-out and that full analyses could show that cement augmentation around pedicle screws increased fixation stiffness in both pull-out and bending, especially for low-density bone.

Cite this article: Bone Joint Res 2021;10(12):797–806.

Mechanical metric for skeletal biomechanics derived from spectral analysis of stiffness matrix

A new metric for the quantitative and qualitative evaluation of bone stiffness is introduced. It is based on the spectral decomposition of stiffness matrix computed with finite element method. The here proposed metric is defined as an amplitude rescaled eigenvalues of stiffness matrix. The metric contains unique information on the principal stiffness of bone and reflects both bone shape and material properties. The metric was compared with anthropometrical measures and was tested for sex sensitivity on pelvis bone. Further, the smallest stiffness of pelvis was computed under a certain loading condition and analyzed with respect to sex and direction. The metric complements anthropometrical measures and provides a unique information about the smallest bone stiffness independent from the loading configuration and can be easily computed by state-of-the-art subject specified finite element algorithms.

Finite element models can reproduce the effect of nucleotomy on the multi-axial compliance of human intervertebral discs

Finite element (FE) models can unravel the link between intervertebral disc (IVD) degeneration and its mechanical behaviour. Nucleotomy may provide the data required for model verification. Three human IVDs were scanned with MRI and tested in multiple loading scenarios, prior and post nucleotomy. The resulting data was used to generate, calibrate, and verify the FE models. Nucleotomy increased the experimental range of motion by 26%, a result reproduced by the FE simulation within a 5% error. This work demonstrates the ability of FE models to reproduce the mechanical compliance of human IVDs prior and post nucleotomy.

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.