Computed tomography-based finite element analysis to assess fracture risk and osteoporosis treatment

Reference value of testicular temperature measured by finite element analysis after first staged inguinal orchidopexy in children with abdominal testis and short spermatic cord

PURPOSE

This study aims to build a 3D reconstruction computed simulation model and to establish a regression equation for detecting the testis's temperature by its location after first staged open orchidopexy in children with abdominal undescended testis (UDT) and short spermatic cords.

METHODS

In this cross-sectional study, we enrolled 31 children with abdominal UDT and short spermatic cords who underwent first staged orchiopexy between 2017 and 2020. Using ultrasonography to obtain the testis's location distance from the skin surface (X1), external iliac vessel (X2), and internal inguinal ring (X3), we input the data into a 3D reconstruction computed simulation along with COMSOL to calculate the testicular temperature. We also used multivariate regression to establish the testicular temperature regression equation from the gathered data.

RESULT

The mean age of the participants was 4.47 ± 1.21 years. The mean size of the operated testis was 0.39 ± 0.13 cc. The mean distance of the testis from X1, X2, and X3 was 3.27 ± 1.25 mm, 21.06 ± 6.42 mm, and 27.19 ± 10.09 mm, respectively. The testicular temperature regression equation derived from testis location was calculated by the formula: 34.57 + 0.0236 X12 - 0.0105 X2 - 0.0018 X3. The concordance for testis temperature calculated via the computational method and regression equation was 83%.

CONCLUSIONS

The current study provided a reference value for the testicular temperature of children with abdominal UDT and short spermatic cords after the first stage of orchiopexy. A testicular temperature regression equation can be established based on the testis location, which will provide relevant information for the testicular development assessment, disease diagnosis, and follow-up, and possibly determination of the time of the second stage of orchiopexy.

Cost-Effectiveness Analysis of CT-Based Finite Element Modeling for Osteoporosis Screening in Secondary Fracture Prevention: An Early Health Technology Assessment in the Netherlands

Objective. To conduct cost-utility analyses for Computed Tomography To Strength (CT2S), a novel osteoporosis screening service, compared with dual-energy X-ray absorptiometry (DXA), treat all without screening, and no screening methods for Dutch postmenopausal women referred to fracture liaison service (FLS). CT2S uses CT scans to generate femur models and simulate sideways fall scenarios for bone strength assessment. Methods. Early health technology assessment (HTA) was adopted to evaluate CT2S as a novel osteoporosis screening tool for secondary fracture prevention. We constructed a 2-dimensional simulation model considering 4 strategies (no screening, treat all without screening, DXA, CT2S) together with screening intervals (5 y, 2 y), treatments (oral alendronate, zoledronic acid), and discount rate scenarios among Dutch women in 3 age groups (60s, 70s, and 80s). Strategy comparisons were based on incremental cost-effectiveness ratios (ICERs), considering an ICER below €20,000 per QALY gained as cost-effective in the Netherlands. Results. Under the base-case scenario, CT2S versus DXA had estimated ICERs of €41,200 and €14,083 per QALY gained for the 60s and 70s age groups, respectively. For the 80s age group, CT2S was more effective and less costly than DXA. Changing treatment from weekly oral alendronate to annual zoledronic acid substantially decreased CT2S versus DXA ICERs across all age groups. Setting the screening interval to 2 y increased CT2S versus DXA ICERs to €100,333, €55,571, and €15,750 per QALY gained for the 60s, 70s, and 80s age groups, respectively. In all simulated populations and scenarios, CT2S was cost-effective (in some cases dominant) compared with the treat all strategy and cost-saving (more effective and less costly) compared with no screening. Conclusion. CT2S was estimated to be potentially cost-effective in the 70s and 80s age groups considering the willingness-to-pay threshold of the Netherlands. This early HTA suggests CT2S as a potential novel osteoporosis screening tool for secondary fracture prevention.

Highlights

For postmenopausal Dutch women who have been referred to the FLS, direct access to CT2S may be cost-effective compared with DXA for age groups 70s and 80s, when considering the ICER threshold of the Netherlands. This study positions CT2S as a potential novel osteoporosis-screening tool for secondary fracture prevention in the clinical setting.A shorter screening interval of 2 y increases the effectiveness of both screening strategies, but the ICER of CT2S compared with DXA also increased substantially, which made CT2S no longer cost-effective for the 70s age group; however, it remains cost-effective for individuals in their 80s.Annual zoledronic acid treatment with better adherence may contribute to a lower cost-effectiveness ratio when comparing CT2S to DXA screening and the treat all strategies for all age groups.

Analysis of stress and stabilization in adolescent with osteoporotic idiopathic scoliosis: finite element method

Objective: To explore the effect of osteoporosis on the stress, stability, and lumbar intervertebral disc of AIS lumbar vertebrae by finite element method. Better understand the biomechanical characteristics of osteoporotic scoliosis.Methods: Based on the CT images of normal lumbar vertebrae and lumbar vertebrae with AIS, the finite element models were established to simulate the estimated osteoporosis by changing the Young's modulus of cortical bone, cancellous bone, and endplate. Four finite element models of normal lumbar, osteoporotic lumbar, normal AIS lumbar and osteoporotic AIS lumbar were established, and the same load and boundary conditions were applied respectively. The displacement, stress, and intervertebral disc strain of the four models were compared to explore the effect of osteoporosis on the stability and injury risk of AIS.Results: After suffering from osteoporosis, under the same load, the displacement of lumbar spine increases, the stability decreases, and the stability of AIS lumbar spine decrease more obviously, especially under extension load. Suffering from osteoporosis will increase the stress of lumbar spine, AIS lumbar spine increases more obviously, and the stress is more concentrated, Osteoporotic lumbar spine only affects the strain of intervertebral disc when AIS lumbar spine bends on the concave side, resulting in greater strain behind the concave side of intervertebral disc.Conclusions: AIS patients with OP have lower lumbar stability, a higher risk of fracture of lumbar vertebrae, and spinal nerves are more likely to be compressed by intervertebral discs. OP can aggravate the scoliosis of lumbar vertebrae.

The Effect of Concave-Side Intertransverse Ligament Laxity on the Stress of AIS Lumbar Spine Based on Finite Element Method

(1) Background: Scoliosis has the mechanical characteristic of asymmetric stress distribution, which is one of the reasons for the aggravation of scoliosis. Bracing therapy is the best treatment for AIS, but it is difficult and costly to operate. Is it possible to reduce pressure in the concave side by relaxing the ITL in the concave side of scoliosis, so as to improve the abnormal stress distribution of scoliosis? In this paper, a finite element method was used to simulate the effect of the relaxation of concave-side ITL on the stress of a lumbar spine with scoliosis, which provides some guidance for the treatment of scoliosis. (2) Methods: Using CT images of a patient with scoliosis whose Cobb Angle was 43° and Lordosis Angle was 45, a scoliosis lumbar was established, and Young's modulus of the ITL of the concave-side lumbar spine was reduced by 95% to simulate ligament relaxation. By comparing the stress condition of the model vertebral body with no ligament relaxation, the effect of concave-side ITL relaxation on the mechanical characteristics of scoliosis lumbar spine was explored. (3) Results: An effective and complete model of the lumbar spine was established. The concave ITL relaxed, which only had a great impact on the bending loads. After the ligament was relaxed, the stability of the spine was reduced. Stress concentration on the concave side of vertebrae and the IVD was aggravated. Under loads on the convex side, the maximum stress on the vertebral body and the IVD increased significantly, making lumbar vertebrae more vulnerable to injury. (4) Conclusions: Laxity of the ITL on the concave side of the AIS lumbar only affects the bending load. Laxity of the concave-side ligament will reduce the stability of the lumbar, aggravate the uneven stress distribution of scoliotic lumbar vertebrae, increase the risk of IVD injury, and be unfavorable for the scoliotic lumbar spine. Relaxation of the concave ITL alone is not an effective way to treat scoliosis.

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.

Optimization of Spinal Reconstructions for Thoracolumbar Burst Fractures to Prevent Proximal Junctional Complications: A Finite Element Study

The management strategies of thoracolumbar (TL) burst fractures include posterior, anterior, and combined approaches. However, the rigid constructs pose a risk of proximal junctional failure. In this study, we aim to systemically evaluate the biomechanical performance of different TL reconstruction constructs using finite element analysis. Furthermore, we investigate the motion and the stress on the proximal junctional level adjacent to the constructs. We used a T10-L3 finite element model and simulated L1 burst fracture. Reconstruction with posterior instrumentation (PI) alone (U2L2 and U1L1+(intermediate screw) and three-column spinal reconstruction (TCSR) constructs (U1L1+PMMA and U1L1+Cage) were compared. Long-segment PI resulted in greater global motion reduction compared to constructs with short-segment PI. TCSR constructs provided better stabilization in L1 compared to PI alone. Decreased intradiscal and intravertebral pressure in the proximal level were observed in U1L1+IS, U1L1+PMMA, and U1L1+Cage compared to U2L2. The stress and strain energy of the pedicle screws decreased when anterior reconstruction was performed in addition to PI. We showed that TCSR with anterior reconstruction and SSPI provided sufficient immobilization while offering additional advantages in the preservation of physiological motion, the decreased burden on the proximal junctional level, and lower risk of implant failure.

Comparison of the susceptibility to implant failure in the lateral, posterior, and transforaminal lumbar interbody fusion: A finite element analysis

OBJECTIVE

There are several techniques for lumbar interbody fusion, and implant failure following lumbar interbody fusion can be troublesome. This study aimed to compare the stress in posterior implant and peri-screw vertebral bodies among lateral lumbar interbody fusion (LLIF), posterior lumbar interbody fusion (PLIF), and transforaminal lumbar interbody fusion (TLIF) and to select the technique that is least likely to cause implant failure.

METHODS

We created an intact L3-L5 model and simulated the LLIF, PLIF, and TLIF techniques at L4-L5 using finite element methods. All models at the lower portion of L5 were fixed and imposed a preload of 400 N and a moment of 7.5 Nm on the upper portion of L3 to simulate flexion, extension, lateral bending, and axial rotation. We investigated the peak stresses and stress concentration in the posterior implant and peri-screw vertebral bodies for the LLIF, PLIF, and TLIF techniques.

RESULTS

The extension, flexion, bending, and rotation peak stresses and stress concentration in the posterior implant, and the peri-screw vertebral bodies, were the lowest in LLIF, followed by PLIF and TLIF, respectively.

CONCLUSIONS

It was found that implant failure was least likely to occur in LLIF, followed by PLIF and TLIF, respectively. Hence, surgeons should be aware of these factors when selecting an appropriate surgical technique and be careful for implant failure during postoperative follow-up.

Biomechanical Study of the Osteoporotic Spine Fracture: Optical Approach

Background and objectives: Osteoporotic spine fractures represent a significant factor for decreasing quality of life in the elderly female population. Understanding the mechanisms involved in producing these fractures can improve their prevention and treatment. This study presents a biomechanical method to produce a vertebral fracture, conducted on a human spine segment, observing the displacements and strains in the intervertebral disc, endplate, and vertebral body. Materials and Methods: We performed two tests, one corresponding to an extension loading, and the second to an axial loading. Results: The maximum displacement in the target vertebral body presented higher values in the case of the extension as compared to the axial strain where it mainly occurred after the fracture was produced. The strains occurred simultaneously on both discs. In the case of the axial strain, due to the occurrence of the fracture, the maximum value was recorded in the spine body, while in the case of the extensions, it occurred in the neural part of the upper disc. The advantage of this method was that the entire study was an experiment, using optical methods, increasing the precision of the material data input. Conclusions: The research method allowed recording in real time of a larger amount of data from the different components of the spine segment. If there was an extension component of the compression force at the moment of the initial loading, part of this load was absorbed by the posterior column with higher mechanical resistance. After the maximum capacity of the absorption was reached, in both situations the behavior was similar.

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.

Testosterone therapy and bone quality in men with diabetes and hypogonadism: Study design and protocol

Context

Type 2 diabetes mellitus (T2D) is often accompanied by male hypogonadism and both conditions are associated with increased risk for fractures. Testosterone (T) has been shown to improve the bone health of hypogonadal men but has not been tested in patients who also have T2D in addition to low T. To date, there is no treatment that is specifically recommended for bone disease among patients with T2D. This study will evaluate the effect of T therapy on the bone health of male veterans with low T who also have T2D.

Methods

This is a randomized double-blind placebo-controlled trial of 166 male veterans 35–65 years old, with T2D and hypogonadism, randomized to either T gel 1.62% or placebo for 12 months. We will evaluate the effect of T therapy on the following primary outcomes:1) changes in bone strength as measured by microfinite elements analysis (μFEA) using high-resolution peripheral quantitative computer tomography, 2) changes in bone turnover markers, and 3) changes in circulating osteoblast progenitors (COP) and osteoclast precursors cells.

Discussion

We anticipate that T therapy will result in improvement in bone strength owing to improvement in bone remodeling through an increase in osteoblastic differentiation and proliferation in patients with hypogonadism and T2D.

Computer-aided diagnosis systems for osteoporosis detection: a comprehensive survey

Computer-aided diagnosis (CAD) has revolutionized the field of medical diagnosis. They assist in improving the treatment potentials and intensify the survival frequency by early diagnosing the diseases in an efficient, timely, and cost-effective way. The automatic segmentation has led the radiologist to successfully segment the region of interest to improve the diagnosis of diseases from medical images which is not so efficiently possible by manual segmentation. The aim of this paper is to survey the vision-based CAD systems especially focusing on the segmentation techniques for the pathological bone disease known as osteoporosis. Osteoporosis is the state of the bones where the mineral density of bones decreases and they become porous, making the bones easily susceptible to fractures by small injury or a fall. The article covers the image acquisition techniques for acquiring the medical images for osteoporosis diagnosis. The article also discusses the advanced machine learning paradigms employed in segmentation for osteoporosis disease. Other image processing steps in osteoporosis like feature extraction and classification are also briefly described. Finally, the paper gives the future directions to improve the osteoporosis diagnosis and presents the proposed architecture. Graphical abstract.

Neuro-musculoskeletal flexible multibody simulation yields a framework for efficient bone failure risk assessment

Fragility fractures are a major socioeconomic problem. A non-invasive, computationally-efficient method for the identification of fracture risk scenarios under the representation of neuro-musculoskeletal dynamics does not exist. We introduce a computational workflow that integrates modally-reduced, quantitative CT-based finite-element models into neuro-musculoskeletal flexible multibody simulation (NfMBS) for early bone fracture risk assessment. Our workflow quantifies the bone strength via the osteogenic stresses and strains that arise due to the physiological-like loading of the bone under the representation of patient-specific neuro-musculoskeletal dynamics. This allows for non-invasive, computationally-efficient dynamic analysis over the enormous parameter space of fracture risk scenarios, while requiring only sparse clinical data. Experimental validation on a fresh human femur specimen together with femur strength computations that were consistent with literature findings provide confidence in the workflow: The simulation of an entire squat took only 38 s CPU-time. Owing to the loss (16% cortical, 33% trabecular) of bone mineral density (BMD), the strain measure that is associated with bone fracture increased by 31.4%; and yielded an elevated risk of a femoral hip fracture. Our novel workflow could offer clinicians with decision-making guidance by enabling the first combined in-silico analysis tool using NfMBS and BMD measurements for optimized bone fracture risk assessment.

The use of subject-specific Finite Element analysis of L1-L4 vertebra to screening osteoporosis in postmenopausal women

The aim of this work is to study Computed Tomography (CT)-based Finite Element (FE) modelling of lumbar vertebra to classify osteoporotic patients using dual-energy x-ray absorptiometry (DXA) as reference. Cohort comprised 15 postmenopausal female patients with CT-DXA pairs. Sixty subject-specific CT-based FE models were constructed of L1 to L4 vertebrae. Load-displacement in each node of the models was analyzed. Mean kurtosis values of displacement distributions between osteoporosis and healthy groups at L1 level showed a statistically significant difference (p<;0.005). The CT-based FE model of L1 vertebra with patient specific material properties can be used to identify osteoporosis patients.

Advances in imaging approaches to fracture risk evaluation

Fragility fractures are a growing problem worldwide, and current methods for diagnosing osteoporosis do not always identify individuals who require treatment to prevent a fracture and may misidentify those not a risk. Traditionally, fracture risk is assessed using dual-energy X-ray absorptiometry, which provides measurements of areal bone mineral density at sites prone to fracture. Recent advances in imaging show promise in adding new information that could improve the prediction of fracture risk in the clinic. As reviewed herein, advances in quantitative computed tomography (QCT) predict hip and vertebral body strength; high-resolution HR-peripheral QCT (HR-pQCT) and micromagnetic resonance imaging assess the microarchitecture of trabecular bone; quantitative ultrasound measures the modulus or tissue stiffness of cortical bone; and quantitative ultrashort echo-time MRI methods quantify the concentrations of bound water and pore water in cortical bone, which reflect a variety of mechanical properties of bone. Each of these technologies provides unique characteristics of bone and may improve fracture risk diagnoses and reduce prevalence of fractures by helping to guide treatment decisions.

Effects of dose reduction on bone strength prediction using finite element analysis

This study aimed to evaluate the effect of dose reduction, by means of tube exposure reduction, on bone strength prediction from finite-element (FE) analysis. Fresh thoracic mid-vertebrae specimens (n = 11) were imaged, using multi-detector computed tomography (MDCT), at different intensities of X-ray tube exposures (80, 150, 220 and 500 mAs). Bone mineral density (BMD) was estimated from the mid-slice of each specimen from MDCT images. Differences in image quality and geometry of each specimen were measured. FE analysis was performed on all specimens to predict fracture load. Paired t-tests were used to compare the results obtained, using the highest CT dose (500 mAs) as reference. Dose reduction had no significant impact on FE-predicted fracture loads, with significant correlations obtained with reference to 500 mAs, for 80 mAs (R²  = 0.997, p < 0.001), 150 mAs (R² = 0.998, p < 0.001) and 220 mAs (R² = 0.987, p < 0.001). There were no significant differences in volume quantification between the different doses examined. CT imaging radiation dose could be reduced substantially to 64% with no impact on strength estimates obtained from FE analysis. Reduced CT dose will enable early diagnosis and advanced monitoring of osteoporosis and associated fracture risk.