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

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.

Single-level subject-specific finite element model can predict fracture outcomes in three-level spine segments under different loading rates

Osteoporosis-related vertebral compression fracture can occur under normal physiological activities. Bone metastasis is another source of vertebral fracture. Different loading rates, either high-energy traumas such as falls or low-energy traumas under normal physiological activities, can result in different fracture outcomes. The aim of the current study was to develop a quantitative computed tomography-based finite element analysis (QCT/FEA) technique for single vertebral bodies to predict fracture strength of three-level spine segments. Developed QCT/FEA technique was also used to characterize vertebral elastic moduli at two loading rates of 5 mm/min, representing a physiologic loading condition, and 12000 mm/min, representing a high-energy trauma. To this end, a cohort of human spine segments divided into three groups of intact, defect, and augmented were mechanically tested to fracture; then, experimental stiffness and fracture strength values were measured. Outcomes of this study showed no significant difference between the elastic modulus equations at the two testing speeds. Areal bone mineral density measured by dual x-ray absorptiometry (DXA/BMD) explained only 53% variability (R² = 0.53) in fracture strength outcomes. However, QCT/FEA could explain 70% of the variability (R² = 0.70) in experimentally measured fracture strength values. Adding disk degeneration grading, testing speed, and sex to QCT/FEA-estimated fracture strength values further increased the performance of our statistical model by 14% (adjusted R² of 0.84 between the prediction and experimental fracture forces). In summary, our results indicated that a single-vertebra model, which is computationally less expensive and more time efficient, is capable of estimating fracture outcomes with acceptable performance (range: 70-84%).

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.

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.

Effect of metastatic lesion size and location on the load-bearing capacity of vertebrae using an optimized ash density-modulus equation

About 1.8 million new cancer cases are estimated in the US in 2019 from which 50-85% might metastasize to the thoracic and lumbar spines. Subject-specific quantitative computed tomography-based finite element analysis (QCT/FEA) is a promising used tool to predict vertebral fracture properties. The aims of this study were twofold: First, to develop an optimized equation for the elastic modulus accounting for all input parameters in FE modeling of fracture properties. Second, to assess the effect of lesion size and location on the predicted fracture loads. An inverse QCT/FEA method was implemented to determine optimal coefficients for the modulus equation as a function of ash density. Lesions of 16 and 20 mm were then virtually located at the center, off-centered, anterior, and posterior regions of the vertebrae. A total of 6426 QCT/FEA models were run to optimize the coefficients and evaluate the effect of lesions on fracture properties. QCT/FEA predicted stiffness showed high correlations (50%) with the experimentally measured values. Compared to a 16 mm lesion size, a 20 mm lesion had a reduction in failure load of 55%, 57%, 52%, and 44% at the center, off-centered, anterior cortex, and pedicle, respectively (p < 0.001). Lesions affecting mostly trabecular bone showed the largest reduction in predicted failure loads (about 55%), and females presented weaker outcomes than males. An optimal elastic modulus equation resulted in accurate vertebral stiffness predictions. A deterioration of the trabecular bone due to the presence of a lesion highly affected the predicted fracture loads, and this reduction was significantly higher in females compared to males.

Dual-energy estimates of volumetric bone mineral densities in the lumbar spine using quantitative computed tomography better correlate with fracture properties when compared to single-energy BMD outcomes

It is estimated that over 200 million people worldwide are affected by osteoporosis. Vertebral fracture risk prediction using dual energy x-ray absorptiometry (DXA) is confounded by limitations of the technology, such as 2D measurements of bone mineral density (BMD), inability to measure bone distribution and heterogeneity, and potential overestimations of BMD due to degenerative diseases. To overcome these shortcomings, single energy (SE) quantitative computed tomography (QCT) imaging estimates of Hounsfield units (HU) and volumetric BMD have been implemented as alternative methodologies for assessing fracture risk. However, marrow fat within the vertebrae can highly affect the vBMD and fracture properties estimations. To address this issue, 54 vertebrae were dissected from nine cadaveric spines and scanned using SE-QCT (120kVp) and dual energy (DE)-QCT (80/140 kVp), with the latter accounting for marrow fat within the vertebrae. The vertebrae were then scanned using DXA and subjected to mechanical testing to obtain fracture properties. aBMD outcomes from DXA showed a better correlation with DE-QCT vBMD versus SE outcomes [DE: aBMD vs. vBMD (R²: 0.61); SE: aBMD vs. vBMD (R²: 0.27)]. SE-QCT underestimated vertebral vBMD by -56% (p<0.0001) when compared to DE-QCT. vBMD estimates from SE-QCT could predict 45% and 37% of the vertebral failure loads and stiffness, respectively, compared to 67% and 46% from DE-QCT. DE-QCT vBMD outcomes highly correlated with fracture properties of vertebrae as compared to SE-QCT metrics. As DE scanning has the ability to correct for the effects of bone marrow fat, estimated vBMD from SE-QCT were significantly underestimated compared to DE-QCT. Dual energy CT scanning has the potential to more accurately predict vertebral failure and aid the clinician in the evaluation of appropriate interventions. Future studies should consider implementing DE-QCT in their fracture assessment.

Instability and impending instability in patients with vertebral metastatic disease

Metastatic disease commonly involves the spine with an increasing incidence due to a worldwide rise of cancer incidence and a longer survival of patients with osseous metastases. Metastases compromise the mechanical integrity of the vertebra and make it susceptible to fracture. Patients with pathological vertebral fracture often become symptomatic, with mechanical pain generally due to intervertebral instability, and may develop spinal cord compression and neurological deficits. Advances in imaging, radiotherapy, as well as in spinal surgery techniques, have allowed the evolution from conventional palliative external beam radiotherapy to modern stereotactic radiosurgery and from traditional open surgery to less-invasive, and sometimes prophylactic stabilization surgical treatments. It is therefore clear that fracture risk prediction, and maintenance or restoration of intervertebral stability, are important objectives in the management of these patients. Correlation between imaging findings and clinical manifestations is crucial, and a common knowledge base for treatment team members rather than a compartmentalized view is very important. This article reviews the literature on the imaging and clinical diagnosis of intervertebral instability and impending instability in the setting of spine metastatic disease, including the spinal instability neoplastic score, which is a reliable tool for diagnosing unstable or potentially unstable metastatic spinal lesions, and on the different elements considered for treatment.

Variation in human vertebral body strength for vertebral body samples from different locations in segments L1-L5

BACKGROUND

The human spine, in particular the lumbar spine, is subject to significant compressive and bending stresses, which affect the structure of the bone tissue of the vertebrae. The more heterogeneous the structure of the spongy bone tissue, the less resistant the whole vertebral body. It is therefore necessary to establish variations in bone strength parameters within one particular vertebral body.

METHODS

The research material comprised human L1-L5 lumbar vertebrae sampled from 15 donors aged 29-35. A total of 975 samples prepared from the collected material were subjected to compressive and bending strength tests. The samples for the tests were collected from carefully selected locations in order to discover the strength properties of various parts of the vertebral body.

FINDINGS

In the case of sample 2 (located in the posterior part of the vertebra, at mid-height) the stress values were the lowest and there were statistically significant differences compared to other samples. Moreover the value of compressive force in this case was lower for vertebrae with higher numbers. Top and bottom samples demonstrated statistically significant higher mean values of destructive stress. In terms of the bending strength test, the mean value of destructive stress in all lumbar vertebrae for all samples increased for vertebrae with higher numbers.

INTERPRETATION

The spongy tissue in healthy vertebral bodies has a very heterogeneous structure. This may be due to the presence of the nutrient canal and the arc structure allowing more springy movement and improved transfer of loads by the vertebral body.

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.