Effect of the intervertebral disc on vertebral bone strength prediction: a finite-element study

Geometric determinants of the mechanical behavior of image-based finite element models of the intervertebral disc

The intervertebral disc is an important structure for load transfer through the spine. Its injury and degeneration have been linked to pain and spinal fractures. Disc injury and spine fractures are associated with high stresses; however, these stresses cannot be measured, necessitating the use of finite element (FE) models. These models should include the disc's complex structure, as changes in disc geometry have been linked to altered mechanical behavior. However, image-based models using disc-specific structures have yet to be established. This study describes a multiphasic FE modeling approach for noninvasive estimates of subject-specific intervertebral disc mechanical behavior based on medical imaging. The models (n = 22) were used to study the influence of disc geometry on the predicted global mechanical response (moments and forces), internal local disc stresses, and tractions at the interface between the disc and the bone. Disc geometry was found to have a strong influence on the predicted moments and forces on the disc (R2 = 0.69-0.93), while assumptions regarding the side curvature (bulge) of the disc had only a minor effect. Strong variability in the predicted internal disc stresses and tractions was observed between the models (mean absolute differences of 5.1%-27.7%). Disc height had a systematic influence on the internal disc stresses and tractions at the disc-to-bone interface. The influence of disc geometry on mechanics highlights the importance of disc-specific modeling to estimate disc injury risk, loading on the adjacent vertebral bodies, and the mechanical environment present in disc tissues.

Preoperative Prediction of New Vertebral Fractures after Vertebral Augmentation with a Radiomics Nomogram

The occurrence of new vertebral fractures (NVFs) after vertebral augmentation (VA) procedures is common in patients with osteoporotic vertebral compression fractures (OVCFs), leading to painful experiences and financial burdens. We aim to develop a radiomics nomogram for the preoperative prediction of NVFs after VA. Data from center 1 (training set: n = 153; internal validation set: n = 66) and center 2 (external validation set: n = 44) were retrospectively collected. Radiomics features were extracted from MRI images and radiomics scores (radscores) were constructed for each level-specific vertebra based on least absolute shrinkage and selection operator (LASSO). The radiomics nomogram, integrating radiomics signature with presence of intravertebral cleft and number of previous vertebral fractures, was developed by multivariable logistic regression analysis. The predictive performance of the vertebrae was level-specific based on radscores and was generally superior to clinical variables. RadscoreL2 had the optimal discrimination (AUC ≥ 0.751). The nomogram provided good predictive performance (AUC ≥ 0.834), favorable calibration, and large clinical net benefits in each set. It was used successfully to categorize patients into high- or low-risk subgroups. As a noninvasive preoperative prediction tool, the MRI-based radiomics nomogram holds great promise for individualized prediction of NVFs following VA.

CT-Based Automatic Spine Segmentation Using Patch-Based Deep Learning

CT vertebral segmentation plays an essential role in various clinical applications, such as computer-assisted surgical interventions, assessment of spinal abnormalities, and vertebral compression fractures. Automatic CT vertebral segmentation is challenging due to the overlapping shadows of thoracoabdominal structures such as the lungs, bony structures such as the ribs, and other issues such as ambiguous object borders, complicated spine architecture, patient variability, and fluctuations in image contrast. Deep learning is an emerging technique for disease diagnosis in the medical field. This study proposes a patch-based deep learning approach to extract the discriminative features from unlabeled data using a stacked sparse autoencoder (SSAE). 2D slices from a CT volume are divided into overlapping patches fed into the model for training. A random under sampling (RUS)-module is applied to balance the training data by selecting a subset of the majority class. SSAE uses pixel intensities alone to learn high-level features to recognize distinctive features from image patches. Each image is subjected to a sliding window operation to express image patches using autoencoder high-level features, which are then fed into a sigmoid layer to classify whether each patch is a vertebra or not. We validate our approach on three diverse publicly available datasets: VerSe, CSI-Seg, and the Lumbar CT dataset. Our proposed method outperformed other models after configuration optimization by achieving 89.9% in precision, 90.2% in recall, 98.9% in accuracy, 90.4% in F-score, 82.6% in intersection over union (IoU), and 90.2% in Dice coefficient (DC). The results of this study demonstrate that our model’s performance consistency using a variety of validation strategies is flexible, fast, and generalizable, making it suited for clinical application.

Development and model form assessment of an automatic subject-specific vertebra reconstruction method

BACKGROUND

Current spine models for analog bench models, surgical navigation and training platforms are conventionally based on 3D models from anatomical human body polygon database or from time-consuming manual-labelled data. This work proposed a workflow of quick and accurate subject-specific vertebra reconstruction method and quantified the reconstructed model accuracy and model form errors.

METHODS

Four different neural networks were customized for vertebra segmentation. To validate the workflow in clinical applications, an excised human lumbar vertebra was scanned via CT and reconstructed into 3D CAD models using four refined networks. A reverse engineering solution was proposed to obtain the high-precision geometry of the excised vertebra as gold standard. The 3D model evaluation metrics and a finite element analysis (FEA) method were designed to reflect the model accuracy and model form errors.

RESULTS

The automatic segmentation networks achieved the best Dice score of 94.20% in validation datasets. The accuracy of reconstructed models was quantified with the best 3D Dice index of 92.80%, 3D IoU of 86.56%, Hausdorff distance of 1.60 mm, and the heatmaps and histograms were used for error visualization. The FEA results showed the impact of different geometries and reflected partial surface accuracy of the reconstructed vertebra under biomechanical loads with the closest percentage error of 4.2710% compared to the gold standard model.

CONCLUSIONS

In this work, a workflow of automatic subject-specific vertebra reconstruction method was proposed while the errors in geometry and FEA were quantified. Such errors should be considered when leveraging subject-specific modelling towards the development and improvement of treatments.

A Review of CT-Based Fracture Risk Assessment with Finite Element Modeling and Machine Learning

PURPOSE OF REVIEW

We reviewed advances over the past 3 years in assessment of fracture risk based on CT scans, considering methods that use finite element models, machine learning, or a combination of both.

RECENT FINDINGS

Several studies have demonstrated that CT-based assessment of fracture risk, using finite element modeling or biomarkers derived from machine learning, is equivalent to currently used clinical tools. Phantomless calibration of CT scans for bone mineral density enables accurate measurements from routinely taken scans. This opportunistic use of CT scans for fracture risk assessment is facilitated by high-quality automated segmentation with deep learning, enabling workflows that do not require user intervention. Modeling of more realistic and diverse loading conditions, as well as improved modeling of fracture mechanisms, has shown promise to enhance our understanding of fracture processes and improve the assessment of fracture risk beyond the performance of current clinical tools. CT-based screening for fracture risk is effective and, by analyzing scans that were taken for other indications, could be used to expand the pool of people screened, therefore improving fracture prevention. Finite element modeling and machine learning both provide valuable tools for fracture risk assessment. Future approaches should focus on including more loading-related aspects of fracture risk.

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.

A Three-Dimensional Cement Quantification Method for Decision Prediction of Vertebral Recompression after Vertebroplasty

Objective

Proposing parameters to quantify cement distribution and increasing accuracy for decision prediction of vertebroplasty postoperative complication.

Methods

Finite element analysis was used to biomechanically assess vertebral mechanics (n = 51) after percutaneous vertebroplasty (PVP) or kyphoplasty (PKP). The vertebral space was divided into 27 portions. The numbers of cement occupied portions and numbers of cement-endplate contact portions were defined as overall distribution number (oDN) and overall endplate contact number (oEP), respectively. And cement distribution was parametrized by oDN and oEP. The determination coefficients of vertebral mechanics and parameters (R 2) can validate the correlation of proposed parameters with vertebral mechanics.

Results

oDN and oEP were mainly correlated with failure load (R 2 = 0.729) and stiffness (R 2 = 0.684), respectively. oDN, oEP, failure load, and stiffness had obvious difference between the PVP group and the PKP group (P < 0.05). The regional endplate contact number in the front column is most correlated with vertebral stiffness (R 2 = 0.59) among all regional parameters. Cement volume and volume fraction are not dominant factors of vertebral augmentation, and they are not suitable for postoperative fracture risk prediction.

Conclusions

Proposed parameters with high correlation on vertebral mechanics are promising for clinical utility. The oDN and oEP can strongly affect augmented vertebral mechanics thus is suitable for postoperative fracture risk prediction. The parameters are beneficial for decision-making process of revision surgery necessity. Parametrized methods are also favorable for surgeon's preoperative planning. The methods can be inspirational for clinical image recognition development and auxiliary diagnosis.

Imaging of the Osteoporotic Spine - Quantitative Approaches in Diagnostics and for the Prediction of the Individual Fracture Risk

Osteoporosis is a highly prevalent systemic skeletal disease that is characterized by low bone mass and microarchitectural bone deterioration. It predisposes to fragility fractures that can occur at various sites of the skeleton, but vertebral fractures (VFs) have been shown to be particularly common. Prevention strategies and timely intervention depend on reliable diagnosis and prediction of the individual fracture risk, and dual-energy X-ray absorptiometry (DXA) has been the reference standard for decades. Yet, DXA has its inherent limitations, and other techniques have shown potential as viable add-on or even stand-alone options. Specifically, three-dimensional (3 D) imaging modalities, such as computed tomography (CT) and magnetic resonance imaging (MRI), are playing an increasing role. For CT, recent advances in medical image analysis now allow automatic vertebral segmentation and value extraction from single vertebral bodies using a deep-learning-based architecture that can be implemented in clinical practice. Regarding MRI, a variety of methods have been developed over recent years, including magnetic resonance spectroscopy (MRS) and chemical shift encoding-based water-fat MRI (CSE-MRI) that enable the extraction of a vertebral body's proton density fat fraction (PDFF) as a promising surrogate biomarker of bone health. Yet, imaging data from CT or MRI may be more efficiently used when combined with advanced analysis techniques such as texture analysis (TA; to provide spatially resolved assessments of vertebral body composition) or finite element analysis (FEA; to provide estimates of bone strength) to further improve fracture prediction. However, distinct and experimentally validated diagnostic criteria for osteoporosis based on CT- and MRI-derived measures have not yet been achieved, limiting broad transfer to clinical practice for these novel approaches. KEY POINTS::   · DXA is the reference standard for diagnosis and fracture prediction in osteoporosis, but it has important limitations.. · CT- and MRI-based methods are increasingly used as (opportunistic) approaches.. · For CT, particularly deep-learning-based automatic vertebral segmentation and value extraction seem promising.. · For MRI, multiple techniques including spectroscopy and chemical shift imaging are available to extract fat fractions.. · Texture and finite element analyses can provide additional measures for vertebral body composition and bone strength.. CITATION FORMAT: · Sollmann N, Kirschke JS, Kronthaler S et al. Imaging of the Osteoporotic Spine - Quantitative Approaches in Diagnostics and for the Prediction of the Individual Fracture Risk. Fortschr Röntgenstr 2022; DOI: 10.1055/a-1770-4626.

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.

Finite Element Method for the Evaluation of the Human Spine: A Literature Overview

The finite element method (FEM) represents a computer simulation method, originally used in civil engineering, which dates back to the early 1940s. Applications of FEM have also been used in numerous medical areas and in orthopedic surgery. Computing technology has improved over the years and as a result, more complex problems, such as those involving the spine, can be analyzed. The spine is a complex anatomical structure that maintains the erect posture and supports considerable loads. Applications of FEM in the spine have contributed to the understanding of bone biomechanics, both in healthy and abnormal conditions, such as scoliosis, fractures (trauma), degenerative disc disease and osteoporosis. However, since FEM is only a digital simulation of the real condition, it will never exactly simulate in vivo results. In particular, when it concerns biomechanics, there are many features that are difficult to represent in a FEM. More FEM studies and spine research are required in order to examine interpersonal spine stiffness, young spine biomechanics and model accuracy. In the future, patient-specific models will be used for better patient evaluations as well as for better pre- and inter-operative planning.

VerSe: A Vertebrae labelling and segmentation benchmark for multi-detector CT images

Vertebral labelling and segmentation are two fundamental tasks in an automated spine processing pipeline. Reliable and accurate processing of spine images is expected to benefit clinical decision support systems for diagnosis, surgery planning, and population-based analysis of spine and bone health. However, designing automated algorithms for spine processing is challenging predominantly due to considerable variations in anatomy and acquisition protocols and due to a severe shortage of publicly available data. Addressing these limitations, the Large Scale Vertebrae Segmentation Challenge (VerSe) was organised in conjunction with the International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI) in 2019 and 2020, with a call for algorithms tackling the labelling and segmentation of vertebrae. Two datasets containing a total of 374 multi-detector CT scans from 355 patients were prepared and 4505 vertebrae have individually been annotated at voxel level by a human-machine hybrid algorithm (https://osf.io/nqjyw/, https://osf.io/t98fz/). A total of 25 algorithms were benchmarked on these datasets. In this work, we present the results of this evaluation and further investigate the performance variation at the vertebra level, scan level, and different fields of view. We also evaluate the generalisability of the approaches to an implicit domain shift in data by evaluating the top-performing algorithms of one challenge iteration on data from the other iteration. The principal takeaway from VerSe: the performance of an algorithm in labelling and segmenting a spine scan hinges on its ability to correctly identify vertebrae in cases of rare anatomical variations. The VerSe content and code can be accessed at: https://github.com/anjany/verse.

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