The impact of metastasis on the mineral phase of vertebral bone tissue

Metastatic spine disease alters spinal load-to-strength ratios in patients compared to healthy individuals

Pathologic vertebral fractures (PVF) are common and serious complications in patients with metastatic lesions affecting the spine. Accurate assessment of cancer patients' PVF risk is an unmet clinical need. Load-to-strength ratios (LSRs) evaluated in vivo by estimating vertebral loading from biomechanical modeling and strength from computed tomography imaging (CT) have been associated with osteoporotic vertebral fractures in older adults. Here, for the first time, we investigate LSRs of thoracic and lumbar vertebrae of 135 spine metastases patients compared to LSRs of 246 healthy adults, comparable by age and sex, from the Framingham Heart Study under four loading tasks. Findings include: (1) Osteolytic vertebrae have higher LSRs than osteosclerotic and mixed vertebrae; (2). In patients' vertebrae without CT observed metastases, LSRs were greater than healthy controls. (3) LSRs depend on the spinal region (Thoracic, Thoracolumbar, Lumbar). These findings suggest that LSRs may contribute to identifying patients at risk of incident PVF in metastatic spine disease patients. The lesion-mediated difference suggests that risk thresholds should be established based on spinal region, simulated task, and metastatic lesion type.

Bone mechanical properties were altered in a mouse model of multiple myeloma bone disease

Multiple myeloma bone disease (MMBD) is characterized by the growth of malignant plasma cells in bone marrow, leading to an imbalance in bone (re)modeling favoring excessive resorption. Loss of bone mass and altered microstructure characterize MMBD in humans and preclinical animal models, although, no study to date has examined bone composition or material properties. We hypothesized that MMBD alters bone composition, mineral crystal properties and mechanical properties in the MOPC315.BM.Luc model after intra-tibial injection of myeloma cells and three weeks of daily in vivo tibial loading. Decreased cortical bone elastic modulus and hardness measured by nanoindentation of tibiae were observed in MM-injected mice compared to PBS-injected mice, whereas cortical bone composition, mineral crystal properties measured by Fourier-transform infrared imaging or small angle X-ray scattering, respectively remained unchanged. However, MM-injected mice had thinner cancellous bone mineral particles compared to PBS-injected mice. Mechanical loading did not lead to altered cortical bone composition, mineral structure, or mechanical properties in the context of MM. Unexpectedly, we observed the intra-tibial injection itself altered the material composition of bone, manifested by increased matrix mineralization and crystal size of the hydroxyapatite crystals in the bone matrix. In conclusion, our data suggest that mechanical stimuli can be used as an adjuvant bone anabolic therapy in patients with MMBD to rebuild bone with unaltered composition and mineral structure to reduce subsequent fracture risk.

Dual-energy CT in musculoskeletal imaging: technical considerations and clinical applications

Dual-energy CT stands out as a robust and innovative imaging modality, which has shown impressive advancements and increasing applications in musculoskeletal imaging. It allows to obtain detailed images with novel insights that were once the exclusive prerogative of magnetic resonance imaging. Attenuation data obtained by using different energy spectra enable to provide unique information about tissue characterization in addition to the well-established strengths of CT in the evaluation of bony structures. To understand clearly the potential of this imaging modality, radiologists must be aware of the technical complexity of this imaging tool, the different ways to acquire images and the several algorithms that can be applied in daily clinical practice and for research. Concerning musculoskeletal imaging, dual-energy CT has gained more and more space for evaluating crystal arthropathy, bone marrow edema, and soft tissue structures, including tendons and ligaments. This article aims to analyze and discuss the role of dual-energy CT in musculoskeletal imaging, exploring technical aspects, applications and clinical implications and possible perspectives of this technique.

Treatment affects load to failure and microdamage accumulation in healthy and osteolytic rat vertebrae

Bone turnover and microdamage are impacted by the presence of skeletal metastases which can contribute to increased fracture risk. Treatments for metastatic disease may further impact bone quality. This exploratory study aimed to establish an initial understanding of microdamage accumulation and load to failure in healthy and osteolytic rat vertebrae following focal and systemic cancer treatment (docetaxel (DTX), stereotactic body radiotherapy (SBRT), or zoledronic acid (ZA)). Osteolytic spine metastases were developed in 6-week-old athymic female rats via intracardiac injection of HeLa human cervical cancer cells (day 0). Additional rats served as healthy controls. Rats were either untreated, received SBRT to the T10-L6 vertebrae on day 14 (15 Gy, two fractions), DTX on day 7 or 14, or ZA on day 7. Rats were euthanized on day 21. Tumor burden was assessed with bioluminescence images acquired on day 14 and 21, histology of the excised T11 and L5 vertebrae, and ex-vivo μCT images of the T13-L4. Microstructural parameters (bone volume/total volume, trabecular number, spacing, thickness, and bone mineral density) were measured from L2 vertebrae. Load to failure was measured with axial compressive loading of the L1-L3 motion segments. Microdamage accumulation was labeled in T13 vertebrae with BaSO4 staining and was visualized with high resolution μCT imaging. Microdamage volume fraction was defined as the ratio of BaSO4 to bone volume. DTX administered on day 7 reduced tumor growth significantly (p < 0.05). Microdamage accumulation was found to be increased by the presence of metastases but was reduced by all treatments with ZA showing the largest improvement in HeLa cell injected rats. Load to failure was decreased in untreated and SBRT HeLa cell injected rats compared to healthy controls (p < 0.01). There was a moderate negative correlation between load to failure and microdamage volume fraction in vertebrae from rats injected with HeLa cells (R = -0.35, p = 0.031). Strong correlations were also found between microstructural parameters and load to failure and microdamage accumulation. Several factors, including the presence of osteolytic lesions and use of cancer therapies, influence microdamage accumulation and load to failure in rat vertebrae. Understanding the impact of these treatments on fracture risk of metastatic vertebrae is important to improve management of patients with spinal metastases.

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.

Temporal and spatial changes in bone mineral content and mechanical properties during breast-cancer bone metastases

Cancer cells favour migration and metastasis to bone tissue for 70–80 % of advanced breast cancer patients and it has been proposed that bone tissue provides attractive physical properties that facilitate tumour invasion, resulting in osteolytic and or osteoblastic metastasis. However, it is not yet known how specific bone tissue composition is associated with tumour invasion. In particular, how compositional and nano-mechanical properties of bone tissue evolve during metastasis, and where in the bone they arise, may affect the overall aggressiveness of tumour invasion, but this is not well understood. The objective of this study is to develop an advanced understanding of temporal and spatial changes in nano-mechanical properties and composition of bone tissue during metastasis. Primary mammary tumours were induced by inoculation of immune-competent BALB/c mice with 4T1 breast cancer cells in the mammary fat pad local to the right femur. Microcomputed tomography and nanoindentation were conducted to quantify cortical and trabecular bone matrix mineralisation and nano-mechanical properties. Analysis was performed in proximal and distal femur regions (spatial analysis) of tumour-adjacent (ipsilateral) and contralateral femurs after 3 weeks and 6 weeks of tumour and metastasis development (temporal analysis). By 3 weeks post-inoculation there was no significant difference in bone volume fraction or nano-mechanical properties of bone tissue between the metastatic femora and healthy controls. However, early osteolysis was indicated by trabecular thinning in the distal and proximal trabecular compartment of tumour-bearing femora. Moreover, cortical thickness was significantly increased in the distal region, and the mean mineral density was significantly higher in cortical and trabecular bone tissue in both proximal and distal regions, of ipsilateral (tumour-bearing) femurs compared to healthy controls. By 6 weeks post-inoculation, overt osteolytic lesions were identified in all ipsilateral metastatic femora, but also in two of four contralateral femora of tumour-bearing mice. Bone volume fraction, cortical area, cortical and trabecular thickness were all significantly decreased in metastatic femora (both ipsilateral and contralateral). Trabecular bone tissue stiffness in the proximal femur decreased in the ipsilateral femurs compared to contralateral and control sites. Temporal and spatial analysis of bone nano-mechanical properties and mineralisation during breast cancer invasion reveals changes in bone tissue composition prior to and following overt metastatic osteolysis, local and distant from the primary tumour site. These changes may alter the mechanical environment of both the bone and tumour cells, and thereby play a role in perpetuating the cancer vicious cycle during breast cancer metastasis to bone tissue.

•

Temporal and spatial analyses of bone tissue properties following breast cancer metastasis

•

Trabecular thinning initiated by 3 weeks but overt osteolysis not evident until 6 weeks.

•

Increased bone mineralisation and distal cortical thickness by 3-weeks post-inoculation

Evaluation of Load-To-Strength Ratios in Metastatic Vertebrae and Comparison With Age- and Sex-Matched Healthy Individuals

Vertebrae containing osteolytic and osteosclerotic bone metastases undergo pathologic vertebral fracture (PVF) when the lesioned vertebrae fail to carry daily loads. We hypothesize that task-specific spinal loading patterns amplify the risk of PVF, with a higher degree of risk in osteolytic than in osteosclerotic vertebrae. To test this hypothesis, we obtained clinical CT images of 11 cadaveric spines with bone metastases, estimated the individual vertebral strength from the CT data, and created spine-specific musculoskeletal models from the CT data. We established a musculoskeletal model for each spine to compute vertebral loading for natural standing, natural standing + weights, forward flexion + weights, and lateral bending + weights and derived the individual vertebral load-to-strength ratio (LSR). For each activity, we compared the metastatic spines' predicted LSRs with the normative LSRs generated from a population-based sample of 250 men and women of comparable ages. Bone metastases classification significantly affected the CT-estimated vertebral strength (Kruskal-Wallis, p < 0.0001). Post-test analysis showed that the estimated vertebral strength of osteosclerotic and mixed metastases vertebrae was significantly higher than that of osteolytic vertebrae (p = 0.0016 and p = 0.0003) or vertebrae without radiographic evidence of bone metastasis (p = 0.0010 and p = 0.0003). Compared with the median (50%) LSRs of the normative dataset, osteolytic vertebrae had higher median (50%) LSRs under natural standing (p = 0.0375), natural standing + weights (p = 0.0118), and lateral bending + weights (p = 0.0111). Surprisingly, vertebrae showing minimal radiographic evidence of bone metastasis presented significantly higher median (50%) LSRs under natural standing (p < 0.0001) and lateral bending + weights (p = 0.0009) than the normative dataset. Osteosclerotic vertebrae had lower median (50%) LSRs under natural standing (p < 0.0001), natural standing + weights (p = 0.0005), forward flexion + weights (p < 0.0001), and lateral bending + weights (p = 0.0002), a trend shared by vertebrae with mixed lesions. This study is the first to apply musculoskeletal modeling to estimate individual vertebral loading in pathologic spines and highlights the role of task-specific loading in augmenting PVF risk associated with specific bone metastatic types. Our finding of high LSRs in vertebrae without radiologically observed bone metastasis highlights that patients with metastatic spine disease could be at an increased risk of vertebral fractures even at levels where lesions have not been identified radiologically.

Establishment and Image based evaluation of a New Preclinical Rat Model of Osteoblastic Bone Metastases

Bone remodeling is disrupted in the presence of metastases and can present as osteolytic, osteoblastic or a mixture of the two. Established rat models of osteolytic and mixed metastases have been identified changes in structural and tissue-level properties of bone. The aim of this work was to establish a preclinical rat model of osteoblastic metastases and characterize bone quality changes through image-based evaluation. Female athymic rats (n = 22) were inoculated with human breast cancer cells ZR-75-1 and tumor development tracked over 3-4 months with bioluminescence and in-vivo µCT imaging. Bone tissue-level stereological features were quantified on ex-vivo µCT imaging. Histopathology verified the presence of osteoblastic bone. Bone mineral density distribution was assessed via backscattered electron microscopy. Newly formed osteoblastic bone was associated with reduced mineral content and increased heterogeneity leading to an overall degraded bone quality. Characterizing changes in osteoblastic bone properties is relevant to pre-clinical therapeutic testing and treatment planning.

Toward an artificial intelligence-assisted framework for reconstructing the digital twin of vertebra and predicting its fracture response

This article presents an effort toward building an artificial intelligence (AI) assisted framework, coined ReconGAN, for creating a realistic digital twin of the human vertebra and predicting the risk of vertebral fracture (VF). ReconGAN consists of a deep convolutional generative adversarial network (DCGAN), image-processing steps, and finite element (FE) based shape optimization to reconstruct the vertebra model. This DCGAN model is trained using a set of quantitative micro-computed tomography (micro-QCT) images of the trabecular bone obtained from cadaveric samples. The quality of synthetic trabecular models generated using DCGAN are verified by comparing a set of its statistical microstructural descriptors with those of the imaging data. The synthesized trabecular microstructure is then infused into the vertebra cortical shell extracted from the patient's diagnostic CT scans using an FE-based shape optimization approach to achieve a smooth transition between trabecular to cortical regions. The final geometrical model of the vertebra is converted into a high-fidelity FE model to simulate the VF response using a continuum damage model under compression and flexion loading conditions. A feasibility study is presented to demonstrate the applicability of digital twins generated using this AI-assisted framework to predict the risk of VF in a cancer patient with spinal metastasis.

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).

Assessment of the efficiency of different chemical treatments and ultrasonic cleaning for defatting of cancellous bone samples

Our study aimed to asses the defatting efficiency of different methods, which are commonly used and easily available in the laboratory in order to find a method that is effective, convenient, safe, and economical. Cylindrical cancellous bone specimens were obtained from fresh-frozen human cadaver femoral condyles, cut into multiple small specimens (Ø8 × 2 mm), and assigned to two groups that were treated with either chemical solvent soaking (Solvent group) or ultrasonic cleaning (Ultrasound group). Each group was divided into several subgroups based on different treatments. Digital photographs were taken of each specimen. The difference of material density (Δρ_b), apparent density (Δρ_app), and porosity (ΔP) before and after treatment were used as evaluation indicators. For the solvent group, in Δρ_b, only the combination of 99% ethanol and detergent solution showed a significant difference before and after treatment (P = 0.00). There was no significant difference in ΔP among acetone, the mixture of 99% ethanol and acetone, and the combination of 99% ethanol and detergent solution (P = 0.93). For the ultrasound group, the median of all subgroups in Δρ_app and ΔP were all lower than the solvent group. The combination of 99% ethanol and detergent solution (v/v = 1:20), as well as the mixture of 99% ethanol and acetone (v/v = 1:1), seem to be the optimal defatting methods for 2 mm thick cancellous bone slices due to their effectiveness, availability, low-cost and safety. Chemical soaking for 24 h is more effective than ultrasonic cleaning with 99% ethanol or acetone for 20 or 40 min.

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.

Impact of radiofrequency ablation (RFA) on bone quality in a murine model of bone metastases

Thermal therapies such as radiofrequency ablation (RFA) are gaining widespread clinical adoption in the local treatment of skeletal metastases. RFA has been shown to successfully destroy tumor cells, yet the impact of RFA on the quality of the surrounding bone has not been well characterized. RFA treatment was performed on femora of rats with bone metastases (osteolytic and osteoblastic) and healthy age matched rats. Histopathology, second harmonic generation imaging and backscatter electron imaging were used to characterize changes in the structure, organic and mineral components of the bone after RFA. RFA treatment was shown to be effective in targeting tumor cells and promoting subsequent new bone formation without impacting the surrounding bone negatively. Mineralization profiles of metastatic models were significantly improved post-RFA treatment with respect to mineral content and homogeneity, suggesting a positive impact of RFA treatment on the quality of cancer involved bone. Evaluating the impact of RFA on bone quality is important in directing the growth of this minimally invasive therapeutic approach with respect to fracture risk assessment, patient selection, and multimodal treatment planning.

Large Lytic Defects Produce Kinematic Instability and Loss of Compressive Strength in Human Spines: An in Vitro Study

BACKGROUND

In patients with spinal metastases, kinematic instability is postulated to be a predictor of pathologic vertebral fractures. However, the relationship between this kinematic instability and the loss of spinal strength remains unknown.

METHODS

Twenty-four 3-level thoracic and lumbar segments from 8 cadaver spines from female donors aged 47 to 69 years were kinematically assessed in axial compression (180 N) and axial compression with a flexion or extension moment (7.5 Nm). Two patterns of lytic defects were mechanically simulated: (1) a vertebral body defect, corresponding to Taneichi model C (n = 13); and (2) the model-C defect plus destruction of the ipsilateral pedicle and facet joint, corresponding to Taneichi model E (n = 11). The kinematic response was retested, and compression strength was measured. Two-way repeated-measures analysis of variance was used to test the effect of each model on the kinematic response of the segment. Multivariable linear regression was used to test the association between the kinematic parameters and compressive strength of the segment.

RESULTS

Under a flexion moment, and for both models C and E, the lesioned spines exhibited greater flexion range of motion (ROM) and axial translation than the control spines. Both models C and E caused lower extension ROM and greater axial, sagittal, and transverse translation under an extension moment compared with the control spines. Two-way repeated-measures analysis revealed that model E, compared with model C, caused significantly greater changes in extension and torsional ROM under an extension moment, and greater sagittal translation under a flexion moment. For both models C and E, greater differences in flexion ROM and sagittal translation under a flexion moment, and greater differences in extension ROM and in axial and transverse translation under an extension moment, were associated with lower compressive strength of the lesioned spines.

CONCLUSIONS

Critical spinal lytic defects result in kinematic abnormalities and lower the compressive strength of the spine.

CLINICAL RELEVANCE

This experimental study demonstrates that lytic foci degrade the kinematic stability and compressive strength of the spine. Understanding the mechanisms for this degradation will help to guide treatment decisions that address inferred instability and fracture risk in patients with metastatic spinal disease.

High resolution 3D structures of mineralized tissues in health and disease

A thorough knowledge of the structures of healthy mineralized tissues, such as bone or cartilage, is key to understanding the pathological changes occurring during disease. Such knowledge enables the underlying mechanisms that are responsible for pathology to be pinpointed. One high-resolution 3D method in particular - focused ion beam-scanning electron microscopy (FIB-SEM) - has fundamentally changed our understanding of healthy vertebrate mineralized tissues. FIB-SEM can be used to study demineralized matrix, the hydrated components of tissue (including cells) using cryo-fixation and even untreated mineralized tissue. The latter requires minimal sample preparation, making it possible to study enough samples to carry out studies capable of detecting statistically significant differences - a pre-requisite for the study of pathological tissues. Here, we present an imaging and characterization strategy for tissue structures at different length scales, describe new insights obtained on healthy mineralized tissues using FIB-SEM, and suggest future research directions for both healthy and diseased mineralized tissues.

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.

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.

CT based quantitative measures of the stability of fractured metastatically involved vertebrae treated with spine stereotactic body radiotherapy

Mechanical instability secondary to vertebral metastases can lead to pathologic vertebral compression fracture (VCF) mechanical pain, neurological compromise, and the need for surgical stabilization. Stereotactic body radiation therapy (SBRT) as a treatment for spinal metastases is effective for pain and local tumor control, it has been associated with an increased risk of VCF. This study quantified computed tomography (CT) based stability measures in metastatic vertebrae with VCF treated with spine SBRT. It was hypothesized that semi-automated quantification of VCF based on CT metrics would be related to clinical outcomes. 128 SBRT treated spinal metastases patients were identified from a prospective database. Of these, 18 vertebral segments were identified with a VCF post-SBRT. A semi-automated system for quantifying VCF was developed based on CT imaging before and after SBRT. The system identified and segmented SBRT treated vertebral bodies, calculated stability metrics at single time points and changes over time. In the vertebrae that developed a new (n = 7) or progressive (n = 11) VCF following SBRT, the median time to VCF/VCF progression was 1.74 months (range 0.53-7.79 months). Fractured thoracolumbar vertebrae that went on to be stabilized (cemented and/or instrumented), had greater fractured vertebral body volume progression over time (12%) compared to those not stabilized (0.4%, p < 0.05). Neither the spinal instability neoplastic score (SINS) or any single timepoint stability metrics in post-hoc analyses correlated with future stabilization. This pilot study presents a quantitative semi-automated method assessing fractured thoracolumbar vertebrae based on CT. Increased fractured vertebral body volume progression post-SBRT was shown to predict those patients who were subsequently stabilized, motivating study of methods that assess temporal radiological changes toward augmenting existing clinical management in the metastatic spine.

The effects of metastatic lesion on the structural determinants of bone: Current clinical and experimental approaches

Metastatic bone disease is incurable with an associated increase in skeletal-related events, particularly a 17-50% risk of pathologic fractures. Current surgical and oncological treatments are palliative, do not reduce overall mortality, and therefore optimal management of adults at risk of pathologic fractures presents an unmet medical need. Plain radiography lacks specificity and may result in unnecessary prophylactic fixation. Radionuclide imaging techniques primarily supply information on the metabolic activity of the tumor or the bone itself. Magnetic resonance imaging and computed tomography provide excellent anatomical and structural information but do not quantitatively assess bone matrix. Research has now shifted to developing unbiased data-driven tools that can predict risk of impending fractures and guide individualized treatment decisions. This review discusses the state-of-the-art in clinical and experimental approaches for prediction of pathologic fractures with bone metastases. Alterations in bone matrix quality are associated with an age-related increase in skeletal fragility but the impact of metastases on the intrinsic material properties of bone is unclear. Engineering-based analyses are non-invasive with the capability to evaluate oncological treatments and predict failure due to the progression of metastasis. The combination of these approaches may improve our understanding of the underlying deterioration in mechanical performance.

Diagnosis and Treatment of Bone Metastases in Breast Cancer: Radiotherapy, Local Approach and Systemic Therapy in a Guide for Clinicians

The standard care for metastatic breast cancer (MBC) is systemic therapies with imbrication of focal treatment for symptoms. Recently, thanks to implementation of radiological and metabolic exams and development of new target therapies, oligometastatic and oligoprogressive settings are even more common-paving the way to a paradigm change of focal treatments role. In fact, according to immunophenotype, radiotherapy can be considered with radical intent in these settings of patients. The aim of this literature review is to analyze available clinical data on prognosis of bone metastases from breast cancer and benefits of available treatments for developing a practical guide for clinicians.

Mechanical behavior of metastatic femurs through patient-specific computational models accounting for bone-metastasis interaction

This paper proposes a computational model based on a finite-element formulation for describing the mechanical behavior of femurs affected by metastatic lesions. A novel geometric/constitutive description is introduced by modelling healthy bone and metastases via a linearly poroelastic constitutive approach. A Gaussian-shaped graded transition of material properties between healthy and metastatic tissues is prescribed, in order to account for the bone-metastasis interaction. Loading-induced failure processes are simulated by implementing a progressive damage procedure, formulated via a quasi-static displacement-driven incremental approach, and considering both a stress- and a strain-based failure criterion. By addressing a real clinical case, left and right patient-specific femur models are geometrically reconstructed via an ad-hoc imaging procedure and embedding multiple distributions of metastatic lesions along femurs. Significant differences in fracture loads, fracture mechanisms, and damage patterns, are highlighted by comparing the proposed constitutive description with a purely elastic formulation, where the metastasis is treated as a pseudo-healthy tissue or as a void region. Proposed constitutive description allows to capture stress/strain localization mechanisms within the metastatic tissue, revealing the model capability in describing possible strain-induced mechano-biological stimuli driving onset and evolution of the lesion. The proposed approach opens towards the definition of effective computational strategies for supporting clinical decision and treatments regarding metastatic femurs, contributing also to overcome some limitations of actual standards and procedures.

Elevated Microdamage Spatially Correlates with Stress in Metastatic Vertebrae

Metastasis of cancer to the spine impacts bone quality. This study aims to characterize vertebral microdamage secondary to metastatic disease considering the pattern of damage and its relationship to stress and strain under load. Osteolytic and mixed osteolytic/osteoblastic vertebral metastases were produced in athymic rats via HeLa cervical or canine Ace-1 prostate cancer cell inoculation, respectively. After 21 days, excised motion segments (T12-L2) were µCT scanned, stained with BaSO₄ and re-imaged. T13-L2 motion segments were loaded in axial compression to induce microdamage, re-stained and re-imaged. L1 (loaded) and T12 (unloaded) vertebrae were fixed, sample blocks cut, polished and BSE imaged. µFE models were generated of all L1 vertebrae with displacement boundary conditions applied based on the loaded µCT images. µCT stereological analysis, BSE analysis and µFE derived von Mises stress and principal strains were quantitatively compared (ANOVA), spatial correlations determined and patterns of microdamage assessed qualitatively. BaSO₄ identified microdamage was found to be spatially correlated with regions of high stress in µFEA. Load-induced microdamage was shown to be elevated in the presence of osteolytic and mixed metastatic disease, with diffuse, crossed hatched areas of microdamage present in addition to linear microdamage and microfractures in metastatic tissue, suggesting diminished bone quality.

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.

The Size of Simulated Lytic Metastases Affects the Strain Distribution on the Anterior Surface of the Vertebra

Metastatic lesions of the vertebra are associated with risk of fracture, which can be disabling and life-threatening. In the literature, attempts are found to identify the parameters that reduce the strength of a metastatic vertebra leading to spine instability. However, a number of controversial issues remain. Our aim was to quantify how the strain distribution in the vertebral body is affected by the presence and by the size of a simulated metastatic defect. Five cadaveric thoracic spine segments were subjected to non-destructive presso-flexion while intact, and after simulation of metastases of increasing size. For the largest defect, the specimens were eventually tested to failure. The full-field strain distribution in the elastic range was measured with digital image correlation (DIC) on the anterior surface of the vertebral body. The mean strain in the vertebra remained similar to the intact when the defects were smaller than 30% of the vertebral volume. The mean strains became significantly larger than in the intact for larger defects. The map of strain and its statistical distribution indicated a rather uniform condition in the intact vertebra and with defects smaller than 30%. Conversely, the strain distribution became significantly different from the intact for defects larger than 30%. A strain peak appeared in the region of the simulated metastasis, where fracture initiated during the final destructive test. This is a first step in understanding how the features of metastasis influence the vertebral strain and for the construction of a mechanistic model to predicted fracture.

Multiscale characterization of the mineral phase at skeletal sites of breast cancer metastasis

Skeletal metastases, the leading cause of death in advanced breast cancer patients, depend on tumor cell interactions with the mineralized bone extracellular matrix. Bone mineral is largely composed of hydroxyapatite (HA) nanocrystals with physicochemical properties that vary significantly by anatomical location, age, and pathology. However, it remains unclear whether bone regions typically targeted by metastatic breast cancer feature distinct HA materials properties. Here we combined high-resolution X-ray scattering analysis with large-area Raman imaging, backscattered electron microscopy, histopathology, and microcomputed tomography to characterize HA in mouse models of advanced breast cancer in relevant skeletal locations. The proximal tibial metaphysis served as a common metastatic site in our studies; we identified that in disease-free bones this skeletal region contained smaller and less-oriented HA nanocrystals relative to ones that constitute the diaphysis. We further observed that osteolytic bone metastasis led to a decrease in HA nanocrystal size and perfection in remnant metaphyseal trabecular bone. Interestingly, in a model of localized breast cancer, metaphyseal HA nanocrystals were also smaller and less perfect than in corresponding bone in disease-free controls. Collectively, these results suggest that skeletal sites prone to tumor cell dissemination contain less-mature HA (i.e., smaller, less-perfect, and less-oriented crystals) and that primary tumors can further increase HA immaturity even before secondary tumor formation, mimicking alterations present during tibial metastasis. Engineered tumor models recapitulating these spatiotemporal dynamics will permit assessing the functional relevance of the detected changes to the progression and treatment of breast cancer bone metastasis.

Collagen fibril organization within rat vertebral bone modified with metastatic involvement

Metastatic involvement diminishes the mechanical integrity of vertebral bone, however its specific impact on the structural characteristics of a primary constituent of bone tissue, the collagen-I fibril matrix, has not been adequately characterized. Female athymic rats were inoculated with HeLa or Ace-1 cancer cells lines producing osteolytic or mixed (osteolytic & osteoblastic) metastases respectively. A maximum of 21days was allowed between inoculation and rat sacrifice for vertebrae extraction. Linear polarization-in, polarization-out (PIPO) second harmonic generation (SHG) and transmission electron microscopy (TEM) imaging was utilized to assess the impact of metastatic involvement on collagen fibril organization. Increased observations of deviations in the typical plywood motif or a parallel packing structure and an increased average measured susceptibility ratio (related to relative degree of in-plane vs. out-plane fibrils in the analyzed tissue area) in bone adjacent to metastatic involvement was indicative of change in fibrilar organization compared to healthy controls. In particular, collagen-I fibrils in tumour-induced osteoblastic bone growth showed no adherence to the plywood motif or parallel packing structure seen in healthy lamellar bone, exhibiting a much higher susceptibility ratio and degree of fibril disorder. Negative correlations were established between measured susceptibility ratios and the hardness and modulus of metastatic bone tissue assessed in a previous study. Characterizing modifications in tissue level properties is key in defining bone quality in the presence of metastatic disease and their potential impact on material behaviour.