Quantification of the effect of osteolytic metastases on bone strain within whole vertebrae using image registration

Influence of Metastatic Bone Lesion Type and Tumor Origin on Human Vertebral Bone Architecture, Matrix Quality, and Mechanical Properties

Metastatic spine disease is incurable, causing increased vertebral fracture risk and severe patient morbidity. Here, we demonstrate that osteolytic, osteosclerotic, and mixed bone metastasis uniquely modify human vertebral bone architecture and quality, affecting vertebral strength and stiffness. Multivariable analysis showed bone metastasis type dominates vertebral strength and stiffness changes, with neither age nor gender having an independent effect. In osteolytic vertebrae, bone architecture rarefaction, lower tissue mineral content and connectivity, and accumulation of advanced glycation end-products (AGEs) affected low vertebral strength and stiffness. In osteosclerotic vertebrae, high trabecular number and thickness but low AGEs, suggesting a high degree of bone remodeling, yielded high vertebral strength. Our study found that bone metastasis from prostate and breast primary cancers differentially impacted the osteosclerotic bone microenvironment, yielding altered bone architecture and accumulation of AGEs. These findings indicate that therapeutic approaches should target the restoration of bone structural integrity. © 2022 American Society for Bone and Mineral Research (ASBMR).

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.

MR-compatible loading device for assessment of heel pad internal tissue displacements under shearing load

In the last decade, the role of shearing loads has been increasingly suspected to play a determinant impact in the formation of deep pressure ulcers. In vivo observations of such deformations are complex to obtain. Previous studies only provide global measurements of such deformations without getting the quantitative values of the loads that generate these deformations. To study the role that shearing loads have in the etiology of heel pressure ulcers, an MR-compatible device for the application of shearing and normal loads was designed. Magnetic resonance imaging is a key feature that allows to monitor deformations of soft tissues after loading in a non-invasive way. Measuring applied forces in an MR-environment is challenging due to the impossibility to use magnetic materials. In our device, forces are applied through the compression of springs made of polylactide. Shearing and normal loads were applied on the plantar skin of the human heel through a flat plate while acquiring MR images. The device materials did not introduce any imaging artifact and allowed for high quality MR deformation measurements of the internal components of the heel. The obtained subject-specific results are an original data set that can be used in validations for Finite Element analysis and therefore contribute to a better understanding of the factors involved in pressure ulcer development.

Improved estimates of strength and stiffness in pathologic vertebrae with bone metastases using CT-derived bone density compared with radiographic bone lesion quality classification

OBJECTIVE

The aim of this study was to compare the ability of 1) CT-derived bone lesion quality (classification of vertebral bone metastases [BM]) and 2) computed CT-measured volumetric bone mineral density (vBMD) for evaluating the strength and stiffness of cadaver vertebrae from donors with metastatic spinal disease.

METHODS

Forty-five thoracic and lumbar vertebrae were obtained from cadaver spines of 11 donors with breast, esophageal, kidney, lung, or prostate cancer. Each vertebra was imaged using microCT (21.4 μm), vBMD, and bone volume to total volume were computed, and compressive strength and stiffness experimentally measured. The microCT images were reconstructed at 1-mm voxel size to simulate axial and sagittal clinical CT images. Five expert clinicians blindly classified the images according to bone lesion quality (osteolytic, osteoblastic, mixed, or healthy). Fleiss' kappa test was used to test agreement among 5 clinical raters for classifying bone lesion quality. Kruskal-Wallis ANOVA was used to test the difference in vertebral strength and stiffness based on bone lesion quality. Multivariable regression analysis was used to test the independent contribution of bone lesion quality, computed vBMD, age, gender, and race for predicting vertebral strength and stiffness.

RESULTS

A low interrater agreement was found for bone lesion quality (κ = 0.19). Although the osteoblastic vertebrae showed significantly higher strength than osteolytic vertebrae (p = 0.0148), the multivariable analysis showed that bone lesion quality explained 19% of the variability in vertebral strength and 13% in vertebral stiffness. The computed vBMD explained 75% of vertebral strength (p < 0.0001) and 48% of stiffness (p < 0.0001) variability. The type of BM affected vBMD-based estimates of vertebral strength, explaining 75% of strength variability in osteoblastic vertebrae (R2 = 0.75, p < 0.0001) but only 41% in vertebrae with mixed bone metastasis (R2 = 0.41, p = 0.0168), and 39% in osteolytic vertebrae (R2 = 0.39, p = 0.0381). For vertebral stiffness, vBMD was only associated with that of osteoblastic vertebrae (R2 = 0.44, p = 0.0024). Age and race inconsistently affected the model's strength and stiffness predictions.

CONCLUSIONS

Pathologic vertebral fracture occurs when the metastatic lesion degrades vertebral strength, rendering it unable to carry daily loads. This study demonstrated the limitation of qualitative clinical classification of bone lesion quality for predicting pathologic vertebral strength and stiffness. Computed CT-derived vBMD more reliably estimated vertebral strength and stiffness. Replacing the qualitative clinical classification with computed vBMD estimates may improve the prediction of vertebral fracture risk.

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.

Biomechanical Properties of Metastatically Involved Osteolytic Bone

PURPOSE OF REVIEW

Skeletal metastasis involves the uncoupling of physiologic bone remodeling resulting in abnormal bone turnover and radical changes in bony architecture, density, and quality. Bone strength assessment and fracture risk prediction are critical in clinical treatment decision-making. This review focuses on bone tissue and structural mechanisms altered by osteolytic metastasis and the resulting changes to its material and mechanical behavior.

RECENT FINDINGS

Both organic and mineral phases of bone tissue are altered by osteolytic metastatic disease, with diminished bone quality evident at multiple length-scales. The mechanical performance of bone with osteolytic lesions is influenced by a combination of tissue-level and structural changes. This review considers the effects of osteolytic metastasis on bone biomechanics demonstrating its negative impact at tissue and structural levels. Future studies need to assess the cumulative impact of cancer treatments on metastatically involved bone quality, and its utility in directing multimodal treatment planning.

Effect of size and location of simulated lytic lesions on the structural properties of human vertebral bodies, a micro-finite element study

Currently, the Spinal Instability Neoplastic Score system is used in clinics to evaluate the risk of fracture in patients with spinal metastases. This method, however, does not always provide a clear guideline due to the complexity in accounting for the effect of metastatic lesions on vertebral stability. The aim of this study was to use a validated micro Finite Element (microFE) modelling approach to analyse the effect of the size and location of lytic metastases on the mechanical properties of human vertebral bodies. Micro Computed Tomography based microFE models were generated with and without lytic lesions simulated as holes within a human vertebral body. Single and multiple lytic lesions were simulated with four different sizes and in five different locations. Bone was assumed homogenous, isotropic and linear elastic, and each vertebra was loaded in axial compression. It was observed that the size of lytic lesions was linearly related with the reduction in structural properties of the vertebral body (reduction of stiffness between 3% and 30% for lesion volume between 4% and 35%). The location of lytic lesions did not show a clear effect on predicted structural properties. Single or multiple lesions with the same volume provided similar results. Locally, there was a homogeneous distribution of axial principal strains among the models with and without lytic lesions. This study highlights the potential of microFE models to study the effect of lesions on the mechanical properties of the human vertebral body.

•

MicroFE models can show the effect of lytic lesions on vertebral properties.

•

The size of the lesions was more critical than the location of the lesions.

•

Lesions affecting the cortical shell had a larger effect on the local strains.

•

Multiple lesions showed a similar effect to single lesions.

Mechanical behavior of metastatic vertebrae are influenced by tissue architecture, mineral content, and organic feature alterations

Diminished vertebral mechanical behavior with metastatic involvement is typically attributed to modified architecture and trabecular bone content. Previous work has identified organic and mineral phase bone quality changes in the presence of metastases, yet limited work exists on the potential influence of such tissue level modifications on vertebral mechanical characteristics. This work seeks to determine correlations between features of bone (structural and tissue level) and mechanical behavior in metastatically involved vertebral bone. It is hypothesized that tissue level properties (mineral and organic) will improve these correlations beyond architectural properties and BMD alone. Twenty-four female athymic rats were inoculated with HeLa or Ace-1 cancer cells lines producing osteolytic (N = 8) or mixed (osteolytic/osteoblastic, N = 7) metastases, respectively. Twenty-one days post-inoculation L1-L3 pathologic vertebral motion segments were excised and μCT imaged. 3D morphometric parameters and axial rigidity of the L2 vertebrae were quantified. Sequential loading and μCT imaging measured progression of failure, stiffness and peak force. Relationships between mechanical testing (whole bone and tissue-level) and tissue-level material property modifications with metastatic involvement were evaluated utilizing linear regression models. Osteolytic involvement reduced vertebral trabecular bone volume, structure, CT-derived axial rigidity, stiffness and failure force compared to healthy controls (N = 9). Mixed metastases demonstrated similar trends. Previously assessed collagen cross-linking and proline-based residues were correlated to mechanical behavior and improved the predictive ability of the regression models. Similarly, collagen organization improved predictive regression models for metastatic bone hardness. This work highlights the importance of both bone content/architecture and organic tissue-level features in characterizing metastatic vertebral mechanics. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:3013-3022, 2018.

Micro Finite Element models of the vertebral body: Validation of local displacement predictions

The estimation of local and structural mechanical properties of bones with micro Finite Element (microFE) models based on Micro Computed Tomography images depends on the quality bone geometry is captured, reconstructed and modelled. The aim of this study was to validate microFE models predictions of local displacements for vertebral bodies and to evaluate the effect of the elastic tissue modulus on model’s predictions of axial forces. Four porcine thoracic vertebrae were axially compressed in situ, in a step-wise fashion and scanned at approximately 39μm resolution in preloaded and loaded conditions. A global digital volume correlation (DVC) approach was used to compute the full-field displacements. Homogeneous, isotropic and linear elastic microFE models were generated with boundary conditions assigned from the interpolated displacement field measured from the DVC. Measured and predicted local displacements were compared for the cortical and trabecular compartments in the middle of the specimens. Models were run with two different tissue moduli defined from microindentation data (12.0GPa) and a back-calculation procedure (4.6GPa). The predicted sum of axial reaction forces was compared to the experimental values for each specimen. MicroFE models predicted more than 87% of the variation in the displacement measurements (R² = 0.87–0.99). However, model predictions of axial forces were largely overestimated (80–369%) for a tissue modulus of 12.0GPa, whereas differences in the range 10–80% were found for a back-calculated tissue modulus. The specimen with the lowest density showed a large number of elements strained beyond yield and the highest predictive errors. This study shows that the simplest microFE models can accurately predict quantitatively the local displacements and qualitatively the strain distribution within the vertebral body, independently from the considered bone types.

Application of digital volume correlation to study the efficacy of prophylactic vertebral augmentation

BACKGROUND

Prophylactic augmentation is meant to reinforce the vertebral body, but in some cases it is suspected to actually weaken it. Past studies only investigated structural failure and the surface strain distribution. To elucidate the failure mechanism of the augmented vertebra, more information is needed about the internal strain distribution. This study aims to measure, for the first time, the full-field three-dimensional strain distribution inside augmented vertebrae in the elastic regime and to failure.

METHODS

Eight porcine vertebrae were prophylactically-augmented using two augmentation materials. They were scanned with a micro-computed tomography scanner (38.8μm voxel resolution) while undeformed, and loaded at 5%, 10%, and 15% compressions. Internal strains (axial, antero-posterior and lateral-lateral components) were computed using digital volume correlation.

FINDINGS

For both augmentation materials, the highest strains were measured in the regions adjacent to the injected cement mass, whereas the cement-interdigitated-bone was less strained. While this was already visible in the elastic regime (5%), it was a predictor of the localization of failure, which became visible at higher degrees of compression (10% and 15%), when failure propagated across the trabecular bone. Localization of high strains and failure was consistent between specimens, but different between the cement types.

INTERPRETATION

This study indicated the potential of digital volume correlation in measuring the internal strain (elastic regime) and failure in augmented vertebrae. While the cement-interdigitated region becomes stiffer (less strained), the adjacent non-augmented trabecular bone is affected by the stress concentration induced by the cement mass. This approach can help establish better criteria to improve vertebroplasty.

Mechanical assessment of the effects of metastatic lytic defect on the structural response of human thoracolumbar spine

To investigate the effects of a clinical lytic defect on the structural response of human thoracolumbar functional spinal unit. A novel CT-compatible mechanical test system was used to image the deformation of a T12-L1 motion segment and measure the change in strain response under compressive loads ranging from 50 to 750 N. A lytic lesion (LM) with cortex involvement (33% by volume) was introduced to the upper vertebral body and the CT experiments were repeated. Finite element models, established from the CT volumes, were used to investigate the defect's effects on the structural response and the state of principal and shear stresses within the affected and adjacent vertebrae. The lytic lesion resulted in severe loss of the vertebral structural competence, resulting in significant, non-linear, and asymmetric increase in the experimentally measured strains and computed stresses within both vertebrae (p < 0.01). At the cortex, the tensile strains were significantly increased, while compressive strains significantly decreased, (p < 0.05). Both the vertebral bone and cortex regions adjacent to the defect showed significant increase in computed compressive, tensile, and shear stresses (p < 0.01). Changes in stress and strain distribution within the affected and adjacent vertebral bone and the experimentally observed bulging and buckling of the vertebral cortices suggested that initiation of catastrophic vertebral failure may occur under load magnitudes encountered in daily living. Although the effect of LM on the global deformation of the spine was well-predicted, our results show that FE predictions of local strain changes must be carefully assessed for clinical relevance. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1808-1819, 2016.

Post-euthanasia micro-computed tomography-based strain analysis is able to represent quasi-static in vivo behavior of whole vertebrae

Three-dimensional image-based strain measurement in whole bones allows representation of physiological, albeit quasi-static, loading conditions. However, such work to date has been limited to specimens postmortem. The main purpose of this study is to verify the efficacy of deformable image registration of post-euthanasia strain to characterize the in vivo mechanical behavior of rat vertebrae. A micro-computed tomography-compatible custom loading device was used to apply 75 N load to a three-level caudal motion segment of a healthy rat. Loaded and unloaded micro-computed tomography scans were acquired in vivo and post-sacrifice. A micro-computed tomography-based deformable image registration algorithm was used to calculate vertebral strains live and post-euthanasia. No significant difference was found in the in vivo strains (-0.011 ± 0.001) and ex vivo strains (-0.012 ± 0.001) obtained from the comparisons of loaded and unloaded images (p = 0.3). Comparisons between unloaded-unloaded and loaded-loaded scans yielded significantly lower axial strains, representing the error of the method. Qualitatively, high strains were observed adjacent to growth plate regions in evaluating the loaded-unloaded images. Strain patterns in the loaded-loaded and unloaded-unloaded scans were inconsistent as would be expected in representing noise. Overall, live and dead loaded to unloaded comparisons yielded similar strain patterns and magnitudes. Point-wise differences in axial strain fields also supported this observation. This study demonstrated a proof of concept, suggesting that post-euthanasia micro-computed tomography-based strain analysis is able to represent the in vivo quasi-static behavior of rat tail vertebrae.

μFE models can represent microdamaged regions of healthy and metastatically involved whole vertebrae identified through histology and contrast enhanced μCT imaging

Micro-damage formation within the skeleton is an important stimulant for bone remodeling, however abnormal build-up of micro-damage can lead to skeletal fragility. In this study, µCT imaging based micro finite element (μFE) models were used to evaluate tissue level damage criteria in whole healthy and metastatically-involved vertebrae. T13-L2 spinal segments were excised from osteolytic (n=3) and healthy (n=3) female athymic rnu/rnu rats. Osteolytic metastasis was generated by intercardiac injection of HeLa cancer cells. Micro-mechanical axial loading was applied to the spinal motion segments under μCT imaging. Vertebral samples underwent BaSO4 staining and sequential calcein/fuchsin staining to identify load induced micro-damage. μCT imaging was used generate specimen specific μFE models of the healthy and osteolytic whole rat vertebrae. Model boundary conditions were generated through deformable image registration of loaded and unloaded scans. Elevated stresses and strains were detected in regions of micro-damage identified through histological and BaSO4 staining within healthy and osteolytic vertebral models, as compared to undamaged regions. Additionally, damaged regions of metastatic vertebrae experienced significantly higher local stresses and strains than those in the damaged regions of healthy specimens. Areas identified by BaSO4 staining, however, yielded lower levels of stress and strain in damaged and undamaged regions of healthy and metastatic vertebrae as compared to fuschin staining. The multimodal (experimental, image-based and computational) techniques used in this study demonstrated the ability of local stresses and strains computed through µFE analysis to identify trabecular micro-damage, that can be applied to biomechanical analyses of healthy and diseased whole bones.

Effect of the metastatic defect on the structural response and failure process of human vertebrae: an experimental study

BACKGROUND

Pathologic vertebral fractures are associated with intractable pain, loss of function and high morbidity in patients with metastatic spine disease. However, the failure mechanisms of vertebrae with lytic defects and the failed vertebrae's ability to retain load carrying capacity remain unclear.

METHODS

Eighteen human thoracic and lumbar vertebrae with simulated uncontained bone defects were tested under compression-bending loads to failure. Failure was defined as 50% reduction in vertebral body height. The vertebrae were allowed to recover under load and re-tested to failure using the initial criteria. Repeated measure ANOVA was used to test for changes in strength and stiffness parameters.

FINDINGS

Vertebral failure occurred via buckling and fracture of the cortex around the defect, followed by collapse of the defect region. Compared to the intact vertebrae, the failed vertebrae exhibited a significant loss in compressive strength (59%, p<0.001), stiffness (53%, p<0.05) and flexion (70%, p<0.01) strength. Significant reduction in anterior-posterior shear (strength (63%, p<0.01) and stiffness (67%, p<0.01)) and lateral bending strength (134%, p<0.05) were similarly recorded. In the intact vertebrae, apart from flexion strength (r(2)=0.63), both compressive and anterior-posterior shear strengths were weakly correlated with their stiffness parameters (r(2)=0.24 and r(2)=0.31). By contrast, in the failed vertebrae, these parameters were strongly correlated, (r(2)=0.91, r(2)=0.86, and r(2)=0.92, p<0.001 respectively).

INTERPRETATION

Failure of the vertebral cortex at the defect site dominated the initiation and progression of vertebral failure with the vertebrae failing via a consolidation process of the vertebral bone. Once failed, the vertebrae showed remarkable loss of load carrying capacity.

EXPERIMENTAL METHODS FOR THE BIOMECHANICAL INVESTIGATION OF THE HUMAN SPINE: A REVIEW

In vitro mechanical testing of spinal specimens is extremely important to better understand the biomechanics of the healthy and diseased spine, fracture, and to test/optimize surgical treatment. While spinal testing has extensively been carried out in the past four decades, testing methods are quite diverse. This paper aims to provide a critical overview of the in vitro methods for mechanical testing the human spine at different scales. Specimens of different type are used, according to the aim of the study: spine segments (two or more adjacent vertebrae) are used both to investigate the spine kinematics, and the mechanical properties of the spine components (vertebrae, ligaments, discs); single vertebrae (whole vertebra, isolated vertebral body, or vertebral body without endplates) are used to investigate the structural properties of the vertebra itself; core specimens are extracted to test the mechanical properties of the trabecular bone at the tissue-level; mechanical properties of spine soft tissue (discs, ligaments, spinal cord) are measured on isolated elements, or on tissue specimens. Identification of consistent reference frames is still a debated issue. Testing conditions feature different pre-conditioning and loading rates, depending on the simulated action. Tissue specimen preservation is a very critical issue, affecting test results. Animal models are often used as a surrogate. However, because of different structure and anatomy, extreme caution is required when extrapolating to the human spine. In vitro loading conditions should be based on reliable in vivo data. Because of the high complexity of the spine, such information (either through instrumented implants or through numerical modeling) is currently unsatisfactory. Because of the increasing ability of computational models in predicting biomechanical properties of musculoskeletal structures, a synergy is possible (and desirable) between in vitro experiments and numerical modeling. Future perspectives in spine testing include integration of mechanical and structural properties at different dimensional scales (from the whole-body-level down to the tissue-level) so that organ-level models (which are used to predict the most relevant phenomena such as fracture) include information from all dimensional scales.