Noninvasive imaging predicts failure load of the spine with simulated osteolytic defects

Computed tomography-based structural rigidity analysis can assess tumor- and treatment-induced changes in rat bones with metastatic lesions

BACKGROUND

Breast cancer (BrCa) is a predominant malignancy, with metastasis occurring in one in eight patients, nearly half of which target the bone, leading to serious complications such as pain, fractures, and compromised mobility. Structural rigidity, crucial for bone strength, becomes compromised with osteolytic lesions, highlighting the vulnerability and increased fracture risk in affected areas. Historically, two-dimensional radiographs have been employed to predict these fracture risks; however, their limitations in capturing the three-dimensional structural and material changes in bone have raised concerns. Recent advances in CT-based Structural Rigidity Analysis (CTRA), offer a promising, more accurate non-invasive 3D approach. This study aims to assess the efficacy of CTRA in monitoring osteolytic lesions' progression and response to therapy, suggesting its potential superiority over existing methodologies in guiding treatment strategies.

METHODS

Twenty-seven female nude rats underwent femoral intra-medullary inoculation with MDA-MB-231 human breast cancer cells or saline control. They were divided into Control, Cancer Control, Ibandronate, and Paclitaxel groups. Osteolytic progression was monitored weekly using biplanar radiography, quantitative computed tomography (QCT), and dual-energy X-ray absorptiometry (DEXA). CTRA was employed to predict fracture risk, normalized using the contralateral femur. Statistical analyses, including Kruskal-Wallis and ANOVA, assessed differences in outcomes among groups and over time.

RESULTS

Biplanar radiographs showed treatment benefits over time; however, only certain time-specific differences between the Control and other treatment groups were discernible. Notably, observer subjectivity in X-ray scoring became evident, with significant inter-operator variations. DEXA measurements for metaphyseal Bone Mineral Content (BMC) did not exhibit notable differences between groups. Although diaphyseal BMC highlighted some variance, it did not reveal significant differences between treatments at specific time points, suggesting a limited ability for DEXA to differentiate between treatment effects. In contrast, the CTRA consistently demonstrated variations across different treatments, effectively capturing bone rigidity changes over time, and the axial- (EA), bending- (EI), and torsional rigidity (GJ) outcomes from the CTRA method successfully distinguished differences among treatments at specific time points.

CONCLUSION

Traditional approaches, such as biplanar radiographs and DEXA, have exhibited inherent limitations, notably observer bias and time-specific inefficacies. Our study accentuates the capability of CTRA in capturing real-time, progressive changes in bone structure, with the potential to predict fractures more accurately and provide a more objective analysis. Ultimately, this innovative approach may bridge the existing gap in clinical guidelines, ushering in enhanced Clinical Decision Support Tool (CDST) for both surgical and non-surgical treatments.

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.

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.

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.

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.

Fracture Risk of Long Bone Metastases: A Review of Current and New Decision-Making Tools for Prophylactic Surgery

Simple Summary

Long bone metastases are frequently a pivotal point in the oncological history of patients. Weakening of the bone results in pathologic fractures that not only compromise patient function but also their survival. Therefore, the main issue for tumor boards remains timely assessment of the risk of fracture, as this is a key consideration in providing preventive surgery while also avoiding overtreatment. As the Mirels scoring system takes into account both the radiological and the clinical criteria, it has been used worldwide since the 1990s. However, due to increasing concern regarding the lack of accuracy, new thresholds have been defined for the identification of impending fractures that require prophylactic surgery, on the basis of axial cortical involvement and biomechanical models involving quantitative computed tomography. The aim of this review is to establish a state-of-the-art of the risk assessment of long bone metastases fractures, from simple radiologic scores to more complex multidimensional bone models, in order to define new decision-making tools.

Abstract

Long bone pathological fractures very much reflect bone metastases morbidity in many types of cancer. Bearing in mind that they not only compromise patient function but also survival, identifying impending fractures before the actual event is one of the main concerns for tumor boards. Indeed, timely prophylactic surgery has been demonstrated to increase patient quality of life as well as survival. However, early surgery for long bone metastases remains controversial as the current fracture risk assessment tools lack accuracy. This review first focuses on the gold standard Mirels rating system. It then explores other unique imaging thresholds such as axial or circumferential cortical involvement and the merits of nuclear imaging tools. To overcome the lack of specificity, other fracture prediction strategies have focused on biomechanical models based on quantitative computed tomography (CT): computed tomography rigidity analysis (CT-RA) and finite element analysis (CT-FEA). Despite their higher specificities in impending fracture assessment, their limited availability, along with a need for standardization, have limited their use in everyday practice. Currently, the prediction of long bone pathologic fractures is a multifactorial process. In this regard, machine learning could potentially be of value by taking into account clinical survival prediction as well as clinical and improved CT-RA/FEA data.

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.

CT-based structural analyses of vertebral fractures with polymeric augmentation: A study of cadaveric three-level spine segments

Pathologic vertebral fractures due to metastasis can occur under normal physiologic activities, leading to pain and neurologic deficit. Prophylactic vertebroplasty is a technique used to augment vertebral strength and reduce the risk of fracture. Currently, no technique is available to objectively assess vertebral fracture risk in metastatically-involved vertebral bodies. The aim of the current study was to develop an image-based computational technique to estimate fracture force outcomes during bending. To this end, mechanical testing was performed on intact, simulated defect, PMMA-augmented, and PPF-augmented 3-level spine segments from both sexes under a compression/flexion-type loading condition. The augmentation performance of poly(methyl methacrylate) (PMMA) and poly(propylene fumarate) (PPF) were also evaluated and compared. Cylindrical defects were created in 3-level spine segments with attached posterior elements and ligaments. Using CT images of each segment, a rigidity analysis technique was developed and used for predicting fracture forces during bending. On average, PPF strengthened the segments by about 630 N, resulting in fracture forces similar to those observed in the intact and PMMA-augmented groups. Female spines fractured at about 1150 N smaller force than did male spines. Rigidity analysis, along with age, explained 66% variability in experimental outcomes. This number increased to 74% when vertebral size and age were added to the rigidity analysis as explanatory variables. Both PPF and PMMA similarly increased fracture strength to the level of intact specimens. The results suggest that PPF can be a suitable candidate for augmentation purposes and rigidity analysis can be a promising predicting tool for vertebral fracture forces.

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.

Density-Dependent Material and Failure Criteria Equations Highly Affect the Accuracy and Precision of QCT/FEA-Based Predictions of Osteoporotic Vertebral Fracture Properties

About 700,000 vertebral fractures occur in the US as a result of bone loss. Quantitative computed tomography (QCT)-based finite element analysis (FEA) is a promising tool for fracture risk prediction that is becoming attractive in the clinical setting. The goals of this study were (1) to perform individual and pooled specimen optimization using inverse QCT/FEA modeling to obtain ash density-elastic modulus equations incorporating the whole vertebral body and accounting for all variables used during FE modeling, and (2) to determine the effect of material equations and failure criteria on the accuracy and precision of mechanical properties. Fifty-four (54) human vertebrae were used to optimize material equations based on experimental outcomes and, together with a previously proposed material equation, were implemented in our models using three different failure criteria to obtain fracture loads. Our robust QCT/FEA approach predicted 78% of the failure loads. Material equations resulted in poor accuracy in the predicted stiffness, yet yielded good precision and, more importantly, strong correlations with fracture loads. Both material and fracture criterion equations are equally important in estimating accurate and precise QCT/FEA predictions. Results suggest that both elastic modulus and fracture criterion equations should be validated against experimental outcomes to better explain the response of the tissue under various conditions.

Simulated lesions representative of metastatic disease predict proximal femur failure strength more accurately than idealized lesions

Metastatic disease in bone is characterized by highly amorphous and variable lesion geometry, with increased fracture risk. Assumptions of idealized lesion geometry made in previous finite element (FE) studies of metastatic disease in the proximal femur may not sufficiently capture effects of local stress/strain concentrations on predicted failure strength. The goal of this study was to develop and validate a FE failure model of the proximal femur incorporating artificial defects representative of physiologic metastatic disease. Data from 11 cadaveric femur specimens were randomly divided into either a training set (n = 5) or a test set (n = 6). Clinically representative artificial defects were created, and the femurs were loaded to failure under offset torsion. Voxel-based FE models replicating the experimental setup were created from the training set pre-fracture computed tomography data. Failure loads from the linear model with maximum principal strain failure criterion correlated best with the experimental data (R² = 0.86, p = 0.024). The developed model was found to be reliable when applied to the test dataset with a relatively low RMSE of 46.9 N, mean absolute percent error of 12.7 ± 17.1%, and cross-validation R² = 0.88 (p < 0.001). Models simulating realistic lesion geometry explained an additional 26% of the variance in experimental failure load compared to idealized lesion models (R² = 0.62, p = 0.062). Our validated automated FE model representative of physiologic metastatic disease may improve clinical fracture risk prediction and facilitate research studies of fracture risk during functional activities and with treatment interventions.

Mechanical testing setups affect spine segment fracture outcomes

The purpose of the work presented here was to establish an experimental testing configuration that would generate a bending compression fracture in a laboratory setting. To this end, we designed and fabricated a fixture to accommodate a three level spine segment and to be able to perform mechanical testing by applying an off-centric compressive loading to create a flexion-type motion. Forces and moments occurring during testing were measured with a six-channel load cell. The initial testing configuration (Fixture A) included plates connected to the superior potted vertebral body and to the ball-socket joint of the testing system ram. Surprisingly, while all cadaveric specimens underwent a similar off-centric compressive loading, most of the specimens showed extension outcomes as opposed to the intended pure-flexion motion. The extension was due to fixture size and weight; by applying an off-centric load directly on the top plate, unintended large shear forces were generated. To resolve the issue, several modifications were made to the original fixture configuration. These modifications included the removal of the superior plates and the implementation of wedges at the superior surface of the fixture (Fixture B). A synthetic sample was used during this modification phase to minimize the number of cadaveric specimens while optimizing the process. The best outcomes were consistently observed when a 15°-wedge was used to provide flexion-type loading. Cadaveric specimens were then experimentally tested to fracture using the modified testing configuration (Fixture B). A comparison between both fixtures, A and B, revealed that almost all biomechanical parameters, including force, moment, and displacement data, were affected by the testing setup. These results suggest that fixture design and implementation for testing is of extreme importance, and can influence the fracture properties and affect the intended motion.

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 application of finite element modelling based on clinical pQCT for classification of fracture status

Fracture risk assessment using dual-energy X-ray absorptiometry (DXA) frequently fails to diagnose osteoporosis amongst individuals who later experience fragility fractures. Hence, more reliable techniques that improve the prediction of fracture risk are needed. In this study, we evaluated a finite element (FE) modelling framework based on clinical peripheral quantitative computed tomography (pQCT) imaging of the tibial epiphysis and diaphysis to predict the stiffness at these locations in compression, shear, torsion and bending. The ability of these properties to identify a group of women who had recently sustained a low-trauma fracture from an age- and weight-matched control group was determined and compared to clinical pQCT and DXA properties and structural properties based on composite beam theory. The predicted stiffnesses derived from the FE models and composite beam theory were significantly different (p < 0.05) between the control and fracture groups, whereas no meaningful differences were observed using DXA and for the stress-strain indices (SSIs) derived using pQCT. The diagnostic performance of each property was assessed by the odds ratio (OR) and the area under the receiver operating curve (AUC), and both were greatest for the FE-predicted shear stiffness (OR 16.09, 95% CI 2.52-102.56, p = 0.003) (AUC: 0.80, 95% CI 0.67-0.93). The clinical pQCT variable total density (ρ_tot) and a number of structural and FE-predicted variables had a similar probability of correct classification between the control and fracture groups (i.e. ORs and AUCs with mean values greater than 5.00 and 0.80, respectively). In general, the diagnostic characteristics were lower for variables derived using DXA and for the SSIs (i.e. ORs and AUCs with mean values of 1.65-2.98 and 0.64-0.71, respectively). For all properties considered, the trabecular-dominant tibial epiphysis exhibited enhanced classification characteristics, as compared to the cortical-dominant tibial diaphysis. The results of this study demonstrate that bone properties may be derived using FE modelling that have the potential to enhance fracture risk assessment using conventional pQCT or DXA instruments in clinical settings.

Assessment of vertebral wedge strength using cancellous textural properties derived from digital tomosynthesis and density properties from dual energy X-ray absorptiometry and high resolution computed tomography

The purpose of this study was to examine the potential of digital tomosynthesis (DTS) derived cancellous bone textural measures to predict vertebral strength under conditions simulating a wedge fracture. 40 vertebral bodies (T6, T8, T11, and L3 levels) from 5 male and 5 female cadaveric donors were utilized. The specimens were scanned using dual energy X-ray absorptiometry (DXA) and high resolution computed tomography (HRCT) to obtain measures of bone mineral density (BMD) and content (BMC), and DTS to obtain measures of bone texture. Using a custom loading apparatus designed to deliver a nonuniform displacement resulting in a wedge deformity similar to those observed clinically, the specimens were loaded to fracture and their fracture strength was recorded. Mixed model regressions were used to determine the associations between wedge strength and DTS derived textural variables, alone and in the presence of BMD or BMC information. DTS derived fractal, lacunarity and mean intercept length variables correlated with wedge strength, and individually explained up to 53% variability. DTS derived textural variables, notably fractal dimension and lacunarity, contributed to multiple regression models of wedge strength independently from BMC and BMD. The model from a scan orientation transverse to the spine axis and in the anterior-posterior view resulted in highest explanatory capability (R²_adj = 0.91), with a scan orientation parallel to the spine axis and in the lateral view offering an alternative (R²_adj = 0.88). In conclusion, DTS can be used to examine cancellous texture relevant to vertebral wedge strength, and potentially complement BMD in assessment of vertebral fracture risk.

Female Human Spines with Simulated Osteolytic Defects: CT-based Structural Analysis of Vertebral Body Strength

Purpose To evaluate a CT structural analysis protocol (SAP) for estimating the strength of human female cadaveric spines with lytic lesions. Materials and Methods Osteolytic foci was created in the middle vertebra of 44 thoracic and lumbar three-level segments from 11 female cadavers (age range, 50-70 years). The segments underwent CT by using standard clinical protocol and their failure strength was assessed at CT SAP. The spines were mechanically tested to failure in pure axial compression or in compression with torsion. The relationships of defect size, bone mineral density, and predicted failure load (at CT SAP) with measured vertebral strength were assessed with linear regression. Analysis of variance and Tukey test were used to evaluate the effect of region and mechanical test on spine strength. Results With axial compression, CT SAP predictions of vertebral strength correlated with the thoracic (r = 0.84; P < .001) and lumbar (r = 0.85; P < .001) segment-measured strength. Bone mineral density correlated with the lumbar (r = 0.64; P = .003) and thoracic (r, 0.51; P = .050) strength. At compression with torsion, CT SAP predictions of strength were moderately correlated with vertebral strength (r = 0.66; P = .018). At compression with torsion, bone mineral density was not correlated with spinal strength (thoracic and lumbar: r = 0.31 and r = 0.26, respectively; P = .539 and .610, respectively). The lytic focus size (range, 28%-41%) was not associated with vertebral strength. Conclusion CT SAP assessment of strength in vertebrae with lytic lesions correlated with the measured strength of female vertebral bodies. ^© RSNA, 2018 Online supplemental material is available for this article.

The trabecular effect: A population-based longitudinal study on age and sex differences in bone mineral density and vertebral load bearing capacity

BACKGROUND

Approximately 16-24% of postmenopausal women are affected by vertebral fractures, negatively affecting their quality of life. Trabecular and cortical bones in vertebrae decline differently with age, thus having a distinct impact on vertebral failure loads. The purpose of this study was to investigate the effect of trabecular and cortical volumetric bone mineral density loss over time on estimated failure loads; and to evaluate the effect of sex and age.

METHOD

Fracture properties from a cohort of 82 patients were evaluated for L₁-L₃ vertebrae at baseline and 6th year using an image-based method that implements axial rigidity analysis. Cortical and trabecular volumetric bone mineral density were obtained, as well as their individual contribution to total failure load. Regression analyses were performed to determine the effect of age and sex on volumetric bone mineral density and failure loads.

FINDINGS

Decline in trabecular and cortical volumetric bone mineral density, and failure load was sex-dependent (p ≤ 0.0095). Cortical and trabecular volumetric bone mineral density reduced 2.08 (g/cm³)/year and 2.02 (g/cm³)/year, respectively. A 1012 N difference in failure load, ~70% attributed to trabecular bone, was found between men and women of similar age. Over 6 years, this difference increased by 287 N. Areal bone mineral density measured by dual X-ray absorptiometry explained ~60% of the vertebral failure load.

INTERPRETATION

Trabecular bone has a significantly greater effect than cortical bone on the structural integrity and load bearing capacity of vertebrae. This might lead to a higher incidence of fragility fractures in osteoporotic women. Our non-invasive, quantitative computed tomography image-based approach may improve prevention, monitoring, and management of fractures.

The Role of the Size and Location of the Tumors and of the Vertebral Anatomy in Determining the Structural Stability of the Metastatically Involved Spine: a Finite Element Study

Vertebral fractures associated with the loss of structural integrity of neoplastic vertebrae are common, and determined to the deterioration of the bone quality in the lesion area. The prediction of the fracture risk in metastatically involved spines can guide in deciding if preventive solutions, such as medical prophylaxis, bracing, or surgery are indicated for the patient. In this study, finite element models of 22 thoracolumbar vertebrae were built based on CT scans of three spines, covering a wide spectrum of possible clinical scenarios in terms of age, bone quality and degenerative features, taking into account the local material properties of bone tissue. Simulations were performed in order to investigate the effect of the size and location of the tumoral lesion, the bone quality and the vertebral level in determining the structural stability of the neoplastic vertebrae. Tumors with random size and positions were added to the models, for a total of 660 simulations in which a compressive load was simulated. Results highlighted the fundamental role of the tumor size, whereas the other parameters had a lower, but non-negligible impact on the axial collapse of the vertebra, the vertebral bulge in the transverse plane and the canal narrowing under the application of the load. All the considered parameters are radiologically measurable, and can therefore be translated in a straightforward way to the clinical practice to support decisions about preventive treatment of metastatic fractures.

Predicting pathological fracture in femoral metastases using a clinical CT scan based algorithm: A case-control study

BACKGROUND

We assessed whether there was a difference in attenuation measurements (in Hounsfield units - HU) and geometric distribution of HU between femora with metastatic lesions that fracture, and metastatic lesions that did not fracture nor underwent prophylactic fixation.

METHODS

Nine patients with femoral metastases who underwent CT and developed a pathological fracture were matched to controls. All femora were delineated in axial CT slices using a region of interest (ROI) tool; the HU within these ROIs were used to calculate: (1) the cumulative HU of the affected over the nonaffected side per slice and presented as a percentage, and (2) the cumulative HU accounting for geometric distribution (polar moment of HU). We repeated the analyses including cortical bone only (HU of 600 and above).

RESULTS

CT-based calculations did not differ between patients with a lesion that fractured and those that did not fracture nor underwent prophylactic fixation when analyzing all tissue. However, when including cortical bone only, the pathological fracture group had a lower cumulative HU value compared to the no fracture and no fixation group for the weakest cross-sectional CT image (pathological fracture group, mean: 71, SD: 23 and no fracture and no prophylactic fixation group, mean: 85, SD: 18, p = 0.042) and the complete lesion analysis (pathological fracture group, mean: 78, SD: 21 and no fracture and no prophylactic fixation group, mean: 92, SD: 15, p = 0.032).

CONCLUSION

The demonstrated CT-based algorithms can be useful for predicting pathological fractures in metastatic lesions.

A pediatric animal model to evaluate the effects of disuse on musculoskeletal growth and development

Prolonged immobilization in hospitalized children can lead to fragility fractures and muscle contractures and atrophy. The purpose of this study was to develop a lower-extremity disuse rabbit model with musculoskeletal changes similar to those observed in children subjected to prolonged immobilization. Six-week-old rabbits were randomly assigned to control (CTRL, n=4) or bilateral sciatic and femoral neurectomy (bSFN, n=4) groups. Trans-axial helical CT scans of each rabbit׳s hind limbs were acquired after eight weeks. The rabbits were then euthanized and the tibiae and calcanea were harvested from each rabbit. μCT imaging was performed on the tibiae and calcanea mid-diaphysis. Four-point bending, gas pycnometry, and ashing were then performed on each tibia. All comparisons reflect the differences between the bSFN and CTRL rabbits. Significant decreases in tibiae bone mineral density (≥9.41%, p≤0.006), axial rigidity (≥50.47%, p≤0.02), and soft tissue mass (55.25%, p=0.006) were observed from the trans-axial helical CT scans. The μCT results indicated significant detriments in tibia and calcaneus cortical thickness and bone volume fraction (p≤0.011). Significant changes in stiffness, yield load, ultimate load, and ultimate displacement (≥30.05%, p≤0.025) were observed from mechanical testing. These data indicate that limb disuse at a time of rapid musculoskeletal growth severely impairs muscle and bone development, reflecting the musculoskeletal complications observed in children with chronic medical conditions causing immobilization. Interventions to reduce these musculoskeletal complications in children are urgently needed. This disuse rabbit model will be useful in pre-clinical studies evaluating novel interventions for improving pediatric musculoskeletal health.

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.

Curved Beam Computed Tomography based Structural Rigidity Analysis of Bones with Simulated Lytic Defect: A Comparative Study with Finite Element Analysis

In this paper, a CT based structural rigidity analysis (CTRA) method that incorporates bone intrinsic local curvature is introduced to assess the compressive failure load of human femur with simulated lytic defects. The proposed CTRA is based on a three dimensional curved beam theory to obtain critical stresses within the human femur model. To test the proposed method, ten human cadaveric femurs with and without simulated defects were mechanically tested under axial compression to failure. Quantitative computed tomography images were acquired from the samples, and CTRA and finite element analysis were performed to obtain the failure load as well as rigidities in both straight and curved cross sections. Experimental results were compared to the results obtained from FEA and CTRA. The failure loads predicated by curved beam CTRA and FEA are in agreement with experimental results. The results also show that the proposed method is an efficient and reliable method to find both the location and magnitude of failure load. Moreover, the results show that the proposed curved CTRA outperforms the regular straight beam CTRA, which ignores the bone intrinsic curvature and can be used as a useful tool in clinical practices.

Noninvasive Failure Load Prediction of Vertebrae with Simulated Lytic Defects and Biomaterial Augmentation

The spine is the most common site for secondary bone metastases, and clinical management for fractures is based on size and geometry of the defect. About 75% of the bone needs to be damaged before lesions are detectable, so clinical tools should measure changes in both geometry and material properties. We have developed an automated, user-friendly, Spine Cancer Assessment (SCA) image-based analysis method that builds on a platform designed for clinical practice providing failure characteristics of vertebrae. The objectives of this study were to (1) validate SCA predictions with experimental failure load outcomes; (2) evaluate the planning capabilities for prophylactic vertebroplasty procedures; and (3) investigate the effect of computed tomography (CT) protocols on predicted failure loads. Twenty-one vertebrae were randomly divided into two groups: (1) simulated defect without treatment (negative control) [n = 9] and (2) with treatment [n = 12]. Defects were created and a polymeric biomaterial was injected into the vertebrae in the treated-defect group. Spines were scanned, reconstructed with two algorithms, and analyzed for fracture loads. To virtually plan for prophylactic intervention, vertebrae with empty lesions were simulated to be augmented with either poly(methyl methacrylate) (PMMA) or a novel bone replacement copolymer, poly(propylene fumarate-co-caprolactone) [P(PF-co-CL)]. Axial rigidities were calculated from the CT images. Failure loads, determined from the cross section with the lowest axial rigidity, were compared with experimental values. Predicted loads correlated well with experimental outcomes (R(2) = 0.73, p < 0.0001). Predictions from negative control specimens highly correlated with measured values (R(2) = 0.90, p < 0.0001). Although a similar correlation was obtained using both algorithms, the smooth reconstruction (B30) tended to underestimate predicted failure loads by ∼50% compared with the ∼10% underestimate of the sharp reconstruction (B70). Percent increase in failure loads after virtual vertebroplasty showed a higher increase in samples with PMMA compared with those with copolymer. The SCA method developed in this study calculated failure loads from quantitative computed tomography scans in vertebrae with simulated metastatic lytic defects, with or without treatment, facilitating clinical applicability and providing more reliable guidelines for physicians to select appropriate treatment options. Furthermore, the analysis could accommodate augmentation planning procedures that aimed to determine the optimum material that would increase vertebral body failure load.

Osteolytic changes around biodegradable cement restrictors in hip surgery

Background and purpose - Biodegradable cement restrictors are widely used in hip arthroplasty. Like others, we observed osteolytic reactions associated with a specific cement restrictor (SynPlug; made of PolyActive) and reviewed our patients. Patients and methods - We identified 703 patients with suitable radiographs from our database (2007 to 2012) who underwent cemented hip arthroplasty and received a SynPlug biodegradable cement restrictor. We reviewed all available radiographs to determine the incidence, severity, and progression of osteolysis. Mean postoperative follow-up was 1.8 (1-7) years Results - 1 year after implantation, the femoral cortex showed thinning by 12% in the anterior-posterior view and by 8% in the axial view. This had increased to 14% and 12%, respectively, at the latest available follow-up postoperatively (at a mean of 4 years). Cortical thinning of less than 10% was found in 37% of patients, but cortical thinning of 10-30% was found in 56% of patients. In the remaining 7%, a reduction of more than 30% of the original cortical thickness was observed. Interpretation - Osteolytic changes associated with the SynPlug biodegradable bone restrictors are inconsistent and highly variable. While some patients showed increased weakening of the femoral cortex with the potential risk of periprosthetic fracture, in others the degree of osteolysis only increased slightly or stabilized after 2 or more years. Any cortical bone loss after total hip replacement should be avoided, so the use of PolyActive biodegradable cement restrictors should be discontinued. Patients with a PolyActive cement restrictor in place should be followed up closely after surgery.

Bone mechanical properties and changes with osteoporosis

This review will define the role of collagen and within-bone heterogeneity and elaborate the importance of trabecular and cortical architecture with regard to their effect on the mechanical strength of bone. For each of these factors, the changes seen with osteoporosis and ageing will be described and how they can compromise strength and eventually lead to bone fragility.

ASSOCIATION BETWEEN COMPUTED TOMOGRAPHIC CHARACTERISTICS AND FRACTURES FOLLOWING STEREOTACTIC RADIOSURGERY IN DOGS WITH APPENDICULAR OSTEOSARCOMA

The objective of this observational, descriptive, retrospective study was to report CT characteristics associated with fractures following stereotactic radiosurgery in canine patients with appendicular osteosarcoma. Medical records (1999 and 2012) of dogs that had a diagnosis of appendicular osteosarcoma and undergone stereotactic radiosurgery were reviewed. Dogs were included in the study if they had undergone stereotactic radiosurgery for an aggressive bone lesion with follow-up information regarding fracture status, toxicity, and date and cause of death. Computed tomography details, staging, chemotherapy, toxicity, fracture status and survival data were recorded. Overall median survival time (MST) and fracture rates of treated dogs were calculated. CT characteristics were evaluated for association with time to fracture. Forty-six dogs met inclusion criteria. The median overall survival time was 9.7 months (95% CI: 6.9-14.3 months). The fracture-free rates at 3, 6, and 9 months were 73%, 44%, and 38% (95% CI: 60-86%, 29-60%, and 22-54%), respectively. The region of bone affected was significantly associated with time to fracture. The median time to fracture was 4.2 months in dogs with subchondral bone involvement and 16.3 months in dogs without subchondral bone involvement (P-value = 0.027, log-rank test). Acute and late skin effects were present in 58% and 16% of patients, respectively. Findings demonstrated a need for improved patient selection for this procedure, which can be aided by CT-based prognostic factors to predict the likelihood of fracture.

Augmentation of failed human vertebrae with critical un-contained lytic defect restores their structural competence under functional loading: An experimental study

BACKGROUND

Lytic spinal lesions reduce vertebral strength and may result in their fracture. Vertebral augmentation is employed clinically to provide mechanical stability and pain relief for vertebrae with lytic lesions. However, little is known about its efficacy in strengthening fractured vertebrae containing lytic metastasis.

METHODS

Eighteen unembalmed human lumbar vertebrae, having simulated uncontained lytic defects and tested to failure in a prior study, were augmented using a transpedicular approach and re-tested to failure using a wedge fracture model. Axial and moment based strength and stiffness parameters were used to quantify the effect of augmentation on the structural response of the failed vertebrae. Effects of cement volume, bone mineral density and vertebral geometry on the change in structural response were investigated.

FINDINGS

Augmentation increased the failed lytic vertebral strength [compression: 85% (P<0.001), flexion: 80% (P<0.001), anterior-posterior shear: 95%, P<0.001)] and stiffness [(40% (P<0.05), 53% (P<0.05), 45% (P<0.05)]. Cement volume correlated with the compressive strength (r(2)=0.47, P<0.05) and anterior-posterior shear strength (r(2)=0.52, P<0.05) and stiffness (r(2)=0.45, P<0.05). Neither the geometry of the failed vertebrae nor its pre-fracture bone mineral density correlated with the volume of cement.

INTERPRETATION

Vertebral augmentation is effective in bolstering the failed lytic vertebrae compressive and axial structural competence, showing strength estimates up to 50-90% of historical values of osteoporotic vertebrae without lytic defects. This modest increase suggests that lytic vertebrae undergo a high degree of structural damage at failure, with strength only partially restored by vertebral augmentation. The positive effect of cement volume is self-limiting due to extravasation.

Treatment Planning and Fracture Prediction in Patients with Skeletal Metastasis with CT-Based Rigidity Analysis

PURPOSE

Pathologic fractures could be prevented if reliable methods of fracture risk assessment were available. A multicenter prospective study was conducted to identify significant predictors of physicians' treatment plan for skeletal metastasis based on clinical fracture risk assessments and the proposed CT-based Rigidity Analysis (CTRA).

EXPERIMENTAL DESIGN

Orthopedic oncologists selected a treatment plan for 124 patients with 149 metastatic lesions based on the Mirels method. Then, CTRA was performed, and the results were provided to the physicians, who were asked to reassess their treatment plan. The pre- and post-CTRA treatment plans were compared to identify cases in which the treatment plan was changed based on the CTRA report. Patients were followed for a 4-month period to establish the incidence of pathologic fractures.

RESULTS

Pain, lesion type, and lesion size were significant predictors of the pre-CTRA plan. After providing the CTRA results, physicians changed their plan for 36 patients. CTRA results, pain, and primary source of metastasis were significant predictors of the post-CTRA plan. Follow-up of patients who did not undergo fixation resulted in 7 fractures; CTRA predicted these fractures with 100% sensitivity and 90% specificity, whereas the Mirels method was 71% sensitive and 50% specific.

CONCLUSIONS

Lesion type and size and pain level influenced the physicians' plans for the management of metastatic lesions. Physicians' treatment plans and fracture risk predictions were significantly influenced by the availability of CTRA results. Due to its high sensitivity and specificity, CTRA could potentially be used as a screening method for pathologic fractures.

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.

Computed tomography-based rigidity analysis: a review of the approach in preclinical and clinical studies

The assessment of fracture risk in patients afflicted with osseous neoplasms has long presented a problem for orthopedic oncologists. These patients are at risk for developing pathologic fractures through lytic defects in the appendicular and axial skeleton with devastating consequences on their quality of life. Lesions with a high risk of fracture may require prophylactic surgical stabilization, whereas low-risk lesions can be treated conservatively. Therefore, effective prevention of pathologic fractures depends on accurate assessment of fracture risk and is a critical step to avoid debilitating complications. Given the complex nature of osseous neoplasms, treatment requires a multidisciplinary approach; yet, little consensus regarding fracture risk assessment exists among physicians involved in the care of these patients. In order to improve the overall standard of care, specific criteria must be adopted to formulate consistent and accurate fracture risk predictions. However, clinicians make subjective assessments about fracture risk on plain radiographs using guidelines now recognized to be inaccurate. Osseous neoplasms alter both the material and geometric properties of bone; failure to account for changes in both of these parameters limits the accuracy of current fracture risk assessments. Rigidity, the capacity to resist deformation upon loading, is a structural property that integrates both the material and geometric properties of bone. Therefore, rigidity can be used as a mechanical assay of the changes induced by lytic lesions to the structural competency of bone. Using this principle, computed tomography (CT)-based structural rigidity analysis (CTRA) was developed and validated in a series of preclinical and clinical studies.

Finite element analysis and CT-based structural rigidity analysis to assess failure load in bones with simulated lytic defects

There is an urgent need to improve the prediction of fracture risk for cancer patients with bone metastases. Pathological fractures that result from these tumors frequently occur in the femur. It is extremely difficult to determine the fracture risk even for experienced physicians. Although evolving, fracture risk assessment is still based on inaccurate predictors estimated from previous retrospective studies. As a result, many patients are surgically over-treated, whereas other patients may fracture their bones against expectations. We mechanically tested ten pairs of human cadaveric femurs to failure, where one of each pair had an artificial defect simulating typical metastatic lesions. Prior to testing, finite element (FE) models were generated and computed tomography rigidity analysis (CTRA) was performed to obtain axial and bending rigidity measurements. We compared the two techniques on their capacity to assess femoral failure load by using linear regression techniques, Student's t-tests, the Bland-Altman methodology and Kendall rank correlation coefficients. The simulated FE failure loads and CTRA predictions showed good correlation with values obtained from the experimental mechanical testing. Kendall rank correlation coefficients between the FE rankings and the CTRA rankings showed moderate to good correlations. No significant differences in prediction accuracy were found between the two methods. Non-invasive fracture risk assessment techniques currently developed both correlated well with actual failure loads in mechanical testing suggesting that both methods could be further developed into a tool that can be used in clinical practice. The results in this study showed slight differences between the methods, yet validation in prospective patient studies should confirm these preliminary findings.

Fracture risk assessment in older adults using a combination of selected quantitative computed tomography bone measures: a subanalysis of the Age, Gene/Environment Susceptibility-Reykjavik Study

Bone mineral density (BMD) and geometric bone measures are individually associated with prevalent osteoporotic fractures. Whether an aggregate of these measures would better associate with fractures has not been examined. We examined relationships between self-reported fractures and selected bone measures acquired by quantitative computerized tomography (QCT), a composite bone score, and QCT-acquired dual-energy X-ray absorptiometry-like total femur BMD in 2110 men and 2682 women in the Age, Gene/Environment Susceptibility-Reykjavik Study. The combined bone score was generated by summing gender-specific Z-scores for 4 QCT measures: vertebral trabecular BMD, femur neck cortical thickness, femur neck trabecular BMD, and femur neck minimal cross-sectional area. Except for the latter measure, lower scores for QCT measures, singly and combined, showed positive (p < 0.05) associations with fractures. Results remained the same in stratified models for participants not taking bone-promoting medication. In women on bone-promoting medication, greater femur neck cortical thickness and trabecular BMD were significantly associated with fracture status. However, the association between fracture and combined bone score was not stronger than the associations between fracture and individual measures or total femur BMD. Thus, the selected measures did not all similarly associate with fracture status and did not appear to have an additive effect on fracture status.

The intravertebral distribution of bone density: correspondence to intervertebral disc health and implications for vertebral strength

This study's goal was to determine associations among the intravertebral heterogeneity in bone density, bone strength, and intervertebral disc (IVD) health. Results indicated that predictions of vertebral strength can benefit from considering the magnitude of the density heterogeneity and the congruence between the spatial distribution of density and IVD health.

INTRODUCTION

This study aims to determine associations among the intravertebral heterogeneity in bone density, bone strength, and IVD health

METHODS

Regional measurements of bone density were performed throughout 30 L1 vertebral bodies using micro-computed tomography (μCT) and quantitative computed tomography (QCT). The magnitude of the intravertebral heterogeneity in density was defined as the interquartile range and quartile coefficient of variation in regional densities. The spatial distribution of density was quantified using ratios of regional densities representing different anatomical zones (e.g., anterior to posterior regional densities). Cluster analysis was used to identify groups of vertebrae with similar spatial distributions of density. Vertebral strength was measured in compression. IVD health was assessed using two scoring systems.

RESULTS

QCT- and μCT-based measures of the magnitude of the intravertebral heterogeneity in density were strongly correlated with each other (p < 0.005). Accounting for the interquartile range in regional densities improved predictions of vertebral strength as compared to predictions based only on mean density (R (2) = 0.59 vs. 0.43; F-test p-value = 0.018). Specifically, after adjustment for mean density, vertebral bodies with greater heterogeneity in density exhibited higher strength. No single spatial distribution of density was associated with high vertebral strength. Analyses of IVD scores suggested that the health of the adjacent IVDs may modulate the effect of a particular spatial distribution of density on vertebral strength.

CONCLUSIONS

Noninvasive measurements of the intravertebral distribution of bone density, in conjunction with assessments of IVD health, can aid in predictions of bone strength and in elucidating biomechanical mechanisms of vertebral fracture.

Vertebral compression fracture after stereotactic body radiotherapy for spinal metastases

The use of stereotactic body radiotherapy for metastatic spinal tumours is increasing. Serious adverse events for this treatment include vertebral compression fracture (VCF) and radiation myelopathy. Although VCF is a fairly low-risk adverse event (approximately 5% risk) after conventional radiotherapy, crude risk estimates for VCF after spinal SBRT range from 11% to 39%. In this Review, we summarise the evidence and predictive factors for VCF induced by spinal SBRT, review the pathophysiology of VCF in the metastatic spine, and discuss strategies used to prevent and manage this potentially disabling complication.

Finite element analysis and computed tomography based structural rigidity analysis of rat tibia with simulated lytic defects

Finite element analysis (FEA), CT based structural rigidity analysis (CTRA) and mechanical testing is performed to assess the compressive failure load of rat tibia with simulated lytic defects. Twenty rat tibia were randomly assigned to four equal groups (n=5). Three of the groups included a simulated defect at various locations: anterior bone surface (Group 1), posterior bone surface (Group 2) and through bone defect (Group 3). The fourth group was a control group with no defect (Group 4). Microcomputed tomography was used to assess bone structural rigidity properties and to provide 3D model data for generation of the finite element models for each specimen. Compressive failure load calculated using CT derived rigidity parameters (FCTRA) was well correlated to failure load recorded in mechanical testing (R(2)=0.96). The relationships between mechanical testing failure load and the axial rigidity (R(2)=0.61), bending rigidity (R(2)=0.71) and FEA calculated failure loads, considering bone as an elastic isotropic (R(2)=0.75) and elastic transversely isotropic (R(2)=0.90) are also well correlated. CTRA stress, calculated adjacent to the defect, were also shown to be well correlated with yield stresses calculated using the minimum density at the weakest cross section (R(2)=0.72). No statistically significant relationship between apparent density and mechanical testing failure load was found (P=0.37). In summary, the results of this study indicate that CTRA analysis of bone strength correlates well with both FEA and results obtained from compression experiments. In addition there exist a good correlation between structural rigidity parameters and experimental failure loads. In contrast, there was no correlation between average bone density and failure load.

The effect of intravertebral heterogeneity in microstructure on vertebral strength and failure patterns

The goal of this study was to determine the influence of intravertebral heterogeneity in microstructure on vertebral failure. Results show that noninvasive assessments of the intravertebral heterogeneity in density improve predictions of vertebral strength and that local variations in microstructure are associated with locations of failure in the vertebral body.

INTRODUCTION

The overall goal of this study was to determine the influence of intravertebral heterogeneity in microstructure on vertebral failure.

METHODS

Trabecular density and microarchitecture were quantified for 32 thoracic vertebrae using micro-computed tomography (μCT)-based analyses of 4.81 mm, contiguous cubes throughout the centrum. Intravertebral heterogeneity in density was defined as the interquartile range and quartile coefficient of variation of the cube densities. The vertebrae were compressed to failure to measure stiffness, strength, and toughness. Pre- and post-compression μCT images were analyzed using digital volume correlation to quantify failure patterns in the vertebrae, as defined by the distributions of residual strain.

RESULTS

Failure patterns consisted of large deformations in the midtransverse plane with concomitant endplate biconcavity and were linked to the intravertebral distribution of bone tissue. Low values of connectivity density and trabecular number, and high values of trabecular separation, were associated with high strains. However, local microstructural properties were not the sole determinants of failure. For instance, the midtransverse plane experienced the highest strain (p < 0.008) yet had the highest density, lowest structure model index, and lowest anisotropy (p < 0.013). Accounting for the intravertebral heterogeneity in density improved predictions of strength and stiffness as compared to predictions based only on mean density (strength: R(2) = 0.75 vs. 0.61, p < 0.001; stiffness: R(2) = 0.44 vs. 0.26, p = 0.001).

CONCLUSIONS

Local variations in microstructure are associated with failure patterns in the vertebra. Noninvasive assessments of the intravertebral heterogeneity in density--which are feasible in clinical settings--can improve predictions of vertebral strength and stiffness.

Can micro-imaging based analysis methods quantify structural integrity of rat vertebrae with and without metastatic involvement?

This study compares the ability of μCT image-based registration, 2D structural rigidity analyses and multimodal continuum-level finite element (FE) modeling in evaluating the mechanical stability of healthy, osteolytic, and mixed osteolytic/osteoblastic metastatically involved rat vertebrae. μMR and μCT images (loaded and unloaded) were acquired of lumbar spinal motion segments from 15rnu/rnu rats (five per group). Strains were calculated based on image registration of the loaded and unloaded μCT images and via analysis of FE models created from the μCT and μMR data. Predicted yield load was also calculated through 2D structural rigidity analysis of the axial unloaded μCT slices. Measures from the three techniques were compared to experimental yield loads. The ability of these methods to predict experimental yield loads were evaluated and image registration and FE calculated strains were directly compared. Quantitatively for all samples, only limited weak correlations were found between the image-based measures and experimental yield load. In comparison to the experimental yield load, we observed a trend toward a weak negative correlation with median strain calculated using the image-based strain measurement algorithm (r=-0.405, p=0.067), weak significant correlations (p<0.05) with FE based median and 10th percentile strain values (r=-0.454, -0.637, respectively), and a trend toward a weak significant correlation with FE based mean strain (r=-0.366, p=0.09). Individual group analyses, however, yielded more and stronger correlations with experimental results. Considering the image-based strain measurement algorithm we observed moderate significant correlations with experimental yield load (p<0.05) in the osteolytic group for mean and median strain values (r=-0.840, -0.832, respectively), and in the healthy group for median strain values (r=-0.809). Considering the rigidity-based predicted yield load, we observed a strong significant correlation with the experimental yield load in the mixed osteolytic/osteoblastic group (r=0.946) and trend toward a moderate correlation with the experimental yield load in the osteolytic group (r=0.788). Qualitatively, strain patterns in the vertebral bodies generated using image registration and FEA were well matched, yet quantitatively a significant correlation was found only between mean strains in the healthy group (r=0.934). Large structural differences in metastatic vertebrae and the complexity of motion segment loading may have led to varied modes of failure. Improvements in load characterization, material properties assignments and resolution are necessary to yield a more generalized ability for image-based registration, structural rigidity and FE methods to accurately represent stability in healthy and pathologic scenarios.

Assessment of axial bone rigidity in rats with metabolic diseases using CT-based structural rigidity analysis

Objectives

This study aims to assess the correlation of CT-based structural rigidity analysis with mechanically determined axial rigidity in normal and metabolically diseased rat bone.

Methods

A total of 30 rats were divided equally into normal, ovariectomized, and partially nephrectomized groups. Cortical and trabecular bone segments from each animal underwent micro-CT to assess their average and minimum axial rigidities using structural rigidity analysis. Following imaging, all specimens were subjected to uniaxial compression and assessment of mechanically-derived axial rigidity.

Results

The average structural rigidity-based axial rigidity was well correlated with the average mechanically-derived axial rigidity results (R2 = 0.74). This correlation improved significantly (p < 0.0001) when the CT-based Structural Rigidity Analysis (CTRA) minimum axial rigidity was correlated to the mechanically-derived minimum axial rigidity results (R2 = 0.84). Tests of slopes in the mixed model regression analysis indicated a significantly steeper slope for the average axial rigidity compared with the minimum axial rigidity (p = 0.028) and a significant difference in the intercepts (p = 0.022). The CTRA average and minimum axial rigidities were correlated with the mechanically-derived average and minimum axial rigidities using paired t-test analysis (p = 0.37 and p = 0.18, respectively).

Conclusions

In summary, the results of this study suggest that structural rigidity analysis of micro-CT data can be used to accurately and quantitatively measure the axial rigidity of bones with metabolic pathologies in an experimental rat model. It appears that minimum axial rigidity is a better model for measuring bone rigidity than average axial rigidity.

Assessment of axial bone rigidity in rats with metabolic diseases using CT-based structural rigidity analysis

Objectives

This study aims to assess the correlation of CT-based structural rigidity analysis with mechanically determined axial rigidity in normal and metabolically diseased rat bone.

Methods

A total of 30 rats were divided equally into normal, ovariectomized, and partially nephrectomized groups. Cortical and trabecular bone segments from each animal underwent micro-CT to assess their average and minimum axial rigidities using structural rigidity analysis. Following imaging, all specimens were subjected to uniaxial compression and assessment of mechanically-derived axial rigidity.

Results

The average structural rigidity-based axial rigidity was well correlated with the average mechanically-derived axial rigidity results (R² = 0.74). This correlation improved significantly (p < 0.0001) when the CT-based Structural Rigidity Analysis (CTRA) minimum axial rigidity was correlated to the mechanically-derived minimum axial rigidity results (R² = 0.84). Tests of slopes in the mixed model regression analysis indicated a significantly steeper slope for the average axial rigidity compared with the minimum axial rigidity (p = 0.028) and a significant difference in the intercepts (p = 0.022). The CTRA average and minimum axial rigidities were correlated with the mechanically-derived average and minimum axial rigidities using paired t-test analysis (p = 0.37 and p = 0.18, respectively).

Conclusions

In summary, the results of this study suggest that structural rigidity analysis of micro-CT data can be used to accurately and quantitatively measure the axial rigidity of bones with metabolic pathologies in an experimental rat model. It appears that minimum axial rigidity is a better model for measuring bone rigidity than average axial rigidity.

Development of quantitative computed-tomography-based strength indicators for the identification of low bone-strength individuals in a clinical environment

The aim of this study was to develop quantitative computed-tomography (QCT)-based bone-strength indicators that highly correlate with finite-element (FE)-based strength. Transaxial QCT scans were obtained from 36 major, cadaveric, long bones (humerus, radius, femur and tibia) from 4 females and 2 males, 53 to 86 years old. These images were used to construct the FE models and to develop the QCT-based bone strength indicators under every-day, simplified loading conditions. We have evaluated the performance of area-weighted (AW), density-weighted (DW) and modulus-weighted (MW) rigidity measures as well as popular strength indicators like section modulus (Z) and stress-strain index (SSI). We have also developed a novel strength metric, the centroid deviation, which analyzes the spatial distribution of the centroids along the length of the bone. The correlation results show that the MW polar moment of inertia and the MW moment of inertia are the two top-performers for all bones and loading conditions (average r>0.89). The MW centroid deviations correlated highly with the estimated load to fracture for all bones under compression (r>0.83), except for the humerus (r=0.67). Consistently DW or MW rigidity measures produced a statistically significant improvement in capturing bone strength compared to AW rigidity measures. As expected, MW rigidity measures showed a higher correlation with the FE-based fracture load than the DW rigidity measures; however, the improvement was not statistically significant. Through this study we present a short-list of useful QCT-based strength parameters that correlate well with FE-based fracture load. Although a few parameters perform reasonably well across most bones and loading conditions, a judicious assessment of bone strength should include multiple parameters evaluated at multiple critical locations in the long bones, with attention to the type of loading and bone type.