Scaling VOI size in 3D μCT studies of trabecular bone: a test of the over-sampling hypothesis

Technical note: Does scan resolution or downsampling impact the analysis of trabecular bone architecture?

The "gold standard" for the assessment of trabecular bone structure is high-resolution micro-CT. In this technical note, we test the influence of initial scan resolution and post hoc downsampling on the quantitative and qualitative analysis of trabecular bone in a Gorilla tibia. We analyzed trabecular morphology in the right distal tibia of one Gorilla gorilla individual to investigate the impact of variation in voxel size on measured trabecular variables. For each version of the micro-CT volume, trabecular bone was segmented using the medical image analysis method. Holistic morphometric analysis was then used to analyze bone volume (BV/TV), anisotropy (DA), trabecular thickness (Tb.Th), spacing (Tb.Sp), and number (Tb.N). Increasing voxel size during initial scanning was found to have a strong impact on DA and Tb.Th measures, while BV/TV, Tb.Sp, and Tb.N were found to be less sensitive to variations in initial scan resolution. All tested parameters were not substantially influenced by downsampling up to 90 μm resolution. Color maps of BV/TV and DA also retained their distribution up to 90 μm. This study is the first to examine the effect of variation in micro-CT voxel size on the analysis of trabecular bone structure using whole epiphysis approaches. Our results indicate that microstructural variables may be measured for most trabecular parameters up to a voxel size of 90 μm for both scan and downsampled resolutions. Moreover, if only BV/TV, Tb.Sp or Tb.N is measured, even larger voxel sizes might be used without substantially affecting the results.

Comparing lateral plantar process trabecular structure to other regions of the human calcaneus

Investigating skeletal adaptations to bipedalism informs our understanding of form-function relationships. The calcaneus is an important skeletal element to study because it is a weight-bearing bone with a critical locomotor role. Although other calcaneal regions have been well studied, we lack a clear understanding of the functional role of the lateral plantar process (LPP). The LPP is a bony protuberance on the inferolateral portion of the calcaneus thought to aid the tuberosity in transmission of ground reaction forces during heel-strike. Here, we analyze LPP internal trabecular structure relative to other calcaneal regions to investigate its potential functional affinities. Human calcanei (n = 20) were micro-CT scanned, and weighted spherical harmonic analysis outputs were used to position 251 volumes of interest (VOI) within each bone. Trabecular thickness (Tb.Th), spacing (Tb.Sp), degree of anisotropy (DA), and bone volume fraction (BV/TV) were calculated for each VOI. Similarities in BV/TV and DA (p = 0.2741) between the LPP and inferior tuberosity support suggestions that the LPP is a weight-bearing structure that may transmit forces in a similar direction. The LPP significantly differs from the inferior tuberosity in Tb.Th and Tb.Sp (p < 0.05). Relatively thinner, more closely spaced trabeculae in the LPP may serve to increase internal surface area to compensate for its relatively small size compared to the tuberosity. Significant differences in all parameters between LPP and joint articular surfaces indicate that trabecular morphology is differently adapted for the transmission of forces associated with body mass through joints.

Definition of the Region of Interest for the Assessment of Alveolar Bone Repair Using Micro-computed Tomography

The objective of this study was to evaluate the influence of the extraction socket (distal or lingual root) and the type of region of interest (ROI) definition (manual or predefined) on the assessment of alveolar repair following tooth extraction using micro-computed tomography (micro-CT). The software package used for scanning, reconstruction, reorientation, and analysis of images (NRecon®, DataViewer®, CT-Analyzer®) was acquired through Bruker <  https://www.bruker.com  > . The sample comprised the micro-CT volumes of seven Wistar rat mandibles, in which the right first molar was extracted. The reconstructed images were analyzed using the extraction sockets, i.e., the distal and intermediate lingual root and the method of ROI definition: manual (MA), central round (CR), and peripheral round (PR). The bone volume fraction (BV/TV) values obtained were analyzed by two-way ANOVA with Tukey's post hoc test (α = 5%). The distal extraction socket resulted in significantly lower BV/TV values than the intermediate lingual socket for MA (P = 0.001), CR (P < 0.001), and PR (P < 0.001). Regarding the ROI, when evaluating the distal extraction socket, the BV/TV was significantly higher (P < 0.001) for MA than for CR and PR, with a lower BV/TV for CR. However, no significant difference was observed for MA (P = 0.855), CR (P = 0.769), or PR (P = 0.453) in the intermediate lingual extraction socket. The bone neoformation outcome (BV/TV) for alveolar bone repair after tooth extraction is significantly influenced by the ROI and the extraction socket. Using the predefined method with a standardized ROI in the central region of the distal extraction socket resulted in the assessment of bone volume, demonstrating the most critical region of the bone neoformation process.

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.

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.

Differing effects of size and lifestyle on bone structure in mammals

BACKGROUND

Mammals are a highly diverse group, with body mass ranging from 2 g to 170 t, and encompassing species with terrestrial, aquatic, aerial, and subterranean lifestyles. The skeleton is involved in most aspects of vertebrate life history, but while previous macroevolutionary analyses have shown that structural, phylogenetic, and functional factors influence the gross morphology of skeletal elements, their inner structure has received comparatively little attention. Here we analysed bone structure of the humerus and mid-lumbar vertebrae across mammals and their correlations with different lifestyles and body size.

RESULTS

We acquired bone structure parameters in appendicular and axial elements (humerus and mid-lumbar vertebra) from 190 species across therian mammals (placentals + marsupials). Our sample captures all transitions to aerial, fully aquatic, and subterranean lifestyles in extant therian clades. We found that mammalian bone structure is highly disparate and we show that the investigated vertebral structure parameters mostly correlate with body size, but not lifestyle, while the opposite is true for humeral parameters. The latter also show a high degree of convergence among the clades that have acquired specialised (non-terrestrial) lifestyles.

CONCLUSIONS

In light of phylogenetic, size, and functional factors, the distribution of each investigated structural parameter reveals patterns explaining the construction of appendicular and axial skeletal elements in mammalian species spanning most of the extant diversity of the clade in terms of body size and lifestyle. These patterns should be further investigated with analyses focused on specific lifestyle transitions that would ideally include key fossils.

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.

Trabecular organization of the proximal femur in Paranthropus robustus: Implications for the assessment of its hip joint loading conditions

Reconstruction of the locomotor repertoire of the australopiths (Australopithecus and Paranthropus) has progressively integrated information from the mechanosensitive internal structure of the appendicular skeleton. Recent investigations showed that the arrangement of the trabecular network at the femoral head center is biomechanically compatible with the pattern of cortical bone distribution across the neck, both suggesting a full commitment to bipedalism in australopiths, but associated with a slightly altered gait kinematics compared to Homo involving more lateral deviation of the body center of mass over the stance limb. To provide a global picture in Paranthropus robustus of the trabecular architecture of the proximal femur across the head, neck and greater trochanter compartments, we applied techniques of virtual imaging to the variably preserved Early Pleistocene specimens SK 82, SK 97, SK 3121, SKW 19 and SWT1/LB-2 from the cave site of Swartkrans, South Africa. We also assessed the coherence between the structural signals from the center of the head and those from the trabecular network of the inferolateral portion of the head and the inferior margin of the neck, sampling the so-called vertical bundle, which in humans represents the principal compressive system of the joint. Our analyses show a functionally related trabecular organization in Pa. robustus that closely resembles the extant human condition, but which also includes some specificities in local textural arrangement. The network of the inferolateral portion of the head shows a humanlike degree of anisotropy and a bone volume fraction intermediate between the extant human and the African ape patterns. These results suggest slight differences in gait kinematics between Pa. robustus and extant humans. The neck portion of the vertical bundle revealed a less biomechanically sensitive signal. Future investigations on the australopith hip joint loading environment should more carefully investigate the trabecular structure of the trochanteric region and possible structural covariation between cortical bone distribution across the neck and site-specific trabecular properties of the arcuate bundle.

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.

Internal architecture of the mandibular condyle of rabbits is related to dietary resistance during growth

Although there is considerable evidence that bone responds to the loading environment in which it develops, few analyses have examined phenotypic plasticity or bone functional adaptation in the masticatory apparatus. Prior work suggests that masticatory morphology is sensitive to differences in food mechanical properties during development; however, the importance of the timing/duration of loading and variation in naturalistic diets is less clear. Here, we examined microstructural and macrostructural differences in the mandibular condyle in four groups of white rabbits (Oryctolagus cuniculus) raised for a year on diets that varied in mechanical properties and timing of the introduction of mechanically challenging foods, simulating seasonal variation in diet. We employed sliding semilandmarks to locate multiple volumes of interest deep to the mandibular condyle articular surface, and compared bone volume fraction, trabecular thickness and spacing, and condylar size/shape among experimental groups. The results reveal a shared pattern of bony architecture across the articular surface of all treatment groups, while also demonstrating significant among-group differences. Rabbits raised on mechanically challenging diets have significantly increased bone volume fraction relative to controls fed a less challenging diet. The post-weaning timing of the introduction of mechanically challenging foods also influences architectural properties, suggesting that bone plasticity can extend well into adulthood and that bony responses to changes in loading may be rapid. These findings demonstrate that bony architecture of the mandibular condyle in rabbits responds to variation in mechanical loading during an organism's lifetime and has the potential to track dietary variation within and among species.

Ontogenetic changes to bone microstructure in an archaeologically derived sample of human ribs

There is considerable variation in the gross morphology and tissue properties among the bones of human infants, children, adolescents, and adults. Using 18 known-age individuals (n_female  = 8, n_male  = 9, n_unknown  = 1; birth to 21 years old), from a well-documented cemetery collection, Spitalfields Christ Church, London, UK, this study explores growth-related changes in cortical and trabecular bone microstructure. Micro-CT scans of mid-shaft middle thoracic ribs are used for quantitative analysis. Results are then compared to previously quantified conventional histomorphometry of the same sample. Total area (Tt.Ar), cortical area (Ct.Ar), cortical thickness (Ct.Th), and the major (Maj.Dm) and minor (Min.Dm) diameters of the rib demonstrate positive correlations with age. Pore density (Po.Dn) increases, but age-related changes to cortical porosity (Ct.Po) appear to be non-linear. Trabecular thickness (Tb.th) and trabecular separation (Tb.Sp) increase with age, whereas trabecular bone pattern factor (Tb.Pf), structural model index (SMI), and connectivity density (Conn.D) decrease with age. Sex-based differences were not identified for any of the variables included in this study. Some samples display clear evidence of diagenetic alteration without corresponding changes in radiopacity, which compromises the reliability of bone mineral density (BMD) data in the study of past populations. Cortical porosity data are not correlated with two-dimensional measures of osteon population density (OPD). This suggests that unfilled resorption spaces contribute more significantly to cortical porosity than do the Haversian canals of secondary osteons. Continued research using complementary imaging techniques and a wide array of histological variables will increase our understanding of age- and sex-specific ontogenetic patterns within and among human populations.

Trabecular architecture in the humeral metaphyses of non-avian reptiles (Crocodylia, Squamata and Testudines): Lifestyle, allometry and phylogeny

The lifestyle of extinct tetrapods is often difficult to assess when clear morphological adaptations such as swimming paddles are absent. According to the hypothesis of bone functional adaptation, the architecture of trabecular bone adapts sensitively to physiological loadings. Previous studies have already shown a clear relation between trabecular architecture and locomotor behavior, mainly in mammals and birds. However, a link between trabecular architecture and lifestyle has rarely been examined. Here, we analyzed trabecular architecture of different clades of reptiles characterized by a wide range of lifestyles (aquatic, amphibious, generalist terrestrial, fossorial, and climbing). Humeri of squamates, turtles, and crocodylians have been scanned with microcomputed tomography. We selected spherical volumes of interest centered in the proximal metaphyses and measured trabecular spacing, thickness and number, degree of anisotropy, average branch length, bone volume fraction, bone surface density, and connectivity density. Only bone volume fraction showed a significant phylogenetic signal and its significant difference between squamates and other reptiles could be linked to their physiologies. We found negative allometric relationships for trabecular thickness and spacing, positive allometries for connectivity density and trabecular number and no dependence with size for degree of anisotropy and bone volume fraction. The different lifestyles are well separated in the morphological space using linear discriminant analyses, but a cross-validation procedure indicated a limited predictive ability of the model. The trabecular bone anisotropy has shown a gradient in turtles and in squamates: higher values in amphibious than terrestrial taxa. These allometric scalings, previously emphasized in mammals and birds, seem to be valid for all amniotes. Discriminant analysis has offered, to some extent, a distinction of lifestyles, which however remains difficult to strictly discriminate. Trabecular architecture seems to be a promising tool to infer lifestyle of extinct tetrapods, especially those involved in the terrestrialization.

Trabecular Bone Structural Variation in the Proximal Sacrum Among Primates

The sacrum occupies a functionally important anatomical position as part of the pelvic girdle and vertebral column. Sacral orientation and external morphology in modern humans are distinct from those in other primates and compatible with the demands of habitual bipedal locomotion. Among nonhuman primates, however, how sacral anatomy relates to positional behaviors is less clear. As an alternative to evaluation of the sacrum's external morphology, this study assesses if the sacrum's internal morphology (i.e., trabecular bone) differs among extant primates. The primary hypothesis tested is that trabecular bone parameters with established functional relevance will differ in the first sacral vertebra (S1) among extant primates that vary in positional behaviors. Results for analyses of individual variables demonstrate that bone volume fraction, degree of anisotropy, trabecular number, and size-corrected trabecular thickness differ among primates grouped by positional behaviors to some extent, but not always in ways consistent with functional expectations. When examined as a suite, these trabecular parameters distinguish obligate bipeds from other positional behavior groups; and, the latter three trabecular bone variables further distinguish knuckle-walking terrestrial quadrupeds from manual suspensor-brachiators, vertical clingers and leapers, and arboreal quadrupeds, as well as between arboreal and terrestrial quadrupeds. As in other regions of the skeleton in modern humans, trabecular bone in S1 exhibits distinctively low bone volume fraction. Results from this study of extant primate S1 trabecular bone structural variation provide a functional context for interpretations concerning the positional behaviors of extinct primates based on internal sacral morphology. Anat Rec, 302:1354-1371, 2019. © 2018 Wiley Periodicals, Inc.

Trabecular architecture and joint loading of the proximal humerus in extant hominoids, Ateles, and Australopithecus africanus

OBJECTIVES

Several studies have investigated potential functional signals in the trabecular structure of the primate proximal humerus but with varied success. Here, we apply for the first time a "whole-epiphyses" approach to analysing trabecular bone in the humeral head with the aim of providing a more holistic interpretation of trabecular variation in relation to habitual locomotor or manipulative behaviors in several extant primates and Australopithecus africanus.

MATERIALS AND METHODS

We use a "whole-epiphysis" methodology in comparison to the traditional volume of interest (VOI) approach to investigate variation in trabecular structure and joint loading in the proximal humerus of extant hominoids, Ateles and A. africanus (StW 328).

RESULTS

There are important differences in the quantification of trabecular parameters using a "whole-epiphysis" versus a VOI-based approach. Variation in trabecular structure across knuckle-walking African apes, suspensory taxa, and modern humans was generally consistent with predictions of load magnitude and inferred joint posture during habitual behaviors. Higher relative trabecular bone volume and more isotropic trabeculae in StW 328 suggest A. africanus may have still used its forelimbs for arboreal locomotion.

DISCUSSION

A whole-epiphysis approach to analysing trabecular structure of the proximal humerus can help distinguish functional signals of joint loading across extant primates and can provide novel insight into habitual behaviors of fossil hominins.

Trabecular bone patterning across the human hand

Hand bone morphology is regularly used to link particular hominin species with behaviors relevant to cognitive/technological progress. Debates about the functional significance of differing hominin hand bone morphologies tend to rely on establishing phylogenetic relationships and/or inferring behavior from epigenetic variation arising from mechanical loading and adaptive bone modeling. Most research focuses on variation in cortical bone structure, but additional information about hand function may be provided through the analysis of internal trabecular structure. While primate hand bone trabecular structure is known to vary in ways that are consistent with expected joint loading differences during manipulation and locomotion, no study exists that has documented this variation across the numerous bones of the hand. We quantify the trabecular structure in 22 bones of the human hand (early/extant modern Homo sapiens) and compare structural variation between two groups associated with post-agricultural/industrial (post-Neolithic) and foraging/hunter-gatherer (forager) subsistence strategies. We (1) establish trabecular bone volume fraction (BV/TV), modulus (E), degree of anisotropy (DA), mean trabecular thickness (Tb.Th) and spacing (Tb.Sp); (2) visualize the average distribution of site-specific BV/TV for each bone; and (3) examine if the variation in trabecular structure is consistent with expected joint loading differences among the regions of the hand and between the groups. Results indicate similar distributions of trabecular bone in both groups, with those of the forager sample presenting higher BV/TV, E, and lower DA, suggesting greater and more variable loading during manipulation. We find indications of higher loading along the ulnar side of the forager sample hand, with high site-specific BV/TV distributions among the carpals that are suggestive of high loading while the wrist moves through the 'dart-thrower's' motion. These results support the use of trabecular structure to infer behavior and have direct implications for refining our understanding of human hand evolution and fossil hominin hand use.

Trabecular bone patterning in the hominoid distal femur

BACKGROUND

In addition to external bone shape and cortical bone thickness and distribution, the distribution and orientation of internal trabecular bone across individuals and species has yielded important functional information on how bone adapts in response to load. In particular, trabecular bone analysis has played a key role in studies of human and nonhuman primate locomotion and has shown that species with different locomotor repertoires display distinct trabecular architecture in various regions of the skeleton. In this study, we analyse trabecular structure throughout the distal femur of extant hominoids and test for differences due to locomotor loading regime.

METHODS

Micro-computed tomography scans of Homo sapiens (n = 11), Pan troglodytes (n = 18), Gorilla gorilla (n = 14) and Pongo sp. (n = 7) were used to investigate trabecular structure throughout the distal epiphysis of the femur. We predicted that bone volume fraction (BV/TV) in the medial and lateral condyles in Homo would be distally concentrated and more anisotropic due to a habitual extended knee posture at the point of peak ground reaction force during bipedal locomotion, whereas great apes would show more posteriorly concentrated BV/TV and greater isotropy due to a flexed knee posture and more variable hindlimb use during locomotion.

RESULTS

Results indicate some significant differences between taxa, with the most prominent being higher BV/TV in the posterosuperior region of the condyles in Pan and higher BV/TV and anisotropy in the posteroinferior region in Homo. Furthermore, trabecular number, spacing and thickness differ significantly, mainly separating Gorilla from the other apes.

DISCUSSION

The trabecular architecture of the distal femur holds a functional signal linked to habitual behaviour; however, there was more similarity across taxa and greater intraspecific variability than expected. Specifically, there was a large degree of overlap in trabecular structure across the sample, and Homo was not as distinct as predicted. Nonetheless, this study offers a comparative sample of trabecular structure in the hominoid distal femur and can contribute to future studies of locomotion in extinct taxa.

Systemic patterns of trabecular bone across the human and chimpanzee skeleton

Aspects of trabecular bone architecture are thought to reflect regional loading of the skeleton, and thus differ between primate taxa with different locomotor and postural modes. However, there are several systemic factors that affect bone structure that could contribute to, or be the primary factor determining, interspecific differences in bone structure. These systemic factors include differences in genetic regulation, sensitivity to loading, hormone levels, diet, and activity levels. Improved understanding of inter-/intraspecific variability, and variability across the skeleton of an individual, is required to interpret properly potential functional signals present within trabecular structure. Using a whole-region method of analysis, we investigated trabecular structure throughout the skeleton of humans and chimpanzees. Trabecular bone volume fraction (BV/TV), degree of anisotropy (DA) and trabecular thickness (Tb.Th) were quantified from high resolution micro-computed tomographic scans of the humeral and femoral head, third metacarpal and third metatarsal head, distal tibia, talus and first thoracic vertebra. We found that BV/TV is, in most anatomical sites, significantly higher in chimpanzees than in humans, suggesting a systemic difference in trabecular structure unrelated to local loading regime. Differences in BV/TV between the forelimb and hindlimb did not clearly reflect differences in locomotor loading in the study taxa. There were no clear systemic differences between the taxa in DA and, as such, this parameter might reflect function and relate to differences in joint loading. This systemic approach reveals both the pattern of variability across the skeleton and between taxa, and helps identify those features of trabecular structure that may relate to joint function.

Trabecular bone microarchitecture analysis, a way for an early detection of genetic dwarfism? Case study of a dwarf mother's offspring

A 66 year-old woman with a disproportionate dwarfism and who bore seven children was discovered at the Middenbeemster archaeological site (The Netherlands). Three are perinates and show no macroscopic or radiological evidence for a FGFR3 mutation causing hypo-or achondroplasia. This mutation induces dysfunction of the growth cartilage, leading to abnormalities in the development of trabecular bone. Because the mutation is autosomal dominant, these perinates have a 50% risk of having been affected. This study determines whether trabecular bone microarchitecture (TBMA) analysis is useful for detecting genetic dwarfism. Proximal metaphyses of humeri were μCT-scanned with a resolution of 7-12 μm. Three volumes of interest were segmented from each bone with TIVMI© software. The TBMA was quantified in BoneJ© using six parameters on which a multivariate analysis was then performed. Two of the Middenbeemster perinates show a quantitatively different TBMA organization. These results and the family's medical history suggest a diagnosis of genetic dwarfism for this two perinates. This study provides evidence to support the efficacy of μCT for diagnosing early-stage bone disease.

A novel use of 3D printing model demonstrates the effects of deteriorated trabecular bone structure on bone stiffness and strength

Trabecular bone structure is crucial to normal mechanical behavior of bones. Studies have shown that osteoporosis negatively affects trabecular bone structure, mainly by reducing bone volume fraction (BV/TV) and thus increasing fracture risk. One major limitation in assessing and quantifying the effect of this structural deterioration is that no two trabecular structures are identical. Thus, when we compare a group of healthy bones against a different group of bones that experienced resorption (i.e. decreased BV/TV) we only discover an "average" mechanical effect. It is impossible to quantify the mechanical effect of individual structural deterioration for each sample, simply because we never have the same sample in both states (intact and deteriorated structure). 3D printing is a new technology that can assist in overcoming this issue. Here we report a preliminary study that compares a healthy 3D printed trabecular bone model with the same model after bone resorption was simulated. Since the deteriorated structural bone model is derived from the healthy one, it is possible to directly estimate (percentage wise) the decrease of tissue stiffness and strength as a result of bone resorption for this specific structure. Our results demonstrate that a relatively small decrease in BV/TV (about 8%) leads to a dramatic decrease in structural strength (24%) and structural stiffness (17%), (P < 0.01). Structural strength decreased from an average of 9.14 ± 2.85MPa to 6.97 ± 2.44MPa, while structural stiffness decreased from an average of 282.5 ± 63.4N/mm to 233.8 ± 51.2N/mm. This study demonstrates that 3D printing is a novel and valuable tool for quantifying the effect of structural deterioration on the mechanical properties of trabecular bone. In the future, this approach may help us attain better personal fracture risk assessments by CT scanning, 3D printing and mechanically testing individual bone replicas from patients suffering excessive bone resorption.

A robust methodology for the quantitative assessment of the rat jawbone microstructure

Micro-computed tomography can be applied for the assessment of the micro-architectural characteristics of the cortical and trabecular bones in either physiological or disease conditions. However, reports often lack a detailed description of the methodological steps used to analyse these images, such as the volumes of interest, the algorithms used for image filtration, the approach used for image segmentation, and the bone parameters quantified, thereby making it difficult to compare or reproduce the studies. This study addresses this critical need and aims to provide standardized assessment and consistent parameter reporting related to quantitative jawbone image analysis. Various regions of the rat jawbones were screened for their potential for standardized micro-computed tomography analysis. Furthermore, the volumes of interest that were anticipated to be most susceptible to bone structural changes in response to experimental interventions were defined. In the mandible, two volumes of interest were selected, namely, the condyle and the trabecular bone surrounding the three molars. In the maxilla, the maxillary tuberosity region and the inter-radicular septum of the second molar were considered as volumes of interest. The presented protocol provides a standardized and reproducible methodology for the analysis of relevant jawbone volumes of interest and is intended to ensure global, accurate, and consistent reporting of its morphometry. Furthermore, the proposed methodology has potential, as a variety of rodent animal models would benefit from its implementation.

Trabecular and cortical bone structure of the talus and distal tibia in Pan and Homo

OBJECTIVES

Internal bone structure, both cortical and trabecular bone, remodels in response to loading and may provide important information regarding behavior. The foot is well suited to analysis of internal bone structure because it experiences the initial substrate reaction forces, due to its proximity to the substrate. Moreover, as humans and apes differ in loading of the foot, this region is relevant to questions concerning arboreal locomotion and bipedality in the hominoid fossil record.

MATERIALS AND METHODS

We apply a whole-bone/epiphysis approach to analyze trabecular and cortical bone in the distal tibia and talus of Pan troglodytes and Homo sapiens. We quantify bone volume fraction (BV/TV), degree of anisotropy (DA), trabecular thickness (Tb.Th), bone surface to volume ratio (BS/BV), and cortical thickness and investigate the distribution of BV/TV and cortical thickness throughout the bone/epiphysis.

RESULTS

We find that Pan has a greater BV/TV, a lower BS/BV and thicker cortices than Homo in both the talus and distal tibia. The trabecular structure of the talus is more divergent than the tibia, having thicker, less uniformly aligned trabeculae in Pan compared to Homo. Differences in dorsiflexion at the talocrural joint and in degree of mobility at the talonavicular joint are reflected in the distribution of cortical and trabecular bone.

DISCUSSION

Overall, quantified trabecular parameters represent overall differences in bone strength between the two species, however, DA may be directly related to joint loading. Cortical and trabecular bone distributions correlate with habitual joint positions adopted by each species, and thus have potential for interpreting joint position in fossil hominoids.

Trabecular mapping: Leveraging geometric morphometrics for analyses of trabecular structure

OBJECTIVES

Trabecular microstructure of limb bone epiphyses has been used to elucidate the relationship between skeletal form and behavior among mammals. Such studies have often relied on the analysis of a single volume of interest (VOI). Here we present a method for evaluating variation in bone microstructure across articular surfaces by leveraging sliding semilandmarks.

METHODS

Two samples were used to demonstrate the proposed methodology and test the hypothesis that microstructural variables are homogeneously distributed: tali from two ape genera (Pan and Pongo, n = 9) and modern human distal femora (n = 10). Sliding semilandmarks were distributed across articular surfaces and used to locate the position of multiple VOIs immediately deep to the cortical shell. Trabecular bone properties were quantified using the BoneJ plugin for ImageJ. Nonparametric MANOVA tests were used to make group comparisons and differences were explored using principal components analysis and visualized using color maps.

RESULTS

Tests reveal that trabecular parameters are not distributed homogeneously and identify differences between chimpanzee and orangutan tali with regards to trabecular spacing and degree of anisotropy, with chimpanzee tali being more anisotropic and having more uniformly spaced trabeculae. Human males and females differed in the pattern of trabecular spacing with males having more uniform trabecular spacing across the joint surface.

CONCLUSIONS

The proposed procedure quantifies variation in trabecular bone parameters across joint surfaces and allows for meaningful statistical comparisons between groups of interest. Consequently it holds promise to help elucidate links between trabecular bone structure and animal behavior.

Comparative analysis of trabecular bone structure and orientation in South African hominin tali

Tali of several hominin taxa are preserved in the fossil record and studies of the external morphology of these often show a mosaic of human-like and ape-like features. This has contributed to a growing recognition of variability characterizing locomotor kinematics of Australopithecus. In contrast, locomotor kinematics of another Plio-Pleistocene hominin, Paranthropus, are substantially less well-documented, in part, because of the paucity of postcranial fossils securely attributed to the genus. Since the talus transmits locomotor-based loads through the ankle and its internal structure is hypothesized to reflect accommodation to such loads, it is a cornerstone structure for reconstructing locomotor kinematics. Here we quantify and characterize trabecular bone morphology within tali attributed to Australopithecus africanus (StW 102, StW 363, StW 486) and Paranthropus robustus (TM 1517), making quantitative comparisons to modern humans, extant non-human apes, baboons, and a hominin talus attributed to Paranthropus boisei (KNM-ER 1464). Using high-resolution images of fossil tali (25 μm voxels), nine trabecular bone subregions of interest beneath the articular surface of the talar trochlea were segmented to quantify localized patterns in distribution and primary strut orientation. It was found that trabecular strut orientation and shape, in some cases, can discriminate amongst species characterized by different locomotor foot kinematics. Discriminant function analyses using standard trabecular bone structural properties align TM 1517 with Pan and Gorilla, while other hominin tali structurally most resemble those of baboons. In primary strut orientation, Paranthropus tali (KNM-ER 1464 and TM 1517) resemble the human condition in the anterior-medial subregion, where strut orientation appears positioned to distribute compressive loads medially and distally toward the talar head. In A. africanus tali (particularly StW 486), primary strut orientation in this region resembles that of apes. These results suggest that Paranthropus may have had a human-like medial weight shift during the last half of stance phase but Australopithecus did not.

The effects of cracks on the quantification of the cancellous bone fabric tensor in fossil and archaeological specimens: a simulation study

Cancellous bone is very sensitive to its prevailing mechanical environment, and study of its architecture has previously aided interpretations of locomotor biomechanics in extinct animals or archaeological populations. However, quantification of architectural features may be compromised by poor preservation in fossil and archaeological specimens, such as post mortem cracking or fracturing. In this study, the effects of post mortem cracks on the quantification of cancellous bone fabric were investigated through the simulation of cracks in otherwise undamaged modern bone samples. The effect on both scalar (degree of fabric anisotropy, fabric elongation index) and vector (principal fabric directions) variables was assessed through comparing the results of architectural analyses of cracked vs. non-cracked samples. Error was found to decrease as the relative size of the crack decreased, and as the orientation of the crack approached the orientation of the primary fabric direction. However, even in the best-case scenario simulated, error remained substantial, with at least 18% of simulations showing a > 10% error when scalar variables were considered, and at least 6.7% of simulations showing a > 10° error when vector variables were considered. As a 10% (scalar) or 10° (vector) difference is probably too large for reliable interpretation of a fossil or archaeological specimen, these results suggest that cracks should be avoided if possible when analysing cancellous bone architecture in such specimens.

Trabecular bone of precocials at birth; Are they prepared to run for the wolf(f)?

Bone is a dynamic tissue adapting to loading according to “Wolff's law of bone adaptation.” During very early life, however, such a mechanism may not be adequate enough to adapt to the dramatic change in environmental challenges in precocial species. Their neonates are required to stand and walk within hours after birth, in contrast to altricial animals that have much more time to adapt from the intrauterine environment to the outside world. In this study, trabecular bone parameters of the talus and sagittal ridge of the tibia from stillborn but full‐term precocials (calves and foals) were analyzed by micro‐CT imaging in order to identify possible anticipatory mechanisms to loading. Calculated average bone volume fraction in the Shetland pony (49–74%) was significantly higher compared to Warmblood foals (28–51%). Bovine trabecular bone was characterized by a low average bone volume fraction (22–28%), however, more directional anisotropy was found. It is concluded that anticipatory strategies in skeletal development exist in precocial species, which differ per species and are most likely related to anatomical differences in joint geometry and related loading patterns. The underlying regulatory mechanisms are still unknown, but they may be based on a genetic blueprint for the development of bone. More knowledge, both about a possible blueprint and its regulation, will be helpful in understanding developmental bone and joint diseases. J. Morphol. 277:948–956, 2016. © 2016 Wiley Periodicals, Inc.

A review of trabecular bone functional adaptation: what have we learned from trabecular analyses in extant hominoids and what can we apply to fossils?

Many of the unresolved debates in palaeoanthropology regarding evolution of particular locomotor or manipulative behaviours are founded in differing opinions about the functional significance of the preserved external fossil morphology. However, the plasticity of internal bone morphology, and particularly trabecular bone, allowing it to respond to mechanical loading during life means that it can reveal greater insight into how a bone or joint was used during an individual's lifetime. Analyses of trabecular bone have been commonplace for several decades in a human clinical context. In contrast, the study of trabecular bone as a method for reconstructing joint position, joint loading and ultimately behaviour in extant and fossil non-human primates is comparatively new. Since the initial 2D studies in the late 1970s and 3D analyses in the 1990 s, the utility of trabecular bone to reconstruct behaviour in primates has grown to incorporate experimental studies, expanded taxonomic samples and skeletal elements, and improved methodologies. However, this work, in conjunction with research on humans and non-primate mammals, has also revealed the substantial complexity inherent in making functional inferences from variation in trabecular architecture. This review addresses the current understanding of trabecular bone functional adaptation, how it has been applied to hominoids, as well as other primates and, ultimately, how this can be used to better interpret fossil hominoid and hominin morphology. Because the fossil record constrains us to interpreting function largely from bony morphology alone, and typically from isolated bones, analyses of trabecular structure, ideally in conjunction with that of cortical structure and external morphology, can offer the best resource for reconstructing behaviour in the past.

Trabecular architecture of the manual elements reflects locomotor patterns in primates

The morphology of trabecular bone has proven sensitive to loading patterns in the long bones and metacarpal heads of primates. It is expected that we should also see differences in the manual digits of primates that practice different methods of locomotion. Primate proximal and middle phalanges are load-bearing elements that are held in different postures and experience different mechanical strains during suspension, quadrupedalism, and knuckle walking. Micro CT scans of the middle phalanx, proximal phalanx and the metacarpal head of the third ray were used to examine the pattern of trabecular orientation in Pan, Gorilla, Pongo, Hylobates and Macaca. Several zones, i.e., the proximal ends of both phalanges and the metacarpal heads, were capable of distinguishing between knuckle-walking, quadrupedal, and suspensory primates. Orientation and shape seem to be the primary distinguishing factors but differences in bone volume, isotropy index, and degree of anisotropy were seen across included taxa. Suspensory primates show primarily proximodistal alignment in all zones, and quadrupeds more palmar-dorsal orientation in several zones. Knuckle walkers are characterized by having proximodistal alignment in the proximal ends of the phalanges and a palmar-dorsal alignment in the distal ends and metacarpal heads. These structural differences may be used to infer locmotor propensities of extinct primate taxa.

Trabecular bone structure in the primate wrist

Trabecular (or cancellous) bone has been shown to respond to mechanical loading throughout ontogeny and thus can provide unique insight into skeletal function and locomotion in comparative studies of living and fossil mammalian morphology. Trabecular bone of the hand may be particularly functionally informative because the hand has more direct contact with the substrate compared with the remainder of the forelimb during locomotion in quadrupedal mammals. This study investigates the trabecular structure within the wrist across a sample of haplorhine primates that vary in locomotor behaviour (and thus hand use) and body size. High-resolution microtomographic scans were collected of the lunate, scaphoid, and capitate in 41 individuals and eight genera (Homo, Gorilla, Pan, Papio, Pongo, Symphalangus, Hylobates, and Ateles). We predicted that particular trabecular parameters would 1) vary across suspensory, quadrupedal, and bipedal primates based on differences in hand use and load, and 2) scale with carpal size following similar allometric patterns found previously in other skeletal elements across a larger sample of mammals and primates. Analyses of variance (trabecular parameters analysed separately) and principal component analyses (trabecular parameters analysed together) revealed no clear functional signal in the trabecular structure of any of the three wrist bones. Instead, there was a large degree of variation within suspensory and quadrupedal locomotor groups, as well as high intrageneric variation within some taxa, particularly Pongo and Gorilla. However, as predicted, Homo sapiens, which rarely use their hands for locomotion and weight support, were unique in showing lower relative bone volume (BV/TV) compared with all other taxa. Furthermore, parameters used to quantify trabecular structure within the wrist scale with size generally following similar allometric patterns found in trabeculae of other mammalian skeletal elements. We discuss the challenges associated with quantifying and interpreting trabecular bone within the wrist.

Trabecular bone structure correlates with hand posture and use in hominoids

Bone is capable of adapting during life in response to stress. Therefore, variation in locomotor and manipulative behaviours across extant hominoids may be reflected in differences in trabecular bone structure. The hand is a promising region for trabecular analysis, as it is the direct contact between the individual and the environment and joint positions at peak loading vary amongst extant hominoids. Building upon traditional volume of interest-based analyses, we apply a whole-epiphysis analytical approach using high-resolution microtomographic scans of the hominoid third metacarpal to investigate whether trabecular structure reflects differences in hand posture and loading in knuckle-walking (Gorilla, Pan), suspensory (Pongo, Hylobates and Symphalangus) and manipulative (Homo) taxa. Additionally, a comparative phylogenetic method was used to analyse rates of evolutionary changes in trabecular parameters. Results demonstrate that trabecular bone volume distribution and regions of greatest stiffness (i.e., Young's modulus) correspond with predicted loading of the hand in each behavioural category. In suspensory and manipulative taxa, regions of high bone volume and greatest stiffness are concentrated on the palmar or distopalmar regions of the metacarpal head, whereas knuckle-walking taxa show greater bone volume and stiffness throughout the head, and particularly in the dorsal region; patterns that correspond with the highest predicted joint reaction forces. Trabecular structure in knuckle-walking taxa is characterised by high bone volume fraction and a high degree of anisotropy in contrast to the suspensory brachiators. Humans, in which the hand is used primarily for manipulation, have a low bone volume fraction and a variable degree of anisotropy. Finally, when trabecular parameters are mapped onto a molecular-based phylogeny, we show that the rates of change in trabecular structure vary across the hominoid clade. Our results support a link between inferred behaviour and trabecular structure in extant hominoids that can be informative for reconstructing behaviour in fossil primates.

Finite element micro-modelling of a human ankle bone reveals the importance of the trabecular network to mechanical performance: new methods for the generation and comparison of 3D models

Most modelling of whole bones does not incorporate trabecular geometry and treats bone as a solid non-porous structure. Some studies have modelled trabecular networks in isolation. One study has modelled the performance of whole human bones incorporating trabeculae, although this required considerable computer resources and purpose-written code. The difference between mechanical behaviour in models that incorporate trabecular geometry and non-porous models has not been explored. The ability to easily model trabecular networks may shed light on the mechanical consequences of bone loss in osteoporosis and remodelling after implant insertion. Here we present a Finite Element Analysis (FEA) of a human ankle bone that includes trabecular network geometry. We compare results from this model with results from non-porous models and introduce protocols achievable on desktop computers using widely available softwares. Our findings show that models including trabecular geometry are considerably stiffer than non-porous whole bone models wherein the non-cortical component has the same mass as the trabecular network, suggesting inclusion of trabecular geometry is desirable. We further present new methods for the construction and analysis of 3D models permitting: (1) construction of multi-property, non-porous models wherein cortical layer thickness can be manipulated; (2) maintenance of the same triangle network for the outer cortical bone surface in both 3D reconstruction and non-porous models allowing exact replication of load and restraint cases; and (3) creation of an internal landmark point grid allowing direct comparison between 3D FE Models (FEMs).

Methodological considerations for analyzing trabecular architecture: an example from the primate hand

Micro-computed tomographic analyses of trabecular bone architecture have been used to clarify the link between positional behavior and skeletal anatomy in primates. However, there are methodological decisions associated with quantifying and comparing trabecular anatomy across taxa that vary greatly in body size and morphology that can affect characterizations of trabecular architecture, such as choice of the volume of interest (VOI) size and location. The potential effects of these decisions may be amplified in small, irregular-shaped bones of the hands and feet that have more complex external morphology and more heterogeneous trabecular structure compared to, for example, the spherical epiphysis of the femoral head. In this study we investigate the effects of changes in VOI size and location on standard trabecular parameters in two bones of the hand, the capitate and third metacarpal, in a diverse sample of nonhuman primates that vary greatly in morphology, body mass and positional behavior. Results demonstrate that changes in VOI location and, to a lesser extent, changes in VOI size had a dramatic affect on many trabecular parameters, especially trabecular connectivity and structure (rods vs. plates), degree of anisotropy, and the primary orientation of the trabeculae. Although previous research has shown that some trabecular parameters are susceptible to slight variations in methodology (e.g. VOI location, scan resolution), this study provides a quantification of these effects in hand bones of a diverse sample of primates. An a priori understanding of the inherent biases created by the choice of VOI size and particularly location is critical to robust trabecular analysis and functional interpretation, especially in small bones with complex arthroses.