Trabecular mapping: Leveraging geometric morphometrics for analyses of trabecular structure

Trabecular bone variation in the gorilla calcaneus

OBJECTIVES

Calcaneal external shape differs among nonhuman primates relative to locomotion. Such relationships between whole-bone calcaneal trabecular structure and locomotion, however, have yet to be studied. Here we analyze calcaneal trabecular architecture in Gorilla gorilla gorilla, Gorilla beringei beringei, and G. b. graueri to investigate general trends and fine-grained differences among gorilla taxa relative to locomotion.

MATERIALS AND METHODS

Calcanei were micro-CT scanned. A three-dimensional geometric morphometric sliding semilandmark analysis was carried out and the final landmark configurations used to position 156 volumes of interest. Trabecular thickness (Tb.Th), trabecular spacing (Tb.Sp), and bone volume fraction (BV/TV) were calculated using the BoneJ plugin for ImageJ and MATLAB. Non-parametric MANOVAs were run to test for significant differences among taxa in parameter raw values and z-scores. Parameter distributions were visualized using color maps and summarized using principal components analysis.

RESULTS

There are no significant differences in raw BV/TV or Tb.Th among gorillas, however G. b. beringei significantly differs in z-scores for both parameters (p = <0.0271). All three taxa exhibit relatively lower BV/TV and Tb.Th in the posterior half of the calcaneus. This gradation is exacerbated in G. b. beringei. G. b. graueri significantly differs from other taxa in Tb.Sp z-scores (p < 0.001) indicating a different spacing distribution.

DISCUSSION

Relatively higher Tb.Th and BV/TV in the anterior calcaneus among gorillas likely reflects higher forces associated with body mass (transmitted through the subtalar joint) relative to forces transferred through the posterior calcaneus. The different Tb.Sp pattern in G. b. graueri may reflect proposed differences in foot positioning during locomotion.

Trabecular architecture of the distal femur in extant hominids

Extant great apes are characterized by a wide range of locomotor, postural and manipulative behaviours that each require the limbs to be used in different ways. In addition to external bone morphology, comparative investigation of trabecular bone, which (re-)models to reflect loads incurred during life, can provide novel insights into bone functional adaptation. Here, we use canonical holistic morphometric analysis (cHMA) to analyse the trabecular morphology in the distal femoral epiphysis of Homo sapiens (n = 26), Gorilla gorilla (n = 14), Pan troglodytes (n = 15) and Pongo sp. (n = 9). We test two predictions: (1) that differing locomotor behaviours will be reflected in differing trabecular architecture of the distal femur across Homo, Pan, Gorilla and Pongo; (2) that trabecular architecture will significantly differ between male and female Gorilla due to their different levels of arboreality but not between male and female Pan or Homo based on previous studies of locomotor behaviours. Results indicate that trabecular architecture differs among extant great apes based on their locomotor repertoires. The relative bone volume and degree of anisotropy patterns found reflect habitual use of extended knee postures during bipedalism in Homo, and habitual use of flexed knee posture during terrestrial and arboreal locomotion in Pan and Gorilla. Trabecular architecture in Pongo is consistent with a highly mobile knee joint that may vary in posture from extension to full flexion. Within Gorilla, trabecular architecture suggests a different loading of knee in extension/flexion between females and males, but no sex differences were found in Pan or Homo, supporting our predictions. Inter- and intra-specific variation in trabecular architecture of distal femur provides a comparative context to interpret knee postures and, in turn, locomotor behaviours in fossil hominins.

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.

Effects of obesity on talar micro- and macro-morphology

The addition of information regarding obesity status to the forensic anthropological biological profile could significantly contribute to the identification of human skeletal remains since over 40% of the U.S. adult population is currently obese. This study examines the differences in talar shape and trabecular bone structure between obese and non-obese individuals. A sample of 20 obese and 20 non-obese divided evenly by sex was selected from the Texas State University Donated Skeletal Collection. Tali were imaged using x-ray computed tomography (voxel size: 28-38.7 μm). Image stacks were processed to produce binary images as well as trabecular thickness and spacing maps. Landmark-based geometric morphometric analyses were conducted to quantify shape variation. Shape coordinates were used to locate 100 geometrically homologous volumes of interest within each talus. Bone volume fraction, trabecular thickness, and trabecular spacing were extracted at each volume of interest. Within each sex, a one-way ANCOVA was used to determine if significant differences exist between obese and non-obese individuals in trabecular bone after controlling for age. The size of the talus as well as subtle aspects of shape were found to distinguish the sexes. The results further indicate that bone volume fraction significantly differs between obese and non-obese males. In females, bone volume fraction is correlated with age but does not differ between obese and non-obese. The study demonstrates that bone microstructure is a promising approach to estimating body mass or body mass index category but age effects diminish the potential for the talus to be used alone.

A computational framework for canonical holistic morphometric analysis of trabecular bone

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

Integrative Approach Uncovers New Patterns of Ecomorphological Convergence in Slow Arboreal Xenarthrans

Identifying ecomorphological convergence examples is a central focus in evolutionary biology. In xenarthrans, slow arboreality independently arose at least three times, in the two genera of ‘tree sloths’, Bradypus and Choloepus, and the silky anteater, Cyclopes. This specialized locomotor ecology is expectedly reflected by distinctive morpho-functional convergences. Cyclopes, although sharing several ecological features with ‘tree sloths’, do not fully mirror the latter in their outstandingly similar suspensory slow arboreal locomotion. We hypothesized that the morphology of Cyclopes is closer to ‘tree sloths’ than to anteaters, but yet distinct, entailing that slow arboreal xenarthrans evolved through ‘incomplete’ convergence. In a multivariate trait space, slow arboreal xenarthrans are hence expected to depart from their sister taxa evolving toward the same area, but not showing extensive phenotypical overlap, due to the distinct position of Cyclopes. Conversely, a pattern of ‘complete’ convergence (i.e., widely overlapping morphologies) is hypothesized for ‘tree sloths’. Through phylogenetic comparative methods, we quantified humeral and femoral convergence in slow arboreal xenarthrans, including a sample of extant and extinct non-slow arboreal xenarthrans. Through 3D geometric morphometrics, cross-sectional properties (CSP) and trabecular architecture, we integratively quantified external shape, diaphyseal anatomy and internal epiphyseal structure. Several traits converged in slow arboreal xenarthrans, especially those pertaining to CSP. Phylomorphospaces and quantitative convergence analyses substantiated the expected patterns of ‘incomplete’ and ‘complete’ convergence for slow arboreal xenarthrans and ‘tree sloths’, respectively. This work, highlighting previously unidentified convergence patterns, emphasizes the value of an integrative multi-pronged quantitative approach to cope with complex mechanisms underlying ecomorphological convergence.

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.

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.

Trabecular variation in the first metacarpal and manipulation in hominids

OBJECTIVES

The dexterity of fossil hominins is often inferred by assessing the comparative manual anatomy and behaviors of extant hominids, with a focus on the thumb. The aim of this study is to test whether trabecular structure is consistent with what is currently known about habitually loaded thumb postures across extant hominids.

MATERIALS AND METHODS

We analyze first metacarpal (Mc1) subarticular trabecular architecture in humans (Homo sapiens, n = 10), bonobos (Pan paniscus, n = 10), chimpanzees (Pan troglodytes, n = 11), as well as for the first time, gorillas (Gorilla gorilla gorilla, n = 10) and orangutans (Pongo sp., n = 1, Pongo abelii, n = 3 and Pongo pygmaeus, n = 5). Using a combination of subarticular and whole-epiphysis approaches, we test for significant differences in relative trabecular bone volume (RBV/TV) and degree of anisotropy (DA) between species.

RESULTS

Humans have significantly greater RBV/TV on the radiopalmar aspects of both the proximal and distal Mc1 subarticular surfaces and greater DA throughout the Mc1 head than other hominids. Nonhuman great apes have greatest RBV/TV on the ulnar aspect of the Mc1 head and the palmar aspect of the Mc1 base. Gorillas possessed significantly lower DA in the Mc1 head than any other taxon in our sample.

DISCUSSION

These results are consistent with abduction of the thumb during forceful "pad-to-pad" precision grips in humans and, in nonhuman great apes, a habitually adducted thumb that is typically used in precision and power grips. This comparative context will help infer habitual manipulative and locomotor grips in fossil hominins.

Metacarpal trabecular bone varies with distinct hand-positions used in hominid locomotion

Trabecular bone remodels during life in response to loading and thus should, at least in part, reflect potential variation in the magnitude, frequency and direction of joint loading across different hominid species. Here we analyse the trabecular structure across all non-pollical metacarpal distal heads (Mc2-5) in extant great apes, expanding on previous volume of interest and whole-epiphysis analyses that have largely focused on only the first or third metacarpal. Specifically, we employ both a univariate statistical mapping and a multivariate approach to test for both inter-ray and interspecific differences in relative trabecular bone volume fraction (RBV/TV) and degree of anisotropy (DA) in Mc2-5 subchondral trabecular bone. Results demonstrate that whereas DA values only separate Pongo from African apes (Pan troglodytes, Pan paniscus, Gorilla gorilla), RBV/TV distribution varies with the predicted loading of the metacarpophalangeal (McP) joints during locomotor behaviours in each species. Gorilla exhibits a relatively dorsal distribution of RBV/TV consistent with habitual hyper-extension of the McP joints during knuckle-walking, whereas Pongo has a palmar distribution consistent with flexed McP joints used to grasp arboreal substrates. Both Pan species possess a disto-dorsal distribution of RBV/TV, compatible with multiple hand postures associated with a more varied locomotor regime. Further inter-ray comparisons reveal RBV/TV patterns consistent with varied knuckle-walking postures in Pan species in contrast to higher RBV/TV values toward the midline of the hand in Mc2 and Mc5 of Gorilla, consistent with habitual palm-back knuckle-walking. These patterns of trabecular bone distribution and structure reflect different behavioural signals that could be useful for determining the behaviours of fossil hominins.

Trabecular architecture of the great ape and human femoral head

Studies of femoral trabecular structure have shown that the orientation and volume of bone are associated with variation in loading and could be informative about individual joint positioning during locomotion. In this study, we analyse for the first time trabecular bone patterns throughout the femoral head using a whole-epiphysis approach to investigate how potential trabecular variation in humans and great apes relates to differences in locomotor modes. Trabecular architecture was analysed using microCT scans of Pan troglodytes (n = 20), Gorilla gorilla (n = 14), Pongo sp. (n = 5) and Homo sapiens (n = 12) in medtool 4.1. Our results revealed differences in bone volume fraction (BV/TV) distribution patterns, as well as overall trabecular parameters of the femoral head between great apes and humans. Pan and Gorilla showed two regions of high BV/TV in the femoral head, consistent with hip posture and loading during two discrete locomotor modes: knuckle-walking and climbing. Most Pongo specimens also displayed two regions of high BV/TV, but these regions were less discrete and there was more variability across the sample. In contrast, Homo showed only one main region of high BV/TV in the femoral head and had the lowest BV/TV, as well as the most anisotropic trabeculae. The Homo trabecular structure is consistent with stereotypical loading with a more extended hip compared with great apes, which is characteristic of modern human bipedalism. Our results suggest that holistic evaluations of femoral head trabecular architecture can reveal previously undetected patterns linked to locomotor behaviour in extant apes and can provide further insight into hip joint loading in fossil hominins and other primates.

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.