Characterising variability and regional correlations of microstructure and mechanical competence of human tibial trabecular bone: An in-vivo HR-pQCT study

Mechanical impact of regional structural deterioration and tissue-level compensation on proximal femur trabecular bone

Introduction

Osteoporosis-induced changes in bone structure and composition significantly reduce bone strength, particularly in the human proximal femur. This study examines how these changes affect the mechanical performance of trabecular bone to enhance diagnosis, prevention, and treatment strategies.

Methods

A proximal femur sample was scanned using micro-CT at 40 μm resolution. Five regions of interest were selected within the femoral head, femoral neck, and greater trochanter. Structural models simulating various stages of osteoporosis were created using image processing software. Micro-finite element analysis evaluated the mechanical properties of trabecular bone under different conditions of structural deterioration and tissue-level elastic modulus variations. The combined effects of structural deterioration and tissue-level mechanical properties on trabecular bone mechanical performance were further analyzed.

Results

The mechanical performance of trabecular bone generally follows a power-law relationship with its microstructural characteristics. However, in any specific region, the apparent mechanical properties linearly decrease with structural deterioration. The femoral neck and greater trochanter are more sensitive to structural deterioration than the femoral head. A 5% bone mass loss in the femoral head led to a 7% reduction in mechanical performance, while the femoral neck experienced a 12% loss. Increasing tissue-level elastic modulus improved mechanical performance, partially offsetting bone mass reduction effects.

Conclusion

Trabecular bone in low bone mass regions is more affected by bone mass loss. Structural deterioration primarily reduces bone strength, but improvements in tissue-level properties can mitigate this effect, especially in early osteoporosis. Targeted assessments and interventions are crucial for effective management. Future research should explore heterogeneous deterioration models to better understand osteoporosis progression.

Comprehensively characterizing heterogeneous and transversely isotropic properties of femur cortical bones

Comprehensive characterization of the transversely isotropic mechanical properties of long bones along both the longitudinal and circumferential gradients is crucial for developing accurate mathematical models and studying bone biomechanics. In addition, mechanical testing to derive elastic, plastic, and failure properties of bones is essential for modeling plastic deformation and failure of bones. To achieve these, we machined a total of 336 cortical specimens, including 168 transverse and 168 longitudinal specimens, from four different quadrants of seven different sections of 3 bovine femurs. We conducted three-point bending tests of these specimens at a loading rate of 0.02 mm/s. Young's modulus, yield stress, tangential modulus, and effective plastic strain for each specimen were derived from correction equations based on classical beam theory. Our statistical analysis reveals that the longitudinal gradient has a significant effect on the Young's modulus, yield stress, and tangential modulus of both longitudinal and transverse specimens, whereas the circumferential gradient significantly influences the Young's modulus, yield stress, and tangential modulus of transverse specimens only. The differences in Young's modulus and yield stress between longitudinal specimens from different sections are greater than 40%, whereas those between transverse specimens are approximately 30%. The Young's modulus and yield stress of transverse specimens in the anterior quadrant were 18.81%/15.46% and 18.34%/14.88% higher than those in the posterior and lateral quadrants, respectively. There is no significant interaction between the longitudinal gradient and the circumferential gradient. Considering the transverse isotropy, it is crucial to consider loading direction when investigating the impact of circumferential gradients in the anterior, lateral, medial, and posterior directions. Our findings indicate that the conventional assumption of homogeneity in simulating the cortical bone of long bones may have limitations, and researchers should consider the anatomical position and loading direction of femur specimens for precise prediction of mechanical responses.

Fracture in porous bone analysed with a numerical phase-field dynamical model

A dynamic phase-field fracture finite element model is applied to discretized high-resolution three-dimensional computed tomography images of human trabecular bone to analyse rapid bone fracture. The model is contrasted to quasi-static experimental results and a quasi-static phase-field finite element model. The experiment revealed complex stepwise crack evolution with multiple crack fronts, and crack arrests, as the global tensile displacement load was incrementally increased. The quasi-static phase-field fracture model captures the fractures in the experiment reasonably well, and the dynamic model converges towards the quasi-static model when mechanically loaded at low rates. At higher load rates, i.e., at larger impulses, inertia effects significantly contribute to an increased initial global stiffness, higher peak forces and a larger number of cracks spread over a larger volume. Since the fracture process clearly is different at large impulses compared to small impulses, it is concluded that dynamic fracture models are necessary when simulating rapid bone fracture.

Multiscale Femoral Neck Imaging and Multimodal Trabeculae Quality Characterization in an Osteoporotic Bone Sample

Although multiple structural, mechanical, and molecular factors are definitely involved in osteoporosis, the assessment of subregional bone mineral density remains the most commonly used diagnostic index. In this study, we characterized bone quality in the femoral neck of one osteoporotic patients as compared to an age-matched control subject, and so used a multiscale and multimodal approach including X-ray computed microtomography at different spatial resolutions (pixel size: 51.0, 4.95 and 0.9 µm), microindentation and Fourier transform infrared spectroscopy. Our results showed abnormalities in the osteocytes lacunae volume (358.08 ± 165.00 for the osteoporotic sample vs. 287.10 ± 160.00 for the control), whereas a statistical difference was found neither for shape nor for density. The osteoporotic femoral head and great trochanter reported reduced elastic modulus (Es) and hardness (H) compared to the control reference (−48% (p < 0.0001) and −34% (p < 0.0001), respectively for Es and H in the femoral head and −29% (p < 0.01) and −22% (p < 0.05), respectively for Es and H in the great trochanter), whereas the corresponding values in the femoral neck were in the same range. The spectral analysis could distinguish neither subregional differences in the osteoporotic sample nor between the osteoporotic and healthy samples. Although, infrared spectroscopic measurements were comparable among subregions, and so regardless of the bone osteoporotic status, the trabecular mechanical properties were comparable only in the femoral neck. These results illustrate that bone remodeling in osteoporosis is a non-uniform process with different rates in different bone anatomical regions, hence showing the interest of a clear analysis of the bone microarchitecture in the case of patients’ osteoporotic evaluation.

The developing juvenile distal tibia: Radiographic identification of distinct ontogenetic phases and structural trajectories

A novel combination of radiographic colour gradient mapping and radiographic absorptiometry was utilised to examine 96 human distal tibiae from 54 individuals ranging in age-at-death from the foetal to 23 years. The purpose of this was to identify previously undocumented changes in the internal organisation during the development of the distal tibia and determine whether these changes could be described as distinct phases. Previous studies have demonstrated a rudimentary structural organisation in other skeletal elements that mirror more mature patterns of bone organisation. Results showed that the perinatal tibia did not exhibit a rudimentary structural pattern similar to the architecture observed within the late adolescent tibia. This lack of early internal organisation is hypothesised to be related to the rudimentary ossification process that is being laid down around a pre-existing vascular template which will be subsequently modified by locomotive forces. Between birth and 2 years of age, the tibia exhibited a period of regression where radiodensity decreased in comparison to the perinatal tibia. This period of regression was postulated to be due to a combination of factors including changing locomotive forces, weaning and growth resulting in a stage of development which is extremely demanding on calcium liberation from the skeleton. After 2 years of age, the distal tibia demonstrated refinement where radiographic trajectories progressively developed into patterns consistent with adult trabecular organisation. These trajectories are linked to the forces associated with the bipedal gait, suggesting a strong influence of biomechanical forces on the development of the distal tibia.

Application of phase-field fracture theories and digital volume correlation to synchrotron X-ray monitored fractures in human trabecular bone: A case study

Fracture processes of trabecular bone have been studied using various approaches over the years. However, reliable methods to analyse fracture at the single trabecula level are limited. In this study, a digital volume correlation (DVC) and a phase-field fracture model are applied and contrasted for human trabecular bone to analyse its failure under global compression at high resolution. A human trabecular bone sample was fractured in situ under synchrotron-based X-ray micro computed tomography (CT). Reconstructed CT data was then used in DVC algorithms to obtain high-resolution displacement fields in the bone at different load steps. A high-resolution specimen-specific structural mesh was discretized from the CT data and used for the phase-field simulation of the fracturing bone. The DVC analysis showed opening mode cracks as well as shear mode cracks. Strains in cracked regions were analysed. The load distribution in the trabecular structure resulted in two completely separated fracture regions in the sample body. A phenomenon that was also captured in the phase-field model. The results encourage us to believe improvements in boundary conditions and material models are worthwhile pursuing. Findings in this study support further development of a phase-field method to analyse fracture in samples with complex morphology, such as trabecular bone, and the capacity of DVC to quantify strains and slowly growing stable fractures during step-wise loading of trabecular bone.

Three-dimensional radiological anatomy of condyle trabecular bone based on a Volume-of-Interest analysis

OBJECTIVES

Three-dimensional radiological anatomic characteristics of condyle trabeculae was obtained quantitatively based on a volume-of-interest (VOI) analysis.

METHODS

Nine human mandibular condyle specimens were scanned by micro-computed tomography (micro-CT). A total of 34 VOIs were selected from each condyle specimen, which were divided into six layers and four parts to analyze the morphological characteristics of trabeculae based on cylindrical VOIs with a diameter and height of 2 mm. One-way analysis of variance was used to compare the regional differences of morphological parameters among each layer and part.

RESULTS

Values for bone mineral density, bone volume/total volume, trabecular thickness, and trabecular bone number were greater in the anterior part compared with the posterior part; and the lateral part was larger than the medial part in the first, second, and third layers, while the medial part was larger in the fourth and fifth layers; these values in the first and sixth layers were much larger, while those in the third and fourth layers were smaller. Bone surface area/bone volume, trabecular spacing, and trabecular bone pattern factor were larger in the posterior part than in the anterior part; and the lateral part was larger than the medial part in the fourth and fifth layers, while the medial part was larger in the first and second layers.

CONCLUSIONS

The morphological distribution of VOIs was anisotropic within trabecular bone of human mandibular condyles. The upper and lower ends of trabecular bone were much more compact, with higher bone density, trabecular thickness, and trabecular number than in the middle layers.

Characterizing the Mechanical Behavior of Bone and Bone Surrogates in Compression Using pQCT

Many axial and appendicular skeleton bones are subjected to repetitive loading during daily activities. Until recently, the structural analysis of fractures has been limited to 2D sections, and the dynamic assessment of fracture progression has not been possible. The structural failure was analyzed using step-wise micro-compression combined with time-lapsed micro-computed tomographic imaging. The structural failure was investigated in four different sample materials (two different bone surrogates, lumbar vertebral bodies from bovine and red deer). The samples were loaded in different force steps based on uniaxial compression tests. The micro-tomography images were used to create three-dimensional models from which various parameters were calculated that provide information about the structure and density of the samples. By superimposing two 3D images and calculating the different surfaces, it was possible to precisely analyze which trabeculae failed in which area and under which load. According to the current state of the art, bone mineral density is usually used as a value for bone quality, but the question can be raised as to whether other values such as trabecular structure, damage accumulation, and bone mineralization can predict structural competence better than bone mineral density alone.

A computational framework for canonical holistic morphometric analysis of trabecular bone

Bone is a remarkable, living tissue that functionally adapts to external loading. Therefore, bone shape and internal structure carry information relevant to many disciplines, including medicine, forensic science, and anthropology. However, morphometric comparisons of homologous regions across different individuals or groups are still challenging. In this study, two methods were combined to quantify such differences: (1) Holistic morphometric analysis (HMA) was used to quantify morphometric values in each bone, (2) which could then be mapped to a volumetric mesh of a canonical bone created by a statistical free-form deformation model (SDM). Required parameters for this canonical holistic morphometric analysis (cHMA) method were identified and the robustness of the method was evaluated. The robustness studies showed that the SDM converged after one to two iterations, had only a marginal bias towards the chosen starting image, and could handle large shape differences seen in bones of different species. Case studies were performed on metacarpal bones and proximal femora of different primate species to confirm prior study results. The differences between species could be visualised and statistically analysed in both case studies. cHMA provides a framework for performing quantitative comparisons of different morphometric quantities across individuals or groups. These comparisons facilitate investigation of the relationship between spatial morphometric variations and function or pathology, or both.

Trabecular bone porosity and pore size distribution in osteoporotic patients - A low field nuclear magnetic resonance and microcomputed tomography investigation

The study of bone morphology is of great importance as bone morphology is influenced by factors such as age and underlying comorbidities and is associated with bone mechanical properties and fracture risk. Standard diagnostic techniques used in bone disease, such as Dual-Energy X-ray absorptiometry and ultrasonography do not provide qualitative and quantitative morphological information. In recent years, techniques such as High Resolution Computed Tomography (HR-CT), micro- CT, Magnetic Resonance Imaging (MRI), and Low Field Nuclear Magnetic Resonance (LF-NMR) have been developed for the study of bone structure and porosity. Data obtained from these techniques have been used to construct models to predict bone mechanical properties thanks to finite element analysis. Cortical porosity has been extensively studied and successfully correlated with disease progression and mechanical properties. Trabecular porosity and pore size distribution, however, have increasingly been taken into consideration to obtain a comprehensive analysis of bone pathology and mechanic. Therefore, we have decided to evaluate the ability of micro- CT (chosen for its high spatial resolving power) and LF-NMR (chosen to analyze the behavior of water molecules within trabecular bone pores) to characterize the morphology of trabecular bone in osteoporosis. Trabecular bone samples from human femoral heads collected during hip replacement surgery were from osteoporosis (test group) and osteoarthritis (control group) patients. Our data show that both micro- CT and LF-NMR can detect qualitative changes in trabecular bone (i.e., transition from plate-like to rod-like morphology). Micro- CT failed to detect significant differences in trabecular bone morphology parameters between osteoporotic and osteoarthritic specimens, with the exception of Trabecular Number and Connectivity Density, which are markers of osteoporosis progression. In contrast, LF-NMR was able to detect significant differences in porosity and pore size of trabecular bone from osteoporotic versus osteoarthritic (control) samples. However, only the combination of these two techniques allowed the detection of structural morphometric changes (increase in the larger pore fraction and enlargement of the larger pores) in the trabecular bone of osteoporotic specimens compared to osteoarthritic ones. In conclusion, the combined use of LF-NMR and micro- CT provides a valuable tool for characterizing the morphology of trabecular bone and may offer the possibility for a new approach to the study and modeling of bone mechanics in the context of aging and disease.

The plate-to-rod transition in trabecular bone loss is elusive

Changes in trabecular micro-architecture are key to our understanding of osteoporosis. Previous work focusing on structure model index (SMI) measurements have concluded that disease progression entails a shift from plates to rods in trabecular bone, but SMI is heavily biased by bone volume fraction. As an alternative to SMI, we proposed the ellipsoid factor (EF) as a continuous measure of local trabecular shape between plate-like and rod-like extremes. We investigated the relationship between EF distributions, SMI and bone volume fraction of the trabecular geometry in a murine model of disuse osteoporosis as well as from human vertebrae of differing bone volume fraction. We observed a moderate shift in EF median (at later disease stages in mouse tibia) and EF mode (in the vertebral samples with low bone volume fraction) towards a more rod-like geometry, but not in EF maximum and minimum. These results support the notion that the plate to rod transition does not coincide with the onset of bone loss and is considerably more moderate, when it does occur, than SMI suggests. A variety of local shapes not straightforward to categorize as rod or plate exist in all our trabecular bone samples.

High-impact exercise stimulated localised adaptation of microarchitecture across distal tibia in postmenopausal women

We provided evidence that a 6-month regular hopping exercise intervention can increase trabecular number and possibly trabecular volume fraction of the distal tibia. Our novel localised analysis demonstrated region-specific changes, predominantly in the anterior region, in postmenopausal women.

INTRODUCTION

The localisation of bone remodelling and microarchitectural adaptation to exercise loading has not been demonstrated previously in vivo in humans. The aim of this study is to assess the feasibility of using 3D image registration and high-resolution peripheral quantitative computed tomography (HR-pQCT) to investigate the effect of high-impact exercise on human trabecular bone variables and remodelling rate across the distal tibia.

METHODS

Ten postmenopausal women were recruited for 6-month unilateral hopping exercises, with HR-pQCT scans taken of both exercise leg (EL) and control leg (CL) for each participant before and after the intervention. A 3D image registration was used to ensure measurements were taken at the same region. Short-term reproducibility tests were conducted prior to the assessment using identical setup. The results were assessed comparing CL and EL, and interaction (time × leg) using a two-way repeated measures analysis of variance (RM-ANOVA).

RESULTS

Across the whole tibia, we observed significant increases in trabecular number (Tb.N) (+ 4.4%) and trabecular bone formation rate (tBFR) (3.3%), and a non-significant increase in trabecular bone volume fraction (BV/TV) (+ 1%) in the EL. Regional resorption was higher in the CL than the EL, with this difference being statistically significant at the lateral tibia. In the EL, tBFR was significantly higher in the anterior region than the medial but a trabecular bone resorption rate (tBRR) showed no significant regional variation. Conversely in the CL, both tBFR and tBRR were significantly higher in the anterior and lateral than the medial region.

CONCLUSION

We demonstrated that it was possible to detect exercise-related bone adaptation with 3D registration of HR-pQCT scan data. Regular hopping exercise increased Tb.N and possibly BV/TV across the whole distal tibia. A novel finding of the study was that tBFR and tBRR responses to loading were localised: changes were achieved by formation rate exceeding resorption rate in the exercise leg, both globally and at the anterior region where turnover was greatest.

TRIAL REGISTRATION

clinicaltrials.gov : NCT03225703.

Remodelling of trabecular bone in human distal tibia: A model based on an in-vivo HR-pQCT study

An abnormal remodelling process of bones can lead to various bone disorders, such as osteoporosis, making them prone to fracture. Simulations of load-induced remodelling of trabecular bone were used to investigate its response to mechanical signal. However, the role of mechanostat in trabecular-bone remodelling has not yet been investigated in simulations underpinned by a longitudinal in-vivo study in humans. In this work, a finite-element model based on a 6-month longitudinal in-vivo HR-pQCT study was developed and validated to investigate the effect of mechanical stimuli on bone remodelling. The simulated changes in microstructural parameters and density of trabecular bone were compared with respective experimental results. A maximum principal strain (MPS) and a maximum principal strain gradient (∇MPS) were used as mechanical signals to drive a five-stage mechanostat remodelling model, including additional over-strain and damage stages. It was found that the density distribution varied with the studied mechanical signals, along with decreasing with time levels of bone volume fraction BV/TV, trabecular thickness Tb.Th and bone surface area Tb.BS as well as increased trabecular separation Tb.Sp. Among these parameters, BV/TV and Tb.Th together with the bone-remodelling parameters from the MPS model demonstrated a significant correlation with the experimental data. The developed model provides a good foundation for further development and investigation of the relationships between mechanical loading and human-bone microarchitecture.

High resolution imaging in bone tissue research-review

This review article focuses on imaging of bone tissue to understand skeletal health with regards to bone quality. Skeletal fragility fractures are due to bone diseases such as osteoporosis which result in low bone mass and bone mineral density (BMD) leading to high risk of fragility fractures. Recent advances in imaging and analysis technologies have highly benefitted the field of biological sciences. In particular, their application in skeletal health has been of significant importance in understanding bone mechanical behavior (structure and properties) at the tissue level. While synchrotron based microCT technique has remained the gold standard for non-destructive evaluation of structure in material and biological sciences, several lab based microCT systems have been developed to provide high resolution imaging of specimens with greater access, and ease of use in laboratory settings. Lab based microCT scanners are widely used in the bone field as a standard tool to evaluate three-dimensional (3D) morphologies of bone structure at image resolutions appropriate for bone samples from small animals to bone biopsy specimens from humans. Both synchrotron and standard lab based microCT systems provide high resolution imaging ex vivo for a small sized specimen. A few X-ray based systems are also commercially available for in vivo scanning at relatively low image resolutions. Synchrotron-based CT microscopy is being used for various ultra-high-resolution image analyses using complex 3D software. However, the synchrotron-based CT technology is in high demand, allows only limited numbers of specimens, expensive, requires complex additional instrumentation, and is not easily available to researchers as it requires access to a synchrotron source which is always limited. Therefore, desktop laboratory scanners (microXCT, Zeiss/Xradia, Scanco, SkyScan. etc.), mimicking the synchrotron based CT technology or image resolution, have been developed to solve the accessibility issues. These lab based scanners have helped both material science, and the bone field to investigate bone tissue morphologies at submicron mage resolutions. Considerable progress has been made in both in vivo and ex vivo imaging towards providing high resolution images of bone tissue. Both clinical and research imaging technologies will continue to improve and help understand osteoporosis and other related skeletal issues in order to develop targeted treatments for bone fragility. This review summarizes the high resolution imaging work in bone research.

Global and Spatial Compartmental Interrelationships of Bone Density, Microstructure, Geometry and Biomechanics in the Distal Radius in a Colles' Fracture Study Using HR-pQCT

BACKGROUND

Bone parameters derived from HR-pQCT have been investigated on a parameter-by-parameter basis for different clinical conditions. However, little is known regarding the interrelationships of bone parameters and the spatial distribution of these interrelationships. In this work: 1) we investigate compartmental interrelationships of bone parameters; 2) assess the spatial distribution of interrelationships of bone parameters; and 3) compare interrelationships of bone parameters between postmenopausal women with and without a recent Colles' fracture.

METHODS

Images from the unaffected radius in fracture cases (n=84), and from the non-dominant radius of controls (n=98) were obtained using HR-pQCT. Trabecular voxel-based maps of local bone volume fraction (L.Tb.BV/TV), homogenized volumetric bone mineral density (H.Tb.BMD), homogenized μFEA-derived strain energy density (H.Tb.SED), and homogenized inter-trabecular distances (H.Tb.1/N) were generated; as well as surface-based maps of apparent cortical bone thickness (Surf.app.Ct.Th), porosity-weighted cortical bone thickness (Surf.Ct.SIT), mean cortical BMD (Surf.Ct.BMD), and mean cortical SED (Surf.Ct.SED). Anatomical correspondences across the parametric maps in the study were established via spatial normalization to a common template. Mean values of the parametric maps before spatial normalization were used to assess compartmental Spearman's rank partial correlations of bone parameters (e.g., between H.Tb.BMD and L.Tb.BV/TV or between Surf.Ct.BMD and Surf.app.Ct.Th). Spearman's rank partial correlations were also assessed for each voxel and vertex of the spatially normalized parametric maps, thus generating maps of Spearman's rank partial correlation coefficients. Correlations were performed independently within each group, and compared between groups using the Fisher's Z transformation.

RESULTS

All within-group global trabecular and cortical Spearman's rank partial correlations were significant; and the correlations of H.Tb.BMD-L.Tb.BV/TV, H.Tb.BMD-H.Tb.1/N, L.Tb.BV/TV-H.Tb.1/N, Surf.Ct.BMD-Surf.Ct.SED and Surf.Ct.SIT-Surf.Ct.SED were significantly different between controls and fracture cases. The spatial analyses revealed significant heterogeneous voxel- and surface-based correlation coefficient maps across the distal radius for both groups; and the correlation maps of H.Tb.BMD-L.Tb.BV/TV, H.Tb.BMD-H.Tb.1/N, L.Tb.BV/TV-H.Tb.1/N, H.Tb.1/N-H.Tb.SED and Surf.app.Ct.Th - Surf.Ct.SIT yielded small clusters of significant correlation differences between groups.

DISCUSSION

The heterogeneous spatial distribution of interrelationships of bone parameters assessing density, microstructure, geometry and biomechanics, along with their global and local differences between controls and fracture cases, may help us further understand different bone mechanisms of bone fracture.

Simulation analysis of impact damage to the bone tissue surrounding a dental implant

Dental implant may suffer transient external impacts. To simulate the effect of impact forces on bone damage is very important for evaluation of damage and guiding treatment in clinics. In this study, an animal model was established by inserting an implant into the femoral condyle of New Zealand rabbit. Implant with good osseointegration was loaded with impact force. A three-dimensional finite element model was established based on the data of the animal model. Damage process to bone tissue was simulated with Abaqus 6.13 software combining dynamic mechanical properties of the femur. The characteristics of bone damage were analyzed by comparing the results of animal testing with numerical simulation data. After impact, cortical bone around the implant and trabecular at the bottom of the implant were prone to damage. The degree of damage correlated with the direction of loading and the magnitude of the impact. Lateral loading was most likely performed to damage cancellous bone. The stress wave formed by the impact force can damage the implant–bone interface and peri-implant trabeculae. The data from numerical simulations were consistent with data from animal experiments, highlighting the importance of a thorough examination and evaluation based on the patient’s medical history.