Age and gender effects on bone mass density variation: finite elements simulation

Bone density of the cervical, thoracic and lumbar spine measured using Hounsfield units of computed tomography - results of 4350 vertebras

INTRODUCTION

The assessment of bone density has gained significance in recent years due to the aging population. Accurate assessment of bone density is crucial when deciding on the appropriate treatment plan for spinal stabilization surgery. The objective of this work was to determine the trabecular bone density values of the subaxial cervical, thoracic and lumbar spine using Hounsfield units.

MATERIAL AND METHODS

Data from 200 patients who underwent contrast-enhanced polytrauma computed tomography at a maximum care hospital over a two-year period were retrospectively analyzed. HUs were measured with an elliptical measurement field in three different locations within the vertebral body: below the upper plate, in the middle of the vertebral body, and above the base plate. The measured Hounsfield units were converted into bone density values using a validated formula.

RESULTS

The mean age of the patient collective was 47.05 years. Mean spinal bone density values decreased from cranial to caudal (C3: 231.79 mg/cm3; L5: 155.13 mg/cm3; p < 0.001), with the highest values in the upper cervical spine. Bone density values generally decreased with age in all spinal segments. There was a clear decrease in values after age 50 years (p < 0.001).

CONCLUSIONS

In our study, bone density decreased from cranial to caudal with higher values in the cervical spine. These data from the individual spinal segments may be helpful to comprehensively evaluate the status of the spine and to design a better preoperative plan before instrumentation.

Modeling and simulation of bone cells dynamic behavior under the late effect of breast cancer treatments

Breast Cancer (BC) treatments have been proven to interfere with the health of bones. Chemotherapy and endocrinal treatment regimens such as tamoxifen and aromatase inhibitors are frequently prescribed for women with BC. However, these drugs increase bone resorption and reduce the Bone Mineral Density (BMD), thus increasing the risk of bone fracture. In the current study, a mechanobiological bone remodeling model has been developed by coupling cellular activities, mechanical stimuli, and the effect of breast cancer treatments (chemotherapy, tamoxifen, and aromatase inhibitors). This model algorithm has been programmed and implemented on MATLAB software to simulate different treatment scenarios and their effects on bone remodeling and also predict the evolution of Bone Volume fraction (BV/TV) and the associated Bone Density Loss (BDL) over a period of time. The simulation results, achieved from different combinations of Breast Cancer treatments, allow the researchers to predict the intensity of each combination treatment on BV/TV and BMD. The combination of chemotherapy, tamoxifen, and aromatase inhibitors, followed by the combination of chemotherapy and tamoxifen remain the most harmful regimen. This is because they have a strong ability to induce the bone degradation which is represented by a decrease of 13.55% and 11.55% of the BV/TV value, respectively. These results were compared with the experimental studies and clinical observations which showed good agreement. The proposed model can be used by clinicians and physicians to choose the most appropriate combination of treatments, according to the patient's case.

A coupled mathematical model between bone remodeling and tumors: a study of different scenarios using Komarova's model

This paper aims to construct a general framework of coupling tumor-bone remodeling processes in order to produce plausible outcomes of the effects of tumors on the number of osteoclasts, osteoblasts, and the frequency of the bone turnover cycle. In this document, Komarova's model has been extended to include the effect of tumors on the bone remodeling processes. Thus, we explored three alternatives for coupling tumor presence into Komarova's model: first, using a "damage" parameter that depends on the tumor cell concentration. A second model follows the original structure of Komarova, including the tumor presence in those equations powered up to a new parameter, called the paracrine effect of the tumor on osteoclasts and osteoblasts; the last model is replicated from Ayati and collaborators in which the impact of the tumor is included into the paracrine parameters. Through the models, we studied their stability and considered some examples that can reproduce the tumor effects seen in clinic and experimentally. Therefore, this paper has three parts: the exposition of the three models, the results and discussion (where we explore some aspects and examples of the solution of the models), and the conclusion.

Bone density of the axis (C2) measured using Hounsfield units of computed tomography

INTRODUCTION

The assessment of bone density is of great importance nowadays due to the increasing age of patients. Especially in regard to the surgical stabilization of the spine, the assessment of bone density is important for therapeutic decision making. The aim of this work was to record trabecular bone density values using Hounsfield units of the second cervical vertebra.

MATERIAL AND METHODS

The study is a monocentric retrospective data analysis of 198 patients who received contrast-enhanced polytrauma computed tomography in a period of two years at a maximum care hospital. Hounsfield units were measured in three different regions within the C2: dens, transition area between dens and vertebral body and vertebral body. The measured Hounsfield units were converted into bone density values using a validated formula.

RESULTS

A total of 198 patients were included. The median bone density varied in different regions of all measured C2 vertebrae: in the dens axis, C2 transition area between dens and vertebral body, and in the vertebral body bone densities were 302.79 mg/cm³, 160.08 mg/cm³, and 240.31 mg/cm³, respectively. The transition area from dens axis to corpus had statistically significant lower bone density values compared to the other regions (p < 0.001). There was a decrease in bone density values after age 50 years in both men and women (p < 0.001).

CONCLUSIONS

The transitional area from dens axis to corpus showed statistically significant lower bone density values compared to the adjacent regions (p < 0.001). This area seems to be a predilection site for fractures of the 2nd cervical vertebra, which is why special attention should be paid here in radiological diagnostics after a trauma.

Risk modeling of femoral neck fracture based on geometric parameters of the proximal epiphysis

BACKGROUND

Fractures of the proximal femur epiphysis are problematic for state health care because they are associated with severe medical and social problems and high morbidity and mortality rates.

AIM

To model the potential risk of hip fracture via femur geometric parameters.

METHODS

Seventy educational cadaveric femurs from people aged 14 to 80 years, 10 X-ray images from the records of the Human Anatomy Department and 10 X-ray images from the Department of Traumatology, Orthopedics and Disaster Surgery of Sechenov University, were evaluated. The parameters of the fractured bone were measured using images captured with a Canon d60 camera. The projection values of the proximal epiphysis of the cadaveric femurs and geometric parameters of the bones shown in the X-ray images were measured with Autodesk software (AutoCAD 2018). Analysis of the video frames showing bone rotation reveal that the greater trochanter can be inscribed in a parallelepiped, where one of the faces is parallel to the plane of view in the frontal standard projection and is rectangular. The angle of bone rotation obtained by turning the cube corresponded to the angle measured with the second technique. This reliable method of calculating the rotation of the bone relative to the anterior projection was employed in subsequent calculations. The geometric parameters of the femur were measured using X-ray images according to the proposed method.

RESULTS

The geometric parameters of 70 femurs were analyzed, and correlation coefficients were calculated. Our measurement results were compared with those reported by other authors. The potential influence of femur geometry on force distribution in the proximal epiphysis of the femur was described, and a 2-dimensional model of the femur epiphysis associated with minimal neck fracture risk was provided. The assessment of the geometric parameters of the femoral epiphysis indicated the greatest risk of a varus fracture of the neck if the angle of the minimal resistance zone (AMRZ) index > 24° and the neck-shaft angle (NSA) < 127.5°. In contrast, the minimum risk was observed at AMRZ < 14° and NSA > 128.87°.

CONCLUSION

The proposed method provides the potential femur neck fracture risk based on geometric parameters.

Pharmacological History of Missing Subjects: Perspective of a Correction Factor to Aid in the Study of Bone Remains

Simple Summary

The reconstruction of the biological profile of skeletal remains of missing subjects is also based on the analysis of the quality of bone tissue. The density of bone mass is a factor that allows us to inscribe the subject in a specific age group. Bone density varies not only according to age but also by the intake of certain drugs and certain abuse substances. The objective of our study is to propose the introduction of pharmacological history in the profile of missing persons. Information on drugs or abuse substances taken by the missing person is a useful corrective factor for tracing the chronological age of bone remains found, increasing the likelihood of identification. We emphasize the usefulness of this information also for the characterization of bone injuries and for the dating of antemortem fractures, useful elements to trace the cause and dynamics of death. The evaluation of these findings is also based on the characteristics of the bone tissue of skeletal remains which is also affected by any drugs and/or substances of abuse. Therefore, we believe that the pharmacological history of the missing subjects could be a new and interesting tool to help the activity of the forensic anthropologist.

Abstract

In forensic anthropology, bone mineral density and the estimation of the dating of fractures based on the degree of progress of healing processes are important parameters of study on bone remains. With our article we aim, on the one hand, to highlight the importance that these parameters have in the reconstruction of the biological profile of the subject, as well as the time and the cause of death; on the other hand, we aim to limit their variability according to the medical substances and/or abuse assumed during life by the subject. The aim of this article is to encourage the introduction of the pharmacological history of missing persons as a new correction factor for the study of bone remains, possibly based on new scientific studies that allow us to establish with greater specificity the effect that certain pharmacological therapies produce on bone mass and the speed of remodeling.

A Pilot Study to Assess Opportunistic Use of CT-Scan for Osteoporosis Screening in Chronic Pancreatitis

Objectives: CT scans are commonly performed in patients with chronic pancreatitis (CP). Osteopathy and fractures are recognized in CP but no osteoporosis screening guidelines are recommended. "Opportunistic" CT scan-derived bone density thresholds are assessed for identifying osteoporosis in CP. Methods: Retrospective pilot cohort study. CP subjects who had CT scans and dual-energy x-ray absorptiometry (DXA) within 1 year were included. CT-derived bone density was measured at the L1 level. Pearson's correlation was performed between age and CT-derived bone density in Hounsfield unit (HU). Univariate analysis using HU to identify osteoporosis was performed at various thresholds of bone density. The discriminatory ability of the model was evaluated with the area under the receiver operating characteristic (ROC) curve (AUC). Several HU thresholds were tested. Results: Twenty-seven CP subjects were included, of whom 11 had normal bone density, 12 osteopenia, and four osteoporosis on DXA. The mean age was 59.9 years (SD 13.0). There was a negative correlation of age with HU (r = -0.519, p = 0.006). CT-derived bone density predicted DXA-based osteoporosis in the univariable analysis (Odds Ratio (OR) = 0.97 95% Confidence Interval (CI) 0.94-1.00, p = 0.03). HU thresholds were tested. A threshold of 106 HU maximized the accuracy (AUC of 0.870). Conclusions: CT scan may be repurposed for "opportunistic" screening to rule out osteoporosis in CP. A larger study is warranted to confirm these results.

3-D cephalometry of the the orbit regarding endocrine orbitopathy, exophthalmos, and sex

Purpose

This study aimed at evaluating the orbital anatomy of patients concerning the relevance of orbital anatomy in the etiology of EO (endocrine orbitopathy) and exophthalmos utilizing a novel approach regarding three-dimensional measurements. Furthermore, sexual dimorphism in orbital anatomy was analyzed.

Methods

Orbital anatomy of 123 Caucasian patients (52 with EO, 71 without EO) was examined using computed tomographic data and FAT software for 3-D cephalometry. Using 56 anatomical landmarks, 20 angles and 155 distances were measured. MEDAS software was used for performing connected and unconnected t-tests and Spearman´s rank correlation test to evaluate interrelations and differences.

Results

Orbital anatomy was highly symmetrical with a mean side difference of 0.3 mm for distances and 0.6° for angles. There was a small albeit statistically significant difference in 13 out of 155 distances in women and 1 in men concerning patients with and without EO. Two out of 12 angles showed a statistically significant difference between female patients with and without EO. Regarding sex, statistically significant differences occurred in 39 distances, orbit volume, orbit surface, and 2 angles. On average, measurements were larger in men. Concerning globe position within the orbit, larger distances to the orbital apex correlated with larger orbital dimensions whereas the sagittal position of the orbital rim defined Hertel values.

Conclusion

In this study, little difference in orbital anatomy between patients with and without EO was found. Concerning sex, orbital anatomy differed significantly with men presenting larger orbital dimensions. Regarding clinically measured exophthalmos, orbital aperture anatomy is an important factor which has to be considered in distinguishing between true exophthalmos with a larger distance between globe and orbital apex and pseudoexophthalmos were only the orbital rim is retruded. Thus, orbital anatomy may influence therapy regarding timing and surgical procedures as it affects exophthalmos.

Longitudinal Tibia Stress Fracture Risk During High-Volume Training: a Multi-Scale Modeling Pipeline Incorporating Bone Remodeling

Tibia stress fractures are prevalent during high-intensity training, yet a mechanistic model linking longitudinal training intensity, bone health, and long-term injury risk has yet to be demonstrated. The objective of this study was to develop and validate a multi-scale model of gross and tissue level loading on the tibia including bone remodeling on a timescale of week. Peak tensile tibial strain (3517 µstrain) during 4 m/s running was below injury thresholds, and the peak anteromedial tibial strain (1248 µstrain) was 0.17 standard deviations away from the mean of reported literature values. An initial study isolated the effects of cortical density and stiffness on tibial strain during a simulated eight week training period. Tibial strains and cortical microcracking correlated to initial cortical modulus, with all simulations presenting peak anteromedial tensile strains (1047-1600 µstrain) near day 11. Average cortical densities decreased by 7-8 percent of their nominal value by day 11, but the overall density change was &lt;2% by the end of the simulated training period, in line with reported results. This study demonstrates the benefits of multi-scale models for investigating stress fracture risk and indicates that peak tibial strain, and thus injury risk, may increase early in a high intensity training program. Future studies could optimize training volume and recovery time to reduce injury risk during the most vulnerable training periods.

A General Mechano-Pharmaco-Biological Model for Bone Remodeling Including Cortisol Variation

The process of bone remodeling requires a strict coordination of bone resorption and formation in time and space in order to maintain consistent bone quality and quantity. Bone-resorbing osteoclasts and bone-forming osteoblasts are the two major players in the remodeling process. Their coordination is achieved by generating the appropriate number of osteoblasts since osteoblastic-lineage cells govern the bone mass variation and regulate a corresponding number of osteoclasts. Furthermore, diverse hormones, cytokines and growth factors that strongly link osteoblasts to osteoclasts coordinated these two cell populations. The understanding of this complex remodeling process and predicting its evolution is crucial to manage bone strength under physiologic and pathologic conditions. Several mathematical models have been suggested to clarify this remodeling process, from the earliest purely phenomenological to the latest biomechanical and mechanobiological models. In this current article, a general mathematical model is proposed to fill the gaps identified in former bone remodeling models. The proposed model is the result of combining existing bone remodeling models to present an updated model, which also incorporates several important parameters affecting bone remodeling under various physiologic and pathologic conditions. Furthermore, the proposed model can be extended to include additional parameters in the future. These parameters are divided into four groups according to their origin, whether endogenous or exogenous, and the cell population they affect, whether osteoclasts or osteoblasts. The model also enables easy coupling of biological models to pharmacological and/or mechanical models in the future.

Influence of cancellous bone material and dead zone on stress-strain, bone stimulus and bone remodelling around the tibia for total ankle replacement

Extreme bone resorption due to bone remodelling is one of the reasons for ankle component loosening. Finite element (FE) analysis has been effectively used nowadays for pre-clinical analysis of orthopaedic implants. For FE modelling, the selection of bone material and dead zone play a vital role to understand the bone remodelling. This study deals with the effects of different cancellous elastic modulus-density relationships and dead zone on bone remodelling around the tibia owing to total ankle replacement (TAR), using finite element analysis with physiological loading conditions. This study also investigated the bone stimulus distribution in the tibia to identify the initial indication of bone density changes due to bone remodelling. Additionally, the Hoffman failure criterion was used to investigate the chances of implant-bone interface failure due to different cancellous bone material modelling and bone remodelling. The present bone remodelling study consists of three different dead or lazy zones (±0.75, ±0.60 and ±0.35) to examine the influence of the dead zone on bone remodelling. Differences in stress/strain distribution were observed in the tibia bone due to different cancellous bone material modelling. Despite little variations, bone density changes due to bone remodelling were found to be almost similar for two FE models having different cancellous bone material. Similar to these results, the effect of different dead zone on bone density changes due to bone remodelling was found to be minimal. Bone stimulus distribution in the cancellous bone was found to be almost similar for FE models having different cancellous bone material modelling and different dead zones. To understand the stress/strain and interface related failure of the tibial component, cancellous bone material modelling plays a crucial role. However, cancellous bone material modelling and dead zone have minimal influence on bone remodelling around the tibia cancellous bone due to TAR.

NUMERICAL SIMULATION OF APPARENT DENSITY EVOLUTION OF TRABECULAR BONE UNDER FATIGUE LOADING: EFFECT OF BONE INITIAL PROPERTIES

Bone remodeling is a physiological phenomenon coupling resorption and formation processes that are mainly mediated by osteoclasts and osteoblasts, in response to mechanical stimuli transduced by osteocytes to biochemical signals activating the bone multicellular unit. Under normal loading conditions, bone resorption and formation are balanced by a homeostasis process. When bone is subjected to overstress, microdamaging occurs, which induces a modification of the structural integrity and microarchitecture. This has drawn significant attention to the mechanical properties of bone. In this context, the current study has been carried out with the aim of numerically investigating the impact of the mechanical properties on the remodeling process of the trabecular bone under cyclic loading, highlighting the effects of different values of the mineral density and the Young’s modulus. This was performed using a mechanobiological model, coupling mechanical and biological approaches, allowing to numerically simulate the effect of the selected parameters for a 20-year-period of cyclic loading for 2D and 3D models of a human femur head. The current work is an explorative numerical study, and the obtained results revealed the changes in the overall stiffness of the bone according to the mechanical properties.