Three-dimensional methods for quantification of cancellous bone architecture

Effect of trabecular architectures on the mechanical response in osteoporotic and healthy human bone

Research at the mesoscale bone trabeculae arrangement yields intriguing results that, due to their clinical resolution, can be applied in clinical field, contributing significantly to the diagnosis of bone-related diseases. While the literature offers quantitative morphometric parameters for a thorough characterization of the mesoscale bone network, there is a gap in understanding relationships among them, particularly in the context of various bone pathologies. This research aims to bridge these gaps by offering a quantitative evaluation of the interplay among morphometric parameters and mechanical response at mesoscale in osteoporotic and non-osteoporotic bones. Bone mechanical response, dependent on trabecular arrangement, is defined by apparent stiffness, computationally calculated using the Gibson-Ashby model. Key findings indicate that: (i) in addition to bone density, measured using X-ray absorptiometry, trabecular connectivity density, trabecular spacing and degree of anisotropy are crucial parameters for characterize osteoporosis state; (ii) apparent stiffness values exhibit strong correlations with bone density and connectivity density; (iii) connectivity density and degree of anisotropy result the best predictors of mechanical response. Despite the inherent heterogeneity in bone structure, suggesting the potential benefit of a larger sample size in the future, this approach presents a valuable method to enhance discrimination between osteoporotic and non-osteoporotic samples.

Age- and sex-related changes in vertebral trabecular bone architecture in Neolithic and Mediaeval populations from Poland

This paper investigates trabecular bone ontogenetic changes in two different Polish populations, one prehistoric and the other historical. The studied populations are from the Brześć Kujawski region in Kujawy (north-central Poland), one from the Neolithic Period (4500-4000 BC) and one from the Middle Ages (twelfth-sixteenth centuries AD), in total 62 vertebral specimens (32 males, 30 females). Eight morphometric parameters acquired from microCT scan images were analysed. Two-way ANOVA after Box-Cox transformation and multifactorial regression model were calculated. A significant decrease in percentage bone volume fraction (BV/TV; [%]) with age at death was observed in the studied sample; Tb.N (trabecular number) was also significantly decreased with age; trabecular separation (Tb.Sp) increased with advancing age; connectivity density (Conn.D) was negatively correlated with biological age and higher in the Neolithic population. These data are found to be compatible with data from the current biomedical literature, while no loss of horizontal trabeculae was recorded as would be expected based on modern osteoporosis.

Method for accurate removal of trabecular bone samples from a curved articulating surface of the distal femur

BACKGROUND

Knowing the mechanical properties of trabecular bone is critical for many branches of orthopaedic research. Trabecular bone is anisotropic and the principal trabecular direction is usually aligned with the load it transmits. It is therefore critical that the mechanical properties are measured as close as possible to this direction, which is often perpendicular to a curved articulating surface.

METHODS

This study presents a method to extract trabecular bone cores perpendicular to a curved articulating surface of the distal femur. Cutting guides were generated from computed tomography scans of 12 human distal femora and a series of cutting tools were used to release cylindrical bone cores from the femora. The bone cores were then measured to identify the angle between the bone core axis and the principal trabecular axis.

FINDINGS

The method yielded an 83% success rate in core extraction over 10 core locations per distal femur specimen. In the condyles, 97% of extracted cores were aligned with the principal trabecular direction.

INTERPRETATION

This method is a reliable way of extracting trabecular bone specimens perpendicular to a curved articular surface and could be useful across the field of orthopaedic research.

Microcomputed tomographic analysis of the ulnar trochlear notch in medium- and large-breed canine cadavers with and without medial coronoid disease

Medial coronoid disease (MCD) is a common disease often associated with thoracic limb lameness in medium- and large-breed dogs. The term MCD includes subchondral changes of the medial coronoid process (MCP), as well as pathologies of both cartilaginous surfaces. As there are only a few comprehensive and detailed studies on the trabecular structure of the medial coronoid, the goal of this study was to compare the trabecular structure of the ulnar trochlear notch of canine cadavers with and without MCD using different micro-computed tomographic (micro-CT) parameters. Fifty-eight elbow joints from 29 canine cadavers of MCD-predisposed and non-predisposed breeds (control group) were examined radiographically, macroscopically and by microcomputed tomography. The study included elbow joints of eight Labrador Retrievers (21.6-37 kg), seven Golden Retrievers (26.3-42 kg), seven Bernese Mountain dogs (31-47 kg) and seven dogs of non-predisposed breeds (19.7-52 kg) (control group). The final diagnosis of MCD was based on necropsy and micro-computed tomographic examinations. Micro-CT examinations were performed using XtremeCT II (Scanco Medical, Zurich, Switzerland) and the following parameters were examined: bone volume fraction (BV/TV), trabecular thickness (Tb.Th), trabecular separation (Tb.Sp), trabecular number (Tb.N), connectivity density (Conn.D) and degree of anisotropy (DA). Twenty-four elbows of 44 elbows of the predisposed breeds (Labrador Retrievers, Golden Retrievers Bernese Mountain dogs) showed subchondral changes and lesions of the cartilage surfaces. The result of this study is a higher density (BV/TV) of the trabecular bone of the ulnar trochlear notch in elbows affected by MCD compared to the control group. The increased density due to trabecular reconstruction in the ulnar trochlear notch is likely the result of selectively increased loading during life.

The trabecular architecture of the popliteal sesamoid bone (cyamella) from a New Zealand white rabbit (Oryctolagus cuniculus)

Sesamoid bones are ossified structures that are embedded in tendons near articulation. They consist of an inner trabecular bone architecture surrounded by a thin cortical shell. While the formation of sesamoid bones is probably mainly controlled by genetic factors, the proper development and mineralization of a sesamoid bone depends also on mechanical stimulation. While most sesamoid bones are not loaded directly by other bones during locomotion, they still experience forces directed from the tendon in which they are embedded. In cases when the sesamoid bone is experiencing forces only from a single tendon, such as the cyamella in the rabbit, this may give us a tool to study bone functional adaptation in a relatively simple loading setting. This study investigates the internal trabecular architecture of the popliteal sesamoid bone (cyamellae) in New Zealand white (NZW) rabbits (Oryctolagus cuniculus). Five hind limbs of NZW rabbits were micro-computed tomography scanned and the cortical and trabecular architectures of the cyamellae were evaluated. The results revealed that similar to the patella, the cyamella has a thin cortex and a high trabecular bone volume fraction (BV/TV), which is derived mostly from the high trabecular thickness (Tb.Th). Trabecular BV/TV and Tb.Th were not distributed homogeneously, but they were lower at the periphery and higher closer to the proximal and middle of the cyamella, near the musculotendinous junction. The results also demonstrated that trabeculae tend to align along two recognizable orientations, one with the direction of tensile stresses, in line with the popliteal tendon, and the second bridging the narrow space between the cranial and caudal cortical faces of the bone.

An improved trabecular bone model based on Voronoi tessellation

BACKGROUND AND OBJECTIVE

Accurate numerical and physical models of trabecular bone, correctly representing its complexity and variability, could be highly advantageous in the development of e.g. new bone-anchored implants due to the limited availability of real bone. Several Voronoi tessellation-based porous models have been reported in the literature, attempting to mimic the trabecular bone. However, these models have been limited to lattice rod-like structures, which are only structurally representative of very high-porosity trabecular bone. The objective of this study was to provide an improved model, more representative of trabecular bone of different porosity.

METHODS

Boolean operations were utilized to merge scaled Voronoi cells, thereby introducing different structural patterns, controlling porosity and to some extent anisotropy. The mechanical properties of the structures were evaluated using analytical estimations, numerical simulations, and experimental compression tests of 3D-printed versions of the structures. The capacity of the developed models to represent trabecular bone was assessed by comparing some key geometric features with trabecular bone characterized in previous studies.

RESULTS

The models gave the possibility to provide pore interconnectivity at relatively low porosities as well as both plate- and rod-like structures. The mechanical properties of the generated models were predictable with numerical simulations as well as an analytical approach. The permeability was found to be better than Sawbones at the same porosity. The models also showed the capability of matching e.g. some vertebral structures for key geometric features.

CONCLUSIONS

An improved numerical model for mimicking trabecular bone structures was successfully developed using Voronoi tessellation and Boolean operations. This is expected to benefit both computational and experimental studies by providing a more diverse and representative structure of trabecular bone.

Study on the association of the microstructure and bone metabolism in the osteoporotic femoral head

BACKGROUND

We compared the bone microstructure and metabolism of the femoral heads in patients with osteoporosis (OP) and non-OP patients to investigate the pathologic mechanism of OP and guide clinical treatment.

METHODS AND RESULTS

From January 2020 to June 2021, we obtained femoral head samples from 30 patients undergoing hip replacement due to femoral neck fracture. All patients were women aged approximately 67 to 80 years (mean age, 74 years). According to the dual-energy X-ray results, the femoral head samples were divided into the OP (T< - 2.5) and non-OP (T > - 1.5) groups. Microcomputed tomography scanning, bone metrology analysis, hematoxylin and eosin staining, and Masson's trichrome staining were used to compare the local bone trabecular microstructure changes. Quantitative reverse transcription PCR was performed to identify changes in the osteogenesis-related genes and the osteoclast-related genes in specific regions to reflect osteogenic and osteoclastic activities. Femoral heads with OP showed significant changes in the local bone microstructure. Bone density, bone volume fraction, and the number and thickness of the bone trabeculae decreased. Local bone metabolism was imbalanced in the areas with microstructural changes in femoral heads with OP, with increased osteoclast activity and decreased osteoblast activity.

CONCLUSIONS

Deterioration of bone microstructure is closely related to abnormal bone metabolism associated with the activity of osteoblasts and osteoclasts in osteoporotic femoral heads. Promoting bone formation by improving local bone metabolism, enhancing osteogenic activity and inhibiting osteoclast activity may be a promising way of preventing local OP and osteoporotic fractures.

A New Microarchitecture-Based Parameter to Predict the Micromechanical Properties of Bone Allografts

Scaffolds are an essential component of bone tissue engineering. They provide support and create a physiological environment for cells to proliferate and differentiate. Bone allografts extracted from human donors are promising scaffolds due to their mechanical and structural characteristics. Bone microarchitecture is well known to be an important determinant of macroscopic mechanical properties, but its role at the microscopic, i.e., the trabeculae level is still poorly understood. The present study investigated linear correlations between microarchitectural parameters obtained from X-ray computed tomography (micro-CT) images of bone allografts, such as bone volume fraction (BV/TV), degree of anisotropy (DA), or ellipsoid factor (EF), and micromechanical parameters derived from micro-finite element calculations, such as mean axial strain (ε_z) and strain energy density (W_e). DAEF, a new parameter based on a linear combination of the two microarchitectural parameters DA and EF, showed a strong linear correlation with the bone mechanical characteristics at the microscopic scale. Our results concluded that the spatial distribution and the plate-and-rod structure of trabecular bone are the main determinants of the mechanical properties of bone at the microscopic level. The DAEF parameter could, therefore, be used as a tool to predict the level of mechanical stimulation at the local scale, a key parameter to better understand and optimize the mechanism of osteogenesis in bone tissue engineering.

A morphological, topological and mechanical investigation of gyroid, spinodoid and dual-lattice algorithms as structural models of trabecular bone

In this study, we evaluate the performance of three algorithms as computational models of trabecular bone architecture, through systematic evaluation of morphometric, topological, and mechanical properties. Here, we consider the widely-used gyroid lattice structure, the recently-developed spinodoid structure and a structure similar to Voronoi lattices introduced here as the dual-lattice. While all computational models were calibrated to recreate the trabecular tissue volume (e.g. BV/TV), it was found that both the gyroid- and spinodoid-based structures showed substantial differences in many other morphometric and topological parameters and, in turn, showed lower effective mechanical properties compared to trabecular bone. The newly-developed dual-lattice structures better captured both morphometric parameters and mechanical properties, despite certain differences being evident their topological configuration compared to trabecular bone. Still, these computational algorithms provide useful platforms to investigate trabecular bone mechanics and for designing biomimetic structures, which could be produced through additive manufacturing for applications that include bone substitutes, scaffolds and porous implants. Furthermore, the software for the creation of the structures has been added to the open source toolbox GIBBON and is therefore freely available to the community.

Evaluation of chronic renal failure with cone beam computed tomography radiomorphometric indices and fractal analysis in the mandible

OBJECTIVES

The purpose of this study is to evaluate the mandibular bone structure of patients with chronic renal failure (CRF) and compare to control group via the fractal analysis (FA) and radiomorphometric indices in the cone beam computed tomography (CBCT) images.

METHODS

Three observers retrospectively investigated 44 CBCT images to compare patients with CRF to healthy controls. FA was performed in two different areas, volume of interests (VOI) were chosen in the mandibular ramus and corpus. The CT cortical index (CTCI), CT mental index (CTMI), and CT mandibular index (CTI) were performed to determine cortical porosity and thickness. The normality distribution of numerical data was tested using the Kolmogorov-Smirnova and Shapiro-Wilk tests. According to the results, the Mann-Whitney U test and independent group t test were used for parameters. The chi-square test was used to evaluate the distribution of categorical variables by groups.

RESULTS

There were statistically significant differences in VOI1 and VOI2. The fractal dimension (FD) values in VOI1 and VOI2 were significantly lower in study group. There were no significant differences in CTCI, CTMI and CTI measurements between both groups.

CONCLUSIONS

CRF is a prevalent cause of radiographic abnormalities in jawbones. The FD values in trabecular bone decreased in study group, although there were no significant differences in the radiomorphometric indices. FA in CBCT images could be useful for a three-dimensional evaluation of trabecular bone structure.

Morphologies in-between: The impact of the first steps on the human talus

OBJECTIVE

The development of bipedalism is a very complex activity that contributes to shaping the anatomy of the foot. The talus, which starts ossifying in utero, may account for the developing stages from the late gestational phase onwards. Here, we explore the early development of the talus in both its internal and external morphology to broaden the knowledge of the anatomical changes that occur during early development.

MATERIALS AND METHODS

The sample consists of high-resolution microCT scans of 28 modern juvenile tali (from 36 prenatal weeks to 2 years), from a broad chronological range from the Late Roman period to the 20th century. We applied geometric morphometric and whole-bone trabecular analysis to investigate the early talar morphological changes.

RESULTS

In the youngest group (<6 postnatal months), the immature external shell is accompanied by an isotropic internal structure, with thin and densely packed trabeculae. After the initial attempts of locomotion, bone volume fraction decreases, while anisotropy and trabecular thickness increase. These internal changes correspond to the maturation of the external shell, which is now more defined and shows the development of the articular surfaces.

DISCUSSION

The internal and external morphology of the human talus reflects the diverse load on the foot during the initial phases of the bipedal locomotion, with the youngest group potentially reflecting the lack of readiness of the human talus to bear forces and perform bipedal walking. These results highlight the link between mechanical loading and bone development in the human talus during the acquisition of bipedalism, providing new insight into the early phases of talar development.

Human talar ontogeny: Insights from morphological and trabecular changes during postnatal growth

Objectives

The study of the development of human bipedalism can provide a unique perspective on the evolution of morphology and behavior across species. To generate new knowledge of these mechanisms, we analyze changes in both internal and external morphology of the growing human talus in a sample of modern human juveniles using an innovative approach.

Materials and Methods

The sample consists of high‐resolution microCT scans of 70 modern juvenile tali, aged between 8 postnatal weeks and 10 years old, from a broad chronological range from Middle/Late Neolithic, that is, between 4800 and 4500 BCE, to the 20th century. We applied geometric morphometric and whole‐bone trabecular analysis (bone volume fraction, degree of anisotropy, trabecular number, thickness, and spacing) to all specimens to identify changes in the external and internal morphology during growth. Morphometric maps were also generated.

Results

During the first year of life, the talus has an immature and globular shape, with a dense, compact, and rather isotropic trabecular architecture, with numerous trabeculae packed closely together. This pattern changes while children acquire a more mature gait, and the talus tends to have a lower bone volume fraction, a higher anisotropy, and a more mature shape.

Discussion

The changes in talar internal and external morphologies reflect the different loading patterns experienced during growth, gradually shifting from an “unspecialized” morphology to a more complex one, following the development of bipedal gait. Our research shows that talar plasticity, even though genetically driven, may show mechanical influences and contribute to tracking the main locomotor milestones.

Ontogenetic Patterning of Human Subchondral Bone Microarchitecture in the Proximal Tibia

High-resolution computed tomography images were acquired for 31 proximal human tibiae, age 8 to 37.5 years, from Norris Farms #36 cemetery site (A.D. 1300). Morphometric analysis of subchondral cortical and trabecular bone architecture was performed between and within the tibial condyles. Kruskal−Wallis and Wilcoxon signed-rank tests were used to examine the association between region, age, body mass, and each morphometric parameter. The findings indicate that age-related changes in mechanical loading have varied effects on subchondral bone morphology. With age, trabecular microstructure increased in bone volume fraction (p = 0.033) and degree of anisotropy (p = 0.012), and decreased in connectivity density (p = 0.001). In the subchondral cortical plate, there was an increase in thickness (p < 0.001). When comparing condylar regions, only degree of anisotropy differed (p = 0.004) between the medial and lateral condyles. Trabeculae in the medial condyle were more anisotropic than in the lateral region. This research represents an innovative approach to quantifying both cortical and trabecular subchondral bone microarchitecture in archaeological remains.

Nonlinear micro finite element models based on digital volume correlation measurements predict early microdamage in newly formed bone

Bone regeneration in critical-sized defects is a clinical challenge, with biomaterials under constant development aiming at enhancing the natural bone healing process. The delivery of bone morphogenetic proteins (BMPs) in appropriate carriers represents a promising strategy for bone defect treatment but optimisation of the spatial-temporal release is still needed for the regeneration of bone with biological, structural, and mechanical properties comparable to the native tissue. Nonlinear micro finite element (μFE) models can address some of these challenges by providing a tool able to predict the biomechanical strength and microdamage onset in newly formed bone when subjected to physiological or supraphysiological loads. Yet, these models need to be validated against experimental data. In this study, experimental local displacements in newly formed bone induced by osteoinductive biomaterials subjected to in situ X-ray computed tomography compression in the apparent elastic regime and measured using digital volume correlation (DVC) were used to validate μFE models. Displacement predictions from homogeneous linear μFE models were highly correlated to DVC-measured local displacements, while tissue heterogeneity capturing mineralisation differences showed negligible effects. Nonlinear μFE models improved the correlation and showed that tissue microdamage occurs at low apparent strains. Microdamage seemed to occur next to large cavities or in biomaterial-induced thin trabeculae, independent of the mineralisation. While localisation of plastic strain accumulation was similar, the amount of damage accumulated in these locations was slightly higher when including material heterogeneity. These results demonstrate the ability of the nonlinear μFE model to capture local microdamage in newly formed bone tissue and can be exploited to improve the current understanding of healing bone and mechanical competence. This will ultimately aid the development of BMPs delivery systems for bone defect treatment able to regenerate bone with optimal biological, mechanical, and structural properties.

Efficacy of bioreactor-activated bone substitute with bone marrow nuclear cells on fusion rate and fusion mass microarchitecture in sheep

Bioreactors have been used for bone graft engineering in pre-clinical investigations over the past 15 years. The ability of bioreactor-incubated bone marrow nuclear cells (BMNCs) to enhance bone-forming potential varies significantly, and the three-dimensional (3D) distribution of BMNCs within the scaffold is largely unknown. The aims of this study were (1) to investigate the efficacy of a carbonated hydroxyapatite (CHA) with/without BMNCs on spine fusion rate and fusion mass microarchitecture using a highly challenging two-level posterolateral spine fusion without instrumentation; and (2) to evaluate 3D distribution of BMNCs within scaffolds characterized by immunohistochemistry. Fusion rate and fusion mass were quantified by micro-CT, microarchitectural analysis, and histology. While the homogenous 3D distribution of BMNCs was not observed, BMNCs were found to migrate towards a substitute core. In the autograft group, the healing rate was 83.3%, irrespective of the presence of BMNCs. In the CHA group, also 83.3% was fused in the presence of BMNCs, and 66.7% fused without BMNCs. A significant decrease in the fusion mass porosity (p = .001) of the CHA group suggested the deposition of mineralized bone. The autograft group revealed more bone, thicker trabeculae, and better trabecular orientation but less connection compared to the CHA group. Immunohistochemistry confirmed the ability of bioreactors to incubate a large-sized substitute coated with viable BMNCs with the potential for proliferation and differentiation. These findings suggested that a bioreactor-activated substitute is comparable to autograft on spine fusion and that new functional bone regeneration could be achieved by a combination of BMNCs, biomaterials, and bioreactors.

Micro-CT-Based Bone Microarchitecture Analysis of the Murine Skull

X-ray micro-computed tomography (micro-CT) imaging has important applications in microarchitecture analysis of cortical and trabecular bone structure. While standardized protocols exist for micro-CT-based microarchitecture assessment of long bones, specific protocols need to be developed for different types of skull bones taking into account differences in embryogenesis, organization, development, and growth compared to the rest of the body. This chapter describes the general principles of bone microarchitecture analysis of murine craniofacial skeleton to accommodate for morphological variations in different regions of interest.

Degenerations in Global Morphometry of Cancellous Bone in Rheumatoid Arthritis, Osteoarthritis and Osteoporosis of Femoral Heads are Similar but More Severe than in Ageing Controls

We have recently revealed significant differences in microarchitectural properties (i.e. global and local morphometries) and mechanical properties between rheumatoid arthritis (RA), osteoarthritis (OA) and osteoporosis (OP) in cancellous bones. This study compared these properties with those of ageing controls by matching bone volume fraction (BV/TV), the most important determinant for bones' mechanical properties, to investigate whether these bones have similar properties and degenerative potentials. RA, OA and OP femoral heads were harvested from patients undergoing total hip replacement surgery. The selected patients were matched by similar cancellous bone BV/TV, with seven patients in each group. Four samples were prepared from each femoral head and scanned with micro-CT to quantify microarchitectural properties and compression tested to determine mechanical properties. In terms of global morphometry, no significant differences were observed between these diseased bones. In terms of local morphometry, the number of plates in the RA group was significantly greater than that of the OP and control groups. Plate volume density in the RA group was significantly greater than in the control group. Interestingly, the ultimate stresses in the three diseased groups were 77% to 195% lower than in the control group (p < 0.001). Degenerations of global morphometry of cancellous bones in these diseased femoral heads are similar but more severe than in ageing controls matched by BV/TV, as evidenced by pronounced reduction in bone strength. This phenomenon suggests that some local morphometric parameters, along with other factors, such as abnormal collagen, mineralisation, erosion and microdamage, may contribute to further compromising mechanical properties.

Determination of anisotropic elastic parameters from morphological parameters of cancellous bone for osteoporotic lumbar spine

In biomechanics, large finite element models with macroscopic representation of several bones or joints are necessary to analyze implant failure mechanisms. In order to handle large simulation models of human bone, it is crucial to homogenize the trabecular structure regarding the mechanical behavior without losing information about the realistic material properties. Accordingly, morphology and fabric measurements of 60 vertebral cancellous bone samples from three osteoporotic lumbar spines were performed on the basis of X-ray microtomography (μCT) images to determine anisotropic elastic parameters as a function of bone density in the area of pedicle screw anchorage. The fabric tensor was mapped in cubic bone volumes by a 3D mean-intercept-length method. Fabric measurements resulted in a high degree of anisotropy (DA = 0.554). For the Young’s and shear moduli as a function of bone volume fraction (BV/TV, bone volume/total volume), an individually fit function was determined and high correlations were found (97.3 ≤ R² ≤ 99.1,p < 0.005). The results suggest that the mathematical formulation for the relationship between anisotropic elastic constants and BV/TV is applicable to current μCT data of cancellous bone in the osteoporotic lumbar spine. In combination with the obtained results and findings, the developed routine allows determination of elastic constants of osteoporotic lumbar spine. Based on this, the elastic constants determined using homogenization theory can enable efficient investigation of human bone using finite element analysis (FEA).

Graphical Abstract

Computed Tomography as a Characterization Tool for Engineered Scaffolds with Biomedical Applications

The ever-growing field of materials with applications in the biomedical field holds great promise regarding the design and fabrication of devices with specific characteristics, especially scaffolds with personalized geometry and architecture. The continuous technological development pushes the limits of innovation in obtaining adequate scaffolds and establishing their characteristics and performance. To this end, computed tomography (CT) proved to be a reliable, nondestructive, high-performance machine, enabling visualization and structure analysis at submicronic resolutions. CT allows both qualitative and quantitative data of the 3D model, offering an overall image of its specific architectural features and reliable numerical data for rigorous analyses. The precise engineering of scaffolds consists in the fabrication of objects with well-defined morphometric parameters (e.g., shape, porosity, wall thickness) and in their performance validation through thorough control over their behavior (in situ visualization, degradation, new tissue formation, wear, etc.). This review is focused on the use of CT in biomaterial science with the aim of qualitatively and quantitatively assessing the scaffolds’ features and monitoring their behavior following in vivo or in vitro experiments. Furthermore, the paper presents the benefits and limitations regarding the employment of this technique when engineering materials with applications in the biomedical field.

Effect of captivity on the vertebral bone microstructure of xenarthran mammals

Captive specimens in museum collections facilitate study of rare taxa, but the lifestyles, diets, and lifespans of captive animals differ from their wild counterparts. Trabecular bone architecture adapts to in vivo forces, and may reflect interspecific variation in ecology and behavior as well as intraspecific variation between captive and wild specimens. We compared trunk vertebrae bone microstructure in captive and wild xenarthran mammals to test the effects of ecology and captivity. We collected μCT scans of the last six presacral vertebrae in 13 fossorial, terrestrial, and suspensorial xenarthran species (body mass: 120 g to 35 kg). For each vertebra, we measured centrum length; bone volume fraction (BV.TV); trabecular number and mean thickness (Tb.Th); global compactness (GC); cross-sectional area; mean intercept length; star length distribution; and connectivity and connectivity density. Wild specimens have more robust trabeculae, but this varies with species, ecology, and pathology. Wild specimens of fossorial taxa (Dasypus) have more robust trabeculae than captives, but there is no clear difference in bone microstructure between wild and captive specimens of suspensorial taxa (Bradypus, Choloepus), suggesting that locomotor ecology influences the degree to which captivity affects bone microstructure. Captive Tamandua and Myrmecophaga have higher BV.TV, Tb.Th, and GC than their wild counterparts due to captivity-caused bone pathologies. Our results add to the understanding of variation in mammalian bone microstructure, suggest caution when including captive specimens in bone microstructure research, and indicate the need to better replicate the habitats, diets, and behavior of animals in captivity.

The local and global geometry of trabecular bone

The organization and shape of the microstructural elements of trabecular bone govern its physical properties, are implicated in bone disease, and serve as blueprints for biomaterial design. To devise fundamental structure-property relationships and design truly bone-mimicking biomaterials, it is essential to characterize trabecular bone structure from the perspective of geometry, the mathematical study of shape. Using micro-CT images from 70 donors at five different sites, we analyze the local and global geometry of human trabecular bone in detail, respectively by quantifying surface curvatures and Minkowski functionals. We find that curvature density maps provide distinct and sensitive shape fingerprints for bone from different sites. Contrary to a common assumption, these curvature maps also show that bone morphology does not approximate a minimal surface but exhibits a much more intricate curvature landscape. At the global (or integral) perspective, our Minkowski analysis illustrates that trabecular bone exhibits other types of anisotropy/ellipticity beyond interfacial orientation, and that anisotropy varies substantially within the trabecular structure. Moreover, we show that the Minkowski functionals unify several traditional morphometric indices. Our geometric approach to trabecular morphometry provides a fundamental language of shape that could be useful for bone failure prediction, understanding geometry-driven tissue growth, and the design of bone-mimicking tissue scaffolds. STATEMENT OF SIGNIFICANCE: The architecture of trabecular bone is key in determining bone properties, and is often a starting point for the design of bone-substitutes. Despite the substantial history of bone morphometry, a fundamental characterization of trabecular bone geometry is still lacking. Therefore, we introduce a robust framework to quantify local and global trabecular bone geometry, which we apply to hundreds of micro-CT scans. Our approach relies on quantifying surface curvatures and Minkowski functionals, which are the most fundamental local and global shape quantifiers. Our results show that these shape metrics are sensitive to differences between bone types and unify traditional metrics within a single mathematical framework. This geometrical framework could also be useful to design bone-mimicking scaffolds and understand geometry-driven tissue growth.

Characterization of Ultralow Density Cellular Solids: Lessons from 30 years of Bone Biomechanics Research

Advances in additive manufacturing techniques have enabled the development of micro-architectured materials displaying a combination of low-density and lightweight structures with high specific strength and toughness. The mechanical performance of micro-architectured materials can be assessed using standard techniques; however, when studying low- and ultralow density micro-architectured materials, standard characterization techniques can be subject to experimental artifacts. Additionally, quantitative assessment and comparisons of microarchitectures with distinct lattice patterns is not always straightforward. Cancellous bone is a natural, ultralow density (porosity often exceeding 90%), irregular, cellular solid that has been thoroughly characterized in terms of micro-architecture and mechanical performance over the past 30 years. However, most the literature on cancellous bone mechanical properties and micro-structure-function relationships is in the medical literature and is not immediately accessible to materials designers. Here we provide a brief review of state-of-the-art approaches for characterizing the micro-architecture and mechanical performance of ultralow density cancellous bone, including methods of addressing experimental artifacts during mechanical characterization of ultralow density cellular solids, methods of quantifying microarchitecture, and currently understood structure-function relationships.

Porous Alumina Ceramics with Multimodal Pore Size Distributions

Pore networks with multimodal pore size distributions combining advantages from isotropic and anisotropic shaped pores of different sizes are highly attractive to optimize the physical properties of porous ceramics. Multimodal porous Al₂O₃ ceramics were manufactured using pyrolyzed cellulose fibers (l = 150 µm, d = 8 µm) and two types of isotropic phenolic resin spheres (d = 30 and 300 µm) as sacrificial templates. The sacrificial templates were homogeneously distributed in the Al₂O₃ matrix, compacted by uniaxial pressing and extracted by a burnout and sintering process up to 1700 °C in air. The amount of sacrificial templates was varied up to a volume content of 67 Vol% to form pore networks with porosities of 0-60 Vol%. The mechanical and thermal properties were measured by 4-point-bending and laser flash analysis (LFA) resulting in bending strengths of 173 MPa to 14 MPa and heat conductivities of 22.5 Wm^-1K^-1 to 4.6 Wm^-1K^-1. Based on µCT-measurements, the representative volume-of-interest (VOI) of the samples digital twin was determined for further analysis. The interconnectivity, tortuosity, permeability, the local and global stress distribution as well as strut and cell size distribution were evaluated on the digital twin's VOI. Based on the experimental and simulation results, the samples pore network can be tailored by changing the fiber to sphere ratio and the overall sacrificial template volume. The presence pore formers significantly influenced the mechanical and thermal properties, resulting in higher strengths for samples containing fibrous templates and lower heat conductivities for samples containing spherical templates.

The Functional and Allometric Implications of Hipbone Trabecular Microarchitecture in a Sample of Eutherian and Metatherian Mammals

The pelvis plays an active role in weight bearing and countering the ground reaction forces incurred by the hindlimbs thus making it a critical component of the locomotor skeleton. Accordingly, this anatomical region is theoretically ideal for inferring locomotor behavior from both external skeletal morphology and trabecular microarchitecture, with the latter possibly offering nuanced insights into the mechanical loading environment given its increased plasticity and higher turnover rate. However, trabecular microarchitecture is also known to be influenced by a variety of factors including body size, sex, age, genetic regulation, diet and activity level, that collectively hinder the ability to generate consistent functional inferences. In this study, a comparative sample of mammals (42 species spanning four orders) of varying sizes, yet comparable locomotor repertoires, were evaluated to determine the effects of body size, phylogeny and locomotion on hipbone trabecular microarchitecture. This study found a weak functional signal detected in differences in bone volume fraction and the degree of anisotropy across certain pre-assigned locomotor categories, while confirming previously recognized allometric scaling trends reported for other mammalian samples based on the femur. Within primates, a more anisotropic pattern was observed for quadrupedal species attributed to their repetitive loading regimes and stereotypical limb excursions, while isotropic values were revealed for taxa utilizing more varied arboreal repertoires. Humans, despite a frequent and predictable loading environment associated with their use of bipedalism, showed relatively isotropic values. This study highlights the confounding factors that influence trabecular microarchitecture and consequently limit its utility as a method for investigating locomotor adaptation.

Trabecular bone properties in the ilium of the Middle Paleolithic/Middle Stone Age Border Cave 3 Homo sapiens infant and the onset of independent gait

The Border Cave 3 (BC3) infant skeleton has been understudied, despite its importance as an example of a well-preserved and fairly complete immature skeleton of early Homo sapiens which potentially provides a rare window into various aspects of ontogenetic development, including locomotor activity (e.g., timing of gait events). Trabecular structure in the BC3 ilium was evaluated to investigate whether it matches that of an equivalently aged infant from a postindustrialized society. Microcomputed tomography (μCT) scans were acquired from the BC3 infant and from an ontogenetic series of 25 postindustrial infants that were divided into three age classes (ACs) ranging from neonates to toddlers (<36 months). All ilia were qualitatively compared and then digitally subdivided into 10 volumes of interest (VOIs) based on anatomical reference points. The VOIs were quantified and ontogenetic differences in trabecular structure were statistically evaluated. Across the comparative ontogenetic series, trabecular architectural properties overlapped in all regions. However, trabecular thickness increased significantly after the first year of life. The BC3 infant demonstrated generally similar trabecular structure to that observed in the age-equivalent postindustrial infants (AC2), including relatively strong development of the trabecular chiasma qualitatively. However, some interesting distinctions were observed in BC3, such as low strut thickness compared with infants from the postindustrial sample, that bear further exploration in future studies. Evaluation of only one individual from the Middle Stone Age (MSA), coupled with the relatively small comparative sample, limit our ability to distinguish more meaningful biological differences in trabecular structure throughout ontogeny from idiosyncratic characteristics. Nonetheless, results of this study extend ongoing research on infant locomotor and morphological development to archeological populations in the Middle Stone Age. Further cross-cultural studies consisting of larger comparative postindustrial samples may provide additional information on trabecular structure in the infant ilium during this important developmental timeframe.

Novel strategies for the characterization of cancellous bone morphology: Virtual isolation and analysis

OBJECTIVES

The advent of micro-computed tomography (μCT) made cancellous bone more accessible than ever before. Nevertheless, the characterization of cancellous bone is made difficult by its inherent complexity and the difficulties in defining homology across datasets. Here we propose novel virtual methodological approaches to overcome those issues and complement existing methods.

MATERIALS AND METHODS

We present a protocol for the isolation of the whole cancellous region within a μCT scanned bone. This method overcomes the subsampling issues and allows studying cancellous bone as a single unit. We test the protocol on a set of primate bones. In addition, we describe a set of morphological indices calculated on the topological skeleton of the cancellous bone: node density, node connectivity, trabecular angle, trabecular tortuosity, and fractal dimension. The usage of the indices is shown on a small comparative sample of primate femoral heads.

RESULTS

The isolation protocol proves reliable in isolating cancellous structures from several different bones, regardless of their shape. The indices seem to detect some functional differences, although further testing on comparative samples is needed to clarify their potential for the study of cancellous architecture.

CONCLUSIONS

The approaches presented overcome some of the difficulties of trabecular bone studies. The methods presented here represent an alternative or supporting method to the existing tools available to address the biomechanics of cancellous bone.

Comparative microcomputed tomographic structural analysis of the trabecular and cortical bone architecture of radius and ulna in toy dog breeds

Radius and ulna fractures can be detected in dogs with an incidence of 18.0% and occur more frequently in toy breeds with increased complication rates. The reasons for the predisposition and the increased complication rate of antebrachial fractures in toy breeds are not yet fully understood. The aim of this study was to research the associated mechanisms by investigating the trabecular structure parameters (BV/TV (%), Tb.N (mm^-1), Tb.Th (mm), Tb.Sp (mm), Conn. D (mm^-3), DA) and the cortical bone density (Mean Density (mg HA/ccm)) for toy breeds and small dogs of other breeds using microcomputed tomography and to identify the differences between both groups. A total of 70 forearms, originating from toy breeds (n = 36) and small dogs of other breeds (n = 34) were evaluated by microcomputed tomography. The results of this study showed a lower bone volume fraction (BV/TV), more (Tb.N) but thinner trabeculae (Tb.Th) with greater separation (Tb.Sp) and decreased diaphyseal cortical density (Mean Density) in toy breed dogs. This could thus contribute to the increased risk of antebrachial fractures in toy breeds. Within the group of toy breeds inhomogeneous results were observed, whereby the values of the Chihuahuas deviated. This breed had significantly thinner trabeculae (Tb.Th) with a greater trabecular separation (Tb.Sp), lower bone volume fraction (BV/TV) and significantly lower diaphyseal cortical density (Mean Density). These parameters could therefore indicate an increased antebrachial fracture and complication risk, especially in Chihuahuas.

Selective Voronoi tessellation as a method to design anisotropic and biomimetic implants

The geometry of a metallic scaffold is important for the success of bone implants, where the introduction of porosity can reduce stress shielding effects and allow for bone tissue integration. In this work, porous scaffolds were designed to closely mimic the natural structure of trabecular bone using selective Voronoi tessellation with preferential seeding. A workflow to generate these structures is introduced, where voided regions of seeds in the starting volume create preferential texture during polyhedral expansion, resulting in modified strut orientation in the implant. Anisotropy was digitally characterized by mean-intercept length and star volume distribution measurements to determine similarity to trabecular orientation. This work demonstrates that selective Voronoi tessellation is an effective method to generate biomimetic porous scaffolds with increased anisotropy and tunable strut architecture in three dimensions as a suitable alternative to patient-derived bone geometries.

Comparison of Trabecular Bone Mineral Density Measurement Using Hounsfield Unit and Trabecular Microstructure in Orthodontic Patients Using Cone-Beam Computed Tomography

The aim of this study was to measure the bone mineral density of specific regions of maxilla, mandible, and first cervical vertebra using the Hounsfield unit and trabecular microstructure pattern analysis and to compare the two methods. In this study, cone-beam computed tomography (CBCT) images were obtained from 58 patients. Trabecular thickness, trabecular number, trabecular separation, and bone volume fraction were measured in 484 regions for trabecular microstructure parameters and Hounsfield unit was measured for the grayscale value. There was no difference in bone mineral density between the right and left side in every site and between males and females. Trabecular thickness and trabecular number were high in the order of anterior base of the maxilla, mandibular body, first cervical vertebra, and mandibular condyle. Bone volume fraction and Hounsfield unit were high in the order of anterior base of the maxilla, mandibular body, mandibular condyle, and first cervical vertebra (p < 0.05). Trabecular thickness, trabecular number, and bone volume fraction was positively correlated to the Hounsfield unit, and trabecular separation was negatively correlated to the Hounsfield unit (p < 0.005). This study suggests that it is possible to compare the bone mineral density of trabecular bone in various sites using the Hounsfield unit and trabecular microstructure pattern analysis.

Assessing patterns of variation in BV/TV in the calcaneus and C2 vertebra of Gorilla gorilla, Pan troglodytes, and populations of Homo sapiens from the Pleistocene and Holocene that differ in physical activity levels

OBJECTIVES

Because trabecular bone volume fraction (BV/TV) is influenced by variations in physical activity recent declines in BV/TV in humans are often attributed to modern sedentary lifestyles. This study tests the hypothesis that presumed variations in mechanical loading between groups can predict the observed BV/TV patterns in humans, chimpanzees and gorillas in two bones: the calcaneus which experiences high and well characterized impact forces, and the C2 vertebrae which experiences reduced locomotor forces.

MATERIALS AND METHODS

BV/TV and other structural variables were quantified from high-resolution microCT scans in gorillas, chimpanzees, and four Homo sapiens populations: Pleistocene, semi-sedentary Natufians; Holocene hunter-gatherers from Point Hope, Alaska; Holocene nomadic pastoralists from medieval Europe; and modern, sedentary Americans.

RESULTS

In the calcaneal tuberosity, Natufian BV/TV was 36, 46, and 46% greater than Alaskans (p = .02), Europeans (p = .005) and modern Americans (p = .002), respectively, but not significantly different from apes. BV/TV was not significantly different between modern Americans and Alaskans or Europeans. In the C2, Natufian BV/TV was 53 and 25% greater than in the Alaskan (p = .0001) and European (p = .048) populations.

DISCUSSION

These results suggest that phenomena other than or in addition to variations in physical activity are needed to explain BV/TV patterns observed in H. sapiens, and point to a systemic decline in H. sapiens BV/TV after the Pleistocene.

Discrimination of Low-Energy Acetabular Fractures from Controls Using Computed Tomography-Based Bone Characteristics

The incidence of low-energy acetabular fractures has increased. However, the structural factors for these fractures remain unclear. The objective of this study was to extract trabecular bone architecture and proximal femur geometry (PFG) measures from clinical computed tomography (CT) images to (1) identify possible structural risk factors of acetabular fractures, and (2) to discriminate fracture cases from controls using machine learning methods. CT images of 107 acetabular fracture subjects (25 females, 82 males) and 107 age-gender matched controls were examined. Three volumes of interest, one at the acetabulum and two at the femoral head, were extracted to calculate bone volume fraction (BV/TV), gray-level co-occurrence matrix and histogram of the gray values (GV). The PFG was defined by neck shaft angle and femoral neck axis length. Relationships between the variables were assessed by statistical mean comparisons and correlation analyses. Bayesian logistic regression and Elastic net machine learning models were implemented for classification. We found lower BV/TV at the femoral head (0.51 vs. 0.55, p = 0.012) and lower mean GV at both the acetabulum (98.81 vs. 115.33, p < 0.001) and femoral head (150.63 vs. 163.47, p = 0.005) of fracture subjects when compared to their matched controls. The trabeculae within the femoral heads of the acetabular fracture sides differed in structure, density and texture from the corresponding control sides of the fracture subjects. Moreover, the PFG and trabecular architectural variables, alone and in combination, were able to discriminate fracture cases from controls (area under the receiver operating characteristics curve 0.70 to 0.79). In conclusion, lower density in the acetabulum and femoral head with abnormal trabecular structure and texture at the femoral head, appear to be risk factors for low-energy acetabular fractures.

Numerical analysis of hip fracture due to a sideways fall

The primary purpose of this paper is to outline a methodology for evaluating the likelihood of cortical bone fracture in the proximal femur in the event of a sideways fall. The approach includes conducting finite element (FE) analysis in which the cortical bone is treated as an anisotropic material, and the admissibility of the stress field is validated both in tension and compression regime. In assessing the onset of fracture, two methodologies are used, namely the Critical Plane approach and the Microstructure Tensor approach. The former is employed in the tension regime, while the latter governs the conditions at failure in compression. The propagation of localized damage is modeled using a constitutive law with embedded discontinuity (CLED). In this approach, the localized deformation is described by a homogenization procedure in which the average properties of cortical tissue intercepted by a macrocrack are established. The key material properties governing the conditions at failure are specified from a series of independent material tests conducted on cortical bone samples tested at different orientations relative to the loading direction. The numerical analysis deals with simulations of experiments involving the sideways fall, and the results are compared with the experimental data. This includes both the evolution of fracture pattern and the local load-displacement characteristics. The proposed approach is numerically efficient, and the results do not display a pathological mesh-dependency. Also, in contrast to the XFEM approach, the analysis does not require any extra degrees of freedom.

Mechanical and morphological characterization of PMMA/bone composites in human femoral heads

PMMA bone cement has gained an important place in a variety of orthopaedic applications in the femur. However, appropriate data on the mechanical properties of bone-cement composites from the human femur are lacking. Therefore, the goal of this study was to determine the morphological and quasi-static compressive properties of proximal femoral bone-cement composites. Thirty trabecular bone specimens were extracted from fifteen pairs of human femoral heads using specimen-specific cutting guides to ensure an accurate alignment with the main trabecular direction (MTD). One specimen from each pair was augmented with PMMA bone cement, while the other one was left untreated. Specimens were scanned with μCT to determine morphological parameters and tested in quasi-static compression until failure. We found that the long axis of the specimens was highly aligned with the MTD (mean error < 5°). A higher compressive modulus and ultimate strength were observed for the bone-cement composite specimens (E = 5.7 ± 0.4 GPa; σ_u = 77.9 ± 5.1 MPa) compared to the bone only specimens (E = 2.9 ± 0.7 GPa; σ_u = 19.0 ± 5.8 MPa). Furthermore, the composites had a higher modulus, but lower strength than cement itself (E = 5.0 ± 0.3 GPa; σ_u = 85.9 ± 2.7 MPa) and the composite modulus was significantly correlated with the bone volume fraction (BV/TV). These results are in contrast to previous findings on human vertebral bone, where the composite was more compliant than cement and no correlation was found between BV/TV and the composite modulus. Thus, properties of bone-cement composites cannot simply be applied across different anatomical sites; the site-specific differences in bone density and trabecular alignment should be taken into account. Collectively, the present results suggest that at low BV/TV, cement dominates the composite properties, while at high BV/TV, the contribution of bone becomes apparent, revealing a positive relationship between BV/TV and the on-axis modulus.

Porous polylactic acid scaffolds for bone regeneration: A study of additively manufactured triply periodic minimal surfaces and their osteogenic potential

Three different triply periodic minimal surfaces (TPMS) with three levels of porosity within those of cancellous bone were investigated as potential bone scaffolds. TPMS have emerged as potential designs to resemble the complex mechanical and mass transport properties of bone. Diamond, Schwarz, and Gyroid structures were 3D printed in polylactic acid, a resorbable medical grade material. The 3D printed structures were investigated for printing feasibility, and assessed by morphometric studies. Mechanical properties and permeability investigations resulted in similar values to cancellous bone. The morphometric analyses showed three different patterns of pore distribution: mono-, bi-, and multimodal pores. Subsequently, biological activity investigated with pre-osteoblastic cell lines showed no signs of cytotoxicity, and the scaffolds supported cell proliferation up to 3 weeks. Cell differentiation investigated by alkaline phosphatase showed an improvement for higher porosities and multimodal pore distributions, suggesting a higher dependency on pore distribution and size than the level of interconnectivity.

Misalignment Error in Cancellous Bone Apparent Elastic Modulus Depends on Bone Volume Fraction and Degree of Anisotropy

Cancellous bone is an anisotropic structure with architectural and mechanical properties that vary due to both skeletal site and disease state. This anisotropy means that, in order to accurately and consistently measure the mechanical properties of cancellous bone, experiments should be performed along the primary mechanical axis (PMA), that is, the orientation in which the mechanical properties are at their maximum value. Unfortunately, some degree of misalignment will always be present, and the magnitude of the resulting error is expected to be architecture dependent. The goal of this work is to quantify the dependence of the misalignment error, expressed in terms of change in apparent elastic modulus (ΔE), on both the bone volume fraction (BV/TV) and the degree of anisotropy (DA). Finite element method (FEM) models of bovine cancellous bone from five different skeletal sites were created at 5 deg and 20 deg from the PMA determined for each region. An additional set of models was created using image dilation/erosion steps in order to control for BV/TV and better isolate the effect of DA. Misalignment error was found to increase with increasing DA and decreasing BV/TV. At 5 deg misaligned from the PMA, error is relatively low (&lt;5%) in all cases but increases to 8–24% error at 20 deg. These results suggest that great care is needed to avoid introducing misalignment error into experimental studies, particularly when studying regions with high anisotropy and/or low bone volume fraction, such as vertebral or osteoporotic bone.

Ultrasounds could be considered as a future tool for probing growing bone properties

Juvenile bone growth is well described (physiological and anatomical) but there are still lacks of knowledge on intrinsic material properties. Our group has already published, on different samples, several studies on the assessment of intrinsic material properties of juvenile bone compared to material properties of adult bone. The purpose of this study was finally to combine different experimental modalities available (ultrasonic measurement, micro-Computed Tomography analysis, mechanical compression tests and biochemical measurements) applied on small cubic bone samples in order to gain insight into the multiparametric evaluation of bone quality. Differences were found between juvenile and adult groups in term of architectural parameters (Porosity Separation), Tissue Mineral Density (TMD), diagonal stiffness coefficients (C₃₃, C_44, C_55, C₆₆) and ratio between immature and mature cross-links (CX). Diagonal stiffness coefficients are more representative of the microstructural and biochemical parameters of child bone than of adult bone. We also found that compression modulus E was highly correlated with several microstructure parameters and CX in children group while it was not at all correlated in the adult group. Similar results were found for the CX which was linked to several microstructure parameters (TMD and E) only in the juvenile group. To our knowledge, this is the first time that, on a same sample, ultrasonic measurements have been combined with the assessment of mechanical and biochemical properties. It appears that ultrasonic measurements can provide relevant indicators of child bone quality (microstructural and biochemical parameters) which is promising for clinical application since, B-mode ultrasound is the preferred first-line modality over other more constraining imaging modalities (radiation, parent–child accessibility and access to the patient's bed) for pediatric patients.

Is there any predictive bone parameter for implant stability in 2-dimensional and 3-dimensional radiologic images?

OBJECTIVES

This ex vivo study aimed to compare radiomorphometric parameters between 2-dimensional (2-D) and 3-dimensional (3-D) radiographs and evaluate the influence of preoperative radiologic bone parameters on the clinical outcomes of implant stability.

STUDY DESIGN

Implant recipient sites in fresh bovine blocks were evaluated on panoramic radiographs for gray value (GV), fractal dimension (FD), number of connected trabeculae (Co), and density of connected trabeculae (CoD). Cone beam computed tomography (CBCT) scans were evaluated for trabecular thickness (TbTh), cortical thickness (CTh), degree of anisotropy (DA), FD, and Co. Insertion torque (IT) and implant stability quotient (ISQ) were measured.

RESULTS

GV was significantly correlated with all parameters in 2-D and 3-D images except FD in 2-D and Co in 3-D, and with all surgical parameters (P ≤ .029). Co and CoD values on panoramic radiographs had significant correlation with TbTh, CTh, and DA values on CBCT images (P < .001). All 2-D parameters and TbTh and CTh in the CBCT data were significantly correlated with IT only (P ≤ .047). Only GV was correlated with ISQ measurements (P ≤ .029).

CONCLUSIONS

GV, Co, and CoD values on panoramic radiographs reflect the architecture of trabecular bone and the thickness of cortical bone, and might help predict implant stability in clinical situations.

Sub-trabecular strain evolution in human trabecular bone

To comprehend the most detrimental characteristics behind bone fractures, it is key to understand the material and tissue level strain limits and their relation to failure sites. The aim of this study was to investigate the three-dimensional strain distribution and its evolution during loading at the sub-trabecular level in trabecular bone tissue. Human cadaver trabecular bone samples were compressed in situ until failure, while imaging with high-resolution synchrotron radiation X-ray tomography. Digital volume correlation was used to determine the strains inside the trabeculae. Regions without emerging damage were compared to those about to crack. Local strains in close vicinity of developing cracks were higher than previously reported for a whole trabecular structure and similar to those reported for single isolated trabeculae. Early literature on bone fracture strain thresholds at the tissue level seem to underestimate the maximum strain magnitudes in trabecular bone. Furthermore, we found lower strain levels and a reduced ability to capture detailed crack-paths with increased image voxel size. This highlights the dependence between the observed strain levels and the voxel size and that high-resolution is needed to investigate behavior of individual trabeculae. Furthermore, low trabecular thickness appears to be one predictor of developing cracks. In summary, this study investigated the local strains in whole trabecular structure at sub-trabecular resolution in human bone and confirmed the high strain magnitudes reported for single trabeculae under loading and, importantly extends its translation to the whole trabecular structure.

Using point clouds to investigate the relationship between trabecular bone phenotype and behavior: An example utilizing the human calcaneus

OBJECTIVES

The objective of this study is to demonstrate a new method for analyzing trabecular bone volume fraction and degree of anisotropy in three dimensions.

METHODS

We use a combination of automatic mesh registration, point-cloud correspondence registration, and P-value corrected univariate statistical tests to compare bone volume fraction and degree of anisotropy on a point by point basis across the entire calcaneus of two human groups with different subsistence strategies.

RESULTS

We found that the patterns of high and low bone volume fraction and degree of anisotropy distribution between the Black Earth (hunter-gatherers) and Norris Farms (mixed-strategy agriculturalists) are very similar, but differ in magnitude. The hunter-gatherers exhibit higher levels of bone volume fraction and less anisotropic trabecular bone organization. Additionally, patterns of bone volume fraction and degree of anisotropy in the calcaneus correspond well with biomechanical expectations of relative forces experienced during walking and running.

CONCLUSIONS

We conclude that comparing site-specific, localized differences in trabecular bone variables such as bone volume fraction and degree of anisotropy in three-dimensions is a powerful analytical tool. This method makes it possible to determine where similarities and differences between groups are located within the whole skeletal element of interest. The visualization of multiple variables also provides a way for researchers to see how the trabecular bone variables interact within the morphology, and allows for a more nuanced understanding of how they relate to one another and the broader mechanical environment.

Influence of nano-hydroxyapatite coating implants on gene expression of osteogenic markers and micro-CT parameters. An in vivo study in diabetic rats

This study evaluated the response of a nano-hydroxyapatite coating implant through gene expression analysis (runt-related transcription factor 2 (Runx2), alkaline phosphatase (Alp), osteopontin (Opn), osteocalcin (Oc), receptor activator of nuclear factor-kappa B (Rank), receptor activator of nuclear factor-kappa B ligand (Rank-L), and osteoprotegerin (Opg)). Three-dimensional evaluation (percent bone volume (BV/TV); percent intersection surface (BIC); bone surface/volume ratio (BS/BV); and total porosity (To.Po)) were also analyzed. Mini implants were surgically placed in tibias of both healthy and diabetic rats. The animals were euthanized at 7 and 30 days. Evaluating all factors the relative expression of Rank showed that NANO surface presented the best results at 7 days (diabetic rats). Furthermore the levels of Runx2, Alp, Oc, and Opn suggest an increase in osteoblasts proliferation, especially in early stages of osseointegration. %BIC in healthy and diabetic (7 days) depicted statistically significant differences for NANO group. BV/TV, BS/BV and To.Po demonstrated higher values for NANO group in all evaluated time point and irrespective of systemic condition, but BS/BV 30 days (healthy rat) and 7 and 30 days (diabetic rat). Microtomographic and gene expression analyses have shown the benefits of nano-hydroxyapatite coated implants in promoting new bone formation in diabetic rats.

Investigation of the "Superior Facet Rule" Using 3D-Printed Thoracic Vertebrae With Simulated Corticocancellous Interface

BACKGROUND

Pedicle screw placement is the most common method of fixation in the thoracic spine. Use of the "superior facet rule" allows the operator to locate the borders of the pedicle reliably using posterior landmarks alone. This study investigated the ability of 3-dimensionally (3D)-printed thoracic vertebrae, made from combined thermoplastic polymers, to demonstrate pedicle screw cannulation accurately using the superior facet as a reliable landmark.

METHODS

An anonymized computed tomography scan of the thoracic spine was obtained. The T1-T12 thoracic vertebrae were anatomically segmented and 3D-printed. The pedicle diameters and distance from the midpoint of the superior facet to the ventral lamina were recorded. A total of 120 thoracic pedicles in 60 thoracic vertebral models were instrumented using a freehand technique based only on posterior landmarks. The vertebral models were then coronally cut and examined for medial or lateral violations of the pedicle after screw placement.

RESULTS

A total of 120 pedicle screws were placed successfully within the 3D-printed thoracic vertebral models. Average measurements fell within 1 standard deviation of previous population studies. There were no pedicle wall violations using standard posterior element landmarks for instrumentation. There were 3 lateral violations of the vertebral body wall during screw placement, all attributable to the insertion technique.

CONCLUSIONS

3D-printed thoracic vertebral models using combined thermoplastic polymers can accurately demonstrate the anatomical ultrastructure and posterior element relationships of the superior facet rule for safe thoracic pedicle screw placement. This method of vertebral model prototyping could prove useful for surgical education and demonstrating spinal anatomy.

Comparison of different microCT-based morphology assessment tools using human trabecular bone

MicroCT-based morphological parameters are often used to quantify the structural properties of trabecular bone. Various software tools are available for calculating these parameters. Studies that examine the comparability of their results are rare. Four different software tools were used to analyse a set of 701 microCT images from human trabecular bone samples. Bone volume to total volume (BV/TV), bone surface (BS), trabecular thickness (Tb. Th.) and degree of anisotropy (DA) were evaluated. BV/TV shows very low difference (-0.18 ± 0.15%). The difference in BS could be reduced below 5% if artificial cut surfaces are not included. Tb. Th. and Tb. Sp. show differences of maximal -12% although the same theoretical background is used. DA is most critical with differences from 4.75 ± 3.70% (medtool vs. Scanco), over -38.61 ± 13.15% (BoneJ vs. Scanco), up to 80.52 ± 50.04% (medtool vs. BoneJ). Quantitative results should be considered with caution, especially when comparing different studies. Introducing standardization procedures and the disclosure of underlying algorithms and their respective implementations could improve this issue.

Evaluation of Structural Anisotropy in a Porous Titanium Medium Mimicking Trabecular Bone Structure Using Mode-Converted Ultrasonic Scattering

The mode-converted (longitudinal to transverse, L-T) ultrasonic scattering method was utilized to characterize the structural anisotropy of a phantom mimicking the structural properties of trabecular bone. The sample was fabricated using metal additive manufacturing from high-resolution computed tomography (CT) images of a sample of trabecular horse bone with strong anisotropy. Two focused transducers were used to perform the L-T ultrasonic measurements. A normal incidence transducer was used to transmit longitudinal ultrasonic waves into the sample, while the scattered transverse signals were received by an oblique incidence transducer. At multiple locations on the sample, four L-T measurements were performed by collecting ultrasonic scattering from four directions. The amplitude of the root mean square (rms) of the collected ultrasonic scattering signals was calculated for each L-T measurement. The ratios of rms amplitudes for L-T measurements in different directions were calculated to characterize the anisotropy of sample. The results show that the amplitude of L-T converted scattering is highly dependent on the direction of microstructural anisotropy. A strong anisotropy of the microstructure was observed, which coincides with simulation results previously published on the same structure as well as with the anisotropy estimated from the CT images. These results suggest the potential of mode-converted ultrasonic scattering methods to assess the anisotropy of materials with porous, complex structures, including trabecular bone.

Relationships Among Bone Morphological Parameters and Mechanical Properties of Cadaveric Human Vertebral Cancellous Bone

Mechanical properties and morphological features of the vertebral cancellous bone are related to resistance to fracture and capability of withstanding surgical treatments. In particular, vertebral strength is related to its elastic properties, whereas the ease of fluid motion, related to the success of incorporation orthopedic materials (eg, bone cement), is regulated by the hydraulic permeability (K). It has been shown that both elastic modulus and permeability of a material are affected by its morphology. The objective of this study was to establish relations between local values of K and the aggregate modulus (H), and parameters descriptive of the bone morphology. We hypothesized that multivariate statistical models, by including the contribution of several morphology parameters at once, would provide a strong correlation with K and H of the vertebral cancellous bone. Hence, μCT scans of human lumbar vertebra were used to determine a set of bone morphology descriptors. Subsequently, indentation tests on the bone samples were conducted to determine local values of K and H. Finally, a multivariate approach supported by principal component analysis was adopted to develop predictive statistical models of bone permeability and aggregate modulus as a function of bone morphology descriptors. It was found that linear combinations of bone volume fraction, trabecular thickness, trabecular spacing, structure model index, connectivity density, and degree of anisotropy provide a strong correlation (R 2 ~ 76%) with K and a weaker correlation (R 2 ~ 47%) with H. The results of this study can be exploited in computational mechanics frameworks for investigating the potential mechanical behavior of human vertebra and to develop strategies to treat or prevent pathological conditions such as osteoporosis, age-related bone loss, and vertebral compression fractures. © 2020 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.