Heterogeneous Strain Distribution in the Subchondral Bone of Human Osteoarthritic Femoral Heads, Measured with Digital Volume Correlation

A New Microarchitecture-Based Parameter to Predict the Micromechanical Properties of Bone Allografts

Scaffolds are an essential component of bone tissue engineering. They provide support and create a physiological environment for cells to proliferate and differentiate. Bone allografts extracted from human donors are promising scaffolds due to their mechanical and structural characteristics. Bone microarchitecture is well known to be an important determinant of macroscopic mechanical properties, but its role at the microscopic, i.e., the trabeculae level is still poorly understood. The present study investigated linear correlations between microarchitectural parameters obtained from X-ray computed tomography (micro-CT) images of bone allografts, such as bone volume fraction (BV/TV), degree of anisotropy (DA), or ellipsoid factor (EF), and micromechanical parameters derived from micro-finite element calculations, such as mean axial strain (ε_z) and strain energy density (W_e). DAEF, a new parameter based on a linear combination of the two microarchitectural parameters DA and EF, showed a strong linear correlation with the bone mechanical characteristics at the microscopic scale. Our results concluded that the spatial distribution and the plate-and-rod structure of trabecular bone are the main determinants of the mechanical properties of bone at the microscopic level. The DAEF parameter could, therefore, be used as a tool to predict the level of mechanical stimulation at the local scale, a key parameter to better understand and optimize the mechanism of osteogenesis in bone tissue engineering.

Bone metastases do not affect the measurement uncertainties of a global digital volume correlation algorithm

Introduction: Measurement uncertainties of Digital Volume Correlation (DVC) are influenced by several factors, like input images quality, correlation algorithm, bone type, etc. However, it is still unknown if highly heterogeneous trabecular microstructures, typical of lytic and blastic metastases, affect the precision of DVC measurements.

Methods: Fifteen metastatic and nine healthy vertebral bodies were scanned twice in zero-strain conditions with a micro-computed tomography (isotropic voxel size = 39 μm). The bone microstructural parameters (Bone Volume Fraction, Structure Thickness, Structure Separation, Structure Number) were calculated. Displacements and strains were evaluated through a global DVC approach (BoneDVC). The relationship between the standard deviation of the error (SDER) and the microstructural parameters was investigated in the entire vertebrae. To evaluate to what extent the measurement uncertainty is influenced by the microstructure, similar relationships were assessed within sub-regions of interest.

Results: Higher variability in the SDER was found for metastatic vertebrae compared to the healthy ones (range 91-1030 με versus 222–599 με). A weak correlation was found between the SDER and the Structure Separation in metastatic vertebrae and in the sub-regions of interest, highlighting that the heterogenous trabecular microstructure only weakly affects the measurement uncertainties of BoneDVC. No correlation was found for the other microstructural parameters. The spatial distribution of the strain measurement uncertainties seemed to be associated with regions with reduced greyscale gradient variation in the microCT images.

Discussion: Measurement uncertainties cannot be taken for granted but need to be assessed in each single application of the DVC to consider the minimum unavoidable measurement uncertainty when interpreting the results.

Comparison of linear and nonlinear stepwise μFE displacement predictions to digital volume correlation measurements of trabecular bone biopsies

Digital volume correlation (DVC) enables to evaluate the ability of μFE models in predicting experimental results on the mesoscale. In this study predicted displacement fields of three different linear and materially nonlinear μFE simulation methods were compared to DVC measured displacement fields at specific load steps in the elastic regime (Step_El) and after yield (Step_Ult). Five human trabecular bone biopsies from a previous study were compressed in several displacement steps until failure. At every compression step, μCT images (resolution: 36 μm) were recorded. A global DVC algorithm was applied to compute the displacement fields at all loading steps. The unloaded 3D images were then used to generate homogeneous, isotropic, linear and materially nonlinear μFE models. Three different μFE simulation methods were used: linear (L), nonlinear (NL), and nonlinear stepwise (NLS). Regarding L and NL, the boundary conditions were derived from the interpolated displacement fields at Step_El and Step_Ult, while for the NLS method nonlinear changes of the boundary conditions of the experiments were captured using the DVC displacement field of every available load step until Step_El and Step_Ult. The predicted displacement fields of all μFE simulation methods were in good agreement with the DVC measured displacement fields (individual specimens: R²>0.83 at Step_El and R²>0.59 at Step_Ult; pooled data: R²>0.97 at Step_El and R²>0.92 at Step_Ult). At Step_El, all three simulation methods showed similar intercepts, slopes, and coefficients of determination while the nonlinear μFE models improved the prediction of the displacement fields slightly in all Cartesian directions at Step_Ult (individual specimens: L: R²>0.59 and NL, NLS: R²>0.68; pooled data: L: R²>0.92 and NL, NLS: R²>0.94). Damaged/overstrained elements in L, NL, and NLS occurred at similar locations but the number of overstrained elements was overestimated when using the L simulation method. Considering the increased solving time of the nonlinear μFE models as well as the acceptable performance in displacement prediction of the linear μFE models, one can conclude that for similar use cases linear μFE models represent the best compromise between computational effort and accuracy of the displacement field predictions.

Full-field experimental analysis of the influence of microstructural parameters on the mechanical properties of humeral head trabecular bone

Understanding the mechanical properties of trabecular bone within the metaphysis of the proximal humerus is becoming increasingly important for the design of humeral head joint replacement components that prioritize bone preservation. The aim of this study was to perform full-field mechanical testing methods on isolated trabecular bone cores from the humeral head to experimentally measure the local magnitude of strain before macroscopic failure and to characterize the ultimate strength of each core. Isolated cubic trabecular bone cores were extracted from the center of humeral head osteotomies retrieved from (1) patients with end-stage osteoarthritis (OA) undergoing total shoulder arthroplasty (TSA) and (2) normal nonpathologic cadaveric humeral heads. A custom computed tomography (CT)-compatible loading device was used to perform compressive mechanical testing. For 10 of the OA specimens, stepwise loading was performed directly within a microCT scanner and digital volume correlation (DVC) was used to measure full-field strains throughout the trabecular structure. A higher variability in ultimate strength was measured for the trabecular cores retrieved from OA humeral heads (range: 2.8-7.6 MPa) compared to the normal cadaveric humeral heads (range: 2.2-5.4 MPa), but no statistically significant difference between the groups was found (p = 0.06). Ultimate strength was strongly correlated with bone volume fraction (OA r²  = 0.72; normal r²  = 0.76) and bone mineral content (OA r²  = 0.79; normal r²  = 0.77). At the trabecular level, 95th percentile of third principal strains, measured at a subvolume size of 152 µm, exceeded 19,000 µε for each of the 10 specimens (range: -19,551 to -36,535 µε) before macroscopic failure of the cores occured. No strong linear correlations (r²  ≥ 0.50) were found between the median or 95th percentile of DVC third principal strain and the corresponding morphometric parameters of each individual bone core. The results of this study indicate that bone volume fraction and bone mineral content heavily influence the apparent ultimate strength of trabecular bone cores collected from OA patients undergoing TSA. Clinical significance: The strong correlations observed within this study further emphasize the importance of considering bone mineral content or bone volume fraction measurements in assessing the localized risk of trabecular bone fracture for orthopedic applications.

A practical guide for in situ mechanical testing of musculoskeletal tissues using synchrotron tomography

Musculoskeletal tissues are complex hierarchical materials where mechanical response is linked to structural and material properties at different dimensional levels. Therefore, high-resolution three-dimensional tomography is very useful for assessing tissue properties at different scales. In particular, Synchrotron Radiation micro-Computed Tomography (SR-microCT) has been used in several applications to analyze the structure of bone and biomaterials. In the past decade the development of digital volume correlation (DVC) algorithms applied to SR-microCT images and its combination with in situ mechanical testing (four-dimensional imaging) have allowed researchers to visualise, for the first time, the deformation of musculoskeletal tissues and their interaction with biomaterials under different loading scenarios. However, there are several experimental challenges that make these measurements difficult and at high risk of failure. Challenges relate to sample preparation, imaging parameters, loading setup, accumulated tissue damage for multiple tomographic acquisitions, reconstruction methods and data processing. Considering that access to SR-microCT facilities is usually associated with bidding processes and long waiting times, the failure of these experiments could notably slow down the advancement of this research area and reduce its impact. Many of the experimental failures can be avoided with increased experience in performing the tests and better guidelines for preparation and execution of these complex experiments; publication of negative results could help interested researchers to avoid recurring mistakes. Therefore, the goal of this article is to highlight the potential and pitfalls in the design and execution of in situ SR-microCT experiments, involving multiple scans, of musculoskeletal tissues for the assessment of their structural and/or mechanical properties. The advice and guidelines that follow should improve the success rate of this type of experiment, allowing the community to reach higher impact more efficiently.

Characterizing the Mechanical Behavior of Bone and Bone Surrogates in Compression Using pQCT

Many axial and appendicular skeleton bones are subjected to repetitive loading during daily activities. Until recently, the structural analysis of fractures has been limited to 2D sections, and the dynamic assessment of fracture progression has not been possible. The structural failure was analyzed using step-wise micro-compression combined with time-lapsed micro-computed tomographic imaging. The structural failure was investigated in four different sample materials (two different bone surrogates, lumbar vertebral bodies from bovine and red deer). The samples were loaded in different force steps based on uniaxial compression tests. The micro-tomography images were used to create three-dimensional models from which various parameters were calculated that provide information about the structure and density of the samples. By superimposing two 3D images and calculating the different surfaces, it was possible to precisely analyze which trabeculae failed in which area and under which load. According to the current state of the art, bone mineral density is usually used as a value for bone quality, but the question can be raised as to whether other values such as trabecular structure, damage accumulation, and bone mineralization can predict structural competence better than bone mineral density alone.

Digital volume correlation for the characterization of musculoskeletal tissues: Current challenges and future developments

Biological tissues are complex hierarchical materials, difficult to characterise due to the challenges associated to the separation of scale and heterogeneity of the mechanical properties at different dimensional levels. The Digital Volume Correlation approach is the only image-based experimental approach that can accurately measure internal strain field within biological tissues under complex loading scenarios. In this minireview examples of DVC applications to study the deformation of musculoskeletal tissues at different dimensional scales are reported, highlighting the potential and challenges of this relatively new technique. The manuscript aims at reporting the wide breath of DVC applications in the past 2 decades and discuss future perspective for this unique technique, including fast analysis, applications on soft tissues, high precision approaches, and clinical applications.

MicroFE models of porcine vertebrae with induced bone focal lesions: Validation of predicted displacements with digital volume correlation

The evaluation of the local mechanical behavior as a result of metastatic lesions is fundamental for the characterization of the mechanical competence of metastatic vertebrae. Micro finite element (microFE) models have the potential of addressing this challenge through laboratory studies but their predictions of local deformation due to the complexity of the bone structure compromized by the lesion must be validated against experiments. In this study, the displacements predicted by homogeneous, linear and isotropic microFE models of vertebrae were validated against experimental Digital Volume Correlation (DVC) measurements. Porcine spine segments, with and without mechanically induced focal lesions, were tested in compression within a micro computed tomography (microCT) scanner. The displacement within the bone were measured with an optimized global DVC approach (BoneDVC). MicroFE models of the intact and lesioned vertebrae, including or excluding the growth plates, were developed from the microCT images. The microFE and DVC boundary conditions were matched. The displacements measured by the DVC and predicted by the microFE along each Cartesian direction were compared. The results showed an excellent agreement between the measured and predicted displacements, both for intact and metastatic vertebrae, in the middle of the vertebra, in those cases where the structure was not loaded beyond yield (0.69 < R² < 1.00). Models with growth plates showed the worst correlations (0.02 < R² < 0.99), while a clear improvement was observed if the growth plates were excluded (0.56 < R² < 1.00). In conclusion, these simplified models can predict complex displacement fields in the elastic regime with high reliability, more complex non-linear models should be implemented to predict regions with high deformation, when the bone is loaded beyond yield.

3D strain analysis of trabecular bone within the osteoarthritic humeral head subjected to stepwise compressive loads

Understanding the local mechanical properties of trabecular bone at the humeral head-neck junction is essential for the safe design of stemless humeral head implants. Recent advancements in mechanical testing coupled with volumetric imaging have allowed for the ability to quantify full-field strain distributions throughout trabecular bone. Within this study, digital volume correlation (DVC) was applied to micro-computed tomography images to investigate the local load carrying characteristics of trabecular bone within osteoarthritic (OA) humeral heads subjected to stepwise loading. A multi-pegged indenter was used to transfer loads from a custom-fabricated loading apparatus to trabecular bone on the resection surface of OA humeral head osteotomies retrieved from patients undergoing total shoulder arthroplasty (TSA). In regions of trabecular bone that eventually fractured, third principal strains were significantly higher (95^th percentile third principal strain = -12,558 μstrain, p < 0.001) compared to regions that did not fracture (95^th percentile third principal strain = -7,806 μstrain). As well, bone volume fraction (p = 0.012), trabecular separation (p = 0.014), and trabecular number (p = 0.007) were found to influence the likelihood of trabecular bone fracture. Collectively, this work has led to a deeper understanding of the local load carrying characteristics of trabecular bone specific to patients receiving TSA for osteoarthritis.

A novel specimen shape for measurement of linear strain fields by means of digital image correlation

Strains on the surface of engineering structures or biological tissues are non-homogeneous. These strain fields can be captured by means of Digital Image Correlation (DIC). However, DIC strain field measurements are prone to noise and filtering of these fields influences measured strain gradients. This study aims to design a novel tensile test specimen showing two linear gradients, to measure full-field linear strain measurements on the surface of test specimens, and to investigate the accuracy of DIC strain measurements globally (full-field) and locally (strain gauges' positions), with and without filtering of the DIC strain fields. Three materials were employed for this study: aluminium, polymer, and bovine bone. Normalized strain gradients were introduced that are load independent and evaluated at two local positions showing 3.6 and 6.9% strain change per mm. Such levels are typically found in human bones. At these two positions, two strain gauges were applied to check the experimental strain magnitudes. A third strain gauge was applied to measure the strain in a neutral position showing no gradient. The accuracy of the DIC field measurement was evaluated at two deformation stages (at [Formula: see text] 500 and 1750 μstrain) using the root mean square error (RMSE). The RMSE over the two linear strain fields was less than 500 μstrain for both deformation stages and all materials. Gaussian low-pass filter (LPF) reduced the DIC noise between 25% and 64% on average. As well, filtering improved the accuracy of the local normalized strain gradients measurements with relative difference less than 20% and 12% for the high- and low-gradient, respectively. In summary, a novel specimen shape and methodological approach are presented which are useful for evaluating and improving the accuracy of the DIC measurement where non-homogeneous strain fields are expected such as on bone tissue due to their hierarchical structure.

A novel approach to evaluate the effects of artificial bone focal lesion on the three-dimensional strain distributions within the vertebral body

The spine is the first site for incidence of bone metastasis. Thus, the vertebrae have a high potential risk of being weakened by metastatic tissues. The evaluation of strength of the bone affected by the presence of metastases is fundamental to assess the fracture risk. This work proposes a robust method to evaluate the variations of strain distributions due to artificial lesions within the vertebral body, based on in situ mechanical testing and digital volume correlation. Five porcine vertebrae were tested in compression up to 6500N inside a micro computed tomography scanner. For each specimen, images were acquired before and after the application of the load, before and after the introduction of the artificial lesions. Principal strains were computed within the bone by means of digital volume correlation (DVC). All intact specimens showed a consistent strain distribution, with peak minimum principal strain in the range -1.8% to -0.7% in the middle of the vertebra, demonstrating the robustness of the method. Similar distributions of strains were found for the intact vertebrae in the different regions. The artificial lesion generally doubled the strain in the middle portion of the specimen, probably due to stress concentrations close to the defect. In conclusion, a robust method to evaluate the redistribution of the strain due to artificial lesions within the vertebral body was developed and will be used in the future to improve current clinical assessment of fracture risk in metastatic spines.