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

Computational homogenization of histological microstructures in human prostate tissue: Heterogeneity, anisotropy and tension-compression asymmetry

Human prostatic tissue exhibits complex mechanical behaviour due to its multiphasic, heterogeneous nature, with hierarchical microstructures involving epithelial compartments, acinar lumens and stromal tissue all interconnected in complex networks. This study aims to establish a computational homogenization framework for quantifying the mechanical behaviour of prostate tissue, considering its multiphasic heterogeneous microstructures and the mechanical characteristics of tissue constituents. Representative tissue microstructure models were reconstructed from high-resolution histology images. Parametric studies on the mechanical properties of the tissue constituents, particularly the fibre-reinforced hyper-elasticity of the stromal tissue, were carried out to investigate their effects on the apparent tissue properties. These were then benchmarked against the experimental data reported in literature. Results showed significant anisotropy, both structural and mechanical, and tension-compression asymmetry of the apparent behaviours of the prostatic tissue. Strong correlation with the key microstructural indices such as area fractions of tissue constituents and the tissue fabric tensor was also observed. The correlation between the stromal tissue orientation and the principal directions of the apparent properties suggested an essential role of stromal tissue in determining the directions of anisotropy and the compression-tension asymmetry characteristics in normal human prostatic tissue. This work presented a homogenization and histology-based computational approach to characterize the apparent mechanical behaviours of human prostatic or other similar glandular tissues, with the ultimate aim of assessing how pathological conditions (e.g., prostate cancer and benign prostatic hyperplasia) could affect the tissue mechanical properties in a future study.

The role of bone metastases on the mechanical competence of human vertebrae

Spine is the most common site for bone metastases. The evaluation of the mechanical competence and failure location in metastatic vertebrae is a biomechanical and clinical challenge. Little is known about the failure behaviour of vertebrae with metastatic lesions. The aim of this study was to use combined micro-Computed Tomography (microCT) and time-lapsed mechanical testing to reveal the failure location in metastatic vertebrae. Fifteen spine segments, each including a metastatic and a radiologically healthy vertebra, were tested in compression up to failure within a microCT. Volumetric strains were measured using Digital Volume Correlation. The images of undeformed and deformed specimens were overlapped to identify the failure location. Vertebrae with lytic metastases experienced the largest average compressive strains (median ± standard deviation: -8506 ± 4748microstrain), followed by the vertebrae with mixed metastases (-7035 ± 15605microstrain), the radiologically healthy vertebrae (-5743 ± 5697microstrain), and the vertebrae with blastic metastases (-3150 ± 4641microstrain). Strain peaks were localised within and nearby the lytic lesions or around the blastic tissue. Failure between the endplate and the metastasis was identified in vertebrae with lytic metastases, whereas failure localised around the metastasis in vertebrae with blastic lesions. This study showed for the first time the role of metastases on the vertebral internal deformations. While lytic lesions lead to failure of the metastatic vertebra, vertebrae with blastic metastases are more likely to induce failure in the adjacent vertebrae. Nevertheless, every metastatic lesion affects the vertebral deformation differently, making it essential to assess how the lesion affects the bone microstructure. These results suggest that the properties of the lesion (type, size, location within the vertebral body) should be considered when developing clinical tools to predict the risk of fracture in patients with metastatic lesions.

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.

Osteolytic vs. Osteoblastic Metastatic Lesion: Computational Modeling of the Mechanical Behavior in the Human Vertebra after Screws Fixation Procedure

Metastatic lesions compromise the mechanical integrity of vertebrae, increasing the fracture risk. Screw fixation is usually performed to guarantee spinal stability and prevent dramatic fracture events. Accordingly, predicting the overall mechanical response in such conditions is critical to planning and optimizing surgical treatment. This work proposes an image-based finite element computational approach describing the mechanical behavior of a patient-specific instrumented metastatic vertebra by assessing the effect of lesion size, location, type, and shape on the fracture load and fracture patterns under physiological loading conditions. A specific constitutive model for metastasis is integrated to account for the effect of the diseased tissue on the bone material properties. Computational results demonstrate that size, location, and type of metastasis significantly affect the overall vertebral mechanical response and suggest a better way to account for these parameters in estimating the fracture risk. Combining multiple osteolytic lesions to account for the irregular shape of the overall metastatic tissue does not significantly affect the vertebra fracture load. In addition, the combination of loading mode and metastasis type is shown for the first time as a critical modeling parameter in determining fracture risk. The proposed computational approach moves toward defining a clinically integrated tool to improve the management of metastatic vertebrae and quantitatively evaluate fracture risk.

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.

Finite Element Analysis of Osteoporotic and Osteoblastic Vertebrae and Its Association With the Proton Density Fat Fraction From Chemical Shift Encoding-Based Water-Fat MRI - A Preliminary Study

PURPOSE

Osteoporosis is prevalent and entails alterations of vertebral bone and marrow. Yet, the spine is also a common site of metastatic spread. Parameters that can be non-invasively measured and could capture these alterations are the volumetric bone mineral density (vBMD), proton density fat fraction (PDFF) as an estimate of relative fat content, and failure displacement and load from finite element analysis (FEA) for assessment of bone strength. This study's purpose was to investigate if osteoporotic and osteoblastic metastatic changes in lumbar vertebrae can be differentiated based on the abovementioned parameters (vBMD, PDFF, and measures from FEA), and how these parameters correlate with each other.

MATERIALS AND METHODS

Seven patients (3 females, median age: 77.5 years) who received 3-Tesla magnetic resonance imaging (MRI) and multi-detector computed tomography (CT) of the lumbar spine and were diagnosed with either osteoporosis (4 patients) or diffuse osteoblastic metastases (3 patients) were included. Chemical shift encoding-based water-fat MRI (CSE-MRI) was used to extract the PDFF, while vBMD was extracted after automated vertebral body segmentation using CT. Segmentation masks were used for FEA-based failure displacement and failure load calculations. Failure displacement, failure load, and PDFF were compared between patients with osteoporotic vertebrae versus patients with osteoblastic metastases, considering non-fractured vertebrae (L1-L4). Associations between those parameters were assessed using Spearman correlation.

RESULTS

Median vBMD was 59.3 mg/cm³ in osteoporotic patients. Median PDFF was lower in the metastatic compared to the osteoporotic patients (11.9% vs. 43.8%, p=0.032). Median failure displacement and failure load were significantly higher in metastatic compared to osteoporotic patients (0.874 mm vs. 0.348 mm, 29,589 N vs. 3,095 N, p=0.034 each). A strong correlation was noted between PDFF and failure displacement (rho -0.679, p=0.094). A very strong correlation was noted between PDFF and failure load (rho -0.893, p=0.007).

CONCLUSION

PDFF as well as failure displacement and load allowed to distinguish osteoporotic from diffuse osteoblastic vertebrae. Our findings further show strong associations between PDFF and failure displacement and load, thus may indicate complimentary pathophysiological associations derived from two non-invasive techniques (CSE-MRI and CT) that inherently measure different properties of vertebral bone and marrow.

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.

Type, size, and position of metastatic lesions explain the deformation of the vertebrae under complex loading conditions

BACKGROUND

Bone metastases may lead to spine instability and increase the risk of fracture. Scoring systems are available to assess critical metastases, but they lack specificity, and provide uncertain indications over a wide range, where most cases fall. The aim of this work was to use a novel biomechanical approach to evaluate the effect of lesion type, size, and location on the deformation of the metastatic vertebra.

METHOD

Vertebrae with metastases were identified from 16 human spines from a donation programme. The size and position of the metastases, and the Spine Instability Neoplastic Score (SINS) were evaluated from clinical Quantitative Computed Tomography images. Thirty-five spine segments consisting of metastatic vertebrae and adjacent healthy controls were biomechanically tested in four different loading conditions. The strain distribution over the entire vertebral bodies was measured with Digital Image Correlation. Correlations between the features of the metastasis (type, size, position and SINS) and the deformation of the metastatic vertebrae were statistically explored.

RESULTS

The metastatic type (lytic, blastic, mixed) characterizes the vertebral behaviour (Kruskal-Wallis, p = 0.04). In fact, the lytic metastases showed more critical deformation compared to the control vertebrae (average: 2-fold increase, with peaks of 14-fold increase). By contrast, the vertebrae with mixed or blastic metastases did not show a clear trend, with deformations similar or lower than the controls. Once the position of the lytic lesion with respect to the loading direction was taken into account, the size of the lesion was significantly correlated with the perturbation to the strain distribution (r² = 0.72, p < 0.001). Conversely, the SINS poorly correlated with the mechanical evidence, and only in case of lytic lesions (r² = 0.25, p < 0.0001).

CONCLUSION

These results highlight the relevance of the size and location of the lytic lesion, which are marginally considered in the current clinical scoring systems, in driving the spinal biomechanical instability. The strong correlation with the biomechanical evidence indicates that these parameters are representative of the mechanical competence of the vertebra. The improved explanatory power compared to the SINS suggests including them in future guidelines for the clinical practice.

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.

Conventional finite element models estimate the strength of metastatic human vertebrae despite alterations of the bone's tissue and structure

INTRODUCTION

Pathologic vertebral fractures are a major clinical concern in the management of cancer patients with metastatic spine disease. These fractures are a direct consequence of the effect of bone metastases on the anatomy and structure of the vertebral bone. The goals of this study were twofold. First, we evaluated the effect of lytic, blastic and mixed (both lytic and blastic) metastases on the bone structure, on its material properties, and on the overall vertebral strength. Second, we tested the ability of bone mineral content (BMC) measurements and standard FE methodologies to predict the strength of real metastatic vertebral bodies.

METHODS

Fifty-seven vertebral bodies from eleven cadaver spines containing lytic, blastic, and mixed metastatic lesions from donors with breast, esophageal, kidney, lung, or prostate cancer were scanned using micro-computed tomography (μCT). Based on radiographic review, twelve vertebrae were selected for nanoindentation testing, while the remaining forty-five vertebrae were used for assessing their compressive strength. The μCT reconstruction was exploited to measure the vertebral BMC and to establish two finite element models. 1) a micro finite element (μFE) model derived at an image resolution of 24.5 μm and 2) homogenized FE (hFE) model derived at a resolution of 0.98 mm. Statistical analyses were conducted to measure the effect of the bone metastases on BV/TV, indentation modulus (E_it), ratio of plastic/total work (W_Pl/W_tot), and in vitro vertebral strength (F_exp). The predictive value of BMC, μFE stiffness, and hFE strength were evaluated against the in vitro measurements.

RESULTS

Blastic vertebral bodies exhibit significantly higher BV/TV compared to the mixed (p = 0.0205) and lytic (p = 0.0216) vertebral bodies. No significant differences were found between lytic and mixed vertebrae (p = 0.7584). Blastic bone tissue exhibited a 5.8% lower median E_it (p< 0.001) and a 3.3% lower median W_pl/W_tot (p<0.001) compared to non-involved bone tissue. No significant differences were measured between lytic and non-involved bone tissues. F_exp ranged from 1.9 to 13.8 kN, was strongly associated with hFE strength (R²=0.78, p< 0.001) and moderately associated with BMC (R²=0.66, p< 0.001) and μFE stiffness (R²=0.66, p< 0.001), independently of the lesion type.

DISCUSSION

Our findings show that tumour-induced osteoblastic metastases lead to slightly, but significantly lower bone tissue properties compared to controls, while osteolytic lesions appear to have a negligible impact. These effects may be attributed to the lower mineralization and woven nature of bone forming in blastic lesions whilst the material properties of bone in osteolytic vertebrae appeared little changed. The moderate association between BMC- and FE-based predictions to fracture strength suggest that vertebral strength is affected by the changes of bone mass induced by the metastatic lesions, rather than altered tissue properties. In a broader context, standard hFE approaches generated from CTs at clinical resolution are robust to the lesion type when predicting vertebral strength. These findings open the door for the development of FE-based prediction tools that overcomes the limitations of BMC in accounting for shape and size of the metastatic lesions. Such tools may help clinicians to decide whether a patient needs the prophylactic fixation of an impending fracture.

Assessing the Mechanical Weakness of Vertebrae Affected by Primary Tumors: A Feasibility Study

Patients spend months between the primary spinal tumor diagnosis and the surgical treatment, due to the need for performing chemotherapy and/or radiotherapy. During this period, they are exposed to an unknown risk of fracture. The aim of this study was to assess if it is possible to measure the mechanical strain in vertebrae affected by primary tumors, so as to open the way to an evidence-based scoring or prediction tool. We performed biomechanical tests on three vertebrae with bone tumor removed from patients. The tests were designed so as not to compromise the standard surgical and diagnostic procedures. Non-destructive mechanical tests in combination with state-of-the-art digital image correlation allowed to measure the distribution of strain on the surface of the vertebra. Our study has shown that the strains in the tumor region is circa 3 times higher than in the healthy bones, with principal strain peaks of 40,000/-20,000 microstrain, indicating a stress concentration potentially triggering vertebral fracture. This study has proven it is possible to analyze the mechanical behavior of primary tumor vertebrae as part of the clinical treatment protocol. This will allow building a tool for quantifying the risk of fracture and improving decision making in spine tumors.

Healing pattern classification for thoracolumbar burst fractures after posterior short-segment fixation

BACKGROUND

Thoracolumbar burst fractures can be treated with posterior short-segment fixation. However, no classification can help to estimate whether the healed vertebral body will have sufficient stability after implant removal. We aimed to develop a Healing Pattern Classification (HPC) to evaluate the stability of the healed vertebra based on cavity size and location.

METHODS

Fifty-two thoracolumbar burst fracture patients treated with posterior short-segmental fixation without fusion and followed up for an average of 3.2 years were retrospectively studied. The HPC was divided into 4 types: type I - no cavity; type II - a small cavity with or without the violation of one endplate; type III - a large cavity with or without the violation of one endplate; and type IV - a burst cavity with the violation of both endplates or the lateral cortical shell. The intraobserver and interobserver intraclass correlation coefficients (ICCs) of the HPC were assessed. The demographic characteristics and clinical outcomes of the cohort were compared between the stable group (types I and II) and the unstable group (types III and IV). Logistic regression was conducted to evaluate risk factors for unstable healing.

RESULTS

The intraobserver and interobserver ICCs of the HPC were 0.86 (95% CI = 0.74-0.90) and 0.77 (95% CI = 0.59-0.86), respectively. While the unstable healing group (types III and IV) accounted for 59.6% of the patients, most of these patients were asymptomatic. The preoperative Load Sharing Classification (LSC) comminution score may predict the occurrence of unstable healing (OR = 8.4, 95% CI = 2.4-29.7).

CONCLUSIONS

A reliable classification for assessing the stability of a healed vertebra was developed. With type I and II healing, the vertebra is considered stable, and the implant can be removed. With type III healing, the vertebra may have healing potential, but the implant should not be removed unless type II healing is achieved. With type IV healing, the vertebra is considered extremely unstable, and instrumentation should be maintained. Assessing the LSC comminution score preoperatively may help to predict unstable healing after surgery.