Patient-Specific Finite Element Modeling of the Whole Lumbar Spine Using Clinical Routine Multi-Detector Computed Tomography (MDCT) Data-A Pilot Study.

(1) Background: To study the feasibility of developing finite element (FE) models of the whole lumbar spine using clinical routine multi-detector computed tomography (MDCT) scans to predict failure load (FL) and range of motion (ROM) parameters. (2) Methods: MDCT scans of 12 subjects (6 healthy controls (HC), mean age ± standard deviation (SD): 62.16 ± 10.24 years, and 6 osteoporotic patients (OP), mean age ± SD: 65.83 ± 11.19 years) were included in the current study. Comprehensive FE models of the lumbar spine (5 vertebrae + 4 intervertebral discs (IVDs) + ligaments) were generated (L1−L5) and simulated. The coefficients of correlation (ρ) were calculated to investigate the relationship between FE-based FL and ROM parameters and bone mineral density (BMD) values of L1−L3 derived from MDCT (BMDQCT-L1-3). Finally, Mann−Whitney U tests were performed to analyze differences in FL and ROM parameters between HC and OP cohorts. (3) Results: Mean FE-based FL value of the HC cohort was significantly higher than that of the OP cohort (1471.50 ± 275.69 N (HC) vs. 763.33 ± 166.70 N (OP), p < 0.01). A strong correlation of 0.8 (p < 0.01) was observed between FE-based FL and BMDQCT-L1-L3 values. However, no significant differences were observed between ROM parameters of HC and OP cohorts (p = 0.69 for flexion; p = 0.69 for extension; p = 0.47 for lateral bending; p = 0.13 for twisting). In addition, no statistically significant correlations were observed between ROM parameters and BMDQCT- L1-3. (4) Conclusions: Clinical routine MDCT data can be used for patient-specific FE modeling of the whole lumbar spine. ROM parameters do not seem to be significantly altered between HC and OP. In contrast, FE-derived FL may help identify patients with increased osteoporotic fracture risk in the future.

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/cm3 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.

Multi-detector computed tomography (MDCT) imaging: association of bone texture parameters with finite element analysis (FEA)-based failure load of single vertebrae and functional spinal units.

BACKGROUND: Osteoporosis is a systemic skeletal disease that is characterized by low bone mass and microarchitectural deterioration, predisposing affected individuals to fragility fractures. Yet, standard measurement of areal bone mineral density (BMD) in dual-energy X-ray absorptiometry (DXA) as the current reference standard has limitations for correctly detecting osteoporosis and fracture risk, with opportunistic osteoporosis screening using computed tomography (CT) showing increasing importance. This study's objective is to compare finite element analysis (FEA)-based vertebral failure load with parameters of texture analysis (TA) derived from multi-detector CT (MDCT). METHODS: MDCT data of seven subjects (mean age: 71.9±7.4 years) were included for FEA and TA. Manual segmentation was performed for the vertebral bodies T11, T12, L1, and L2 and the intervertebral discs (IVDs) T11/12, T12/L1, L1/2, and L2/3. Correlation analyses between FEA-derived failure loads and parameters of TA for the single vertebrae and two functional spinal units (FSUs) were calculated, defining FSU-1 as T11-IVD-T12-IVD-L1 and FSU-2 as T12-IVD-L1-IVD-L2. Furthermore, multivariate regressions were performed to identify the texture parameters that predicted the failure load best. RESULTS: For single vertebrae, the strongest correlations were observed for skewnessglobal, kurtosisglobal, and gray level variance (rho =-0.7668 to -0.7362; P<0.001), while for FSUs, SumAverage, long-run emphasis, long-run low gray-level emphasis, homogeneity, and energy showed the strongest correlations (rho =-0.8187 to 0.8407; P<0.05) to failure loads. SumAverage best predicted the failure load for single vertebrae (R2 adj =0.523, P<0.001). For the two FSUs, kurtosisglobal (FSU-1: R2 adj =0.611, P=0.001) and skewnessglobal (FSU-2: R2 adj =0.579, P=0.002) were the best predictors. CONCLUSIONS: TA using MDCT data of the spine was significantly associated with FEA-derived failure loads of both, single vertebrae and FSUs. Texture parameters predicted failure loads of FSUs as a more realistic in-vivo scenario equally well as compared to single vertebrae analyses. TA may reflect a less complex and time-consuming approach to accurately and non-invasively evaluate vertebral bone strength.

Prediction of incident vertebral fractures in routine MDCT: Comparison of global texture features, 3D finite element parameters and volumetric BMD.

PURPOSE: In this case-control study, we evaluated different quantitative parameters derived from routine multi-detector computed tomography (MDCT) scans with respect to their ability to predict incident osteoporotic vertebral fractures of the thoracolumbar spine. METHODS: 16 patients who received baseline and follow-up contrast-enhanced MDCT and were diagnosed with an incident osteoporotic vertebral fracture at follow-up, and 16 age-, sex-, and follow-up-time-matched controls were included in the study. Vertebrae were labelled and segmented using a fully automated pipeline. Volumetric bone mineral density (vBMD), finite element analysis (FEA)-based failure load (FL) and failure displacement (FD), as well as 24 texture features were extracted from L1 - L3 and averaged. Odds ratios (OR) with 95% confidence intervals (CI), expressed per standard deviation decrease, receiver operating characteristic (ROC) area under the curve (AUC), as well as logistic regression models, including all analyzed parameters as independent variables, were used to assess the prediction of incident vertebral fractures. RESULTS: The texture feature Correlation (AUC = 0.754, p = 0.014; OR = 2.76, CI = 1.16-6.58) and vBMD (AUC = 0.750, p = 0.016; OR = 2.67, CI = 1.12-6.37) classified incident vertebral fractures best, while the best FEA-based parameter FL showed an AUC = 0.719 (p = 0.035). Correlation was the only significant predictor of incident fractures in the logistic regression analysis of all parameters (p = 0.022). CONCLUSION: MDCT-derived FEA parameters and texture features, averaged from L1 - L3, showed only a moderate, but no statistically significant improvement of incident vertebral fracture prediction beyond BMD, supporting the hypothesis that vertebral-specific parameters may be superior for fracture risk assessment.

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.

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 (R2), 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 (R2 = 0.93, RMSCV = 8.6% for 50%; R2 = 0.89, RMSCV = 11.90% for 75%; R2 = 0.86, RMSCV = 11.30% for 90%) and up to 50% dose reduction through tube current reduction method (R2 = 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 (R2 = 0.88, RMSCV = 22.04% for 75%; R2 = 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.

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 (R2) 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%, R2 = 0.96 vs. RMSCV = 20.78%, R2 = 0.79). Seventy-five percent dose reduction achieved through both methods (RMSCV = 22.38%, R2 = 0.80 for sparse sampling; RMSCV = 24.58%, R2 = 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.

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.