Sensitivity of patient-specific vertebral finite element model from low dose imaging to material properties and loading conditions

Effect of postural alignment alteration with age on vertebral strength

EOS biplane radiographs of 117 subjects between 20 and 83 years were analyzed to compute the upper body lever arm over the L1 vertebra and its impact on vertebral strength. Postural sagittal alignment alteration was observed with age and resulted in a greater lever arm causing vertebral strength to decrease.

PURPOSE

The purpose of this study was to analyze the impact of postural alignment changes with age on vertebral strength using finite element analysis and barycentremetry.

METHODS

A total of 117 subjects from 20 to 83 years were divided in three age groups: young (20 to 40 years, 62 subjects), intermediate (40 to 60 years, 26 subjects), and elderly (60 years and over, 29 subjects). EOS biplane radiographs were acquired, allowing 3D reconstruction of the spine and body envelope as well as spinal, pelvic, and sagittal alignment parameter measurements. A barycentremetry method allowed the estimation of the mass and center of mass (CoM) position of the upper body above L1, relatively to the center of the L1 vertebra (lever arm). To investigate the effect of this lever arm, vertebral strength of a generic finite element model (with constant geometry and mechanical properties for all subjects) was successively computed applying the personalized lever arm of each subject.

RESULTS

A combination of an increase in thoracic kyphosis, cervical lordosis, and pelvic tilt with a loss of lumbar lordosis was observed between the young and the older groups. Sagittal alignment parameters indicated a more forward position as age increased. The lever arm of the CoM above L1 varied from an average of 1 mm backward for the young group, to averages of 10 and 24 mm forward, respectively, for the intermediate and elderly group. As a result, vertebral strength decreased from 2527 N for the young group to 1820 N for the elderly group.

CONCLUSION

The global sagittal alignment modifications observed with age were consistent with the literature. Posture alteration with age reduced vertebral strength significantly in this simplified loading model. Postural alignment seems essential to be considered in the evaluation of osteoporotic patients.

A Review on Multiscale Bone Damage: From the Clinical to the Research Perspective

The investigation of bone damage processes is a crucial point to understand the mechanisms of age-related bone fractures. In order to reduce their impact, early diagnosis is key. The intricate architecture of bone and the complexity of multiscale damage processes make fracture prediction an ambitious goal. This review, supported by a detailed analysis of bone damage physical principles, aims at presenting a critical overview of how multiscale imaging techniques could be used to implement reliable and validated numerical tools for the study and prediction of bone fractures. While macro- and meso-scale imaging find applications in clinical practice, micro- and nano-scale imaging are commonly used only for research purposes, with the objective to extract fragility indexes. Those images are used as a source for multiscale computational damage models. As an example, micro-computed tomography (micro-CT) images in combination with micro-finite element models could shed some light on the comprehension of the interaction between micro-cracks and micro-scale bone features. As future insights, the actual state of technology suggests that these models could be a potential substitute for invasive clinical practice for the prediction of age-related bone fractures. However, the translation to clinical practice requires experimental validation, which is still in progress.

ASSESSMENT OF PARAMETER UNCERTAINTY IN RIGID MUSCULOSKELETAL SIMULATION USING A PROBABILISTIC APPROACH

Experimental investigation coupled with numerical simulations is commonly used for solving multi-physical problems. In the field of biomechanics, in which the aim is to understand the mechanics of living system, the main difficulties are to provide experimental data reflecting the multi-physical behavior of the system of interest. These experimental data are used as input data for numerical simulations to quantify output responses through physical and/or biological laws expressed by constitutive mathematical equations. However, uncertainties on the experimentally available data exist from factors such as human variability and differences in protocols parameters and techniques. Thus, the true values of these data could never be experimentally measured. The objective of this study was to develop a modeling workflow to assess and account for the parameter uncertainty in rigid musculoskeletal simulation. A generic musculoskeletal model was used. Data uncertainties of the right thigh mass, physiological cross-sectional area (pCSA) and muscle tension coefficient of the rectus femoris were accounted to estimate their effect on the joint moment and muscle force computing, respectively. A guideline was developed to fuse data from multiple sources into a sample variation space leading to establish input data distribution. Uncertainty propagation was performed using Monte Carlo and most probable point methods. A high degree of sensitivity of 0.98 was noted for the effect of thigh mass uncertainty on the hip joint moment using inverse dynamics method. A strong deviation of rectus femoris muscle force (around 260 N) was found under effect of pCSA and muscle tension coefficient on the force estimation using static optimization method. Accounting parameter uncertainty into rigid musculoskeletal simulation plays an essential role in the evaluation of the confidence in the model outputs. Thus, simulation outcome may be computed and represented in a more reliable manner with a global range of plausible values.

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

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

Effects of densitometry, material mapping and load estimation uncertainties on the accuracy of patient-specific finite-element models of the scapula

Patient-specific biomechanical models including patient-specific finite-element (FE) models are considered potentially important tools for providing personalized healthcare to patients with musculoskeletal diseases. A multi-step procedure is often needed to generate a patient-specific FE model. As all involved steps are associated with certain levels of uncertainty, it is important to study how the uncertainties of individual components propagate to final simulation results. In this study, we considered a specific case of this problem where the uncertainties of the involved steps were known and the aim was to determine the uncertainty of the predicted strain distribution. The effects of uncertainties of three important components of patient-specific models, including bone density, musculoskeletal loads and the parameters of the material mapping relationship on the predicted strain distributions, were studied. It was found that the number of uncertain components and the level of their uncertainty determine the uncertainty of simulation results. The 'average' uncertainty values were found to be relatively small even for high levels of uncertainty in the components of the model. The 'maximum' uncertainty values were, however, quite high and occurred in the areas of the scapula that are of the greatest clinical relevance. In addition, the uncertainty of the simulation result was found to be dependent on the type of movement analysed, with abduction movements presenting consistently lower uncertainty values than flexion movements.

Prediction of the vertebral strength using a finite element model derived from low-dose biplanar imaging: benefits of subject-specific material properties

STUDY DESIGN

A finite element analysis on osteoporotic vertebrae.

OBJECTIVE

This study aims to validate subject-specific finite element models (FEMs) derived from a low-dose imaging system (EOS, Biospace Med, France) for the prediction of vertebral strength. The vertebrae are submitted to an eccentric compression force leading to compression and anterior bending.

SUMMARY OF BACKGROUND DATA

Given the aging population, osteoporosis and vertebral fractures are a major public health issue. A low bone mineral density (BMD) does not always explain incident fractures, and multifactorial analyses are required. In this context, FEMs based on quantitative computed tomography (QCT) have been proposed to predict vertebral strength in vitro or quantify effects of treatments. However, the clinical use of such a model for the in vivo follow-up of the whole spine is limited by the high-radiation dose induced by QCT and the lying position, which does not allow postural assessment with the same modality.

METHODS

Fourteen vertebrae were modeled using a parametric meshing method. The mesh was subject-specific using geometric parameters computed on the 3-dimensional (3D) reconstructions obtained from the EOS biplanar radiographs. The contribution of cortical bone was taken into account by modeling a cortico-cancellous shell whose properties were derived from experimental data. The effect of subject-specific bone Young's moduli derived from EOS vertebral areal BMD was quantified. The 3D position of the point-of-load application and the 3D orientation of the force was faithfully reproduced in the model to compare the predicted strength and experimental strength under the same loading conditions.

RESULTS

The relative error of prediction decreased from 43% to 16% (2.5 times) when subject-specific mechanical properties, derived from EOS areal BMD, were implemented in the FEM compared with averaged material properties. The resulting subject-specific FEMs predicted vertebral strength with a level of significance close to the QCT-based models (r adjusted = 0.79, root mean square error = 367 N).

CONCLUSION

This work underlines the potential of low-dose biplanar x-ray devices to make subject-specific FEMs for prediction of vertebral strength.

A framework for the study of multiple realizations: the importance of levels of analysis

The brain may undergo functional reorganizations. Selective loss of sensory input or training within a restricted part of a modality cause "shifts" within for instance somatotopic or tonotopic maps. Cross-modal plasticity occurs when input within a modality is absent - e.g., in the congenitally blind. Reorganizations are also found in functional recovery after brain injury. Focusing on such reorganizations, it may be studied whether a cognitive or conscious process can exclusively be mediated by one neural substrate - or may be associated with multiple neural representations. This is typically known as the problem of multiple realization - an essentially empirical issue with wide theoretical implications. This issue may appear to have a simple solution. When, for instance, the symptoms associated with brain injury disappear and the recovery is associated with increased activities within spared regions of the brain, it is tempting to conclude that the processes originally associated with the injured part of the brain are now mediated by an alternative neural substrate. Such a conclusion is, however, not a simple matter. Without a more thorough analysis, it cannot be concluded that a functional recovery of for instance language or attention is necessarily associated with a novel representation of the processes lost to injury. Alternatively, for instance, the recovery may reflect that apparently similar surface phenomena are obtained via dissimilar cognitive mechanisms. In this paper we propose a theoretical framework, which we believe can guide the design and interpretations of studies of post-traumatic recovery. It is essential to distinguish between a number of levels of analysis - including a differentiation between the surface phenomena and the underlying information processing - when addressing, for instance, whether a pre-traumatic and post-traumatically recovered cognitive or conscious process are actually the same. We propose a (somewhat preliminary) system of levels of analysis, which can be applied to such studies.

A framework for the study of multiple realizations: the importance of levels of analysis

The brain may undergo functional reorganizations. Selective loss of sensory input or training within a restricted part of a modality cause "shifts" within for instance somatotopic or tonotopic maps. Cross-modal plasticity occurs when input within a modality is absent - e.g., in the congenitally blind. Reorganizations are also found in functional recovery after brain injury. Focusing on such reorganizations, it may be studied whether a cognitive or conscious process can exclusively be mediated by one neural substrate - or may be associated with multiple neural representations. This is typically known as the problem of multiple realization - an essentially empirical issue with wide theoretical implications. This issue may appear to have a simple solution. When, for instance, the symptoms associated with brain injury disappear and the recovery is associated with increased activities within spared regions of the brain, it is tempting to conclude that the processes originally associated with the injured part of the brain are now mediated by an alternative neural substrate. Such a conclusion is, however, not a simple matter. Without a more thorough analysis, it cannot be concluded that a functional recovery of for instance language or attention is necessarily associated with a novel representation of the processes lost to injury. Alternatively, for instance, the recovery may reflect that apparently similar surface phenomena are obtained via dissimilar cognitive mechanisms. In this paper we propose a theoretical framework, which we believe can guide the design and interpretations of studies of post-traumatic recovery. It is essential to distinguish between a number of levels of analysis - including a differentiation between the surface phenomena and the underlying information processing - when addressing, for instance, whether a pre-traumatic and post-traumatically recovered cognitive or conscious process are actually the same. We propose a (somewhat preliminary) system of levels of analysis, which can be applied to such studies.