MDCT-based Finite Element Analysis of Vertebral Fracture Risk: What Dose is Needed?

Characterization of mechanical stiffness using additive manufacturing and finite element analysis: potential tool for bone health assessment

BACKGROUND

Bone health and fracture risk are known to be correlated with stiffness. Both micro-finite element analysis (μFEA) and mechanical testing of additive manufactured phantoms are useful approaches for estimating mechanical properties of trabecular bone-like structures. However, it is unclear if measurements from the two approaches are consistent. The purpose of this work is to evaluate the agreement between stiffness measurements obtained from mechanical testing of additive manufactured trabecular bone phantoms and μFEA modeling. Agreement between the two methods would suggest 3D printing is a viable method for validation of μFEA modeling.

METHODS

A set of 20 lumbar vertebrae regions of interests were segmented and the corresponding trabecular bone phantoms were produced using selective laser sintering. The phantoms were mechanically tested in uniaxial compression to derive their stiffness values. The stiffness values were also derived from in silico simulation, where linear elastic μFEA was applied to simulate the same compression and boundary conditions. Bland-Altman analysis was used to evaluate agreement between the mechanical testing and μFEA simulation values. Additionally, we evaluated the fidelity of the 3D printed phantoms as well as the repeatability of the 3D printing and mechanical testing process.

RESULTS

We observed good agreement between the mechanically tested stiffness and μFEA stiffness, with R2 of 0.84 and normalized root mean square deviation of 8.1%. We demonstrate that the overall trabecular bone structures are printed in high fidelity (Dice score of 0.97 (95% CI, [0.96,0.98]) and that mechanical testing is repeatable (coefficient of variation less than 5% for stiffness values from testing of duplicated phantoms). However, we noticed some defects in the resin microstructure of the 3D printed phantoms, which may account for the discrepancy between the stiffness values from simulation and mechanical testing.

CONCLUSION

Overall, the level of agreement achieved between the mechanical stiffness and μFEA indicates that our μFEA methods may be acceptable for assessing bone mechanics of complex trabecular structures as part of an analysis of overall bone health.

Computed Tomography of the Spine

The introduction of the first whole-body CT scanner in 1974 marked the beginning of cross-sectional spine imaging. In the last decades, the technological advancement, increasing availability and clinical success of CT led to a rapidly growing number of CT examinations, also of the spine. After initially being primarily used for trauma evaluation, new indications continued to emerge, such as assessment of vertebral fractures or degenerative spine disease, preoperative and postoperative evaluation, or CT-guided interventions at the spine; however, improvements in patient management and clinical outcomes come along with higher radiation exposure, which increases the risk for secondary malignancies. Therefore, technical developments in CT acquisition and reconstruction must always include efforts to reduce the radiation dose. But how exactly can the dose be reduced? What amount of dose reduction can be achieved without compromising the clinical value of spinal CT examinations and what can be expected from the rising stars in CT technology: artificial intelligence and photon counting CT? In this article, we try to answer these questions by systematically reviewing dose reduction techniques with respect to the major clinical indications of spinal CT. Furthermore, we take a concise look on the dose reduction potential of future developments in CT hardware and software.

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

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

MDCT-Based Finite Element Analysis for the Prediction of Functional Spine Unit Strength-An In Vitro Study

(1) Objective: This study aimed to analyze the effect of ligaments on the strength of functional spine unit (FSU) assessed by finite element (FE) analysis of anatomical models developed from multi-detector computed tomography (MDCT) data. (2) Methods: MDCT scans for cadaveric specimens were acquired from 16 donors (7 males, mean age of 84.29 ± 6.06 years and 9 females, mean age of 81.00 ± 11.52 years). Two sets of FSU models (three vertebrae + two disks), one with and another without (w/o) ligaments, were generated. The vertebrae were segmented semi-automatically, intervertebral disks (IVD) were generated manually, and ligaments were modeled based on the anatomical location. FE-predicted failure loads of FSU models (with and w/o ligaments) were compared with the experimental failure loads obtained from the uniaxial biomechanical test of specimens. (3) Results: The mean and standard deviation of the experimental failure load of FSU specimens was 3513 ± 1029 N, whereas of FE-based failure loads were 2942 ± 943 N and 2537 ± 929 N for FSU models with ligaments and without ligament attachments, respectively. A good correlation (ρ = 0.79, and ρ = 0.75) was observed between the experimental and FE-based failure loads for the FSU model with and with ligaments, respectively. (4) Conclusions: The FE-based FSU model can be used to determine bone strength, and the ligaments seem to have an effect on the model accuracy for the failure load calculation; further studies are needed to understand the contribution of ligaments.

MDCT-Based Finite Element Analyses: Are Measurements at the Lumbar Spine Associated with the Biomechanical Strength of Functional Spinal Units of Incidental Osteoporotic Fractures along the Thoracolumbar Spine?

Assessment of osteoporosis-associated fracture risk during clinical routine is based on the evaluation of clinical risk factors and T-scores, as derived from measurements of areal bone mineral density (aBMD). However, these parameters are limited in their ability to identify patients at high fracture risk. Finite element models (FEMs) have shown to improve bone strength prediction beyond aBMD. This study aims to investigate whether FEM measurements at the lumbar spine can predict the biomechanical strength of functional spinal units (FSUs) with incidental osteoporotic vertebral fractures (VFs) along the thoracolumbar spine. Multi-detector computed tomography (MDCT) data of 11 patients (5 females and 6 males, median age: 67 years) who underwent MDCT twice (median interval between baseline and follow-up MDCT: 18 months) and sustained an incidental osteoporotic VF between baseline and follow-up scanning were used. Based on baseline MDCT data, two FSUs consisting of vertebral bodies and intervertebral discs (IVDs) were modeled: one standardly capturing L1-IVD-L2-IVD-L3 (FSU_L1-L3) and one modeling the incidentally fractured vertebral body at the center of the FSU (FSU_F). Furthermore, volumetric BMD (vBMD) derived from MDCT, FEM-based displacement, and FEM-based load of the single vertebrae L1 to L3 were determined. Statistically significant correlations (adjusted for a BMD ratio of fracture/L1-L3 segments) were revealed between the FSU_F and mean load of L1 to L3 (r = 0.814, p = 0.004) and the mean vBMD of L1 to L3 (r = 0.745, p = 0.013), whereas there was no statistically significant association between the FSU_F and FSU_L1-L3 or between FSU_F and the mean displacement of L1 to L3 (p > 0.05). In conclusion, FEM measurements of single vertebrae at the lumbar spine may be able to predict the biomechanical strength of incidentally fractured vertebral segments along the thoracolumbar spine, while FSUs seem to predict only segment-specific fracture risk.

Prediction of Incidental Osteoporotic Fractures at Vertebral-Specific Level Using 3D Non-Linear Finite Element Parameters Derived from Routine Abdominal MDCT

To investigate whether finite element (FE) analysis of the spine in routine thoracic/abdominal multi-detector computed tomography (MDCT) can predict incidental osteoporotic fractures at vertebral-specific level; Baseline routine thoracic/abdominal MDCT scans of 16 subjects (8(m), mean age: 66.1 ± 8.2 years and 8(f), mean age: 64.3 ± 9.5 years) who sustained incidental osteoporotic vertebral fractures as confirmed in follow-up MDCTs were included in the current study. Thoracic and lumbar vertebrae (T5-L5) were automatically segmented, and bone mineral density (BMD), finite element (FE)-based failure-load, and failure-displacement were determined. These values of individual vertebrae were normalized globally (g), by dividing the absolute value with the average of L1-3 and locally by dividing the absolute value with the average of T5-12 and L1-5 for thoracic and lumbar vertebrae, respectively. Mean-BMD of L1-3 was determined as reference. Receiver operating characteristics (ROC) and area under the curve (AUC) were calculated for different normalized FE (Kload, Kdisplacement,K(load)g, and K(displacement)g) and BMD (KBMD, and K(BMD)g) ratio parameter combinations for identifying incidental fractures. Kload, K(load)g, KBMD, and K(BMD)g showed significantly higher discriminative power compared to standard mean BMD of L1-3 (BMDStandard) (AUC = 0.67 for Kload; 0.64 for K(load)g; 0.64 for KBMD; 0.61 for K(BMD)g vs. 0.54 for BMDStandard). The combination of Kload, Kdisplacement, and KBMD increased the AUC further up to 0.77 (p < 0.001). The combination of FE with BMD measurements derived from routine thoracic/abdominal MDCT allowed an improved prediction of incidental fractures at vertebral-specific level.

Predicting Vertebral Bone Strength Using Finite Element Analysis for Opportunistic Osteoporosis Screening in Routine Multidetector Computed Tomography Scans-A Feasibility Study

Purpose

To investigate the feasibility of using routine clinical multidetector computed tomography (MDCT) scans for conducting finite element (FE) analysis to predict vertebral bone strength for opportunistic osteoporosis screening.

Methods

Routine abdominal MDCT with and without intravenous contrast medium (IVCM) of seven subjects (five male; two female; mean age: 71.86 ± 7.40 years) without any bone disease were used. FE analysis was performed on individual vertebrae (T11, T12, L1, and L2) including the posterior elements to investigate the effect of IVCM and slice thickness (1 and 3 mm) on vertebral bone strength. Another subset of data from subjects with vs. without osteoporotic vertebral fractures (n = 9 age and gender-matched pairs) was analyzed for investigating the ability of FE-analysis to differentiate the two cohorts. Bland-Altman plots, box plots, and coefficient of correlation (R2) were calculated to determine the variations in FE-predicted failure loads for different conditions.

Results

The FE-predicted failure loads obtained from routine MDCT scans were strongly correlated with those from without IVCM (R2 = 0.91 for 1mm; R2 = 0.92 for 3mm slice thickness, respectively) and different slice thicknesses (R2 = 0.93 for 1mm vs. 3mm with IVCM). Furthermore, a good correlation was observed for 3mm slice thickness with IVCM vs. 1mm without IVCM (R2 = 0.87). Significant difference between FE-predicted failure loads of healthy and fractured patients was observed (4,705 ± 1,238 vs. 4,010 ± 1,297 N; p=0.026).

Conclusion

Routine clinical MDCT scans could be reliably used for assessment of fracture risk based on FE analysis and may be beneficial for patients who are at increased risk for osteoporotic fractures.

Low-dose and sparse sampling MDCT-based femoral bone strength prediction using finite element analysis

This study aims to evaluate the impact of dose reduction through tube current and sparse sampling on multi-detector computed tomography (MDCT)-based femoral bone strength prediction using finite element (FE) analysis. FE-predicted femoral failure load obtained from MDCT scan data was not significantly affected by 50% dose reductions through sparse sampling. Further decrease in dose through sparse sampling (25% of original projections) and virtually reduced tube current (50% and 25% of the original dose) showed significant effects on the FE-predicted failure load results.

PURPOSE

To investigate the effect of virtually reduced tube current and sparse sampling on multi-detector computed tomography (MDCT)-based femoral bone strength prediction using finite element (FE) analysis.

METHODS

Routine MDCT data covering the proximal femur of 21 subjects (17 males; 4 females; mean age, 71.0 ± 8.8 years) without any bone diseases aside from osteoporosis were included in this study. Fifty percent and 75% dose reductions were achieved by virtually reducing tube current and by applying a sparse sampling strategy from the raw image data. Images were then reconstructed with a statistically iterative reconstruction algorithm. FE analysis was performed on all reconstructed images and the failure load was calculated. The root mean square coefficient of variation (RMSCV) and coefficient of correlation (R²) were calculated to determine the variation in the FE-predicted failure load data for dose reductions, using original-dose MDCT scan as the standard of reference.

RESULTS

Fifty percent dose reduction through sparse sampling showed lower RMSCV and higher correlations when compared with virtually reduced tube current method (RMSCV = 5.70%, R² = 0.96 vs. RMSCV = 20.78%, R² = 0.79). Seventy-five percent dose reduction achieved through both methods (RMSCV = 22.38%, R² = 0.80 for sparse sampling; RMSCV = 24.58%, R² = 0.73 for reduced tube current) could not predict the failure load accurately.

CONCLUSION

Our simulations indicate that up to 50% reduction in radiation dose through sparse sampling can be used for FE-based prediction of femoral failure load. Sparse-sampled MDCT may allow fracture risk prediction and treatment monitoring in osteoporosis with less radiation exposure in the future.

X-ray-based quantitative osteoporosis imaging at the spine

Osteoporosis is a metabolic bone disease with a high prevalence that affects the population worldwide, particularly the elderly. It is often due to fractures associated with bone fragility that the diagnosis of osteoporosis becomes clinically evident. However, early diagnosis would be necessary to initiate therapy and to prevent occurrence of further fractures, thus reducing morbidity and mortality. X-ray-based imaging plays a key role for fracture risk assessment and monitoring of osteoporosis. Whereas over decades dual-energy X-ray absorptiometry (DXA) has been the main method used and still reflects the reference standard, another modality reemerges with quantitative computed tomography (QCT) because of its three-dimensional advantages and the opportunistic exploitation of routine CT scans. Against this background, this article intends to review and evaluate recent advances in the field of X-ray-based quantitative imaging of osteoporosis at the spine. First, standard DXA with the recent addition of trabecular bone score (TBS) is presented. Secondly, standard QCT, dual-energy BMD quantification, and opportunistic BMD screening in non-dedicated CT exams are discussed. Lastly, finite element analysis and microstructural parameter analysis are reviewed.

Finite Element Analysis-Based Vertebral Bone Strength Prediction Using MDCT Data: How Low Can We Go?

Objective: To study the impact of dose reduction in MDCT images through tube current reduction or sparse sampling on the vertebral bone strength prediction using finite element (FE) analysis for fracture risk assessment. Methods: Routine MDCT data covering lumbar vertebrae of 12 subjects (six male; six female; 74.70 ± 9.13 years old) were included in this study. Sparsely sampled and virtually reduced tube current-based MDCT images were computed using statistical iterative reconstruction (SIR) with reduced dose levels at 50, 25, and 10% of the tube current and original projections, respectively. Subject-specific static non-linear FE analyses were performed on vertebra models (L1, L2, and L3) 3-D-reconstructed from those dose-reduced MDCT images to predict bone strength. Coefficient of correlation (R²), Bland-Altman plots, and root mean square coefficient of variation (RMSCV) were calculated to find the variation in the FE-predicted strength at different dose levels, using high-intensity dose-based strength as the reference. Results: FE-predicted failure loads were not significantly affected by up to 90% dose reduction through sparse sampling (R² = 0.93, RMSCV = 8.6% for 50%; R² = 0.89, RMSCV = 11.90% for 75%; R² = 0.86, RMSCV = 11.30% for 90%) and up to 50% dose reduction through tube current reduction method (R² = 0.96, RMSCV = 12.06%). However, further reduction in dose with the tube current reduction method affected the ability to predict the failure load accurately (R² = 0.88, RMSCV = 22.04% for 75%; R² = 0.43, RMSCV = 54.18% for 90%). Conclusion: Results from this study suggest that a 50% radiation dose reduction through reduced tube current and a 90% radiation dose reduction through sparse sampling can be used to predict vertebral bone strength. Our findings suggest that the sparse sampling-based method performs better than the tube current-reduction method in generating images required for FE-based bone strength prediction models.