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

The predictive ability of a QCT-FE model of the proximal femoral stiffness under multiple load cases is strongly influenced by experimental uncertainties

Despite significant improvements in terms of the predictive ability of Quantitative Computed Tomography based Finite Element (QCT-FE) models in estimating femoral strength (fracture load and stiffness), no substantial clinical adoption of this method has taken place to date. Narrowing the wide variability of FE results by standardizing the methodology and validation protocols, as well as reducing the uncertainties in the FEA process have been proposed as routes towards improved reliability. The aim of this study was to: First, validate a QCT-FE model of proximal femoral stiffness in multiple stance load cases, and second, using a parametric approach, determine the influence of select experimental and modeling parameters on the predictive ability of our model. Ten fresh frozen human femoral samples were tested in neutral stance, 15° adducted and 15° abducted load cases. Voxel-based linear-elastic QCT-FE models of the samples were generated to predict the models' stiffness values in all load cases. The base FE models were validated against the experimental results using linear regression. Thirty six deviated models were created using the minimum and maximum values of experiment-based "plausible range" for 18 parameters in 4 categories of embedding, loading, material, and segmentation. The predictive ability of the models were compared in terms of the coefficient of determination (R²) of the linear regression between the measured and predicted stiffness values in all load cases. Our model was capable of capturing 90% of the variation in the experimental stiffness of the samples in neutral stance position (R² = 0.9, concordance correlation coefficient (CCC) = 0.93, percent root mean squared error (RMSE%) = 8.4%, slope and intercept not significantly different from unity and zero, respectively). Embedding and loading categories strongly affected the predictive ability of the models with an average percent difference in R² of 4.36% ± 2.77 and 2.96% ± 1.69 for the stance-neutral load case, respectively. The performance of the models were significantly different in adducted and abducted load cases with their R² dropping to 71% and 70%, respectively. Similarly, off-axes load cases were affected by the parameters differently compared to the neutral load case, with the loading parameter category imposing more than 10% difference on their R², larger than all other categories. We also showed that automatically selecting the best performing plausible value for each parameter and each sample would result in a perfectly linear correlation (R²> 0.99) between the "tuned" model's predicted stiffness and experimental results. Based on our results, high sensitivity of the model performance to experimental parameters requires extra diligence in modeling the embedding geometry and the loading angles since these sources of uncertainty could dwarf the effects of material modeling and image processing parameters. The results of this study could help in improving the robustness of the QCT-FE models of proximal femur by limiting the uncertainties in the experimental and modeling steps.

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.

Fracture Risk Evaluation of Bone Metastases: A Burning Issue

Major progress has been achieved to treat cancer patients and survival has improved considerably, even for stage-IV bone metastatic patients. Locomotive health has become a crucial issue for patient autonomy and quality of life. The centerpiece of the reflection lies in the fracture risk evaluation of bone metastasis to guide physician decision regarding physical activity, antiresorptive agent prescription, and local intervention by radiotherapy, surgery, and interventional radiology. A key mandatory step, since bone metastases may be asymptomatic and disseminated throughout the skeleton, is to identify the bone metastasis location by cartography, especially within weight-bearing bones. For every location, the fracture risk evaluation relies on qualitative approaches using imagery and scores such as Mirels and spinal instability neoplastic score (SINS). This approach, however, has important limitations and there is a need to develop new tools for bone metastatic and myeloma fracture risk evaluation. Personalized numerical simulation qCT-based imaging constitutes one of these emerging tools to assess bone tumoral strength and estimate the femoral and vertebral fracture risk. The next generation of numerical simulation and artificial intelligence will take into account multiple loadings to integrate movement and obtain conditions even closer to real-life, in order to guide patient rehabilitation and activity within a personalized-medicine approach.

Fracture Risk of Long Bone Metastases: A Review of Current and New Decision-Making Tools for Prophylactic Surgery

Simple Summary

Long bone metastases are frequently a pivotal point in the oncological history of patients. Weakening of the bone results in pathologic fractures that not only compromise patient function but also their survival. Therefore, the main issue for tumor boards remains timely assessment of the risk of fracture, as this is a key consideration in providing preventive surgery while also avoiding overtreatment. As the Mirels scoring system takes into account both the radiological and the clinical criteria, it has been used worldwide since the 1990s. However, due to increasing concern regarding the lack of accuracy, new thresholds have been defined for the identification of impending fractures that require prophylactic surgery, on the basis of axial cortical involvement and biomechanical models involving quantitative computed tomography. The aim of this review is to establish a state-of-the-art of the risk assessment of long bone metastases fractures, from simple radiologic scores to more complex multidimensional bone models, in order to define new decision-making tools.

Abstract

Long bone pathological fractures very much reflect bone metastases morbidity in many types of cancer. Bearing in mind that they not only compromise patient function but also survival, identifying impending fractures before the actual event is one of the main concerns for tumor boards. Indeed, timely prophylactic surgery has been demonstrated to increase patient quality of life as well as survival. However, early surgery for long bone metastases remains controversial as the current fracture risk assessment tools lack accuracy. This review first focuses on the gold standard Mirels rating system. It then explores other unique imaging thresholds such as axial or circumferential cortical involvement and the merits of nuclear imaging tools. To overcome the lack of specificity, other fracture prediction strategies have focused on biomechanical models based on quantitative computed tomography (CT): computed tomography rigidity analysis (CT-RA) and finite element analysis (CT-FEA). Despite their higher specificities in impending fracture assessment, their limited availability, along with a need for standardization, have limited their use in everyday practice. Currently, the prediction of long bone pathologic fractures is a multifactorial process. In this regard, machine learning could potentially be of value by taking into account clinical survival prediction as well as clinical and improved CT-RA/FEA data.

Femoral strength can be predicted from 2D projections using a 3D statistical deformation and texture model with finite element analysis

Ultimate force of the proximal human femur can be predicted using Finite Element Analysis (FEA), but the models rely on 3D computed tomography images. Landmark-based statistical appearance models (SAM) and B-Spline transformation-based statistical deformation models (SDM) have been used to estimate 3D images from 2D projections, which facilitates model generation and reduces the radiation dose. However, there is no literature on the accuracy of SDM-based FEA models of bones with respect to experimental results. In this study, a methodology for an enhanced SDM with textural information is presented. The statistical deformation and texture models (SDTMs) are based on a set of 37 quantitative CT (QCT) images. They were used to estimate 3D images from two or one projections of the set in a leave-one-out setup. These estimations where then used to create FEA models. The ultimate force predicted by FEA models estimated from two or one projection using the SDTMs were compared to the experimental ultimate force from a previous study on the same femora and to the results of standard QCT-based FEA models. High correlations between predictions and experimental measurements were found for FEA models reconstructed from 2D projections with R²=0.835 when based on two projections and R²=0.724 when using one projection. The correlations were comparable to those reached with standard QCT-based FE-models with the experimental results (R²=0.795). This study shows the high potential of SDTM-based 3D image reconstruction and FEA modelling from 2D projections to predict femoral ultimate force.

Post-operative fracture risk assessment following tumor curettage in the distal femur: a hybrid in vitro and in silico biomechanical approach

The distal femur is the predominant site for benign bone tumours and a common site for fracture following tumour removal or cementation. However, the lack of conclusive assessment criterion for post-operative fracture risk and appropriate devices for cement augmentation are serious concerns. Hence, a validated biomechanical tool was developed to assess bone strength, depending on the size and location of artificially created tumorous defects in the distal femora. The mechanics of the bone–cement interface was investigated to determine the main causes of reconstruction failure. Based on quantitative-CT images, non-linear and heterogeneous finite element (FE) models of human cadaveric distal femora with simulated tumourous defects were created and validated using in vitro mechanical tests from 14 cadaveric samples. Statistical analyses demonstrated a strong linear relationship (R² = 0.95, slope = 1.12) with no significant difference between bone strengths predicted by in silico analyses and in vitro tests (P = 0.174). FE analyses showed little reduction in bone strength until the defect was 35% or more of epiphyseal volume, and reduction in bone strength was less pronounced for laterally located defects than medial side defects. Moreover, the proximal end of the cortical window and the most interior wall of the bone–cement interface were the most vulnerable sites for reconstruction failure.

Image-based finite-element modeling of the human femur

Fracture is considered a critical clinical endpoint in skeletal pathologies including osteoporosis and bone metastases. However, current clinical guidelines are limited with respect to identifying cases at high risk of fracture, as they do not account for many mechanical determinants that contribute to bone fracture. Improving fracture risk assessment is an important area of research with clear clinical relevance. Patient-specific numerical musculoskeletal models generated from diagnostic images are widely used in biomechanics research and may provide the foundation for clinical tools used to quantify fracture risk. However, prior to clinical translation, in vitro validation of predictions generated from such numerical models is necessary. Despite adopting radically different models, in vitro validation of image-based finite element (FE) models of the proximal femur (predicting strains and failure loads) have shown very similar, encouraging levels of accuracy. The accuracy of such in vitro models has motivated their application to clinical studies of osteoporotic and metastatic fractures. Such models have demonstrated promising but heterogeneous results, which may be explained by the lack of a uniform strategy with respect to FE modeling of the human femur. This review aims to critically discuss the state of the art of image-based femoral FE modeling strategies, highlighting principal features and differences among current approaches. Quantitative results are also reported with respect to the level of accuracy achieved from in vitro evaluations and clinical applications and are used to motivate the adoption of a standardized approach/workflow for image-based FE modeling of the femur.

Empirical relationships between bone density and ultimate strength: A literature review

INTRODUCTION

Ultimate strength-density relationships for bone have been reported with widely varying results. Reliable bone strength predictions are crucial for many applications that aim to assess bone failure. Bone density and bone morphology have been proposed to explain most of the variance in measured bone strength. If this holds true, it could lead to the derivation of a single ultimate strength-density-morphology relationship for all anatomical sites.

METHODS

All relevant literature was reviewed. Ultimate strength-density relationships derived from mechanical testing of human bone tissue were included. The reported relationships were translated to ultimate strength-apparent density relationships and normalized with respect to strain rate. Results were grouped based on bone tissue type (cancellous or cortical), anatomical site, and loading mode (tension vs. compression). When possible, the relationships were compared to existing ultimate strength-density-morphology relationships.

RESULTS

Relationships that considered bone density and morphology covered the full spectrum of eight-fold inter-study difference in reported compressive ultimate strength-density relationships for trabecular bone. This was true for studies that tested specimens in different loading direction and tissue from different anatomical sites. Sparse data was found for ultimate strength-density relationships in tension and for cortical bone properties transverse to the main loading axis of the bone.

CONCLUSIONS

Ultimate strength-density-morphology relationships could explain measured strength across anatomical sites and loading directions. We recommend testing of bone specimens in other directions than along the main trabecular alignment and to include bone morphology in studies that investigate bone material properties. The lack of tensile strength data did not allow for drawing conclusions on ultimate strength-density-morphology relationships. Further studies are needed. Ideally, these studies would investigate both tensile and compressive strength-density relationships, including morphology, to close this gap and lead to more accurate evaluation of bone failure.

Patient-Specific Bone Multiscale Modelling, Fracture Simulation and Risk Analysis-A Survey

This paper provides a starting point for researchers and practitioners from biology, medicine, physics and engineering who can benefit from an up-to-date literature survey on patient-specific bone fracture modelling, simulation and risk analysis. This survey hints at a framework for devising realistic patient-specific bone fracture simulations. This paper has 18 sections: Section 1 presents the main interested parties; Section 2 explains the organzation of the text; Section 3 motivates further work on patient-specific bone fracture simulation; Section 4 motivates this survey; Section 5 concerns the collection of bibliographical references; Section 6 motivates the physico-mathematical approach to bone fracture; Section 7 presents the modelling of bone as a continuum; Section 8 categorizes the surveyed literature into a continuum mechanics framework; Section 9 concerns the computational modelling of bone geometry; Section 10 concerns the estimation of bone mechanical properties; Section 11 concerns the selection of boundary conditions representative of bone trauma; Section 12 concerns bone fracture simulation; Section 13 presents the multiscale structure of bone; Section 14 concerns the multiscale mathematical modelling of bone; Section 15 concerns the experimental validation of bone fracture simulations; Section 16 concerns bone fracture risk assessment. Lastly, glossaries for symbols, acronyms, and physico-mathematical terms are provided.

Radiofrequency echographic multi-spectrometry for the in-vivo assessment of bone strength: state of the art-outcomes of an expert consensus meeting organized by the European Society for Clinical and Economic Aspects of Osteoporosis, Osteoarthritis and Musculoskeletal Diseases (ESCEO)

PURPOSE

The purpose of this paper was to review the available approaches for bone strength assessment, osteoporosis diagnosis and fracture risk prediction, and to provide insights into radiofrequency echographic multi spectrometry (REMS), a non-ionizing axial skeleton technique.

METHODS

A working group convened by the European Society for Clinical and Economic Aspects of Osteoporosis and Osteoarthritis met to review the current image-based methods for bone strength assessment and fracture risk estimation, and to discuss the clinical perspectives of REMS.

RESULTS

Areal bone mineral density (BMD) measured by dual-energy X-ray absorptiometry (DXA) is the consolidated indicator for osteoporosis diagnosis and fracture risk assessment. A more reliable fracture risk estimation would actually require an improved assessment of bone strength, integrating also bone quality information. Several different approaches have been proposed, including additional DXA-based parameters, quantitative computed tomography, and quantitative ultrasound. Although each of them showed a somewhat improved clinical performance, none satisfied all the requirements for a widespread routine employment, which was typically hindered by unclear clinical usefulness, radiation doses, limited accessibility, or inapplicability to spine and hip, therefore leaving several clinical needs still unmet. REMS is a clinically available technology for osteoporosis diagnosis and fracture risk assessment through the estimation of BMD on the axial skeleton reference sites. Its automatic processing of unfiltered ultrasound signals provides accurate BMD values in view of fracture risk assessment.

CONCLUSIONS

New approaches for improved bone strength and fracture risk estimations are needed for a better management of osteoporotic patients. In this context, REMS represents a valuable approach for osteoporosis diagnosis and fracture risk prediction.

A multiscale model to predict current absolute risk of femoral fracture in a postmenopausal population

Osteoporotic hip fractures are a major healthcare problem. Fall severity and bone strength are important risk factors of hip fracture. This study aims to obtain a mechanistic explanation for fracture risk in dependence of these risk factors. A novel modelling approach is developed that combines models at different scales to overcome the challenge of a large space–time domain of interest and considers the variability of impact forces between potential falls in a subject. The multiscale model and its component models are verified with respect to numerical approximations made therein, the propagation of measurement uncertainties of model inputs is quantified, and model predictions are validated against experimental and clinical data. The main results are model predicted absolute risk of current fracture (ARF0) that ranged from 1.93 to 81.6% (median 36.1%) for subjects in a retrospective cohort of 98 postmenopausal British women (49 fracture cases and 49 controls); ARF0 was computed up to a precision of 1.92 percentage points (pp) due to numerical approximations made in the model; ARF0 possessed an uncertainty of 4.00 pp due to uncertainties in measuring model inputs; ARF0 classified observed fracture status in the above cohort with AUC = 0.852 (95% CI 0.753–0.918), 77.6% specificity (95% CI 63.4–86.5%) and 81.6% sensitivity (95% CI 68.3–91.1%). These results demonstrate that ARF0 can be computed using the model with sufficient precision to distinguish between subjects and that the novel mechanism of fracture risk determination based on fall dynamics, hip impact and bone strength can be considered validated.